pyfits と pyds9 で楽々生活

とりあえず ds9 に array を表示する．

import ds9, pyfits

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

activeds9list = ds9.ds9_targets()
if activeds9list is not None:
	view=ds9.ds9(activeds9list[0])
else:
	view=ds9.ds9()

def showregion(array):
     view.set('array [xdim=%d,ydim=%d,BITPIX=-32]' % \
          (len(array[0]),len(array)) , array[0:len(array),0:len(array)])

使い方

>>> showregion(hdulist[0].data)
>>> showregion(hdulist[0].data[0:100,0:100])