TP partie 2. Systèmes de coordonnées

TP Utilisation en ligne de l'outil cs2cs

Dans un tableur, ouvrez le fichier Donnees_Climato_Dept34.xlsx dans le répertoire data/Occitanie/CD34. Ce fichier contient la position des stations météo départementales de l'Hérault, géoréférencées dans le système de coordonnées Lambert 93 (EPSG:2154).

Avec l'outil en ligne MyGeoData CS2CS (https://mygeodata.cloud/cs2cs/) , convertissez en WGS84 (EPSG:4326) les coordonnées des 10 premières lignes du fichier.

TP Identification du système de coordonnées

Ouvrez Quantum GIS 3.16 Desktop. Dans un nouveau projet, ajoutez les couches suivantes, dans cet ordre :

• Départements (DEPARTEMENT_CARTO.shp dans data/IGN_ETALAB/ADMIN-EXPRESS-COG)
• Lignes SNCF (formes-des-lignes-du-rfn.shp dans data/SNCF)

Identifiez le système de coordonnées du projet.

Identifiez le système de coordonnées de chaque couche.

Ajoutez la couche des gares géoréférencées à partir du fichier referentiel-gares-voyageurs.csv : vous devez utiliser les colonnes Longitude et Latitude en WGS84 (EPSG:4326).

Enregistrer cette couche de points dans un nouveau shapefile TP2/gares_l93.shp en Lambert 93 (clic-droit sur la couche de points + Exporter / Sauvegarder les entités sous).

TP Projection à la volée, calcul de longueur

Exemple avec un fichier GPX

Dans un nouveau projet, ouvrir le fichier monthaut2.gpx dans le répertoire GPS (waypoints, tracks et track_points). Les fichiers GPX sont en WGS84.

Superposer la trace issue du fichier GPS monthaut2.gpx et la carte topo de la zone concernée : dalle Scan25 de Saint-Guilhem-Le-Désert (répertoire : SCAN25, fichier SC25_TOPO_0740_6300_L93_E100.jp2 et SC25_TOPO_0740_6300_L93_E100.jp2).

Définir le SCR du projet : EPSG:2154 (Lambert 93).

Cliquez sur la trace GPS avec l'outil Identifier : dans (Dérivé), relevez la longueur sur l'ellipsoïde.

Avec l'outil "Ajouter des attributs de géométrie" (dans la boîte à outils de traitement : Géométrie vectorielle), mesurez la longueur totale de la trace en mètres avec le SCR du projet. Spécifiez comme nom de couche en sortie : ajouter_longueur_rando.shp

Dans la couche obtenue, relevez la valeur du champ length : il s'agit de la longueur de la trace projetée en Lambert 93.

Indicatrices de Tissot, Pseudo-Mercator et extension QuickMapServices

Dans un nouveau projet, ouvrez ces couches incluses dans le fichier GeoPackage data/NaturalEarth/NaturalEarth2020.gpkg :

• ne_50m_admin_0_countries
• ne_50m_coastline
• ne_50m_geographic_lines
• ne_50m_graticules_30

Utilisez l'outil Identifier sur la couche ne_50m_geographic_lines et la longueur indiquée dans (Dérivé) pour relever :

• la longueur de l'Equateur
• la longueur des Tropiques
• la longueur des Cercles Polaires

Superposer les indicatrices de Tissot et le quadrillage à 30° (fichier ne_50m_graticules_30.shp et tissot.shp.

Utiliser le mécanisme de projection à la volée pour afficher successivement la carte du monde avec les 3 projections suivantes :

1. World Robinson / ESPG:54030
2. World Winkel-Tripel / EPSG:54042
3. Pseudo-mercator / EPSG:3857

Laisser le projet en Pseudo-Mercator. Il s'agit d'un système de coordonnées couramment utilisé en webmapping ; avec l'extension "QuickMapServices" (cf TP1), ajoutez le fond de carte OpenStreetMap.

Avec l'outil "Ajouter des attributs de géométrie" (cf. TP précédent), calculer la surface des indicatrices le système Pseudo-Mercator (SCR du projet) pour quantifier la distortion de surface des cercles.

TP Transformation de données raster dans un nouveau système de coordonnées

Ouvrez une des dalles NASADEM du Liban (WGS84) et convertissez-la dans le système de coordonnées UTM 36N (EPSG:32636) avec une résolution de 30 mètres par pixel. Pour cela, utilisez l'outil "Projection(Warp)" accessible depuis la boîte à outils GDAL > Projections