最后,我找到了一种使用以下代码可视化德国所有邮政编码区域的好方法:
import json
import matplotlib.pyplot as plt
from descartes import PolygonPatch
jfile = open('postal_code_areas.geojson', 'r')
jdata = json.load(jfile)
jfile.close()
fig = plt.figure()
ax = fig.gca()
for poly in jdata['features']:
poly = poly['geometry']
ax.add_patch(PolygonPatch(poly, fc='#6699cc', ec='#000000', alpha=0.5, zorder=2 ))
ax.axis('scaled')
plt.show()
不幸的是,结果图像并不完全符合我的要求:
首先,我不明白为什么这看起来如此压缩?
我需要这种图像的可移植格式(如PNG),我需要能够缩放才能调查地图的某些部分(尚未应用任何颜色,例如出生率或无论每个邮政编码区域)。 这样做的正确方法是什么?