@jul está en lo cierto al decir que "necesita calcular una superficie de tendencia y luego restarla de su DEM inicial para obtener la 'desredagada'", pero parece que en este caso se necesitan procedimientos más simples para "preservar las caídas". Si la "superficie de tendencia" sigue muy de cerca el DEM original, entonces los residuos no retendrán las características locales de la superficie. Por lo tanto, entre las técnicas que se deben evitar están todas las locales (splines, filtros y, especialmente, kriging) y las que se prefieren son globales.
Un enfoque simple, robusto y directo es ajustar un avión al DEM en las cercanías del río. Esto no requiere ninguna tecnología sofisticada o cálculos pesados, porque (de acuerdo con Euclides) un plano está determinado por tres puntos (no colineales) en el espacio. En consecuencia, seleccione un punto (x1 ', y1', z1 ') = (coordenadas, elevación) en la cabecera del río, otro punto (x2', y2 ', z2') en el extremo aguas abajo y un tercer punto en (x0 ', y0', z0 ') de su elección lejos del segmento de línea que conecta los dos primeros puntos. (Estas coordenadas se indican con números primos porque pronto los cambiaremos). Este último punto no¡tiene que corresponder a un punto en o incluso cerca de la superficie del suelo! De hecho, una buena opción inicial podría ser establecer su elevación al promedio de las elevaciones aguas arriba y aguas abajo, z0 '= (z1' + z2 ') / 2.
Los cálculos se facilitan adoptando el punto (x0 ', y0', z0 ') como el origen de un sistema de coordenadas local. En estas coordenadas, los otros dos puntos están en
(x1,y1,z1) = (x1'-x0',y1'-y0',z1'-z0')
(x2,y2,z2) = (x2'-x0',y2'-y0',z2'-z0').
Cualquier ubicación arbitraria, en (x ', y') en el sistema de coordenadas original, tiene coordenadas (x, y) = (x'-x0 ', y'-y0') en este nuevo sistema. Debido a que cualquier plano que pase por el origen (0,0,0) debe tener una ecuación de la forma z = a * x + b * y, esto reduce el problema a lo siguiente:
Encuentre una ecuación en la forma z = a * x + b * y para el plano que pasa por los puntos (0,0,0), (x1, y1, z1) y (x2, y2, z2).
La solución única es calcular
u = z1 y2 - z2 y1
v = x1 z2 - x2 z1
w = x1 y2 - x2 y1
en que terminos
a = u/w, b = v/w.
Habiendo encontrado estos dos números a y b, y recordando las dos coordenadas originales x0 'e y0', un cálculo ráster de la forma
[DEM] - a * ([X'] - x0') - b * ([Y'] - y0')
elimina la "inclinación" del DEM . En esta expresión, [X '] se refiere a la cuadrícula de coordenadas x en las coordenadas originales e [Y'] se refiere a la cuadrícula de coordenadas y original. Se garantiza que el DEM resultante tendrá la misma elevación (es decir, z0 ') en cada uno de los tres puntos que eligió originalmente; ¡lo que hace en otros lugares depende del DEM mismo!
(Espero que algunos lectores aprecien cómo este enfoque evita todas las referencias a la trigonometría o la maquinaria de mínimos cuadrados. :-)
He encontrado que la tendencia con el análisis de tendencias no es satisfactoria para los ríos. Puede funcionar en regiones como el medio oeste, pero en el oeste una sola superficie lineal o polinómica nunca parece aproximarse muy bien al gradiente del río debido a la compleja topografía. Lo que he usado en su lugar es un promedio ponderado donde comparo la elevación de un píxel de tierras altas con un promedio ponderado del río usando una densidad de kernel (ver artículo en ArcUser de invierno de 2010). Esto produce un resultado más "localizado". Todavía es posible obtener algunos valores negativos, pero he descubierto que esos valores parecen estar confinados al canal del río y son más comunes en áreas de gradiente empinado. Más recientemente, he experimentado con enfoques de distancia de ruta para mapear "altura de inundación", pero esto no ganó '
fuente
Necesita calcular una superficie de tendencia, y luego restarla de su DEM inicial para obtener la "destendada".
Para calcular una superficie de tendencia, se puede utilizar cualquier función de simplificación de superficie (filtrado, suavizado de spline, kriging, etc.). Esta elección depende de cómo desee que se vea esta superficie de tendencia.
Aquí hay un ejemplo sobre cómo usar kriging en arcGIS 10.0 para eso.
fuente