Function b0_angle_from_affine 

pub fn b0_angle_from_affine(affine: &[f64; 16]) -> f64

Angle (degrees) between B0 and the voxel +z axis — the tilt that matters for the dipole kernel. Zero for an axial acquisition; equals the scanner’s slice tilt for a simple oblique.