Antes de comenzar cualquier análisis, recomendaría aplicar un filtro a sus datos para limpiar el efecto "sal y pimienta". Cualquier algoritmo tendrá problemas con el patrón estructural actual de sus datos. Una mayoría focal simple probablemente arrojaría resultados indeseables. Un método más robusto es aplicar un enfoque de tamiz, donde se puede especificar una unidad de mapeo mínimo. Esto se puede hacer a través de la función gdal_sieve.py en GDAL, la raster > analysis > sievefunción en QGIS o la función de tamizado en el ArcGIS Gradient Metrics Toolbox .
Esto parece algo que podría abordarse con los operadores de Morfología matemática (por ejemplo, extraer caminos de las imágenes ). Me imagino que un operador de dilatación seguido de un operador de cierre aclararía los pasillos. A continuación, puede aplicar un operador de apertura para eliminar los corredores y diferenciar los rásteres para extraer los corredores identificados como objetos separados. Estos tipos de funciones de descomposición de imágenes están algo automatizadas en el software MSPA y GUIDOS pero, una vez más, se verían notablemente afectadas por las discontinuidades en sus datos.
Hay un complemento QGIS para MSPA , así como funciones disponibles en GRASS (disponible a través de la GUI QGIS). Uno de los problemas con MSPA y GUIDOS es que tiene un tamaño limitado de imagen. Desafortunadamente, en el software ESRI, los operadores morfológicos solo están disponibles en la extensión ArcScan . Con un poco de excavación, encontrará otras opciones de software, así como métodos para definir operadores morfológicos a través del álgebra ráster con matrices de kernel personalizadas.
Otro enfoque serían los métodos de filtrado de detección de bordes, como un operador de kernel Sobal . Hay una función sobal en ArcGIS Gradient Metrics Toolbox , así como en el paquete spaceialEco R. La ventaja de la implementación de R es que puede devolver la función de gradiente del operador, mientras que la implementación de ArcGIS solo devuelve la función de primer orden (que puede ser todo lo que necesita). Creo que la caja de herramientas Orfeo (disponible como un complemento QGIS) tiene una opción simple en la función EdgeExtraction.
