TRANSVERSE MERCATOR AND LAMBERT CONFORMAL …

1 TRANSVERSE MERCATOR AND LAMBERT CONFORMAL CONIC. MAP PROJECTION FUNCTIONS. Alan Vonderohe (February 2020). Contents 1. Background. 2. General Notation and Definitions. 3. Ellipsoid Constants. 4. TRANSVERSE MERCATOR Projections (After Stem (1989)). a. Projection-Specific Notation. b. Projection Constants. c. Direct Transformation ( , to N,E). d. Inverse Transformation (N,E to , ). e. Scale and Convergence. 5. LAMBERT CONFORMAL Conic Projections. a. Secant (Two Standard Parallels) (After Stem (1989)). i. Projection-Specific Notation. ii. Projection Constants. iii. Direct Transformation ( , to N,E). iv. Inverse Transformation (N,E to , ). b. Non-Intersecting (Central Parallel and Its Scale Factor) (After Bomford (1985)). i. Projection-Specific Notation.

2 Ii. Projection Constants. iii. Direct Transformation ( , to N,E). iv. Inverse Transformation (N,E to , ). v. Scale and Convergence. 6. Linear Distortion at , . 7. Geocentric Coordinates. a. Direct Transformation ( , ,h to X,Y,Z). b. Inverse Transformation (X,Y,Z to , ,h). 8. List of References. 1. 1. Background. A map projection is a mathematical surface with functional relationships between points having geodetic coordinates (latitude and longitude) on a reference ellipsoid and corresponding points with two-dimensional rectangular coordinates (northing and easting) on the map projection surface. A map projection is described by mathematical transformations between the two types of coordinates. A direct transformation computes northing and easting (N,E) from latitude and longitude ( , ): N = f1( , , ellipsoid parameters, map projection parameters).

3 E = f2( , , ellipsoid parameters, map projection parameters). An inverse transformation computes , from N,E: = g1(N,E, ellipsoid parameters, map projection parameters). = g2(N,E, ellipsoid parameters, map projection parameters). In the equations above, ellipsoid parameters are two descriptors that define the size and shape of the reference ellipsoid. In this document, those two parameters are the semi-major axis (a) and the semi-minor axis (b). Map projection parameters are descriptors that define the size and shape of the map projection surface and its location and orientation with respect to the reference ellipsoid. Map projection parameters also define the location of the rectangular coordinate origin and its false northing and false easting values.

4 On any map projection, each point has a scale factor and a convergence. On CONFORMAL map projections, such as TRANSVERSE MERCATOR and LAMBERT CONFORMAL conic, the scale factor is the same in all directions at any given point but is variable from point to point. Convergence, sometimes referred to as the mapping angle , also varies from point to point and is the angle between geodetic north and grid north. It is defined as a geodetic azimuth minus the projection of that azimuth on the map projection coordinate grid. Convergence at a point is not the difference between geodetic azimuth and grid azimuth. Such a difference depends not only upon convergence but also upon the arc-to-chord or second term correction. Arc-to-chord corrections, at any point, vary with distance and direction to an arbitrary second point, whose coordinates must be specified.

5 This 2. document presents methods for computing scale and convergence. It does not address methods for computing arc-to-chord corrections. For many applications arc-to-chord corrections are negligible. The scale factor (denoted k) and convergence (denoted ) are found by functions: k = h1( , , a, b, map projection parameters). = h2( , , a, b, map projection parameters). Finally, any given point on a map projection has a linear distortion that is the ratio of a very small distance on Earth's surface to the corresponding very small distance on the map projection surface. Linear distortion accounts for the separation between the two surfaces and expresses the error to be encountered if ground distances are used for grid distances: LD = i( , a, b, k, ellipsoid height).

6 This document presents each of the functions described above for TRANSVERSE MERCATOR projections and two types of LAMBERT CONFORMAL conic projections. 2. General Notation and Definitions. Geodetic latitude, positive north. Geodetic longitude, positive east (0 to 360 ). N Northing coordinate on the projection. E Easting coordinate on the projection. o Longitude of the central meridian and the coordinate origin. Eo False easting of the coordinate origin. o Latitude of the coordinate origin ( TRANSVERSE MERCATOR and non-intersecting LAMBERT CONFORMAL conic). Also, latitude of the central parallel ( LAMBERT CONFORMAL conic). ko Scale factor along the central meridian ( TRANSVERSE MERCATOR ) or central parallel ( LAMBERT CONFORMAL conic). k Scale factor. Convergence.

7 A Semi-major axis of the reference ellipsoid. b Semi-minor axis of the reference ellipsoid. e First eccentricity of the reference ellipsoid. e' Second eccentricity of the reference ellipsoid. Radius of curvature in the prime vertical. 3. Radius of curvature in the meridian. RG Gaussian or geometric mean radius of curvature. Ng Geoid height. h Ellipsoid height. H Orthometric height. 3. Ellipsoid Constants. a 2 b2. e2 =. a2. e2. e '2 =. 1 e2. a b n=. a+b 4. TRANSVERSE MERCATOR Projections (After Stem (1989)). TRANSVERSE MERCATOR projections are based upon right cylinders whose axes lie in the equatorial plane and pass through the center of the reference ellipsoid. The selected right cylinder can be secant or tangent to the reference ellipsoid. It can also not intersect the reference ellipsoid at all.

8 Projection-Specific Notation. NOTE: Projection parameters are o, o, ko, No, Eo. Rectifying latitude. S Meridional distance from the equator, multiplied by ko. o Rectifying latitude at o. So Meridional distance from the equator to o, multiplied by ko. r Radius of the rectifying sphere. No False northing of the coordinate origin. 4. Projection Constants. 9n 2 225n 4 . r = a (1 n)(1 n 2 ) 1 + + . 4 64 . 3n 9n3. u2 = +. 2 16. 15n 15n 4. 2. u4 = . 16 32. 3. 35n u6 = . 48. 315n 4. u8 =. 512. U 0 = 2(u2 2u4 + 3u6 4u8 ). U 2 = 8(u4 4u6 + 10u8 ). U 4 = 32(u6 6u8 ). U 6 = 128u8. 3n 27n3. v2 = . 2 32. 2. 21n 55n 4. v4 = . 16 32. 3. 151n v6 =. 96. 1097n 4. v8 =. 512. V0 = 2(v2 2v4 + 3v6 4v8 ). V2 = 8(v4 4v6 + 10v8 ). V4 = 32(v6 6v8 ). V6 = 128v8. o = o + sin o cos o (U 0 + U 2 cos2 o + U 4 cos4 o + U 6 cos6 o ).

9 So = ko o r Direct Transformation ( , to N,E). 5. t = tan . 2 = e'2 cos2 . ko a R=. 1 e 2 sin 2 . A1 = R. Rt A2 =. 2. 1 t2 + 2. A3 =. 6. 5 t 2 + 2 (9 + 4 2 ). A4 =. 12. 5 18t 2 + t 4 + 2 (14 58t 2 ). A5 =. 120. 61 58t + t 4 + 2 (270 330t 2 ). 2. A6 =. 360. 61 479t + 179t 4 t 6. 2. A7 =. 5040.. N = S So + N o + A2 L2 1 + L2 ( A4 + A6 L2 ) .. E = Eo + A1L 1 + L2 A3 + L2 ( A5 + A7 L2 ) . Inverse Transformation (N,E to , ). N N o + So =. ko r f = + sin cos (V0 + V2 cos 2 + V4 cos 4 + V6 cos6 ). t f = tan f 2f = e' 2 cos 2 f ko a Rf =. 1 e 2 sin 2 f E Eo Q=. Rf t f (1 + 2f ). B2 =. 2. 6. (1 + 2t 2f + 2f ). B3 =. 6. (5 + 3t 2f + 2f (1 9t 2f ) 4 4f ). B4 =. 12. 5 + 28t + 24t 4f + 2f (6 + 8t 2f ). 2. B5 =. f 120. 61 + 90t + 45t 4f + 2f (46 252t 2f 90t 4f ).

10 2. B6 =. f 360. (61 + 662t + 1320t 4f + 720t 6f ). 2. B7 =. f 5040. Scale and Convergence. Convergence ( ) in radians. t = tan . 2 = e' 2 cos2 . 1. C5 = (2 t 2 ). 15.. = C1 L 1 + L2 (C3 + C5 L2 ) . Scale Factor (k). 1. F2 = (1 + 2 ). 2. F4 =. 1. 12.. 5 4t 2 + 2 (9 24t 2 ) .. k = k o 1 + F2 L2 (1 + F4 L2 ) . 7. 5. LAMBERT CONFORMAL Conic Projections. LAMBERT CONFORMAL conic projections are based upon right circular cones whose axes coincide with the minor axis the reference ellipsoid. The selected right circular cone can be secant or tangent to the reference ellipsoid. It can also not intersect the reference ellipsoid at all. Secant (Two Standard Parallels) (After Stem (1989)). Projection-Specific Notation. NOTE: Projection parameters are N, S, b, o, Nb, Eo.