El cálculo de "pendiente relativa" por Mora-Vahrson método en ArcGIS

Question

Soy un veterano de lurker que finalmente se decidió a registrar en el foro (de ninguna manera impulsado por mi necesidad de tener un cierto una pregunta, por supuesto...).

Fui contactado por un colega pidiendo ayuda calcuating una entrada para un deslizamiento de tierra modelo que está trabajando. Específicamente, el modelo de la siguiente (la de Mora-Vahrson de papel/modelo) pide la pendiente de entrada que no sea en grados o porcentajes, pero "pendiente relativa." Esto se obtiene mediante la fórmula:

Sr= [(ElevMAX-ElevMIN)/(Área)],

donde el el min y max elevaciones son el mínimo y máximo de puntos dentro de un área limitada (más comúnmente 1 km cuadrados); estas búsquedas/cálculos se llevan a cabo a través de la deseada región de estudio, sin que se superpongan.

He intentado varias formas de calcular esta utilizando mis habilidades (limitado; soy un arqueólogo con experiencia en SIG a quien se le pidió ayuda...) y el programa (ArcGIS). Estoy usando un 90m resolución del DEM SRTM de la zona y en la actualidad se establecieron en el uso de la Cuadra herramienta de Estadísticas para calcular el MIN y el MAX) para para una determinada área de 1000 x de 1000m, y, a continuación, entrar en esas resultante rásteres en la Trama de la Calculadora como parte de la Mora-Vahrson ecuación), y dividiendo por 1 (como el área del Bloque de cálculo Stat, expresada en kilómetros cuadrados).

Lo que me gustaría calibre de los más experimentados usuarios de SIG es si o no el Bloque de Estadísticas de proceso que he descrito parece cumplir los requisitos para el modelo de la pendiente relativa de los componentes. Si alguien ha tenido experiencia calcular la pendiente relativa o un proceso similar, yo estaría muy interesado en escuchar acerca de él. Gracias.

Answer 1

El uso de focales estadísticas en lugar de bloque de estadística: cuando se utiliza rectangular barrios este produce los mismos resultados en los centros de las manzanas, pero focal estadísticas son calculadas con el movimiento (superposición) de windows, creando una representación de una superficie de relación de pendientes. Por otra parte, el coordinador de estadísticas puede ser calculada con la más natural de los barrios, tales como círculos, en lugar de la artificial isothetic rectángulos requiere de bloque de estadísticas.

ArcGIS es focal estadísticas comando calcular un intervalo de estadística: que el numerador. Para compensar los efectos de borde, es aconsejable no asumir el movimiento de las ventanas constante de la zona (debido a que sus áreas de disminuir hacia los bordes). En su lugar, crear una unidad de red (usando Not(IsNull(...)) , por ejemplo) y calcular el centro de la suma utilizando el mismo barrio de la especificación. Multiplicar la focal de la suma por el cuadrado de la cellsize (en millones de kilómetros cuadrados, si eso es lo que necesitas) da el "área central" de la cuadrícula. Este cálculo se compensa no sólo para efectos de borde, sino también para cualquier falta (NoData) los valores dentro de la red propia. Este es el denominador: dividir el rango focal por esto para obtener la pendiente relativa.

Si, al final, usted realmente necesita que no se solapan los barrios (aunque no puedo ver el motivo esencial por lo que es necesario), puede submuestra resultante de la cuadrícula a lo largo de cualquier matriz de puntos, como en los centros de los que originalmente la intención de bloques. Quizás sorprendentemente, esta no puede ser menos eficiente que la realización de bloque de estadísticas directamente: el tiempo para llevar a cabo (simple) de la cuadrícula de operaciones tiende a ser más limitado por la sobrecarga de datos de e/S que por los cálculos de sí mismos. Tener una red de resultados puede, al menos, proporcionar una mejor y más informativo mapa del cálculo que el bloque de estadísticas de la cuadrícula.