Mosaic raster files using GDAL Utilities

Description

Mosaic raster files using GDAL Utilities

Usage

1
2
mosaic_rasters(gdalfile, dst_dataset, output.vrt = NULL,
  output_Raster = FALSE, trim_margins = NULL, verbose = FALSE, ...)

Arguments

gdalfile

Character. Input files (as a character vector) or a wildcard search term (e.g. "*.tif")

dst_dataset

Character. The destination file name.

output.vrt

Character. Output VRT file. If NULL a temporary .vrt file will be created.

output_Raster

Logical. Return output dst_dataset as a RasterBrick?

trim_margins

Numeric. Pre-crop the input tiles by a fixed number of pixels before mosaicking. Can be a single value or four values representing the left, top, right, and bottom margins, respectively.

verbose

Logical. Enable verbose execution? Default is FALSE.

...

Parameters to pass to gdalbuildvrt or gdal_translate.

Details

This function mosaics a set of input rasters (gdalfile) using parameters found in gdalbuildvrt and subsequently exports the mosaic to an output file (dst_dataset) using parameters found in gdal_translate. The user can choose to preserve the intermediate output.vrt file, but in general this is not needed.

Value

Either a list of NULLs or a list of RasterBricks depending on whether output_Raster is set to TRUE.

Author(s)

Jonathan A. Greenberg (gdalUtils@estarcion.net)

See Also

gdalbuildvrt, gdal_translate

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# We'll pre-check to make sure there is a valid GDAL install
# and that raster and rgdal are also installed.
# Note this isn't strictly neccessary, as executing the function will
# force a search for a valid GDAL install.
gdal_setInstallation()
valid_install <- !is.null(getOption("gdalUtils_gdalPath"))
if(require(raster) && require(rgdal) && valid_install)
{
layer1 <- system.file("external/tahoe_lidar_bareearth.tif", package="gdalUtils")
layer2 <- system.file("external/tahoe_lidar_highesthit.tif", package="gdalUtils")
mosaic_rasters(gdalfile=c(layer1,layer2),dst_dataset="test_mosaic.envi",separate=TRUE,of="ENVI",
		verbose=TRUE)
gdalinfo("test_mosaic.envi")
}