This function identifies movement track segments whose attributes suggest corridor use behavior
1 2 3 4
numeric between 0 and 1, defines the proportion of speeds which are high enough to be a valid corridor point (default value selects speeds that are greater than 75 % of all speeds).
numeric between 0 and 1, defines the proportion of the circular variances that is low enough to be a valid corridor point. Low values of the circular variance indicate that the segments are (near) parallel (default value selects variances that are lower than 25 % of all variances).
logical, if TRUE the track is plotted together with dots that indicate corridor points when a
The corridor function uses the attributes of a movement step (segment) to identify movement steps that exhibit corridor use behavior. For each segment, the speed and the azimuth are calculated and assigned to the segment midpoint.
A circular buffer is created around the midpoint of each segment whose radius is equal to half the segment length. The segment azimuth (180 >= azimuth > -180) is converted into a new unit, the 'pseudo-azimuth' (0 <= 360), removing the directional information.
Subsequent, the circular variance of the pseudo-azimuths of all segment midpoints that fall within the circular buffer is calculated. Low values of the circular variance indicate that the segments are (near) parallel.
Next, it is determined whether a segment's speed is higher than
speedProp (by default the upper 25% speeds) and its circular variance is lower than
circProp (by default the lower 25% of all variances).
Segment midpoints that meet both of these requirements are considered as a 'corridor' point, all others are considered 'non-corridor' points. Finally, a corridor point is determined to be within a true corridor if within its circular buffer there are more 'corridor' points than 'non-corridor' points.
The function returns a
moveBurst object or a list of moveBurst objects (if a MoveStack is supplied). The MoveBurst date.frame stores the following information:
- segment midpoint
- circular variance
The object is bursted by the factor that indicates whether the segment belongs to a corridor segment or not, and is specified in the "burstId" slot.
The default values for the
circProp can be changed by the users discretion using the according argument.
Marco Smolla & Anne Scharf
LaPoint, S., Gallery, P., Wikelski, M. and Kays, R. (2013), Animal Behavior, Cost-based Corridor Models, and Real Corridors. Landscape Ecology. doi:10.1007/s10980-013-9910-0.
1 2 3 4 5 6 7 8 9 10
## with a move object data(leroy) tmp <- corridor(leroy, plot=TRUE) plot(tmp, type="l", col=c("red","black")[c(tmp@burstId,NA)]) ## with a moveStack object data(fishers) stacktmp <- corridor(fishers[c(1:400,sum(n.locs(fishers))-(400:1)),]) plot(stacktmp[], col=c("red","black")[stacktmp[]@burstId]) lines(stacktmp[], col=c("red","black")[c(stacktmp[]@burstId,NA)])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.