Function pad_centred 

fn pad_centred<T: Copy + Default>(
    x: &[T],
    dims: (usize, usize, usize),
    new: (usize, usize, usize),
    pad: [usize; 3],
) -> Vec<T>

Pad with pad zero voxels on each side (STI’s padarray(x, padsize)).