if (FALSE) {
### Set-up parameters of interest
coll <- "MOD11A1.061"
bands <- c("LST_Day_1km", "LST_Night_1km")
time_range <- as.Date(c("2017-01-01", "2017-01-30"))
roi <- sf::st_as_sf(
data.frame(
id = "roi_test",
geom = "POLYGON ((-5.82 9.54, -5.42 9.55, -5.41 8.84, -5.81 8.84, -5.82 9.54))"
),
wkt = "geom", crs = 4326
)
### Download and import the data
modis_ts <- mf_modisfast(
collection = coll,
variables = bands,
roi = roi,
time_range = time_range,
earthdata_username = "earthdata_un",
earthdata_password = "earthdata_pw"
)
### Plot the data
terra::plot(modis_ts)
}
Run the code above in your browser using DataLab