Tutoriel: Introduction à GDAL
3. Reprojection de données SIG
3.1. Reprojecter des données raster
GDAL permet de changer le système de coordonnées d'un raster en utilisant la syntaxe suivante:
gdalwarp
-t_srs EPSG:... <entrée> <sortie>
L'argument -t_srs spécifie
le système de coordonnées cible. Si le
système de coordonnées source est inconnu, il doit être spécifié avec
l'argument -s_srs . EPSG:...
spécifie
le code EPSG de la projection. <entrée> et <sortie> sont
respectivement les données d'entrée et de sortie.
Nous
allons maintenant reprojeter un Modèle Numérique d'Elévation (DEM) acquis par Shuttle
Radar Topography Mission (SRTM). Vous
pouvez télécharger des DEM pour votre propre zone d'intérêt à partir de USGS Earth Explorer. Ici,
nous utiliserons les données du cours.
Afin de reprojeter le MNE de WGS-84 lat/lon à Amersfoort/RD New nous utilisons cette commande:
gdalwarp -t_srs EPSG:XXXXX srtm_37_02.tif dem_rd.tif
XXXXX par le
code EPSG approprié pour Amersfoort/RD New (voir une de vos réponses
précédentes en utilisant ogrinfo).2. Exécutez la commande:
gdalwarp -t_srs EPSG:28992 srtm_37_02.tif dem_rd.tif <ENTER>
3. Visualisez le résultat dans QGIS.
Pour plus d'informations sur la commande gdalwarp, cliquez ici.