RGISTools (version 0.9.7)

lsMosaic: Mosaic a set of Landsat-7 or Landsat-8 images

Description

lsMosaic merges the Landsat-7 or Landsat-8 imagery that covers a region of interest on the same dates.

Usage

lsMosaic(src, AppRoot, extent = NULL, out.name = "outfile",
  verbose = FALSE, gutils = FALSE, overwrite = FALSE, ...)

Arguments

src

the path to the folder with the Landsat-7 or Landsa-8 images in GTiff format.

AppRoot

the directory to save the mosaicked images.

extent

an axtent, Raster*, or Spatial* object representing the region of interest.

out.name

the name of the folder that stores the outputs. By default, <U+201C>outfile<U+201D> is assigned.

verbose

logical argument. If TRUE, the function prints the running steps and warnings.

gutils

logical argument. If TRUE, the function uses `GDAL' utilities for mosaicking.

overwrite

logical argument. If TRUE, overwrites the existing images with the same name.

...

arguments for nested functions:

• pathrow a list of vectors with the path and row numbers of the tiles concerning the region of interest.

• bFilter a vector with the bands to be mosaicked. If not supplied, all bands are mosaicked.

• dayFilter a vector with the capturing dates being considered for mosaicking. If not supplied, all dates are mosaicked.

Details

The function mosaics the imagery in the src folder. The folder can hold GTiff images from several tiles, dates and bands. When only a subset of bands or dates has to be mosaicked, the band names or dates should be provided through the argument bFilter or dateFilter. Band names are defined by the letter <U+201C>b<U+201D> and the two-digit band number (e.g., <U+2018>b01<U+2019>). The dates must be provided as a Date class object. Once mosaicked, the images can be cropped to fit the extent (optional). The extent can be defined in any coordinate reference system, since lsMosaic automatically reprojects the extent to match the projection of the image. The outputs are placed in the AppRoot directory, under the folder named as out.name. If no name is provided, the folder is named <U+201C>outfile<U+201D>. To use `gutils = TRUE', a proper installation of `GDAL' and the `gdalUtils' library is required. This method is faster than native `R' functions.

Examples

# NOT RUN {
# load a spatial polygon object of Navarre
data(ex.navarre)
# main output directory
src <- paste0(tempdir(),"/Path_for_downloading_folder")
print(src)
# download Landsat-8 images
lsDownload(satellite = "ls8",
           username = "username",
           password = "password",
           startDate = as.Date("01-01-2018", "%d-%m-%Y"),
           endDate = as.Date("20-01-2018", "%d-%m-%Y"),
           extent = ex.navarre,
           untar = TRUE,
           AppRoot = src)
# folder with the Landsat-8 untared images
ls8.src <- file.path(src, "Landsat8")
tif.src <- file.path(ls8.src,"untar")
# mosaic the Landsat-8 images
lsMosaic(src = tif.src,
         AppRoot = ls8.src,
         out.name = "Navarre")

lsMosaic(src = tif.src,
         AppRoot = ls8.src,
         out.name = "Navarre",
         extent = ex.navarre,
         gutils = TRUE, # using gdalUtils
         overwrite = TRUE) # overwrite
# }