Calculating NDVI from MODIS rasters using GDAL and Python