import pylab as pl

#help(pylab)
import pyfits

#import matplotlib as mp 

hdulist = pyfits.open('ngc6302.fits')

myhead = hdulist[0].header

scidata = hdulist[0].data

hdulist.close()

myhead.ascardlist() 


pl.imshow(scidata)