Saya sedang memeriksa cara terbaik untuk memotong linestrings dengan poin.
Skenarionya adalah: banyak jalan, perlu bagian dipotong oleh titik potong, semacam ini:
saya mendapatkan
tabel linestrings (penuh tanpa poin)
tabel poin st_intersection
Saya perlu memotong segmen linestring independen oleh tabel titik persimpangan.
Saya menggunakan fungsi PostGIS, dan menemukan beberapa pendekatan, tetapi masing-masing dari mereka telah memberi saya semacam masalah.
Inilah yang sudah saya uji:
1
Tabel baris: 1 baris, st_memunion dari 1200 baris Tabel poin: 1700 baris (poin)
Yang buruk: benar-benar membutuhkan banyak waktu dan memori memerah. Tidak dapat membuat beberapa tabel sekaligus karena memori tidak bisa menanganinya. Dan hasilnya kotor dan berantakan. Alih-alih memberi saya nomor baris yang benar, dan saya perlu membersihkannya nanti (dijelaskan dengan baik di sini, Pisahkan garis pada titik persimpangan )
CREATE TABLE lines_with_messy_result AS (
SELECT
((ST_DUMP(ST_SPLIT(a.geom,b.ix))).geom) as geom
FROM st_union_lines a
INNER JOIN lots_of_points b
ON ST_INTERSECTS(a.geom, b.ix)
);
--then need to clean this up
create table lines_segments_cleaned as (
SELECT DISTINCT ON (ST_AsBinary(geom))
geom
FROM
lines_with_messy_result
);
sumber cara / pendekatan ini: /programming/25753348/how-do-i-divide-city-streets-by-intersection-using-postgis
2
Tabel garis / titik yang sama. Masih berantakan hasilnya dan perlu dibersihkan ini. Masih banyak waktu untuk menyelesaikan kueri.
--messy table
Create table messy_lines_segments as
Select
row_number() over() as id,
(st_dump(st_split(input.geom, blade.ix))).geom as geom
from st_union_lines input
inner join lots_of_points blade on st_intersects(input.geom, blade.ix);
--then cleaning the messy table
delete from messy_lines_segments a
where exists
(
select 1 from messy_lines_segments b where a.id != b.id and st_coveredby(b.geom,a.geom)
);
sumber cara / pendekatan ini: Membagi garis pada titik persimpangan
3
Saya juga menemukan fungsi ini: https://github.com/Remi-C/PPPP_utilities/blob/master/postgis/rc_Split_Line_By_Points.sql
yang memiliki hal baik yang tidak meninggalkan messy_result itu maka saya perlu membersihkannya. Tapi Anda perlu st_memunion dari kedua sisi (tabel garis dan tabel poin)
Ini semacam:
create table osm.calles_cortadas_03segmentos_sanluis as (
SELECT result.geom
FROM
osm.calles_cortadas_01uniones_sanluis AS line,
osm.calles_cortadas_00intersecciones_sanluis AS point,
rc_split_line_by_points(
input_line:=line.geom
,input_points:=point.ix
,tolerance:=4
) AS result
);
Tapi ini juga jam yang sangat lama untuk mendapatkan hasil. Dan saya juga mencoba dengan tabel yang lebih panjang (garis 10k, poin 14k) dan saya hanya mendapatkan masalah memori.
Saya juga mencoba Esri's ArcGIS dengan hasil yang juga buruk ...
Jadi, apa cara terbaik untuk menyelesaikan ini dengan fungsi geom PostGIS?
Maksudku, tanpa melangkah ke topologi.
Atau apa rekomendasi terbaik Anda?