;
; nim array light transfer tests
;
;goto,point1

intframe=fltarr(256,256,10)
noiseframe=fltarr(256,256,10)
x=[1,2,3,6,10,18,30,56,100,178]

for i=0,9 do begin
 print,i
; dark=fltarr(256,256)
 temp=fltarr(256,256,32)
 for j=0,31 do begin
;  dark=dark+float(swap_endian(rfits('c:\data\02sep25\light_xfer_'+string(x(i),format='(i3.3)')+'.d.i.'+string(j+1,format='(i3.3)'))))
  temp(*,*,j)=(swap_endian(rfits('/usr2/mpenn/02sep25/lamp2/light_xfer_'+string(x(i),format='(i3.3)')+'.o.i.'+string(j+1,format='(i3.3)'))))
 endfor
; dark=dark/32.
; temp2=temp
; for k=0,31 do temp2(*,*,k)=temp2(*,*,k)-mean(temp2(*,*,k))
 for k=0,255 do begin
  for kk=0,255 do begin
   intframe(k,kk,i)=mean(temp(k,kk,*))
   noiseframe(k,kk,i)=stdev(temp(k,kk,*))
  endfor
 endfor
endfor

point1:

xx=alog(intframe+1.e-5)/alog(10.)
yy=alog(noiseframe+1.e-5)/alog(10.)


xbin=findgen(35)/10.+0.05
ybin=fltarr(35)
yerr=fltarr(35)

for i=0,34 do begin
 b=where(xx gt xbin(i)-0.05 and xx le xbin(i),numb)
 if (numb gt 3) then begin
  ybin(i)=median(yy(b))
  yerr(i)=stdev(yy(b))
 endif else begin
  ybin(i)=0.
  yerr(i)=0.
 endelse
endfor

plot,xbin,ybin,psym=2,xtit='log!D10!N Intensity',ytit='log!D10!N Noise',chars=1.6,$
     xran=[0,4],xstyl=1,yran=[-0.5,1.0],ystyl=1
errplot,xbin,ybin-yerr,ybin+yerr
plots,[0,4],[0.096,0.096],linestyl=2
plots,[1.56,4],[-0.5,0.72],linestyl=2


;plot,[xx],[yy],psym=3,xtit='intensity',ytit='noise',tic=1,yran=[-.5,1.5],ysty=1,xran=[0,3.7],xsty=1

the_end:
end