Làm thế nào để tăng tốc độ phân vùng không gian trong postgis?

Tôi có một loạt các đa giác chồng chéo và tôi đang cố gắng phân vùng không gian để tránh có các đa giác chồng chéo. Tôi nghĩ vấn đề của tôi khá đơn giản. Sử dụng một số sản phẩm ESRI và http://arcscripts.esri.com/details.asp?dbid=16700 đồng nghiệp của tôi đã tính toán nó trong 48 giây.

Tôi đang cố gắng thực hiện với postgis bằng cách sử dụng http://s3.opengeo.org/postgis-power.pdf#page=24 (đoán chi tiết, sử dụng http://trac.osgeo.org/postgis/wiki/UsersWikiSimplifyPreserveTopology như một nguồn cảm hứng) nhưng nó chậm đến mức tôi không thể làm điều đó với hơn 10 polys (tôi đã có 800 trong số chúng để chia tách). Phần chậm là ST_Union, tôi đã thử nhiều thứ, nhưng không có nơi nào thành công, đây là tình trạng hiện tại:

select geom from
(select st_linemerge(st_union(geom)) as geom from
    (select st_exteriorring((st_dumprings((st_dump(t.geom)).geom)).geom) as geom from
        (SELECT geometry AS geom, id
               FROM tt
              WHERE campaign_id = 204
              ORDER BY id limit 200) t) t2) t3

điều này đã được tính toán trong 26 phút (linTablege () thực tế là không). Các polys là MultiPolygons trong trường hợp st_dump làm phiền bạn.

Bạn có bất cứ lời khuyên? St_union () của linework là phần rất chậm.

Cảm ơn,

Nico.

PS: đây là một số con số: 852 đa cực, dẫn đến 14880 đa giác, dẫn đến 21467 linestrings có tổng số 315513 đỉnh.

Câu trả lời này có thể không giúp @nraynaud trực tiếp, nhưng hy vọng sẽ làm sáng tỏ chủ đề này.

Có một vấn đề tương tự trong spatiaLite <4.0 do vấn đề với GEOS. Xem liên kết này để thảo luận về vấn đề này.

Cách giải quyết là thay thế hàm ST_Union () bằng ST_UnaryUnion (ST_Collect ()). Thật không may, ST_UnaryUnion không có sẵn cho đến postGIS 2.0 (theo như tôi có thể nói.)

Phiên bản nào của PostGIS bạn đang sử dụng? Liên kết chậm hơn nhiều nếu bạn đang sử dụng PostGIS <1.4 hoặc GEOS <3.2. Liên minh nhanh hơn nhiều đã được giới thiệu trong 1.4, nhưng cũng yêu cầu GEOS 3.2+. Vì vậy, trước tiên nếu bạn đang sử dụng thấp hơn 1,4, tôi sẽ nâng cấp lên ít nhất 1,5.

SELECT postgis_full_version();

Để kiểm tra.

Cũng là ý định của bạn để giữ các cạnh ban đầu của đa giác. Nếu bạn chỉ muốn giải thể các khu vực chồng chéo,

SELECT ST_Union(geom) FROM tt WHERE campaign_id = 204;

Sẽ làm các mẹo.