使用python

时间:2018-01-09 12:27:19

标签: python shapefile pyshp

我有一个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没有深入探讨这些事情。

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

2 个答案:

答案 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