如何以编程方式检查形状的数量=表记录的数量？

我有大约1000个shapefile已损坏（请参阅附件的错误消息）。shapefile是从eCognition Developer 8生成的。有一个脚本工具似乎可以在将shapefile识别为损坏后对其进行修复。

编辑：

我想创建一个快速脚本来遍历我所有的shapefile，并检查形状的数量是否与表记录匹配。我可以使用以下方法对表记录进行计数：

# Name: fcCount.py
# Purpose: calculate the number of features in a featureclass

# Import system modules
import arcpy
from arcpy import env

env.workspace = "C:/data"
Sample = "MyShp.shp"
result_dbf = int(arcpy.GetCount_management(Sample).getOutput(0)) 
print result_dbf

我最终想创建某种逻辑检查，例如：

if result_dbf = result_shp:
    pass
else:
    print "There is a problem with" + str(Sample)

如何在不访问.dbf文件的情况下直接计算形状？或者换句话说，以编程方式检查形状的数量是否与表记录的数量匹配的最佳方法是什么？

那使用pyshp呢？我用pip安装了它，下面尝试的几乎是自述文件：

>>> import shapefile
>>> sf = shapefile.Reader("/Users/chad/CoalOutcrops.shp")
>>> shapes = sf.shapes()
>>> len(shapes)
33732
>>> records = sf.records()
>>> len(records)
33732
>>>

不幸的是（或者幸运的是？）我没有任何可测试的shapefile来测试是否为否。形状可以！=否。记录。

请稍等，由于Kirk在下面的评论中的想法，我现在有了一个shapefile。我备份了dbf，制作了整个shapefile的副本，删除了一些功能，然后将备份的dbf重命名回原始文件，瞧瞧，形状的数量<记录的数量：

>>> sf = shapefile.Reader("/Users/chad/CoalOutcrops.shp")
>>> records = sf.records()
>>> len(records)
33732
>>> shapes = sf.shapes()
>>> len(shapes)
33721
>>>

从您的疑问中看来，您似乎真正想要做的只是确定shapefile是否存在问题（在这种情况下，记录不匹配）。如果您只需要确定有问题的对象，则实际上不需要统计DBF和Shapefile中的记录来确定是否有错误。原因如下：

如果您尝试在具有不同记录计数的shapefile上运行GetCount函数，它将失败并显示以下错误：

错误000229：无法打开。执行失败（GetCount）。

由于GetCount函数在这种情况下会失败，并且您要做的就是识别错误的shapefile，因此可以在代码中使用try / except子句来捕获此文件，而不是之前尝试使用的if / else。

我随意添加“列出FeatureClasses”代码并进行循环，以便您可以测试工作区中的所有FC，而无需手动测试每个FC。

# Import system modules
import arcpy
from arcpy import env

env.workspace = "C:/data"

fcList = arcpy.ListFeatureClasses()

for fc in fcList:
    try:
        result_dbf = int(arcpy.GetCount_management(fc).getOutput(0))
        print fc + ": " + str(result_dbf) + " records"
    except:
        print "There is a problem with: " + str(fc)

shapefile格式已记录。我猜想shp文件中的记录数与dbf文件中的记录数不对应。

shp文件格式在此处记录。因此，您可以编写一个程序来计算形状的数量。dbf格式在许多地方都有记录，您应该能够找到用于计数行的样本，例如here。

附加的脚本循环遍历目录，并检查形状的数量是否与每个shapefile的记录数量相匹配。

import arcpy, os, shapefile
from arcpy import env

env.workspace = r"C:\path\to\shapefiles"
Dir = env.workspace

fclist = arcpy.ListFeatureClasses()

for fc in fclist:

    myfc = os.path.join(Dir, fc)
    sf = shapefile.Reader(str(myfc))
    shapes = sf.shapes()
    shape_total = len(shapes)
    records = sf.records()
    record_total = len(records)

    if shape_total != record_total:
        print "There is a problem with " + str(fc)
    else:
        print str(fc) + " passed"

使用检查几何图形应该可以帮助您完成第一步。
Onus
Repair Geometry将允许您选择要修复的问题的顺序和优先级。
这是其他一些旧版本的链接。当您运行shapefile checker时，您是否需要重新构建dbf？
这是创建要匹配的记录的步骤。发生了两件事之一导致该错误。

1. shp具有一个已被另一个软件/进程删除/删除的对象（空间）。
2. dbf包含引用空几何的记录。
   有几件事会导致此。
   shx实际上是两者之间的索引。
   不计算dbf记录而计算形状只是解决方案的一半。

查看关于shapefile的维基百科文章，.shx文件应在.shp文件上包含索引，而不是.dbf文件上的索引。因此，可能有必要检查.shx和.shp是否合适。

可以打开不带.dbf的shapefile（这意味着您没有属性表），但是索引损坏将生成错误消息。