Use astropy.io.fits indexing convention (y, x rather than x, y)