Le portail francophone de la géomatique

st_perimeter ne fonctionne qu'avec des géométries de type surface.
Or, le résultat de l'intersection de votre requête ne sera jamais de type
surface,  étant donné que celle ci n'est effectuée que sur des géométries
ne partageant que leur limites, les géométries seront de type point ou ligne.

Si ce que vous voulez récupérer est effectivement la longueur des limites partagées
entre vos géométries, je pense que vous devriez remplacer st_perimeter par st_length.
(st_length renvoie 0 quand il est utilisé avec un (multi)point comme paramètre).

Dernière modification par tumasgiu (Mon 05 September 2016 19:12)

Bonjour michcarto,

tu peux aussi essayer ceci :

st_length(st_collection_extract(st_collect(ST_intersection(geom1987,geom)),2))/10000

On force la création d'une collection à partir de l'intersection, et on en extrait que les lignes (2).

Pour info st_length renvoie 0 pour les géométries de type (multi)point/polygones.
Pour les geometrycollection, elle est capable de retrouver les lignes
mêmes dans des collections imbriquées.
Exemple :

SELECT st_length(
              st_geomfromtext(
                     'GEOMETRYCOLLECTION(
                        POLYGON((0 0, 0 1, 1 1, 1 0, 0 0)), 
                        GEOMETRYCOLLECTION(
                           LINESTRING(0 0, 0 1), 
                           POLYGON((0 0, 0 1, 1 1, 1 0, 0 0))), 
                           MULTILINESTRING((0 0, 0 1), (0 0, 0 1)))'
                       )
                     ); 
--  trois lignes dans la collection, renvoie 3.0

En tout cas c'est le comportement de st_length sous Postgis.

Dernière modification par tumasgiu (Mon 12 September 2016 22:27)