10.2 Map Projections: Universal Transverse Mercator (UTM) and Philippine Transverse Mercator (PTM) Grid Systems

10.2 Map Projections: Universal Transverse Mercator (UTM) and Philippine Transverse Mercator (PTM) Grid Systems

Mapping a three-dimensional, doubly curved reference ellipsoid onto a two-dimensional flat plane cannot be accomplished without spatial distortion. A map projection is a mathematically defined systematic transformation that maps geographic coordinates (latitude $\phi$, longitude $\lambda$) onto plane Cartesian coordinates (Easting $E$, Northing $N$).

$$\begin{pmatrix} E \\ N \end{pmatrix} = f(\phi, \lambda)$$

1. Distortion Characteristics and Projection Typology

According to Carl Friedrich Gauss's Theorema Egregium, a curved surface cannot be flattened into a plane without stretching, tearing, or compressing the surface geometry. Cartographers characterize projection distortion using Tissot's Indicatrix—an infinitely small circle on the ellipsoid that distorts into an ellipse on the map projection plane.

   Ellipsoid Surface              Conformal Projection            Equal-Area Projection
     (Unit Circle)                     (Enlarged)                      (Flattened)
       ┌───────┐                        ┌─────────┐                     ┌───────────┐
       │   O   │       ─────────►       │    O    │     ─────────►      │    ( )    │
       └───────┘                        └─────────┘                     └───────────┘
     a = b = 1.0                      a = b = 1.05                     a = 1.25, b = 0.8
    Area = 3.1415                    Area = 3.4636                    Area = 3.1415
 (Shape Preserved)                 (Shape Preserved)                (Area Preserved)

Primary Distortion Classes

Map projections are categorized into four major functional classes based on which geometric property they preserve:

1. Conformal (Orthomorphic) Projections: Preserve local angles, differential shapes, and direction relationships. The scale factor at any given point is independent of direction ($k_x = k_y$). Tissot's indicatrix remains a circle, although its size varies across the map. Conformal projections are mandatory for geodetic surveying, navigation, and engineering mapping (e.g., Mercator, Transverse Mercator, Lambert Conformal Conic).
2. Equivalent (Equal-Area) Projections: Preserve area proportions across the entire map sheet. An area of $1\text{ cm}^2$ anywhere on the map represents the exact same ground area. However, local shapes and angles are severely distorted (e.g., Albers Equal-Area, Gall-Peters, Lambert Azimuthal Equal-Area).
3. Equidistant Projections: Maintain true scale along specified lines (such as standard parallels or radial lines from a central point). No projection can maintain true distance scale in all directions across the entire map sheet.
4. Azimuthal (True Direction) Projections: Maintain correct bearings and directions from a chosen central point to all other points on the projection plane.

---

2. Transverse Mercator (TM) Projection Mechanics

The Transverse Mercator projection is a secant, transverse, cylindrical conformal projection. Unlike the standard Mercator projection (where the cylinder is tangent or secant along the Equator), the cylinder axis of the Transverse Mercator lies in the equatorial plane, perpendicular to the Earth's axis of rotation. The cylinder intersects the ellipsoid along a central meridian and two parallel lines of secancy.

Key geometric characteristics of the Transverse Mercator projection include:

• The Central Meridian ($\lambda_0$) projects as a straight vertical line representing Grid North.
• The Equator projects as a straight horizontal line representing the Easting baseline.
• All other meridians and parallels project as complex curved lines intersecting at $90^\circ$ angles.
• The point scale factor ($k$) equals a predefined central meridian scale factor ($k_0 < 1.0$) at the central meridian, increases to $k = 1.0$ along standard secancy lines, and grows greater than $1.0$ as distance from the central meridian increases.

---

3. Universal Transverse Mercator (UTM) Grid System

The Universal Transverse Mercator (UTM) system is an international metric grid framework covering the globe from $80^\circ\text{S}$ latitude to $84^\circ\text{N}$ latitude.

UTM System Parameters

• Zone Width: The globe is divided into 60 longitudinal zones, each spanning $6^\circ$ of longitude.
• Zone Numbering: Zones are numbered 1 to 60 starting at longitude $180^\circ\text{W}$ and proceeding eastward.
• Central Meridian Scale Factor ($k_0$): $0.9996$ (representing a secant cylinder scale reduction of $1/2,500$ along the central meridian).
• False Easting ($E_0$): Assigned a value of $500,000.00\text{ meters}$ at the central meridian to eliminate negative Easting values.
• False Northing ($N_0$): Assigned $0.00\text{ meters}$ for the Northern Hemisphere and $10,000,000.00\text{ meters}$ for the Southern Hemisphere.
• Philippine UTM Coverage: The Philippines spans UTM Zone 51N (Central Meridian $123^\circ\text{E}$) and UTM Zone 52N (Central Meridian $129^\circ\text{E}$).

---

4. Philippine Transverse Mercator (PTM) Grid System

Because the $6^\circ$ zone width of UTM causes excessive scale distortion ($k$ reaching up to $1.00030$ at zone edges) for high-precision cadastral surveys, the Philippine Government established the Philippine Transverse Mercator (PTM) grid system under DENR Administrative Regulations (Manual for Land Surveys of the Philippines / DAO 2007-29).

PTM Grid Specifications

• Zone Width: Divided into 5 narrow zones, each spanning $2^\circ$ of longitude ($1^\circ$ on either side of the central meridian).
• Scale Factor at Central Meridian ($k_0$): $0.99995$ (a scale reduction of only $1/20,000$, yielding far higher precision across the zone).
• False Easting ($E_0$): $500,000.00\text{ meters}$ at the central meridian.
• False Northing ($N_0$): $0.00\text{ meters}$ at the Equator.
• Reference Ellipsoid / Datum: Historical surveys used the Clarke 1866 ellipsoid (Luzon Datum 1911). Modern geodetic operations use the Philippine Reference System of 1992 (PRS92), which retains the Clarke 1866 ellipsoid and is connected to WGS84 through a 7-parameter datum transformation.

Complete PTM Zone Layout

Structural Comparison: UTM vs. PTM

---

5. Grid Corrections and Numerical Calculations

To transform terrestrial field measurements into PTM grid coordinates, Geodetic Engineers apply two fundamental corrections: Grid Convergence and Grid Scale Factor.

Grid Convergence ($\gamma$)

Grid Convergence is the angular difference between True Geographic North and PTM Grid North at a given point. It is approximated by the mathematical expression:

$$\gamma \approx (\lambda - \lambda_0) \sin \phi$$

where $\lambda$ is point longitude, $\lambda_0$ is zone central meridian longitude, and $\phi$ is point latitude.

Worked Example 3: Grid Convergence Calculation

Calculate the grid convergence angle $\gamma$ for a survey station located in Quezon City with coordinates Latitude $\phi = 14^\circ 39'\text{ N}$ ($14.65^\circ$) and Longitude $\lambda = 121^\circ 03'\text{ E}$ ($121.05^\circ$) within PTM Zone III (Central Meridian $\lambda_0 = 121^\circ 00'\text{ E}$).

$$\text{Solution:}$$ $$\Delta \lambda = \lambda - \lambda_0 = 121.05^\circ - 121.00^\circ = +0.05^\circ = +3' = 180''$$ $$\gamma \approx (+0.05^\circ) \times \sin(14.65^\circ) = +0.05^\circ \times 0.25293 = +0.0126465^\circ$$ $$\gamma \text{ in arcseconds} = +0.0126465^\circ \times 3600''/^\circ = +45.53'' \text{ (East of True North)}$$

Point Grid Scale Factor ($k$) and Elevation Factor ($h_{\text{factor}}$)

The point scale factor $k$ at an Easting coordinate $E$ is calculated using the formula:

$$k = k_0 \left[ 1 + \frac{(E - E_0)^2}{2 R_m^2} \right]$$

where $k_0 = 0.99995$, $E_0 = 500,000\text{ m}$, and $R_m \approx 6,371,000\text{ m}$ is the mean radius of curvature of the ellipsoid.

The Elevation Reduction Factor ($h_{\text{factor}}$) reduces ground distance at mean elevation $H$ to the reference ellipsoid:

$$h_{\text{factor}} = \frac{R_m}{R_m + H}$$

The Combined Scale Factor (CSF) is the product of the Point Scale Factor and the Elevation Factor:

$$\text{CSF} = k \times h_{\text{factor}}$$ $$\text{Grid Distance} = \text{Ground Distance} \times \text{CSF}$$ $$\text{Ground Distance} = \frac{\text{Grid Distance}}{\text{CSF}}$$

Worked Example 4: Grid Distance Transformation

A boundary baseline measured with a Total Station on a plateau at an average elevation $H = 500\text{ meters}$ yields a horizontal ground distance of $1,250.00\text{ meters}$. The midpoint of the baseline has a PTM Easting of $E = 560,000\text{ meters}$ in Zone III. Calculate the corresponding PTM Grid Distance.

$$\text{Solution:}$$ $$\text{1. Calculate Point Scale Factor } k\text{:}$$ $$E - E_0 = 560,000 - 500,000 = 60,000\text{ m}$$ $$k = 0.99995 \left[ 1 + \frac{(60,000)^2}{2 \times (6,371,000)^2} \right] = 0.99995 \left[ 1 + \frac{3.6 \times 10^9}{8.118 \times 10^{13}} \right]$$ $$k = 0.99995 \times [1 + 0.000044346] = 0.99995 \times 1.000044346 = 0.99999435$$ $$\text{2. Calculate Elevation Factor } h_{\text{factor}}\text{:}$$ $$h_{\text{factor}} = \frac{6,371,000}{6,371,000 + 500} = \frac{6,371,000}{6,371,500} = 0.99992152$$ $$\text{3. Calculate Combined Scale Factor (CSF):}$$ $$\text{CSF} = k \times h_{\text{factor}} = 0.99999435 \times 0.99992152 = 0.99991587$$ $$\text{4. Calculate PTM Grid Distance:}$$ $$\text{Grid Distance} = 1,250.00\text{ m} \times 0.99991587 = 1,249.895\text{ m} \approx 1,249.90\text{ meters}$$