Le portail francophone de la géomatique

Bonjour,

Pour ce genre de problématique (si je la comprends bien, ce qui n'est pas sûr), les fonctions de référencement linéaires marchent bien (https://postgis.net/docs/reference.html … eferencing).

Elles permettent de localiser un point par rapport à une ligne et de découper une ligne entre deux points de localisation (=locus).
Ces fonctions sont très (très) performantes par rapport à des intersections classiques, notamment entre point et ligne (pb de précision des points par rapport aux lignes).

Ici, pour découper des routes par rapport à des parcelles, on peut calculer les points d'intersection entre les routes et les contours des parcelles, puis calculer les locus de ces points d'intersection par rapport à chaque route, puis enfin découper les routes entre les locus consécutifs trouvés: le découpage est ultra rapide.

Les données: Tables route (routes bdtopo) et parc (parcelles cadastrales) sur Paris.
Les parcelles sont fusionnées et dumpées pour extraire des polygones simples (plus de frontieres entre parcelles contigües) => table parcfusion:

drop table if exists parcfusion;
CREATE TABLE parcfusion as (
  select (st_dump(st_union(geom))).path, (st_dump(st_union(geom))).geom
  from parc
);

alter table parcfusion add column id SERIAL PRIMARY KEY ;
create index parcfusion_geom_gist on parcfusion USING GIST (geom);
VACUUM ANALYSE parcfusion;

Ensuite, dans la CTE inter, on calcule l'intersection entre les routes et les contours externes des parcelles => on obtient des points ou des multipoints qu'on dump.

Dans la CTE locus, on calcule le locus (=localisation du point sur la ligne, entre 0 et 1) entre chaque route et les points d'intersection.
On en profite pour rajouter artificiellement les locus 0 et 1 (début et fin de ligne) pour chaque route. On a alors tout le temps, pour les routes, une sequence de valeurs entre 0 et 1 qui sont les points de découpage de la route.

Dans locus_id, on génère un id croissant continu sur les locus, apres avoir ordonné comme il faut dans l'etape précédente: par id de route, puis locus.

Dans route_cut, on fait un self-join sur la table locus_id histoire de pouvoir mettre sur une meme colonne une route, son locus n et le locus n+1: on découpe la route entre ses deux locus. On obtient alors des segments à l'intérieur ou à l'extérieur des parcelles. Dans une étape finale (pas faite ici), on peut retrouver les segments contenus dans les parcelles (le plus efficace pour cela est de comparer les longueurs d'intersection parcelle/troncon vs troncon).

Le processus tourne en qq secondes sur Paris: on fait une seule intersection entre les routes et les parcelles (apres fusion de celles-ci, qui peut prendre un peu de temps).

with inter as (
    SELECT
      r.gid,
      p.id,
      (st_dump(st_intersection(r.geom, st_boundary(p.geom)))).geom AS geom,
      r.geom as rgeom
    FROM route r
      JOIN parcfusion p ON st_intersects(r.geom, p.geom)
), locus as (
  select gid, st_linelocatepoint(rgeom, geom) as locus,
    rgeom
  from inter union all
    select gid, 0 as locus, rgeom
  from inter UNION ALL
    select gid, 1 as locus, rgeom
  from inter
  order by gid, locus
), locus_id as (
    SELECT
      row_number()
      OVER () AS id,
      *
    FROM locus
  where locus between 0.001 and 0.999
), route_cut as (
    SELECT l1.gid as gid1, st_linesubstring(l1.rgeom, l1.locus, l2.locus) as geom
    FROM locus_id l1
      JOIN locus_id l2 ON l1.gid = l2.gid AND l2.id = l1.id + 1
) select * from route_cut;

Quelques captures en attachement pour montrer le process: les parcelle initiales et les routes, puis les parcelles fusionnées, enfin les troncons de route a cheval sur des parcelles, découpés par les points d'intersection route/parcelle.

Nicolas

merci

Oui les versions recentes de PG et postgis permettent de faire des traitements de fou