Good example in GDAL/Python: http://www.gis.usu.edu/~chrisg/python/2009/lectures/ospy_slides5.pdf Script for GDAL: http://www.jeremymsmith.us/davidson/NDVI.py Remember, NDVI is: Infrared - Visible / Infrared + Visible