Problemas con la vectorización de la acumulación de flujo diluido mediante r.to.vect en QGIS

Question

Estoy tratando de producir una red fluvial shapefile a partir de un raster de acumulación de flujo utilizando el programa r.to.vect y he superado algunas dificultades relacionadas con la profundidad de bits y la trama adelgazada.

n.b. No tengo el DEM original del que se derivó el raster de acumulación de flujo, por lo que no puedo basarme en el Channel Network and Drainage Basins función.

Mi objetivo es obtener una red fluvial con el valor 'upcells' del antiguo conjunto de datos de redes fluviales HydroSHEDS pero mi resultado presenta algunos problemas que no he encontrado la forma de resolver de forma automática:

Aquí tenemos una trama booleana adelgazada y el vector lineal resultante.

1. El primer problema es la ambigüedad entre celdas adyacentes (círculo rojo) que hace que la r.to.vect crear líneas redundantes;
2. El segundo problema se produce cuando hay un ángulo recto entre celdas, pero sólo en la dirección horizontal/vertical (círculo azul). En este caso, el r.to.vect divide las líneas cercanas a estos ángulos rectos, como vemos en la etiqueta 'id' (el elemento 758 se divide del elemento 755, del mismo modo el elemento 759 se divide del elemento 764, pero esto no ocurre para el elemento 764).

La segunda cuestión afecta al siguiente paso de mi flujo de trabajo cuando voy a utilizar v.rast.stats para asignar a cada característica el valor máximo del ráster del mapa de acumulación de flujo. Este es el resultado:

Los vectores lineales se gradúan ahora en función del valor de acumulación de caudal (etiquetas rojas), pero el v.rast.stats no se pueden calcular estadísticas para algunos de estos segmentos cortos generados por la segunda emisión (elementos 759 y 794; el elemento 799 no es relevante porque es redundante)

He probado algunas configuraciones opcionales del r.to.vect pero no parece capaz de gestionar estos dos problemas.

¿Hay alguna forma de superarlas de manera automática/semiautomática?

Answer 1

Ya que has etiquetado la pregunta también con GRASS, me gustaría señalar que ninguno de estos problemas se produce cuando se delinean redes de arroyos utilizando la función r.watershed en GRASS. Ese módulo utiliza la dirección de flujo múltiple (MFD) para determinar la dirección del flujo, evitando así las divisiones en el círculo rojo. Y nunca he visto que un curso de agua se divida por la mitad como se muestra en el círculo azul.

A falta del DEM original, se podría intentar tamponar los flujos a una distancia suficiente para eliminar la ambigüedad "roja", y luego "desamponar" (tamponar con negativo la misma distancia) para restaurar las líneas.