volume: retrieve a set of slices along any dimension

It would be convenient to have a util function that can retrieve a set of slices from (axis, indices)