使用matplotlib和geopandas在具有不同投影的地图上注释标签时出现问题

Question

我使用的默认贴图如下：

world = geopandas.read_file(gpd.datasets.get_path('naturalearth_lowres'))

我可以从另一个GeoDataFrame（这里称为sample_gdf）成功地将标签注释到此地图，下面是循环示例：

for idx, row in sample_gdf.iterrows():
     plt.annotate(text=row['country_name'], # e.g. this column contains the names of each countries
                  xy=(row['longitude'], row['latitude']), # e.g. these columns are showing the coordinates of middle points of each countries
                  horizontalalignment='center')

这就是epsg=4326时的样子当我想改变地图的投影时，问题就开始了。上面变量world的默认CRS是epsg： 4326。一旦我像这样改变投影：

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
world = world.to_crs(epsg=3035)

专注于欧洲，我的标签不再出现在正确的位置。一周来，我一直在寻找解决这个问题的建议，但目前找不到任何解决方案。谢谢你的帮助。这是epsg=3035时的情况，标签显示在左下角。

Answer 1

所有这些都归结为您用来定位注释的['longitude']、['latitude']列中存储的内容。您需要确保这些都在正确的投影中，因为坐标需要反映绘图上使用的实际单位。

重新投影后，获取新坐标并使用这些坐标。

sample_gdf["x"] = sample_gdf.centroid.x
sample_gdf["y"] = sample_gdf.centroid.y

sample_gdf.plot()

for idx, row in sample_gdf.iterrows():
     plt.annotate(text=row['country_name'], # e.g. this column contains the names of each countries
                  xy=(row['x'], row['y']), # e.g. these columns are showing the coordinates of middle points of each countries
                  horizontalalignment='center')