Eigenvalues and Eigenvectors of a Set of Vectors

Decomposition of Orientation Tensor Eigenvectors and Eigenvalues

Usage

or_eigen(x, scaled = FALSE)

Arguments

x

numeric. Can be three element vector, three column array, or an object of class "line" or "plane"

scaled

logical. Whether the Eigenvectors should be scaled by the Eigenvalues (only effective if x is in Cartesian coordinates).

Value

list containing

values

Eigenvalues

vectors

Eigenvectors in coordinate system of x

See also

ortensor(), eigen()

Examples

set.seed(1)
mu <- rvmf(n = 1) |> vec2line()
x <- rfb(100, mu = mu, k = 1, A = diag(c(10, 0, 0)))
x_eigen <- or_eigen(x)
x_eigen
#> eigen() decomposition
#> $values
#> [1] 60.136563 35.136468  4.726968
#> 
#> $vectors
#>        azimuth   plunge
#> [1,] 127.45652 25.49194
#> [2,] 282.15637 62.19332
#> [3,]  32.45254 10.36689
#> attr(,"class")
#> [1] "line"
#> 
stereoplot()
stereo_point(x, col = "grey")
stereo_point(mu, lab = "mu", col = 4)
stereo_point(x_eigen$vectors, col = c(1, 2, 3), lab = c("E1", "E2", "E3"))