Estoy tratando de evaluar el cambio de la cubierta terrestre que se ha producido en Asia oriental en los últimos 18 años sólo en los inviernos (octubre-enero) para lo cual quiero utilizar el conjunto de datos MOD13Q1 de modis de 2000-2018. Estoy utilizando la función mt_subset para descargar los datos.
Ya he descargado una muestra de datos de los Países Bajos, pero no consigo entender cómo obtener el conjunto de datos que necesito a través de esta función, que debe contener todos los productos de octubre a enero (2000-2018) y la capa debe ser 250m NDVI. Este sitio web proporciona información sobre la biblioteca para la descarga de productos https://cran.r-project.org/web/packages/MODISTools/vignettes/modistools-vignette.html
Así es como estoy descargando los datos.
VI <- mt_subset(product = "MOD13Q1",
site_id = "nl_gelderland_loobos",
band = "250m_16_days_NDVI",
start = "2000-01-01",
km_lr = 2,
km_ab = 2,
site_name = "testsite",
internal = TRUE,
progress = TRUE)
Lo ideal es que la ubicación sea la que se muestra en la imagen (este de China) y que las fechas de inicio y fin sean sólo para los inviernos, de modo que sólo descargue productos para esa zona y fechas concretas. ¿Cómo debo modificar los argumentos de esta función para cumplir esta tarea?