Table of Contents

Qu'est-ce qu'un MNT ?

Il s'agit d'un fichier RASTER, c'est-à-dire défini cellule à cellule. Chacune de ces cellules (les célèbres pixels) a pour valeur l'altitude moyenne de la zone et le fichier peut, dans certains formats, se présenter sous la forme d'une matrice humainement lisible (mais pas humainement interprétable !!). Un MNT peut donc être plus ou moins détaillé, plus ou moins fin. Le logiciel de SIG (dans notre cas, GRASS !!!) va ensuite interpréter ces valeurs avec des niveaux de couleur ou alors visualiser cette zone en 3D…

Où trouver des MNT ?

Plusieurs sources existent, pour disposer de MNT gratuits (c'est quand même pas mal…) : le SRTM (Shuttle Radar Topography Mission) a eu lieu en Février 2000, à bord de la navette spatiale Endeavour. Ces données couvrent PRATIQUEMENT l'ensemble de la surface de la Terre, avec une précision de 3 secondes d'arc, soit environ 90 mètres (c'est déjà pas mal !). Les données sont disponibles via un serveur ftp, par tuiles mesurant 1 degrés de côté : ftp://e0srp01u.ecs.nasa.gov/srtm/
Il s'agit de l'adresse actuelle (en décembre 2006), mais elle est susceptible de changer fréquemment… J'ai récupéré sur mon ordi les tuiles qui concernent la France, et je peux vous les envoyer si vous ne parvenez pas à les récupérer…

Si vous disposez des courbes de niveaux digitalisées, ou si vous voulez travailler à une résolution très fine (par exemple sur une micro-région), vous pouvez créer votre MNT par vous-même…

Les données du SRTM

Elles sont de bonnes qualités, avec cependant quelques lacunes parfois… Un script permet de les importer directement dans GRASS :
ATTENTION !! Les fichiers présentent la double extension .hgt.zip . VOUS NE DEVEZ PAS LES DEZIPPER !!! Si une telle manipulation se fait automatiquement, faites attention à conserver le fichier zippé !
Premièrement, pour pouvoir les importer, il faut que vous soyez dans une location en PROJECTION LATITUDE-LONGITUDE selon l'ellipsoïde WGS84… Si elle n'existe pas, créez la avec une résolution de 3 secondes d'arc, ou moins…
Si vous avez votre location adéquate, allons-y :
Dans l'interface graphique, sélectionnez :

File -> Import -> Raster map -> SRTM hgt files

Vous n'avez à préciser que le chemin d'accès au fichier récupéré et le nom du fichier de sortie (elev par exemple…).
En ligne de commande, vous tapez :

r.in.srtm input=xxx output=xxx 

Le traitement prend un certain temps.

Une fois qu'elles sont traitées, vous devez combler les vides (présence de nuages lors des mesures du SRTM…). Une commande redoutable existe pour celà : il s'agit du rééchantillonnage par interpolation de splines régularisées avec tension.
Sous ce nom barbare, se cache un outil très puissant qui fait des calculs très précis et de bonne qualité, y compris si l'on souhaite affiner la résolution (à 50 mètres par exemple).
Dans l'interface graphique, sélectionnez :

Raster -> Interpolate surfaces -> Fill NULL cells by interpolation using regularized spline tension

ou dans le terminal :

r.fillnulls input=xxx output=yyy tension=40. smooth=0.1

Dans l'interface graphique, les paramètres tension et smooth sont fixés avec les paramètres ci-dessus (soit 40 et 0.1), qui conviennent bien. input est le nom du raster d'origine et output celui du raster terminal.
Cette opération est TRES longue et TRES gourmande en puissance. Un conseil : pendant le traitement, ne faites rien sur votre ordinateur !!

Par la suite, vous pouvez augmenter la résolution en utilisant un outil tout à fait comparable… Dans l'interface graphique :

Raster -> Develop map -> Resample (change resolution) using regularized spline tension

ou dans le Terminal :

r.resamp.rst input=yyy ew_res=0.000003 ns_res=0.000003 elev=elev_fine overlap=3 zmult=1.0 tension=40. 

Les éléments que vous devez précisez sont :
* input, le nom du fichier de départ ;
* ew_res et ns_res, les résolutions en sortie ;
* elev qui est le nom du fichier de sortie.
Vous pouvez générer d'autres fichiers : pente (slope), aspect (aspect)… mais l'élévation est celle qui nous convient pour l'instant !

Faire un MNT à partir de courbes de niveau