2

我有一个笔记本(github 链接),我在其中使用 geopandas 绘制带有不同国家颜色的地图。根据绘图的顺序,有时它不尊重我指定的 figsize() 。我在 Ubuntu 20.04 和 Firefox 中本地运行的 jupyter 以及在 Chromium 中运行的 Binder 和 Colab 中反复看到这种行为。

有人可以帮助我了解发生了什么吗?这是一个错误还是我控制 geopandas/matplotlib 错误?

import matplotlib.pyplot as plt
import geopandas as gpd

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
sixc = world[ world['continent'] != 'Antarctica' ]
asia = world[ world['continent'] == 'Asia' ]
noam = world[ world['continent'] == 'North America']
swed = world[ world['iso_a3'] == 'SWE' ] 

# This works, makes a 2x1 landscapey aspect
axes = sixc.plot(figsize=(8,4), color='lightgrey')
asia.plot(ax=axes, color='green')
noam.plot(ax=axes, color='purple')

# Plotting swed at the end breaks figsize, makes it squareish
axes = sixc.plot(figsize=(8,4), color='lightgrey')
asia.plot(ax=axes, color='green')
noam.plot(ax=axes, color='purple')
swed.plot(ax=axes, color='yellow')

# Plotting swed in the middle makes it ok again
axes = sixc.plot(figsize=(8,4), color='lightgrey')
asia.plot(ax=axes, color='green')
swed.plot(ax=axes, color='yellow')  
noam.plot(ax=axes, color='purple')

(另外,对于我的一些学生(但不是全部!),那个破碎的情节也没有浅灰色背景国家。这可能是导致方面问题的任何结果/副作用吗?)

4

2 回答 2

4

geopandas您使用的 v 0.8.1 中,您使用default的绘图命令的纵横比为auto. 因此,您获得的输出图将具有不可预测的纵横比值。尝试

print(axes.get_aspect()) 

对于每个情节。在以前版本的 geopandas 中,会得到一个equal输出,并且绘图是正确的。但在您的情况下,您将获得不代表相等方面的值。

为了简单地解决您的问题,您可以添加这行代码:

 axes.set_aspect('equal')

在最后一个情节声明之后。

于 2021-01-09T04:44:11.803 回答
1

这是关于纵横比如何工作的一个很好的答案matplotlib;然而,尽管这些概念很有帮助,但它GeoSeries.plot是建立在 matplotlib 之上的,因此就解决问题的代码而言,这个答案可能无济于事。

最简单的方法是传递aspect='equal'给每个链接GeoSeries.plot默认是'auto'链接中的“方面”)

至于您的其他问题,我无法重现,我会让这些学生升级到最新版本或使用不同的 IDE。这在 Jupyter Notebook 中对我有用。

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
sixc = world[ world['continent'] != 'Antarctica' ]
asia = world[ world['continent'] == 'Asia' ]
noam = world[ world['continent'] == 'North America']
swed = world[ world['iso_a3'] == 'SWE' ] 

# This works, makes a 2x1 landscapey aspect
axes = sixc.plot(figsize=(8,4), color='lightgrey', aspect='equal')
asia.plot(ax=axes, color='green', aspect='equal')
noam.plot(ax=axes, color='purple', aspect='equal')

# Plotting swed at the end breaks figsize, makes it squareish
axes = sixc.plot(figsize=(8,4), color='lightgrey', aspect='equal')
asia.plot(ax=axes, color='green', aspect='equal')
noam.plot(ax=axes, color='purple', aspect='equal')
swed.plot(ax=axes, color='yellow', aspect='equal')

# Plotting swed in the middle makes it ok again
axes = sixc.plot(figsize=(8,4), color='lightgrey', aspect='equal')
asia.plot(ax=axes, color='green', aspect='equal')
swed.plot(ax=axes, color='yellow', aspect='equal')
noam.plot(ax=axes, color='purple', aspect='equal')

在此处输入图像描述

此外,作为旁注,数据集与代码有一些不规则99ISO。我相信它会在更新中得到修复(参见https://github.com/geopandas/geopandas/issues/1041)....您可以使用以下代码手动修复它们:

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
world.loc[world['name'] == 'France', 'iso_a3'] = 'FRA'
world.loc[world['name'] == 'Norway', 'iso_a3'] = 'NOR'
world.loc[world['name'] == 'Somaliland', 'iso_a3'] = 'SOM'
world.loc[world['name'] == 'Kosovo', 'iso_a3'] = 'RKS'
于 2021-01-09T02:03:33.117 回答