¿Cálculo del área total que se superpone a un intervalo de color DEM raster específico cuyos límites están delimitados por otro vector?

Question

Sé que esto suena bastante difícil, pero quizá no lo sea para ti o al menos no debería serlo si lo explico correctamente.

(Estoy usando QGIS 2.14) Me gustaría saber para calcular el área total que se superpone a un color DEM raster específico (intervalo de altitud especificado por el usuario) cuyos límites están confinados por otro vector.

Así que básicamente hay un raster DEM como mapa de fondo, he asignado un color específico a un intervalo de altitud. Apliqué un vector que creé para mostrar dónde están los límites de un área determinada. Ahora, dentro de esta área, apliqué otra trama que se solapa parcialmente con el intervalo de altitud que apliqué en el DEM de fondo.

Me gustaría saber cuál es el área, dentro de los límites del vector, que se solapa con mi intervalo de altitud en el mapa DEM.

EDIT: He intentado ejecutar calculadora raster en las intersecciones a continuación, las estadísticas zonales comprobación de "recuento" y "suma", pero me sale una enorme tabla de 8000 + filas y dos columnas, sin embargo, cada célula por lo general tienen un valor alrededor de 0,4 o 1,x, así que no sé cómo interpretar esto o si debo de alguna manera la suma de todas esas células juntos para obtener el número de píxeles que caen en mi área de interés (superposición de rango de altitud, sombreado y dentro de vector bounday).

Answer 1

Para poner las cosas claras, supongo que quieres calcular el área que responde a las siguientes condiciones:

• situado dentro de su intervalo de altitud DEM
• situado dentro de los límites de su vector
• situado donde su segundo raster (llamémosle raster2) tiene valores distintos de "nodata" (esta es la parte poco clara de su mensaje, siéntase libre de corregirme si lo he entendido mal).

Hay varias formas de hacerlo; ésta es la que creo que es la más sencilla (llamemos z1-z2 a tu intervalo de altitud):

1. En primer lugar, sus datos deben tener un proyección equiárea si quieres calcular áreas significativas. Por lo tanto, si aún no es el caso, elija una proyección que se ajuste a la región que está mapeando y, a continuación, vuelva a proyectar sus capas raster y vectoriales en ella (utilizando "Guardar como" para las capas vectoriales y "Raster > Proyecciones > Deformar" para las raster).

2. utilice la calculadora de raster para crear una máscara que represente las partes de raster2 situadas dentro de su intervalo DEM. La fórmula será algo como ("raster2@1" != -9999) * ("dem@1" >= z1) * ("dem@1" <= z2) . Debería obtener un raster con 1s en su intervalo de altitud y 0s o "nodata" en todos los demás lugares. (Nota "raster2@1" != -9999 debería estar bien, pero dependiendo de los valores de raster2 puede que tengas que ajustarlo).

3. instale el Plugin de estadísticas zonales . Este plugin calcula estadísticas sobre una trama, dentro de polígonos. Una vez instalado, se encontrará en el menú "raster". Elija la máscara creada anteriormente como "capa rasterizada", y la capa vectorial que contiene sus límites como "capa de polígonos". Si esta capa vectorial contiene muchos polígonos, puede que prefiera exportar el polígono que representa sus límites como una capa vectorial separada (le llevará menos tiempo). Si todo va bien, obtendrá campos adicionales en su capa vectorial; entre ellos habrá una columna "count" que representa el número de píxeles de máscara dentro de cada polígono de la capa vectorial. Sólo tienes que multiplicar el número de esta columna por el área de un píxel (tamaño de celda) -que depende de la proyección- para obtener el área final.