This program produces the Gauss-Kruger (constant meridional scale) Transverse Mercator Projection which is used to construct the U.S. Army's Universal Transverse Mercator (UTM) Grid System. The method is capable of mapping the entire northern hemisphere of the earth (and, by symmetry of the projection, the entire earth) accurately with respect to a single principal meridian, and is therefore mathematically insensitive to proximity either to the pole or the equator, or to the departure of the meridian from the central meridian. This program could be useful to any map-making agency.
The program overcomes the limitations of the "series" method (Thomas, 1952) presently used to compute the UTM Grid, specifically its complicated derivation, non-convergence near the pole, lack of rigorous error analysis, and difficulty of obtaining increased accuracy. The method is based on the principle that the parametric colatitude of a point is the amplitude of the Elliptic Integral of the 2nd Kind, and this (irreducible) integral is the desired projection. Thus, a specification of the colatitude leads, most directly (and with strongest motivation) to a formulation in terms of amplitude. The most difficult problem to be solved was setting up the method so that the Elliptic Integral of the 2nd Kind could be used elsewhere than on the principal meridian.
The point to be mapped is specified in conventional geographic coordinates (geodetic latitude and longitudinal departure from the principal meridian). Using the colatitude (complement of latitude) and the longitude (departure), the initial step is to map the point to the North Polar Stereographic Projection. The closed-form, analytic function that coincides with the North Polar Stereographic Projection of the spheroid along the principal meridian is put into a Newton-Raphson iteration that solves for the tangent of one half the parametric colatitude, generalized to the complex plane. Because the parametric colatitude is the amplitude of the (irreducible) Incomplete Elliptic Integral of the 2nd Kind, the value for the tangent of one half the amplitude of the Elliptic Integral of the 2nd Kind is now known. The elliptic integral may now be computed by any desired method, and the result will be the Gauss-Kruger Transverse Mercator Projection. This result is a consequence of the fact that these steps produce a computation of real distance along the image (in the plane) of the principal meridian, and an analytic continuation of the distance at points that don't lie on the principal meridian. The elliptic-integral method used by this program is one of the "transformations of the elliptic integral" (similar to Landen's Transformation), appearing in standard handbooks of mathematical functions. Only elementary transcendental functions are utilized. The program output is the conventional (as used by the mapping agencies) cartesian coordinates, in meters, of the Transverse Mercator projection. The origin is at the intersection of the principal meridian and the equator.