Hòa tan đa giác trong QGIS (GEOS) với dung sai

Tôi có một lớp vectơ đa giác đại diện cho khu vực ngập trong một trận lụt. Lớp này đã được tạo ra theo cách mà một số khu vực được tạo thành từ nhiều đa giác chạm nhỏ hơn (xem ảnh chụp màn hình bên dưới, được gắn nhãn "dữ liệu gốc").

Tôi muốn hòa tan các đa giác nơi chúng chạm vào. Nếu tôi sử dụng công cụ xử lý địa lý hòa tan trong fTools của QGIS, sử dụng thư viện GEOS, tôi sẽ nhận được một số kết quả không mong muốn (xem "hòa tan với QGIS fTools"). Có vẻ như GEOS không sử dụng ngưỡng trong việc xác định xem đa giác có chạm vào nhau không, như đã đề cập trong các vấn đề này:

http://trac.osgeo.org/qgis/ticket/3126

http://trac.osgeo.org/qgis/ticket/3184

http://hub.qgis.org/issues/2806

Cách tiếp cận tương tự trong ArcGIS cho kết quả như mong đợi (xem "hòa tan với ArcGIS 10.1").

Rõ ràng đây là một "tính năng" của GEOS, chứ không phải là một lỗi. Có cách nào để sửa dữ liệu tôi phải tạo ra kết quả như mong đợi không? Một sự mất mát nhỏ của độ chính xác (một vài cm?) Có thể được chấp nhận đối với tôi. Lý tưởng nhất là giải pháp tôi có thể viết vào plugin QGIS, như một bước xử lý trước.

Câu hỏi liên quan (không có câu trả lời): Làm cách nào để đặt giá trị dung sai của công cụ Geoprocessing trong QGIS?

Liên kết với dữ liệu mẫu dưới dạng shapefile: https://www.dropbox.com/s/c8ns7wwfv289ad0/sample.zip

Trong trường hợp này, có vẻ như bạn cần trợ giúp công cụ hòa tan bằng cách sửa cấu trúc liên kết trước. Đây là cách bạn có thể làm điều này trong GRASS.

Khi bạn tải shapefile, đây là cách GRASS nhìn thấy nó:

Mỗi chữ thập đánh dấu một khu vực. Bạn có thể thấy rằng có một số đường chéo trên các đường biên: Đây là các đa giác cúi do độ lệch nhỏ giữa các đường biên.

1. Loại bỏ các đa giác mảnh được hình thành bởi các khoảng trống giữa các đa giác thực:

   v.clean input=sample output=tempA --o tool=rmarea thresh=0.001

2. Chụp các đường biên:

   v.clean input=tempA output=tempB --o tool=snap thresh=0.001

3. Phá vỡ các đường ranh giới tại các giao lộ mới:

   v.clean input=tempB output=tempA --o tool=break thresh=0.001

4. Xóa các đường biên bây giờ trùng lặp:

   v.clean input=tempA output=tempB --o tool=rmdupl thresh=0.001

5. Hòa tan các đa giác

   v.dissolve input=tempB output=cleaned --o column=ID

Et voilà:

Sử dụng lệnh GRASS GIS v.dissolve (trong plugin Sextante -> Lệnh GRASS -> Vector)

Nếu tôi nhanh chóng sao chép ví dụ của bạn:

Kết quả với v.dissolve:

Giải trình

Vấn đề rất dễ hiểu. Một shapefile không có cấu trúc liên kết, nếu hai khu vực có chung đường viền chung thì đường viền sẽ được số hóa hai lần và cũng được lưu trữ trùng lặp.

Bạn thấy rõ điều này khi bạn cố gắng chỉnh sửa shapefile:

GRASS GIS (như ArcGIS) là một gis tô pô. Khi bạn nhập shapefile trong GRASS, lớp kết quả có cấu trúc liên kết: đường viền chung giữa hai khu vực tồn tại một lần và được chia sẻ.

Vì lý do này, khi bạn muốn hòa tan một shapefile ( hòa tan hai ranh giới ), bạn có thể gặp phải các vấn đề (trong QGIS, GvSIG, OpenJUMP hoặc ArcView 3.x cũ), xem ArcUserI: Hiểu về cấu trúc liên kết và Shapefiles .

Trong GRASS, v.dissolve được thiết kế để xóa tan ranh giới chung giữa các khu vực có cùng số danh mục