Help for GTPROJ

PURPOSE

   "gtproj" converts map coordinates from one map projection to
another using the EPOSC map conversion routines.  The input can
either by (x,y) or (x,y,z) columns in an IBIS file.  If two
columns are specified (the x,y case), then an all-zero z column
is generated for input to the EPOSC routines.  Similarly, if
no z column is provided for output, then the EPOSC z output is
ignored.  The name of the IBIS file is the first element of the
INP parameter.

   The map coordinate specifications for both input and output
must be in the headers of two VICAR files which are the second
and third elements of the INP parameter.  These files must have
GeoTIFF state labels with sufficient parameters to specify a
mapping.

If zero or one datum is given, then the conversion will assume
same datum for input and output.  If both input and output have
a datum, then datum conversion will take place as well.
 
The columns in the IBIS file for coordinate I/O are given
by the parameters INCOL and OUTCOL.

The EP_PAR parameter is added for special use of the EPOSC parameters
in code testing situations.  To obtain a listing of these, 
call Gerald I. Evenden at (508)563-6766 or write to same at:



        P.O. Box 1027

        N. Falmouth, MA 02556-1027

        USA


or e-mail  gie@kai.er.usgs.gov

The 'TEC keyword is provided to use an alternative mapping package
from US ARMY TEC to calculate the mapping.  This can be used to test
or confirm the EPOSC mapping routine, or to replace it if it doesn't 
work.  The TEC package had fourteen map projection routines that 
matched up to the eposc package so those have been activated:

Cassini
Equidistant Cylindrical (Plate Caree)
van der Grinten (I)
Lambert Conformal Conic
Mercator
Miller Cylindrical
Orthographic
Polyconic
Sinusoidal
Transverse Mercator
Universal Polar Stereographic
Albers Equal Area
UTM

Spheroids that have been activated have the codes:

Airy                           AA 6377563.396 6356256.9090 299.324964600

Modified Airy                  AM 6377340.189 6356034.4480 299.324964600

Australian National            AN 6378160.000 6356774.7190 298.250000000

Bessel 1841                    BR 6377397.155 6356078.9630 299.152812800

Bessel 1841(Namibia)           BN 6377483.865 6356165.3830 299.152812800

Clarke 1866                    CC 6378206.400 6356583.8000 294.978698200

Clarke 1880                    CD 6378249.145 6356514.8700 293.465000000

Everest                        EA 6377276.345 6356075.4130 300.801700000

Everest (E. Malasia, Brunei)   EB 6377298.556 6356097.5500 300.801700000

Everest 1956 (India)           EC 6377301.243 6356100.2280 300.801700000

Everest 1969 (West Malasia)    ED 6377295.664 6356094.6680 300.801700000

Everest 1948(W.Mals. & Sing.)  EE 6377304.063 6356103.0390 300.801700000

Mod. Fischer 1960(South Asia)  FA 6378155.000 6356773.3200 298.300000000

GRS 80                         RF 6378137.000 6356752.3141 298.257222101

Helmert 1906                   HE 6378200.000 6356818.1700 298.300000000

Hough                          HO 6378270.000 6356794.3430 297.000000000

Indonesian 1974                ID 6378160.000 6356774.5040 298.247000000

International                  IN 6378388.000 6356911.9460 297.000000000

Krassovsky                     KA 6378245.000 6356863.0190 298.300000000

South American 1969            SA 6378160.000 6356774.7190 298.250000000

WGS 72                         WD 6378135.000 6356750.5200 298.260000000

WGS 84                         WE 6378137.000 6356752.3142 298.257223563

Other spheroids can be input by using their a,b values

EXAMPLE

gtproj (myibisfile,img_lambert,img_albers) incol=(2,3,4) outcol=(5,6,7)


RESTRICTIONS

None.  The IBIS column length is dynamically allocated. 

WRITTEN BY:            Albert L. Zobrist	Feb 2000

COGNIZANT PROGRAMMER:  Barbara McGuffie

Revisions:
2008-01-11 WLB Switched to USES_ANSI_C AND LIB_CARTO; misc cleanup  
2016-01-08 WLB Migrated to MIPL; replaced printfs with zifmessage calls
2017-08-08 WLB Removed DEBIG flag.

PARAMETERS: