Astropy currently only supports tiling along the first two dimensions but we can try and generalize this - I did it in #49 in the following commits:
e0533e5
242f7d6
but we decided to probably remove it for the first PR. We can re-do this as a follow-up PR.