如何讓使用python繪制中國地圖并給特定地區(qū)上色？

在使用python Basemap繪制地圖并給特定地區(qū)染色的過程中，有些地區(qū)無法上色。請問各位大神有無解決方法？
圖片描述

我是準備給14個地區(qū)染色，但是能成功上色的只有6個省和直轄市。

請問這是Basemap畫地圖的bug嗎？

在繪制的過程中參考的是
鏈接描述

大神的教程。

我是準備給14個地區(qū)染色，但是能成功上色的只有6個省和直轄市。

請問這是Basemap畫地圖的bug嗎？



代碼如下：

plt.figure(figsize=(16,10))

map = Basemap(llcrnrlon=77, llcrnrlat=14, urcrnrlon=140, urcrnrlat=51,projection='lcc',lat_1=33, lat_2=45, lon_0=100)

map.drawcoastlines()

map.readshapefile('CHN_adm1', 'states', drawbounds=True)

#map.drawcountries()





statenames=[]

for shapedict in map.states_info:

    statename = shapedict['NL_NAME_1']

    p = statename.split('|')

    if len(p) > 1:

        s = p[1]

    else:

        s = p[0]

    statenames.append(s)



ax = plt.gca()



seg = map.states[statenames.index('四川')]

seg1 = map.states[statenames.index('湖南')]

seg2 = map.states[statenames.index('江蘇')]

seg3 = map.states[statenames.index('上海')]

seg4 = map.states[statenames.index('安徽')]

seg5 = map.states[statenames.index('浙江')]

seg6 = map.states[statenames.index('北京')]

seg7 = map.states[statenames.index('山東')]

seg8 = map.states[statenames.index('江西')]

seg9 = map.states[statenames.index('廣西')]

seg10 = map.states[statenames.index('湖北')]





poly = Polygon(seg,facecolor = 'crimson',edgecolor = 'r')

poly1 = Polygon(seg1,facecolor = 'crimson',edgecolor = 'r')

poly2 = Polygon(seg2,facecolor = 'crimson',edgecolor = 'r')

poly3 = Polygon(seg3,facecolor = 'violet',edgecolor = 'r')

poly4 = Polygon(seg4,facecolor = 'violet',edgecolor = 'r')

poly5 = Polygon(seg5,facecolor = 'violet',edgecolor = 'r')

poly6 = Polygon(seg6,facecolor = 'gold',edgecolor = 'r')

poly7 = Polygon(seg7,facecolor = 'gold',edgecolor = 'r')

poly8 = Polygon(seg8,facecolor = 'coral',edgecolor = 'r')

poly9 = Polygon(seg9,facecolor = 'coral',edgecolor = 'r')

poly10 = Polygon(seg10,facecolor = 'coral',edgecolor = 'r')





ax.add_patch(poly)

ax.add_patch(poly1)

ax.add_patch(poly2)

ax.add_patch(poly3)

ax.add_patch(poly4)

ax.add_patch(poly5)

ax.add_patch(poly6)

ax.add_patch(poly7)

ax.add_patch(poly8)

ax.add_patch(poly9)

ax.add_patch(poly10)



plt.show()

回答

編輯回答

夏木

這里有一篇現(xiàn)成的文章 -- https://segmentfault.com/a/11...

2017年7月5日 21:17

編輯回答

心沉

import matplotlib.pyplot as plt
from matplotlib.patches import Polygon
from mpl_toolkits.basemap import Basemap

plt.figure(figsize=(16,10))
map = Basemap(llcrnrlon=77, llcrnrlat=14, urcrnrlon=140, urcrnrlat=51,projection='lcc',lat_1=33, lat_2=45, lon_0=100)
map.drawcoastlines()
map.readshapefile('map/CHN_adm1', 'states', drawbounds=True)
ax = plt.gca()

for i, shapedict in enumerate(map.states_info):
    p = shapedict['NL_NAME_1'].split('|')
    s = p[1] if len(p) > 1 else p[0]
    if s == '山東':
        ax.add_patch(Polygon(map.states[i], facecolor = 'gold', edgecolor = 'r'))
    elif s == '浙江':
        ax.add_patch(Polygon(map.states[i], facecolor = 'violet', edgecolor = 'r'))

只畫山東和浙江的例子

圖片描述

2018年8月16日 06:27