En esta práctica se explica como obtener datos satelitales y almacenarlos en el directorio de trabajo. Estos datos serán utilizados posteriormente en las próximas prácticas. Aunque las imagenes pueden obtenerse por descarga directa en este caso utilizaremos código para almanecar las imágenes.
A continuación el código de descarga:
dir.create('data', showWarnings = FALSE)
if (!file.exists('data/rs/samples.rds')) {
download.file('https://biogeo.ucdavis.edu/data/rspatial/rsdata.zip', dest = 'data/rsdata.zip')
unzip('data/rsdata.zip', exdir='data')
}Es importante verificar y recordar que se creó una carpeta llamada data y otra carpeta dentro de esta llamada rs donde se almacenaron los datos. Estas carpetas se ubicaron por defecto en el directorio de trabajo.
Aunque no está incluido en el taller podemos verificar que existen archivos en estas carpetas y cuales son de la siguiente forma, utilizanto la funcion *list.files:
list.files(path = "data/rs/")## [1] "centralvalley-2001LE7.tif" "centralvalley-2011LT5.tif"
## [3] "cropped-landsat.tif" "LC08_044034_20170614_B1.tif"
## [5] "LC08_044034_20170614_B10.tif" "LC08_044034_20170614_B11.tif"
## [7] "LC08_044034_20170614_B2.tif" "LC08_044034_20170614_B3.tif"
## [9] "LC08_044034_20170614_B4.tif" "LC08_044034_20170614_B5.tif"
## [11] "LC08_044034_20170614_B6.tif" "LC08_044034_20170614_B7.tif"
## [13] "LC08_044034_20170614_B8.tif" "LC08_044034_20170614_B9.tif"
## [15] "nlcd-L1.tif" "nlcd-L2.tif"
## [17] "samples.rds"De esta forma verificamos que se descargaron correctamente 17 archivos, 11 de ellos corresponden a las bandas de Landsat 8 y los restantes a imagenes “.tiff”.