(1)los puntos para extrar la informacion (en coord. geograficas decimales, en el ejemplo se tomaron 2) :
lat <- c(-0.5,-1.3)
long <- c(0.1,-78)
(2)las fechas de incio y fin (en el ejemplo se puso del 2013 al 2014):
modis.subset$start.date <- rep(2013, nrow(modis.subset))
modis.subset$end.date <- rep(2014, nrow(modis.subset))
(3)el directorio de salida se debe reemplezar al empleado:
SaveDir = "C:/DATA_BASE/MODIS"
(4)la linea define todos los productos del MODIS EVI disponibles
GetDates(Product = "MOD13Q1", Lat = modis.subset$lat[3], Long = modis.subset$long[3])
(5)para revisar los resultados especificar el directorio y el archivo *.asc creado
#############################################################################
#load packages
library (MODISTools)
#centers points
lat <- c(-62.1,-0.5,-1.3)
long <- c(-60.1,0,-78)
#create data frame
modis.subset <- data.frame(lat,long)
#define the years
modis.subset$start.date <- rep(2013, nrow(modis.subset))
modis.subset$end.date <- rep(2014, nrow(modis.subset))
#check products
GetDates(Product = "MOD13Q1", Lat = modis.subset$lat[3], Long = modis.subset$long[3])
#download
MODISSubsets(LoadDat = modis.subset[3,],Products = "MOD13Q1",
Bands = c("250m_16_days_EVI", "250m_16_days_pixel_reliability"), Size = c(1,1),
SaveDir = "C:/DATA_BASE/MODIS",StartDate = TRUE)
#review results
data <- read.csv("C:/DATA_BASE/MODIS/Lat-1.30000Lon-78.00000Start2013-01-01End2014-12-31_MOD13Q1.asc",header = FALSE, as.is = TRUE)