我有一个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没有深入探讨这些事情。
我如何实现目标?如何在保存文件之前检查删除是否成功?
答案 0 :(得分: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>]
答案 1 :(得分:0)
我对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