stack: Create a RasterStack object

In raster: Geographic Data Analysis and Modeling

Description

A RasterStack is a collection of RasterLayer objects with the same spatial extent and resolution. A RasterStack can be created from RasterLayer objects, or from raster files, or both. It can also be created from a SpatialPixelsDataFrame or a SpatialGridDataFrame object.

Usage

## S4 method for signature 'character'
stack(x, ..., bands=NULL, varname="", native=FALSE, RAT=TRUE, quick=FALSE)

## S4 method for signature 'Raster'
stack(x, ..., layers=NULL)

## S4 method for signature 'missing'
stack(x)

## S4 method for signature 'list'
stack(x, bands=NULL, native=FALSE, RAT=TRUE, ...)

Arguments

x	filename (character), Raster* object, missing (to create an empty RasterStack), SpatialGrid*, SpatialPixels*, or list (of filenames and/or Raster* objects). If x is a list, additional arguments ... are ignored
bands	integer. which bands (layers) of the file should be used (default is all layers)
layers	integer (or character with layer names) indicating which layers of a RasterBrick should be used (default is all layers)
native	logical. If TRUE native drivers are used instead of gdal drivers (where available, such as for BIL and Arc-ASCII files)
RAT	logical. If TRUE a raster attribute table is created for files that have one
quick	logical. If TRUE the extent and resolution of the objects are not compared. This speeds up the creation of the RasteStack but should be use with great caution. Only use this option when you are absolutely sure that all the data in all the files are aligned, and you need to create RasterStack for many (>100) files
varname	character. To select the variable of interest in a NetCDF file (see raster)
...	additional filenames or Raster* objects

Value

RasterStack

See Also

addLayer, dropLayer, raster, brick

Examples

# file with one layer
fn <- system.file("external/test.grd", package="raster")
s <- stack(fn, fn)
r <- raster(fn)
s <- stack(r, fn) 
nlayers(s)

# file with three layers
slogo <- stack(system.file("external/rlogo.grd", package="raster")) 
nlayers(slogo)
slogo

Example output

Loading required package: sp
[1] 2
[1] 3
class       : RasterStack 
dimensions  : 77, 101, 7777, 3  (nrow, ncol, ncell, nlayers)
resolution  : 1, 1  (x, y)
extent      : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=merc +datum=WGS84 
names       : red, green, blue 
min values  :   0,     0,    0 
max values  : 255,   255,  255 

raster documentation built on Jan. 5, 2021, 3:01 a.m.