Le portail francophone de la géomatique

Bonsoir,

Le résultat de unnest_st_clusterIntersecting(... renvoie déja une géométrie: vous avez une geometryCollection contenant tous les objets qui s'intersectent entre eux.

Si vous etes interessés par la partie surfacique de ce résultat, il faut utiliser st_collectionExtract(geom, 3) pour renvoyer les polygones et multipolygones contenus dans la geometryCollection:
(https://postgis.net/docs/ST_CollectionExtract.html)

SELECT ST_CollectionExtract(unnest(st_clusterIntersecting(geom))), 3) FROM ma_table

Nicolas

de peu

Bonjour, après plusieurs tests, j'ai explosé le st_clusterIntersecting avec un st_dump. Le résultat me paraît incertain. En effet le résultat ne me donne pas les géométries qui s'intersectent. Pour 2 parcelles qui se croisent dans ma couche, je devrais avoir un polygone se limitant au croisement. Le résultat me donne un polygone englobant l'intersection + parcelle et l'autre la parcelle sans l'intersection.

Cordialement

Bonjour,

Pour le suivi, j'ai réussi à trouver une solution à mon problème avec cette requête ;

        -- gestion des intersections
    WITH
    req_inter AS
    (
    SELECT
    row_number() over() AS gid, 
        -- je garde uniquement ici les polygones d'intersections (avec un buffer pour éliminer les artéfacts) 
        st_buffer(st_union(ST_CollectionExtract(St_intersection(p1.geom,p2.geom),3)),0.5) as geom
    FROM
     [ma_table_parcelle] p1, [ma_table_parcelle]  p2
    WHERE
        - - condition que le résultat du croisement soit une géométrie
    ST_Overlaps(p1.geom,p2.geom) is true
    AND p1.insee = p2.insee AND p1.insee = 'nnnnn'
    GROUP BY p1.insee, p2.insee
    )
    SELECT
    row_number() over() as gid, 
        p.*
        -- je supprime de ma table parcelle les intersections calculées ci-avant
    st_difference(p.geom,i.geom) as geom
    FROM
    req_inter i, [ma_table_parcelle]  p
    WHERE p.insee = 'nnnnn'

Ce résultat de parcelles sans intersection, va me permettre de déterminer les zonages de notre PLU sans remonter des erreurs dues au chevauchement présent dans les données de la DGFIP (en attendant le RPCU !)

Cordialement