Créer une carte de gravité des dommages causés par le feu
L’analyse des données raster et d’image donne des informations précieuses pour la prise de décision dans le cadre des applications de gestion des urgences. L’imagerie représente les informations actuelles d’un emplacement ou d’une situation et peut être analysée parallèlement à des images historiques ou à d’autres informations spatiales. Le traitement des données raster peut s’avérer ardu car les fichiers sont habituellement complexes et volumineux. Par ailleurs, il est souvent nécessaire de partager et d’analyser à l’échelle de l’entreprise des données provenant de différentes organisations et de différentes sources. La conception élastique du traitement distribué de l’analyse raster offre un moyen efficace de partager et de traiter rapidement de gros volumes de données raster afin de prendre en charge les applications soumises à des contraintes de temps. Une fois votre chaîne de traitement optimisée, vous pouvez l’enregistrer, la partager avec les membres de votre organisation et exécuter le traitement dans votre déploiement d’analyse raster depuis ArcGIS Pro ou Map Viewer Classic dans ArcGIS Enterprise.
Dans ce module, vous allez utiliser des chaînes de fonctions raster pour dériver une carte déterminant la gravité des dommages causés par le feu, cela étant essentiel pour savoir quelles zones présentent le plus de dommages par le feu.
Remarque :
Ce didacticiel utilise Map Viewer Classic. Map Viewer Classic est le prédécesseur de Map Viewer. Certaines fonctionnalités ne sont pas encore disponibles dans Map Viewer. Il est recommandé d’utiliser Map Viewer Classic pour le processus suivant jusqu’à ce que ce dernier soit pris en charge dans une version ultérieure de Map Viewer.
Télécharger les données et définir le portail actif
En premier lieu, vous allez télécharger les données du didacticiel, puis allez définir le portail actif dans ArcGIS Pro. Cette étape de la configuration transmet les informations de licence de votre compte à ArcGIS Pro et vous permet de partager les données traitées dans ArcGIS Enterprise.
- 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 Settings (Paramètres).
- Cliquez sur Portals (Portails).
- Si votre portail ne figure pas dans la liste, cliquez sur Add Portal (Ajouter un portail) et saisissez l’URL du portail.
Lorsqu’il est détecté, le portail est ajouté à la liste.
- Cliquez avec le bouton droit sur votre portail et choisissez Sign in (Se connecter).
Assurez-vous de vous connecter avec un compte qui dispose des autorisations permettant d’utiliser la fonctionnalité d’analyse raster.
- Entrez votre nom d’utilisateur et votre mot de passe, puis cliquez sur Sign in (Se connecter).
- Cliquez avec le bouton droit sur votre portail et sélectionnez Set As Active Portal (Définir comme portail actif).
Une coche verte apparaît en regard de la connexion à votre portail.
- Cliquez sur la flèche de retour.
Maintenant que vous êtes connecté à Portal, vous allez ouvrir un modèle de projet pour commencer l’analyse.
- Cliquez sur Start with another template (Démarrer avec un autre modèle).
- Accédez à C:\LandslideData sur votre machine Enterprise. Choisissez le fichier du modèle de projet, Landslide_Risk_Project.aptx, puis cliquez sur OK.
Conseil :
Si vous ne parvenez pas à trouver ce fichier, vous pouvez télécharger et décompresser le fichier LandslideData.zip sur votre machine locale.
- 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).
Vous avez téléchargé les données, défini le portail actif, ouvert un modèle de projet et créé une connexion aux dossiers. À présent, vous pouvez procéder à 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 les imageries Landsat 8 acquises avant et après les incendies autour de Santa Rosa, en Californie, en 2017.
Pour commencer, vous allez ajouter les imageries 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 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 sur l’unité C de votre machine Enterprise. Sélectionnez-le et cliquez sur OK (OK).
- 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 srband5.
- Pour le vert, sélectionnez srband4 et pour le bleu, srband3.
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 srband5.
- Associez la pastille de couleur verte à srband4.
- Associez la pastille de couleur bleue à srband3.
À présent, vous allez utiliser l’outil 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).
- 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 réglé sur Raster Layer (Couche raster), puis cliquez sur Create new layer (Créer une 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.
- Dans la barre d’outils Accès rapide, cliquez sur le bouton Save (Enregistrer).
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 Remap (Classification) à 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.
- Pour Description (Description), saisissez Modèle de fonction raster permettant de déduire la pente d’un MNA en entrée.
- Cliquez sur OK.
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).
Conseil :
Il peut être nécessaire d’effacer le contenu de la barre de recherche pour pouvoir voir le modèle de fonction 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. Veillez à ce que Output Layer Type (Type de couche en sortie) ait pour valeur 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 plusieurs modèles de fonction raster et critères 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 qui sera traitée dans votre déploiement ArcGIS Enterprise.
- 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).
Le modèle de fonction raster de risque de glissement de terrain est complet.
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 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.
- Pour Description (Description), saisissez Modèle de fonction raster permettant de calculer le risque de glissement de terrain en fonction de la gravité des incendies, 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), choisissez Web Image Layer (Couche Web d’imagerie), puis cliquez sur Next (Suivant).
Remarque :
Si vous ne voyez pas d’option permettant d’associer votre sortie à une couche Web d’imagerie, demandez à votre administrateur de s’assurer que le serveur d’analyse raster est configuré correctement et est opérationnel.
La fenêtre Output Generation (Génération en sortie), dans laquelle vous pouvez définir les propriétés de la couche Web que vous allez créer dans Enterprise, s’ouvre.
- Pour Name (Nom), entrez Landslide_Risk (Risque_glissement_de_terrain). Dans Description, saisissez This web layer shows a landslide risk estimate for Sonoma County, CA (Cette couche web montre une estimation des risques de glissement de terrain pour le comté de Sonoma, en Californie).
- Pour Tags (Balises), entrez Sonoma County, landslide, risk, wildfire (comté de Sonoma, glissement de terrain, risque, incendie). Appuyez sur la touche Entrée.
Quatre balises distinctes sont ajoutées.
- Cliquez sur Run (Exécuter).
La chaîne de traitement est soumise à votre déploiement Enterprise en vue d’un traitement distribué. Un message vous informe que le processus a été soumis.
- Cliquez sur la notification Raster Functions (Fonctions raster) pour afficher le statut de votre modèle de fonction raster.
La fenêtre History (Historique), dans laquelle l’onglet Raster Functions (Fonctions raster) est sélectionné, s’ouvre. Toutes les fonctions raster que vous avez utilisées précédemment sont répertoriées ici.
Remarque :
Si la notification Raster Functions (Fonctions raster) disparaît, une autre méthode permet d’accéder au statut du traitement raster. Sur le ruban, cliquez sur l’onglet Analysis (Analyse). Dans le groupe Géotraitement, cliquez sur Historique. Dans la fenêtre History (Historique), cliquez sur l’onglet Raster Functions (Fonctions raster).
Vous verrez une coche verte en regard de la fonction raster une fois que le processus a abouti.
Vous avez créé la couche des risques de glissement de terrain en utilisant des fonctions raster, puis l’avez partagée en tant que couche d’imagerie Web sur votre portail. Ensuite, vous allez utiliser cette couche dans Map Viewer Classic pour procéder à une analyse raster.
Synthétiser le risque de glissement de terrain selon le bassin du bassin versant
Partagez la couche de bassins sur votre portail
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. Vous allez d’abord publier la couche des bassins versants dans votre Portal. Vous allez ensuite utiliser le jeu de données des risques de glissement de terrain généré dans la section précédente pour synthétiser le risque par sous-bassin.
- Dans la fenêtre Catalog (Catalogue), cliquez sur l’onglet Project (Projet), puis accédez à la connexion au dossier LandslideData.
- Cliquez avec le bouton droit sur Basins.tif et choisissez Share As Web Layer (Partager en tant que couche web).
La fenêtre Share As Web Layer (Partager en tant que couche Web) s’ouvre. Dans cette fenêtre, vous pouvez nommer et configurer ce jeu de données raster lorsque vous le publiez dans votre Portal.
- Sous Item Details (Détails de l’élément), définissez les valeurs suivantes :
- Pour Name (Nom), entrez Basins (Bassins).
- Pour Summary (Résumé), saisissez Sous-bassin du bassin versant dans les comtés Napa et Sonoma, en Californie.
- Pour Tags (Balises), entrez watershed, basin, California (bassin versant, bassin, Californie), puis appuyez sur Entrée.
Vous allez maintenant vous assurer que ce jeu de données raster est référencé depuis le dossier LandslideData et n’est pas enregistré en tant que copie dans le Portal.
- Sous Layer and Data Type (Type de données et de couche), assurez-vous que l’option Reference registered data (Référencer les données inscrites) est sélectionnée.
- Sous Location (Localisation), pour Server and Folder (Serveur et dossier), choisissez l’emplacement de stockage de la couche Basins (Bassins).
Remarque :
L’emplacement de votre serveur dépend de votre configuration et peut être différent de l’exemple présenté dans l’image.
- Cliquez sur Analyze (Analyser).
Il est fort probable qu’un message d’erreur s’affiche et indique que la source de données n’est pas inscrite auprès du serveur. C’est normal. En l’absence d’erreur ou d’avertissement, vous pouvez passer à l’étape 12, car un datastore est déjà inscrit.
Remarque :
Lorsque vous publiez des services Web dans Enterprise et choisissez de référencer des données inscrites, la source de données doit être inscrite auprès de ArcGIS Server. Cette inscription permet au serveur d’accéder à vos données et de les utiliser comme source pour les couches Web. La création d’un data store est essentielle pour que le serveur puisse accéder à vos données. Un data store peut correspondre à n’importe quel emplacement (bases de données d’entreprise, dossiers, stockages Cloud ou bases de données NoSQL) hébergeant les données que vous voulez utiliser. Une fois les données inscrites auprès du serveur, le service Web publié établit une connexion directe à la source de données. Cette connexion garantit que les services Web référencent les données du data store sans les dupliquer.
En règle générale, vous pouvez passer à la publication lorsqu’un avertissement s’affiche, auquel cas, il est recommandé d’inscrire les données auprès du serveur pour être sûr de partager une référence aux données et non pas une copie.
- Cliquez avec le bouton droit de la souris sur l’avertissement, puis sélectionnez Register Data Source With Server (Inscrire la source de données auprès du serveur).
- Dans la fenêtre Add data store (Ajouter un data store), pour Title (Titre), saisissez BasinsDS.
- Pour Tags (Balises), saisissez basins (bassins) et data store.
- Pour Select the servers to which you want to add your data store (Sélectionnez les serveurs auxquels vous souhaitez ajouter votre data store), cochez la case correspondant à votre serveur, puis cliquez sur Validate (Valider).
Remarque :
Pour en savoir plus sur l’inscription et la gestion des datas stores, reportez-vous à la rubrique Gérer les data stores inscrits.
- Cliquez sur Create (Créer).
Dès que vous avez créé le data store, les données sont inscrites auprès du serveur et une coche verte apparaît en regard de l’avertissement. Vous pouvez à présent publier la couche Web.
- Cliquez sur Publish (Publier).
- Cliquez sur le lien Manage the web layer (Gérer la couche Web).
La page des détails de l’élément pour votre couche Basins (Bassins) apparaît.
- Sur le ruban, cliquez sur l’onglet Content (Contenu).
Lors du partage de la couche Web, vous avez créé deux éléments dans votre portail : une couche d’imagerie et un data store. Le data store a été créé lors de l’inscription de la source de données auprès du serveur et était exigé pour le partage sur le portail.
Effectuer une analyse raster
Vous avez partagé la couche Basins (Bassins) sur votre portail et pouvez à présent l’inclure dans votre analyse. Vous allez utiliser un outil d’analyse raster appelé Summarize Within (Synthétiser - À l’intérieur) pour synthétiser le risque de glissement de terrain en fonction des caractéristiques du sous-bassin du bassin versant.
- Cliquez sur la couche Basins (Bassins) pour retourner à la page de détails de ses éléments.
- Cliquez sur la flèche orientée vers le bas en regard de Open in Map Viewer (Ouvrir dans Map Viewer) et sélectionnez Open in Map Viewer Classic (Ouvrir dans Map Viewer Classic).
Map Viewer Classic s’ouvre.
À présent, vous allez ajouter le jeu de données Landslide_Risk (Risque_glissement_de_terrain) à la carte.
- Sur le ruban, cliquez sur Add (Ajouter) et choisissez Search for Layers (Rechercher des couches).
- Vérifiez que la recherche est définie sur My Content (Mon contenu). Recherchez la couche Landslide_Risk (Risque_glissement_de_terrain), puis cliquez sur le bouton Add (Ajouter).
La couche Landslide_Risk (Risque_glissement_de_terrain) est ajoutée à la carte.
Remarque :
Si la couche de risques de glissement de terrain ne peut pas être chargée sur la carte, vous devrez peut-être attendre la fin du traitement pendant quelques minutes, puis recommencer.
Maintenant que vous disposez d’une couche pour les bassins et d’une couche pour les risques, vous allez calculer le risque moyen par bassin à l’aide de l’outil Synthétiser un raster à l’intérieur.
- Sur le ruban, cliquez sur Analysis (Analyse) et choisissez Raster Analysis (Analyse raster).
- Développez Summarize Data (Synthétiser les données) et cliquez sur l’outil Synthétiser un raster à l’intérieur.
- Pour l’outil Summarize Raster Within (Synthétiser un raster à l’intérieur), renseignez les paramètres suivants :
- Pour Choose area layer to summarize a raster layer within defined boundaries (Choisir une couche surfacique pour synthétiser une couche raster au sein des limites définies), choisissez Basins (Bassins).
- Pour Choose field to define the boundaries (Choisir un champ pour définir les limites), choisissez Value (Valeur).
- Pour Choose raster layer to summarize (Choisir la couche raster à synthétiser), choisissez Landslide_Risk (Risque_glissement_de_terrain).
- Pour Choose statistic to calculate (Choisir les statistiques à calculer), choisissez Average (Moyenne).
- Pour Ignore missing values in calculation (Ignorer les valeurs manquantes dans le calcul), sélectionnez la case à cocher Ignore (Ignorer).
- Pour Process as multidimensional (Traiter comme multidimensionnel), désélectionnez la case à cocher All slices (Toutes les tranches).
- Pour Result layer name (Nom de la couche de résultat), entrez Risk per Basin (Risque par bassin).
- Désélectionnez la case à cocher Use current map extent (Utiliser l’étendue actuelle de la carte).
- Cliquez sur Run Analysis (Exécuter l’analyse).
Une fois l’analyse terminée, un jeu de données illustrant le risque moyen par bassin versant est ajouté à la carte.
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. Cependant, il est difficile de distinguer les résultats actuels avec cette symbologie. Après avoir enregistré la carte, vous allez mettre à jour la symbologie. Ainsi, il sera plus facile de comprendre les résultats de votre analyse.
- Sur le ruban, cliquez sur Save (Enregistrer) et choisissez Save As(Enregistrer sous).
La fenêtre Save Map (Enregistrer la carte) apparaît.
- Dans la fenêtre Save Map (Enregistrer la carte), définissez les valeurs suivantes :
- Pour Title (Titre), entrez Landslide Risk per Basin (Sonoma and Napa) (Risque de glissement de terrain par bassin [Sonoma et Napa]).
- Pour Tags (Balises), entrez landslide, risk, wildfire (glissement de terrain, risque, incendie), puis appuyez sur Entrée.
- Pour Summary (Résumé), saisissez Risque de glissement de terrain par bassin dans les comtés de Sonoma et Napa, en fonction de la gravité des incendies, de la pente et de l’occupation du sol.
- Cliquez sur Enregistrer la carte.
Votre carte Web est désormais enregistrée dans votre contenu et vous pouvez y accéder pour une analyse plus approfondie ou l’inclure dans d’autres applications.
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 des couches de critères, partagé des couches d’imagerie sur votre Portal, puis synthétisé le risque selon le sous-bassin du bassin versant à l’aide d’un outil d’analyse raster dans Map Viewer Classic. 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.