Cómo funciona De topo a ráster

La herramienta De topo a ráster es un método de interpolación diseñado específicamente para crear modelos digitales de elevación (DEM) hidrológicamente correctos. Está basado en el programa ANUDEM desarrollado por Michael Hutchinson (1988, 1989, 1996, 2000, 2011). Vea Hutchinson and Dowling (1991), ANU Fenner School of Environment and Society y Geoscience Australia (2008) para obtener aplicaciones de ANUDEM para producción de DEM en todo el continente. Las aplicaciones de DEM para modelado ambiental se tratan en Hutchinson and Gallant (2000) y Hutchinson (2008). Otros desarrollos de ANUDEM se tratan en Hutchinson et al. (2009, 2011). La versión actual de ANUDEM utilizada en ArcGIS es 5.3.

De topo a ráster interpola los valores de elevación para un ráster mientras se imponen restricciones que aseguran:

Como tal, es el único interpolador de ArcGIS diseñado específicamente para trabajar de forma inteligente con entradas de curvas de nivel.

La herramienta De topo a ráster por archivo es útil para ejecutar la herramienta De topo a ráster varias veces, ya que es a menudo más fácil cambiar una entrada única en el archivo de parámetros y volver a ejecutar la herramienta que volver a completar el cuadro de diálogo de la herramienta todas las veces.

El proceso de interpolación

El procedimiento de interpolación se diseñó para aprovechar los tipos de datos de entrada comúnmente disponibles y las características conocidas de las superficies de elevación. Este método utiliza una técnica de interpolación de diferencia finita iterativa. Se optimiza para tener la eficacia computacional de los métodos de interpolación local, como la interpolación de distancia inversa ponderada (IDW), sin perder la continuidad de la superficie de los métodos de interpolación global, como Kriging y Spline. Es, esencialmente, una técnica de spline de lámina delgada discretizada (Wahba, 1990) a la cual se le modificó la penalización por rugosidad para permitir que el DEM ajustado siga los cambios abruptos en el terreno, como arroyos, crestas y acantilados.

El agua es la principal fuerza erosiva que determina la forma general de la mayoría de los paisajes. Es por esto que la mayoría de los paisajes tienen varias cumbres (máximos locales) y pocos sumideros (mínimos locales), lo que resulta en un patrón de drenaje conectado. De topo a ráster utiliza este conocimiento de las superficies e impone restricciones al proceso de interpolación que resulta en una estructura de drenaje conectado y la correcta representación de crestas y arroyos. Esta condición de drenaje impuesto produce superficies con mayor precisión y menos datos de entrada. La cantidad de datos de entrada puede ser de hasta un orden de magnitud menor que lo que se requiere normalmente para describir de manera adecuada una superficie con curvas de nivel digitalizadas, lo que minimiza aún más el coste de obtener DEM confiables. La condición de drenaje global también elimina prácticamente cualquier necesidad de edición o postproceso para quitar los sumideros espurios de la superficie generada.

El programa actúa de manera conservadora en la eliminación de sumideros y no impone condiciones de drenaje en las ubicaciones que podrían contradecir los datos de elevación de entrada. Estas ubicaciones suelen aparecer en el archivo de diagnóstico como sumideros. Utilice esta información para corregir los errores en los datos, en particular cuando procesa datasets de gran tamaño.

El proceso de aplicación de drenaje

El propósito del proceso de aplicación de drenaje es quitar todos los puntos de sumideros en el DEM de salida que no se identificaron como sumideros en el dataset de entidades de sumidero de entrada. El programa presupone que todos los sumideros no identificados son errores, ya que, por lo general, es poco común encontrarlos en paisajes naturales (Goodchild y Mark, 1987).

El algoritmo de aplicación de drenaje intenta eliminar los sumideros espurios modificando el DEM, infiriendo las líneas de drenaje a través del punto de collado más bajo en el área de drenaje que circunda cada sumidero espurio. No intenta eliminar los sumideros reales suministrados por la función Sumidero. Como la eliminación de sumideros está sujeta a la tolerancia de elevación, el programa es conservador al intentar eliminar los sumideros espurios. Es decir, no elimina los sumideros espurios que puedan contradecir los datos de elevación de entrada por más del valor de Tolerancia 1.

La aplicación de drenaje también se puede complementar con la incorporación de los datos de línea de corriente. Esto es útil cuando se requiere una ubicación más precisa de los arroyos. Distributarios de arroyos se modelan permitiendo a cada celda tener hasta dos direcciones de aguas abajo.

Es posible desactivar la aplicación de drenaje, en cuyo caso, se ignora el proceso de eliminación de sumideros. Esto puede ser útil si tiene datos de curvas de nivel de algo diferente a la elevación (por ejemplo, la temperatura) para los que desea crear una superficie.

Uso de los datos de curvas de nivel

Las curvas de nivel eran originalmente el método más común para el almacenamiento y la presentación de la información de elevación. Lamentablemente, este método también es el más difícil de utilizar de manera adecuada con las técnicas de interpolación general. La desventaja radica en el submuestreo de la información entre las curvas de nivel, en especial en las áreas de relieve bajo.

Al comienzo del proceso de interpolación, De topo a ráster emplea información inherente a las curvas de nivel para construir un modelo de drenaje generalizado inicial. Esto se realiza mediante la identificación de puntos de curvatura local máxima en cada curva de nivel. Una red de arroyos y crestas curvilineales que interseca esos puntos se deriva utilizando la cuadrícula de elevación inicial (Hutchinson, 1988). Las ubicaciones de esas líneas se actualizan iterativamente a medida que las elevaciones de DEM se actualizan iterativamente. Esta información se utiliza para garantizar propiedades hidrogeomorfológicas adecuadas del DEM de salida y también se pueden utilizar para verificar la precisión del DEM de salida.

Los puntos de datos de curvas de nivel también se utilizan en la interpolación de valores de elevación en cada celda. Todos los datos de curvas de nivel se leen y generalizan. Un máximo de 100 puntos de datos se leen desde las curvas de nivel dentro de cada celda, con el valor de la elevación promedio que se utiliza como el único punto de datos de elevación para cada celda que interseca los datos de línea de curvas de nivel. En cada resolución DEM, se utiliza solo un punto crítico para cada celda. Por esta razón, es redundante tener una densidad de curvas de nivel con varias curvas de nivel que cruzan las celdas de salida.

Una vez que se determina la morfología general de la superficie, los datos de curvas de nivel también se utilizan en la interpolación de valores de elevación en cada celda.

Cuando estos datos se usan para interpolar información de elevación, se leen y generalizan todos los datos de curvas de nivel. Se lee un máximo de 50 puntos de datos de estas curvas de nivel dentro de cada celda. Después de la resolución final, se utiliza solo un punto crítico para cada celda. Por esta razón, es redundante tener una densidad de curvas de nivel con varias curvas de nivel que cruzan las celdas de salida.

Uso de datos de lago

Los polígonos de lago en versiones anteriores de Topo a ráster eran simples máscaras que establecían la elevación de cada superficie de lago a la elevación mínima de todos los valores del DEM vecinos del lago. Se actualizó el algoritmo de los límites del lago para habilitar la determinación automática de las alturas del lago que son totalmente compatibles con líneas de arroyos de conexión y los valores de elevación vecinos.

El método revisado de los límites del lago también trata cada límite del lago como una curva de nivel con elevación desconocida y calcula iterativamente la elevación de esta curva de nivel desde los valores de celda en el límite del lago. Al mismo tiempo, la elevación de cada límite del lago se establece para ajustarse a las elevaciones de cualquier lago de aguas arriba o aguas abajo. La elevación de cada límite del lago también se establece para que sea coherente con los valores de DEM vecinos. Los valores de celda inmediatamente fuera del lago se establecen para estar sobre la elevación del límite del lago y los valores de celda en el interior del lago se establecen para estar por debajo de la elevación del límite del lago.

Los límites del lago pueden incluir islas dentro de lagos y lagos dentro de islas. Todos los valores de DEM que se encuentran dentro de los lagos, según lo determinan los polígonos del límite del lago, se establecen a la altura estimada del DEM en el límite del lago.

Uso de datos de acantilados

Las líneas de acantilado permiten un corte completo continuo entre los valores de celda vecinos en cada lado de las líneas de datos de acantilados, según están codificados en el ráster de salida. Las líneas de acantilados se deben proporcionar como líneas dirigidas, con el lado bajo de cada línea de acantilado a la izquierda y el lado alto de la línea de acantilado a la derecha.. Esto permite la eliminación de puntos de datos de elevación que se encuentran en el lado equivocado de los acantilados, según están codificados en el ráster, así como una mejor ubicación de los acantilados en relación con los arroyos.

También se ha descubierto que los pequeños cambios en la posición de los arroyos y los acantilados a medida que se incorporan en el ráster pueden dar lugar a falsas interacciones entre estos datos. Por lo tanto, se ha desarrollado un método automatizado para realizar pequeños ajustes en la ubicación de las líneas de arroyos y acantilados en el ráster para minimizar estas falsas interacciones.

Uso de datos de líneas de costa

Las celdas en el DEM de salida final que se encuentran fuera de los polígonos especificados por esta clase de entidad poligonal se establece en un valor especial determinado internamente y que es menor que el límite de altura mínima especificado por el usuario. El resultado de esto es que un polígono costero completo se puede utilizar como entrada y automáticamente se recorta a la extensión del procesamiento.

Interpolación de varias resoluciones

El programa utiliza un método de interpolación de varias resoluciones, que comienza con un ráster grueso y sigue con la resolución más fina especificada por el usuario. En cada resolución, se aplican las condiciones de drenaje, se lleva a cabo la interpolación y se registra la cantidad de sumideros restantes en el archivo de diagnóstico de salida.

Procesar datos de arroyos

La herramienta De topo a ráster requiere que los datos de la red de arroyos tengan todos los arcos apuntando hacia abajo y que no haya polígonos (lagos) en la red.

Los datos de arroyos deben estar compuestos por arcos individuales en un patrón dendrítico, con los bancos de arroyos paralelos, polígonos de lago, etc. que se limpien a través de la edición interactiva. Al editar polígonos de lago fuera de la red, se debe ubicar un arco individual desde el principio al final del área contenida. El arco debe seguir la ruta de un lecho fluvial histórico si se conoce o existe uno. Si se conoce la elevación del lago, el polígono del lago y su elevación se pueden utilizar como entrada de CURVA DE NIVEL.

Para visualizar la dirección de las secciones de línea, cambie la Simbología a la opción Flecha en el extremo. Con esta opción se dibujarán las secciones de línea con un símbolo de flecha que muestra las direcciones de las líneas.

Crear rásteres adyacentes y realizar mosaicos con ellos

A veces es necesario crear DEM a partir de teselas adyacentes de datos de entrada. Generalmente, esto ocurre cuando las entidades de entrada se derivan de una serie de hojas de mapa o cuando, debido a limitaciones de memoria, los datos de entrada deben procesarse en varias partes.

El proceso de interpolación utiliza datos de entrada de las áreas circundantes para definir la morfología y el drenaje de la superficie, e interpolar los valores de salida. Sin embargo, los valores de celda de los bordes de cualquier DEM de salida no son tan fiables como los del área central ya que están interpolados con la mitad de la información.

Para realizar las predicciones más precisas en los bordes del área de interés, la extensión de los datasets de entrada debe ser mayor que el área de interés. El parámetro Margen en celdas proporciona un método para acortar los bordes de los DEM de salida según una distancia especificada por el usuario. Los bordes de las áreas superpuestas deben tener al menos 20 celdas de ancho.

Debe haber un poco de superposición de los datos de entrada en las áreas adyacentes cuando se combinan varios DEM de salida en un ráster único. Sin esta superposición, es posible que los bordes de los DEM fusionados no sean suaves. Las extensiones de los datasets de entrada de cada una de las interpolaciones deben tener un área aún más grande que si solo se fuera a hacer una interpolación para una interpolación simple, de forma de asegurar que los bordes se puedan prever de manera tan precisa como sea posible.

Después de crear los DEM, se pueden combinar mejor utilizando la herramienta de geoprocesamiento Mosaico con las opciones Combinar o Valor medio. Esta función proporciona opciones para administrar áreas superpuestas a fin de suavizar la transición entre los datasets.

Evaluar la salida

Cada superficie creada se debe evaluar para asegurar que los datos y parámetros suministrados al programa resulten en una representación realista de la superficie. Hay varias formas de evaluar la calidad de una superficie de salida, según el tipo de entrada disponible para crear la superficie.

La evaluación más común es crear curvas de nivel a partir de la superficie nueva con la herramienta Curva de nivel y compararlas con los datos de curvas de nivel de entrada. Es mejor crear estas curvas de nivel nuevas a la mitad del intervalo de la curva de nivel original para examinar los resultados entre éstas. Dibujar las curvas de nivel originales y las nuevas una sobre la otra puede ayudar a identificar los errores de interpolación.

Otro método de comparación visual es comparar la cubierta opcional de drenaje de salida con arroyos y crestas conocidos. La clase de entidad de drenaje contiene los arroyos y las crestas que se generaron mediante el programa durante el proceso de aplicación de drenaje. Estos arroyos y crestas deben coincidir con los arroyos y las crestas conocidos del área. Si se utilizó una clase de entidad de arroyo como entrada, los arroyos de salida deben superponerse casi perfectamente con los arroyos de entrada, aunque pueden ser levemente más generalizados.

Un método común para evaluar la calidad de una superficie generada es retener un porcentaje de los datos de entrada del proceso de interpolación. Después de generar la superficie, la altura de estos puntos conocidos se puede sustraer de la superficie generada para examinar con qué precisión la superficie nueva representa la superficie real. Estas diferencias se pueden usar para calcular una medida de error para la superficie, como el error cuadrático medio (RMS).

De topo a ráster tiene un conjunto completo de procedimientos para evaluar la calidad del DEM ajustado, optimizar la resolución del DEM y detectar errores en los datos de entrada.

Influencia de las curvas de nivel

Existe una influencia menor en el algoritmo de interpolación que hace que las curvas de nivel de entrada tengan un efecto más fuerte en la superficie de salida en la curva de nivel. Esta influencia puede resultar en un aplanamiento leve de la superficie de salida cuando cruza la curva de nivel. Esto puede generar resultados confusos al calcular la curvatura del perfil de la superficie de salida pero no es perceptible de otro modo.

Causas posibles de problemas con De topo a ráster

Si encuentra algún problema al ejecutar De topo a ráster, revise los siguientes puntos para obtener explicaciones y soluciones para los problemas más comunes.

Referencias

ANU Fenner School of Environment and Society and Geoscience Australia, 2008. GEODATA 9 Second DEM and D8 Digital Elevation Model and Flow Direction Grid, User Guide. Geoscience Australia, 43 pp. Consulte: http://www.ga.gov.au/image_cache/GA11644.pdf.

Goodchild, M. F. y D. M. Mark. 1987. The fractal nature of geographic phenomena. Annals of Association of American Geographers. 77 (2): 265–278.

Hutchinson, M. F. 1988. Calculation of hydrologically sound digital elevation models. Documento presentado en el Tercer simposio internacional sobre manejo de datos espaciales en Sydney, Australia.

Hutchinson, M. F. 1989. A new procedure for gridding elevation and stream line data with automatic removal of spurious pits. Journal of Hydrology, 106: 211–232.

Hutchinson, M. F. y T. I. Dowling. 1991. A continental hydrological assessment of a new grid-based digital elevation model of Australia. Hydrological Processes 5: 45–58.

Hutchinson, M. F. 1993. Development of a continent-wide DEM with applications to terrain and climate analysis. En Environmental Modeling with GIS, ed. M. F. Goodchild et al., 392–399. Nueva York: Oxford University Press.

Hutchinson, M. F. 1996. A locally adaptive approach to the interpolation of digital elevation models. En Procedimientos, Tercer conferencia/taller internacional sobre la integración del modelado ambiental con SIG. Santa Barbara, CA: National Center for Geographic Information and Analysis. Consulte: http://www.ncgia.ucsb.edu/conf/SANTA_FE_CD-ROM/sf_papers/hutchinson_michael_dem/local.html.

Hutchinson, M.F. 2000. 2000. Optimising the degree of data smoothing for locally adaptive finite element bivariate smoothing splines. ANZIAM Journal 42(E): C774–C796.

Hutchinson, M.F. and Gallant, J.C. 2000. Digital elevation models and representation of terrain shape. In: J.P. Wilson and J.C. Gallant (eds) Terrain Analysis. Wiley, New York, pp. 29–50.

Hutchinson, M.F. 2008. Adding the Z-dimension. In: J.P. Wilson and A.S. Fotheringham (eds), Handbook of Geographic Information Science, Blackwell, pp 144–168.

Hutchinson, M.F., Stein, J.A., Stein, J.L. and Xu, T. 2009. Locally adaptive gridding of noisy high resolution topographic data. In Anderssen, R.S., R.D. Braddock and L.T.H. Newham (eds) 18th World IMACS Congress. Modelling and Simulation Society of Australia and New Zealand and International Association for Mathematics and Computers in Simulation, July 2009, pp. 2493–2499. Consulte: http://www.mssanz.org.au/modsim09/F13/hutchinson.pdf.

Hutchinson, M.F., Xu, T. and Stein, J.A. 2011. Recent Progress in the ANUDEM Elevation Gridding Procedure. En: Geomorphometry 2011, edited by T. Hengel, I.S. Evans, J.P. Wilson and M. Gould, pp. 19–22. Redlands, California, USA. Consulte: http://geomorphometry.org/HutchinsonXu2011.

Wahba, G. 1990. Spline models for Observational data. Documento presentado en la Serie de conferencias regionales de CBMS-NSF sobre matemáticas aplicadas. Philadelphia: Soc. Ind. Appl. Maths.

Temas relacionados

5/10/2014