TauDEM: valores idénticos de acumulación de caudal a lo largo del arroyo

Question

Hace poco descubrí TauDEM y la experiencia en general es fantástica.

Sin embargo, he encontrado algo extraño, cuando se calcula la acumulación de flujo (ESRI ArcGIS término) o D8 área de contribución (TauDEM término). Algunas células vecinas, muy pocos, sin embargo, a lo largo de la línea de flujo de alta almacenar el mismo valor:

lo que, en mi opinión, contradice la definición de acumulación de flujo. Aún más extraño, esto está sucediendo en una sola ruta de flujo, a menos que no he podido encontrar otros duplicados.

Publico esta pregunta porque el autor de software David Tarboton sugiere hacerlo en su sitio . Así que esperemos que pueda comentar esta cuestión (menor) en beneficio de múltiples usuarios.

Answer 1

Lo más probable es que esto se deba al redondeo numérico y a la precisión. Tal y como se implementan actualmente en TauDEM, los resultados de aread8 se guardan en una variable "float" en C++ que tiene unos 7 dígitos decimales de precisión (línea 152 en https://github.com/dtarb/TauDEM/blob/Develop/src/aread8.cpp ). Concretamente, si se añade el número 1 a 26648268 los efectos numéricos no producirán ningún cambio.

En la medida en que esto sea un problema, podría solucionarse cambiando a doble a costa de utilizar más memoria. También tendríamos que asegurarnos de que la rejilla donde se escriben los resultados admite doble. Otra solución podría ser que la versión no ponderada utilizara una representación entera para el resultado del área contribuyente. La decisión (tomada por mí hace mucho tiempo en el momento de la implementación) de utilizar float fue para que el mismo código manejara el área de contribución ponderada y no ponderada, reconociendo que los pesos pueden ser números de coma flotante. El problema reportado es una consecuencia de esta decisión cuando el área contribuyente se hace muy grande y en retrospectiva se podría argumentar que el uso de doble o entero para el área contribuyente no ponderada habría sido mejor.

La representación numérica y la precisión podrían cambiarse si existe una preocupación suficiente de que esto es un problema, o alguien con el tiempo y la energía para trabajar a través del código para cambiarlo hace los cambios y los ofrece de nuevo a TauDEM como una solicitud de pull.