Le portail francophone de la géomatique

bonjour.

merci de ta réponse. je ne suis pas expert non plus sur postgis, mais on va s'en sortir
j'ai ajouté un polygone avec INSERT INTO polygons_department (id_po_de,   id_de,   polygone) VALUES (154,64,ST_GeomFromText( 'POLYGON((-0.278513523192217 42.9392355563593, ... ))',  -1 ) );

je crois que la projection correspond au second paramètre de ST_GeomFromText (pour moi -1), que faut-il mettre à la place ?
               
merci d'avance

bonjour,

Il faut mettre 4326 à la place de -1 dans ton cas.

st_area renvoie une surface en unité de la projection, donc il ne renverra pas de km²

Il y a peut être plus simple mais tu peux toujours reprojeter en lambert 2 étendue et calculer la surface.

st_area(transform(the_geom,27572))


Y.

Faut etre sur du SRID avant de mettre 4326.
Il existe des dizaines de SRS geographiques (lat/long).

D'ou viennent donc ces données ?

Bonjour,

Si le fichier suit bien le standard KML il doit être en 4326.

Y.

Bonsoir,

Mettre a jour le SRID des objets geometriques ne suffit pas. Ils sont toujours en lat/lon et postgis calcule toujours l'aire en degres carrés.
Pour calculer une aire en m2, il faut reprojeter les données dans un systeme de projection metrique, comme l'indiquait Yves.
Le systeme de projection choisi doit etre en accord avec la zone geographique des données. Pour la france, vous avez le choix entre Lambert 2 Etendu; SRID 27582, ou mieux, Lambert 93, srid 2154.

quelque chose comme:

select st_area(st_transform(polygone, 27582)) as aire from polygons_departement;

Nicolas

PS: pensez aussi a mettre la table geometry_columns a jour, apres avoir changé le SRID de la colonne polygone.
Cette table doit contenir une entrée décrivant la table polygons_departement.
Certains logiciels se servent de cette table pour lister les données geographiques presentes dans une base postgis.

hi hi! ca à l'air de fonctionner )
pour info, voici ce que j'ai passé pour remplir mes spatial_ref_sys :
psql toupil < /usr/share/postgresql-8.4/contrib/postgis-1.5/spatial_ref_sys.sql

merci à tous! je suis impressionné de ce qu'il est possible de faire avec postgis !!!

Hmm... /contrib/postgis-1.5/spatial_ref_sys.sql ?

Vous avez postgis 1.5 alors, finalement ? ou plusieurs installations ?. Si c'est le cas, il faut mieux installer les bons fichiers SQL sur la version de postgis courante.

Si vous avez 1.5, le nouveau type geography est a essayer, dans votre cas

Nico

Bonjour,

Petite question sur ce vieux post, en lien avec l'utilisation de la fonction st_area dans QGIS.

Après pas mal de recherches, l'utilisation de la fonction st_area, y compris avec l'outil transform, me renvoie toujours la même chose que l'utilisation de l'outil area($geometry) et donc, si j'ai bien compris, la surface "géométrique", sans prise en compte du système de projection et de l'ellipsoïde que je souhaiterais utiliser... ce qui me pose un souci de cohérence avec des surfaces calculées "en dur" en utilisant l'outil $area dans la table attributaire.

Un peu de contexte est sans doute nécessaire !
Jusque là, notre équipe calculait des surfaces "en dur" dans la table attributaire via la fonction $area, sur des entités contenues dans des couches projetées en Lambert 93/EPSG 2154, dans des projets paramétrés dans leurs propriétés pour calculer les surfaces en utilisant l'ellipsoïde GRS 80 (EPSG 7019) avec des m² comme unité.

Je cherche à calculer de manière dynamique ces surfaces grâce à une couche virtuelle.

J'ai essayé différentes syntaxes :
select st_area(COUCHE.geometry) as surf_sig from COUCHE

> Cela me retourne la même surface que l'outil area($geometry) en dur, assez logiquement si j'ai bien compris, la fonction renvoyant aux entités géométriques (et géographiques ?)

select st_area(COUCHE.geometry, true) as surf_sig from COUCHE

> Cela ne me retourne rien dans ma couche virtuelle

select st_area(st_transform(COUCHE.geometry, 2154)) as surf_sig_2154 from COUCHE

> Cela ne me retourne rien dans ma couche virtuelle

J'ai sans doute un problème de syntaxe et/ou de référence au système de référence spatiales mais que je n'arrive pas à identifier.

En vous remerciant pour votre aide ... Claire

Bonjour,
Peut-être que cela est à mettre en rapport avec la différence entre '$area' et 'area' ?
Cf. ici, je cite :