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 ##GDAL commands## splice out a single band and save as independent geotiff: gdal_translate -of GTiff -b 2 morflex-freeflight-geo.tif morflex-freeflight-green.tif