ipython -pylab
help(pylab)
import pyfits
hdulist = pyfits.open('ngc6302.fits')
myhead = hdulist[0].header
scidata = hdulist[0].data
hdulist.close()
myhead.ascardlist()
help(imshow)
imshow(scidata)