Python-Modul zum Löschen von SHP-Funktionen (ohne installiertes Desktop-GIS)

Ich habe eine Software (nicht Arc), die jede Nacht auf einem Nicht-Produktions-Computer ausgeführt wird, der Shapefile-Features aus einer externen Datenbank aktualisiert. In regelmäßigen Abständen möchte ich alle Features im Shapefile löschen (nicht die Datei selbst, die verbleiben muss) und die Software das Shapefile von Grund auf neu erstellen lassen. Ich möchte diesen Prozess automatisieren.

Ich habe derzeit keine GIS-Software auf diesem Computer installiert. Ich hatte gehofft, in Python eine Routine erstellen zu können , mit der die Features automatisch gelöscht werden, ähnlich wie mit dem Geoverarbeitungswerkzeug "Features löschen" von Arc .

Gibt es Python-Module, mit denen ich das machen könnte? Am liebsten Open Source? Ich habe mir Shapely und PyShp angesehen , aber nichts gefunden, das es mir erlaubt, die Features in der Masse zu löschen, oder das mit einer WHERE-Klausel übereinstimmt. Sie können Features schreiben und analysieren, haben jedoch nirgendwo DELETE FEATURES-Funktionen gesehen.

Mir muss auf jeden Fall etwas Einfaches fehlen ...

BEARBEITEN: Ich habe 35 Ordner (verschiedene geografische Ausmaße, alle in ihrer eigenen Projektion) mit 35-65 Shapefiles, mit denen ungefähr 1000 Shapefiles bearbeitet werden können.

Sie können die GDAL / OGR-Python-API verwenden . Der Code sieht folgendermaßen aus :

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)

Die Befehlszeile ogr2ogr mit einer where-Klausel, die garantiert leere Ergebnisse erzeugt, ist eine schnelle und einfache Methode:

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

Die Übersichtsseite für Python und OGR (und GDAL) lautet http://trac.osgeo.org/gdal/wiki/GdalOgrInPython

Sie können dies in pyshp tun. Es ist einfach, aber nicht offensichtlich, da ich mir diesen Anwendungsfall nie vorgestellt habe. Bei automatisierten Update-Anwendungen ist dies jedoch sinnvoll. Ich habe die folgenden 6 Codezeilen getestet und es hat super funktioniert:

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

Sie haben jetzt ein Shapefile über das Original geschrieben, das die richtigen Header und die ursprünglichen DBF-Felder enthält. Es wird sicher in GIS-Software- und Shapefile-Bibliotheken geöffnet, hat jedoch keine Features oder DBF-Datensätze.

Die Lambda-Funktion überträgt den ursprünglichen Begrenzungsrahmen als Platzhalter. Sie können beliebige Gleitkommawerte in ein Array von [xmin, ymin, xmax, ymax] einfügen. Beispiel:

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

Das Ändern von DBF-Feldern ist ebenfalls einfach und in den pyshp-Dokumenten dokumentiert.

Ich hoffe, das hilft.

Speichern Sie eine Kopie des leeren Shapefiles und überschreiben Sie das gewünschte Shapefile.