Help for GTPWARP

PURPOSE
     GTPWARP acts like GTWARP or GTWARPXD in that it converts an input
     image to an output image, matching the mapping of a reference
     image.  Consider the following four cases:
     
        input                 reference               use
       
     1. GeoTIFF label exact   GeoTIFF label exact     GTWARP
     2. GeoTIFF label approx  GeoTIFF label exact     GTPWARP or GTWARPXD
     3. GeoTIFF label exact   GeoTIFF label approx    GTPWARP or GTWARPXD
     4. GeoTIFF label approx  GeoTIFF label approx    GTPWARP or GTWARPXD
     
     For cases 2 through 4, GTPWARP is highly automatic, but will
     depend upon PICMTCH5 being able to do a good enough correlation
     to give the desired result.  GTWARPXD is more general purpose,
     but the user has to come up with a tiepoint data set.
     
     GTPWARP reads the GeoTIFF label of a "master" or "reference" image,
     then reads the GeoTIFF label of the input image.  It calculates
     a set of GEOMV parameters from these.  Then it applies PICMTCH5 to
     the input and reference to get a mesh of matching points that
     specify the deviation of the mappings due to inaccuracy in the
     mappings.  Both the deviations and the mapping are input to routine
     GTWARPXD to warp the input to the output giving an accurate match
     to the "reference" image.
     
     1.  The keyword 'coverinp causes the procedure to map all of the
     input image to an output image that will just contain it.   The
     keyword gorewid allows the user to specify a zero-filled margin
     to make the output a little larger (or negative trims smaller).
	
     2.  The keyword 'coverref causes the procedure to map the input
     image to an output image that exactly matches the "reference" image.
     Some parts of the input can be lost, or huge areas of zero-fill
     could be added to make the images match.  The gorewid keyword has
     no effect in this case.

    Note that if the resolution of the reference image is much lower than
    the input image, you should specify ref_lowres=y. This will
    swap the order of images in picmtch5 which leads to better matches.
	   
CALL
     gtpwarp INPUT OUTPUT REF '(QUALIFIERS) PARAMS
  WHERE:
     INPUT          is the input data set (must have a GeoTIFF label).
     OUTPUT         is the output data set (will have a GeoTIFF label).
     REF            is the reference data set (must have a GeoTIFF label).
     QUALIFIERS     consist of any of the following keywords:
          COVERINP         smallest output that covers the input geographic area
          COVERREF         output exactly matches the reference image.
     PARAMS         interp,nah,nav,gorewid, fftsize.

  
OPERATION

The program calls PICMTCH5 to create a deviation data set
PERFORMANCE

see PICMTCH5 and GTWARPXD.

Restrictions
------------

The input and reference images must have GeoTIFF labels.


Original Programmer: A. L. Zobrist, 10 Dec, 2000
Last Update:         A. L. Zobrist, 4 Sep, 2007     Additional parameters for Tom
Current Cognizant Programmer: A. L. Zobrist

HISTORY
  2021-05-24 WLB - Migrated from AFIDS


PARAMETERS:


INP

Input file name with GeoTIFF label

OUT

Output file name

REF

Reference file name with GeoTIFF label

REF_RES

Reference file name with GeoTIFF label to resample to. This can be ref at a different resolution (e.g., using gtsize), useful when we want to match using one resolution but resample the input file to a different resolution.

TYPREF

'COVERINP - output minimally covers the input data 'COVERREF - output matches the ref image exactly

NAH

Number of grid cells horiz.

NAV

Number of grid cells vert.

INTERP

interpolation options Valid: NOIN,ZNOIN,BILIN

GOREWID

Added width of gores for 'coverinp case only

FFTSIZE

FFT window size

FFTGRID

how many fft points, nah x nav

MAGNIFY

enlarge footprint of fft window by this factor

POLYFIT

if given, overrides triangu- lation surface fit with poly- nomial fit (see tieconv)

ACCPLOT

if given, the accuracy plot image will have this name, otherwise no plot output for no plot calculation

ACCOFF

set to "y" to turn off accuracy

GRIDOUT

name of file to save warp grid

MAGMIN

minimum magnifier as corr- elation proceeds

FFTHALF

0 = no halving of fftsize near edge 1 = single halving of fftsize near edge 2 = double halving of fftsize near edge

AUTOFIT

Sets the FITMIN threshold auto- matically using the first n correlations (see help 2)

REDO

At the end of processing, this number of points will be redone

ITIE

Control points in first image, only if map labels not used

REFTIE

Control points in second image, only if map labels not used

SEED

seed for randomizing picmtch5 grid

WMASK

If n skip using the watermask in accck

REF_LOWRES

If y, then switch the order of matching in picmtch5 to give better results

THR_RES

Prevents offset correlations (further than this value from the predicted location) from entering into the model fit of image 1 to image 2

See Examples:


Cognizant Programmer: