Skip to contents

The UD represents the minimum area in which an animal has some specified probability of being located (Cumulative Distribution Function).

Usage

# S4 method for class '.UD'
getVolumeUD(x, ...)

Arguments

x

a DBBMM, DBBMMStack, dynBGB, .UD or .UDStack object

...

when x is of class DBBMM or dynBGB, several objects of class DBBMM or dynBGB can be added (see 'Examples')

Note

To obtain this modified UD raster from a DBBMMBurstStack object, transform the object with the UDStack function into a '.UDStack' class, and than use the getVolumeUD function upon the obtained object.

Value

'raster' or 'rasterStack'
If several objects are provided, a list of rasters is returned

Author

Marco Smolla & Anne Scharf

Examples

data(leroydbbmm)
data(dbbmmstack)
data(leroydbgb)
getVolumeUD(leroydbbmm) # for a DBBMM object
#> class      : RasterLayer 
#> dimensions : 42, 45, 1890  (nrow, ncol, ncell)
#> resolution : 759.9677, 759.9677  (x, y)
#> extent     : -17098.3, 17100.24, -15958.9, 15959.75  (xmin, xmax, ymin, ymax)
#> crs        : +proj=aeqd +lat_0=42.73884025 +lon_0=-73.8871629 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source     : memory
#> names      : layer 
#> values     : 1.110223e-16, 1  (min, max)
#> 
getVolumeUD(dbbmmstack) # for a DBBMMStack object 
#> class      : RasterStack 
#> dimensions : 45, 21, 945, 2  (nrow, ncol, ncell, nlayers)
#> resolution : 617.1116, 617.1116  (x, y)
#> extent     : -6478.702, 6480.641, -13884.56, 13885.46  (xmin, xmax, ymin, ymax)
#> crs        : +proj=aeqd +lat_0=42.778423 +lon_0=-73.8871629 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> names      :        Leroy,      Ricky.T 
#> min values : 1.110223e-16, 0.000000e+00 
#> max values :            1,            1 
#> 
getVolumeUD(leroydbgb) # for a dynBGB object 
#> class      : RasterLayer 
#> dimensions : 174, 102, 17748  (nrow, ncol, ncell)
#> resolution : 20, 20  (x, y)
#> extent     : -1022.437, 1017.563, -1740.689, 1739.311  (xmin, xmax, ymin, ymax)
#> crs        : +proj=aeqd +lat_0=42.76087175 +lon_0=-73.89786045 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source     : memory
#> names      : layer 
#> values     : 0, 1  (min, max)
#> 

getVolumeUD(leroydbbmm, leroydbgb) # for several objects 
#> [[1]]
#> class      : RasterLayer 
#> dimensions : 42, 45, 1890  (nrow, ncol, ncell)
#> resolution : 759.9677, 759.9677  (x, y)
#> extent     : -17098.3, 17100.24, -15958.9, 15959.75  (xmin, xmax, ymin, ymax)
#> crs        : +proj=aeqd +lat_0=42.73884025 +lon_0=-73.8871629 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source     : memory
#> names      : layer 
#> values     : 1.110223e-16, 1  (min, max)
#> 
#> 
#> [[2]]
#> class      : RasterLayer 
#> dimensions : 174, 102, 17748  (nrow, ncol, ncell)
#> resolution : 20, 20  (x, y)
#> extent     : -1022.437, 1017.563, -1740.689, 1739.311  (xmin, xmax, ymin, ymax)
#> crs        : +proj=aeqd +lat_0=42.76087175 +lon_0=-73.89786045 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source     : memory
#> names      : layer 
#> values     : 0, 1  (min, max)
#> 
#> 

plot(getVolumeUD(leroydbbmm))


## e.g. select the raster corresponding to the 95% UD
leroyUD <- getVolumeUD(leroydbbmm)
leroyUD[leroyUD>0.95] <- NA
plot(leroyUD)