3. Reprojeter des données SIG

3.1. Reprojeter des données raster

GDAL a la capacité de modifier un système de coordonnées raster en utilisant la syntaxe suivante:

gdalwarp -t_srs EPSG:... <input> <output>

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. <input> et <output> sont respectivement les données d'entrée et de sortie.

Nous allons maintenant reprojeter un modèle numérique d'élévation (DEM) acquis par la Shuttle Radar Topography Mission (SRTM). Vous pouvez utiliser ce site Web pour télécharger des DEM pour votre domaine d'intérêt: http://srtm.csi.cgiar.org/index.asp

Afin de reprojeter le DEM de WGS-84 lat / lon à Amersfoort / RD New, nous utilisons cette commande:

gdalwarp -t_srs EPSG:XXXXX srtm_37_02.tif dem_rd.tif

1. Remplacez XXXXX par le code EPSG approprié pour Amersfoort / RD New (voir l'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.

Plus d'informations sur la commande gdalwarppeuvent être trouvées ici.