GeoPandas .to_crs() 不能作为函数的一部分工作，但可以在独立程序中使用答案

【问题标题】：GeoPandas .to_crs() won't work as part of a function but does in a standalone programGeoPandas .to_crs() 不能作为函数的一部分工作，但可以在独立程序中使用
【发布时间】：2021-09-24 12:36:06
【问题描述】：

我正在尝试将单个 GeoJSON 层的点从 EPSG:4326 重新投影到 EPSG:3857。当我在一个独立的程序（如下）中单独运行以下代码时，输出符合预期。

geojsonSaveDir = 'path_to_input_file'
firePoints4326 = gpd.read_file(geojsonSaveDir)
firePointsReproject = firePoints4326.copy()
firePointsReproject['geometry'] = firePointsReproject['geometry'].to_crs(3857)
print(firePoints4326['geometry'].head())
print(firePointsReproject['geometry'].head())
firePointsReproject.crs = from_epsg(3857)
reprojectSaveDir = geojsonSaveDir.replace('.geojson', '') + '_3857'
firePointsReproject.to_file(reprojectSaveDir)

以上代码的输出为：

0    POINT (-61.56000 -11.96000)
1    POINT (-61.58000 -11.96000)
2    POINT (-59.21000 -11.66000)
3    POINT (-59.23000 -11.66000)
4    POINT (-51.55000 -11.28000)
Name: geometry, dtype: geometry
1    POINT (-6855054.243 -1341156.488)
2    POINT (-6591227.050 -1307038.377)
3    POINT (-6593453.440 -1307038.377)
4    POINT (-5738519.750 -1263874.866)
Name: geometry, dtype: geometry

但是，当我使用完全相同的代码创建函数时，输出是 EPSG:3857 中的 shapefile，其中所有点的坐标均为 0.00000、0.00000。大多数情况下，这甚至不会显示在 QGIS 中。此函数的输入参数是要重新投影的 GeoJSON 文件的路径。我对输出是 shapefile 很好，但只是无法弄清楚为什么重新投影不会像以前那样工作。

def pointReprojector(geojsonSaveDir):
firePoints4326 = gpd.read_file(geojsonSaveDir)
firePointsReproject = firePoints4326.copy()
firePointsReproject['geometry'] = firePointsReproject['geometry'].to_crs(epsg=3857)
print(firePoints4326['geometry'].head())
print(firePointsReproject['geometry'].head())
firePointsReproject.crs = from_epsg(3857)
reprojectSaveDir = geojsonSaveDir.replace('.geojson', '') + '_3857'
firePointsReproject.to_file(reprojectSaveDir)
return reprojectSaveDir

这个的输出是：

0    POINT (-61.56000 -11.96000)
1    POINT (-61.58000 -11.96000)
2    POINT (-59.21000 -11.66000)
3    POINT (-59.23000 -11.66000)
4    POINT (-51.55000 -11.28000)
Name: geometry, dtype: geometry
0    POINT (0.00000 0.00000)
1    POINT (0.00000 0.00000)
2    POINT (0.00000 0.00000)
3    POINT (0.00000 0.00000)
4    POINT (0.00000 0.00000)
Name: geometry, dtype: geometry

我对 GeoPandas 还很陌生，因此我们将不胜感激。如果需要更多信息，请告诉我。

【问题讨论】：

标签： python geojson geopandas

【解决方案1】：

尝试改用这个：

def pointReprojector(geojsonSaveDir):
    firePoints4326 = gpd.read_file(geojsonSaveDir)
    firePointsReproject = firePoints4326.copy()
    firePointsReproject['geometry'] = firePointsReproject['geometry'].to_crs('epsg:3857')
    print(firePoints4326['geometry'].head())
    print(firePointsReproject['geometry'].head())
    #firePointsReproject.crs = from_epsg(3857)
    reprojectSaveDir = geojsonSaveDir.replace('.geojson', '') + '_3857'
    firePointsReproject.to_file(reprojectSaveDir)
    return reprojectSaveDir

请注意，重新投影后无需设置 CRS。这就是为什么我评论了这条线。另外，我看到您为to_crs 函数使用了两种不同样式的输入（to_crs(3857) 和.to_crs(epsg=3857)）。只要坚持这种输入方式就不会出错：to_crs("epsg:3857")。

我用我下载的 GeoJSON 尝试了上面的代码，结果很好 - 生成的 shapefile 在 QGIS 中非常完美。

【讨论】：