Çok geometriden çözünmüş tampon oluşturma (Paylaşılan öznitelik ve uzamsal kesişime göre birleştirme)

Çok noktalı giriş özelliklerinden çözünmüş tamponlar oluşturmam gerekiyor. Aşağıdaki örnekte, giriş tablosu 4 özellik içermektedir. Özellik #2iki nokta geometrisinden oluşur. Bir tampon oluşturduktan sonra, 4 çokgen geometrisi elde ediyorum:

Sonucu gruplamanın bir yolu var mı? Puan tamponlar #1ve #2eritilir ve tek bir çoklu-çokgen özelliği olmalıdır ( a).

Şimdiye kadar yaptığım:

-- 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;

DÜZENLE:

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

Bazı rasgele noktalarla başlayarak, OP görüntüsünde, ilk ikisinin mekansal olarak kesiştiği yerlere taklit etme girişiminde, 2. ve 3. aynı özellik kimliğine (2) sahip, mekansal olarak kesişmeyen veya sahip olmayan birkaç nokta daha var. aynı öznitelikte, aşağıdaki sorgu 3 küme üretir:

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;

Burada birkaç adım var:

1. kullanmak ST_Unionözelliğe göre ilk gruba, kimliğine göre, gruplama
2. ST_ClusterIntersectingmekansal olarak kesişen aynı gruptaki insanları birleştirmek için kullanın
3. kümelerin her birine bir kimlik ekleyin (tablo multis) - bunu doğrudan ClusterIntersecting'te yapmaya çalışmak, 1
4. 2. adımdaki dökümlenmiş geometrileri, 3. adımdaki kimliğe göre gruplandırın - bu çözülme kısmıdır. Bu, küme A'nızdaki üst üste gelen iki çokgenin üst üste gelmek yerine adım 2'nin sonunda olduğu gibi birbirine birleştirilmesine neden olur.

Oldukça uzun, ama işe yarıyor (ve daha kısa bir yol olduğuna eminim).

QGIS'de WKT aracını kullanmak (ve düzenleme araçlarıyla ne kadar korkunç olduğumu keşfetmek), a olarak etiketlenmiş kümeyi görebileceğiniz aşağıdaki gibi kümeler üretir - yani tek renk.

Bir ST_AsText finaline ST_UNION (d.geom) koyarsanız, sonuçları doğrudan görebilirsiniz.

Yorumlarda daha fazla bilgi DÜZENLEYİN: Noktalarla başlarken tamponu orijinal çözümüme dahil etmeniz gerekecek - bu da şemayı taklit etmek için başlangıçta geçici CTE'ye koydum. Tamponları CTE sendikalarına eklemek daha kolay olacaktır, böylece tüm geometrileri bir kerede yapabilirsiniz. Dolayısıyla, örnek olarak 1000 arabellek mesafesi kullanıldığında, aşağıdakiler artık beklendiği gibi 3 küme döndürmektedir.

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;

Bunu yapmanın bir yolu, ST_Uniontüm tamponlara birlikte, ST_Dumpsonuçtaki poligonun bileşenlerini elde etme ve ST_Intersectsher bir kümede kaç tane / hangi noktayı oluşturduğunu bulmak için giriş noktalarına geri birleşmektir .

Bu, çağrı yapmadan önce noktaları bir araya getirerek birleştirme gerektirmeden yapılabilir ST_Buffer. Aynı çözünmüş tampon içinde iki noktanın bulunabilmesi için, bunlardan daha kısa bir mesafe arasındaki sıçramalarla erişilebilmelidir eps. Bu, yalnızca aşağıdakiler kullanılarak çözülebilen minimum bağlantı kümeleme sorunudur 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;

PostGIS arabellekleri mükemmel daireler olmadığından ve 1000m aralıkla iki nokta iki 500m arabellekle bağlanamayabileceğinden, bunun arabellek ilkiyle aynı sonucu vermeyeceğini unutmayın.

Bu cevaba göre alt sorgunuzda ST_DUMP yapmak istiyorsunuz.

Bunun gibi bir şey:

-- 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;

Bunun nedeni ST_UNION, tüm özelliklerin çözünmüş bir ST_DUMPçokgenini döndürür ve bunu ayrı ayrı çokgen özelliklerine (çözülmüş olan) ayırır.