Given a 4x4 orthonormal matrix with a determinant of 1, when interpreting this matrix as a rotation basis in 4D space, there exists a 4D geometric algebra rotor with scalar and bivector components (1 scalar number, 6 bivector numbers), and possibly a pseudoscalar, that describes the rotation of this matrix relative to a 4x4 identity matrix. How can this be computed?
The Wikipedia article about 4D rotations is completely missing information on geometric algebra rotors: https://en.wikipedia.org/wiki/Rotations_in_4-dimensional_Euclidean_space
The Wikipedia article about Quaternions and spatial rotation explains how to convert from a matrix, but does not explain how this generalizes to higher dimensions with rotors: https://en.wikipedia.org/wiki/Quaternions_and_spatial_rotation
This question is related to Recover a 4D Rotor in geometric algebra from rotated vectors, since a 4x4 matrix can be thought of as a set of 4 rotated column vectors. However, an operation on 2 vectors would only allow for getting the rotation from the identity to one of the columns, while this question is asking for a general solution of the entire matrix at once.
This question is related to Decomposition of a single 4D rotation, however this question is not looking for a solution involving Euler angles, rather geometric algebra rotors and bivectors. I actually already have an algorithm for converting a 4D basis into 4D Euler angles, but then I also have no way to convert that to a 4D rotor.