1

我正在尝试通过以下方式处理底图中的 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可能不干净，但我不知道如何解决这个问题。任何想法？提前谢谢了。

4

3 回答 3

0

我刚刚用 Qgis 创建的 shapefile 出现了这个错误，并试图在 python 中加载。我将它跟踪到正在读取属性表的 mpl_toolskits.basemap 库中。我并没有真正进一步调试它，但我有一堆可以删除的属性（因为我不再使用它们了）。完成此操作后，我能够在 python 中成功加载 shapefile。

于 2019-02-14T19:15:28.730 回答
0

我有同样的错误。它是由错误的 shapefile 引起的。我尝试了 QGIS Goeprocessing 工具“检查有效性”和“多部分到单部分”。我以为我过去使用这些方法解决了它，但这次没有任何帮助。

我通过从头开始绘制一个 shapefile 并替换损坏的文件来解决它。有趣的是，geodjango 完美地读取了相同的 shapefile，只有底图导致了错误。

于 2019-01-10T17:57:53.240 回答
0

好的，不用担心，我只是用从官方网站下载的另一个 shp 文件解决了这个问题。

于 2017-06-23T20:54:05.177 回答