Query of Amazon AWS OLI-Landsat 8 cloud service
oli.asw(path, row, dates, cloud.cover = 10, processing)
data.frame object with:
entityId - Granule ID
L = Landsat
X = Sensor
SS = Satellite
PPP = WRS path
RRR = WRS row
YYYYMMDD = Acquisition date
yyyymmdd = Processing date
CC = Collection number
TX = Collection category
acquisitionDate - POSIXct YYYY-MM-DD (eg., 2015-01-02)
cloudCover -
processingLevel - USGS processing level
path - Landsat path
row - Landsat row
landsat path
landsat row
dates, single or start-stop range in YYYY-MM-DD format
percent cloud cover
processing level ("L1GT" or "L1T")
Jeffrey S. Evans <jeffrey_evans@tnc.org>
if (FALSE) {
# Query path 126, row 59, 2013-04-15 to 2017-03-09, <20% cloud cover
( p126r59.oli <- oli.asw(path=126, row=59, dates = c("2013-04-15", "2017-03-09"),
cloud.cover = 20) )
# Download images from query
bands <- c("_B1.TIF", "_B2.TIF", "_B3.TIF", "_B4.TIF", "_B5.TIF",
"_B6.TIF","_B7.TIF", "_B8.TIF", "_B9.TIF", "_B10.TIF",
"_B11.TIF", "_BQA.TIF","_MTL.txt")
for(i in 1:length(p126r59.oli$download_url)) {
oli.url <- gsub("/index.html","",p126r59.oli$download_url[i])
all.bands <- paste(oli.url, paste0(unlist(strsplit(oli.url, "/"))[8], bands), sep="/")
for(j in all.bands) {
try(utils::download.file(url=j, destfile=basename(j), mode = "wb"))
}
}
}
Run the code above in your browser using DataLab