# NOT RUN {
# load a spatial polygon object of Navarre
data(ex.navarre)
# set the download folder
s.start <- Sys.time()
wdir <- file.path(tempdir(),"Path_for_downloading_folder")
print(wdir)
# download the images
modDownSearch(product = "MOD09GA",
startDate = as.Date("30-07-2018", "%d-%m-%Y"),
endDate = as.Date("06-08-2018", "%d-%m-%Y"),
username = "username",
password = "password",
AppRoot = wdir,
extract.tif = TRUE,
collection = 6,
extent = ex.navarre)
# set folder path where MOD09GA images will be saved
wdir.mod <- file.path(wdir,"Modis","MOD09GA")
# set the tif folder path
wdir.mod.tif <- file.path(wdir.mod,"tif")
# mosaic and cut navarre region
modMosaic(wdir.mod.tif,
AppRoot = wdir.mod,
out.name = "Navarre",
extent = ex.navarre)
# change src to navarre folder
wdir.mod.navarre <- file.path(wdir.mod,"Navarre")
# calculate NDVI from navarre folder
modFolderToVar(wdir.mod.navarre,
fun = varNDVI,
AppRoot = dirname(wdir.mod.navarre),
overwrite = TRUE)
# change src TS_sample
wdir.mod.ndvi <- file.path(dirname(wdir.mod.navarre),"NDVI")
# create the Rdata
tiles.mod.ndvi<-genSaveTSRData(wdir.mod.ndvi, ts.name = "ModisNDVI")
# remove values out of 0-1 range
tiles.mod.ndvi.lim <- clamp(tiles.mod.ndvi,lower=0,upper=1)
# plot the ndvi images
spplot(tiles.mod.ndvi.lim)
s.end <- Sys.time()
s.end - s.start
# }
Run the code above in your browser using DataLab