Utilization distribution (UD)
getVolumeUD.Rd
The UD represents the minimum area in which an animal has some specified probability of being located (Cumulative Distribution Function).
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.
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)