Le portail francophone de la géomatique

Pierre Louis · Mon 27 January 2020 21:42

Bonjour,

Je m'adresse à vous pour vous demandez confirmation : la valeur de la surface calculée avec la projection mondiale "Pseudo-Mercator" (epsg 3857) est-elle acceptable pour des précisions au kilomètre près ?
En effet travaillant sur plusieurs polygones de tout le globe je me suis enquiquiné à calculer chaque polygone selon la projection locale et je vois en comparant avec le 3857 que les différences n’excèdent pas le km².
Cependant quand j'utilise l'extension "calculate geometry" les surfaces sont complètement différentes et donne la valeur selon les caractéristiques de la projection (Suède = 2 163 011 km² au lieu de 446 070 km² - sources : gadm). Je pensais donc que la fonction "$area" convertissait automatiquement au système "local" selon la position x y du polygone.

merci bien,

A bientôt

Pierre Louis

dominique.lys · Tue 28 January 2020 07:26

Les distorsion dépasseront facilement le km, elles dépendent de la latitude, on peut se faire une idée de la valeur réelle en multipliant une distance mesurée par cos(lat). Voir aussi les indicatrices de Tissot : https://en.wikipedia.org/wiki/Mercator_ … Distortion

La fonction $area s'appuie désormais sur l'ellipsoide pour éviter les surprises, cf. la doc de la fonction :

Returns the area of the current feature. The area calculated by this function respects both the current project's ellipsoid setting and area unit settings. For example, if an ellipsoid has been set for the project then the calculated area will be ellipsoidal, and if no ellipsoid is set then the calculated area will be planimetric.

ChristopheV · Tue 28 January 2020 08:39

Bonjour,

Pierre Louis faudrait commencer par travailler en CC42 localement

T'as pris de mauvaises habitudes Cours Général Leclerc ?

(pas pu m'empêcher)

Pierre Louis · Tue 28 January 2020 20:46

Merci beaucoup pour les indications et informations.
Je m'excuse à l'avance si je vulgarise jusqu'aux frontières de l'erreur en parlant des projections, je ne suis pas spécialiste.

Je poursuis donc le job avec le $area sans avoir le souci alors.
Si je comprends bien "calculate geometry" repose donc sur la projection directement, sur le plan quoi. Que je sois donc en RGF93 ou WGS84 j'aurai le même résultat étant données que les systèmes se basent sur l’ellipsoïde IAG-GRS80 ?

Cordialement,

ps : CC 42 c quoi ça encore ??
A bientôt j'espère sur le beau cours général Leclerc, haut lieu symbolique de bien des choses !

ChristopheV · Wed 29 January 2020 01:34

Salut PL, bonsoir à tous

Avec geometry tu travailles sur le globe en coordonnées lat long et les surfaces ne sont pas planes mais courbes.
Avec st_area() classique tu travailles sur le plan (à plat ) donc là tu fais de la géométrie classique avec des polygones 2D.

Or faire une carte (un truc plan) qui représente l'ensemble d'une sphère forcément ça déforme un peu non ? Donc pour ce faire suivant les besoins ie privilégier les zones que l'ont veut voir le moins déformées, on a inventer des méthodes pour passer des coordonnées de la sphère (coordonnées géographiques) aux coordonnées planes. Et ça se nomme une projection.

Ensuite la sphère c'en est pas une c'est un œuf, un ellipsoïde. Pareil il en existe plein.

Donc si tu veux comparer des surfaces de façon homogène sur l'ensemble du globe il te faut indiquer l'ellipsoïde choisi comme cité par dominique.lys et faire le calcul en coordonnées géographiques et pas planimétriques.

Et WGS84 ou RGF93 ça le fait bien.

Dernière modification par ChristopheV (Wed 29 January 2020 01:36)

Jean-Michel · Wed 29 January 2020 14:29

Bonjour,

Pour CC42, Christophe y tient beaucoup !
Pour le quoitesesquecé, voir par ici notamment : https://www.ideobfc.fr/upload/gedit/1/D … _zones.pdf
Géomatiquement

Pierre Louis · Wed 29 January 2020 14:52

slt Christophe et toutes et tous,

Oui j'ai suivi ce que dominique.lys disait à propos de la fonction $area dans Qgis qui calcule directement sur l'ellipsoïde.
Je vais donc transformer toutes mes données en epsg 4326 et arrêter avec mon 3857don-bi..
Je te rassure j'utilise autre chose lorsque je suis en Corse

Mais c'est le plugin "calculate geometry" qui m'a aussi embrouillé... Je vais lire la notice pour apprendre plus sur son fonctionnement..
Me suis arraché les cheveux...

merci bien en tout cas je suis éclairé et rassuré.

A + tout le monde !

Pierre Louis · Wed 29 January 2020 16:15

Bonjour Jean-Michel,
C'est que christophe me taquinait sur le CC42 par rapport au L93 etc.

merci en tout cas.

a +

Bruno · Thu 30 January 2020 08:36

Bonjour les amis corses et bretons!
Et finalement ce sujet a-t-il à voir avec QGIS?

Bonne journée!

fbecir · Tue 18 August 2020 18:32

Bonjour

Il y a deux choses à considérer :
1. Le mode de calcul de QGIS. Que fait exactement QGIS lors de son calcul ? Corrige-t-il automatiquement l'altération linéaire de la projection ?
2. La projection à considérer. Une projection que vous qualifiez de "locale" ne vous garantit pas un calcul juste des surfaces. En effet, la plupart des projections utilisées sont conformes, donc préservent les angles mais pas les surfaces. Il vous faut une projection équivalente qui conserve les surfaces (voir https://fr.wikipedia.org/wiki/Projectio … raphique).

Si vous travaillez en Europe, il y a justement une projection équivalente qui a été définie dans le cadre d'INSPIRE : ETRS89-LAEA. En utilisant une telle projection, vous vous évitez les questions existentielles sur le mode de calcul de votre SIG (que ce soit QGIS ou un autre).

Cordialement

Patrice · Wed 19 August 2020 10:27

Hello

Par exemple dans AutoCAD MAP ou AutoCAD CIVIL, sur la commande clavier : MAPCSASSIGN < Entree >
( qui permet d assigner un System de Projection a un DWG ) ,
cette projection speciale se nomme : ETRF89.Europe/EqArea

Je dois avouer que je ne l'ai jamais utilisee ! Merci pour l Info !!

LA SANTE, GeoBye, Pat

SANTANNA · Wed 19 August 2020 12:48

Bonjour,
Euh... Le message de OwenJoseph est une copie du premier message de ce fil, de Pierre Louis. "Spammeur"?

Patrice · Wed 19 August 2020 13:04

Hello Mr Santanna

Bien vu ! Je m en suis occupe ...

LA SANTE (Stay Safe), GeoBye, Pat

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 January 2020 21:42

projection et valeurs surfaciques avec Qgis

#2 Tue 28 January 2020 07:26

Re: projection et valeurs surfaciques avec Qgis

#3 Tue 28 January 2020 08:39

Re: projection et valeurs surfaciques avec Qgis

#4 Tue 28 January 2020 20:46

Re: projection et valeurs surfaciques avec Qgis

#5 Wed 29 January 2020 01:34

Re: projection et valeurs surfaciques avec Qgis

#6 Wed 29 January 2020 14:29

Re: projection et valeurs surfaciques avec Qgis

#7 Wed 29 January 2020 14:52

Re: projection et valeurs surfaciques avec Qgis

#8 Wed 29 January 2020 16:15

Re: projection et valeurs surfaciques avec Qgis

#9 Thu 30 January 2020 08:36

Re: projection et valeurs surfaciques avec Qgis

#10 Tue 18 August 2020 18:32

Re: projection et valeurs surfaciques avec Qgis

#11 Wed 19 August 2020 10:27

Re: projection et valeurs surfaciques avec Qgis

#12 Wed 19 August 2020 12:48

Re: projection et valeurs surfaciques avec Qgis

#13 Wed 19 August 2020 13:04

Re: projection et valeurs surfaciques avec Qgis

Pied de page des forums