Python 在 basemap 上绘制 shapefile

Question

我正在使用 python 3.6 绘制一个包含美国大型分水岭的 shapefile。我遇到的问题是该 shapefile 有许多组件，例如 .dbf 文件、.prj 文件...等等。我不确定是否必须单独读取所有文件然后进行绘图，或者是否有一个命令可以让我一次读取所有文件然后进行绘图。

这是迄今为止我的代码:

from mpl_toolkits.basemap import Basemap

import matplotlib.pyplot as plt

mydbf =('filename.dbf')
myprj =('filename.prj')
myqpj =('filename.qpj')
myshp =('filename.shp')
myshx =('filename.shx')


map= Basemap(projection='cyl',
            lon_0=180,
            lat_0=0,
            resolution='l')
map.drawmapboundary(fill_color='aqua')
map.fillcontinents(color='#ddaa66',lake_color='aqua')
map.drawcoastlines()

map.readshapefile(mydbf,myprj,myqpj,myshp,myshx,'Watersheds')
plt.show()

运行此代码时收到的错误消息是:OSError:无法找到文件名.dbf.shp

Answer 1

在 basemap tutorial在 readthedocs 上它清楚地说明了 readshapefile 函数

The first parameter shapefile name must go without the shp extension. The library assumes that all shp, sbf and shx files will exist with this given name

因此您的电话应该是

map.readshapefile("filename",'Watersheds')