La calculadora de rasterización de QGIS, ¿se ejecuta como un proceso por lotes?

Question

Tengo una serie de archivos ráster (150 más o menos) que necesito multiplicar por un ráster de referencia. Desgraciadamente, la serie de rásters tiene un tamaño de celda diferente (30m×30m, es decir, 0,00027777778 grados × 0,000277778 grados) que el archivo de línea de base (,1 grado × 0,1 grado).

Hice una prueba utilizando la calculadora de rasterización de qgis (que está en la ruta de menú: raster->calculadora de rasterización...) y funcionó bien, por lo que puedo decir. Es decir, QGIS fue capaz de reconciliar los diferentes tamaños de celdas y completar la acción como se esperaba. Sin embargo, cuando fui a realizar el proceso por lotes, descubrí que no es posible - al menos no desde la GUI.

Cuando probé los algoritmos SAGA, GRASS y GDAL_CALC parece que los diferentes tamaños de las celdas están presentando problemas. Como nota, el archivo de la línea base tiene actualmente una extensión mundial, mientras que los otros archivos son subconjuntos más pequeños.

¿Existe una forma de llamar a la calculadora de rasterización de QGIS como proceso por lotes o desde un lenguaje de scripting?

Ejemplo gdal_calculate entrada:

python gdal_calculate.py --outfile= C:/Users/asd/Documents/MangWork/TestImage/gdalcalcout.tif \
                         --calc "(han*adj)" \
                         --han=C:/Users/asd/Documents/MangWork/TestImage/cl_MEX_30N_110W.tif \
                         --adj=C:/Users/asd/Documents/MangWork/TestImage/Correction_Raster_Comp.tif \
                         --cellsize=MINOF

Answer 1

¿Por qué no cambias el tamaño de las celdas del archivo base? Si las resoluciones son las mismas, podrás utilizar correctamente GRASS (y si es la misma región) o GDAL_CALC. De hecho, es lo que probablemente hace QGIS antes de ejecutar la multiplicación.

De lo contrario, deberías considerar el uso de la consola de Python con la clase QqgRasterCalculator: http://qgis.org/api/classQgsRasterCalculator.html .

Answer 2

Este es también mi problema al restar dos imágenes rasterizadas con diferentes dimensiones. También he intentado gdal_calc He probado la calculadora de rasterización de QGIS, pero como has dicho no es posible con el procesamiento por lotes. También he buscado en todas partes y he preguntado en todas partes, pero sólo dicen que restar imágenes con dos dimensiones diferentes no es posible. No hasta que encontré gdal_calculate.py . Pude restar dos imágenes rasterizadas con diferentes dimensiones usando esta utilidad gdal. Y uno de sus parámetros opcionales es el --cellsize que creo que es lo que necesitas?