cmnd = 'ls -1 OPALS/opal_S_p?' ; cmnd = 'ls -1 OPALS/opal_M_p?' cmnd = 'ls -1 OPALS/opal25_M_p?' ; cmnd = 'ls -1 OPALS/opal_L_p?' if(n_elements(p1) eq 0 ) then begin ;{ p1=fltarr(2034,2034) p2=fltarr(2034,2034) p3=fltarr(2034,2034) p4=fltarr(2034,2034) window,xsize=339*3, ysize=339 endif ;} spawn,cmnd,files openr,1,files(0) & readu,1,p1 & close,1 openr,1,files(1) & readu,1,p2 & close,1 openr,1,files(2) & readu,1,p3 & close,1 openr,1,files(3) & readu,1,p4 & close,1 Tp1 = total(p1(50:50+1933,50:50+1933)) / (1934.*1934.) Tp2 = total(p2(50:50+1933,50:50+1933)) / (1934.*1934.) Tp3 = total(p3(50:50+1933,50:50+1933)) / (1934.*1934.) Tp4 = total(p4(50:50+1933,50:50+1933)) / (1934.*1934.) avgTp = ( Tp1 + Tp2 + Tp3 + Tp4 ) *0.25 Tp1 = Tp1 / avgTp Tp2 = Tp2 / avgTp Tp3 = Tp3 / avgTp Tp4 = Tp4 / avgTp TotI = p1 + p2 + p3 + p4 Q = p1 - p2 - p3 + p4 U = p1 + p2 - p3 - p4 donorm = 1 if ( donorm ) then begin ;{ Q = p1/Tp1 - p2/Tp2 - p3/Tp3 + p4/Tp4 U = p1/Tp1 + p2/Tp2 - p3/Tp3 - p4/Tp4 endif ;} Gain = total(TotI(50:50+1933,50:50+1933)) / (1934.*1934.) X21 = (total(Q (50:50+1933,50:50+1933)) / (1934.*1934.) ) / Gain X31 = (total(U (50:50+1933,50:50+1933)) / (1934.*1934.) ) / Gain print,cmnd,' ',Gain,X21,X31 mn=Gain-Gain*.1 mx=Gain+Gain*.1 smI=rebin(TotI,339,339) & tv,bytscl(smI, mn, mx), 0 ,0 smQ=rebin(Q ,339,339) & tv,bytscl(smQ, -150, 150),339 ,0 smU=rebin(U ,339,339) & tv,bytscl(smU, -150, 150),339*2,0 smX21= smQ / smI smX31= smU / smI x21mn = min(smX21(20:338-20,20:338-20),max=x21mx) x31mn = min(smX31(20:338-20,20:338-20),max=x31mx) tv,bytscl(smX21, x21mn , x21mx ),339 ,0 tv,bytscl(smX31, x31mn , x31mx ),339*2,0 print,x21mn,x21mx,' X21 mn,mx', x21mx - x21mn print,x31mn,x31mx,' X31 mn,mx', x31mx - x31mn end