Function crop_c_order 

fn crop_c_order<T: Copy + Default>(
    data: &[T],
    dims: [usize; 3],
    bounds: [(usize, usize); 3],
    n_ch: usize,
) -> Vec<T>

Crop a C-order [d0, d1, d2] volume of n_ch-element voxels to bounds.