[AIT logo]

Institut für Astronomie und Astrophysik

Abteilung Astronomie

Waldhäuser Str. 64, D-72076 Tübingen, Germany
[Uni logo]


CORREL_OPTIMIZE Source code in correl_optimize.pro

CORREL_OPTIMIZE

Name
       CORREL_OPTIMIZE
Purpose
       Find the optimal (x,y) pixel offset of image_B relative to image_A
Explanation
       Optimal offset is computed by means of maximizing the correlation
       function of the two images.
Calling Sequence
       CORREL_OPTIMIZE, image_A, image_B, xoffset_optimum, yoffset_optimum
               [ XOFF_INIT=, YOFF_INIT=, MAGNIFICATION=, /PRINT, /NUMPIX,
                 /MONITOR, PLATEAU_THRESH=  ]
Input Parameters
       image_A, image_B = the two images of interest.
Optional Input Keywords
       XOFF_INIT = initial X pixel offset of image_B relative to image_A,
       YOFF_INIT = Y pixel offset, (default offsets are 0 and 0).
       MAGNIFICATION = option to determine offsets up to fractional pixels,
                       (example: MAG=2 means 1/2 pixel accuracy, default=1).
       /NUMPIX: sqrt( sqrt( # pixels )) used as correlation weighting factor.
       /MONITOR causes the progress of computation to be briefly printed.
       /PRINT causes the results of analysis to be printed.
       PLATEAU_THRESH = threshold used for detecting plateaus in
               the cross-correlation matrix near maximum, (default=0.01),
               used only if MAGNIFICATION > 1.    Decrease this value for
               high signal-to-noise data
Output Parameters
       xoffset_optimum = optimal X pixel offset of image_B relative to image_A.
       yoffset_optimum = optimal Y pixel offset.
Procedures Used
       function  correl_images( image_A, image_B )
       pro  corrmat_analyze
Procedure
       The combination of function correl_images( image_A, image_B ) and
       corrmat_analyze of the result is used to obtain the (x,y) offset
       yielding maximal correlation. The combination is first executed at
       large REDUCTION factors to speed up computation, then zooming in
       recursively on the optimal (x,y) offset by factors of 2.
       Finally, the MAGNIFICATION option (if specified)
       is executed to determine the (x,y) offset up to fractional pixels.
Revision History
       Written, July,1991, Frank Varosi, STX @ NASA/GSFC
       Added PLATEAU_THRESH keyword  June 1997,  Wayne Landsman  STX
       Converted to IDL V5.0   W. Landsman   September 1997

Last modified by pro2html on 2001 April 26 at 03:13 UTC

[Home Page] [Software, Documentation] [IDL Documentation] [Quick Reference] [Feedback]

Jörn Wilms (wilms@astro.uni-tuebingen.de)
Updated automatically