¿Cómo comprobar mediante programación si el número de formas es igual al número de registros de la tabla?

Question

Tengo un puñado de aproximadamente 1000 shapefiles que están corruptos (ver mensaje de error adjunto). Los archivos shape se generaron con eCognition Developer 8. Hay un herramienta de script que parece reparar el shapefile una vez identificado como corrupto.

Edita:

Quiero crear un script rápido para recorrer todos mis shapefiles y comprobar si el número de shapes coincide con los registros de la tabla. Puedo contar los registros de la tabla utilizando el siguiente:

# 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

En última instancia, me gustaría crear algún tipo de comprobación lógica como:

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

¿Cómo puedo contar las formas directamente sin acceder al archivo .dbf? O, en otras palabras, ¿cuál es la mejor manera de comprobar mediante programación si el número de formas coincide con el número de registros de la tabla?

Answer 1

¿Qué le parece utilizar pyshp ? Lo instalé con pip y lo que he probado a continuación es más o menos directamente de la LÉAME :

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

Por desgracia (¿o por suerte?), no dispongo de ningún archivo shape para comprobar si el número de formas puede ser igual al número de registros.

Espera un momento, ahora sí tengo un shapefile de la hostia gracias a la idea de Kirk en los comentarios de abajo. Hice una copia de seguridad de la dbf, hizo una copia de todo el shapefile, eliminado algunas características, a continuación, cambiar el nombre de la copia de seguridad dbf de nuevo a la original, y he aquí, el número de formas < número de registros:

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

Answer 2

Por el sonido de su pregunta, parece que todo lo que realmente quiere hacer es determinar si un shapefile tiene o no problemas (en este caso, registros no coincidentes). Si sólo necesita identificar los que tienen problemas, en realidad no es necesario contar los registros en el DBF y Shapefile para determinar si está en error. He aquí por qué:

Si intenta ejecutar la función GetCount en un shapefile que tiene diferentes recuentos de registros, fallará con el error:

ERROR 000229 : No se puede abrir . Error al ejecutar (GetCount).

Dado que la función GetCount falla en este escenario, y todo lo que quieres hacer es identificar los shapefiles en error, puedes atrapar esto con una cláusula try/except en tu código, en lugar del if/else que estabas intentando usar previamente.

Me he tomado la libertad de añadir el código y el bucle "List FeatureClasses" para que pueda probar todas las FC de su espacio de trabajo sin tener que probar manualmente cada una de ellas.

# 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)

Answer 3

El formato shapefile está documentado. Supongo que el número de registros del archivo shp no se corresponde con el número de registros del archivo dbf.

El formato de archivo shp es documentado aquí . Así que usted podría escribir un programa para contar el número de formas. El formato dbf está documentado en muchos sitios y deberías poder encontrar ejemplos para contar filas, por ejemplo aquí .

Answer 4

El script adjunto recorre un directorio y comprueba si el número de formas coincide con el número de registros de cada 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"

Answer 5

Observando el artículo de la wikipedia sobre shapefiles , el archivo .shx debe contener un índice en el archivo .shp, no en el archivo .dbf. Así que puede ser necesario comprobar si .shx y .shp encajan.

Es posible abrir un shapefile sin .dbf (lo que significa que no tiene tabla de atributos), pero un índice roto generará un mensaje de error.