Mark Litwintschik

```sql D select ST_GEOMETRYTYPE(geom), count(*) from roads2 group by 1 order by 1; ┌───────────────────────┬──────────────┐ │ st_geometrytype(geom) │ count_star() │ │ "geometry_type" │ int64 │ ├───────────────────────┼──────────────┤ │ LINESTRING │ 4628534 │...

```sql D select ST_IsValid(geom), count(*) from roads2 group by 1 order by 1; ┌──────────────────┬──────────────┐ │ st_isvalid(geom) │ count_star() │ │ boolean │ int64 │ ├──────────────────┼──────────────┤ │ true │ 4628534 │...

Just a heads up, I tried running the steps again on the latest master branch using the steps from the blog post on my MBP and it ran out of...

The script is meant to be run from Ubuntu for Windows, not any of the other Windows shells. I'll research pure python geojson to gpkg libraries like https://github.com/geomet/geomet which could...

Seems that fastparquet doesn't spot the issue. ```python df = pd.read_parquet("sample.parquet", engine='pyarrow') # Exception df = pd.read_parquet("sample.parquet", engine='fastparquet') # Runs fine ```

Understood. I'll look into that. What do you mean by RPC?

I asked Umbra for advice on a set of gdal parameters https://twitter.com/marklit82/status/1809151161825427959

A single call to ``gdalwarp`` fixed it. ```bash $ wget https://umbra-open-data-catalog.s3.amazonaws.com/sar-data/tasks/ad%20hoc/Haneda_Japan/0028c361-193c-49b7-ad59-808033d5dc74/2024-02-27-01-18-51_UMBRA-07/2024-02-27-01-18-51_UMBRA-07_CSI.tif $ gdalwarp -t_srs "EPSG:4326" 2024-02-27-01-18-51_UMBRA-07_CSI.tif warped.tif $ gdal_retile.py \ -s_srs "EPSG:4326" \ -ps 4096 4096 \ -targetDir ./ \...