Le portail francophone de la géomatique

Bonjour à tous, ceci n'est pas une question mais la solution au problème que j'ai rencontré.

J'imagine que pas mal de personnes étaient au courant de la manipulation à effectuer et cela va leur paraître simpliste, mais n'étant pas un grand spécialiste du SIG et n'ayant trouvé aucun post à ce sujet, je poste ici la réponse au cas où elle pourrait aider des personnes rencontrant le même soucis.

Je m'explique. Je disposais d'un MNT au format Mapinfo .grd . A son ouverture sous QGIS ou MapInfo, aucun soucis, il s'affichait en tant que Raster avec des couleurs différentes selon l'altitude. Cependant mon intérêt était de pouvoir extraire les courbes de niveau afin de les afficher par la suite pour réaliser des exports avec ortho + cadastres + courbe topo (je travaille en BE d'urbanisme).

Ne disposant pas de Vertical Mapper, je me tourne donc vers ce merveilleux QGIS et j'utilise l'outil GDAL Raster/Extraction/Création de contours. Par défaut cet outil utilise les valeurs de la bande 1 du MNT pour créer les courbes topo. Je charge la nouvelle couche et m'aperçois que le shp créé ne correspond pas du tout (courbes manquantes, pas les bonnes valeurs d'altitude). En utilisant l'outil "identifier les entités", je remarque que c'est à la bande 4 que sont associées les bonnes valeurs d'altitude mais comment l'indiquer au script (c'était mon problème). En regardant dans d'autres script j'ai trouvé la solution, la voici :

On utilise toujours le même outil GDAL création de contours, on renseigne l'ensemble des paramètres (fichiers de sortie, nom d'attribut) et là on clique sur le stylo pour éditer le script. Il suffit de rajouter à la fin de celui-ci

'espace'-b'espace''numéro de la bande'

dans mon cas

 -b 4

(voir le fichier joint)
En espérant que cela puisse aider quelques personnes

Merci à tous de faire vivre ce forum qui est très utile !

Thomas

Dernière modification par Kinad Vandar (Thu 15 November 2012 14:54)