"Cannot infer original voxel size for reslicing... "
"set mri_resolution to boolean value or apply morph first.")
// Now deal with the fact that we may have multiple sub-volumes
inuse = [morph[k]["inuse"] for k in range(len(morph))]
src_shape = [morph[k]["shape"] for k in range(len(morph))]
assert len(set(map(tuple, src_shape))) == 1
morph = BunchConst(src_data=_get_src_data(morph, mri_resolution)[0])