python - shapefile和matplotlib : plot polygon collection of shapefile coordinates

我正在尝试使用python中的matplotlib在世界 map 上绘制国家的填充多边形。

我有一个shapefile，其中包含每个国家/地区的国家/地区边界坐标。现在，我想将这些坐标(针对每个国家)转换为带有matplotlib的多边形。不使用 basemap 。不幸的是，这些部分交叉或重叠。有没有动工，也许是使用点到点的距离。或者重新排序？

最佳答案

哈!
我发现了，如何..我完全忽略了sf.shapes [i] .parts信息!然后归结为:

#   -- import --
import shapefile
import matplotlib.pyplot as plt
import matplotlib.patches as patches
from matplotlib.patches import Polygon
from matplotlib.collections import PatchCollection
#   -- input --
sf = shapefile.Reader("./shapefiles/world_countries_boundary_file_world_2002")
recs    = sf.records()
shapes  = sf.shapes()
Nshp    = len(shapes)
cns     = []
for nshp in xrange(Nshp):
    cns.append(recs[nshp][1])
cns = array(cns)
cm    = get_cmap('Dark2')
cccol = cm(1.*arange(Nshp)/Nshp)
#   -- plot --
fig     = plt.figure()
ax      = fig.add_subplot(111)
for nshp in xrange(Nshp):
    ptchs   = []
    pts     = array(shapes[nshp].points)
    prt     = shapes[nshp].parts
    par     = list(prt) + [pts.shape[0]]
    for pij in xrange(len(prt)):
     ptchs.append(Polygon(pts[par[pij]:par[pij+1]]))
    ax.add_collection(PatchCollection(ptchs,facecolor=cccol[nshp,:],edgecolor='k', linewidths=.1))
ax.set_xlim(-180,+180)
ax.set_ylim(-90,90)
fig.savefig('test.png')

然后它将如下所示: