pro triplets,Log

db='/DB'

; PC Version
; Last Modified 2005.01.12 csee
; For individual images...to do the whole Log use tripletfiles.pro
;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

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

;____________________________________________________________________________
;check to see if Bushroe's processor has been run on this log.

;file_ck=findfile(log+'\DB2')
;if file_ck(0) ne '' then db='\DB2'
;;if file_ck(0) ne '' then begin      ;i.e. /DB2 does exist
 ;;ans='n'
 ;;read,'Do you want to use the images in /DB2 (Bushroe Processed)? (y or n): ',ans
 ;;if ans eq 'y' then begin & print,'using /DB2' & db='/DB2' & endif else print,'using /DB'
;;endif

cd,Log+db+'/Image'

;____________________________________________________________________________
;List Files...

files=t_getfil('*')
;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
first_file=fix(strmid(files(0),4,3))       ;the number of the first file
last_file=fix(strmid(files(nfiles-1),4,3))  ;the number of the last file

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

;____________________________________________________________________________
;Initalize...

auto=2   ;for auto=0:manual_stretch  auto=1:tvscl,min-to-max  auto=2:tvscl,min-to-2.5xaverage
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.5 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(last_file)+'? (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,'Old Min = '+sstr(min)+'; New Min? ',min
  print,'Min set to: ',sstr(min)
  file_no=ilast
 endif

 if file_no eq -3 then begin       ;set maximum of stretch
  auto=0
  read,'Old Max = '+sstr(max)+'; New Max? ',max
  print,'Max set to: ',sstr(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.5 x average
  auto=2
  print,'Set to Autoscale, min to 2.5xAve.'
  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,'Images: '+Images+'  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 last_file 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...

   file_name=t_getfil('*'+sstr(file_no)+'I'+'*')
   file_name=file_name(0)
   if file_name eq '' then begin & print,'Invalid File No.' & goto,menu & endif
   file=where(files eq file_name)
   file=file(0)
   filetime=strmid(files(file),13,13)

   triplet_files=t_getfil('*'+filetime+'*') & siz=size(triplet_files) & nfiles2=siz(1)
   imno=intarr(3)      ;actual image numbers from file names
   for i=0,nfiles2-1 do imno(i)=fix(strmid(triplet_files(i),4,3))
   if nfiles2 ne 3 then print,'Not a complete triplet',imno

   array_sort,triplet_files     ;sorts files in assending order

   p0=intarr(176,256)      ;MRI frame 176x256 or 254
   p1=intarr(128,256)      ;SLI frame 128x256 or 254
   p2=intarr(160,256)      ;HRI frame 160x256 or 254

   for i=0,nfiles2-1 do begin
    d_read,triplet_files(i),h,p     ;read in each image in order
    type=d_value(h,96)       ;21=MRI,22=SLI, 23=HRI
     case type of
      21: p0=p
      22: p1=p
      23: p2=p
      else:print,'type of ',triplet_files(i),' is unknown'
     endcase
    endfor

   ;;d_read,triplet_files(0),h,p0   ;MRI 176x256 or 254
   ;;d_read,triplet_files(1),h,p1   ;SLI 128x256 or 254
   ;;d_read,triplet_files(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

;____________________________________________________________________________
;condition data...

;remove hot pixels & cosmic rays (less than 15 pixels greater that 6 sigma)
index=where(p0 gt mean(p0)+6.*stdev(p0))
if index(0) ne -1 then begin
 siz=size(index)
 if siz(1) le 15 then p0(index)=Max(p0(where(p0 le mean(p0)+6.*stdev(p0))))
endif

index=where(p1 gt mean(p1)+6.*stdev(p1))
if index(0) ne -1 then begin
 siz=size(index)
 if siz(1) le 15 then p1(index)=Max(p1(where(p1 le mean(p1)+6.*stdev(p1))))
endif

index=where(p2 gt mean(p2)+6.*stdev(p2))
if index(0) ne -1 then begin
 siz=size(index)
 if siz(1) le 15 then p2(index)=Max(p2(where(p2 le mean(p2)+6.*stdev(p2))))
endif


;scale each from 0 min to 4095 max (skip if all zeros = missing image)...
if mean(p0) ne 0 then p0=(P0-min(p0))*4095./(max(p0)-min(p0))
if mean(p1) ne 0 then p1=(p1-min(p1))*4095./(max(p1)-min(p1))
if mean(p2) ne 0 then p2=(p2-min(p2))*4095./(max(p2)-min(p2))

;scale so that their means are equal...
;equilize images to the mean of the brightest...
q=max([mean(p0),mean(p1),mean(p2)])
if mean(p0) ne 0 then  p0=p0*q/mean(p0)
if mean(p1) ne 0 then  p1=p1*q/mean(p1)
if mean(p2) ne 0 then  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

;;print,'p0 ',min(p0),mean(p0),max(p0),stdev(p0),(max(p0)-mean(p0))/stdev(p0)
;;print,'p1 ',min(p1),mean(p1),max(p1),stdev(p1),(max(p1)-mean(p1))/stdev(p1)
;;print,'p2 ',min(p2),mean(p2),max(p2),stdev(p2),(max(p2)-mean(p2))/stdev(p2)
;;; index=where(p0 ge mean(p0)+6*stdev(p0))

;____________________________________________________________________________
;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(where(p eq 0))=min([p0,p1,p2])    ;set all 'empty' pixels to the min of the data

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

 if auto eq 2 then begin & min=min(p) & max=2.5*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

 Images=sstr(imno(0))+', '+sstr(imno(1))+' & '+sstr(imno(2))
 print,string(10b),'Images: ',Images,'  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'


;____________________________________________________________________________
;Done!

continue:

ilast=file_no
goto,again
exit:
cd,cdir
end