Làm thế nào để kiểm tra theo chương trình nếu số lượng hình = số bản ghi bảng?

Tôi có một số khoảng 1000 shapefile bị hỏng (xem thông báo lỗi đính kèm). Các shapefile được tạo từ eCognition Developer 8. Có một công cụ tập lệnh dường như sửa chữa shapefile một khi nó được xác định là bị hỏng.

Biên tập:

Tôi muốn tạo một tập lệnh nhanh để lặp qua tất cả các shapefile của mình và kiểm tra xem số lượng hình có khớp với các bản ghi bảng không. Tôi có thể đếm các bản ghi bảng bằng cách sử dụng như sau:

# 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

Cuối cùng tôi muốn tạo ra một số loại kiểm tra logic như:

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

Làm cách nào để đếm hình trực tiếp mà không cần truy cập tệp .dbf? Hay nói cách khác, cách tốt nhất để kiểm tra theo chương trình nếu số lượng hình phù hợp với số lượng bản ghi bảng là gì?

Còn việc sử dụng pyshp thì sao? Tôi đã cài đặt nó với pip và những gì tôi đã thử bên dưới khá nhiều so với README :

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

Thật không may (hoặc có thể may mắn thay?) Tôi không có bất kỳ shapefile jack-up nào để kiểm tra xem có không. của hình dạng có thể! = không. của hồ sơ.

Đợi một chút, bây giờ tôi đã có một shapefile jacked nhờ ý tưởng của Kirk trong các ý kiến ​​dưới đây. Tôi đã sao lưu dbf, tạo một bản sao của toàn bộ shapefile, xóa một số tính năng, sau đó đổi tên dbf sao lưu trở lại ban đầu, và lo và kìa, số lượng hình dạng <số lượng bản ghi:

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

Từ âm thanh của câu hỏi của bạn, có vẻ như tất cả những gì bạn thực sự muốn làm là xác định xem một shapefile có vấn đề với nó hay không (trong trường hợp này là các bản ghi không khớp). Nếu bạn chỉ cần xác định những người có vấn đề, bạn thực sự không cần phải đếm các bản ghi trong DBF và Shapefile để xác định xem nó có bị lỗi hay không. Đây là lý do tại sao:

Nếu bạn cố chạy chức năng GetCount trên một shapefile có số lượng bản ghi khác nhau, nó sẽ thất bại với lỗi:

ERROR 000.229 : Không thể mở. Không thể thực thi (GetCount).

Vì hàm GetCount không thành công trong kịch bản này và tất cả những gì bạn muốn làm là xác định các shapefile bị lỗi, bạn có thể bắt gặp điều này với mệnh đề try / trừ trong mã của bạn, thay vì if / other trước đây bạn đang cố sử dụng.

Tôi đã tự do thêm mã và vòng lặp "Danh sách tính năng" để bạn có thể kiểm tra tất cả các FC trong không gian làm việc của mình mà không phải kiểm tra thủ công từng cái.

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

Các định dạng shapefile được ghi lại. Tôi đoán số lượng bản ghi trong tệp shp không tương ứng với số lượng bản ghi trong tệp dbf.

Định dạng tệp shp được ghi lại ở đây . Vì vậy, bạn có thể viết một chương trình để đếm số lượng hình dạng. Định dạng dbf được ghi lại ở nhiều nơi và bạn sẽ có thể tìm thấy các mẫu để đếm hàng, ví dụ ở đây .

Tập lệnh đính kèm lặp qua một thư mục và kiểm tra xem số lượng hình có khớp với số lượng bản ghi cho mỗi shapefile không.

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"

Sử dụng hình học kiểm tra sẽ giúp bạn vượt qua bước đầu tiên.
Onus
Repair Geometry sẽ cho phép bạn chọn thứ tự và mức độ ưu tiên của vấn đề bạn muốn sửa chữa.
đây là một số liên kết phiên bản cũ hơn . Khi bạn chạy trình kiểm tra shapefile, sau đó bạn kết thúc với việc xây dựng lại dbf?
Đó là bước tạo ra các hồ sơ để phù hợp. Một trong hai điều đã xảy ra để gây ra lỗi.

1. Shp có một đối tượng (không gian) đã bị xóa / bỏ bởi một phần mềm / quy trình khác.
2. Dff có một bản ghi tham chiếu hình học null.
   Một số điều có thể gây ra điều này.
   Các shx thực sự là chỉ số giữa hai.
   Đếm các hình dạng mà không đếm các bản ghi dbf chỉ là một nửa của giải pháp.

Nhìn vào bài viết wikipedia về shapefiles , tệp .shx nên chứa một chỉ mục trên tệp .shp, không phải trên tệp .dbf. Vì vậy, có thể cần phải kiểm tra xem .shx và .shp có khớp với nhau không.

Có thể mở một shapefile mà không cần .dbf (có nghĩa là bạn không có bảng thuộc tính), nhưng một chỉ mục bị hỏng sẽ tạo ra một thông báo lỗi.