Stavo controllando il modo migliore per tagliare i linestring per punti.

Lo scenario è: molte strade, hanno bisogno dei segmenti tagliati da punti di intersezione, questo è il seguente:

ho ottenuto

• tabella di stringhe di linea (interamente non suturata da punti)

• tabella dei punti st_intersection

Devo tagliare i segmenti di stringa indipendenti dalla tabella dei punti di intersezione.

Sto usando le funzioni di PostGIS e ho trovato diversi approcci, ma ognuno di loro mi sta dando qualche tipo di problema.

Ecco cosa ho già testato:

1

Tabella delle linee: 1 riga, st_memunion di 1200 linee Tabella dei punti: 1700 righe (punti)

Cosa c'è di male: richiede davvero molto tempo e memoria. Impossibile creare più tabelle contemporaneamente perché la memoria non è in grado di gestirla. E il risultato è sporco e disordinato. Invece di darmi il numero di riga corretto, e ho bisogno di ripulirlo in seguito (ben spiegato qui Dividi linee ai punti di intersezione )

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
);

fonte di questo modo / approccio: /programming/25753348/how-do-i-divide-city-streets-by-intersection-using-postgis

2

Stessa tabella di linee / punti. Risultati ancora disordinati e necessità di ripulire. Ancora molto tempo per terminare la query.

--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)
);

fonte di questo modo / approccio: linee di divisione in punti di intersezione

3

Ho trovato anche questa funzione: https://github.com/Remi-C/PPPP_utilities/blob/master/postgis/rc_Split_Line_By_Points.sql

che ha il bello che non lascia un messy_result che poi ho bisogno di ripulire. Ma hai bisogno di st_memunion da entrambi i lati (tabella delle linee e tabella dei punti)

È una specie di:

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
);

Ma sono anche ore lunghissime per ottenere i risultati. E ho anche provato con tabelle più lunghe (10k righe, 14k punti) e ho solo problemi di memoria.

Ho anche provato ArcGIS di Esri con risultati anche cattivi ...

Allora, qual è il modo migliore per farlo con le funzioni geom di PostGIS?

Voglio dire, senza entrare nella topologia.

O qual è la tua migliore raccomandazione?

Questo è il modo!

Ok, ho ricevuto un ottimo feedback da Remi-C e ora funziona come un incantesimo:

La migliore soluzione non topologica di sempre .. funziona DAVVERO in modo semplice e veloce (credetemi, ho provato molti modi per farlo):

--please add this function:
https://github.com/Remi-C/PPPP_utilities/blob/master/postgis/rc_split_multi.sql

-- please create a universal sequence for unique id multipurpose
CREATE SEQUENCE select_id
  INCREMENT 1
  MINVALUE 1
  MAXVALUE 9223372036854775807
  START 1
  CACHE 1
  CYCLE;



-- lines (non st_unioned linestrings, with gist index) : linetable
-- points (non st_unioned points, with gist index) : pointtable

-- first, lines with cuts (pointtable)
create table segment_lines_00 as (
WITH points_intersecting_lines AS( 
   SELECT lines.id AS lines_id, lines.geom AS line_geom,  ST_Collect(points.ix) AS blade
   FROM 
   linetable as lines, 
   pointtable as points
   WHERE st_dwithin(lines.geom, points.ix, 4) = true
   GROUP BY lines.id, lines.geom
)
SELECT lines_id, rc_Split_multi(line_geom, blade, 4)
FROM points_intersecting_lines
);
CREATE INDEX ON segment_lines_00 USING gist(rc_split_multi);


-- then, segments cutted by points
create table segment_lines_01 as (
select 
(ST_Dump(rc_split_multi)).geom as geom,  nextval('SELECT_id'::regclass) AS gid  from segment_lines_00
);
CREATE INDEX ON segment_lines_01 USING gist(geom);

Questo è tutto!.

Aggiungo anche, in aggiunta, il mio metodo, senza usare st_split ():

Per ogni linea controllo se ci sono alcuni punti di intersezione.

Se sì, il mio tentatore sarà simile a questo:

[line_id, fraction]
[1      , 0       ]
[1      , 0.3     ]
[1      , 1       ]
[2      , 0       ]
[2      , 0.88    ]
[2      , 1       ]
...

Una tabella che contiene l'ID linea e la frazione della lunghezza della linea in cui un punto interseca la linea.

Quindi unisco la frazione per coppia per creare le nuove linee tagliate

Recquirement:

• La tabella "punti" contiene gli endpoint di ciascuna riga.
• I punti intersecano perfettamente le linee.

La domanda:

DROP TABLE IF EXISTS temptable;

CREATE TABLE temptable as 
(
    SELECT ST_LineLocatePoint(a.geom,b.geom) as fraction,a.id, a.geom as shape
    FROM lines a, points b
    WHERE st_intersects(a.geom,b.geom)
    UNION
    SELECT 1 as fraction ,a.id, a.geom as shape -- It avoid to get an incomplete line when the original line is circular
    FROM lines a
    ORDER BY a.id, fraction -- important to keep the same order
);

ALTER TABLE temptable  ADD COLUMN gid SERIAL PRIMARY KEY;

DROP TABLE IF EXISTS LINE_SPLIT;
CREATE TABLE LINE_SPLIT AS
(
SELECT b.fend,a.*, ST_LineSubstring(a.shape,a.fstart,b.fend) as geom
FROM
    (
    SELECT *, fraction as fstart FROM temptable
    LIMIT (SELECT COUNT(*) FROM temptable)-1
    ) a,
    (
    SELECT fraction as fend,gid-1 as gid FROM temptable
    OFFSET 1
    ) b
WHERE a.gid = b.gid
AND a.fstart < b.fend
);

ALTER TABLE LINE_SPLIT DROP COLUMN IF EXISTS shape;
DROP TABLE IF EXISTS temptable;