GeoPandas – Using Specific Transform on ‘to_crs’ Function

geopandaspyprojpython

I have a GeoDataFrame with CRS epsg:4326 and want to convert to epsg:32066

Such as:

gdf.to_crs(epsg=32066)

I know that GeoPandas is selecting the optimal transformation and it works.
But I want to specify which Transformation code it should use. Is this possible?

I have attempted to do this with Pyproj by specifying the transformation:

transform_step1 = Transformer.from_pipeline("EPSG:15851")
transform_step1.transform(-84.81629793076491, 24.233538817965595, direction="INVERSE")

But I want to know if it is possible to specify the transformation code in the .to_crs() method.

Best Answer

It is not currently possible, but has been discussed in the past: https://github.com/geopandas/geopandas/issues/1175

Related Solutions

GeoPandas CRS – to_crs init proj4 String Works but New Style Doesn’t

I reported this as a bug here: https://github.com/pyproj4/pyproj/issues/549

In the meantime, you can keep using the deprecated "{'init':'epsg:4326'}" syntax a bit longer (and ignore the warning for now).
Alternatively, you can use the following workaround to read the CRS directly from the .prj file (avoiding going through GDAL):

gdf = geopandas.read_file('note/map/bnd_emd_2012/BND_ADM_EMD_2012.shp')    

with open("note/map/bnd_emd_2012/BND_ADM_EMD_2012.prj") as prj_file:
    crs = prj_file.read()

gdf.crs = crs

gdf2 = gdf.to_crs('epsg:4326')

Pyproj – How to Use Specific Transformation in Pyproj With Python

>>> from pyproj import Transformer
>>> tr = Transformer.from_pipeline("EPSG:15851")
>>> tr
<Transformation Transformer: pipeline>
Description: NAD27 to WGS 84 (79)
Area of Use:
- name: United States (USA) - CONUS including EEZ -onshore and offshore - Alabama; Arizona; Arkansas; California; Colorado; Connecticut; Delaware; Florida; Georgia; Idaho; Illinois; Indiana; Iowa; Kansas; Kentucky; Louisiana; Maine; Maryland; Massachusetts; Michigan; Minnesota; Mississippi; Missouri; Montana; Nebraska; Nevada; New Hampshire; New Jersey; New Mexico; New York; North Carolina; North Dakota; Ohio; Oklahoma; Oregon; Pennsylvania; Rhode Island; South Carolina; South Dakota; Tennessee; Texas; Utah; Vermont; Virginia; Washington; West Virginia; Wisconsin; Wyoming. US Gulf of Mexico (GoM) OCS.
- bounds: (-129.17, 23.81, -65.69, 49.38)

Best Answer

Related Solutions

GeoPandas CRS – to_crs init proj4 String Works but New Style Doesn’t

Pyproj – How to Use Specific Transformation in Pyproj With Python

Related Question