Question

我有一个ESRI shapefile .shp（包含所有相关文件，如.shx，.dbf等等）我想编辑 - 我需要删除第一条记录并保存文件。

为此，我安装了pyshp并尝试解析和编辑shapefile。这就是我的尝试：

import shapefile
e = shapefile.Editor('location/of/my/shp')
e.shapes()
# example output
>>> [<shapefile._Shape instance at 0x7fc5e18d93f8>,
     <shapefile._Shape instance at 0x7fc5e18d9440>,
     <shapefile._Shape instance at 0x7fc5e18d9488>,
     <shapefile._Shape instance at 0x7fc5e18d94d0>,
     <shapefile._Shape instance at 0x7fc5e18d9518>]

从这里我想删除第一个条目<shapefile._Shape instance at 0x7fc5e18d93f8>，然后保存文件：

e.delete(0) # I tried e.delete(shape=0) too
e.save()

但是，新保存的文件中仍然可以使用该记录。

不幸的是，documentation没有深入探讨这些事情。

我如何实现目标？如何在保存文件之前检查删除是否成功？

Answer 1

完全遵循您描述的程序似乎对我来说效果很好。我首先打开一个shapefile：

>>> e = shapefile.Editor('example')

此文件有三种形状：

>>> e.shapes()
[<shapefile._Shape instance at 0x7f6cb5f67dd0>, <shapefile._Shape instance at 0x7f6cb5f67f38>, <shapefile._Shape instance at 0x7f6cb5f6e050>]

我删除第一个形状并保存文件：

>>> e.delete(0)
>>> e.save('example')

现在我重新打开文件：

>>> e = shapefile.Editor('example')

我可以看到它现在只有两种形状：

>>> e.shapes()
[<shapefile._Shape instance at 0x7f6cb5f6e518>, <shapefile._Shape instance at 0x7f6cb5f6e560>]

Answer 2

我对pyshp并不熟悉，但可以使用ogr轻松解决这个问题，它允许使用矢量数据并成为gdal库的一部分。

from osgeo import ogr

fn = r"file.shp"  # The path of your shapefile
ds = ogr.Open(fn, True)  # True allows to edit the shapefile
lyr = ds.GetLayer()

print("Features before: {}".format(lyr.GetFeatureCount()))
lyr.DeleteFeature(0)  # Deletes the first feature in the Layer

# Repack and recompute extent
# This is not mandatory but it organizes the FID's (so they start at 0 again and not 1)
# and recalculates the spatial extent.
ds.ExecuteSQL('REPACK ' + lyr.GetName())
ds.ExecuteSQL('RECOMPUTE EXTENT ON ' + lyr.GetName())

print("Features after: {}".format(lyr.GetFeatureCount()))

del ds

使用python

2 个答案: