Remuestrear binario de trama para dar *proporción* dentro de la nueva ventana de la celda

Question

Me gustaría ir de un binario trama de bosque/no-bosque a 30m de resolución de una trama en 240m de la resolución, con el valor de cada celda, siendo la proporción de células boscosas - es decir, los valores de los datos para la 240m células que van desde 0 (donde todas las células de los 30m de la trama fueron cero/no-bosque) a través de 0.5 (donde la mitad de la 30m células fueron bosque, la mitad eran no-bosque) a 1 (donde todas las células de los 30m de la trama eran boscosas).

Aunque se ha sugerido a no utilizar la interpolación bilineal sobre datos discretos, como lo que puedo decir que el resultado binario con 0/1 datos sería dar un valor promedio (es decir, una proporción entre el 1 y el 0). Es esta una forma sensata de hacerlo, o es que hay un método mejor?

Puedo usar el Arco, QGIS y Idrisi. Gracias por tu ayuda.

Answer 1

El uso de Bloques de estadísticas.

Esto funciona como Focal de estadísticas mediante el cálculo de un estadístico resumen (como el que significa que usted desee) dentro de un determinado barrio de las células (como un 8 por 8 metros, donde 8 = 240 m / 30 m), excepto que realiza esta solo para regular la subdivisión de la red, en lugar de con un conjunto de la superposición de los barrios, uno en cada celda.

Usted también puede hacer uso de Focales estadísticas , si usted realmente quiere: después de la computación focal media de más de 8 por 8 plazas, volver a muestrear a una de 240 m de red con vecino más cercano de remuestreo. Cuando las rejillas son registrados a la una de la otra (es decir, tienen el mismo origen) esto debería dar el mismo resultado como block statistics. (No voy a garantizar que: algunas de las arbitrarias decisiones tienen que ser tomadas cuando la nueva célula de centros de coincidir con el viejo de las esquinas de las celdas, a medida que va aquí--y, si los diferentes comités codificados los dos procedimientos, que pueden haber tomado decisiones diferentes: poco en ArcGIS es verdaderamente consistente, me temo.)

Otro método es crear una red de zonas, una zona por la plaza donde un promedio deseado, y realizar un zonal resumen como una cuadrícula. Las zonas pueden ser calculada matemáticamente de las cuadrículas de fila y columna de las coordenadas (por medio de la floor o int función, por reclasificación, o uniéndose a una mesa apropiada a la tabla de atributos).

Voy a cerrar diciendo que la interpolación bilineal, aunque, de hecho, el rendimiento de los valores en el rango 0..1, no es lo que usted desea: trabaja por encontrar en la mayoría de los cuatro originales (30 m) de la cuadrícula de celdas que rodean el centro de un nuevo (240 m) de la célula y la interpolación sólo sus valores. Como tal, va a pasar por alto los otros 8*8 - 4 = 60 células originales que caen dentro de cada bloque nuevo. Me ilustrar la interpolación bilineal en http://www.quantdec.com/SYSEN597/GTKAV/section9/map_algebra.htm: la discusión comienza cerca de la mitad de la página.

Answer 2

En ArcGIS whne remuestrea de datos utilizando bilineal de remuestreo que solo se fija en los valores del centro de cuatro celdas (remuestrear la documentación). Como tal, el uso de este método que voy a perder los datos si no compensar la pérdida de datos.

Dado que sabemos que para el cambiadas las células usted está después de una proporción de las células que están forestadas, podemos pensar que es la suma de los 30m células dividido por 64 (hay 64 30m células en el 240m bloque).

Esto significa que si podemos crear una nueva trama con el centro de valores como la suma de los valores que rodean a 30m de resolución, la reducción de la resolución mediante el remuestreo con el vecino más cercano o interpolación bilineal nos dará 240 millones de células que son la suma de los 30m células que cubren. Podemos hacer esto con la focal de la suma de la herramienta a través de la 30m de la trama.

Por último una vez que tenemos una 240m focal de la suma de la trama, dividir por 64 para obtener la proporción de respuesta.

En Idrisi no estoy seguro de que la imagen algoritmos de toma de muestras, de la misma forma en QGIS, pero me imagino que hay algo similar. Ciertamente, en QGIS puede procesar la trama en Python usando scipy ndimage o similar.