我正在尝试通过以下方式处理底图中的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,我得到了错误:
Traceback (most recent call last):
File "maps.py", line 104, in <module>
zorder=2)
File "/usr/lib/python2.7/dist-packages/mpl_toolkits/basemap/__init__.py", line 2146, in readshapefile
for shprec in shf.shapeRecords():
File "/usr/lib/python2.7/dist-packages/mpl_toolkits/basemap/shapefile.py", line 543, in shapeRecords
for rec in zip(self.shapes(), self.records())]
File "/usr/lib/python2.7/dist-packages/mpl_toolkits/basemap/shapefile.py", line 515, in records
r = self.__record()
File "/usr/lib/python2.7/dist-packages/mpl_toolkits/basemap/shapefile.py", line 480, in __record
value = int(value)
ValueError: invalid literal for int() with base 10: '***********'我认为问题在于shapefile可能并不干净,但我不知道如何解决这个问题。有什么想法吗?在此之前,非常感谢您。
发布于 2017-06-24 04:54:05
好的,不用担心,我只是用从官方网站下载的另一个shp文件解决了这个问题。
发布于 2019-01-11 01:57:53
我也犯了同样的错误。它是由错误的shapefile引起的。我尝试了QGIS Goeprocessing工具“检查有效性”和“多部分到单个部分”。我以为我以前用这些方法解决了这个问题,但这次什么都没有用。
我通过从头开始绘制一个shapefile并替换损坏的shapefile解决了这个问题。有趣的是,同样的shapefile被geodjango完美地读取了,只有底图导致了错误。
发布于 2019-02-15 03:15:29
我用Qgis创建了一个shapefile,并试图在python中加载,但我遇到了这个错误。我将其跟踪到正在读取属性表的mpl_toolskits.basemap库中。我真的没有进一步调试,但我有一堆属性可以删除(因为我不再使用它们)。这样做之后,我就可以成功地在python中加载shapefile了。
https://stackoverflow.com/questions/44727931
复制相似问题