pro delta_image,file1,file2

; created 8/5/97, csee
; This program shows the difference betweent two images (file1-file2)

set_plot,'x'
d_read,file1,h1,p1
d_read,file2,h2,p2

n1=name(file1)
n2=name(file2)
t1=ccdtemp(file1)
t2=ccdtemp(file2)
e1=etime(file1)
e2=etime(file2)
l1=lamps(file1)
l2=lamps(file2)
mt1=mtime(file1)
mt2=mtime(file2)
mdl1=model(file1)
mdl2=model(file2)

pmax1=max(p1)
pmin1=min(p1)
pave1=mean(p1)
pmax2=max(p2)
pmin2=min(p2)
pave2=mean(p2)

p=p1-p2
pmax=max(p)
pmin=min(p)
pave=mean(p) 

print,'                                        exposure     ccdtemp    mission time  max.',$ 
  '    min.     average'

title1='Min., Max., Ave.: '+sstr(pmin1)+', '+sstr(pmax1)+', '+sstr(pave1)+ $
 ':  '+sstr(pmin2)+', '+sstr(pmax2)+', '+sstr(pave2)+ $
 ':  '+sstr(pmin)+', '+sstr(pmax)+', '+sstr(pave)
window,0,title=title1
print,n1,' from ',mdl1,', lamps= ',l1,e1,t1,mt1,pmax1,pmin1,pave1
tvscl,p1,10,10
;wait,1


;title2='Second Image: '+n2+' from '+mdl2+', lamps= '+l2+' exp= '+sstr(e2)+ $
;  ', temp= '+sstr(t2)+', mission time= '+sstr(mt2)+', Max.= '+sstr(pmax2)+', Min.= '+$
;  sstr(pmin2)+', Ave.= '+sstr(pave2)
;window,1,title=title2
print,n2,' from ',mdl2,', lamps= ',l2,e2,t2,mt2,pmax2,pmin2,pave2
tvscl,p2,210,10
;wait,1


;title3='Image1 - Image2:  Max.= '+sstr(pmax)+', Min.= '+$
;  sstr(pmin)+', Ave.= '+sstr(pave)
;window,2,title=title3
print,'Difference:                                                              ', $
  pmax,pmin,pave
tvscl,p,410,10

print,''
ans=''
read,'Do you want a printout (y or n)? ',ans
if ans eq 'y' then begin
 set_plot,'ps'
 device,filename = '~/idl.ps',ysize=7.5,xsize=10.0,/palatino,/landscape,/inches,bits=8

 close,3 
 openw,3,'temp3.txt'
 printf,3,format='(//)'
 printf,3,'File1= ',file1
 printf,3,'File2= ',file2
 printf,3,format='(T39,"exposure",2x,"ccdtemp",2x,"mission t",4x,"Max.",7x,"Min.",5x,"Ave.")'
 printf,3,n1,mdl1,l1,e1,t1,mt1,pmax1,pmin1,pave1,format='(a," from ",a,", lamps= ",a,T36,6f10.2)'
 printf,3,n2,mdl2,l2,e2,t2,mt2,pmax2,pmin2,pave2,format='(a," from ",a,", lamps= ",a,T36,6f10.2)'
 printf,3,pmax,pmin,pave,format='("Difference: ",T66,3f10.2)' 
 printf,3,'' 
 close,3
 spawn,'lp -ofp14  temp3.txt'
 spawn,'rm temp3.txt'
 
 tvscl,p1
 ;;xyouts,n1+' from '+mdl1+', lamps= '+l1+' exp= '+sstr(e1)+ $
 ;; ', temp= '+sstr(t1)+', mission time= '+sstr(mt1)+', Max.= '+sstr(pmax1)+', Min.= '+$
 ;; sstr(pmin1)+', Ave.= '+sstr(pave1), charsize=.8, alignment=
 erase
 
 tvscl,p2
 ;;xyouts,n2+' from '+mdl2+', lamps= '+l2+' exp= '+sstr(e2)+ $
 ;; ', temp= '+sstr(t2)+', mission time= '+sstr(mt2)+', Max.= '+sstr(pmax2)+', Min.= '+$
 ;; sstr(pmin2)+', Ave.= '+sstr(pave2), charsize=.8
 erase
 
 tvscl,p
;; xyouts,'Image1 - Image2:  Max.= '+sstr(pmax)+', Min.= '+$
 ;; sstr(pmin)+', Ave.= '+sstr(pave)+',  File1= '+file1+',  File2= '+file2, charsize=.8
 erase

 ;!p.title='First Image'
 ;show3,p1

 ;!p.title='Second Image'
 ;show3,p2

 !p.title='First Image - Second Image'
 show3,p

 device,/close_file
 spawn,'lp ~/idl.ps'
 set_plot,'x'

endif

exit:
;stop
end