pro jpegtriplets

;This program makes a jpeg file of the image that you select
; (along with its corrispoding triplet mates). 
; The output file is: ~/triplet.jpg


cd,'',current=cdir & cd,cdir
if n_params() eq 0 then log = cdir  ;current log is the default
cd,'DB/Image'

print,''
print,'The Test Log is: ',log

;____________________________________________________________________________
;List Files...

spawn,'lsf',files,/noshell
s=size(files)
nfiles=fix(s(1))
numfiles=sstr(nfiles)				;numfiles=short string (nfiles)
if files(0) eq '' then begin & print,'No files found' & goto,exit & endif
array_sort,files	;sorts files in assending order

print,'The number of Images in this log is: ',sstr(nfiles),' (1 thru ',sstr(nfiles),')'
window,1,xsize=250, ysize=850

;____________________________________________________________________________
;Initalize...

auto=1		;for auto=1 tvscl
min=0		;the tvscl min setting
max=4095	;the tvscl max setting
ilast=0		;ilast is the previous value of i
space=1.5	;x spacing in inches
sf=.75		;scale factor (scales the size of the images on the page)
ypos=0.5	;sets the y position for each row
ff=-1		;flag for flat fields, 1=>enabled, -1=>disabled

menu:
print,''
print,'Enter  to...              Enter  to...
print,'  0    Get this menu       -4    Autoscale, min to max
print,' -1    Quit the program    -5    Scale from min to 2 x average
print,' -2    Change minimum      -6    Print current image
print,' -3    Change maximum      -7    Toggle flat field

;____________________________________________________________________________
;begin the loop that does the deed...

again:
file_no=0
read,'Which triplet do you want to display, 1 thru '+sstr(nfiles)+'? (0 for menu): ',file_no
 
if file_no eq 0 then goto,menu	;display menu

if file_no eq -1 then goto, exit	;end of program
 
 if file_no eq -2 then begin		;set minimum of stretch
  auto=0
  read,'Min? ',min
  print,'Min set to: ',min
  file_no=ilast
 endif
 
 if file_no eq -3 then begin		;set maximum of stretch
  auto=0
  read,'Max? ',max
  print,'Max set to: ',max
  file_no=ilast
 endif
 
 if file_no eq -4 then begin		;autoscale to max and min
  auto=1  
  print,'Set to Autoscale, min to max.'
  file_no=ilast
 endif 
 
 if file_no eq -5 then begin		;autoscale, range = minimum to 2 x average
  auto=2  
  print,'Set to Autoscale, min to 2xAve.'
  file_no=ilast
 endif 
 
 if file_no eq -6 then begin		;prints current image
  set_plot,'ps'  
  device,filename = '~/idl.ps',/times,/portrait,xsize=7.5,ysize=9,yoffset=1.0,/inches,bits_per_pixel=8
  ;;device,filename = '~/idl.ps',/times,/portrait,xsize=4.5,ysize=5.4,yoffset=1.0,/inches,bits_per_pixel=8
  if auto eq 1 then tvscl,bytscl(p,min=min(p),max=max(p),top=255) else tvscl,bytscl(p,min=min,max=max,top=255)	;p>min<max
  xyouts,0,-.02,'triplet'+sstr(ilast)+'  Min= '+sstr(min(p))+', Ave= '+sstr(mean(p))+', Max= '+sstr(max(p)),/norm
  xyouts,0,-.04,'Exposure (ms) = '+sstr(d_value(h,36))+'   lamps = '+sstr(d_value(h,102)),/norm
  xyouts,0,1.02,'From Log: '+cdir,/norm
  device,/close_file
  spawn,'lp ~/idl.ps'
  set_plot,'x'		;returns to X windows environment
  file_no=ilast
 endif 
 
 if file_no eq -7 then begin		;use flat fields
  ff=-ff  
  if ff eq 1 then print, 'flat field enabled' else print,'flat field disabled'
  file_no=ilast
 endif 
 
 
 if file_no le 0 or file_no gt nfiles then begin	;for i input not in range
  if file_no le -8 or file_no gt nfiles then print,'Invalid value'
  goto,menu
 endif
 



;____________________________________________________________________________
;find three matching files...
   filetime=strmid(files(file_no-1),13,13)
   ;stop
   triplets=t_getfil('*'+filetime+'*') & siz=size(triplets) & nfiles2=siz(1)	
   if nfiles2 ne 3 then print,'Not a complete triplet'
   
   array_sort,triplets		;sorts files in assending order

   d_read,triplets(0),h,p0	;MRI 176x256 or 254
   d_read,triplets(1),h,p1	;SLI 128x256 or 254
   d_read,triplets(2),h,p2	;HRI 160x256 or 254

;____________________________________________________________________________
;apply flat field?...
 
model=d_value(h,219)	;the model is needed to get the right flat field...

if ff eq 1 then begin
 
 case model of
  ;;'DISR1': model='DISR#1'	;DISR1 Flat field is wrong
  'DISR2': model='DISR#2'
  'DISR3': model='DISR#3'
  else: model='Unknown'
  endcase
  
  if model eq 'Unknown' then begin
   print,'There is no flat field for model: ',d_value(h,219)
   ff=-1
   print,' Flat Field Disaabled'
   goto,menu
  endif
   
  d_read,'/local/Imagers/Flat_Field/'+model+'/dli2_flat_field.disrsoft',hff,pff ;MRI
  P0=p0/pff
  
  d_read,'/local/Imagers/Flat_Field/'+model+'/sli_flat_field.disrsoft',hff,pff	;SLI
  P1=p1/pff

  d_read,'/local/Imagers/Flat_Field/'+model+'/dli1_flat_field.disrsoft',hff,pff	;HRI
  P2=p2/pff

endif
  
;____________________________________________________________________________
;scale so that their means are equal ...
;stop

q=max([mean(p0),mean(p1),mean(p2)])
 p0=p0*q/mean(p0)
 p1=p1*q/mean(p1)
 p2=p2*q/mean(p2)

;;p1=bytscl(p1,min=min(p1),max=2*mean(p1),top=255)	;SLI
;;p2=bytscl(p2,min=min(p2),max=2*mean(p2),top=255)	;HRI
;;p0=bytscl(p0,min=min(p0),max=2*mean(p0),top=255)	;MRI 


;____________________________________________________________________________
;combine, display & create jpeg...

p=intarr(176,768)
s=size(p0)
hight=s(2)
p(24:151,0:hight-1)=p1		;SLI
p(0:175,hight:2*hight-1)=p0	;MRI
p(8:167,2*hight:3*hight-1)=p2	;HRI

 p=rotate(p,2)				;right side up

 if auto eq 2 then begin & min=min(p) & max=2*mean(p) & endif
 if auto eq 1 then begin & min=min(p) & max=max(p) & endif


   ;tvscl,bytscl(p,min=min(p),max=2*mean(p),top=255)
   tvscl,p>min<max
   
 print,'file '+sstr(file_no),'  Min= ',sstr(min(p)),', Ave= ',sstr(mean(p)),', Max= ',sstr(max(p))
 print,'Exposure (ms) = ',sstr(d_value(h,36)),'   lamps = ',sstr(d_value(h,102))
 if ff eq -1 then print,'Model= ',model,'  Flat Field is off' $
  else print,'Model= ',model,'  Flat Field is on'
 
   for n=0,5 do begin
   
   pout=bytscl(p,min=min(p),max=2*(.7+n*.1)*mean(p),top=255)
   outfile="~/triplet1."+sstr(n)+".jpg"
   write_jpeg,outfile,pout,quality=100
   
  pout2=bytscl(p,min=min(p),max=2*(.8+n*.1)*mean(p),top=4095)
   minpt=min(pout2)>0
   maxpt=max(pout2)
   display_in=uint(65535d0*(double(pout2)-double(minpt))/(double(maxpt)-double(minpt)))
   storfile="~/triplet2."+sstr(n)+".png"
   write_png,storfile,display_in
   
  ;;pout2=bytscl(p,min=min(p),max=2*(.8+n*.1)*mean(p),top=4095)
  p4=p
   minpt=min(p4)>0
   maxpt=max(p4)
   display_in=uint(65535d0*(double(p4)-double(minpt))/(double(maxpt)-double(minpt)))
   storfile2="~/triplet4."+sstr(n)+".png"
   write_png,storfile2,display_in
   
  p5=p^.5
   minpt=min(p5)>0
   maxpt=max(p5)
   display_in=uint(65535d0*(double(p5)-double(minpt))/(double(maxpt)-double(minpt)))
   storfile3="~/triplet5."+sstr(n)+".png"
   write_png,storfile3,display_in
   
   proot=(p)^0.5
    pout3=bytscl(proot,min=min(proot),max=2*(.8+n*.1)*mean(proot),top=255)
   outfile="~/triplet3."+sstr(n)+".jpg"
   write_jpeg,outfile,pout,quality=100
   
   endfor
  
  ;stop
 
ilast=file_no 
goto,again
exit:
cd,cdir
end