¿Área de extracción de DEM por encima de donde se encuentra el 60% de la región con ArcGIS Desktop?

9

Tengo un DEM de una región montañosa. Me gustaría extraer las áreas del DEM por encima de las cuales se encuentra el 60% del área total del DEM. Entonces, esencialmente me gustaría extraer las áreas que representan las partes superiores de las montañas.

He oído hablar de esto conocido como la línea hipsométrica (línea H60), pero no he podido localizar ninguna documentación sobre cómo lograrlo.

Mi instinto es calcular de alguna manera las estadísticas del DEM para decirme dónde está el nivel del 60% y luego contornear el DEM en esta elevación. Luego puedo recortar el DEM por este contorno y realizar mis cálculos en los subconjuntos de DEM de salida.

¿Cómo determino este nivel de 60%?

Radar
fuente

Respuestas:

7

Creo que solo está buscando el valor del percentil 60 para su ráster. La forma más fácil de hacer esto probablemente es a través de la funcionalidad del histograma. Haga clic derecho en la capa en la Tabla de contenido, luego vaya a Propiedades de capa> Simbología> Clasificado> Clasificar. Elija Cuantil como Método y use 10 clases. La sexta clase comenzará en el 60% del valor del percentil. Vea la imagen a continuación para ver un ejemplo.

Si desea especificar otros valores de porcentaje, puede presionar el símbolo "%" en el lado derecho de la ventana de Clasificación, al lado de donde dice "Romper valores". ingrese la descripción de la imagen aquí

dmahr
fuente
3
Esto es correcto, excepto que Radar buscará el percentil 40 en lugar del 60. Solo el 40% de la tierra estaría por encima del percentil 60.
nmpeterson el
66
¡Cuidado con los errores de ArcGIS! No he probado Arc 10, pero las versiones anteriores no podían calcular las estadísticas correctamente. ¿Has notado cómo esta captura de pantalla indica que hay "9667 elementos en la clase" entre un "Count:" de 85842? Eso es sustancialmente más del 10% (que es lo que debería ser). No hay nada de malo en usar el histograma para adivinar el percentil 40 (o lo que sea), pero luego debe verificar esta suposición (usando cálculos de Álgebra de mapas) para asegurarse de que sea correcto.
whuber
Oh wow ... gracias por señalar eso Whuber. Nunca supe sobre ese error. Suspiro ... es difícil creer que ArcMap no pueda hacer cosas tan simples de manera confiable.
dmahr
Genial, esto es justo lo que estaba buscando. Me aseguraré de verificar los valores.
Radar
44
@dmahr Si desea que se incluyan 'todas' las muestras, aumente el tamaño de 'Muestreo ...' a cualquier valor superior al 'conteo'; 90,000 deberían funcionar para la información en la captura de pantalla anterior. Por defecto, ArcGIS solo incluirá hasta los primeros 10,000 registros / valores.
SaultDon
3

Un enfoque alternativo para la extracción DEM propuesta, suponiendo que tenga Spatial Analyst, sería:

  1. Use la herramienta 'Slice' con el método "ÁREA IGUAL" y 5 o 10 zonas.
  2. Entonces podrías:

    a) Reclasifique la salida de Slice para crear una máscara para extraer su DEM

    - o -

    b) Úselo con la herramienta 'CON' o la Calculadora ráster para enmascarar su DEM en las 3 o 6 zonas principales (respectivamente, para el número de zonas de Slice definidas)

Además, Slice parece obtener cuantiles algo más precisos que el método de clasificación Cuantil (como lo indicó @whuber el 28/11/2011 20:29)

Joe Champagne
fuente