Skip site navigation (1)Skip section navigation (2)

FreeBSD Manual Pages

  
 
  

home | help
TRANSVERSEMERCATORPROJ(1)    GeographicLib Utilities   TRANSVERSEMERCATORPROJ(1)

NAME
     TransverseMercatorProj -- perform transverse Mercator projection

SYNOPSIS
     TransverseMercatorProj  [ -s | -t ] [ -l lon0 ] [ -k k0 ] [ -r ] [ -e a f ]
     [ -w ] [ -p prec ] [ --comment-delimiter commentdelim ] [ --version | -h  |
     --help ] [ --input-file infile | --input-string instring ] [ --line-separa-
     tor linesep ] [ --output-file outfile ]

DESCRIPTION
     Perform  the transverse Mercator projections.  Convert geodetic coordinates
     to transverse Mercator coordinates.  The central meridian is given by lon0.
     The longitude of origin is the equator.  The scale on the central	meridian
     is  k0.  By default an implementation of the exact transverse Mercator pro-
     jection is used.

     Geodetic coordinates are provided on standard input as a set of lines  con-
     taining  (blank  separated)  latitude and longitude (decimal degrees or de-
     grees, minutes, seconds); for detils on the allowed  formats  for	latitude
     and  longitude,  see the "GEOGRAPHIC COORDINATES" section of GeoConvert(1).
     For each set of geodetic coordinates, the corresponding projected	easting,
     x,  and  northing, y, (meters) are printed on standard output together with
     the meridian convergence gamma (degrees) and scale k.  The meridian conver-
     gence is the bearing of grid north (the y	axis)  measured  clockwise  from
     true north.

OPTIONS
     -s  use the sixth-order Krueger series approximation to the transverse Mer-
	 cator projection instead of the exact projection.

     -t  use  the  exact  algorithm  with the "EXTENDED DOMAIN"; this is the de-
	 fault.

     -l lon0
	 specify the longitude of origin lon0 (degrees, default 0).

     -k k0
	 specify the scale k0 on the central meridian (default 0.9996).

     -r  perform the reverse projection.  x and y are given  on  standard  input
	 and  each line of standard output gives latitude, longitude, gamma, and
	 k.

     -e a f
	 specify the ellipsoid via the equatorial radius, a and the  flattening,
	 f.  Setting f = 0 results in a sphere.  Specify f < 0 for a prolate el-
	 lipsoid.   A  simple  fraction,  e.g., 1/297, is allowed for f.  By de-
	 fault, the WGS84 ellipsoid is used, a = 6378137 m, f = 1/298.257223563.
	 If the exact algorithm is used, f must be positive.

     -w  on input and output, longitude precedes latitude (except that on  input
	 this can be overridden by a hemisphere designator, N, S, E, W).

     -p prec
	 set  the  output  precision to prec (default 6).  prec is the number of
	 digits after the decimal point for lengths (in meters).  For  latitudes
	 and  longitudes  (in  degrees),  the number of digits after the decimal
	 point is prec + 5.  For the convergence (in  degrees)	and  scale,  the
	 number of digits after the decimal point is prec + 6.

     --comment-delimiter commentdelim
	 set the comment delimiter to commentdelim (e.g., "#" or "//").  If set,
	 the  input  lines will be scanned for this delimiter and, if found, the
	 delimiter and the rest of the line will be removed prior to  processing
	 and subsequently appended to the output line (separated by a space).

     --version
	 print version and exit.

     -h  print usage and exit.

     --help
	 print full documentation and exit.

     --input-file infile
	 read  input from the file infile instead of from standard input; a file
	 name of "-" stands for standard input.

     --input-string instring
	 read input from the string instring instead  of  from	standard  input.
	 All  occurrences  of  the  line separator character (default is a semi-
	 colon) in instring are converted to newlines before the reading begins.

     --line-separator linesep
	 set the line separator character to linesep.	By  default  this  is  a
	 semicolon.

     --output-file outfile
	 write	output to the file outfile instead of to standard output; a file
	 name of "-" stands for standard output.

EXTENDED DOMAIN
     The exact transverse Mercator projection has a branch point on the  equator
     at  longitudes (relative to lon0) of +/- (1 - e) 90 = 82.636..., where e is
     the eccentricity of the ellipsoid.  The standard  convention  for	handling
     this  branch  point  is  to map positive (negative) latitudes into positive
     (negative) northings y; i.e., a branch cut is placed on the equator.   With
     the  extended domain, the northern sheet of the projection is extended into
     the south hemisphere by pushing  the  branch  cut	south  from  the  branch
     points.  See the reference below for details.

EXAMPLES
	echo 0 90 | TransverseMercatorProj
	=> 25953592.84 9997964.94 90 18.40
	echo 260e5 100e5 | TransverseMercatorProj -r
	=> -0.02 90.00 90.01 18.48

ERRORS
     An illegal line of input will print an error message to standard output be-
     ginning  with  "ERROR:" and causes TransverseMercatorProj to return an exit
     code of 1.  However, an error does not cause TransverseMercatorProj to ter-
     minate; following lines will be converted.

AUTHOR
     TransverseMercatorProj was written by Charles Karney.

SEE ALSO
     The algorithms for the transverse Mercator projection are described  in  C.
     F.  F. Karney, Transverse Mercator with an accuracy of a few nanometers, J.
     Geodesy	    85(8),	  475-485	 (Aug.	      2011);	     DOI
     <https://doi.org/10.1007/s00190-011-0445-3>;			preprint
     <https://arxiv.org/abs/1002.1417>.  The explanation of the extended  domain
     of the projection with the -t option is given in Section 5 of this paper.

HISTORY
     TransverseMercatorProj  was  added  to  GeographicLib,  <https://geographi-
     clib.sourceforge.io>, in 2009-01.	Prior  to  version  1.9  it  was  called
     TransverseMercatorTest (and its interface was slightly different).

GeographicLib 2.2		   2022-12-13	       TRANSVERSEMERCATORPROJ(1)

Want to link to this manual page? Use this URL:
<https://man.freebsd.org/cgi/man.cgi?query=TransverseMercatorProj&sektion=1&manpath=FreeBSD+Ports+15.1.quarterly>

home | help