Pages: 1
- Sujet précédent - QGIS: Tampons et Systeme de Projections / Echelle mondiale - Sujet suivant
#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
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 )
(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
Pages: 1
- Sujet précédent - QGIS: Tampons et Systeme de Projections / Echelle mondiale - Sujet suivant