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] 55.973192 39.737797  4.289011
#> 
#> $vectors
#>       azimuth    plunge
#> [1,] 128.2222 18.441302
#> [2,] 278.6912 69.030969
#> [3,]  34.9778  9.632663
#> 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"))