Function for_each_slice 

fn for_each_slice(
    data: &mut [f64],
    dims: [usize; 3],
    axis: usize,
    f: impl FnMut(usize, &mut f64),
)

Apply f(slice_index_along_axis, &mut value) to every voxel of a C-order volume.