Le portail francophone de la géomatique

l'image en PJ que je n'ai pas réussi à attacher au précédent message désolé

EDIT : bon j'arrive toujours pas à attacher l'image... tant pis!

Dernière modification par zang (Wed 10 April 2019 13:07)

Merci Jean-Marie,
J'ai juste ajouté une parenthèse pour fermer le ST_Intersection.
Par contre j'ai le message suivant:
[Err] ERREUR:  plus d'une ligne renvoyée par une sous-requête utilisée comme une expression
C'est peut-être pas la meilleure manière d'opérer pour générer ces lignes frontières, si?
Merci!
Franck

Bonjour Franck,
Il faut rajouter une clause WHERE dans la sous-requete afin que ST_intersection ne "reçoive" qu'une seule geometrie

En isolant la segmentisation d'encaissant_v3
dans une CTE ou une sous requête,
je pense que ca devrait passer :

with segments as (
    SELECT
        DISTINCT
                St_makeline(
                  ST_PointN(geom, generate_series(1, ST_NPoints(geom)-1)),
                  ST_PointN(geom, generate_series(2, ST_NPoints(geom) ))
                    ) line
    FROM
        (SELECT 
            (ST_Dump(ST_Boundary(geom))).geom
        FROM 
            encaissant_v3
        ) AS linestrings 
    )

INSERT INTO matable
SELECT 
    DISTINCT inter 
FROM 
    (SELECT 
        st_intersection(geom, line) inter 
    FROM 
        segments
    JOIN
        matable ON geom && line
    ) foo 
WHERE
    not st_isempty(inter)

Dernière modification par tumasgiu (Thu 11 April 2019 10:04)

Quelle version de PostgreSQL utilisez vous ?

j'ai l'impression que c'est

                St_makeline(
                  ST_PointN(geom, generate_series(1, ST_NPoints(geom)-1)),
                  ST_PointN(geom, generate_series(2, ST_NPoints(geom) ))
                    ) line

qui coince.

La requête est tres simple :
le bloc WITH est une autre manière d'écrire des sous-requêtes,
en les écrivant en préambule.
Pour la requête en elle même, plutot que d'utiliser une sous requête
pour générer le second argument de st_intersection,
on croise une sous requête qui est l'ensemble des segments,
avec l'ensemble des polygones, on génère toutes les intersections probables,
on ecarte les intersections nulles et les doublons.

Il faudrait que vous utilisez Postgres >= 10 (je crois)
pour que cette syntaxe fonctionne.
Ce qui coince ce sont les deux appels à generate_series
comme arguments d'une même fonction.
Ceux ci sont censés générer tout les couples d'indices
(Exemple pour un polygone à 4 coté, donc composé de 5 points: (0, 1), (1, 2), (2, 3) , (3, 4) )
pour récupérer les points grâce à la fonction st_pointn.

En déplaçant la création de cet indice dans une requête LATERAL (ou corrélée),
çà devrait fonctionner :

with segments as (
    SELECT
        DISTINCT
                St_makeline(
                  ST_PointN(geom, sp),
                  ST_PointN(geom, sp + 1)
                    ) line
    FROM
        (SELECT 
            (ST_Dump(ST_Boundary(geom))).geom
        FROM 
            encaissant_v3
        ) AS linestrings,
    LATERAL generate_series(1, st_npoint(geom) - 1) sp
    )

INSERT INTO talus (geom)
SELECT 
    DISTINCT inter 
FROM 
    (SELECT 
        st_intersection(geom, line) inter 
    FROM 
        segments
    JOIN
        talus ON geom && line
    ) foo 
WHERE
    not st_isempty(inter)

Si vous voulez faire la MAJ, c'est une très bonne chose
(dans ce cas là, essayez directement la 11 )
mais la requête précédente devrait fonctionner de manière équivalente.