Hey Milton,

yes there is a way. You could read you raster as numpy array and conduct any kind of analyis on it. Of course basic scriting is necessary. Here you'll find a nice compendium on gdal/ogr snippets that you could reuse for your purpose: http://pcjericks.github.io/py-gdalogr-cookbook/. In particular I like also this presentation that helped me some years ago to start with the gdal and python: http://www.gis.usu.edu/~chrisg/python/2008/os5_slides.pdf.

Best regards,

Jewgenij

Thank you Jewgenij

I have worked on Osgearth for developing GIS appilcation. Osgearth uses gdal libraries. But using Osgearth one can do lot many things.

http://osgearth.readthedocs.org/en/latest/user/earthfiles.html

A simple Python raster calculator is actually shipped with GDAL - gdal_calc.py

http://www.gdal.org/gdal_calc.html

Usage is pretty straightforward, you can declare input images to the letters A-Z, declare an output file as well as calculation. Your Example would look something like this:

gdal_calc.py -A imageA.tif -B imageB.tif --outfile=ABmultiplied.tif --calc="A*B"

In the current version of GDAL you can even employ numpy array functions on your input files (e.g. logical_and, logical_or, etc.)

If you find the GDAL handling from Python a bit cumbersome might I suggest taking a look at the rasterio module developed by the Mapbox team: https://github.com/mapbox/rasterio

in a Unix machine (Linux or OSX) you can do  it in the terminal very easily:

gdal_calc.py -A file1.tif --B file2.tif --outfile=fileresult.tif --calc="A*B"