Le portail francophone de la géomatique

sfp · Wed 27 April 2022 11:22

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

Sylvain M. · Wed 27 April 2022 16:30

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. · Thu 28 April 2022 15:45

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. · Thu 28 April 2022 15:51

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. · Thu 28 April 2022 16:10

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

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 Wed 27 April 2022 11:22

QGIS: Tampons et Systeme de Projections / Echelle mondiale

#2 Wed 27 April 2022 16:30

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

#3 Thu 28 April 2022 15:45

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

Code:

#4 Thu 28 April 2022 15:51

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

Code:

Code:

#5 Thu 28 April 2022 16:10

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

Pied de page des forums