Decomposition of Orientation Tensor Eigenvectors and Eigenvalues
See also
ortensor()
, eigen()
Examples
set.seed(20250411)
mu <- rvmf(n = 1)
x <- rfb(100, mu = mu, k = 1, A = diag(c(10, 0, 0)))
x_eigen <- eigen(x)
x_eigen
#> eigen() decomposition
#> $values
#> [1] 51.78685 42.84371 5.36944
#>
#> $vectors
#> Vector (Vec3) object (n = 3):
#> x y z
#> [1,] -0.5187265 0.3210950 0.7923514
#> [2,] 0.8037616 0.4990043 0.3239783
#> [3,] 0.2913590 -0.8049178 0.5169306
#>
plot(x, col = "grey")
points(mu, labels = "mu", col = 4)
#> Warning: "labels" is not a graphical parameter
points(x_eigen$vectors, col = c(1, 2, 3), labels = c("E1", "E2", "E3"))
#> Warning: "labels" is not a graphical parameter