我正在尝试将 shapefile 绘制为 Folium 地图的点,但我正在拉扯头发以使投影正确。数据在 EPSG:28992 系统,RD 新坐标系中。这是一个投影坐标系,因此它以米为单位给出坐标。要在folium 地图上绘制它,我需要将其转换为 EPSG:3857,如果我是正确的,它适用于度数。不幸的是,在用函数转换它之后,我完全是胡说八道to_crs()
。
首先我加载数据:
import geopandas
WTG = geopandas.read_file('WTGs.shp')
print(WTG.head())
显示了这一点:
TEXTSTRING geometry
0 A27-09 POINT (151809.960 484552.520)
哪个是对的。它还具有正确的 CRS,如以下输出所示print(WTG.crs)
:
epsg:28992
现在要将其转换为其他 CRS,我使用以下简单代码:
WTG = WTG.to_crs(epsg=3857)
但是现在检查结果print(WTG.head())
对于几何列来说完全是一派胡言。
TEXTSTRING geometry
0 A27-09 POINT (594489.076 6863453.514)
几何点应该以度为单位。换句话说,坐标应该在 52 和 5 左右。print(WTG.crs)
虽然的输出是正确的:
epsg:3857
你们中有人知道如何解决这个问题吗?还是我错过了什么?提前致谢!