aslopect {simba}R Documentation

Calculate similarity of plots based on slope aspect and inclination

Description

Allows for the comparison of plots regarding the two variables slope aspect and slope inclination at once. To obtain a distance measure integrating aspect and inclination the model of a unit sphere is used and great-circle distances between virtual locations are calculated. For each plot a virtual location on the sphere is defined using the values for aspect as longitude and 90°-inclination as latitude. See details for more...

Usage

aslopect(asp, slo, names=rownames(asp), listout = FALSE)

Arguments

asp Numeric vector with aspect values, given in degree. Expects values between 0° (North) and 180°/-180° (South). Eastward directions count positive, westward directions count negative.
slo Numeric vector with slope inclination values. Expects values between 0°(flat) and 90° (vertical wall)
names Plot names, defaults to the rownames of asp, but a separate vector can be specified. Its length has to match the length of asp and slo.
listout Shall the results be given in list-format (data.frame). Defaults to FALSE which means that a matrix of class dist is returned

Details

To obtain a distance measure integrating aspect and inclination the model of a unit sphere is used and great-circle distances between virtual locations are calculated. For each plot a virtual location on the sphere is defined using the values for aspect as longitude and 90°-inclination as latitude. This means that as long as the inclination is low the virtual points are located in the pole region so that, regardless of aspect, plots with low inclination are rather close to each other regarding these qualities. The idea behind is, that solar radiation, wind or other factors highly depending on aspect and inclination are not really differing between plots with different aspect as long as the slope is low. The longitude values on the unit sphere are derived from the values of slope aspect. The equator of the sphere is thought as the compass circle. The Prime Meridian of the virtual sphere is the great circle through North and South of the compass. As in geographic terms longitude counts positive in Eastern and negative in Western direction. With phi = latitude = 90°-inclination and λ = longitude = aspect the great-circle distance between A and B can be calculated with the following formula.

sim = arccos(sin(phiA)*sin(phiB) + cos(phiA)*cos(phiB)*cos(lambdaB-lambdaA))

Value

Returns a dist object or a data.frame (depending if listout = FALSE or TRUE). As a unit sphere is used, the maximum distance between two inclination/aspect pairs is perimeter/2 of the sphere which is by definition Pi. To scale the possible distances between 0 and 1 the results are divided by Pi. Thus, a great-circle distance of 1 is rather scarce in the real world, however, two vertical rock walls with opposite aspect would share it. If listout = TRUE a data.frame with the following variables returns.

NBX one of the compared sampling units
NBY the other part of the pair
x The returned aslopect value

Author(s)

Gerald Jurasinski

Examples

data(abis)
## identify columns with slope and aspect data
names(abis.env)

## calculate aslopect
abis.aslop <- aslopect(abis.env[,4], abis.env[,5])


[Package simba version 0.2-5 Index]