Persimpangan beberapa lapisan dalam PostGIS


9

Saya perlu menulis prosedur tersimpan PL / pgSQL untuk perhitungan persimpangan beberapa lapisan. Misalnya, untuk tiga lapisan A, B, C fungsi ini harus mengembalikan geometri ABC: masukkan deskripsi gambar di sini

Fungsinya sebagai input id layer yang perlu disilangkan. Adakah yang bisa memberi saya saran untuk membangun fungsi ini? Pikiranku tentang ini aku harus memotong lapisan pertama dan kedua kemudian hasil persimpangan ini berpotongan dengan lapisan ketiga dll.


1
Apakah pertanyaan ini sama dengan yang Anda cari: gis.stackexchange.com/q/83/97
Nathan W

Tidak, ini pertanyaan yang berbeda.
drnextgis

Jawaban:


3

Tolong beri saya komentar tentang solusi yang saya usulkan:

CREATE OR REPLACE FUNCTION fp_intersect(lids varchar)
    RETURNS integer AS
$$
DECLARE
    lid_new  integer;
    lndx     integer := 1;
    lids_arr integer[];
BEGIN

IF regexp_replace(lids, E'\\s+', '', 'g') ~ E'^-?\\d+$' THEN
    RETURN -1;
END IF;

SELECT nextval ('g_layer_lid_seq') INTO lid_new; 

lids_arr := string_to_array(regexp_replace(lids, E'\\s+', '', 'g'), ',');
DROP TABLE IF EXISTS tmp_intersect;
CREATE TEMPORARY TABLE tmp_intersect AS SELECT geom FROM g_lgeom WHERE lid = lids_arr[1];

WHILE lndx < array_length(lids_arr, 1) LOOP
    DROP TABLE IF EXISTS tmp;
    CREATE TEMPORARY TABLE tmp AS SELECT ST_Intersection(geom, g_next) AS geom
    FROM tmp_intersect
    JOIN (
        SELECT geom AS g_next
        FROM g_lgeom
        WHERE lid = lids_arr[lndx+1]
    ) AS _
    ON ST_Intersects(geom, g_next);
    lndx := lndx+1;
    DROP TABLE IF EXISTS tmp_intersect;
    CREATE TEMPORARY TABLE tmp_intersect AS SELECT geom FROM tmp;
END LOOP;

INSERT INTO g_lgeom(lid, geom) SELECT lid_new, (_.p_geom).geom FROM (SELECT ST_Dump(geom) AS p_geom FROM tmp_intersect) AS _;

RETURN lid_new;
END
$$
LANGUAGE plpgsql;
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.