;
;
;

flist=file_search('20050325/speed_test8hz/spd*.fit',count=numfiles)
;flist=file_search('20050325/speed_test10hz/spd*.fit',count=numfiles)

window,0,xs=512,ys=512,retain=2
;
filetimes=fltarr(numfiles)
counts=fltarr(numfiles)
counts2=fltarr(numfiles)
itime=fltarr(numfiles)
othertemp=strarr(numfiles)
first_img=fltarr(1024,1024)
data=fltarr(1024,1024,100)

for i=300,399 do begin
 print,flist(i)
 r=fitshline(file=flist(i),error=e,KEY='ITIME')
 itime(i)=float(strmid(r,11,8))
 r=fitshline(file=flist(i),error=e,KEY='OTHERT')
 othertemp(i)=r
 temp=readfits(flist(i))
 data(*,*,i-300)=temp
 tvhist,rebin(temp,512,512)
 counts(i)=mean(temp)
 counts2(i)=mode(temp)
 if(i eq 0) then first_img=temp
endfor

noise=fltarr(1024,1024)
signal=fltarr(1024,1024)

for i=0,1023 do begin
 for j=0,1023 do begin
  noise(i,j)=stdev(data(i,j,*))
  signal(i,j)=mean(data(i,j,*))
 endfor
endfor
 

the_end:
end