Le portail francophone de la géomatique

Vous travaillez sur pgadmin ET postgis ?
Pourriez-vous decrire un peu plus precisement ce que vous souhaitez faire ? ces points existent ou representent des positions particulieres sur la ligne ?

Les fonctions de referencement lineaire remplissent cette tache, notamment st_line_interpolate_point():
http://www.postgis.org/documentation/ma … eferencing
Ainsi que les fonction d'intersection ou de difference.

Nicolas

Un hack rapide:

select astext(st_snaptogrid(st_difference(
    geomfromtext('LINESTRING(0 0, 2 2)', -1),
    st_expand(geomfromtext('POINT (1 1)', -1)::geometry, 0.00001)
), 0.1));

renvoie:

MULTILINESTRING((0 0,1 1),(1 1,2 2))

Decouper les lignes avec des points legerements etendus (petits rectangles), puis "resnapper" les geometries obtenues sur une grille plus grande pour supprimer les imprecisions.
faire un expand tres petit, en deca de la precisions des données, puis revenir a la precision des données (1 decimale apres la virgule, ici).

Sinon, plus robuste: utiliser st_line_interpolate_point pour calculer le point d'intersection entre une ligne et un point, puis reconstruire les lignes avec ce nouveau point d'intersection.

Nicolas

Quelque chose comme ca:

create table lignes (id serial primary key, geom geometry);
create table points (id serial primary key, geom geometry);

insert into points(geom) values (geomfromtext('POINT (0 0)'));
insert into points(geom) values (geomfromtext('POINT (6 0)'));
insert into points(geom) values (geomfromtext('POINT (3 3)'));
insert into lignes(geom) values (geomfromtext('MULTILINESTRING ((-1 0, 7 0), (3 0, 3 4))'));

select astext(snapToGrid(st_difference(lignes.geom, st_expand(points.geom, 0.001)), 0.1))
from lignes, points
where st_intersects(lignes.geom, points.geom);

create table resultat as
select lignes.id idlignes, points.id idpoints, snapToGrid(st_difference(lignes.geom, st_expand(points.geom, 0.001)), 0.1) geom
from lignes, points
where st_intersects(lignes.geom, points.geom);

avec une clause where un peu plus subtile pour eviter les segments redondants, ou un traitement apres le snaptogrid pour exploser les multilinestring en linestring et virer les doublons.

Nicolas

Bonjour à tous j'ai lu les postes précédents et je me suis renseigné sur le net sur le découpage d'une ligne par un point.

J'arrive bien à découper une ligne avec un point (voir requête jointe) mais je souhaiterais maintenant réaliser une fonction qui permettrait de découper l'ensemble des lignes au niveau des points. Je viens donc par ici pour chercher un peu d'aide sur la réalisation d'une fonction qui réitérerait le découpage en utilisant le gid de chaque vanne (dans la requête jointe je découpe au niveau du gid n°6).

Merci d'avance

DROP TABLE IF EXISTS coupe;
CREATE TABLE coupe(gid int4) with oids;
ALTER TABLE coupe ADD COLUMN the_geom geometry;

INSERT INTO coupe(the_geom)
/*INSERE LA PREMIERE PARTIE DE LA GEOMETRIE*/
(SELECT (ST_SnapToGrid(ST_Line_Substring(ln_geom, 0.0, location), 0))
/*CALCUL DE 0 A LOCATION %, ATTENTION A LA PRECISION*/
FROM
(
SELECT  (SELECT reseau.the_geom FROM reseau WHERE reseau.gid=(SELECT a.gid FROM
(SELECT reseau.gid AS gid, Distance((SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6 ),reseau.the_geom) AS distance FROM reseau
ORDER BY distance ASC limit 1) AS a )) AS ln_geom,
 ST_Line_Locate_Point((SELECT reseau.the_geom FROM reseau WHERE reseau.gid=
/*SELECTION DE LA POLYLIGNE LA PLUS PROCHE DU POINT*/
(SELECT a.gid FROM
(SELECT reseau.gid AS gid, Distance((SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6 ),reseau.the_geom) AS distance FROM reseau
ORDER BY distance ASC limit 1) AS a )),
 
(SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6))AS location
/*ST_LINE_LOCATE_POINT QUI CORRESPOND A LOCATION*/
   ) AS foo);
-------------------------------------------------------------------------------------------------------------------------------------------

INSERT INTO coupe(the_geom)
/*INSERE LA SECONDE PARTIE DE LA GEOMETRIE*/
(SELECT  (ST_SnapToGrid(ST_Line_Substring(ln_geom,location, 1.0), 0) ) 
/*CALCUL DE LOCATION % A 1 DE LA LIGNE, ATTENTION A LA PRECISION*/
FROM
(
SELECT  (SELECT reseau.the_geom FROM reseau WHERE reseau.gid=(SELECT a.gid FROM
(SELECT reseau.gid AS gid, Distance((SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6 ),reseau.the_geom) AS distance FROM reseau
ORDER BY distance ASC limit 1) AS a )) AS ln_geom,
 ST_Line_Locate_Point((SELECT reseau.the_geom FROM reseau WHERE reseau.gid=
/*JE DEBUTE LA SELECTION DE LA POLYLIGNE LA PLUS PROCHE DU POINT*/
 (SELECT a.gid FROM
(SELECT reseau.gid as gid, Distance((SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6 ),reseau.the_geom) AS distance FROM reseau
ORDER BY distance ASC LIMIT 1) AS a )), 

(SELECT vannes.the_geom FROM vannes WHERE vannes.gid =6))AS location
/*ST_LINE_LOCATE_POINT QUI CORRESPOND A LOCATION*/
   ) AS foo);