我正在尝试将形状文件绘制为指向 Folium 地图,但正在拔头发试图使投影正确。数据位于 EPSG:28992 系统(RD 新坐标系(中。这是一个投影坐标系,因此它以米为单位给出坐标。要在叶子地图上绘制它,我需要将其转换为 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
你们中有谁知道如何解决这个问题吗?还是我错过了什么?提前感谢!
EPSG 3857 不是以度为单位,而是以米为单位,因此没有理由假设转换是错误的。如果你想要学位,你需要4326。
该点位于乌得勒支北部,对吗?如果是这样,一切都按预期工作。