前準備

まず用意するものを挙げます。

• ``あすか''の生イメージ (256x256 が望ましい)
• ``あすか''の point spread function のイメージ
• Deconvolution を IDL で行なうための Script

以下に、Deconvolution を IDL で行なうための Script の例を示します。 適当なファイルに書いてください。

     p=readfits(psffile)
     im=readfits(infile,imhdr)
     help, p
     help, im
     s=size(im)
     sp=size(p)
     i=(s(1)-sp(1))/2
     j=(s(2)-sp(2))/2
     img=im(i:sp(1)-1+i,j:sp(2)-1+j)
     help, img
     psf=p/total(p)
     cdpix = sxpar(imhdr, 'CRPIX*')
     sxaddpar, imhdr, 'NAXIS1', sp(1)
     sxaddpar, imhdr, 'NAXIS2', sp(2)
     sxaddpar, imhdr, 'CRPIX1', cdpix(0)-i
     sxaddpar, imhdr, 'CRPIX2', cdpix(1)-i
     print, imhdr
     maxlike2,img,psf,r,/init,niter=niter, header=imhdr,$
     outfile=outfile, chimin=0.0, stall=0.0, display='tvscl,iter'

仮に、これのファイルを deconvolve.idl としておきます。