I've been reading the wikipedia for 4D rotations. I understand the algebraic argument for why any rotation matrix can be decomposed into a right-isoclinic and a left-isoclinic rotation, and why these correspond to right and left multiplication by unit quaternions. But how does one go from a unit quaternion to the planes and the angle of the corresponding isoclinic rotation?
In 3D, I have the nice and easy fact that the imaginary part of a unit quaternion specifies the axis (or plane) of rotation while the real part is the cosine of one half times the rotation angle. Is there an analogously simple way of thinking about it in 4D?
I tried reasoning this out but ran into a wall almost immediately: specifying a unit quaternion is the same as specifying a line through the origin of $\mathbb{R}^4$ as well as a hyperplane orthogonal to this line, but I have no idea how one could go from this information to a pair of orthogonal 2-planes plus an angle.