Decomposition of Orientation Tensor Eigenvectors and Eigenvalues
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)