とりあえず 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(…