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.

# S4 method for missing
sds(x) 

# S4 method for character
sds(x, ids=0)

# 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 SpatRaster objects, 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

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
#> 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