最后,我找到了一种使用以下代码可视化德国所有邮政编码区域的好方法:
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),并且我需要能够缩放以调查地图的某些部分(尚未应用颜色,例如出生率或每个邮政编码区域的任何内容)。这样做的正确方法是什么?