如何手动定义Lisa群集的颜色?

时间:2019-05-26 00:40:47

标签: python maps cluster-analysis geopandas pysal

我正在尝试绘制LISA CLUSTERS的一些地图。 O更改了lisa_cluster函数的code以指定我想要的颜色。我使用了通用的5种颜色列表,并手动对其进行了更改

from matplotlib import patches, colors
import palettable

palettable.colorbrewer.sequential.Greys_5_r.colors = [[60,60,60],[105,105,105],[0,0,255],[255,255,0],[240,240,240]]
paleta = palettable.colorbrewer.sequential.Greys_5_r.mpl_colormap

def lisa_cluster(moran_loc, gdf, p=0.05, ax=None,
                 legend=True, legend_kwds=None, **kwargs):
...
 if ax is None:
        figsize = kwargs.pop('figsize', None)
        fig, ax = plt.subplots(1, figsize=figsize)
    else:
        fig = ax.get_figure()

    gdf.assign(cl=labels).plot(column='cl', categorical=True,
                               k=2, cmap=paleta, linewidth=0.1, ax=ax,
                               edgecolor='white', legend=legend,
                               legend_kwds=legend_kwds, **kwargs)
    ax.set_axis_off()
    ax.set_aspect('equal')
    return fig, ax

所以我希望每个象限中的区域具有以下颜色:

1(HH)-黑色
2(HL)-深灰色
3(LL)-黄色 4(LH)-蓝色
不重要-浅灰色

问题在于颜色正在融合,我不知道为什么。 我用相应象限标记了区域,以显示enter image description here

2003年和2004年还可以。在2002年的地图上,黄色和蓝色(我认为是蓝色和浅灰色)融合了

1 个答案:

答案 0 :(得分:2)

解决方案是更改从功能mask_local_auto导入的对象。 代码是here

颜色在对象colors5中定义:

(...)
#create a mask for local spatial autocorrelation  
cluster = moran_hot_cold_spots(moran_loc, p)

cluster_labels = ['ns', 'HH', 'LH', 'LL', 'HL']
labels = [cluster_labels[i] for i in cluster]

colors5 = {0: 'lightgrey',
           1: '#d7191c',
           2: '#abd9e9',
           3: '#2c7bb6',
           4: '#fdae61'}
colors = [colors5[i] for i in cluster]  # for Bokeh
# for MPL, keeps colors even if clusters are missing:
x = np.array(labels)
y = np.unique(x)
colors5_mpl = {'HH': '#d7191c',
               'LH': '#abd9e9',
               'LL': '#2c7bb6',
               'HL': '#fdae61',
               'ns': 'lightgrey'}
colors5 = [colors5_mpl[i] for i in y]  # for mpl
(...)