我正试图将一个形状文件绘制成一个叶状图,但为了得到正确的投影,我把头发拔了出来。数据位于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
现在,要将其转换为其他CR,我使用以下简单代码:
WTG = WTG.to_crs(epsg=3857)
但是现在用print(WTG.head())
检查结果会使geometry列完全失去意义
TEXTSTRING geometry
0 A27-09 POINT (594489.076 6863453.514)
几何图形点应以度为单位。换句话说,坐标应该在52和5之间。
print(WTG.crs)
的输出是正确的,尽管:
epsg:3857
你们有谁知道怎么解决这个问题吗?还是我遗漏了什么?提前谢谢
EPSG 3857不是以度数为单位,而是以米为单位,因此没有理由假设转换是错误的。如果你想要学位,你需要4326
该点位于乌得勒支北部,对吗?如果是这样的话,一切都按预期进行
相关问题 更多 >
编程相关推荐