Skip to contents

Decomposition of Orientation Tensor Eigenvectors and Eigenvalues

Usage

ot_eigen(x, scaled = FALSE)

Arguments

x

Object of class "Vec3", "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

Examples

set.seed(20250411)
mu <- rvmf(n = 1)
x <- rfb(100, mu = mu, k = 1, A = diag(c(10, 0, 0)))
x_eigen <- ot_eigen(x)
x_eigen
#> eigen() decomposition
#> $values
#> [1] 52.511342 42.798819  4.689839
#> 
#> $vectors
#> Vector (Vec3) object (n = 3):
#>              x          y          z
#> [1,] 0.4918288 -0.2979796 -0.8181153
#> [2,] 0.8109276  0.4988802  0.3058022
#> [3,] 0.3170187 -0.8138346  0.4870035
#> 
plot(x, col = "grey")
points(mu, col = 4)
text(mu, labels = "Mean", col = 4, pos = 4)
points(x_eigen$vectors, col = c(1, 2, 3))
text(x_eigen$vectors, col = c(1, 2, 3), labels = c("E1", "E2", "E3"), pos = 4)