Modulo Python per eliminare le funzionalità SHP (senza Desktop GIS installato)

Ho un software (non Arc) che gira di notte su una macchina non di produzione che aggiorna le funzionalità dello shapefile da un database esterno. Periodicamente vorrei eliminare tutte le funzionalità del file di forma (non il file stesso, che deve rimanere) e lasciare che il software "ricostruisca" il file di forma da zero. Vorrei automatizzare questo processo.

Non ho alcun software GIS attualmente installato su quella macchina. Speravo di poter scrivere una routine in Python che eliminasse automaticamente le funzionalità, proprio come lo strumento di geoprocessing di Elimina funzionalità di Arc .

Ci sono dei moduli Python che mi permetterebbero di farlo? Preferibilmente open-source? Ho guardato Shapely e PyShp ma non ho visto nulla che mi permettesse di eliminare le funzionalità in massa o che corrispondessero a una clausola WHERE. Possono scrivere funzionalità e analizzarle, ma non hanno mai visto le funzioni ELIMINA FUNZIONI da nessuna parte.

Devo sicuramente mancare qualcosa di semplice ...

EDIT: ho 35 cartelle (diverse estensioni geografiche, tutte nella loro proiezione), con 35-65 shapefile con circa 1000 shapefile da gestire.

Puoi usare l' API python GDAL / OGR , il codice sarà così:

from osgeo import ogr

shapefile = ogr.Open( "shapfile.shp",1 )
layer=shapefile.GetLayerByIndex(0)
count=layer.GetFeatureCount()
for feature in range(count):
    layer.DeleteFeature(feature)

La riga di comando ogr2ogr con una clausola where garantita per creare risultati vuoti è un metodo rapido e semplice:

ogr2ogr output.shp input.shp -where "FID < 0"

La pagina panoramica di python e OGR (e GDAL) è http://trac.osgeo.org/gdal/wiki/GdalOgrInPython

Puoi farlo in pyshp. È semplice ma non ovvio perché non ho mai immaginato questo caso d'uso. Ma ha senso per le applicazioni di aggiornamento automatizzato. Ho testato le seguenti 6 righe di codice e ha funzionato alla grande:

import shapefile
r = shapefile.Reader("myshape")
w = shapefile.Writer(r.shapeType)
# This line will give us the same dbf schema
w.fields = r.fields
# Use the original bounding box as a place holder in the header
w.bbox = lambda: r.bbox
w.save("myshape")

Ora hai uno shapefile scritto sopra l'originale che ha le intestazioni corrette e i campi dbf originali. Si aprirà in modo sicuro nel software GIS e nelle librerie shapefile ma non ha funzionalità o record dbf.

La funzione lambda trasferisce il riquadro di delimitazione originale come segnaposto. Puoi mettere tutti i valori float che desideri in una matrice di [xmin, ymin, xmax, ymax]. Esempio:

w.bbox = lambda: [0.0, 0.0, 0.0, 0.0]

Anche cambiare i campi dbf è semplice e documentato nei documenti pyshp.

Spero possa aiutare.

Perché non salvare una copia dello Shapefile vuoto e sovrascrivere lo Shapefile di interesse.