RAPPORTS
Direction rgionale
de lenvironnement, de
l amnagement et du
logement
Rhne-Alpes
CEPE
UDG
Fvrier 2013
Traitements de MNT avec QGIS :
extraction, dcoupage,
assemblage, rchantilonnage
Ministre de lcologie,
du Dveloppement durable des et de l''Energie
http://www.rhone-alpes.developpement-durable.gouv.fr/
Affaire suivie par
Yann CONVERS CEPE/UDG
Tl. 03 26 28 67 89
[email protected]
Version 0.4 du 05/03/2013
Relecture assurance qualit [email protected]
Sommaire
1PRAMBULE........................................................................................................... 4
2DELIMITER LEMPRISE DE LA ZONE................................................................... 4
3LECTURE DES CARACTRISTIQUES D'UN MNT................................................5
4DCOUPE D'UN RASTER DPARTEMENTAL SUR UNE EMPRISE .................6
5RCHANTILLONAGE DU RASTER FINAL ..........................................................7
6ASSEMBLAGE DE PLUSIEURS RASTERS ..........................................................10
7CONVERTIR LE RSULTAT EN FICHIER ASC .....................................................11
8UTILISATION DE GDAL SANS QGIS..................................................................... 13
2/12
1 Prambule
Les manipulations dcrites ci-dessous sappuie sur la bibliothque GDAL-OGR et peuvent tre
excutes en lignes de commandes directement. Toutefois certaines captures dcran
sappuient sur limplmentation de GDAL-OGR dans QGIS version 1.8
2 Delimiter lemprise de la zone
Dans QGIS ouvrir la table contenant les emprises ainsi que les dalles reprsentant le MNT des
dpartements traiter :
Pour mmoire les dalles dpartementales des MNT (.ASC) se situent dans
X:\REF_EXT\BDTOPO\MNE
Vous devez obtenir un rsultat proche de lillustration ci-dessous :
3/12
3 Lecture des caractristiques d'un MNT.
Pour connatre le nombre de lignes et de colonne, le pas (rsolution) et les coordonnes du
coin bas gauche du MNT source ditez-le avec un diteur de texte, vous obtiendrez les valeurs
suivantes :
Lecture dtaille des caractristiques d'un fichier ASC :
ncolls, nrows : nombre respectivement de colonne et de ligne du fichier dpartemental.
Xllcorner,ullcorner : coordonnes X et Y (en loccurrence en lambert93) du coin bas gauche du premier
pixel gnr. (LL = Lower Left?)
Cellsize : rsolution du fichier 1 pixel reprsente un carr de 25 sur 25 m en ralit.
NODATA_value : l'image tant rectangulaire, une partie ce celle-ci peut ne pas contenir d'altitude (pays
frontalier par exemple), auquel cas pour viter la confusion avec le 0 altimtrique, on fixe celle-ci
-9999.
4/12
4 Dcoupe d'un raster dpartemental sur une emprise
partir du menu raster\extraction\dcouper :
1. Slectionnez le MNT sur lequel vous souhaiter faire la dcoupe (fichier source).
2. Indiquez le nom et l'emplacement du fichier final (fichier en sortie) et son format on prfrera les
fichiers ASC si possible.
3. Cochez la case Valeur nulle et saisissez la valeur NODATA contenu dans le fichier source (en
loccurrence ici -9999)
4. Enfin cliquez sur emprise et dessinez un rectangle d'encombrement autour de votre zone d'tude
(cf ci-dessus)
Les coordonnes sont directement extraites dans la bote de dialogue.
Dans l'exemple ci-dessus le rectangle est cheval sur les dpartements 38 et 42 et on ne peut pas
extraire un raster en dehors de son emprise. On va donc limiter l'emprise au bord stricte du raster du
dpartement 38.
Le plus simple consiste lire la coordonne X du bord gauche du raster dans l'entte du fichier
ASC en loccurrence 834987.50 et la substituer dans la boite de dialogue.
5/12
Attention en cas de dcoupe d'autres dpartements il faut penser sauvegarder les
coordonnes prcdentes, car si l'emprise diffre en Y on aura une bande de valeur nulle lors de
l'assemblage des rasters.
Sans autre action de notre part gdal ne tronquera pas les pixels en dehors de l'emprise dfinie et
fournira un nombre entier de ligne colonne.
En ligne de commande, la commande taper est
gdal_translate -a_nodata -9999 -projwin X1 Y1 X2 Y2 -og AAIGrid chemin_depart chemin arrive
5 Rchantillonage du raster final
Si vous disposez sur le territoire d'tude d'un MNT de meilleure prcision (pas) vous pouvez
rchantillonner le MNT en provenance de la Bdtopo, sachant effectivement qu'il ne s'agira que d'une
interpolation de valeur existante... la prcision finale n'en sera donc pas meilleure.
6/12
Le rchantillonnage ne fonctionne qu'avec des TIF aussi ?
Dans lexemple ci-dessous on traitera le cas dun MNT disponible au pas de 5m.
1. Utiliser la commande du menu Raster\Projection\projection (gdalwarp)
1. Spcifier en fichier source lassemblage ralis au point prcdent, puis en sortie lemplacement
et le nom du fichier rchantillonn.
2. Prcisez le systme de projection en entre et en sortie
3. Choisissez la mthode de rchantillonnage, ici bilinaire.
Enfin pour saisir le nouveau pas il faut diter la ligne de commande, le plugin Qgis noffrant pas
nativement cette possibilit :
7/12
1. Cliquez sur le stylo et rajouter l'option -tr 5 5 . Ou la premire valeur reprsente la largeur et
la seconde la hauteur du pixel dans lunit de mesure initiale du raster.
Il n'est pas ncessaire de spcifier un src source et destination si on ne change pas la
projection
En ligne de commande, la commande taper est
gdalwarp -r bilinear -tr 5 5 dstnodata -9999 -of Gtiff fichier_depart.tif fichier_arrivee.tif
8/12
6 Assemblage de plusieurs rasters
1. Ouvrez les deux rasters dpartementaux
2. Puis utiliser la commander Raster\Divers\fusionner
1. Spcifier en sources les deux fichiers rasters crs prcdemment en cliquant sur les noms des
fichiers avec la touche CTRL pour en slectionner plusieurs
2. ne pas oublier de spcifier -9999 en valeur nulle
3.
4. En sortie le format de fichier ASC n'est pas disponible pour la commande gdal_merge,
choisissez le format Gtiff pour crer l'assemblage.
9/12
Vous obtiendrez alors un GTIFF sur l'ensemble du territoire d'tude au pas de 25 m (ou diffrent selon
votre source de donnes)
En ligne de commande, la commande taper est
gdal_merge.bat -n -9999 -of Gtiff -o fichier_destination.tif fichier_depart1.asc fichier_depart2.asc
7 Convertir le rsultat en fichier ASC
1. Utiliser la commande Raster/conversion/Convertir
2. Spcifier le fichier source ;
3. puis le fichier destination en pensant changer le type de fichier dans la boite de dialogue (ASC)
4. Ajouter la valeur -9999 pour le donnes sans valeur (nodata)
Si le Systme de projection est inchang (SCR) inutile de le prciser
En ligne de commande, la commande taper est :
10/12
gdal_translate -a_nodata -9999 -of AAIGrid fichier_depart.tif fichier_arrivee.asc
voil vous avez obtenu un MNT raster au format ASC, du pas souhait et sur le territoire d'tude voulu !
11/12
8 Utilisation de GDAL sans QGIS
(je placerais plutt ce chapitre en fin de document)
Si certains calcultraitementss n'aboutissent pas dans QGIS vous avez la possibilit d'utiliser directement
GDAL en shell (invite de commande) :
1. Sous l'explorateur windows, placez-vous dans le rpertoire dexcution de QGIS (rpertoire
bin ) et copiez le chemin d'accs dans le presse papier ( CTRL+C)
2. Cliquez sur Dmarrer\Excuter puis tapez cmd et validez
3. Tapez cd (aller au rpertoire) et collez le chemin (clic droit dans l'invite puis coller) : pour obtenir
le rsultat suivant (exemple) : cd C:\Program Files\Quantum GIS\bin
4. validez en appuyant sur la touche entre
A partir de cette fentre vous pouvez copier coller toutes les lignes de commande de Gdal qui suivent
dans cette notice et utilisez ventuellement d'autres options disponibles grce GDAL.
12/12
