consulta la trama celular valores por punto

Question

Tengo un binario de trama (0s y 1s) que necesito para consulta, dado un conjunto de coordenadas de un punto. Esto es bastante fácil de hacer con gdal (gdallocationinfo), pero donde es complicado es que también tengo que devolver la suma de los 8 células barrio. Hay una manera fácil de hacer esto sin el procesamiento de todo el mapa?

En la HIERBA de la única manera que se me ocurre para hacer esto es utilizar r.mapcalc de devolución de una suma de la celda que me interesa y, a continuación, cambie el valor de todas las otras células a 0. Esto me permite calcular una suma de agregados para el mapa con r.sum. Pero esto requiere una gran cantidad de procesamiento y si necesito hacer esto para cientos de puntos, la tarea se vuelve enorme.

Ideas?

Answer 1

Aquí es código python gdal para obtener matrices de NxN numpy de las coordenadas del punto:

Answer 2

En GRASS GIS, uso r.neighbors que calcula nuevos valores en función de las categoría los valores asignados a las células alrededor de él y almacena nuevos valores de las celdas en una capa de mapa raster de salida. Cuenta con un operador de suma (suma de focal). A continuación, consulta el mapa resultante.

Answer 3

Si usted está familiarizado con R, se usa:

column <- extract(RASTER, Samplingsites.spdf, buffer=...)
SamplingSites.df <- cbind(SamplingSites.df, column)

Donde la TRAMA es su binario de trama, y el SpatialPointsDataFrame son el conjunto de coordenadas de un punto.

Usted puede hacer esto desde dentro de hierba con spgrass6 para acelerar el proceso.

Si usted desea permanecer en la hierba quiero suponer que el uso de un flujo de trabajo con:

1. r.los vecinos con su conjunto de puntos(se calcula su suma con un movimiento de la ventana, tamaño:3x3 celdas ... así que su punto y 8 puntos adicionales)
2. v. qué.rast para hacer una consulta espacial