PyFITS provides an interface to FITS formatted files in the Python
scripting language. PyFITS is a development project of the Science
Software Branch at the Space Telescope Science Institute. 

This requires numpy.
