Usando Python de código abierto y módulos, ¿cómo puedo realizar la Trama funciones de la Calculadora?

Question

Me gustaría realizar funciones matemáticas similares a las que haría uso de la Calculadora Ráster en ArcMap, sin embargo, me gustaría hacerlo en Python sin utilizar el módulo Arcpy. Tengo el módulo de osgeo, y parece que no debería ser posible, sin embargo, estoy teniendo un tiempo difícil encontrar cualquier útil de la sintaxis. Cualquier ayuda se agradece.

Un ejemplo de la función sería multiplicar cada valor de la celda en un ráster de entrada por 2.

Gracias, -Scott

Answer 1

Sin reinventar la rueda, le sugiero que utilice gdal_calc.py. Ejemplo:

gdal_calc.py -A input.tif --outfile=result.tif --calc="A*2"

Answer 2

Si usted desea hacer esto en python, gdal_calc.py no es de importación amable así que sólo puede ser usado desde la línea de comandos. Una alternativa es el gdal cálculos de paquetes que pueden ser importados y utilizados en secuencias de comandos de python.

from gdal_calculations import Dataset
ds1=Dataset('../testdata/landsat7.tif')
ds2=ds1 * 2
ds2.save('/tmp/test.tif')