Tạo bộ đệm hòa tan từ đa hình học (Liên kết theo thuộc tính chia sẻ và giao điểm không gian)

Tôi đã tạo bộ đệm hòa tan từ các tính năng đầu vào đa điểm. Trong ví dụ dưới đây, bảng đầu vào chứa 4 tính năng. Tính năng #2bao gồm hai hình học điểm. Sau khi tạo bộ đệm, tôi nhận được 4 hình học đa giác:

Có cách nào để nhóm kết quả? Bộ đệm của các điểm #1và #2được hòa tan và phải là một tính năng đa giác đơn ( a).

Những gì tôi đã làm cho đến nay:

-- collect all buffers to a single multi-polygon feature
-- dissolve overlapping polygon geometries
CREATE TABLE public.pg_multibuffer AS SELECT
    row_number() over() AS gid,
    sub_qry.*
FROM (SELECT
    ST_Union(ST_Buffer(geom, 1000, 8))::geometry(MultiPolygon, /*SRID*/) AS geom
FROM
public.multipoints)
AS sub_qry;

BIÊN TẬP:

-- create sample geometries

CREATE TABLE public.multipoints (
gid serial NOT NULL,
geom geometry(MultiPoint, 31256),
CONSTRAINT multipoints_pkey PRIMARY KEY (gid)
);

CREATE INDEX sidx_multipoints_geom
ON public.multipoints
USING gist
(geom);

INSERT INTO public.multipoints (gid, geom) VALUES
(1, ST_SetSRID(ST_GeomFromText('MultiPoint(12370 361685)'), 31256)),
(2, ST_SetSRID(ST_GeomFromText('MultiPoint(13520 360880, 19325 364350)'), 31256)),
(3, ST_SetSRID(ST_GeomFromText('MultiPoint(11785 367775)'), 31256)),
(4, ST_SetSRID(ST_GeomFromText('MultiPoint(19525 356305)'), 31256));

Bắt đầu với một số điểm ngẫu nhiên, trong một nhận thức bắt chước những điểm trong hình ảnh của OP, trong đó hai điểm đầu tiên giao nhau về mặt không gian, sau đó điểm thứ 2 và thứ 3 có cùng thuộc tính id (2), với một vài điểm khác không giao nhau về mặt không gian và không có cùng thuộc tính, truy vấn sau đây tạo ra 3 cụm:

WITH 
  temp (id, geom) AS 
     (VALUES (1, ST_Buffer(ST_Makepoint(0, 0), 2)),
        (2, ST_Buffer(ST_MakePoint(-0.7,0.5), 2)),
        (2, ST_Buffer(ST_MakePoint(10, 10), 2)), 
        (3, ST_Buffer(ST_MakePoint(-2, 12), 2)), 
        (4, ST_Buffer(ST_MakePoint(5, -6), 2))),
 unions(geoms) AS 
      (SELECT ST_Union(geom) FROM temp GROUP BY id),
 clusters(geoms) AS 
      (SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(geoms)), 3) 
         FROM unions),
 multis(id, geoms) AS 
      (SELECT row_number() over() as id, geoms FROM clusters)
 SELECT ST_UNION(d.geom) FROM 
      (SELECT id, (ST_DUMP(geoms)).geom FROM multis) d GROUP BY id;

Có một số bước ở đây:

1. sử dụng ST_Union, nhóm theo id, để nhóm đầu tiên theo thuộc tính
2. sử dụng ST_ClusterIntersectingđể kết hợp những người từ cùng một nhóm giao nhau trong không gian
3. thêm một id vào mỗi cụm (bảng multis) - cố gắng thực hiện điều này trực tiếp trong ClusterIntersecting dẫn đến tất cả các hình học nhận được id là 1
4. Liên kết các hình học bị đổ từ bước 2, nhóm theo id từ bước 3 - đây là phần hòa tan . Điều này làm cho hai đa giác chồng chéo trong cụm A của bạn, được nối với nhau, thay vì chồng chéo, vì chúng ở cuối bước 2.

Khá dài, nhưng nó hoạt động (và, tôi chắc chắn có một cách ngắn hơn).

Sử dụng công cụ WKT trong QGIS, (và khám phá mức độ khủng khiếp của tôi với các công cụ chỉnh sửa) tạo ra các cụm như sau, trong đó bạn có thể thấy cụm được gắn nhãn là a, hoàn toàn cùng nhau - tức là một màu.

Nếu bạn đặt ST_AsText vào vòng chung kết, ST_UNION (d.geom), thì bạn có thể xem kết quả trực tiếp.

EDIT theo dõi thêm thông tin trong các bình luận: Khi bạn bắt đầu với các điểm, bạn sẽ cần kết hợp bộ đệm vào giải pháp ban đầu của tôi - mà tôi đã đưa vào CTE tạm thời khi bắt đầu mô phỏng sơ đồ của bạn. Sẽ dễ dàng hơn để thêm bộ đệm trong các hiệp hội CTE, vì vậy bạn có thể thực hiện tất cả các hình học cùng một lúc. Vì vậy, sử dụng khoảng cách bộ đệm là 1000, làm ví dụ, sau đây trả về 3 cụm, như mong đợi.

WITH temp(id, geom) AS 
  (VALUES 
      (1, ST_SetSRID(ST_GeomFromText('MultiPoint(12370 361685)'), 31256)),   
      (2, ST_SetSRID(ST_GeomFromText('MultiPoint(13520 360880, 19325 364350)'), 31256)),                                                
      (3, ST_SetSRID(ST_GeomFromText('MultiPoint(11785 367775)'), 31256)),
      (4, ST_SetSRID(ST_GeomFromText('MultiPoint(19525 356305)'), 31256))
),                                              
unions(geoms) AS 
  (SELECT st_buffer(ST_Union(geom), 1000) FROM temp GROUP BY id),
clusters(geoms) AS 
  (SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(geoms)), 3) 
     FROM unions),
multis(id, geoms) AS 
  (SELECT row_number() over() as id, geoms FROM clusters)
SELECT id, ST_UNION(d.geom) FROM 
  (SELECT id, (ST_DUMP(geoms)).geom FROM multis) d GROUP BY id;

Một cách để làm điều này là kết hợp ST_Uniontất cả các bộ đệm với nhau, ST_Dumpkết quả để có được các thành phần của đa giác kết quả và nối ST_Intersectslại với các điểm đầu vào để tìm ra có bao nhiêu / điểm tạo thành mỗi cụm.

Điều này có thể được thực hiện mà không yêu cầu tham gia bằng cách nhóm các điểm lại với nhau trước khi gọi ST_Buffer. Để hai điểm được đặt trong cùng một bộ đệm hòa tan, chúng phải có thể truy cập được bằng các bước nhảy giữa các điểm có khoảng cách nhỏ hơn eps. Đây chỉ là một vấn đề phân cụm liên kết tối thiểu, có thể được giải quyết bằng cách sử dụng ST_ClusterDBSCAN:

SELECT
  cluster_id,
  ST_Union(ST_Buffer(geom, 1000)) AS geom,
  count(*)                        AS num_points,
  array_agg(point_id)             AS point_ids
FROM (
  SELECT
    point_id,
    ST_ClusterDBSCAN(geom, eps := 2000, minpoints := 1) OVER() AS cluster_id ,
    geom
  FROM points) sq
 GROUP BY cluster_id;

Lưu ý rằng điều này sẽ không tạo ra kết quả chính xác như phương pháp đệm đầu tiên, bởi vì bộ đệm PostGIS không phải là vòng tròn hoàn hảo và hai điểm cách nhau 1000m có thể không được kết nối bởi hai bộ đệm 500m.

Theo câu trả lời này, bạn muốn thực hiện ST_DUMP trong truy vấn con của mình.

Một cái gì đó như thế này:

-- collect all buffers to a single multi-polygon feature
-- dissolve overlapping polygon geometries
CREATE TABLE public.pg_multibuffer AS SELECT
    row_number() over() AS gid,
    sub_qry.*
FROM (SELECT
    ST_Dump(ST_Union(ST_Buffer(geom, 1000, 8))::geometry(MultiPolygon, /*SRID*/)) AS geom
FROM
public.multipoints)
AS sub_qry;

Lý do tại sao điều đó ST_UNIONtrả về một đa giác hòa tan của tất cả các tính năng và ST_DUMPphá vỡ điều này thành các tính năng đa giác riêng lẻ (đã bị hòa tan).