PostGIS kullanarak daha büyük çokgen oluşturmak için çok sayıda küçük çokgen birleştirmek?

Postgis'te SRID 27700 kullanarak aşağıdaki katman var:

Birleşik Krallık'taki her idari bölge ve (renk gruplarından da görebileceğiniz gibi) her birinin içinde bulunduğu bölgeyi belirten bir metin alanı var.

Yapmak istediğim, belirli bir ilçedeki küçük ilçelerden daha büyük ilçeli çokgenler yapmaktır, bu nedenle, tüm deniz mavisi renk çokgenlerinin üzerindeki resimde EG, buradaki tüm kutupları içeren tek dış halkadan büyük bir çokgen oluşturur. renk, mesela tamamen mor, kahverengi, pembe, gri vb. hepsi bir poligon oluşturmalı.

Ben zaten aşağıdakileri denedim:

insert into parishesmerged (geometry)
select astext(multi(ST_Union(the_geom))) as the_geom from parishes
group by county_name

Fakat daha sonra işlemede büyük problemler yaşadığım kırılmış geometriler oluşturuyor.

Başlıca çıktı alanlarının bulunduğu daha basit bir ilçe düzeyinde harita yapmaya çalışıyorum.

Herhangi bir çözümün Postgis'te de olması gerekmiyor, tam OS4Geo yığınının yüklü olması, QGis'in en son sürümü ve bir çubuğumu sallayabileceğimden daha fazla dosya kullanıyorum.

Sahip olmadığım tek şey ArcGis gibi büyük çocuklar (Etrafta bir yerde uzanmış eski bir Mapinfo olsa da)

Kayıt için oluşturmaya çalıştığım veri kümesi, .NET kullanarak GIS uygulamaları yazmak isteyen .NET programcılarına yönelik yazdığım bir GIS kitabına eşlik etmektir.

Aşağıdaki önerileri denedikten sonra, en iyi sonuç veren, 'Paul Ramsey' çözümü oldu.

Şimdi kitabım için yeterince basit, ancak ilginç bir coğrafi-mekansal SQL göstermeme izin verecek kadar karmaşık, hoş bir basitleştirilmiş ilçe ve ilçe dosyasına sahibim.

Sonunda Paul'un çözümü benim için çalışan olsa da, poligon haritasını basitleştirmek ve karmaşıklığı daha da azaltmak gibi şeylerin diğer cevaplarını da çizdim.

Ancak bunu yaparken gözlemlediğim şey, ST_Collect gerçekten ST_Union'dan daha hızlı olsa da, koşmak için koşmak aynı zamanda kırılmış geometrilerden en çok sorumlu olanıydı. Tahminimce hız artışı çekirdek fonksiyonunda daha az doğruluk pahasına elde ediliyor.

ST_Union işe yarayacaktı ama işin neredeyse kesin olarak temiz değildi. Demek ki küçük şeylerinin sınırları tamamen mükemmel değil. Köşelerin hizaya girme ihtimalini artırmak için onları yavaşça bir ızgaraya bağlayabilirsiniz, ancak bahse girerim hala işe yaramayan birkaç vakanız olacak. Ya toleransların ötesinde olacaklar ya da daha muhtemel olarak, tepe noktalarının eşleştirilmediği yerler olacak, yani bir tarafta bir çizgi ve diğer tarafta bir köşe olacak.

 CREATE TABLE merged AS
 SELECT ST_Union(ST_SnapToGrid(the_geom,0.0001)) 
 FROM parishes
 GROUP BY county_name;

PostGIS 2.0'a sahipseniz, toleranslı bir topoloji yapısı inşa etmek, aradığınız cevaba, eğer biraz şansınız varsa, getirebilir.

Diyelim ki, "... tüm halkaları içeren tek dış halkadan büyük bir çokgen ...". ST_ExteriorRing bunu yapar,

SELECT ST_MakePolygon(ST_ExteriorRing(ST_Union(GEOM)))
FROM GEOMTABLE GROUP BY ATTRCOLUMN

ST_Union () öğesini önerildiği gibi kullanabilir veya ST_Collection () ile test edebilirsiniz.

NOTLAR: Küçük halkalardan veya "kırılmış geometrilerden" kaçınmak için her geom için st_convexhull ve / veya ST_Simplify kullanabilirsiniz,

SELECT ST_MakePolygon(ST_ExteriorRing(ST_union(ST_Simplify(GEOM,0.5))))
FROM GEOMTABLE GROUP BY ATTRCOLUMN

ve geometrilerinizi kontrol edin,

SELECT * FROM (
   SELECT gid, ST_IsValid(geom) as valid, ST_IsSimple(geom) as simple 
   FROM GEOMTABLE) AS t  
WHERE NOT(valid AND simple); 

ST_Collect işlevi PostgreSQL terminolojisinde bir "toplama" işlevidir.

" SELECT ST_Collect(GEOM) FROM GEOMTABLE GROUP BY ATTRCOLUMN", ATTRCOLUM’un her bir değeri için ayrı bir GEOMETRY KOLEKSİYONU döndürecek

http://postgis.net/docs/ST_Collect.html

Not: ST_Collect, ST_Union'dan çok daha hızlıdır

Sorunuzdan Ordnance Survey'in Sınır Çizgisi ürününü kullandığınızı varsayıyorum. Eğer durum buysa, o zaman zaten bir İlçe düzeyinde veri seti içermektedir, bu yüzden onu daha düşük seviye bölge alanlarından kendiniz üretmeye çalışmanıza gerek yoktur.

Sınır Çizgisi'ni kullanmıyorsanız, OS OpenData lisansı altında ücretsiz olduğundan ve doğrudan PostGIS'e yükleyebileceğiniz bir şekil dosyası olarak İlçe düzeyine sahip olmanızı öneririm.