Fonctionnement de l’outil Centile de surface multi-échelle

Disponible avec une licence Spatial Analyst.

L’analyse de la topographie et d’autres surfaces constitue une part importante de nombreuses disciplines, allant de l’hydrologie à l’écologie. Les résultats de ces analyses dépendent souvent de la résolution spatiale des données ou des calculs pour une caractéristique topographique en particulier. Cette dépendance a généré une généralisation de l’approche multi-échelle des analyses, où les calculs sont réalisés pour plusieurs résolutions spatiales. Ces approches multi-échelles peuvent servir à déterminer l’échelle optimale pour caractériser une topographie et mesurer la façon dont les paramètres réagissent aux changements d’échelle.

L’outil Multiscale Surface Percentile (Centile de surface multi-échelle) calcule le centile le plus extrême sur une plage d’échelles spatiales (voisinages dont la taille varie). Le centile le plus éloigné de 50 (par exemple, les valeurs les plus proches de 0 ou de 100) est considéré comme la valeur la plus extrême pour une cellule donnée. Les sorties de cet outil identifient ce centile pour une cellule et son échelle.

Les sorties peuvent servir à interpréter une entité sur un raster de surface en entrée et les échelles associées. L’image ci-dessous montre les résultats de deux échelles différentes pour la même surface en entrée. Une sortie a utilisé une échelle de 9 cellules par 9 cellules, tandis que l’autre sortie a utilisé une échelle de 29 cellules par 29 cellules. Ici, l’échelle inférieure est plus sensible aux variations locales dans le paysage et capture des entités de surface plus petites. De l’autre côté, l’échelle supérieure montre moins de détails en capturant uniquement les entités de surface de plus grande taille.

Illustration d’un centile en sortie avec la même étendue pour deux échelles différentes : une petite échelle (deuxième image) et une plus grande échelle (troisième image).

Mode de calcul de la valeur du centile le plus extrême

Les étapes suivantes fournissent une vue d’ensemble des processus internes utilisés par l’outil :

1. Pour l’analyse, les échelles sont définies à l’aide des paramètres Minimum Neighborhood Distance (Distance de voisinage minimale), Maximum Neighborhood Distance (Distance de voisinage maximale), Base Distance Increment (Incrément de distance de base) et Nonlinearity Factor (Facteur de non-linéarité). Les unités de ces paramètres sont déterminées par le paramètre Distance Units (Unités de distance).

2. Pour chaque cellule, un centile est calculé à chaque échelle identifiée.

3. Les centiles calculés sont comparés sur les différentes échelles. Le centile le plus éloigné de 50 est considéré comme le centile le plus extrême pour une cellule donnée.

Chacune de ces étapes est expliquée plus en détail dans les sections ci-dessous.

Mode d’identification des échelles à analyser

Pour l’analyse, les échelles sont déterminées à l’aide des paramètres facultatifs de l’outil Centile de surface multi-échelle. Les paramètres Minimum Neighborhood Distance (Distance de voisinage minimale) et Maximm Neighborhood Distance (Distance de voisinage maximale) définissent les échelles minimale et maximale pour l’analyse. Les paramètres Base Distance Increment (Incrément de distance de base) et Nonlinearity Factor (Facteur de non-linéarité) déterminent l’augmentation de la distance de voisinage entre la valeur minimale et la valeur maximale.

Chaque échelle est représentée par une valeur de distance de voisinage. L’analyse est effectuée pour plusieurs distances de voisinage en fonction des paramètres en entrée.

Pour une cellule cible donnée, la distance de voisinage est mesurée à partir du centre de la cellule cible vers l’extérieur, ce qui crée un carré de cellules autour de la cellule cible. Par exemple, une distance de voisinage de 30 mètres pour un raster de surface en entrée avec une taille de cellule de 10 mètres génère un voisinage de 7 cellules par 7 cellules, comme illustré dans la figure ci-dessous. Cette valeur, 30 mètres, correspond à l’une des échelles pour lesquelles le centile le plus extrême est calculé.

Illustration de la relation entre la distance de voisinage (ligne orange) et le nombre de cellules dans la fenêtre mobile. Pour une taille de cellule de 10 mètres, une distance de voisinage de 10 mètres utilisera une fenêtre de 3 cellules par 3 cellules (par défaut), une distance de voisinage de 20 mètres utilisera une fenêtre de 5 cellules par 5 cellules et une distance de voisinage de 30 mètres utilisera une fenêtre de 7 cellules par 7 cellules.

La plus petite distance de voisinage autorisée est égale à la taille de cellule du raster en entrée. Il s’agit d’une valeur de 1 cellule, qui crée un voisinage de cellules de 3 par 3. Dans l’exemple ci-dessus, cette distance de voisinage minimale est de 5 mètres.

La distance de voisinage ne peut pas être supérieure au raster de surface en entrée.

Si une distance de voisinage qui ne génère pas un intervalle de la taille de cellule est spécifiée, l’outil arrondit la distance à l’intervalle supérieur suivant de la taille de cellule. Par exemple, dans l’illustration ci-dessus, si une distance de voisinage de 25 mètres est spécifiée, elle est arrondie à l’intervalle supérieur suivant de la taille de cellule, à savoir 30 mètres.

Les calculs débutent avec la valeur du paramètre Minimum Neighborhood Distance (Distance de voisinage minimale). Ensuite, chaque distance de voisinage suivante est calculée.

Voici l’expression à utiliser pour calculer les distances de voisinage suivantes :

\(n_i = n_o + [\Delta n * (i - n_o)]{^p}\)

• où :

  \(n_i\) est la distance de voisinage pour l’étape \(i\).

  \(n_o\) est la distance de voisinage minimale.

  \(\Delta n\) est l’incrément de distance de base.

  \(i\) est l’étape pour laquelle la distance de voisinage est calculée (où la première étape a la valeur \(1 + n_o\)).

  \(p\) est le facteur de non-linéarité.

Chaque nouvelle distance de voisinage identifiée est vérifiée afin de déterminer si elle est inférieure ou égale à la valeur du paramètre Maximum Neighborhood Distance (Distance de voisinage maximale). Si la nouvelle valeur de distance est inférieure ou égale à la valeur maximale, les calculs de la distance de voisinage se poursuivent. Si la nouvelle valeur est supérieure à la valeur maximale, toutes les distances de voisinage ont été identifiées et les calculs du centile débutent.

Reportez-vous à la section Mode de calcul du centile ci-dessous pour plus d’informations sur cette partie de l’analyse.

Impact du facteur de non-linéarité sur les distances de voisinage

Le paramètre Nonlinearity Factor (Facteur de non-linéarité) contrôle le taux d’augmentation de la distance de voisinage. La valeur par défaut est définie sur 1, ce qui génère une augmentation linéaire de la distance de voisinage. Cela signifie que les incréments entre les distances de voisinage seront égaux à la valeur du paramètre Base Distance Increment (Incrément de distance de base).

Lorsque la valeur du paramètre Nonlinearity Factor (Facteur de non-linéarité) dépasse 1, les incréments entre les distances de voisinage changent après le premier. Le premier incrément est égal à la valeur Base Distance Increment (Incrément de distance de base), mais la taille de tous les incréments suivants augmente progressivement.

Lorsque le paramètre Nonlinearity Factor (Facteur de non-linéarité) est définie sur des valeurs supérieures à 1, l’incrément de distance entre les distances de voisinage suivantes après le premier incrément augmente progressivement. Autre conséquence : pour les mêmes valeurs de voisinage minimale et maximale, un facteur de non-linéarité supérieur entraîne dans l’ensemble un nombre réduit de distances de voisinage.

La figure ci-dessous illustre l’effet de trois réglages différents du paramètre Nonlinearity Factor (Facteur de non-linéarité). Dans cet exemple, les paramètres sont 1,0, 1,5 et 2,0. Pour chacun de ces paramètres, les valeurs des autres paramètres restent identiques. La valeur du paramètre Minimum Neighborhood Distance (Distance de voisinage minimale) est 1, la valeur du paramètre Maximum Neighborhood Distance (Distance de voisinage maximale) est 10 et la valeur du paramètre Base Distance Increment (Incrément de distance de base) est 1.

Pour le premier incrément, la distance de voisinage pour les trois paramètres du facteur de non-linéarité a la même valeur, c’est-à-dire 2 cellules. Ensuite, les valeurs pour la distance de voisinage commencent à varier. Les incréments augmentent progressivement pour un facteur de non-linéarité de 1,5 et encore plus pour un facteur de 2,0.

Lorsque la valeur du paramètre Nonlinearity Factor (Facteur de non-linéarité) est 1,0, le nombre total d’incréments est de 9 et chaque incrément est plus élevé que le précédent, de façon linéaire. Avec un facteur d’une valeur de 1,5, on ne compte que 4 incréments, et 3 incréments pour une valeur de 2,0.

L’augmentation du facteur de non-linéarité au-dessus de 1,0 entraîne l’augmentation progressive des tailles d’incrément et par conséquent, l’augmentation plus rapide des distances de voisinage. De plus, le nombre d’échelles pour les mêmes valeurs de distance de voisinage minimale et de distance de voisinage maximale est réduit lorsqu’une valeur de facteur de non-linéarité plus élevée est utilisée.

Le paramètre Nonlinearity Factor (Facteur de non-linéarité) vous permet de personnaliser la densité d’échantillonnage des échelles. Pour une surface d’élévation, un centile d’élévation est souvent plus sensible à la taille de voisinage à de plus petites échelles et moins sensible à de plus grandes échelles. Une valeur de paramètre Nonlinearity Factor (Facteur de non-linéarité) supérieure à 1,0 permet d’augmenter la densité d’échantillonnage pour les plus petites échelles et de la réduire pour les plus grandes échelles. Dans ce cas, il peut toutefois être nécessaire d’augmenter la valeur du paramètre Maximum Neighborhood Distance (Distance de voisinage maximale) pour obtenir le nombre souhaité d’incréments. Dans la plupart des cas, une valeur comprise entre 1,0 et 2,0 est utilisée pour le facteur de non-linéarité.

Mode de calcul du centile

Les centiles sont une mesure statistique qui indique le pourcentage de valeurs dans un jeu de données qui sont inférieures à une valeur donnée. Par exemple, le 80^e centile est une valeur dans laquelle vous trouverez que 80 % des valeurs dans le jeu de données sont inférieures, tandis que les 20 % restantes sont supérieures.

Pour chaque distance de voisinage identifiée pour le calcul et chaque cellule du raster de surface en entrée, l’outil Centile de surface multi-échelle calcule le centile. Les valeurs de centile les plus extrêmes sont identifiées et enregistrées dans la valeur du paramètre Output Percentile Raster (Raster de centile en sortie). Les échelles auxquelles ces centiles ont été trouvés sont enregistrées comme valeurs de cellule dans le paramètre Output Scale Raster (Raster d’échelle en sortie).

L’équation permettant de calculer le centile pour chaque cellule est la suivante :

\(Percentile = count{\large_i}{\large_\in}{\large_C}(z_i < z_0) * (100/n_C)\)

• où :

  \(C\) est le voisinage identifié pour le traitement.

  \(count{\large_i}{\large_\in}{\large_C}\) est le nombre de cellules dans le voisinage \(C\) où (\(z_i < z_0\)) est vrai.

  \(z_i\) est la valeur de la cellule \(i\) dans le voisinage \(C\).

  \(z_0\) est la valeur de la cellule centrale dans le voisinage \(C\).

  \(n_C\) est le nombre de cellules contenues dans le voisinage \(C\).

Voici l’illustration d’un exemple de ce calcul avec un voisinage de 3 cellules par 3 cellules.

L’échelle est de 3 cellules par 3 cellules et comporte 9 cellules au total. À gauche, la valeur de la cellule de traitement et les 8 valeurs de cellule environnantes (en gris) sont identifiées. À droite, les cellules dont les valeurs sont inférieures à celle de la cellule centrale (en vert) et les valeurs qui sont égales ou supérieures à la cellule centrale (en brun) sont affichées.

L’application de la formule ci-dessus à cet exemple génère le résultat suivant :

Percentile = (Count of cell values less than center cell value)
           * 100 / (Count of cells in neighborhood)
        = 5 * 100 / 9
        = 55.5556

Si une valeur de centile est trouvée, elle est comparée au centile le plus extrême précédemment identifié. Si la nouvelle valeur est plus extrême, avec un centile plus éloigné de 50, elle est enregistrée pour cette localisation dans la valeur du paramètre Output Percentile raster (Raster de centile en sortie). La valeur d’échelle est enregistrée dans le paramètre Output Scale Raster (Raster d’échelle en sortie).

Utilisation d’un GPU

Cet outil propose des performances accrues si un certain matériel GPU est installé sur votre système. Reportez-vous à la rubrique Traitement GPU avec Spatial Analyst pour des détails sur la prise en charge, la configuration et l’activation de cette fonctionnalité.

Bibliographie

Newman, Daniel R., John B. Lindsay et Jaclyn Mary Helen Cockburn. 2018. "Evaluating metrics of local topographic position for multiscale geomorphometric analysis." Geomorphology 312, 40 à 50. https://doi.org/10.1016/j.geomorph.2018.04.003

Rubriques connexes

• Vue d’ensemble du jeu d’outils Surface

• Centile de surface multi-échelle

• Fonctionnement de l’outil Multiscale Surface Percentile (Différence de surface multi-échelle)