在 Geopandas 中绘图时管理投影

Managing projections when plotting in 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 正确定义为 'geometry' GeoSeries 的一部分之后的结果。

但是,我无法修改图形纵横比,即使在 plt.subplots 中改变 figsize。我是否遗漏了一些微不足道的东西,或者它可能是 geopandas 的问题?

谢谢

您的源数据集 (region_map) 在 地理坐标系 中显然是 "encoded"(单位:lats 和 lons)。可以安全地假设您的情况是 WGS84(EPSG:4326)。如果你想让你的地块看起来更像它在 Google 地图中的样子,你将不得不 重新投影 它的坐标到许多 投影坐标系之一(单位:米)。您可以使用全球通用的 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 有关管理预测的内容。