Le portail francophone de la géomatique

Le fichier contexte.png montre la situation :
-  il y a 2 polylignes (TRC1 et TRC2)
- il y a une polyligne (1) qui se trouve à proximité de TRC1 et TRC2.

En fait, j'aimerais obtenir le résultat indiqué par result.png à savoir d'obtenir comme résultat le TRC1 et TRC2 découpé (par projection du point début/fin de la polyligne rouge).

Je ne sais pas si c'est très clair. Mais il s'agit de pouvoir mesurer l'impact de la polyligne (1) sur les TRC = la polyligne 1 concerne TRC1 sur x mètres et TRC2 sur y mètres.

Cordialement

j'ai projeté les startpoint/endpoint sur les lignes TRC.

Par contre, je me questionne sur la manière de découpe des lignes avec ces points projetés.

Si j'écris cela, j'ai un message d'erreur m'indiquant qu'il est impossible de "spliter" avec du multipoint

with point as (
    SELECT idcrit, ST_Union(geom) as geom FROM test4 GROUP BY idcrit

) SELECT ST_Split(point.geom, fea_troncon.geom) FROM fea_troncon JOIN point ON ST_DWithin(fea_troncon.geom, point.geom,1); 

Faut-il que je découpe point par point mais comment faire alors pour ne conserver que les segments situés entre les points start/end ?

Cordialement

Je n'arrive pas à rédiger la bonne syntaxe.

Voici ce que j'ai crée (je passe par plusieurs vues pour bien poser les étapes ... mais je pense que le code peut être largement épuré ....)

# Conversion des polylignes de découpe (fea_crit200) en points

CREATE OR REPLACE VIEW GeneratePoint AS
With tmp as
(
  SELECT idcrit, (ST_DumpPoints(geom)).geom as geom FROM fea_crit200)
SELECT idcrit, geom FROM tmp;

# Projection des points précédemment générés sur fea_troncon (polylignes à découper)

CREATE VIEW ProjectPoint as
with mpoint_by_line as( 
    SELECT idcrit, geom FROM GeneratePoint
) SELECT row_number() over() as num, mpoint_by_line.idcrit, ST_ClosestPoint(fea_troncon.geom, mpoint_by_line.geom) As geom FROM fea_troncon JOIN mpoint_by_line ON ST_DWithin(fea_troncon.geom, mpoint_by_line.geom,100);

# Sélection des segments concernés

CREATE OR REPLACE VIEW SplitFea AS
with point as (
    SELECT idcrit, geom FROM ProjectPoint

) SELECT ST_LineSubstring(fea_troncon.geom, 0, ST_LineLocatePoint(fea_troncon.geom, point.geom)) as geom FROM fea_troncon JOIN point ON ST_DWithin(fea_troncon.geom, point.geom,1);


Je n'y arrive pas dans le sens ou ma dernière vue traite point par point alors qu'il me faudrait le point start et le point end que je passerais en paramètre de ST_LineSubstring

Dernière modification par schummi (Wed 11 February 2015 18:05)