Suite

Système de coordonnées de référence avec unités en mètres pour le centre de l'Allemagne ?


J'ai trouvé plusieurs questions concernant la reprojection des données dans un CRS avec une unité de mètres et quel CRS utiliser. Habituellement, la contre-question est « pour quelle région, à quelle échelle » et là les discussions s'arrêtent. J'ai moi-même essayé de trouver un bon CRS, mais comme je ne suis pas un expert en la matière, toute cette affaire de CRS est extrêmement déroutante.

Donc, je cherche un bon CRS avec des unités en mètres pour : l'Allemagne centrale, une ville qui s'appelle Paderborn. La région d'intérêt est assez petite, ca. 15km autour de Paderborn.

Je crée des simulations de bruit pour les parcs éoliens. J'ai donc des données ponctuelles (actuellement dans EPSG:3857) pour les éoliennes. Ensuite, dans un script python, je crée une grille équidistante autour de ces points, calcule les niveaux de bruit sur ces points de grille, utilise des méthodes de contour pour créer des "isophones" (lignes avec un niveau de bruit identique, par exemple 45 dBA) et enregistre ces données sous forme de ligne ou des données de polygone.

En général, tout cela fonctionne très bien et donne de beaux résultats.

Mais la partie vraiment sale est la suivante : dans mon script, je prends les données EPSG:3857 et les transforme en EPSG:4326. Ensuite, j'utilise une transformation haversine auto-codée pour utiliser les données déjà transformées dans "ma grille" avec une origine que je choisis. J'exécute la simulation, je fais l'analyse, j'utilise à nouveau une transformée de haversine inverse autocodée en EPSG:4326 et la transformée à nouveau en EPSG:3857.

Ce n'est probablement pas une bonne idée… Des recommandations ?


je suggère d'utiliserEPSG : 25832 ETRS89 / zone UTM 32N, car les zones Gauss-Krueger DHDN CRS sont basées sur l'ancien ellipsoïde de Bessel, et vous devez appliquer un décalage de référence qui n'est pas constant sur l'ensemble de l'Allemagne.

ETRS89 n'a pas de décalage de référence versEPSG : 4326 WGS84que vous pouvez obtenir à partir de récepteurs GPS, et toutes les informations cadastrales officielles ont déjà été transférées vers ETRS89.

Vous pouvez également configurer une projection de Mercator transversale locale sur un point central de votre zone d'intérêt, mais cela n'améliorerait que peu la précision par rapport à l'UTM.

D'ailleurs,EPSG : 3857n'a pas de compteurs réels comme unités, ils ne correspondent qu'à l'équateur.


J'utiliserais le CRS officiel en Allemagne : http://crs.bkg.bund.de/crseu/crs/eu-description.php?crs_id=Y0RFX0RIRE4gLyBHS18z

Pour la postérité, des projections de Gauss-Krueger comme 31466, 31467, 31468 ou 31469 (une par zone). Votre ville appartient à la zone 2 (ESPG:31467) : http://georepository.com/crs_31467/DHDN-3-degree-Gauss-Kruger-zone-3.html