SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
FICHE 12 : TRAITEMENTS SUR SRTM
Objectifs spécifiques
Visualiser et manipuler une image SRTM
Extraire et mettre en forme des données altimétriques
Extraire et mettre en forme des données hydrologiques
**************************
A partir du CD-Rom fourni, copiez le dossier SRTM sur le bureau
Les données SRTM sont issues de la Shuttle Radar Topography Mission, campagne
d’observation de la terre réalisée en février 2000 par la navette spatiale Endeavour (STS-99) à
une altitude de 233 km en utilisant l'interférométrie radar. Il existe 3 types de résolution pour
les données srtm, 1m pour les Etats-Unis, 30m et 90m pour le reste du monde.
Aperçu de l’image Srtm téléchargé sur Internet1 : le fond de carte représente la grille des
images Landsat disponibles en téléchargement pour le Cameroun (36 images).
1ère étape : Visualisation du srtm et extraction de la zone d’étude
Lancez ArcMap,
Ouvrez un nouveau document ArcMap, nommez-le « Traitements_srtm » et enregistrez le
dans le dossier « SRTM » que vous venez de copier sur le bureau. Ajoutez à la vue courante le
11
Adresse de téléchargement : [Link]
ENSTP - Topographie 1 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
fichier raster « Srtm_région_ydé.tif » situé dans le dossier « SRTM » puis dans le dossier
« Srtm_région_ydé ». Changez la symbologie du fichier en cliquant droit sur la couche dans
la légende gauche, puis en allant dans « Propriétés ». Dans la fenêtre qui s’ouvre
« Propriétés de la couche », puis dans l’onglet « Symbologie », choisissez d’afficher les
valeurs « classé » du srtm. Validez et visualisez le srtm dans la vue courante. Apres cela,
choisissez à nouveau d’afficher les valeurs « étiré », mais avant de validez votre choix,
sélectionnez dans la sous-rubrique ‘Etirer’ le ‘Type’ « Personnalisé » puis cliquez ensuite sur
le bouton « Histogrammes ».
Pour étirer les valeurs de l’histogramme, cliquez sur le graphique à l’endroit indiqué ci-
dessous puis validez vos choix. Visualisez à nouveau le srtm dans la vue courante.
Cliquez ici pour étirer l’histogramme (dans le violet) Résultat obtenu
Résultat : Les valeurs qui étaient concentrés au sein d’un faible intervalle de valeurs (proche
de 1566 sont ensuite dispersés sur l’ensemble de l’intervalle [-32768, 1566]
Aperçus au sein de la vue courante :
Srtm faiblement étiré Srtm fortement étiré Srtm classé (5 classes)
ENSTP - Topographie 2 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
Sous ArcMap,
Ajoutez à la vue courante le fichier vecteur « Centre_dpt.shp » situé dans le dossier
« SRTM », puis dans le dossier « Région Centre ». Sélectionnez le département du
« Mfoundi », créez une couche à partir de l’entité sélectionnée et exportez-là sous le nom de
« [Link] » dans le dossier « Région Centre ».
Lancez ArcToolbox,
Dans le menu « Outils Spatial Analyst » sélectionnez « Extraction » puis l’option « Extraction
par masque ». Dans la fenêtre « Extraction par masque » qui s’ouvre, choisissez comme
« Raster en entrée » le fichier raster « Srtm_région_ydé.tif », choisissez comme « Données
raster ou vecteur de masque en entrée », le fichier vecteur « [Link] ». Nommez et
enregistrez le nouveau « Raster en sortie », sous le nom de « Srtm_Mfoundi » dans le dossier
« SRTM ».
Dans la vue courante, le srtm à été découpé selon la forme du département du Mfoundi.
2ème étape : Extraction des données altimétriques
Les courbes de niveau
Sous « ArcToolbox », dans le menu « Outils Spatial Analyst » sélectionnez « Surface » puis
l’option « Isoligne ». Dans la fenêtre « Isoligne » qui s’ouvre, choisissez comme « Raster en
entrée », le fichier « Srtm_Mfoundi », choisissez comme « Intervalle des isolignes » la valeur
« 25 ». Nommez le fichier en sortie « Mfoundi_Cdn25.shp » et enregistrez le dans le dossier
« Altimétrie » que allez créer maintenant dans le dossier « SRTM ».
Dans la vue courante le fichier « Mfoundi_Cdn30.shp » a dessiné des courbes de niveau
espacées de 30m sur l’ensemble de la surface du Mfoundi.
ENSTP - Topographie 3 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
Créez également un nouveau fichier « Mfoundi_Cdn10.shp » pour obtenir des courbes de
niveau espacées de 10m et enregistrez-le dans le dossier « Altimétrie ».
L’ombrage
Sous « ArcToolbox », dans le menu « Outils Spatial Analyst » sélectionnez « Surface » puis
l’option « Ombrage ». Dans la fenêtre « Ombrage » qui s’ouvre, choisissez comme « Raster
en entrée », le fichier « Srtm_Mfoundi ». Ne modifiez pas les autres configurations
facultatives par défaut. Nommez le fichier raster en sortie « Mf_Ombrage » et enregistrez le
dans le dossier « Altimétrie ».
Créez également un nouveau fichier « Mf_Ombrage2.shp » en modifiant le « Facteur Z »
supprimez la valeur « 1 » inscrite par défaut et tapez « 2 » à sa place. Visualisez la différence
entre les deux rasters d’ombrage que vous avez créé. Les rasters d’ombrage sont le plus
souvent utilisés en fond sur certaines cartes avec un facteur de transparence plus ou moins
élevé pour donner une impression de relief et d’ombrage à celles-ci.
Les pentes
Sous « ArcToolbox », dans le menu « Outils Spatial Analyst » sélectionnez « Surface » puis
l’option « Pente ». Dans la fenêtre « Pente » qui s’ouvre, choisissez comme « Raster en
entrée », le fichier « Srtm_Mfoundi », choisissez comme « Mesure en sortie » des
« Degrés/DEGREE ». Nommez le fichier raster en sortie « Mf_Pente » et enregistrez le dans
le dossier « Altimétrie ».
Dans la vue courante, le fichier « Mf_Pente » a dessiné un raster indiquant le degré des pentes
et leur localisation.
ENSTP - Topographie 4 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
L’orientation des pentes (exposition)
Sous « ArcToolbox », dans le menu « Outils Spatial Analyst » sélectionnez « Surface » puis
l’option « Exposition ». Dans la fenêtre « Exposition » qui s’ouvre, choisissez comme
« Raster en entrée », le fichier « Srtm_Mfoundi ». Nommez le fichier raster en sortie
« Mf_Exposition » et enregistrez le dans le dossier « Altimétrie ».
Dans la vue courante, le fichier « Mf_Exposition » a dessiné un raster indiquant l’orientation
et la direction des pentes.
Le Modèle numérique de terrain
Sous « ArcToolbox », dans le menu « 3D Analyst » sélectionnez « Gestion des données » puis
« TIN », puis l’option « Créer un TIN ». Dans la fenêtre « Créer un TIN » qui s’ouvre,
choisissez comme « Classe d’entités en entrée », le fichier des courbes de niveau «
Mfoundi_Cdn25 ». Enregistrez le TIN sous le nom de « Mf_Tin » dans le dossier
« Altimétrie ».
Dans la vue courante, le fichier « Mf_Tin » a dessiné un raster montrant le relief du Mfoundi
avec une palette par défaut. Ce fichier peut être également visualisé en 3 dimensions sous le
module ArcScène.
ENSTP - Topographie 5 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
2ème étape : Extraction des données hydrologiques
La direction des flux
Sous ArcToolbox,
Dans le menu « Outils Spatial Analyst » sélectionnez « Hydrologie » puis l’option « Direction
de flux ». Dans la fenêtre « Direction de flux » qui s’ouvre, choisissez comme « Raster de
surface en entrée » le fichier raster « Srtm_Mfoundi », Nommez le nouveau fichier « Raster
de direction de flux en sortie », « Mf_Dirflux » et enregistrez le dossier « Hydrologie » que
vous allez créez maintenant dans le dossier « SRTM ».
Dans la vue courante, le fichier « Mf_Dirflux » a généré un raster montrant les éventuelles
directions des flux en cas d’écoulement hydrique du fait de l’orientation des pentes de la
surface analysée.
Les bassins
Sous ArcToolbox,
Dans le menu « Outils Spatial Analyst » sélectionnez « Hydrologie » puis l’option « Bassin ».
Dans la fenêtre « Bassin » qui s’ouvre, choisissez comme « Raster de direction de flux en
entrée » le fichier raster « Mf_Dirflux », Nommez le nouveau fichier « Raster en sortie »,
« Mf_Bassin » et enregistrez le dossier « Hydrologie ».
Dans la vue courante, le fichier « Mf_Bassin » a dessiné un raster découpant la surface
étudiée en plusieurs bassins d’écoulement selon la direction des flux et des pentes.
ENSTP - Topographie 6 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
Les écoulements (accumulation de flux)
Sous ArcToolbox,
Dans le menu « Outils Spatial Analyst » sélectionnez « Hydrologie » puis l’option
« Accumulation de flux ». Dans la fenêtre « Accumulation de flux » qui s’ouvre, choisissez
comme « Raster de direction de flux en entrée » le fichier raster « Mf_Dirflux ». Nommez le
nouveau fichier « Raster en sortie », « Mf_Accuflux » et enregistrez le dossier « Hydrologie ».
Dans la vue courante, le fichier « Mf_Accuflux » a généré un raster qui met en évidence
l’accumulation de faible valeurs d’altitude lieux de prédilection des écoulements de surfaces.
Pour obtenir un fichier vecteur des cours d’eau, nous allons d’abord reclasser ce raster en 2
classes (Eau et Non_eau), puis le convertir en vecteur. Pour connaitre les valeurs de classe, il
est intéressant de consulter l’histogramme de classification du fichier « Mf_Accuflux » dans
l’onglet « Symbologie » en cliquant droit sur la couche dans la légende de gauche et en allant
dans « Propriétés ». Si l’on choisit de visualiser le raster de façon « classé » et non « étiré »,
s’il l’on choisit 2 classes selon la méthode des seuils naturels de Jenks. Il semble que la
valeur 20 partage les pixels représentant l’eau de ceux ne représentant pas l’eau. Nous allons
donc reclasser le raster « Mf_Accuflux » selon les critères suivants : les pixels de 0 à 20
deviendront du « Nodata », les pixels de 20 à 338 prendront la valeur 1 (indiquant que les
pixels en question représente l’eau et les écoulements), pour finir les pixels qui étaient du
« Nodata », resteront du « Nodata ».
La reclassification
Sous ArcToolbox,
Dans le menu « Outils Spatial Analyst » sélectionnez « Reclassement » puis l’option
« Reclassification ». Dans la fenêtre « Reclassification » qui s’ouvre, choisissez comme
« Raster en entrée » le fichier raster « Mf_Accuflux ».
Définissez les nouvelles valeurs pour la reclassification :
De 0 à 20 = Nodata
De 20 à 338 = 1
Nodata = Nodata
Nommez le fichier raster en sortie « Mf_Accuflux2 » et enregistrez le dans le dossier
« Hydrologie ».
ENSTP - Topographie 7 Année 2017 / 2018
SIG – Travaux pratiques, fiche (12) Traitements sur srtm Enseignant : Fredy MBIANDA
Dans la vue courante le nouveau fichier raster « Mf_Accuflux2 » a gardé uniquement les
pixels correspondant à l’eau en leur donnant la valeur 1.
Pour transformer ce fichier raster en format vecteur nous allons utiliser dans ArcToolbox,
puis dans les « Outils Spatial Analyst » puis dans « Hydrologie », l’option « Ecoulement vers
entité ». Dans la fenêtre « Ecoulement vers entité », choisissez comme « Raster d’écoulement
en entrée », le fichier « Mf_Accuflux2 ». Choisissez comme « Raster de direction de flux en
entrée » le fichier « Mf_Dirflux ». Nommez le nouveau fichier vecteur (lignes) en sortie
« Ecoulement_Mfoundi » et enregistrez-le dans le dossier « Hydrologie ».
ENSTP - Topographie 8 Année 2017 / 2018
