_______________

The calling sequence of DESTRE_F program is

IDL> destre_f,sel,wid,nref,reb

where SEL is the array containing file numbers (and it controls which
frames will be processed), WID is the FWHM of the Gaussian window
which determines the "cell size" in destretching, NREF is a number
of frames which are averaged to get the reference image and REB
is a down-rebinning factor to diminish the frame size and speed-up
the calculations. The optimum combination of the parameters is:

         Run no.   wid       nref    reb
         -------------------------------
           1       64 (8")     7      3		; large scale
           2        8 (1")     5      2		; small scale
           3        8          5      1		; residual jittering

At the beginning, the program asks for the first part of the input
filename (default = spotgs) and of the output filename (default =
= spotde). The filenames are constructed in the form spotgs.x
(spotde.x), where x is the file number contained in the 0-column
of sel. You can change the default settings in the "read" and "if"
commands at the beginning of the destre_f.pro source file.
Do not change the unsharp masking parameters.
_______________

In the working directory we need to have the programs:

destre_f.pro
deformo.pro		; called by destre_f
nint.pro			"
smoothe.pro			"
distorti.pro			"
sconvol.pro			"
kconvol.pro			"
odd.pro				"
mean.pro			"

vyber2.pro		; to convert 'log' to 'sel'
cinema3.pro		; to create movie of the spot
cinema4.pro		; to create movie of the pore
animate.pro		; to look at the movie
sta.pro			; to look at the statistics

(all the support programs are present in /home/sobotka)

and the data files:

log			; logfile (SAVE file)
spotgs.*		; spot (SAVE2 files)
poregs.*		; pore (     "     )


How to run the program:
-----------------------

1.  IDL> restore,'log'

2.  IDL> vyber2,log,sel,0.,0.

  Conversion of 'log' to 'sel' array. Selects all files.
  If you want to process only a part of the series, you can use
    IDL> vyber2,log(*,from:to),sel,0.,0.
  where from, to are indices in the log array. 

3.  IDL> destre_f,sel,64,7,3
    "INPUT  filenames (default = spotgs):"  press Enter (default)
    "OUTPUT filenames (default = spotde):"  press Enter (default)
    
  1st destretching run, runtime for 760 frames at Sparc10 = 1h 10m.

4.  IDL> destre_f,sel,8,5,2
    "INPUT  filenames (default = spotgs):"  spotde
    "OUTPUT filenames (default = spotde):"  press Enter (default)
    
  2nd run, runtime for 760 frames at Sparc10 = 45m.
	   
5.  IDL> destre_f,sel,8,5,1
    "INPUT  filenames (default = spotgs):"  spotde
    "OUTPUT filenames (default = spotde):"  press Enter (default)
    
  3rd run, runtime for 760 frames at Sparc10 = 1h 30m.
      
  This way, the spotde.x are always overwritten by the next run
  to save the disk space. You can also assign them different names
  to have all the intermediate results, but it is not so important.
  
6.  IDL> cinema3,sel(*,from:to),mov

  Create a 3-D movie array MOV of a part of the series determined
  by FROM, TO. We know that 200 frames are acceptable for Sparc10.
  Do not forget to CHANGE PATH AND FILENAME BEFORE YOU RUN CINEMA3
  in the cinema3.pro source file!
  
7.  IDL> animate,bytscl(mov)

  Enjoy the movie and check the result...
_________________

I did not check the performance in case of the pore destretching but
it can be assumed that the parameters are the same. The only changes
necessary are in the filenames (default settings in destre_f.pro and
path-and-filename change in cinema4.pro,  which is used for the pore
instead of cinema3.pro).