flow.dispersion_int: Compute flow-based cost or conductance

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/wind_functions2.R

Description

flow.dispersion_int computes movement conductance through a flow either, sea or wind currents. It implements the formula described in Felícisimo et al. 2008:

Usage

1
flow.dispersion_int(stack, fun = cost.FMGS, output = "transitionLayer", ...)

Arguments

stack

RasterStack object with layers obtained from wind2raster function ("rWind" package) with direction and speed flow values.

fun

A function to compute the cost to move between cells. The default is cost.FMGS from Felicísimo et al. (2008), see details.

output

This argument allows to select different kinds of output. "raw" mode creates a matrix (class "dgCMatrix") with transition costs between all cells in the raster. "transitionLayer" creates a TransitionLayer object with conductance values to be used with "gdistance" package.

...

Further arguments passed to or from other methods.

Details

Cost=(1/Speed)*(HorizontalFactor)

being HorizontalFactor a "function that incrementally penalized angular deviations from the wind direction" (Felicísimo et al. 2008).

Value

In "transitionLayer" output, the function returns conductance values (1/cost)to move between all cells in a raster having into account flow speed and direction obtained from wind.fit function("rWind" package). As wind or sea currents implies directionality, flow.dispersion produces an anisotropic conductance matrix (asymmetric). Conductance values are used later to built a TransitionLayer object from "gdistance" package.

In "raw" output, flow.dispersion creates a sparse Matrix with cost values.

Note

Note that for large data sets, it could take a while. For large study areas is strongly advised perform the analysis in a remote computer or a cluster.

Author(s)

Javier Fernández-López; Klaus Schliep; Yurena Arjona

References

Felicísimo, Á. M., Muñoz, J., & González-Solis, J. (2008). Ocean surface winds drive dynamics of transoceanic aerial movements. PLoS One, 3(8), e2928.

Jacob van Etten (2017). R Package gdistance: Distances and Routes on Geographical Grids. Journal of Statistical Software, 76(13), 1-21. doi:10.18637/jss.v076.i13

See Also

wind.dl, wind2raster

Examples

1
2
3
4
5
6
7
8
9
data(wind.data)
wind <- wind2raster(wind.data)
Conductance <- flow.dispersion(wind, type = "passive")
## Not run: 
require(gdistance)
transitionMatrix(Conductance)
image(transitionMatrix(Conductance))

## End(Not run)

rWind documentation built on Oct. 19, 2021, 5:06 p.m.