banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

#1 Wed 27 April 2022 11:22

sfp
Juste Inscrit !
Date d'inscription: 27 Apr 2022
Messages: 2

QGIS: Tampons et Systeme de Projections / Echelle mondiale

Bonjour,

J'ai une série de points GPS repartis un peu partout dans le monde. Je souhaite créer des tampons circulaires de 5 km de rayon autour de ces points. Cependant, en essayant plusieurs systèmes de projection (Lambert, WGS84, WGS84WMSV, Robinson, etc.) j'obtiens souvent des tampons déformés/ovales ou d'un rayon plus petit que celui voulu, notamment aux latitudes et longitudes fortes (par exemple Alaska ou Nouvelle Zélande).

Comment réussir à avoir des tampons circulaires à un rayon toujours constant, peu importe sa localisation dans le monde ? Quel système de projection utiliser ?

Merci d'avance pour votre aide smile

Hors ligne

 

#2 Wed 27 April 2022 16:30

Sylvain M.
Participant assidu
Lieu: Saint-Pierre-des-Nids (53)
Date d'inscription: 8 Sep 2005
Messages: 994

Re: QGIS: Tampons et Systeme de Projections / Echelle mondiale

Personnellement, j'ajouterais un champ "5Km_To_Degree" dans la table attributaire des points, dans lequel je calculerais la valeur de 5 km convertis en degrés, selon ce type de formule :

1 degree° = 60 arc minutes '
1 arc minute ' = 60 arc seconds ''
1 arc second '' at equatorial sea level = 1852m/60 = 30.86666667m

Formlua: cos(degree latitude) * (1852/60)


trouvée sur ce site : https://www.opendem.info/arc2meters.html
(je ne suis pas expert en trigonométrie, mais pour quelqu'un qui s'y connait, c'est la base wink )

(tu peux d'abord passer par un champ intermédiaire "latitude", que tu calcules avec $x sur le point)

Ensuite, tu peux faire un Tampon dont la distance sera contrôlée par ce champ "5Km_To_Degree".


Sylvain M.

Hors ligne

 

#3 Thu 28 April 2022 15:45

Sylvain M.
Participant assidu
Lieu: Saint-Pierre-des-Nids (53)
Date d'inscription: 8 Sep 2005
Messages: 994

Re: QGIS: Tampons et Systeme de Projections / Echelle mondiale

Pour le fun, j'ai tenté de trouver la formule (même si mes cours de trigo remontent au lycée...).
Voici donc une formule dans QGis pour calculer l'équivalent en degrés de 5km selon la latitude* :

Code:

(5000/(COS("lat"*PI()/180)*(1852/60)))/3600

(en ayant préalablement calculé un champ "lat" avec la latitude du point, qui peut par exemple se calculer avec $y par exemple si la couche est en WGS84)


Sylvain M.

Hors ligne

 

#4 Thu 28 April 2022 15:51

Sylvain M.
Participant assidu
Lieu: Saint-Pierre-des-Nids (53)
Date d'inscription: 8 Sep 2005
Messages: 994

Re: QGIS: Tampons et Systeme de Projections / Echelle mondiale

Sans passer par le champ "Lat", si la couche est en WGS84 :

Code:

(5000/(COS($y*PI()/180)*(1852/60)))/3600

Si elle est dans un autre EPSG (en remplaçant XXXX par le code EPSG) :

Code:

(5000/(COS(y(transform($geometry,'EPSG:XXXX','EPSG:4326'))*PI()/180)*(1852/60)))/3600

Sylvain M.

Hors ligne

 

#5 Thu 28 April 2022 16:10

Sylvain M.
Participant assidu
Lieu: Saint-Pierre-des-Nids (53)
Date d'inscription: 8 Sep 2005
Messages: 994

Re: QGIS: Tampons et Systeme de Projections / Echelle mondiale

Et ça peut donc se faire directement depuis l'outil Tampon, donc, sans passer par des champs intermédiaires, en définissant la distance de la zone tampon par la formule :
https://i.ibb.co/9TPxFkP/tampon.png


Sylvain M.

Hors ligne

 

Pied de page des forums

Powered by FluxBB