python - 在Geopandas中绘图时管理投影

我正在使用geopandas绘制意大利地图。

import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize = (20,30))

region_map.plot(ax=ax, color='white', edgecolor='black')
plt.xlim([6,19])
plt.ylim([36,47.7])
plt.tight_layout()
plt.show()

这是在将region_map正确定义为“几何” GeoSeries之后的结果。



但是，即使更改figsize中的plt.subplots，我也无法修改图形长宽比。我是否错过了一些琐碎的事情，或者这可能是Geopandas问题？

谢谢

最佳答案

您的源数据集（region_map）显然是在地理坐标系中“编码”的（单位：经度和纬度）。可以安全地假设这是WGS84（EPSG：4326）。如果您想让自己的绘图看起来更像在Google Maps中，则必须将其坐标重新投影到许多投影坐标系之一（单位：米）中。您可以使用全球公认的WEB MERCATOR（EPSG：3857）。

Geopandas使这个过程变得尽可能容易。您只需要了解我们如何处理计算机科学中的坐标投影以及通过其EPSG代码学习最流行的CRS的基础知识。

import matplotlib.pyplot as plt

#If your source does not have a crs assigned to it, do it like this:
region_map.crs = {"init": "epsg:4326"}

#Now that Geopandas what is the "encoding" of your coordinates, you can perform any coordinate reprojection
region_map = region_map.to_crs(epsg=3857)

fig, ax = plt.subplots(figsize = (20,30))
region_map.plot(ax=ax, color='white', edgecolor='black')

#Keep in mind that these limits are not longer referring to the source data!
# plt.xlim([6,19])
# plt.ylim([36,47.7])
plt.tight_layout()
plt.show()

我强烈建议阅读official GeoPandas docs有关管理预测的信息。

关于python - 在Geopandas中绘图时管理投影，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/54445754/