SIG974

Quelquefois on voudrait récupérer le relief grossier d’un territoire, histoire de faire de belles cartes avec le relief en dégradé de gris. Mais voilà où trouver des données altimétriques gratuites et pas trop imprécises ?

Et puis on va bien avoir un MNT / MNE très précis élaboré par l’IGN mais il ne sera pas gratuit donc il faut se rabattre sur le SRTM, et c’est déjà pas mal !

En effet, depuis plus de 10 ans maintenant la NASA et la NGA (National Geospatial-Intelligence Agency) on monté un programme pour reconstituer le Modèle Numérique de Terrain du globe: La Shuttle Radar Topography Mission (SRTM).

Ils ont fabriqué, à l’aide d’un radar monté sur la navette Endeavour, des milliers de dalles de 5° par 5°.  Et comme les américains sont sympas ils en ont profité pour couvrir d’une seule dalle MNT,  Maurice et La Réunion.

Sur leur site, ils promettent une résolution inférieure à 90m. Alors dans la Beauce c’est pas terrible, mais à La Réunion, 90m de résolution c’est assez intéressant.

Ça donne une image comme celle-ci  (beaucoup de vide, hein !) :

Dalle SRTM Maurice / La Réunion

Et voici la description du MNT sur le site du CGIAR-CSI:

Product:  SRTM 90m DEM version 4
Data File Name: srtm_48_17.zip
Latitude min: 25 S
Latitude max: 20 S
Longitude min: 55 E
Longitude max: 60 E
Central Point : 22.50 S / 57.50 E

Évidemment ces données sont gratuites et utilisables directement dans Quantum GIS. Et si vous avez récupéré l’extension GdalTools vous pourrez créer des courbes de niveau sans problème.

Bien sûr voici tous les liens nécessaires:

• Données STRM sur le site du CGIAR-CSI

• Tout sur le format et la mission SRTM

• Le STRM de Maurice / La Réunion (1 Mo)

• Les courbes de niveau de La Réunion créées avec Quantum GIS (1,5 Mo)

Format:  shape
Équidistance: 50m
Projection RGR92 / UTM40Sud

Voici une astuce bien utile pour les adeptes des couches logiques dans MapInfo, et qui se trouvent un peu frustrés dans QGIS avec leurs 130 dalles BDOrtho dans le gestionnaire de couches !

Oui, c’est possible dans QGIS de faire une couche logique comme sous MapInfo avec le format VRT. Simplement c’est juste un peu plus long au départ. Et merci à Mathieu RAJERISON de m’avoir soufflé cette combine qui prend en tout 5 minutes

Alors on y va:

1. Il faut installer FWTools (2.4.7) : la boîte à outils en lignes de commande (don’t stress qui contient notamment GDAL/OGR.

2. Une fois FWTools installé:

• Ouvrir l’invite de commandes
• Se positionner dans le répertoire C:\FWTools2.4.7(à coup de commandes cd.. et cd FWTools2.4.7)
• Lancer le fichier setfw.bat pour initialiser les variables
• Aller dans le répertoire C:\FWTools2.4.7\bin (cd bin)
• Depuis ce répertoire, taper la commande suivante: gdalbuildvrt OrthoAssemblage.vrt D:/REFERENTIELS/BDORTHO/DALLES/*.tif

Le résultat est un fichier texte (de type xml) nommé OrthoAssemblage.vrt, situé dans le répertoire C:\FWTools2.4.7\bin,. Ce fichier sera le tableau d’assemblage de toutes les dalles de la BDOrtho.

3. Il suffit enfin d’ouvrir dans QGIS ce fichier par le menu Couche > Ajouter Raster > [typeVRT]

Et voilà, le tour est joué ! Une seule couche logique apparaît.

MàJ du 23 juillet 2010

Sinon pour les plus flemmards, ou ceux qui ont pas trop de dalles, une solution encore plus simple:

• Installer FWTools et lancer le fichier setfw.bat d’initialisation des variables de base
• Aller dans QGIS, installer le plugin GdalTools qui créera le menu Raster
• Utiliser alors le menu Raster > Construit un Raster Virtuel

Dans la série on remonte le temps, cette fois-ci on va du côté de Saint-André et de la côte est et jusqu’à Saint-Benoît .

Voici une nouvelle mosaïque géoréférencée avec QGIS, effectuée à partir de d’une trentaine de photographies aériennes prise en 1950 par l’Institut Géographique National.

[Si vous cliquez sur , vous pouvez alterner entre les photos de 1950 et GoogleMap !]

[iframe http://sig974.free.fr/zoom/1950/STANDRE/STANDRE.html 680 600]