Le portail francophone de la géomatique

JeanMich · Mon 27 June 2022 17:05

Bonjour,

Je cherche à calculer des surfaces non planes dans des virtual de layer.

Si je compare cette requète:

select LIBELLE,
SUM(st_area(ZONE_URBA.geometry)) as 'surf'
from ZONE_URBA

aux résultats obtenus avec la calculatrice de champs avec $area, j'ai des différences de surfaces
car je ne sais comment prendre en compte le src dans cette requête.

Dans postgis on pourrait écrire st_area(geom::geography)
mais je ne connais pas la correspondance avec le Sqlite utilisé dans les virtual layer.

J'espère être assez clair, merci de vos conseils.

Dernière modification par JeanMich (Mon 27 June 2022 17:07)

JeanMich · Tue 19 July 2022 16:06

.

Dernière modification par JeanMich (Tue 19 July 2022 16:11)

JeanMich · Tue 19 July 2022 16:08

J'ai essayé

select sum(st_area(ZONE_URBA.geometry, true)) as surf
from ZONE_URBA

mais sans succès.

D'après la doc: http://www.gaia-gis.it/gaia-sins/spatia … .0.html#p4

pour utiliser st_area avec use_ellipsoid il faut le module LWGEOM. C'est peut être ca le problème ?

Dernière modification par JeanMich (Tue 19 July 2022 16:12)

JeanMich · Tue 19 July 2022 16:08

.

Dernière modification par JeanMich (Tue 19 July 2022 16:11)

JeanMich · Tue 13 September 2022 09:55

Bonjour,

Je cherche toujours la réponse. Je ne sais pas si la question est pas clair, ou si il n'y a pas de réponse...

Mathieu Denat · Wed 14 September 2022 16:19

Bonjour,
Sans grande conviction, avez vous essayé:

Code:

select LIBELLE,
SUM(st_area(geometry)) as 'surf'
from ZONE_URBA

Bonne continuation,

Sylvain M. · Wed 14 September 2022 16:35

Ne manque-t-il pas tout simplement un "GROUP BY" ?

Code:

select LIBELLE,
SUM(st_area(ZONE_URBA.geometry)) as 'surf'
from ZONE_URBA
GROUP BY LIBELLE

Sylvain M. · Wed 14 September 2022 16:38

Sinon, si vous travaillez en France métropolitaine, vous pouvez reprojeter vos géométries en Lambert 93 avec cette formule : ST_Transform(ZONE_URBA.geometry, 2154)

Ce qui donnerait donc :

Code:

select LIBELLE,
SUM(st_area(ST_Transform(ZONE_URBA.geometry, 2154))) as 'surf_m2'
from ZONE_URBA
GROUP BY LIBELLE

Miniopterine · Wed 14 September 2022 16:43

JeanMich a écrit:

J'ai essayé

select sum(st_area(ZONE_URBA.geometry, true)) as surf
from ZONE_URBA

mais sans succès.

D'après la doc: http://www.gaia-gis.it/gaia-sins/spatia … .0.html#p4

pour utiliser st_area avec use_ellipsoid il faut le module LWGEOM. C'est peut être ca le problème ?

Ca ressemble à ce pb: https://github.com/qgis/QGIS/issues/41890

JeanMich · Fri 16 September 2022 11:21

Mathieu Denat a écrit:

Bonjour,
Sans grande conviction, avez vous essayé:
Code:
select LIBELLE,
SUM(st_area(geometry)) as 'surf'
from ZONE_URBA
Bonne continuation,

Oui, ca donne la surface cartésienne, donc fausse. Mais merci !

JeanMich · Fri 16 September 2022 11:25

Sylvain M. a écrit:

Sinon, si vous travaillez en France métropolitaine, vous pouvez reprojeter vos géométries en Lambert 93 avec cette formule : ST_Transform(ZONE_URBA.geometry, 2154)

Ce qui donnerait donc :
Code:
select LIBELLE,
SUM(st_area(ST_Transform(ZONE_URBA.geometry, 2154))) as 'surf_m2'
from ZONE_URBA
GROUP BY LIBELLE

J'ai aussi essayé, mais là on reprojette qlq chose de faux, le résultat n'est pas bon. Merci d'avoir essayé.

JeanMich · Fri 16 September 2022 11:33

Miniopterine a écrit:

JeanMich a écrit:
J'ai essayé

select sum(st_area(ZONE_URBA.geometry, true)) as surf
from ZONE_URBA

mais sans succès.

D'après la doc: http://www.gaia-gis.it/gaia-sins/spatia … .0.html#p4

pour utiliser st_area avec use_ellipsoid il faut le module LWGEOM. C'est peut être ca le problème ?

Ca ressemble à ce pb: https://github.com/qgis/QGIS/issues/41890

j'ai essayé avec:

select LIBELLE,
st_area(geometry),
st_area(geometry,true),
st_area(geometry,false)
from ZONE_URBA
GROUP BY LIBELLE

st_area(geometry) = renvoi un résultat faux
st_area(geometry,true), st_area(geometry,false) = renvoi un NULL que je passe par les virtual layer ou le db manager.

Je ne suis pas sûr de tout comprendre dans les réponses de github. Y a t il une solution ?

Merci à tous ceux qui ont essayé

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Mon 27 June 2022 17:05

QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#2 Tue 19 July 2022 16:06

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#3 Tue 19 July 2022 16:08

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#4 Tue 19 July 2022 16:08

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#5 Tue 13 September 2022 09:55

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#6 Wed 14 September 2022 16:19

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Code:

#7 Wed 14 September 2022 16:35

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Code:

#8 Wed 14 September 2022 16:38

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Code:

#9 Wed 14 September 2022 16:43

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

#10 Fri 16 September 2022 11:21

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Code:

#11 Fri 16 September 2022 11:25

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Code:

#12 Fri 16 September 2022 11:33

Re: QGIS 3.16: Calcul de surface non cartesienne avec les Virtual Layer

Pied de page des forums