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] 58.211774 36.621408  5.166818
#> 
#> $vectors
#>       azimuth   plunge
#> [1,] 128.9529 10.92743
#> [2,] 258.1261 73.00448
#> [3,]  36.4276 12.85540
#> 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"))