在Geopandas中更改df CRS后向点添加坐标

Question

基本上，这个想法是在不使用Qgis的情况下自动化一些工作流。

我无法获得与Geopandas中Qgis功能“向点添加坐标”类似的结果，该功能允许您在当前投影中获取点的x、y坐标，并为表创建新属性。

所以我有一套我玩过的积分。原始形状文件的CRS为epsg 2154（Lambert 93）。我需要以与谷歌地图兼容的格式获取纬度和经度。

谷歌使用epsg 3857谷歌地图。

points=pd.DataFrame（{'Id'：['001'，'002'，'003'，'geometry'：['POINT（909149.3986619939 6881986.232659903）'，'POINT（909649.3986619939 6882486.232659903）'，'POINT（909149.3986619939 6882486.232659903）'）

这个想法是切换到epsg 3857（wgs84），并从那里创建lat/long列，这些列将填充wgs84坐标，例如47.357955,1.7317783。

因此，我所做的是明显地改变CRS:pointswgs84=points.to_CRS（espg=3857）

然后

pointswgs84['lon']=pointswgs84.geometry.x

pointswgs84['lat']=pointswgs84.geometry.y

但是我的lat/long列会被对应于原始点数据框的坐标填充：

points=pd.DataFrame（{'Id'：['001'，'002'，'003'，'geometry'：['POINT（909149.3986619939 6881986.232659903）'，'POINT（909649.3986619939 6882486.232659903）'，'long'：['689149.3986619939 6882486.232659903'，'6882486.232659903'，'lat'：['909149.3986619999'，'909649.3986619999'，'909149.3986619999']}）

看起来我在这里遗漏了一些东西，但因为我是Python新手

谢谢你的帮助。

Answer 1

好了，问题解决了。这是一个简单的espg问题。如果这能帮助像我这样的新手，下面是一个解决方案：

即使谷歌使用espg 3857，如果你需要将你的点的X，Y转换为可以利用的东西（如网址参数），你也需要使用espg 4326。

然后你会得到你的长lat正确的方式。如5.859162,49.003710.

因此：

步骤1：将您的原始shapefile转换为espg 4326

points.to_crs（{'init'：'epsg： 4326'}）

步骤2：只需提取x和y并将它们添加到新列中

pointswgs84['lon']=pointswgs84.geometry.x

pointswgs84['lat']=pointswgs84.geometry.y

Answer 2

您可以尝试：

import geopandas as gpd
from shapely.geometry import Point

points = pd.DataFrame({'Id': ['001', '002', '003'],
                       'geometry': [Point(909149.3986619939, 6881986.232659903), Point(909649.3986619939, 6882486.232659903), Point(909149.3986619939, 6882486.232659903)]})

gdf = gpd.GeoDataFrame(points, crs={'init': 'epsg:2154'}, geometry=points.geometry) # https://epsg.io/2154
gdf.geometry = gdf.geometry.to_crs({'init': 'epsg:3857'}) # https://epsg.io/3857

print(gdf)

结果:

    Id                                     geometry
0  001   POINT (652238.9846394559 6275490.88115861)
1  002  POINT (653026.8173626668 6276225.689838496)
2  003  POINT (652266.5262486706 6276253.346635176)