GELE Geodesy — Geodetic and Cartesian CoordinatesSummary

Geodetic latitude φ is the angle from the equator to a point on the ellipsoid along the meridian, measured at the ellipsoid's center. Geodetic longitude λ is the angle from the Greenwich meridian to the meridian through the point. Ellipsoidal height h is the perpendicular distance from the ellipsoid surface to the point. For Philippine work, these are typically referenced to WGS84 (global) or PRS92 (national datum adopted by the Land Registration Authority under RA 4374). Unlike geocentric latitude (used in older systems), geodetic latitude accounts for the Earth's flattening and is what surveyors and GNSS receivers report.

A three-dimensional rectangular coordinate system with origin at the Earth's center of mass. The Z-axis points toward the North Pole, the X-axis toward the intersection of the Greenwich meridian and the equator, and the Y-axis completes a right-handed system (toward 90° East). GNSS satellites transmit their orbits and Earth's shape in this frame; receivers compute positions by trilateration in X, Y, Z space. The system is universal—the same for all points on Earth—and is independent of local projections or datums.

The radius of curvature of the ellipsoid in the prime vertical (east–west) plane. Defined as N = a / √(1 − e² sin² φ), where a is the semi-major axis and e² is the first eccentricity squared. N varies with latitude: at the equator N = a, and at the poles N approaches a/(1−e²) = a(1/(1−f)), where f is flattening. For WGS84: a = 6,378,137 m, e² = 0.00669438. N represents the distance from a point on the ellipsoid (in the prime vertical plane) to the Z-axis. It appears in all forward and inverse transformation formulas and must be recalculated for each latitude iteration in the inverse problem.

Given (φ, λ, h), compute (X, Y, Z) using: X = (N + h) cos φ cos λ, Y = (N + h) cos φ sin λ, Z = [N(1 − e²) + h] sin φ. The factor (N + h) is the distance from the Z-axis to the point, lying in the equatorial plane; it is then projected onto the X and Y axes using cos and sin of the longitude. The Z component uses N(1 − e²) instead of N to account for the ellipsoid's oblateness (flattening); omitting (1 − e²) would incorrectly model a sphere. This transformation is direct (non-iterative) and fast, making it ideal for batch processing of survey data.

Given (X, Y, Z), compute (φ, λ, h). Longitude is immediate: λ = atan2(Y, X), but requires resolving the quadrant (quadrant I: X>0, Y>0 gives λ in [0°, 90°]; quadrant II: X<0, Y>0 gives λ in [90°, 180°], etc.). Latitude and height are coupled through N and typically solved iteratively: (1) Initialize p = √(X² + Y²) and φ₀ = atan(Z / p(1 − e²)). (2) Compute N(φᵢ) and hᵢ = p / cos φᵢ − N. (3) Update φᵢ₊₁ = atan(Z / p · (1 − e² N / (N + h))⁻¹). Repeat steps 2–3 until |φᵢ₊₁ − φᵢ| < ε (e.g., 10⁻¹² radians). Alternatively, closed-form solutions (e.g., Bowring's or Heikkinen's method) avoid iteration but are algebraically complex.

Ellipsoidal height h is the perpendicular distance from the WGS84/PRS92 ellipsoid surface to the point. It is what GNSS receivers measure directly. Orthometric height H is the distance measured vertically (along the direction of gravity) from the geoid (mean sea level surface) to the point. They are related by the geoid undulation N_geoid: H = h − N_geoid. For Philippine land titles and cadastral surveys, orthometric heights (measured by traditional levelling) are often required for land value assessment and flood-risk mapping. GNSS-derived ellipsoidal heights must be converted using a geoid model (e.g., PRS92-specific geoid grid). This distinction is critical when mixing GNSS and levelling data.

The Earth's ellipsoid is defined by its semi-major axis a and either its eccentricity e or flattening f. Eccentricity e = √(1 − b²/a²), where b is the semi-minor axis; e² = 0.00669438 for WGS84. Flattening f = (a − b) / a ≈ 1/298.257 for WGS84. These parameters quantify the ellipsoid's deviation from a sphere and appear in all geodetic calculations. The factor (1 − e²) ≈ 0.9933 in the Z formula reflects the flattening's effect. Understanding the relationship between a, b, e, and f is essential for working with different datums (e.g., PRS92 vs. WGS84) or historical ellipsoids.

The arctangent function atan(Y/X) returns an angle in (−π/2, π/2) radians or (−90°, 90°), ambiguous regarding which of the four quadrants the point lies in. To resolve the quadrant, use the signs of X and Y: Quadrant I (X>0, Y>0): λ = atan(Y/X); Quadrant II (X<0, Y>0): λ = 180° + atan(Y/X); Quadrant III (X<0, Y<0): λ = −180° + atan(Y/X) or 180° − |atan(Y/X)|; Quadrant IV (X>0, Y<0): λ = atan(Y/X) or 360° + atan(Y/X). Most programming languages provide atan2(Y, X) which handles this automatically. For board exams, always explicitly state which quadrant and show the quadrant-adjustment step.

The inverse transformation cannot be solved in a single closed form for φ and h because they are interdependent (N depends on φ, and h depends on N). The standard iterative method is: φᵢ₊₁ = atan(Z / (p(1 − e² N(φᵢ) / (N(φᵢ) + hᵢ)))). Convergence is rapid (typically 2–4 iterations) because the Earth's ellipsoid is only slightly flattened. A stopping criterion is |φᵢ₊₁ − φᵢ| < 1×10⁻¹² radians (about 3 × 10⁻⁸ m on the Earth's surface). For licensing exams, one full iteration is usually sufficient to demonstrate understanding; in practice, modern software uses closed-form solutions or fast iteration with careful stopping rules.

The Philippine Plane Coordinate System (PPCS) and UTM projections are defined by converting from geodetic coordinates to planar (Easting, Northing) coordinates. GNSS receivers output Cartesian coordinates, which are first converted to geodetic (φ, λ, h on PRS92), then to PPCS or UTM for map sheets and cadastral plans. Land titles registered under the Land Registration Authority (RA 4374, RA 8560) reference these planar coordinates, as do engineering designs for roads, dams, and utilities. Field surveyors use GNSS to establish control points in Cartesian coordinates, which are then transformed to geodetic, and finally to the projection used for the survey region (typically PPCS for large projects, UTM for smaller areas or international consistency under PD 1529).

The semi-major axis a and eccentricity e² of the ellipsoid (WGS84 or PRS92) define the shape of the reference surface. These parameters appear in the formula for N, which in turn governs all forward and inverse transformations. Changes in a or e² (due to refined measurements or choice of datum) directly affect N and propagate into X, Y, Z or φ, λ, h. This relationship underscores the importance of specifying the datum clearly in any surveying project.

These are two representations of the same point in space. The forward transformation uses trigonometric projections and the prime-vertical radius to map from the curved ellipsoid to a rectangular frame. The inverse requires both direct calculation (for λ) and iteration (for φ, h) due to the coupling between latitude and height through N. Neither transformation is approximate; they are exact mathematical relationships. Understanding the duality is essential for working with different geodetic systems and software packages.

In the inverse transformation, latitude φ appears in the formula for N, and N appears in the formula for φ. This circular dependence is resolved by iteration: assume an initial φ₀, compute N(φ₀), use N to refine φ, and repeat. Convergence is fast because the ellipsoid's flattening is small (e² ≈ 0.0067). This relationship explains why the inverse problem is more complex than the forward and why understanding iteration is critical for licensure.

Longitude is determined solely from the signs and magnitudes of X and Y. The raw arctangent atan(Y/X) is ambiguous; the quadrant must be resolved using the signs of X and Y (Quadrant I, II, III, or IV). This relationship reminds students that Cartesian coordinates encode full positional information, whereas geographic coordinates require care in interpreting them. Misreading the quadrant is a frequent exam error with serious practical consequences (e.g., a survey boundary placed 180° away).

Ellipsoidal height (from GNSS) and orthometric height (from levelling) are connected by the geoid undulation: H = h − N_geoid. This relationship bridges geodetic theory and practical surveying. In the Philippines, the geoid undulation varies across the archipelago (roughly −10 to +40 m), so h and H can differ significantly. Understanding this relationship is essential for projects requiring orthometric heights (e.g., flood mapping, infrastructure design) and for reconciling GNSS data with traditional levelled networks.

This is the workflow in modern Philippine surveying. GNSS receivers compute X, Y, Z in real time. Field software then converts to φ, λ, h on PRS92, and finally to PPCS or UTM for map sheets and cadastral plans. Each transformation introduces rounding, but when executed correctly, the final planar coordinates are accurate and traceable to the global datum. This relationship illustrates why mastering both forward and inverse transformations is essential for practicing surveyors.

The iterative inverse solution for φ is deemed converged when |φᵢ₊₁ − φᵢ| < threshold (typically 10⁻¹² rad). This criterion corresponds to a linear precision of about 3×10⁻⁸ m on the Earth's surface, well below practical surveying accuracies (which are typically ±1 cm or larger). Thus, 2–4 iterations guarantee submillimeter precision. This relationship shows why high-precision geodesy requires careful attention to numerical methods and stopping criteria.

A surveyor sets up a GNSS receiver on a property corner in Metro Manila. The receiver outputs X = -3,188,054.9 m, Y = 5,305,814.4 m, Z = 1,532,993.9 m (WGS84). To register the boundary with the Land Registration Authority under RA 4374, the position must be expressed in geodetic coordinates (φ, λ, h) and then converted to the Philippine Plane Coordinate System (PPCS Zone III for the Metro). Using the inverse transformation, the surveyor calculates φ ≈ 14°N, λ ≈ 121°E, h ≈ 50 m. This geodetic position is then mapped to PPCS using the UTM/Transverse Mercator projection, yielding planar coordinates (Easting, Northing) that are written into the land title deed. Without correct coordinate transformations, the title is ambiguous and legally unenforceable.

The Philippine Atmospheric, Geophysical and Astronomical Services Administration (PAGASA) conducts a GNSS survey to map flood elevations in low-lying areas. Raw GNSS data are in WGS84 Cartesian (X, Y, Z). The survey office converts to geodetic coordinates (φ, λ, h on PRS92), but map production requires orthometric heights H (elevation above mean sea level for geotechnical and hydrological meaning). Using the PRS92 geoid model, the office computes H = h − N_geoid. A point with h = 15 m and N_geoid = 12 m yields H = 3 m, showing it is 3 m above sea level (flood risk). Without proper height system conversion, the flood maps would be meaningless and infrastructure at the wrong elevation would be designed.

The National Mapping and Resource Information Authority (NAMRIA) produces 1:50,000-scale topographic map sheets covering the Philippines. Surveys are conducted using GNSS receivers (Cartesian coordinates) and differential positioning techniques. The survey data are post-processed to geodetic coordinates on PRS92. For each map sheet (e.g., 'Marikina Valley, 1 of 4'), the geodetic coordinates are converted to Philippine Plane Coordinate System (PPCS) using Transverse Mercator projection formulas. The planar coordinates are then printed on the map grid, along with geographic graticule (latitude/longitude lines). Surveyors and planners use these maps; coordinate transformations ensure consistency between the map grid, GNSS positioning, and land titles. Any error in forward or inverse transformations would misalign field surveys with published maps.

A boundary dispute arises over a 10-hectare agricultural property in Nueva Ecija. The original title (issued in 1985) references hand-surveyed coordinates in a local grid. A modern survey using GNSS generates coordinates in WGS84 Cartesian, which are converted to geodetic (PRS92), then to PPCS. A careful inverse transformation reveals that the old boundary has shifted, possibly due to the choice of ellipsoid in the original survey (before PRS92 was adopted). The surveyors resolve the dispute by demonstrating the mathematical chain: old coordinates → WGS84/PRS92 geodetic → PPCS → field positions. This chain is auditable and legally defensible. Without understanding coordinate transformations, surveyors cannot reconcile old and new data, leaving disputes unresolved.

The Department of Public Works and Highways (DPWH) designs a 50-km expressway. Survey control points are established via GNSS (Cartesian output). Design engineers receive these as geodetic coordinates and project them onto the design datum (PPCS Zone II or UTM). Horizontal alignment is laid out in the planar system, but vertical alignment (grade elevation) requires orthometric heights. A culvert elevation must be set 0.5 m above the 100-year flood level; this is specified in orthometric height. The survey team converts GNSS h to H using a geoid model, sets the culvert at the correct elevation, and construction proceeds. An error in height-system conversion would either allow flooding (if the culvert is too low) or waste material (if too high). Coordinate transformations are thus safety-critical.

The Philippines and Malaysia cooperate on a marine boundary survey in the Celebes Sea. Malaysian data are referenced to local geodetic datum; Philippine data use PRS92. To align the surveys, both are converted to WGS84 Cartesian, compared for consistency, and a datum transformation (shift, rotation, scale) is computed. The transformation involves forward and inverse conversions for hundreds of control points. Only through rigorous understanding of Cartesian and geodetic coordinates can such international efforts succeed. Errors here can lead to maritime boundary disputes with significant economic and political consequences.

A farmer in Bulacan uses RTK-GNSS guidance for precision farming (e.g., variable-rate fertilizer application). The RTK receiver outputs Cartesian coordinates with 2-cm accuracy in real time. Field software converts these to geodetic (PRS92) and then to PPCS for the local farm plot. Grid-based application maps are then generated, aligned with the farm's existing boundary (from a traditional survey). Without correct Cartesian-to-geodetic-to-planar transformations, the GPS-guided equipment would apply fertilizer to the wrong areas, reducing yield or harming adjacent crops. High-accuracy transformations are essential for modern agricultural technology.

PHIVOLCS operates a network of continuous GNSS stations across the Philippines to monitor tectonic motion (Philippines sits on the boundary of four major tectonic plates). Each station outputs X, Y, Z coordinates at regular intervals. Over months, deformation accumulates; for instance, a station near Mount Pinatubo may drift a few centimeters northwestward per year. Scientists convert time-series Cartesian coordinates to geodetic differences (ΔX, ΔY, ΔZ → Δφ, Δλ, Δh) to quantify deformation rates and identify anomalies that might precede volcanic eruption or earthquakes. Accurate transformations are vital for disaster preparedness; a misinterpretation due to transformation error could lead to false alarms or missed warnings.