Skip to contents

Methods to create a SpatRasterDataset. This is an object to hold "sub-datasets", each a SpatRaster that in most cases will have multiple layers.

See describe for getting information about the sub-datasets present in a file.

Usage

# S4 method for missing
sds(x) 

# S4 method for character
sds(x, ids=0, opts=NULL, raw=FALSE)

# S4 method for SpatRaster
sds(x, ...) 

# S4 method for list
sds(x) 

# S4 method for array
sds(x, crs="", extent=NULL)

Arguments

x

character (filename), or SpatRaster, or list of SpatRasters, or missing. If multiple filenames are provided, it is attempted to make SpatRasters from these, and combine them into a SpatRasterDataset

ids

optional. vector of integer subdataset ids. Ignored if the first value is not a positive integer

opts

character. GDAL dataset open options

raw

logical. If TRUE, scale and offset values are ignored

crs

character. Description of the Coordinate Reference System (map projection) in PROJ.4, WKT or authority:code notation. If this argument is missing, and the x coordinates are within -360 .. 360 and the y coordinates are within -90 .. 90, longitude/latitude is assigned

extent

SpatExtent

...

additional SpatRaster objects

Value

SpatRasterDataset

See also

Examples

s <- rast(system.file("ex/logo.tif", package="terra"))   
x <- sds(s, s/2)
names(x) <- c("first", "second")
x
#> class       : SpatRasterDataset 
#> subdatasets : 2 
#> dimensions  : 77, 101 (nrow, ncol)
#> nlyr        : 3, 3 
#> resolution  : 1, 1  (x, y)
#> extent      : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
#> coord. ref. : Cartesian (Meter) 
#> source(s)   : logo.tif, memory 
#> names       : first, second 
length(x)
#> [1] 2

# extract the second SpatRaster
x[2]
#> class       : SpatRaster 
#> dimensions  : 77, 101, 3  (nrow, ncol, nlyr)
#> resolution  : 1, 1  (x, y)
#> extent      : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
#> coord. ref. : Cartesian (Meter) 
#> source(s)   : memory
#> varname     : logo 
#> names       :   red, green,  blue 
#> min values  :   0.0,   0.0,   0.0 
#> max values  : 127.5, 127.5, 127.5 

a <- array(1:9, c(3,3,3,3))
sds(a)
#> class       : SpatRasterDataset 
#> subdatasets : 3 
#> dimensions  : 3, 3 (nrow, ncol)
#> nlyr        : 3, 3, 3 
#> resolution  : 1, 1  (x, y)
#> extent      : 0, 3, 0, 3  (xmin, xmax, ymin, ymax)
#> coord. ref. :  
#> source(s)   : memory