我有一个geopandas dataframe,其Crs为EPSG:4326,我想转换为ESPG:102001的CRC,该CRC的定义如下:
:proj4_102001 = '+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs'
在我的/usr/share/proj目录中,该投影位于esri.extra文件中,但在proj_def.dat文件中不在下面的屏幕截图中看到:
1)GREP命令的结果搜索102001和2)/usr/share/proj目录的LS
根据Geopandas文档,定义投影的正确方法如下:
world = world.to_crs({'init': 'epsg:3395'})
我尝试过" world = world.to_crs({'init':'esri.extra:102001'})"和许多其他变体,但没有任何变化。我尝试在互联网上搜索包括ProJ.4文档在内的解决方案,但尚未找到任何解决方案。如果您知道该怎么做,请告诉我。我正在使用Ubuntu(Ubuntu-18.04),Python(3.7.1),Geopandas(0.4.0)和Proj.4(5.2.0)。
预先感谢您。
tom
您应该能够在to_crs
函数中提供完整的投影字符串。
world = world.to_crs('+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs')