Tách chồng chéo giữa các đa giác và gán cho đa giác gần nhất bằng PostGIS?

Câu hỏi của tôi khá liên quan đến sự trùng lặp Chia không được trả lời giữa các đa giác thành các tập hợp rời rạc trừ khi tôi không sử dụng shapely, tôi đang sử dụng PostGIS 2.4 trên PostgreQuery 10.

Tôi có một bảng chứa nhiều bộ dữ liệu cho các đơn vị hành chính. Các hình học chỉ bao gồm diện tích đất. Mục tiêu của tôi là kéo dài chúng một khoảng cách nhất định xuống nước.

Để mở rộng các đa giác hiện tại vào trong nước, tôi thực hiện một bộ đệm theo sau bằng cách trừ tất cả các hình học khác trong bảng. Điều này cắt bộ đệm trong đó hình học có hàng xóm và giữ lại bộ đệm nơi lớp hiện không có gì. Vì vậy, một cái gì đó như:

SELECT c1.name
     , ST_Difference( ST_Buffer( ST_Union(c1.geom), 250 )
                    , (SELECT ST_Union(c2.geom)
                         FROM foo.county AS c2
                        WHERE c2.name <> c1.name)
                    )
  FROM foo.county AS c1
 GROUP BY c1.name;

Tuy nhiên, điều này tất nhiên mang lại cho tôi bộ đệm chồng chéo trong đó hai đa giác gặp nhau ở mép nước, cũng như trên các eo hẹp:

(Lưu ý rằng mảnh màu tím nhỏ ở bên trái phải được gán hoàn toàn cho đa giác màu đỏ, vì không có điểm nào trong hình học đó gần với màu xanh hơn màu đỏ do bán đảo giống như phụ lục.)

Đây là nơi tôi bị mắc kẹt. Tôi muốn cắt các phần chồng chéo (các phần màu tím) thành hai đa giác rời rạc dọc theo một đường thẳng tương đương với các đa giác không có bộ đệm, và sau đó trừ đi các phần tử từ các hạt được đệm.

Tôi không nghĩ rằng tôi có một trường hợp có ba hoặc nhiều bộ đệm chồng chéo. Tôi sẽ đi qua cây cầu đó nếu và khi tôi đến đó.

Nó có mùi đối với tôi giống như một cái gì đó tương tự nhưng không hoàn toàn giống với việc sử dụng ST_VoronoiPolygonsnhưng xung quanh các đa giác thay vì các điểm.

Điều này có tồn tại hay không - lý tưởng - như các chức năng của PostGIS, như chức năng xử lý địa lý của QGIS hoặc có thể trong một số thư viện python?

Nếu giải pháp cho vấn đề này có liên quan và tôi hiểu nó một cách chính xác, thì giải pháp của nó là như sau: (Dữ liệu ban đầu là trừu tượng, vì tôi hiểu nó nằm liền kề trên bộ đệm với một số mét nhất định, trong trường hợp của tôi là 2 bảng : exper_list3 và exper_list4 xem hình 1):

1) Trong pgAdmin, chạy truy vấn: create table exper_list5 as SELECT st_intersection (exper_list4.geom, exper_list3.geom) AS geom FROM exper_list4, exper_list3 WHERE st_intersects (exper_list4.geom, exper_list3.geom) = true; result: bảng exper_list5 xem hình 2;

2) Trong pgAdmin, hãy chạy truy vấn: create table exper_list6 as SELECT (ST_Dumppoints (geom)). Geom FROM exper_list5; Kết quả là một bảng có dấu chấm, ngoài ra đi đến thuộc tính bảng, tạo trường gid và đặt các hạn chế đối với nó;

3) Trong QGIS, chạy công cụ Đa giác Voronoi trên bảng (lớp vectơ) exper_list6 Kết quả được hiển thị trong Hình 3.

4) Trong pgAdmin, chạy truy vấn: create table exper_list8 as SELECT (ST_Dumppoints (geom)).geom FROM exper_list7; Kết quả được hiển thị trong Hình 4.

5) Chọn từ exper_list8 những điểm nhập đa giác exper_list5 Trong pgAdmin, chạy truy vấn: create table exper_list9 as SELECT st_intersection (exper_list8.geom, exper_list5.geom) AS geom FROM exper_list8, exper_list5 WHERE st_intersects (exper_list8.geom, exper_list5.geom) = true; và sao chép các điểm này vào kết quả bảng exper_list6 hình 5;

6) Trong QGIS Trên bảng exper_list6 chạy công cụ "Delaunay triangulation" và nhận kết quả exper_list10 xem hình 6

7) Trong QGIS ở chế độ chỉnh sửa "Bút chì" trong bảng exper_list10, hãy chọn những hình tam giác chạm vào vùng quan tâm, sử dụng công cụ "Vector" - "Geoprocessing" - "Thống nhất theo tính năng" (kiểm tra "chỉ các đối tượng được chọn" "Phân loại" - bằng tất cả các dấu hiệu) kết quả exper_list11 hình 7;

8) Từ exper_list6 cắt exper_list11 kết quả hình 8 exper_list12;

9) Từ exper_list4 cắt exper_list11 kết quả hình 9 exper_list13;

10) Các hành động tương tự phải được thực hiện với đa giác bên trái (trong trường hợp của tôi, bảng exper_list4)

11) Kết quả chung là hai bảng exper_list13 và exper_list14 xem hình 10.

chúc may mắn, bạn trung thành