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 としておきます。