barf=1 if(barf eq 0) then begin ;{ filenmO='980226_181138_eclp.fts' old = readfits(filenmO,hdr, SILENT=1) filenmN='980226_181151_eclp.fts' ; filenmN='980226_181204_eclp.fts' ; filenmN='980226_181216_eclp.fts' new = readfits(filenmN,hdr, SILENT=1) byteorder,old byteorder,new ;; Convert the signed integers to unsigned integers. imgN = float(new) imgN = (imgN lt 0.) * ( imgN + 2.*32768. ) + (imgN ge 0.)*(imgN) imgO = float(old) imgO = (imgO lt 0.) * ( imgO + 2.*32768. ) + (imgO ge 0.)*(imgO) endif else begin;}{ filenmO='CALIBRATED/mp1' filenmN='CALIBRATED/mp2' ; filenmN='CALIBRATED/mp3' ; filenmN='CALIBRATED/mp4' imgN=fltarr(2034,2034) imgO=fltarr(2034,2034) openr,1,filenmN & readu,1,imgN & close,1 openr,1,filenmO & readu,1,imgO & close,1 endelse;} ;window,xsize=2034,ysize=2034,xpos=-400,ypos=-400 diff = imgN-imgO mn= -5. mx= 5. tv,bytscl(diff, mn,mx) xbeg= 910 & ybeg= 890 xend= 1200 & yend= 1200 leastsq,imgN(xbeg:xend,ybeg:yend),imgO(xbeg:xend,ybeg:yend),aa,bb,RR print,filenmO,filenmN,aa,bb,RR ;window,1 wset,1 plot,imgN(xbeg:xend,ybeg:yend),imgO(xbeg:xend,ybeg:yend),/ynoz,psym=1 wset,0 end