Function resample_to_axial 

pub fn resample_to_axial(
    data: &[f64],
    nx: usize,
    ny: usize,
    nz: usize,
    affine: &[f64; 16],
) -> ResampledVolume

Resample a continuous volume (magnitude, unwrapped field, χ) to a cardinal-aligned grid.

Not for wrapped phase — see the module docs and resample_complex_to_axial.