R scripts to process MODIS data