[gdal-dev] basic raster math

[Jeshua Lacock]

On Feb 14, 2008, at 2:24 PM, Greg Fiske wrote:

> I'm a new python user, so pardon my simple question.  I'd like to
> do some simple raster math within Python.  For example, I'd like
> to read in 2 Geotiffs and add one to the other -- creating a new
> output Geotiff.  I've experimented with the gdal and Numpy
> libraries, but I haven't made much progress on my own.  Can
> somebody send me some example code?

Hi Greg,

Although somewhat off-topic and not too steal thunder from the  
fabulous GDAL, you might want to have a look at GRASS GIS' mapcalc  
command - it lets you use plain math on any rasters and is quite easy  
(dare I say fun) to use.


Cheers,

Jeshua Lacock, Owner
<http://OpenOSX.com>
phone: 877.240.1364