从点生成geopandas数据帧，然后投影gpd，得到错误：无法转换原始几何体。请先在对象上设置crs

Question

我将一个具有lat/lon的数据集转换为geopandas数据帧（我们称之为gpd）。gpd没有CRS，我试图使用以下方法将此gpd投影到“EPSG:3857”：

gpd=对crs的gpd（“EPSG:3857”）

我得到一个错误，说“无法转换原始几何体。请先在对象上设置一个CRS”。这是否意味着，对于任何没有CRS的gpd，我必须首先为其分配CRS，然后将其重新投影到我感兴趣的CRS，例如，从4326到3857：

gpd.crs=“EPSG:4326”

gpd=对crs的gpd（“EPSG:3857”）

我的困惑是CRS和投影之间的区别，为什么我不能在这种情况下直接使用3857。此外，pyproj项目提到重新投影geodata有一些变化，并就如何投影geodata给出了建议，我不清楚细节，有人能就如何适当投影/重新投影geodata给出一些建议吗？谢谢

Answer 1

您需要先定义crs，然后才能对其进行投影。如果.shp文件没有附带的.prj文件，则通常会发生这种情况。

import geopandas as gpd
gdf = gpd.read_file('GSHHS_c_L1.shp')
print(gdf.crs)
None

#this is how we define the projection
gdf.crs = "EPSG:4326"
print(gdf.crs)
EPSG:4326

#In Jupyter if you don't use a print statement, but shift+enter, you'll get this for your crs

<Geographic 2D CRS: EPSG:4326>
Name: WGS 84
Axis Info [ellipsoidal]:
- Lat[north]: Geodetic latitude (degree)
- Lon[east]: Geodetic longitude (degree)
Area of Use:
- name: World
- bounds: (-180.0, -90.0, 180.0, 90.0)
Datum: World Geodetic System 1984
- Ellipsoid: WGS 84
- Prime Meridian: Greenwich

定义了crs之后，就可以对其进行投影了。

#this is how we project the data
gdf = gdf.to_crs("EPSG:3857")
# this will display crs information in jupyter
gdf.crs

<Projected CRS: EPSG:3857>
Name: WGS 84 / Pseudo-Mercator
Axis Info [cartesian]:
- X[east]: Easting (metre)
- Y[north]: Northing (metre)
Area of Use:
- name: World - 85°S to 85°N
- bounds: (-180.0, -85.06, 180.0, 85.06)
Coordinate Operation:
- name: Popular Visualisation Pseudo-Mercator
- method: Popular Visualisation Pseudo Mercator
Datum: World Geodetic System 1984
- Ellipsoid: WGS 84
- Prime Meridian: Greenwich