Convert raster to contour lines
raster2contour.RdThe function converts a raster UD(stack) object to a SpatialLinesDataFrame. This allows to re-project the contours to different projections.
Arguments
- x
a
DBBMM,DBBMMStack,dynBGB,.UDor.UDStackobject- ...
additional arguments, like
levelsandnlevels, that can be passed to 'rasterToContour' function
Details
The contour function creates a shape of the area in which the animal can be found by a certain probability (i.e. the 90% contour describes the area in which the animal can be found with the 90% probability).
One or several probabilities can be set with levels (numeric or vector of values between 0 and 1). If no value is set all contour lines are returned.
You can also use nlevel to set a number of fixed distance levels.
The raster2contour function creates a SpatialLinesDataFrame from the input raster object. This allows to re-project the contours to different projections.
Examples
data(leroydbbmm)
data(leroydbgb)
data(dbbmmstack)
## from a DBBMM object
(cont1 <- raster2contour(leroydbbmm))
#> geometry level
#> C_1 MULTILINESTRING((-1518.966 287.8897 ...)) 0.1
#> C_2 MULTILINESTRING((-1518.966 195.3691 ...)) 0.2
#> C_3 MULTILINESTRING((-2278.934 1786.655 ...)) 0.3
#> C_4 MULTILINESTRING((-2278.934 1624.06 ...)) 0.4
#> C_5 MULTILINESTRING((-2278.934 1461.464 ...)) 0.5
#> C_6 MULTILINESTRING((-2278.934 1298.869 ...)) 0.6
#> C_7 MULTILINESTRING((-3038.901 1809.509 ...)) 0.7
#> C_8 MULTILINESTRING((-3038.901 1434.866 ...)) 0.8
#> C_9 MULTILINESTRING((-3038.901 975.3713 ...)) 0.9
#> C_10 MULTILINESTRING((-6078.772 -2659.461 ...)) 1
plot(cont1)
## from a dynBGB object
(cont2 <- raster2contour(leroydbgb, level=.95))
#> geometry level
#> C_1 MULTILINESTRING((7.56312 -140.8608 ...)) 0.95
plot(cont2)
## from a DBBMMStack object
(cont3 <- raster2contour(dbbmmstack))
#> geometry level
#> C_1 Leroy MULTILINESTRING((-616.1418 -2498.635 ...)) 0.1
#> C_2 Leroy MULTILINESTRING((-2467.477 -3724.124 ...)) 0.2
#> C_3 Leroy MULTILINESTRING((-2467.477 -3803.409 ...)) 0.3
#> C_4 Leroy MULTILINESTRING((-2467.477 -3882.693 ...)) 0.4
#> C_5 Leroy MULTILINESTRING((-3084.588 -2005.183 ...)) 0.5
#> C_6 Leroy MULTILINESTRING((-3084.588 -2201.046 ...)) 0.6
#> C_7 Leroy MULTILINESTRING((-3084.588 -2396.908 ...)) 0.7
#> C_8 Leroy MULTILINESTRING((-3084.588 -2796.651 ...)) 0.8
#> C_9 Leroy MULTILINESTRING((-3084.588 -3338.411 ...)) 0.9
#> C_10 Leroy MULTILINESTRING((-6170.146 -6787.781 ...)) 1
#> C_1 Ricky.T MULTILINESTRING((618.0813 6702.374 ...)) 0.1
#> C_2 Ricky.T MULTILINESTRING((618.0813 6616.074 ...)) 0.2
#> C_3 Ricky.T MULTILINESTRING((618.0813 6529.774 ...)) 0.3
#> C_4 Ricky.T MULTILINESTRING((618.0813 6443.475 ...)) 0.4
#> C_5 Ricky.T MULTILINESTRING((618.0813 6357.175 ...)) 0.5
#> C_6 Ricky.T MULTILINESTRING((618.0813 6270.876 ...)) 0.6
#> C_7 Ricky.T MULTILINESTRING((618.0813 6184.576 ...)) 0.7
#> C_8 Ricky.T MULTILINESTRING((618.0813 5982.249 ...)) 0.8
#> C_9 Ricky.T MULTILINESTRING((618.0813 5759.316 ...)) 0.9
#> C_10 Ricky.T MULTILINESTRING((-2467.477 5554.45 ...)) 1
#> individual.local.identifier
#> C_1 Leroy Leroy
#> C_2 Leroy Leroy
#> C_3 Leroy Leroy
#> C_4 Leroy Leroy
#> C_5 Leroy Leroy
#> C_6 Leroy Leroy
#> C_7 Leroy Leroy
#> C_8 Leroy Leroy
#> C_9 Leroy Leroy
#> C_10 Leroy Leroy
#> C_1 Ricky.T Ricky.T
#> C_2 Ricky.T Ricky.T
#> C_3 Ricky.T Ricky.T
#> C_4 Ricky.T Ricky.T
#> C_5 Ricky.T Ricky.T
#> C_6 Ricky.T Ricky.T
#> C_7 Ricky.T Ricky.T
#> C_8 Ricky.T Ricky.T
#> C_9 Ricky.T Ricky.T
#> C_10 Ricky.T Ricky.T
plot(cont3)
(cont4 <- raster2contour(dbbmmstack, level=c(.5,.95)))
#> geometry level
#> C_1 Leroy MULTILINESTRING((-3084.588 -2005.183 ...)) 0.5
#> C_2 Leroy MULTILINESTRING((-3701.7 -2805.881 ...)) 0.95
#> C_1 Ricky.T MULTILINESTRING((618.0813 6357.175 ...)) 0.5
#> C_2 Ricky.T MULTILINESTRING((618.0813 5647.85 ...)) 0.95
#> individual.local.identifier
#> C_1 Leroy Leroy
#> C_2 Leroy Leroy
#> C_1 Ricky.T Ricky.T
#> C_2 Ricky.T Ricky.T
plot(cont4)
