我正在尝试通过以下方式处理basemap中的shapefile:
shp = fiona.open('Wien_test2.shp')
bds = shp.bounds
shp.close()
extra = 0.01
ll = (float(bds[0]), float(bds[1]))
ur = (float(bds[2]), float(bds[3]))
coords = list(chain(ll, ur))
w, h = coords[2] - coords[0], coords[3] - coords[1]
m = Basemap(
projection='tmerc',
lon_0=16.,
lat_0=48.,
ellps = 'WGS84',
llcrnrlon=coords[0] - extra * w,
llcrnrlat=coords[1] - extra + 0.01 * h,
urcrnrlon=coords[2] + extra * w,
urcrnrlat=coords[3] + extra + 0.01 * h,
lat_ts=0,
resolution='i',
suppress_ticks=True)
m.readshapefile(
'Wien_test2',
'Wien_test2',
color='none',
zorder=2)
虽然代码可以很好地处理另一个shapefile,但是有了这个shapefile,我得到了一个错误:
^{pr2}$我认为问题是形状文件可能不干净,但我不知道如何解决。有什么想法吗?提前致谢。在
我也犯了同样的错误。它是由错误的形状文件引起的。我尝试了QGIS Goeprocessing工具“check validity”和“multiparts to singleparts”。我以为我以前用这些方法解决了这个问题,但这次没有任何帮助。在
我从零开始画了一个形状文件并替换了损坏的形状文件。有趣的是,同一个shapefile被geodjango完美地读取了,只有basemap导致了错误。在
好的,不用担心,我只是用另一个从官方网站下载的shp文件解决了这个问题。在
我刚刚用Qgis创建了一个shapefile,并试图在python中加载它。我查到了mpl_toolskits.basemap正在读取属性表的库。我没有进一步调试,但是我有一堆属性可以删除(因为我不再使用它们了)。在我完成这项工作之后,我能够成功地用python加载shapefile。在
相关问题 更多 >
编程相关推荐