Using this CartesianIndices, you can create a n-dimensional recetengular
shape and iterate over the region. This should be quite a powerful but
easy to use abstraction over the volume.

The provided free functions, can be used to visit neighbouring voxels
from a given position.