转换Geopandas中形状文件的坐标

Question

我是空间分析的新手，但我在任何地方都找不到这个答案。在

我有一份用CRS坐标、经度和纬度表示的邮政编码列表，用OSN坐标表示的伦敦市行政区形状文件，我想把它们映射在一起，但事情就是这样。这是邮政编码的头

london_post_codes.head()
Out[81]:
postcode    latitude    longitude
0   WD6 1GS 51.658021   -0.255663
1   WD17 1LA    51.660366   -0.397525
2   WC2N 6LE    51.509413   -0.121676
3   WC2N 6NA    51.508363   -0.124454
4   WC2N 6ND    51.508216   -0.123829

而这是geopandas中读取的形状文件

^{pr2}$

我们可以看到多边形的坐标和邮政编码的坐标有什么不同。当我把它们画在一起的时候

fig, ax = plt.subplots()
borough.plot(ax = ax)
borough = gpd.read_file('statistical-gis-boundaries-london/ESRI/London_Borough_Excluding_MHW.shp')
london_post_codes.plot(kind='scatter',s=10, x='longitude', y='latitude',ax=ax)

有什么建议吗？在

Answer 1

解决办法很简单，就是改变坐标系。这些代码被称为epsg。横向/纵向CRS的代码为epsg=4326。因此

borough = gpd.read_file('London_Borough_Excluding_MHW.shp')
borough = borough.to_crs(epsg=4326)

然后剩下的。在