Come verificare a livello di codice se il numero di forme = numero di record della tabella?

Ho una manciata di circa 1000 shapefile che sono danneggiati (vedi messaggio di errore allegato). Gli shapefile sono stati generati da eCognition Developer 8. Esiste uno strumento di script che sembra riparare lo shapefile una volta identificato come corrotto.

Modificare:

Voglio creare uno script rapido per scorrere tutti i miei shapefile e verificare se il numero di forme corrisponde ai record della tabella. Posso contare i record della tabella usando quanto segue:

# 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

Alla fine vorrei creare una sorta di controllo logico come:

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

Come posso contare le forme direttamente senza accedere al file .dbf? O, in altre parole, qual è il modo migliore per verificare a livello di codice se il numero di forme corrisponde al numero di record della tabella?

Che dire dell'uso di pyshp ? L'ho installato con pip e quello che ho provato di seguito è praticamente uscito da README :

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

Sfortunatamente (o forse per fortuna?) Non ho alcun shapefile sollevato da testare per vedere se no. di forme può! = no. dei record.

Aspetta solo un minuto, ora ho uno shapefile sollevato grazie all'idea di Kirk nei commenti qui sotto. Ho eseguito il backup del dbf, ho creato una copia dell'intero file di forma, ho eliminato alcune funzionalità, quindi ho rinominato il dbf di cui è stato eseguito il backup con l'originale, ed ecco il numero di forme <numero di record:

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

Dal suono della tua domanda, sembra che tutto ciò che vuoi davvero fare sia determinare se un file di forma ha o meno problemi con esso (in questo caso, record non corrispondenti). Se hai solo bisogno di identificare quelli con problemi, in realtà non devi contare i record nel DBF e Shapefile per determinare se è un errore. Ecco perché:

Se si tenta di eseguire la funzione GetCount su uno shapefile con conteggi di record diversi, non riuscirà con l'errore:

ERRORE 000229 : impossibile aprire. Impossibile eseguire (GetCount).

Poiché la funzione GetCount non riesce in questo scenario e tutto ciò che vuoi fare è identificare gli shapefile in errore, puoi prenderli con una clausola try / tranne nel tuo codice, invece che se / else stavi tentando in precedenza di utilizzare.

Mi sono preso la libertà di aggiungere il codice e il loop "Elenco FeatureClasses" in modo da poter testare tutti gli FC nell'area di lavoro senza doverli testare manualmente.

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

Il formato dello shapefile è documentato. Immagino che il numero di record nel file shp non corrisponda al numero di record nel file dbf.

Il formato del file shp è documentato qui . Quindi potresti scrivere un programma per contare il numero di forme. Il formato dbf è documentato in molti punti e dovresti essere in grado di trovare campioni per contare le righe, ad esempio qui .

Lo script allegato scorre in una directory e controlla se il numero di forme corrisponde al numero di record per ciascun file di forma.

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"

L'uso della geometria di controllo dovrebbe farti passare per il primo passo.
Onus
Repair Geometry ti permetterà di selezionare l'ordine e la priorità del problema per il quale vuoi riparare.
ecco alcuni altri collegamenti alla versione precedente . Quando esegui il controllo dello shapefile, finisci con ricostruire dbf?
Questo è il passaggio che crea i record da abbinare. Una delle due cose si è verificata per causare l'errore.

1. Lo shp ha un oggetto (spaziale) che è stato eliminato / rilasciato da un altro software / processo.
2. Il dbf ha un record che fa riferimento a geometria nulla.
   Diverse cose possono causare questo.
   Lo shx è in realtà l'indice tra i due.
   Contare le forme senza contare i record dbf è solo metà della soluzione.

Guardando l' articolo di Wikipedia sui file shape , il file .shx dovrebbe contenere un indice sul file .shp, non sul file .dbf. Quindi potrebbe essere necessario verificare se .shx e .shp si adattano insieme.

È possibile aprire un file di forma senza .dbf (il che significa che non hai una tabella di attributi), ma un indice rotto genererà un messaggio di errore.