Distinguer les couloirs des fragments dans une carte de l'habitat


12

Je me demandais si quelqu'un avait déjà abordé ce genre de problème:

Je veux délimiter des fragments basés sur des pixels voisins de la même classe (forêt) avec un seuil de distance de 3 pixels (ce qui est biologiquement significatif pour mon espèce d'intérêt).

Ma préoccupation est, voir l'image ci-jointe pour un exemple, que parfois ces fragments sont en fait des couloirs, et souvent les couloirs et les fragments réels sont regroupés en un même fragment en raison de leur proximité.

entrez la description de l'image ici

Je me demande s'il existe un moyen de distinguer les couloirs des fragments en fonction de la forme, du nombre de pixels environnants, etc.?

Par exemple, dans la case suivante, les couloirs possibles sont indiqués par les cases rouges et les fragments par les cases vertes.

J'ai accès à QGIS et R.

entrez la description de l'image ici


pourriez-vous s'il vous plaît ajouter un chiffre avec le type de couloirs que vous souhaitez extraire de l'exemple?
radouxju

Réponses:


4

Avant de commencer une analyse, je vous recommande fortement d'appliquer un filtre à vos données pour nettoyer l'effet "sel et poivre". Tout algorithme aura du mal avec le modèle structurel actuel de vos données. Une simple majorité focale donnerait probablement des résultats indésirables. Une méthode plus robuste applique une approche par tamisage, où une unité de mappage minimale peut être spécifiée. Cela peut être fait via la fonction gdal_sieve.py dans GDAL, la raster > analysis > sievefonction dans QGIS ou la fonction sieve dans ArcGIS Gradient Metrics Toolbox .

Cela ressemble à quelque chose qui pourrait être résolu avec les opérateurs de morphologie mathématique (par exemple, extraire des routes de l'imagerie ). J'imagine qu'un opérateur de dilatation suivi d'un opérateur de fermeture clarifierait les couloirs. Vous pouvez ensuite appliquer un opérateur d'ouverture pour supprimer les couloirs et différencier les rasters pour tirer les couloirs identifiés en tant qu'objets distincts. Ces types de fonctions de décomposition d'images sont quelque peu automatisés dans les logiciels MSPA et GUIDOS mais, encore une fois, seraient notamment affectés par les discontinuités de vos données.

Il existe un plug-in QGIS pour MSPA ainsi que des fonctions disponibles dans GRASS (disponibles via l'interface graphique QGIS). L'un des problèmes avec MSPA et GUIDOS est que la taille de l'image est limitée. Malheureusement, dans le logiciel ESRI, les opérateurs morphologiques ne sont disponibles que dans l' extension ArcScan . Avec quelques fouilles, vous trouverez d'autres options logicielles ainsi que des méthodes pour définir des opérateurs morphologiques via l'algèbre raster avec des matrices de noyau personnalisées.

Une autre approche consisterait à utiliser des méthodes de filtrage de détection de bord, comme un opérateur de noyau Sobal . Il existe une fonction sobal dans ArcGIS Gradient Metrics Toolbox ainsi que dans le package spatialEco R. L'avantage de l'implémentation R est que vous pouvez renvoyer la fonction de gradient de l'opérateur alors que l'implémentation ArcGIS ne renvoie que la fonction de premier ordre (qui peut être tout ce dont vous avez besoin). Je crois que la boîte à outils Orfeo (disponible en tant que module complémentaire QGIS) a une option sobal dans la fonction EdgeExtraction.


peut-être que l'application d'un filtre «saler et poivrer» peut avoir pour effet de supprimer de petites zones d'habitat importantes pour les espèces cibles reliant les fragments (bien que parfois les zones de qualité sous-optimales). Je déciderais d'appliquer un tel filtre en fonction des caractéristiques de dispersion des espèces et de la résolution spatiale des données. Un oiseau (avec une capacité de dispersion généralement plus élevée) peut utiliser ces petites parcelles alors qu'un amphibien (avec une capacité de dispersion comparativement plus faible) ne peut pas. Juste une pensée ..
Kamo

C'est toujours un compromis mais l'unité de cartographie minimale devrait tenir compte de cette hypothèse donnée. Mais, vous devez également garder à l'esprit que l'équilibre est que le modèle présente un certain degré d'incertitude et que vous traitez le résultat comme un absolu. Honnêtement, pour étayer votre affirmation, le paysage doit être traité comme un gradient probabiliste et non comme un processus binomial. Ce type de filtrage est une pratique reconnue depuis longtemps pour représenter des paysages plus fonctionnels. À moins que le modèle n'inclue un terme tenant compte de la structure spatiale dans les estimations, l'incertitude spatiale est une réalité.
Jeffrey Evans

5

Ce n'est pas une solution complète mais, consultez ces outils pour l'analyse de la connectivité (le premier correspond bien à ce que vous recherchez):

Pensez également à former un algorithme d'apprentissage automatique pour classer vos instances (corridor vs fragments). Vous pouvez donner des attributs spatiaux au niveau du patch (par exemple, la taille du patch, le rapport de zone de périmètre, le rapport de cercle) et les entités basées sur la distance (par exemple, la distance aux fragments) à classer. Pour calculer les fonctionnalités de niveau correctif nécessaires à la classification, vous pouvez essayer FragStats ( https://www.umass.edu/landeco/research/fragstats/fragstats.html ).

Vous pouvez également envisager un «système de règles basé sur l'expert» plus simple pour classer chaque instance. Par exemple, les corridors auront un rapport périmètre-surface plus élevé que les fragments d'habitat, etc.

Plus d'informations intéressantes ici pour les analyses de connectivité: http://conservationcorridor.org/corridor-toolbox/programs-and-tools/

Cependant, le fait que certains couloirs ne disposent pas d'une «connectivité pixel complète» sera un problème que vous devrez d'abord régler. Je pense que vous devez définir une sorte de critère de seuil basé sur la distance pour décider si un pixel fait partie d'un couloir ou non.


0

Il semble que cela dépende de la forme de votre fragment. Si la largeur est plus de 2 (ou 3) fois la hauteur (ou l'inverse), vous pouvez l'appeler un couloir?

Êtes-vous arrivé à la délimitation des fragments?

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.