用Basemap创建一些地图后,我变得很热情.我想集成shapefile信息,比如说一个多边形,但是有一些问题.我在这里下载了巴伐利亚村庄的寄宿生:
现在,我想为雷根斯堡集成一个多边形.我可以使用此代码获取信息,但是我有一些问题
#!/usr/bin/env python
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import shapefile
map = Basemap(projection='merc',
resolution='l',
area_thresh=0.01,
llcrnrlon=9.497681, llcrnrlat=47.827908,
urcrnrlon=12.683716, urcrnrlat=50.408517)
map.drawcountries(color="gray")
map.fillcontinents(color='#c8dfb0', lake_color='#53BEFD')
map.drawmapboundary(color='black', linewidth=0.5, fill_color='#53BEFD')
sf = shapefile.Reader("BY_Gemeinden/BY_Gemeinden_WM.shp")
shapes = sf.shapes()
records = sf.records()
for record, shape in zip(records