;
;
;

flist=file_search('20050331/s*test*.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)

for i=0,numfiles-1 do begin
 print,flist(i)
 r=fitshline(file=flist(i),error=e,KEY='TIME_SYS')
 itime(i)=hms2h('0'+strmid(r,11,7));+float(strmid(r,19,3))/(1000.*3600.)
 print,i,itime(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)=temp
 tvhist,rebin(temp,512,512)
 wait,0.1
 counts(i)=mean(temp)
; counts2(i)=mode(temp)
; if(i eq 0) then first_img=temp
endfor

;c=poly_fit(findgen(numfiles),itime,1)
;print,1./(c(1)*3600.),' frames per second'

the_end:
end