使用 Geopandas 绘制缺失值

Question

nan我的 shapefile在某些列（例如 GDP）上有一些缺失值（用 表示）。当绘制而不处理这些缺失值时，图例显示如下：

这不是我想要的。因此，我用字符串“missing”替换缺失值，然后重做绘图。毫不奇怪，我收到了这样的错误消息TypeError: '<' not supported between instances of 'str' and 'float'。

我的问题是： 1. Geopandas 如何处理缺失值？它将缺失值存储在字符串或其他类型的数据中吗？2.如何保留这些缺失值并使用图例标签重新绘制以显示缺失值？

Answer 1

import numpy as np
import matplotlib.pyplot as plt
import geopandas as gpd
import pysal.viz.mapclassify as mc
from matplotlib.colors import rgb2hex
from matplotlib.colors import ListedColormap
plt.style.use('seaborn')

gdf = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
# generate random data
gdf['random'] = np.random.normal(100, 10, len(gdf))
# assign missing values
gdf.loc[np.random.choice(gdf.index, 40), 'random'] = np.nan

这里的基本思想是根据您想要用于数字数据的分类方法（例如分位数、百分位数等）创建类别/字符串列。之后，我们绘制该字符串列，以便我们可以传递自定义的颜色图（用灰色表示缺失值）。

# categorize the numerical column
k = 5
quantiles = mc.Quantiles(gdf.random.dropna(), k=k)
gdf['random_cat'] = quantiles.find_bin(gdf.random).astype('str')

gdf.loc[gdf.random.isnull(), 'random_cat'] = 'No Data'

# add grey to a colormap to represent missing value
cmap = plt.cm.get_cmap('Blues', k)
cmap_list = [rgb2hex(cmap(i)) for i in range(cmap.N)]
cmap_list.append('grey')
cmap_with_grey = ListedColormap(cmap_list)

# plot map
fig, ax = plt.subplots(figsize=(12, 10))
gdf.plot(column='random_cat', edgecolor='k', cmap=cmap_with_grey,
         legend=True, legend_kwds=dict(loc='center left'),
         ax=ax)

# get all upper bounds in the quantiles category
upper_bounds = quantiles.bins
# get and format all bounds
bounds = []
for index, upper_bound in enumerate(upper_bounds):
    if index == 0:
        lower_bound = gdf.random.min()
    else:
        lower_bound = upper_bounds[index-1]

    bound = f'{lower_bound:.2f} - {upper_bound:.2f}'
    bounds.append(bound)

# get all the legend labels
legend_labels = ax.get_legend().get_texts()
# replace the numerical legend labels
for bound, legend_label in zip(bounds, legend_labels):
    legend_label.set_text(bound)

您可能想看看以下帖子：

GeoPandas 中的格式/圆形数字图例标签

以十六进制格式提取 matplotlib 颜色图

python 中的 Matplotlib.colors.ListedColormap

更改主图图例标签文本

更新为 geopandas0.8.1：

missing_kwds您现在可以简单地在绘图函数中传递一个参数：

fig, ax = plt.subplots(figsize=(12, 10))

missing_kwds = dict(color='grey', label='No Data')

gdf.plot(column='random', scheme='Quantiles', k= 5,
         legend=True, legend_kwds=dict(loc='center left'),
         missing_kwds=missing_kwds, ax=ax)