## ----product, message=FALSE---------------------------------------------------
library(terra)
library(luna)

# lists all products that are currently searchable
prod <- getProducts() 
head(prod)

# to find the MODIS products
modis <- getProducts("^MOD|^MYD|^MCD")
head(modis)


## ----prod---------------------------------------------------------------------
product <- "MOD09A1"


## ----dates--------------------------------------------------------------------
start <- "2010-01-01" 
end <- "2010-01-07"


## ----aoi1---------------------------------------------------------------------
ken <- geodata::gadm("Kenya", level=1, path=".") 
ken


## ----aoi2---------------------------------------------------------------------
i <- ken$NAME_1 == "Marsabit"
aoi <- ken[i,]


## ----aoi3---------------------------------------------------------------------
plot(ken, col="light gray")
lines(aoi, col="red", lwd=2)


## ----getfilenames-------------------------------------------------------------
mf <- luna::getModis(product, start, end, aoi=aoi, download = FALSE)
mf


## ----downloaddir--------------------------------------------------------------
datadir <- file.path(dirname(tempdir()), "_modis")
dir.create(datadir, showWarnings=FALSE)


## ----pwd----------------------------------------------------------------------
up <- readRDS("../../../../pwds.rds")
up <- up[up$service == "EOSDIS", ]


## ----download-----------------------------------------------------------------
mf <- luna::getNASA(product, start, end, aoi=aoi, download=TRUE,
                     path=datadir, username=up$user, password=up$pwd)
mf