เข้าร่วมรูปหลายเหลี่ยมขนาดเล็กจำนวนมากเพื่อสร้างรูปหลายเหลี่ยมที่มีขนาดใหญ่ขึ้นโดยใช้ PostGIS

ฉันมีเลเยอร์ต่อไปนี้ใช้ SRID 27700 ใน postgis:

เป็นเขตการปกครองทุกแห่งในสหราชอาณาจักรและ (ตามที่คุณเห็นจากการจัดกลุ่มสี) แต่ละเขตข้อมูลมีเขตข้อมูลข้อความระบุเขตที่พวกเขาอยู่

สิ่งที่ฉันต้องการทำคือการทำให้รูปหลายเหลี่ยมของมณฑลใหญ่ขึ้นจากรูปเล็กในเขตที่กำหนดดังนั้น EG ในภาพด้านบนรูปหลายเหลี่ยมสีนกเป็ดน้ำจะกลายเป็นรูปหลายเหลี่ยมขนาดใหญ่หนึ่งรูปจากวงแหวนด้านนอกเดียวที่มี polys ทั้งหมดในนั้น สีเช่นสีม่วง, น้ำตาล, ชมพู, เทา ฯลฯ อย่างชาญฉลาดควรรวมกันเป็นหนึ่งเดียว

ฉันได้ลองทำสิ่งต่อไปนี้แล้ว:

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

แต่มันก็สร้างรูปทรงที่แตกซึ่งฉันมีปัญหาใหญ่ในการประมวลผลต่อไป

ฉันกำลังพยายามทำแผนที่ระดับเคาน์ตีที่เรียบง่ายขึ้นด้วยพื้นที่เอาท์พุทหลักใน

โซลูชันใด ๆ ก็ไม่จำเป็นต้องอยู่ใน Postgis เช่นกันฉันได้ติดตั้งสแต็ค OS4Geo เต็มเวอร์ชั่นล่าสุดของ QGis และอุปกรณ์อื่น ๆ มากกว่าที่ฉันสามารถจับได้

สิ่งเดียวที่ฉันไม่มีคือชายใหญ่อย่าง ArcGis (แม้ว่าฉันอาจมี Mapinfo แบบเก่าอยู่รอบ ๆ )

สำหรับบันทึกชุดข้อมูลที่ฉันพยายามสร้างคือการรวบรวมหนังสือ GIS ที่ฉันกำลังเขียนโดยมีวัตถุประสงค์เพื่อโปรแกรมเมอร์. NET ที่ต้องการเขียนแอปพลิเคชัน GIS โดยใช้. NET

หลังจากลองทำตามคำแนะนำด้านล่างสิ่งที่ได้ผลดีที่สุดคือทางออกของ 'Paul Ramseys'

ตอนนี้ฉันมีไฟล์มณฑลและเมืองที่เรียบง่ายที่ดีพอสำหรับหนังสือของฉัน แต่ซับซ้อนพอที่จะให้ฉันแสดง SQL ทางภูมิศาสตร์เชิงพื้นที่ที่น่าสนใจ

แม้ว่าทางออกของ Paul ในท้ายที่สุดจะเป็นสิ่งที่เหมาะกับฉัน แต่ฉันก็ตอบคำถามอื่น ๆ เช่นการทำให้แผนที่รูปหลายเหลี่ยมง่ายขึ้นและลดความซับซ้อนลงอีก

ในสิ่งที่ฉันได้สังเกตในขณะที่ทำสิ่งนี้อย่างไรก็ตามในขณะที่ ST_Collect เร็วกว่า ST_Union แน่นอนการวิ่งเพื่อวิ่งก็เป็นสิ่งที่รับผิดชอบต่อรูปทรงเรขาคณิตที่แตกหัก ฉันเดาว่าการเพิ่มความเร็วนั้นทำได้ด้วยความแม่นยำน้อยกว่าในฟังก์ชั่นคอร์

ST_Union จะทำงานได้ แต่งานของคุณเกือบจะไม่สะอาด ดังนั้นขอบเขตของสิ่งเล็ก ๆ น้อย ๆ ของคุณจึงไม่เหมือนกันอย่างสมบูรณ์ คุณสามารถสแน็ปอินไปที่กริดเพื่อลองและเพิ่มอัตราต่อรองที่แนวตั้ง แต่ฉันเดิมพันว่าคุณจะยังมีบางกรณีที่ไม่ได้ผล ไม่ว่าพวกเขาจะอยู่เหนือความอดทนหรือมีโอกาสมากขึ้นที่จะมีสถานที่ที่จุดยอดไม่ได้ถูกจับคู่ดังนั้นจึงมีเส้นข้างหนึ่งและจุดยอดอีกด้านหนึ่ง

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

หากคุณมี PostGIS 2.0 การสร้างโครงสร้างโทโพโลยีด้วยความอดทนอาจทำให้คุณได้คำตอบที่คุณกำลังมองหาหากคุณโชคดี

คุณบอกว่าต้อง "" สร้างรูปหลายเหลี่ยมขนาดใหญ่หนึ่งรูปจากวงแหวนรอบนอกเดียวที่มี polys ทั้งหมด ... " ST_ExternalRing ทำเช่นนี้

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

คุณสามารถใช้ ST_Union () ตามที่แนะนำหรือทดสอบกับ ST_Collection ()

หมายเหตุ: เพื่อหลีกเลี่ยงการเกิดอุบัติเหตุเล็ก ๆ น้อย ๆ หรือ "รูปทรงเรขาคณิตที่แตกหัก" คุณสามารถใช้ st_convexhull และ / หรือ ST_Simplify สำหรับแต่ละ geom

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

และตรวจสอบรูปทรงเรขาคณิตของคุณ

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เป็นฟังก์ชัน "รวม" ในคำศัพท์ของ PostgreSQL

" SELECT ST_Collect(GEOM) FROM GEOMTABLE GROUP BY ATTRCOLUMN" จะส่งกลับ GEOMETRYCOLLECTION แยกต่างหากสำหรับแต่ละค่า ATTRCOLUM ที่แตกต่างกัน

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

หมายเหตุ: ST_Collect เร็วกว่า ST_Union มาก

ฉันสมมติว่าจากคำถามของคุณว่าคุณกำลังใช้ผลิตภัณฑ์เขตแดนของ Ordnance Survey หากเป็นกรณีนี้จะมีชุดข้อมูลระดับเขตอยู่แล้วดังนั้นจึงไม่จำเป็นต้องพยายามสร้างด้วยตนเองจากพื้นที่ตำบลระดับต่ำกว่า

หากคุณไม่ได้ใช้ Boundary-Line ฉันขอแนะนำให้คุณทำเพราะมันฟรีภายใต้ใบอนุญาต OS OpenData และมีระดับ County เป็นไฟล์รูปร่างคุณสามารถโหลดลงใน PostGIS ได้โดยตรง