Le portail francophone de la géomatique

Bonjour,

Regardez du côté de la fonction ST_Intersection(). Cette dernière devrait vous donner le résultat attendu.

https://postgis.net/docs/ST_Intersection.html

Bonjour,

Effectivement st_intersection mais avec quelques clauses WHERE.

St_intersection() uniquement quand st_intersects(a,b) renvoie vrai. et st_intersects(a,b) renvoie vrai si a && b.

Bonjour,

Oui mais : intersection(intersection(a,b),c) = intersection (a,b,c).

Bonjour,

Perso, pour ce genre d'explosion totale, je me sers de st_polygonize et st_union sur les lignes (qui node les lignes):

with tmp as (
    SELECT st_union(st_boundary(geom)) as geom
    FROM testpg
), tmp1 as (
    SELECT st_dump(st_polygonize(geom)) as dmp
    FROM tmp
) select (dmp).path[1] as id, (dmp).geom
from tmp1;

Ca explose vite et bien, sauf pour les gros datasets, a traiter par parties.

Ensuite, des requetes spatiales de contenance, de distance, de % d'intersection permettent de retrouver les liens vers les objets initiaux (complètement perdus avec st_polygonize).

Nico

Une petite précision concerant st_polygonize vs st_buildArea:

st_polygonize() reconstruit tous les polygones possibles depuis des lignes qui s'intersectent.
st_buildArea() reconstruit toutes les surfaces possibles (au sens topologique): si 2 polygones se touchent, ils sont fusionnés par st_buildArea()

Et puis, pour répondre vraiment à la démande, même si je me doute que ce n'est pas ce que vous aviez en tête, une requête pour "exploser" les polygones, c'est à dire les décaler vers l'exterieur du centre de la scène:

WITH tmp AS (
    SELECT st_union(st_boundary(geom)) AS geom
    FROM testpg
), tmp1 AS (
    SELECT st_dump(st_polygonize(geom)) AS dmp
    FROM tmp
), tmp2 AS (
    SELECT
      (dmp).path [1] AS id,
      (dmp).geom
    FROM tmp1
), tmp3 AS (
    SELECT st_centroid(st_extent(geom) :: GEOMETRY) AS geom_center
    FROM tmp2
) SELECT
  t2.id,
  st_translate(
      geom,
      ((1.2 * st_X(st_centroid(geom))) - (st_x(geom_center) * (1.2 - 1))) - st_x(st_centroid(t2.geom)),
      ((1.2 * st_Y(st_centroid(geom))) - (st_Y(geom_center) * (1.2 - 1))) - st_y(st_centroid(t2.geom))
  )
FROM tmp2 t2, tmp3 t3;

Nicolas