1TRANSVERSEMERCATORPROJ(1) GeographicLib Utilities TRANSVERSEMERCATORPROJ(1)
2
6 TransverseMercatorProj -- perform transverse Mercator projection
7
9 TransverseMercatorProj [ -s | -t ] [ -l lon0 ] [ -k k0 ] [ -r ] [ -e a
10 f ] [ -w ] [ -p prec ] [ --comment-delimiter commentdelim ] [ --version
11 | -h | --help ] [ --input-file infile | --input-string instring ] [
12 --line-separator linesep ] [ --output-file outfile ]
13
15 Perform the transverse Mercator projections. Convert geodetic
16 coordinates to transverse Mercator coordinates. The central meridian
17 is given by lon0. The longitude of origin is the equator. The scale
18 on the central meridian is k0. By default an implementation of the
19 exact transverse Mercator projection is used.
20 Geodetic coordinates are provided on standard input as a set of lines
22 containing (blank separated) latitude and longitude (decimal degrees or
23 degrees, minutes, seconds); for detils on the allowed formats for
24 latitude and longitude, see the "GEOGRAPHIC COORDINATES" section of
25 GeoConvert(1). For each set of geodetic coordinates, the corresponding
26 projected easting, x, and northing, y, (meters) are printed on standard
27 output together with the meridian convergence gamma (degrees) and scale
28 k. The meridian convergence is the bearing of grid north (the y axis)
29 measured clockwise from true north.
30
32 -s use the sixth-order Krueger series approximation to the transverse
33 Mercator projection instead of the exact projection.
34 -t use the exact algorithm with the "EXTENDED DOMAIN".
36 -l lon0
38 specify the longitude of origin lon0 (degrees, default 0).
39 -k k0
41 specify the scale k0 on the central meridian (default 0.9996).
42 -r perform the reverse projection. x and y are given on standard
44 input and each line of standard output gives latitude, longitude,
45 gamma, and k.
46 -e a f
48 specify the ellipsoid via the equatorial radius, a and the
49 flattening, f. Setting f = 0 results in a sphere. Specify f < 0
50 for a prolate ellipsoid. A simple fraction, e.g., 1/297, is
51 allowed for f. By default, the WGS84 ellipsoid is used, a =
52 6378137 m, f = 1/298.257223563. If the exact algorithm is used, f
53 must be positive.
54 -w on input and output, longitude precedes latitude (except that on
56 input this can be overridden by a hemisphere designator, N, S, E,
57 W).
58 -p prec
60 set the output precision to prec (default 6). prec is the number
61 of digits after the decimal point for lengths (in meters). For
62 latitudes and longitudes (in degrees), the number of digits after
63 the decimal point is prec + 5. For the convergence (in degrees)
64 and scale, the number of digits after the decimal point is prec +
65 6.
66 --comment-delimiter commentdelim
68 set the comment delimiter to commentdelim (e.g., "#" or "//"). If
69 set, the input lines will be scanned for this delimiter and, if
70 found, the delimiter and the rest of the line will be removed prior
71 to processing and subsequently appended to the output line
72 (separated by a space).
73 --version
75 print version and exit.
76 -h print usage and exit.
78 --help
80 print full documentation and exit.
81 --input-file infile
83 read input from the file infile instead of from standard input; a
84 file name of "-" stands for standard input.
85 --input-string instring
87 read input from the string instring instead of from standard input.
88 All occurrences of the line separator character (default is a
89 semicolon) in instring are converted to newlines before the reading
90 begins.
91 --line-separator linesep
93 set the line separator character to linesep. By default this is a
94 semicolon.
95 --output-file outfile
97 write output to the file outfile instead of to standard output; a
98 file name of "-" stands for standard output.
99
101 The exact transverse Mercator projection has a branch point on the
102 equator at longitudes (relative to lon0) of +/- (1 - e) 90, where e is
103 the eccentricity of the ellipsoid. The standard convention for
104 handling this branch point is to map positive (negative) latitudes into
105 positive (negative) northings y; i.e., a branch cut is placed on the
106 equator. With the extended domain, the northern sheet of the
107 projection is extended into the south hemisphere by pushing the branch
108 cut south from the branch points. See the reference below for details.
109
111 echo 0 90 | TransverseMercatorProj
112 => 25953592.84 9997964.94 90 18.40
113 echo 260e5 100e5 | TransverseMercatorProj -r
114 => -0.02 90.00 90.01 18.48
115
117 An illegal line of input will print an error message to standard output
118 beginning with "ERROR:" and causes TransverseMercatorProj to return an
119 exit code of 1. However, an error does not cause
120 TransverseMercatorProj to terminate; following lines will be converted.
121
123 TransverseMercatorProj was written by Charles Karney.
124
126 The algorithms for the transverse Mercator projection are described in
127 C. F. F. Karney, Transverse Mercator with an accuracy of a few
128 nanometers, J. Geodesy 85(8), 475-485 (Aug. 2011); DOI
129 <https://doi.org/10.1007/s00190-011-0445-3>; preprint
130 <https://arxiv.org/abs/1002.1417>. The explanation of the extended
131 domain of the projection with the -t option is given in Section 5 of
132 this paper.
133
135 TransverseMercatorProj was added to GeographicLib,
136 <https://geographiclib.sourceforge.io>, in 2009-01. Prior to version
137 1.9 it was called TransverseMercatorTest (and its interface was
138 slightly different).
139GeographicLib 1.50.1 2019-12-12 TRANSVERSEMERCATORPROJ(1)