pub fn rotation_3d(radians: f64, ax: Axes) -> Array2<f64>
generates a 3x3 matrix describing a rotation around either the index coordinate axes (X,Y) or in the perpendicular axes to the image (Z)