Créer une carte de gravité des dommages causés par le feu
Télécharger les données et ouvrir un modèle de projet
Vous allez commencer par préparer l’analyse en téléchargeant les données nécessaires et en ouvrant un modèle de projet ArcGIS Pro.
- Téléchargez le fichier LandslideData.zip.
- Sur le lecteur C, créez un dossier nommé LandslideData.
- Accédez au dossier Downloads (Téléchargements) et extrayez le contenu du fichier .ZIP que vous avez téléchargé dans le dossier LandslideData.
- Démarrez ArcGIS Pro. Si vous y êtes invité, connectez-vous via votre compte d’organisation ArcGIS sous licence.
Remarque :
Si vous n’avez pas accès à ArcGIS Pro ou ne disposez pas d’un compte d’organisation ArcGIS, consultez les options disponibles pour accéder aux logiciels.
- Cliquez sur Start with another template (Démarrer avec un autre modèle).
- Accédez à C:\LandslideData, choisissez le fichier du modèle de projet, Landslide_Risk_Project.aptx, puis cliquez sur OK.
- Dans la fenêtre New Project (Nouveau projet), nommez votre projet Landslide_Risk et enregistrez-le à l’emplacement par défaut (dans le dossier Documents) sur votre ordinateur. Cliquez sur OK.
Le nouveau projet s’ouvre et affiche une carte centrée sur Santa Rosa, en Californie.
Vous allez ensuite établir une connexion au dossier LandslideData afin d’accéder rapidement aux données.
- Sur le ruban, cliquez sur l’onglet Insert (Insertion). Dans le groupe Project (Projet), cliquez sur Add Folder (Ajouter un dossier).
La fenêtre Add Folder Connection (Ajouter une connexion aux dossiers) apparaît.
- Accédez au dossier LandslideData que vous avez créé. Sélectionnez-le et cliquez sur OK (OK).
Maintenant que le projet est ouvert et qu’une connexion au dossier est établie pour les données, vous êtes prêt à réaliser l’analyse.
Analyser la gravité des dommages causés par le feu
Les différents types de végétation et de matériaux brûlent à des intensités diverses en fonction de leur composition, de leur densité, de la topographie, du vent, de l’humidité du sol et d’autres facteurs. En général, les intensités les plus élevées découlent d’un pouvoir hydrophobe et d’une possibilité d’érosion supérieurs en cas de pluie. L’intensité des dommages causés par le feu, ou gravité, peut être dérivée d’une image multispectrale comportant des canaux infrarouge à ondes courtes et proche infrarouge, telle qu’une imagerie Landsat 8. Grâce aux fonctions raster dans une chaîne de traitement, vous allez calculer le ratio de dommages causés par le feu en exploitant l’imagerie Landsat 8 acquise avant et après l’incendie qui s’est produit autour de Santa Rosa, en Californie, en 2017.
Pour commencer, vous allez ajouter l’imagerie Landsat 8 avant et après l’incendie, un modèle numérique d’élévation et des jeux de données raster d’occupation du sol.
- Sur le ruban, cliquez sur l'onglet Vue. Dans le groupe Windows (Fenêtres), choisissez Catalog Pane (Fenêtre du catalogue).
La fenêtre Catalogue apparaît.
- Dans l’onglet Project (Projet) de la fenêtre Catalog (Catalogue), développez Folders (Dossiers), puis le dossier LandslideData. Sélectionnez tous les jeux de données à l’intérieur de ce dossier, excepté Basins.tif, puis faites-les glisser sur la carte.
Conseil :
Vous pouvez sélectionner plusieurs éléments simultanément en appuyant sur la touche Ctrl et en cliquant sur les fichiers de votre choix.
Les données apparaissent dans la fenêtre Contents (Contenu) et sur la carte.
Les données comportent les images Landsat 8 précédant (Before_L8.tif) et suivant (After_L8.tif) les incendies d’octobre 2017 dans les comtés de Napa et Sonoma. Il existe également deux couches que vous utiliserez en entrée pour établir votre carte de risque. La couche DEM_30m.tif est un modèle numérique d’élévation affichant l’élévation du terrain. La couche Sonoma_NLCD2011.tif est une partie du jeu de données d’occupation des sols américains qui présente l’utilisation du sol et le type de végétation prédominant.
- Dans la fenêtre Contents (Contenu), faites glisser la couche After_L8.tif en tête de liste et placez la couche Before_L8.tif juste en dessous. Désactivez et réduisez les couches Sonoma_NLCD2011.tif et DEM_30m.tif.
Actuellement, l’imagerie Landsat est affichée sur la carte à l’aide des canaux rouge, vert et bleu. Pour comparer les zones brûlées, vous allez afficher certains des canaux multispectraux de l’imagerie.
- Sous After_L8.tif, cliquez avec le bouton droit de la souris sur la pastille de couleur rouge et choisissez sr_band5.
- Pour le vert, sélectionnez sr_band4 et pour le bleu, sr_band3.
Cette combinaison de canaux affiche les canaux de l’imagerie Landsat 8 en mode de couleur infrarouge. La végétation est indiquée en rouge vif. Les entités non couvertes de végétation telles que les zones nues et les zones urbaines s’affichent dans différentes nuances de gris et de bleu.
- Définissez la couche Before_L8.tif comme suit :
- Associez la pastille de couleur rouge à la valeur sr_band5.
- Associez la pastille de couleur verte à sr_band4.
- Associez la pastille de couleur bleue à la valeur sr_band3 ou vérifiez qu’elle est déjà sélectionnée.
À présent, vous allez utiliser l’outil Swipe (Balayer) pour comparer l’imagerie avant et après l’incendie.
- Dans la fenêtre Contents (Contenu), cliquez sur le paramètre After_L8.tif pour le sélectionner.
- Sur le ruban, cliquez sur l’onglet Raster Layer (Couche raster). Dans le groupe Compare (Comparer), cliquez sur l’outil Swipe (Balayer).
Le pointeur prend la forme d’une flèche dans l’affichage cartographique.
- Cliquez tout en effectuant un balayage de l’image sur la carte pour comparer l’imagerie avant et après l’incendie.
Vous pouvez effectuer un balayage vertical ou horizontal. Notez que de nombreuses zones sont en rouge dans Before_L8.tif et en gris ou vert dans After_L8.tif, pour indiquer la végétation perdue.
- Sur le ruban, cliquez sur l’onglet Map (Carte). Dans le groupe Navigate (Naviguer), cliquez sur Explore (Explorer).
Le pointeur retrouve son apparence habituelle. Puisque toutes les données figurent sur la carte, vous allez utiliser les fonctions raster pour calculer la gravité des dommages causés par le feu.
- Dans le ruban, cliquez sur l’onglet Imagerie. Dans le groupe Analysis (Analyse), choisissez Raster Functions (Fonctions raster).
La fenêtre Fonctions raster apparaît. L’onglet System (Système) présente les catégories de fonctions disponibles pour l’analyse raster. Pour ce didacticiel, deux modèles de fonction raster ont été créés pour vous. Ces modèles de fonctions personnalisées sont répertoriés sous l’onglet Project (Projet).
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur l’onglet Project (Projet).
Deux modèles de fonction raster sont inclus dans ce projet : Landcover_Remap (Classification_occupation_du_sol) et Burn_Severity (Gravité_des_dommages_causés_par_le_feu).
Remarque :
Si vous ne voyez pas les modèles de fonction raster, procédez comme suit :
Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur l’onglet Custom (Personnalisé). En regard de Landslide Grenada (Glissement de terrain à Grenade), cliquez sur le bouton Import functions (Importer des fonctions). Dans la fenêtre Select Processing Templates (Sélectionner les modèles de traitement), parcourez Folders > LandslideData (Dossiers > LandslideData) et cliquez sur Landcover_Remap.rft.xml, puis sur OK. Répétez cette procédure et ajoutez Burn_Severity.rft.xml à partir du même dossier.
- Cliquez avec le bouton droit sur Landcover_Remap (Classification_occupation_du_sol), pointez le curseur de la souris sur Move to (Déplacer vers), puis sur Custom (Personnalisé) et cliquez sur Custom1 (Personnalisé1).
Le fait de déplacer ces fonctions dans une catégorie personnalisée permet d’enregistrer toutes les mises à jour que vous apportez aux modèles de fonction raster dans l’éditeur de fonctions. Les modifications apportées à la catégorie Project (Projet) seront perdues si le projet n’est pas enregistré.
- Cliquez sur l’onglet Project (Projet). Cliquez avec le bouton droit de la souris sur Burn_Severity (Gravité_des_dommages_causés_par_le_feu) et déplacez cet élément dans Custom1 (Personnalisé1).
À présent, vous allez ouvrir et explorer le modèle de fonction raster Burn_Severity (Gravité_des_dommages_causés_par_le_feu).
- Sous l’onglet Custom (Personnalisé), cliquez avec le bouton droit sur le modèle Burn_Severity (Gravité_des_dommages_causés_par_le_feu) et choisissez Edit (Mettre à jour).
L’éditeur de fonctions s’ouvre et affiche la chaîne de traitement.
Les fonctions Band Arithmetic (Arithmétique de canal) convertissent les pixels des images en expressions. Les images post-incendie sont soustraites des images pré-incendie et traitées via une fonction de classification. La fonction de classification organise les valeurs de pixels en cinq catégories selon l’importance des dommages causés par le feu. Les seuils des cinq valeurs de gravité des dommages causés par le feu sont fournis par une étude d’évaluation des paysages (Key et Benson, 2005). La fonction Attribute Table (Table attributaire) dans la chaîne de traitement affecte un dégradé de couleurs à la carte de gravité des dommages causés par le feu. Celle-ci a déjà été créée pour vous.
- Fermez le modèle de fonction raster Burn Severity (Gravité des dommages causés par le feu) dans la fenêtre Function Editor (Éditeur de fonctions).
Maintenant que vous avez compris comment fonctionne cette fonction raster, vous allez l’utiliser pour calculer la gravité des dommages causés par le feu dans la zone de projet.
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur le modèle Burn_Severity (Gravité_des_dommages_causés_par_le_feu).
La fonction raster Burn_Severity Properties (Propriétés de Gravité_des_dommages_causés_par_le_feu) apparaît.
- Dans le champ Pre-Fire Imagery (Imagerie pré-incendie), choisissez Before_L8.tif et pour Post-Fire Imagery (Imagerie post-incendie), choisissez After_L8.tif.
- Assurez-vous que Output Layer Type (Type de couche en sortie) est défini sur Raster Layer (Couche raster), puis cliquez sur Create new layer (Créer une nouvelle couche).
Le traitement peut prendre quelques minutes. Lorsqu’il est terminé, la couche obtenue s’affiche dans la carte et son nom figure dans la fenêtre Contents (Contenu). Par nature, les fonctions raster sont temporaires ; les calculs sont effectués à la volée, ou en temps réel, au fur et à mesure que vous vous déplacez sur la carte et ne sont pas enregistrés automatiquement. La gravité des dommages causés par le feu est calculée de façon dynamique dans l’affichage lorsque vous naviguez dans la couche.
Vous pouvez faire défiler la carte et localiser les zones les plus affectées. Comme indiqué dans la légende, ces zones s’affichent en orange et en rouge sur la carte.
- Dans la barre d’outils Accès rapide, cliquez sur le bouton Save (Enregistrer).
Vous avez créé et exécuté un modèle de fonction raster afin de créer une couche de gravité des dommages causés par le feu.
Maintenant que vous avez découvert les fonctions raster, vous allez apprendre dans la section suivante à créer un modèle de fonction raster afin de créer une carte des indices de pente.
Créer une carte des indices de pente
Générer un modèle de fonction raster de pente
La carte de pente est une couche critique dans la détermination de la stabilité de la pente. La déclivité de la pente est dérivée d’un modèle numérique de terrain (MNE). Plus la pente est raide, plus elle est susceptible de faire l’objet d’un glissement de terrain, notamment en cas de fortes pluies après une disparition soudaine (suite à un incendie, par exemple) de la végétation qui stabilisait le terrain. Vous allez maintenant construire et enregistrer un modèle de fonction raster pour calculer le pourcentage de pente et l’utiliser pour créer un indice d’escarpement.
- Dans le ruban, cliquez sur l’onglet Imagerie. Dans le groupe Analysis (Analyse), cliquez sur Function Editor (Éditeur de fonctions).
La fenêtre Function Editor (Éditeur de fonctions) s’ancre en bas de la fenêtre de la carte.
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur l’onglet System (Système ). Développez le groupe Surface.
- Faites glisser la fonction Slope (Pente) dans la fenêtre Function Editor (Éditeur de fonctions).
- Dans la fenêtre Function Editor (Éditeur de fonctions), sur le ruban, cliquez sur Add Raster Variable (Ajouter une variable raster).
Une zone verte intitulée Raster (Raster) est ajoutée à la fenêtre Function Editor (Éditeur de fonctions).
- Sélectionnez la zone Raster et placez-la à gauche de la zone de fonction Slope (Pente).
La zone Raster (Raster) définira le jeu de données en entrée de la fonction Slope (Pente).
- Pointez le curseur de la souris sur la zone Raster pour afficher le paramètre Out (Sortie). Cliquez sur Out (Sortie) et faites-le glisser pour connecter l’élément raster (paramètre Out) à la fonction Slope (Pente) (paramètre DEM).
- Cliquez avec le bouton droit sur Raster et cliquez sur Rename (Renommer). Saisissez Input DEM (MNE en entrée) et appuyez sur Entrée.
- Double-cliquez sur la fonction Slope (Pente). Dans la fenêtre Slope Properties (Propriétés de la pente), cliquez sur l’onglet Variables. Pour DEM (MNE), sélectionnez le champ IsPublic (EstPublic).
L’option IsPublic permet de modifier les données en entrée dans un traitement d’outil ultérieur.
- Cliquez sur OK.
Votre fonction raster utilisera un modèle numérique d’élévation en entrée et calculera sa pente. La fonction suivante, Remap (Classification), classifiera la pente dans l’une des cinq catégories de déclivité.
- Dans la fenêtre Raster Functions (Fonctions raster), sur l’onglet System (Système), recherchez Remap (Classification).
- Faites glisser la fonction Remap (Classification) dans la fenêtre Function Editor (Éditeur de fonctions) à droite de la fonction Slope (Pente).
- Connectez Slope (Pente) à Remap (Classification) en définissant le paramètre Out (Sortie) de Slope (Pente) sur l’entrée Raster pour Remap (Classification).
Vous allez ensuite définir les entrées de la classification afin d’indexer les pentes (exprimées en degrés) et les classer en cinq catégories.
- Double-cliquez sur la fonction Remap (Classification). Dans la fenêtre Remap Properties (Propriétés de classification), cliquez sur le premier champ sous l’en-tête Minimum et entrez 0. Dans Maximum, tapez 5 et dans Output (Sortie), tapez 1.
- Définissez d’autres catégories avec les valeurs suivantes :
Minimum Maximum Sortie 2
5
15
2
3
15
25
3
4
25
35
4
5
35
91
5
- Cliquez sur l’onglet Général. Pour Output Pixel Type (Type de pixel en sortie), dans le menu déroulant, choisissez 8 Bit Signed (Signé 8 bits).
L’outil Attribute Table (Table attributaire), qui sera le prochain outil à ajouter, n’accepte que les rasters en entrée de 8 bits. Vous allez donc définir le paramètre Output Pixel Type (Type de pixel en sortie) sur 8 Bit Signed (Signé 8 bits).
- Cliquez sur OK.
- Cliquez avec le bouton droit sur la fonction Remap (Classification) et sélectionnez Rename (Renommer). Saisissez Slope Remap (Classification de la pente) et appuyez sur Entrée.
Ainsi, vous pourrez distinguer les fonctions ultérieurement dans ce didacticiel lorsque vous combinerez plusieurs chaînes de fonctions raster.
- Dans la fenêtre Raster Functions (Fonctions raster), recherchez la fonction Attribute Table (Table attributaire) et faites-la glisser dans la fenêtre Function Editor (Éditeur de fonctions) à droite de Slope Remap (Classification de la pente).
- Associez la sortie de la fonction Slope Remap (Classification de la pente) à l’entrée (Raster) de la fonction Attribute Table (Table attributaire).
- Double-cliquez sur la fonction Attribute Table (Table attributaire) et réglez Table Type (Type de table) sur Manual (Manuel).
- Sous la table vide, cliquez sur le bouton Generate (Générer).
- Pour Maximum Value (Valeur maximale), sélectionnez 5. Cliquez sur OK.
Cinq lignes (de 1 à 5) sont ajoutées ainsi qu’une palette par défaut allant du vert au rouge.
- Cliquez sur Class Name (Nom de la classe) pour chaque ligne et affectez les classes de pente comme suit :
Valeur Nom de la classe 1
Plate
2
Faible
3
Modéré
4
Raide
5
Très raide
- Cliquez sur OK.
Maintenant que votre modèle de fonction raster est terminé vous allez l’enregistrer.
Enregistrer et exécuter la fonction Pente
Enregistrez le modèle de fonction raster dans votre projet et exécutez-le pour créer la couche d’indice de pente.
- Dans la fenêtre Function Editor (Éditeur de fonctions), cliquez sur Save As (Enregistrer sous).
La fenêtre Save As (Enregistrer sous) s’ouvre.
- Dans la fenêtre Save As (Enregistrer sous), pour Name (Nom), saisissez Slope_Index (Indice_pente).
- Assurez-vous que Category (Catégorie) est défini sur Custom (Personnalisé) et que Sub-Category (Sous-catégorie) est défini sur Custom1 (Personnalisé1).
- Pour Description, copiez et collez A raster function template to derive slope from an input DEM (Modèle de fonction raster permettant de déduire la pente d’un MNE en entrée).
- Cliquez sur OK et effacez la barre de recherche de la fenêtre Raster Functions (Fonctions raster).
Votre modèle de fonction raster Slope_Index (Indice_pente) apparaît désormais dans la catégorie Custom (Personnalisé) dans la fenêtre Raster Functions (Fonctions raster).
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur Slope_Index (Indice_pente).
- Pour Input DEM (MNE en entrée), choisissez DEM_30m.tif. Assurez-vous que Output Layer Type (Type de couche en sortie) est défini sur Raster Layer (Couche raster).
- Cliquez sur Créer une nouvelle couche.
Lorsque le traitement est terminé, la couche apparaît sur la carte et figure dans la fenêtre Contents (Contenu), sous le nom Slope_Index_DEM_30m.tif.
- Fermez le modèle de fonction raster Slope_Index (Indice_pente) dans la fenêtre Function Editor (Éditeur de fonctions) et enregistrez votre projet.
Maintenant que vous disposez d’une fonction permettant de classifier les pentes, vous allez combiner un certain nombre de modèles de fonction raster afin de créer une carte des risques de glissement de terrain.
Créer une carte des risques de glissement de terrain
Générer la fonction raster de risque de glissement de terrain
Le calcul du risque de glissement de terrain combine deux variables que vous venez d’utiliser : la gravité des dommages causés par le feu et la pente. Il inclut aussi l’occupation du sol, autre facteur déterminant en matière de risque de glissement de terrain. En effet, la végétation stabilise les pentes par l’intermédiaire de ses systèmes racinaires. Mais, les incendies, ou feux de forêt, peuvent faire disparaître une bonne partie de la végétation stabilisatrice. Certains types de végétation, en particulier les espèces constituant le chaparral, se sont adaptés au feu et leurs racines atteignent des profondeurs particulièrement élevées afin d’y résister. L’occupation du sol a déjà fait l’objet d’une indexation préalable en cinq catégories, suivant son effet stabilisateur sur les pentes. Pour calculer le risque de glissement de terrain, vous allez ajouter les trois modèles de fonction raster dans une même chaîne et les traiter.
- Sur le ruban, dans l’onglet Imagery (Imagerie), cliquez sur Function Editor (Éditeur de fonctions).
La fenêtre Function Editor (Éditeur de fonctions) s’ouvre.
- Dans la fenêtre Raster Functions (Fonctions raster), appuyez sur la touche Ctrl et sélectionnez les modèles de fonction raster Landcover_Remap (Classification_occupation_du_sol), Burn_Severity (Gravité_des_dommages_causés_par_le_feu) et Slope_Index (Indice_pente). Faites-les glisser dans la zone d’affichage Function Editor (Éditeur de fonctions).
Par défaut, les modèles de fonction raster sont regroupés. Vous allez les séparer de manière à pouvoir connecter plus facilement leur sortie.
- Cliquez et maintenez enfoncé le bouton de la souris pour tracer un rectangle autour des fonctions Slope_Index (Indice_pente), puis faites glisser le groupe entier de sorte qu’il se trouve sous le modèle de fonction raster Burn_Severity (Gravité_des_dommages_causés_par_le_feu).
Conseil :
Les zones en entrée de couleur verte indiquent le début de chaque chaîne de fonctions. Le modèle de fonction raster Slope_Index (Indice_pente) commence avec le champ vert Input DEM (MNE en entrée).
- Faites glisser le modèle de fonction raster Landcover_Remap (Classification_occupation_du_sol) sous le modèle de fonction raster Slope_Index (Indice_pente).
- Cliquez avec le bouton droit sur chaque fonction Attribute Table (Table attributaire) et renommez chacune d’elles de manière à les associer à la fonction Remap (Classification) correspondante.
Leurs nouveaux noms doivent être Burn Severity Attribute Table (Table attributaire de la gravité des dommages causés par le feu), Slope Attribute Table (Table attributaire de la pente) et Landcover Attribute Table (Table attributaire de l’occupation du sol).
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur l’onglet System (Système) et recherchez la fonction Weighted Overlay (Superposition pondérée). Faites glisser la fenêtre Function Editor (Éditeur de fonctions) à droite des trois autres modèles de fonction raster.
- Connectez les trois sorties Attribute Table (Table attributaire) au paramètre en entrée de la fonction Weighted Overlay (Superposition pondérée).
- Dans la fenêtre Function Editor (Éditeur de fonctions), cliquez sur le bouton Auto Layout (Mise en page automatique).
Les modèles de fonction raster sont présentés dans un format compact.
- Double-cliquez sur la fonction Weighted Overlay (Superposition pondérée).
La fenêtre Weighted Overlay Properties (Propriétés de la superposition pondérée) apparaît. Dans la Weighted Overlay Table (Table de superposition pondérée), vous pouvez affecter des pondérations en pourcentage à chaque raster.
- Dans la Weighted Overlay Table (Table de superposition pondérée), dans la cellule en regard de <Burn Severity Attribute Table.OutputRaster>, saisissez 30. Pour la couche <Slope Remap Attribute Table.OutputRaster> (<Table attributaire de classification de la pente.Raster en sortie>), indiquez 55 pour cent. Pour la couche <Landcover Remap Attribute Table.OutputRaster> (<Table attributaire de classification de l’occupation du sol.Raster en sortie>), indiquez 15 pour cent.
Ces pondérations de danger reposent sur les recherches menées par l’USGS dans le cadre du programme national sur les risques de glissement de terrain.
La table Remap Table (Table de classification) est toujours vide. Les trois couches possédant le même nombre de catégories d’indices, vous allez associer chaque couche avec une relation un vers un.
- Dans Weighted Overlay Table (Table de superposition pondérée), cliquez sur la couche <Burn Severity Attribute Table.OutputRaster> (<Table attributaire de la gravité des dommages causés par le feu.Raster en sortie>). Dans la Remap Table (Table de classification), sous Value (Valeur), cliquez sur NODATA pour mettre à jour le champ attributaire et saisissez 1. Sous Scale (Échelle), cliquez sur NODATA et sélectionnez 1.
- Dans la Remap Table (Table de classification), double-cliquez sur la ligne vide en bas de la colonne Value (Valeur) et saisissez 2. Pour Scale (Échelle), sélectionnez 2.
- Répétez les étapes 1 à 5 pour les trois rasters de la Weighted Overlay Table (Table de superposition pondérée).
- Cliquez sur OK.
- Dans la fenêtre Raster Functions (Fonctions raster), onglet System (Système), recherchez la fonction Attribute Table (Table attributaire) et faites-la glisser vers la Function Editor (Éditeur de fonctions), à droite de la fonction Weighted Overlay (Superposition pondérée).
- Associez la sortie de la fonction Weighted Overlay (Superposition pondérée) à l’entrée de la fonction Attribute Table (Table attributaire).
- Double-cliquez sur la fonction Attribute Table (Table attributaire). Dans la fenêtre Attribute Table Properties (Propriétés de la table attributaire), pour Table Type (Type de table), choisissez Manual (Manuel).
- Sous la table vide, cliquez sur le bouton Generate (Générer) et définissez Maximum Value (Valeur maximale) sur 5. Cliquez sur OK.
Cinq lignes (de 1 à 5) sont ajoutées, ainsi qu’une palette par défaut allant du vert au rouge.
- Cliquez sur Class Name (Nom de la classe) pour chacune des lignes et affectez les valeurs suivantes :
Valeur Nom de la classe 1
Faible
2
Moyenne
3
Elevée
4
Très élevé
5
Extrême
- Cliquez sur OK pour fermer la fenêtre Attribute Table Properties (Propriétés de la table attributaire).
Enregistrer et exécuter la fonction raster de risque de glissement de terrain
Vous allez ensuite enregistrer le modèle de fonction raster et l’exécuter pour créer la couche des risques de glissement de terrain. Ce modèle considère la pente, la gravité des dommages causés par le feu et l’occupation du sol comme les principaux critères pour déterminer le risque de glissement de terrain.
- Dans la fenêtre Function Editor (Éditeur de fonctions), cliquez sur Save As (Enregistrer sous).
- Dans la fenêtre Save As (Enregistrer sous), pour Name (Nom), saisissez Landslide_Risk. Assurez-vous que Category (Catégorie) est défini sur Custom (Personnalisé) et que Sub-Category (Sous-catégorie) est défini sur Custom1 (Personnalisé1).
- Pour Description, copiez et collez Raster Function Template to calculate landslide risk based on wildfire burn severity, slope, and landcover (Modèle de fonction raster permettant de calculer le risque de glissement de terrain en fonction de la gravité des dommages causés par le feu, de la pente et de l’occupation du sol).
- Cliquez sur OK et fermez la fenêtre Function Editor (Éditeur de fonctions).
Votre fonction raster est terminée. À présent, vous allez l’exécuter en utilisant l’analyse raster distribuée.
- Dans la fenêtre Raster Functions (Fonctions raster), effacez votre recherche.
- Si nécessaire, cliquez sur l’onglet Custom (Personnalisé). Cliquez sur Landslide_Risk.
La fonction raster Landslide_Risk (Risque_glissement_de_terrain) apparaît.
- Renseignez les paramètres en sortie comme suit :
- Pour Pre-Fire Imagery (Imagerie pré-incendie), choisissez Before_L8.tif.
- Pour Post-Fire Imagery (Imagerie post-incendie), choisissez After_L8.tif.
- Pour Slope Input DEM (Pente - MNE en entrée), choisissez DEM_30m.tif.
- Pour Landcover Remap Raster (Occupation du sol - Raster de classification), choisissez Sonoma_NLCD2011.tif.
- Pour Output Layer Type (Type de couche en sortie), sélectionnez Raster Layer (Couche raster) et cliquez sur Create new layer (Créer une nouvelle couche).
- Fermez le modèle de fonction raster Slope_Index (Indice_pente) dans la fenêtre Function Editor (Éditeur de fonctions) et enregistrez votre projet.
Vous avez créé la couche des risques de glissement de terrain pour la zone en ajoutant plusieurs critères et fonctions à une fonction. Même si la couche des risques de glissement de terrain est utile par elle-même afin d’identifier les zones à risque, il est toujours possible d’optimiser une analyse.
Vous allez affiner votre analyse en tenant compte des bassins versants au sein de la zone d’étude, car ils ont également un impact sur le risque de glissement de terrain.
Synthétiser le risque de glissement de terrain selon le sous-bassin du bassin versant
Effectuer une analyse raster
La carte des risques de glissement de terrain est utile, mais vous désirez en savoir plus et isoler les zones présentant le risque le plus élevé. Étant donné que le risque de glissement de terrain est impacté par les précipitations et les caractéristiques des bassins versants, vous allez synthétiser les risques en fonction des bassins versants au sein de la zone d’étude à l’aide d’une fonction raster. Vous allez ensuite partager la couche obtenue sur ArcGIS Online, où vous pourrez l’ajouter à une application ou une carte Web et la partager.
- Dans la fenêtre Catalog (Catalogue), cliquez sur l’onglet Project (Projet), puis accédez à la connexion au dossier LandslideData.
- Ajoutez le raster Basins.tif à la carte actuelle.
Il s’agit de la couche qui vous servira à synthétiser le risque de glissement de terrain. Vous allez utiliser une fonction raster nommée Zonal Statistics (Statistiques zonales) pour réaliser cette analyse.
- Dans la fenêtre Raster Functions (Fonctions raster), cliquez sur l’onglet System (Système ). Dans la zone de recherche, saisissez zonal et ouvrez la fonction Zonal Statistics (Statistiques zonales).
- Dans la fenêtre Zonal Statistics Properties (Propriétés des statistiques zonales), renseignez les paramètres suivants :
- Pour Zone Raster (Raster de zones), sélectionnez Basins.tif.
- Pour Value Raster (Raster de valeur), sélectionnez Landslide_Risk.
- Acceptez les autres valeurs par défaut.
- Cliquez sur Create new raster (Créer un nouveau raster).
- Dans la fenêtre Contents (Contenu), renommez la couche en sortie des statistiques zonales en Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]). Désactivez toutes les couches à l'exception de la couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]) et du fond de carte.
Une fois l’analyse terminée, un jeu de données illustrant le risque moyen par bassin versant est ajouté à la carte.
![Couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa])](GUID-133863D8-1DD4-411D-B43A-6A155B32C0EA-web.png "Couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa])")
Les zones en blanc et en gris clair sont plus exposées au risque de glissement de terrain. La carte est générée selon le risque que vous avez calculé à l’aide du modèle de fonction raster faisant intervenir les entrées pondérées de pente, de gravité de dommages causés par le feu et d’occupation du sol.
Symboliser et partager vos résultats d’analyse
Vous avez localisé les zones les plus exposées au risque de glissement de terrain selon les bassins versants. Les résultats actuels sont cependant difficiles à distinguer avec cette symbologie. Vous allez donc la modifier afin de faciliter leur compréhension.
- Dans la fenêtre Contents (Contenu), cliquez avec le bouton droit sur la couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]) et sélectionnez Symbology (Symbologie).
La fenêtre Symbology (Symbologie) apparaît.
- Pour Color scheme (Combinaison de couleurs), cliquez sur la combinaison de couleurs actuelle Black to White (Noir vers blanc).
- Dans le menu qui s’affiche, cliquez sur Show names (Afficher les noms) et sélectionnez la combinaison de couleurs Slope (Pente).
La combinaison de couleurs correspond désormais à la sortie des autres couches que vous avez créées. Vous allez fusionner la couche pour améliorer son affichage.
- Dans la fenêtre Contents (Contenu), vérifiez que la couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]) est sélectionnée.
- Sur l’onglet Raster Layer (Couche raster), dans le groupe Effects (Effets), pour Layer Blend (Fusion des couches), cliquez sur le menu déroulant et choisissez Multiply (Multiplication).
La couche s’assombrit et une transparence lui est appliquée.
- Zoomez sur la couche pour l’agrandir.
L’aspect de la couche s’est amélioré et vous pouvez voir le fond de carte sous-jacent. La carte est maintenant prête à être partagée.
- Dans la fenêtre Contents (Contenu), cliquez avec le bouton droit sur la couche Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]), pointez sur Sharing (Partage) et sélectionnez Share As Web Layer (Partager en tant que couche Web).
- Pour Name (Nom), si des traits de soulignement sont présents, supprimez-les. Pour Layer Type (Type de couche), sélectionnez Tile (Tuile).
- Cliquez sur l’onglet Configuration.
Vous allez modifier certains paramètres afin de contrôler la projection en sortie et le mode de gestion de la mise en cache. Comme vous souhaitez que la projection en sortie soit Web Mercator, vous allez sélectionner la structure de tuilage ArcGIS Online.
Remarque :
Si vous ne configurez pas la structure de tuilage et les options de mise en cache, votre couche sera publiée, mais elle ne sera pas utilisable dans ArcGIS Online.
- Pour Tiling Scheme (Structure de tuilage), cliquez sur le menu déroulant et sélectionnez ArcGIS Online / Bing Maps / Google Maps.
- Développez Options et cliquez sur Cache locally (Mettre en cache en local).
- Cliquez sur Publish (Publier).
- Au bas de la fenêtre Share As Web Layer (Partager en tant que couche Web), cliquez sur le lien Manage web layer (Gérer la couche Web).
La page des détails des éléments de la couche de tuiles s’affiche.
- Cliquez sur Open in Map Viewer (Ouvrir dans Map Viewer).
La couche s’affiche, mais la fusion s’est perdue au cours du partage. Vous pouvez rapidement fusionner la couche dans la carte Web afin d’optimiser son apparence.
- Dans la fenêtre Properties (Propriétés), dans la section Appearance (Apparence), pour Blending (Fusion), cliquez sur le menu déroulant et sélectionnez Multiply (Multiplication).
- Fermez la fenêtre Blending (Fusion) et la fenêtre Properties (Propriétés).
À présent, la couche des glissements de terrain par sous-bassin est symbolisée de la même façon que dans ArcGIS Pro. Vous allez maintenant enregistrer la carte Web afin de la partager ou de l’intégrer dans d’autres applications, par exemple un récit ou une Instant App.
- Dans la barre d’outils (sombre) Contents (Contenu), cliquez sur Save and open (Enregistrer et ouvrir).
- Dans la fenêtre Save map, pour Name (Nom), saisissez Landslide Risk by Subbasin - Napa & Sonoma (Risque de glissement de terrain par sous-bassin [Napa et Sonoma]) et cliquez sur Save (Enregistrer).
La carte Web fait désormais partie de votre contenu. Vous pouvez la partager avec votre organisation ou le public telle qu’elle ou dans le cadre d’une autre application.
Remarque :
Pour en savoir plus sur la création d’applications, suivez ce didacticiel.
Dans ce didacticiel, vous avez créé des modèles de fonction raster pour traiter les couches de critères, synthétisé le risque selon le sous-bassin du bassin versant à l’aide d’un outil d’analyse raster et vous avez partagé la couche obtenue. Vous disposez désormais d’une couche de résultat finale illustrant les zones avec le risque de glissement de terrain le plus élevé que vous pouvez inclure dans d’autres applications ou cartes Web et partager avec votre organisation.
