Mô-đun Python để xóa các tính năng SHP (không cài đặt Desktop Desktop)

Tôi có một phần mềm (không phải Arc) chạy hàng đêm trên một máy không sản xuất, cập nhật các tính năng shapefile từ cơ sở dữ liệu bên ngoài. Định kỳ tôi muốn xóa tất cả các tính năng trong shapefile (không phải chính tệp, phải giữ lại) và để phần mềm "xây dựng lại" shapefile từ đầu. Tôi muốn tự động hóa quá trình này.

Tôi hiện chưa có phần mềm GIS nào được cài đặt trên máy đó. Tôi đã hy vọng rằng tôi có thể tạo kịch bản một thói quen trong Python để tự động xóa các tính năng, giống như công cụ xử lý địa lý Xóa tính năng Xóa của Arc .

Có mô-đun Python nào cho phép tôi làm điều này không? Ưu tiên nguồn mở? Tôi đã xem Shapely và PyShp nhưng không thấy bất cứ điều gì cho phép tôi xóa các tính năng hàng loạt hoặc phù hợp với mệnh đề WHERE. Họ có thể viết các tính năng và phân tích chúng, nhưng chưa thấy các hàm TÍNH NĂNG XÓA ở bất cứ đâu.

Tôi chắc chắn phải thiếu một cái gì đó đơn giản ...

EDIT: Tôi có 35 thư mục (các phạm vi địa lý khác nhau, tất cả đều nằm trong hình chiếu riêng của chúng), với 35-65 shapefile với khoảng 1000 shapefile để xử lý.

Bạn có thể sử dụng API python GDAL / OGR , mã sẽ như thế:

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)

Dòng lệnh ogr2ogr với mệnh đề where được đảm bảo để tạo kết quả trống là một phương pháp nhanh chóng và dễ dàng:

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

Trang tổng quan về python và OGR (và GDAL) là http://trac.osgeo.org/gdal/wiki/GdalOgrInPython

Bạn có thể làm điều này trong pyshp. Nó đơn giản nhưng không rõ ràng vì tôi chưa bao giờ hình dung trường hợp sử dụng này. Nhưng nó có ý nghĩa cho các ứng dụng cập nhật tự động. Tôi đã thử nghiệm 6 dòng mã sau đây và nó hoạt động rất tốt:

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

Bây giờ bạn có một shapefile được viết trên bản gốc có các tiêu đề chính xác và các trường dbf gốc. Nó sẽ mở một cách an toàn trong phần mềm và thư viện shapefile của GIS nhưng không có tính năng hoặc bản ghi dbf.

Hàm lambda chuyển hộp giới hạn ban đầu dưới dạng giữ chỗ. Bạn có thể đặt giá trị float bao giờ bạn muốn trong một mảng [xmin, ymin, xmax, ymax]. Thí dụ:

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

Việc thay đổi các trường dbf cũng đơn giản và được ghi lại trong các tài liệu pyshp.

Mong rằng sẽ giúp.

Tại sao không lưu một bản sao của Shapefile trống và ghi đè lên Shapefile quan tâm của bạn.