A tip on that... I don't know if you use matplotlib.pyplot, but if you do I recommend not using `imshow` - it plots images as if they were Matrices, as in with columns going Left->Right and rows going *down* from Top->Bottom by default.

Instead, I prefer `pcolor`, which produces output that is oriented in a FITS compliant sense.  E.g. if you have something in a galsim `Image` object called `image`, then the code below...
```
import matplotlib.pyplot as plt
plt.pcolor(image.array)
plt.colorbar()
image.write('array.fits')
```
should plot the contents of `image` to screen in a pretty easy way, with a colorbar, and when you load up `array.fits` with DS9 it will be oriented in the same manner.  This is sometimes handy for sanity checking without always having to write to FITS.