;  .com destretch.pro
;   .r  destretch.pro
pro cw4_tio_destretch

  If ( !version.os eq 'Win32' ) Then st = '\'
  If ( !version.os eq 'linux' ) Then st = '/'

  Print, 'Where are the  files after alignment?'

  dir_read = DIALOG_PICKFILE( /DIRECTORY )

   file_mkdir,dir_read+'destretch'
   dir_write = dir_read+'destretch'+st
   list = file_search(dir_read + '*aln.fts',count=nn)
   pos = STRPOS(list[0], 'aln.fts')
   spos = Strpos( list[0], st, /reverse_search )
   sn = pos - spos
   listname = STRMID(list, spos+1, sn-1)

   ;head = strarr(71,nn)
   img0 = readfits(list[0],head0)
   ;head[*,0]=head0
   sz = size(img0)
   img = fltarr(sz[1],sz[2],nn)
   img[*,*,0]=img0
   print,sz
  for i = 1,nn-1 do begin
     img[*,*,i] = readfits(list[i],head0)
     ;head[*,i] = head0
  endfor

   ;movie,bytscl(img,2500,8500),order=0
   ;xtvscl,img[*,*,fix(nn/2)]
   ;box_cursor,x,y,dx,dy
   ;img=img[round(x):round(x)+dx-1,round(y):round(y)+dy-1,*]

   img = do_destretch(img)
    ;movie,bytscl(img,2500,8500),order=0

;for the img not equal bright
  ; print,'Enter y if you do not want to adjust the bright'
   ;  cri = get_kbrd(60)

   ; case strupcase(cri) of
   ;  'Y' : begin
   ;          i = fix(nn/2)
    ;         xtvscl,img[*,*,i]
     ;        box_cursor,x,y,dx,dy
      ;      imgg = img[round(x):round(x)+dx-1,round(y):round(y)+dy-1,i]
       ;     refb = mean(imgg)
        ;      for i = 0,sz(3)-1 do begin
         ;      imgg = img[*,*,i]*refb/mean(img[round(x):round(x)+dx-1,round(y):round(y)+dy-1,i])
          ;     img[*,*,i] = imgg
           ;   endfor
           ;end
     ;else :
    ;endcase

  for i = 0,nn-1 do WRITEFITS, dir_write+listname[i]+'dstr.fts', img[*,*,i]
  ;, head[*,i]




end