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: