gpq
gpq copied to clipboard
Specify the dst crs in convert?
I don't think gpq currently contains a method to specify the target crs. Also I see that by default you use "OGC:CRS84", what is your rationale for that? Why not, for example, use "EPSG:4326"?
I'll add a little bit of context on my use case. So I just used gpq
to convert a 'big' collection of parquet files to geoparquet by simply doing gpq convert non-geo.parquet valid-geo.parquet
in a for loop. Further in my processing chain I load these geoparquet files using GeoPandas
, but I ran into an issue because when the crs == "OGC:CRS84"
it cannot be converted to epgs. Although it's expected behaviour I'm mostly just curious why you use "OGC:CRS84" instead of "EPSG:4326".
gdf = gpd.read_parquet("valid-geo.parquet")
print(gdf.crs.to_epsg()) # None
print(gdf.to_crs(4326).to_epsg()) # 4326
I'll probably change my routines from gdf.crs.to_epsg()
to gdf.crs.to_string()
, but I guess that several others rely on to_epsg() as well when using GeoPandas, so I thought it's worth opening a discussion point here.