Study guide
3.1 Advanced Algebra, Analytic Geometry, and Solid Geometry in Geodetic Calculations
PRC Geodetic Engineering Licensure Examination Exam
Track your progress as you read
Finished sections get a check mark.
1 Chapter 1: Laws, Professional Practice, and Code of Ethics
1.1 Republic Act 8560 & RA 9200 (Philippine Geodetic Engineering Act) and PRC Board Rules 1.2 Code of Ethics for Geodetic Engineers, Professional Responsibilities & Practice Limits 1.3 PRC Licensure Examination Rules, Continuing Professional Development (CPD), and Legal Penalties
2 Chapter 2: Property, Land Registration, and Natural Resources Laws
2.1 Presidential Decree 1529 (Property Registration Decree) & The Torrens System 2.2 Commonwealth Act 141 (Public Land Act), Land Classification, and Public Land Patents 2.3 Comprehensive Agrarian Reform Law (RA 6657) & Environmental/Forestry/Mining Laws (RA 7942, PD 705) 2.4 Civil Code Provisions on Ownership, Property Boundaries, Easements, Obligations, and Contracts 2.5 DENR Administrative Order DAO 2007-29 (Revised Regulations for Land Surveys) & Technical Standards
3 Chapter 3: Mathematics, Geometry, and Trigonometry
3.1 Advanced Algebra, Analytic Geometry, and Solid Geometry in Geodetic Calculations 3.2 Plane and Spherical Trigonometry for Surveying and Astronomical Applications 3.3 Differential and Integral Calculus Applied to Curves, Surfaces, and Surveying Computations
4 Chapter 4: Engineering Sciences, Mechanics, and Economics
4.1 Engineering Mechanics (Statics, Kinematics, Dynamics) for Engineering Applications 4.2 Engineering Economics, Time Value of Money, Cash Flows, Depreciation, and Capital Costing
5 Chapter 5: Theory of Errors and Adjustment of Observations
5.1 Theory of Errors, Precision vs. Accuracy, Random/Systematic Errors, and Error Propagation 5.2 Principles of Least Squares, Most Probable Value, and Weighted Observations 5.3 Least-Squares Adjustment of Level Nets, Triangulation, and Plane Traverses
6 Chapter 6: Plane & Distance/Angle Surveying Fundamentals
6.1 Distance Measurement Methods, Taping Corrections (Temperature, Tension, Sag), and EDM Operations 6.2 Angular Measurements, Transit/Theodolite/Total Station Operations, Bearings, Azimuths, and Magnetic Declination 6.3 Differential, Trigonometric, Barometric, and Profile Leveling Methods
7 Chapter 7: Cadastral, Property, and Special Surveys
7.1 Traverse Computations, Latitudes, Departures, Error of Closure, and Transit/Compass Rule Adjustments 7.2 Area Computations (DMD, DPD, Coordinates) and Subdivision/Consolidation Survey Standards 7.3 Mineral and Mining Surveys, Hydrographic Surveying, Route Alignment, and Construction Layout
8 Chapter 8: Geodetic Surveying and Reference Datums
8.1 Ellipsoids, Geoid Undulation, Geodetic Datums, and the Philippine Reference System of 1992 (PRS92) 8.2 Geodetic Position Computations: Direct and Inverse Problems on the Reference Ellipsoid 8.3 Triangulation, Trilateration, Baseline Measurements, and National Geodetic Control Networks
9 Chapter 9: Geodetic Astronomy and Satellite Positioning (GNSS)
9.1 Celestial Sphere, Equatorial/Horizon Coordinate Systems, and Astronomical Time Systems (UT, ST, LMT) 9.2 Determination of Latitude, Longitude, and Solar/Polaris Azimuth by Astronomical Observations 9.3 Global Navigation Satellite Systems (GNSS/GPS), Carrier-Phase Positioning, RTK, and Static Baselines
10 Chapter 10: Cartography, Map Projections, and Coordinate Systems
10.1 Cartographic Principles, Map Scales, Map Symbolization, and Cartographic Generalization 10.2 Map Projections: Universal Transverse Mercator (UTM) and Philippine Transverse Mercator (PTM) Grid Systems 10.3 Survey Plan Preparation, Technical Description Formatting, and DENR Approval Workflows
11 Chapter 11: Photogrammetry, Remote Sensing, and Geographic Information Systems (GIS)
11.1 Aerial Photogrammetry, Photographic Scale, Relief Displacement, Flight Planning, and Stereoscopy 11.2 Digital Photogrammetry, UAV/Drone Surveying, Orthophoto Generation, and LiDAR Data Processing 11.3 Remote Sensing Spectral Signatures, Multispectral Image Classification, and GIS Spatial Analysis
3.1 Advanced Algebra, Analytic Geometry, and Solid Geometry in Geodetic Calculations
Key Takeaways Explain This Topic AI
Polynomial, quadratic, logarithmic, and exponential functions form the mathematical basis for geodetic coordinate transformations, atmospheric refraction modeling, and height decay functions.
Spatial vector representation of lines and planes in 3D coordinate systems enables direct modeling of total station sightlines, satellite vector geometry, and plane intersections.
Conic sections—specifically ellipses—define reference ellipsoid geometry (semi-major axis a, semi-minor axis b, flattening f, and first eccentricity e^2), satellite Keplerian orbits, and reflector optics.
Conformal coordinate transformations (2D similarity and 3D Helmert 7-parameter transformations) allow seamless conversion of coordinates between local geodetic datums (such as PRS92) and global satellite reference frames (such as WGS84 and ITRF).
Solid geometry formulas for polyhedra, spherical caps, and ellipsoidal segments provide exact volumetric and surface area calculations for geodetic reservoir surveys, spherical earth modeling, and earthwork estimations.
Last updated: July 2026
3.1 Advanced Algebra, Analytic Geometry, and Solid Geometry in Geodetic Calculations
Geodetic engineering calculations rely heavily on rigorous mathematical foundations spanning advanced algebra, spatial analytic geometry, conic sections, coordinate transformations, and 3D solid geometry. Whether computing position vectors from GNSS satellites, transforming local grid coordinates to national reference systems, or modeling the geometric shape of the Earth, geodetic engineers must apply these mathematical principles with exact precision.
Advanced Algebra Applications in Geodetic Engineering
Algebraic equations provide the essential structure for modeling physical and geometric phenomena in surveying and geodesy.
Polynomial and Quadratic Functions
Quadratic equations ( a x 2 + b x + c = 0 ax^2 + bx + c = 0 a x 2 + b x + c = 0 ) appear routinely in curve fitting, intersection of lines and circular arcs, and vertical road curve design. Higher-order polynomials are utilized in empirical models of terrain surfaces and gravity anomalies. The general solution for a quadratic equation is given by: x = − b ± b 2 − 4 a c 2 a x = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a} x = 2 a − b ± b 2 − 4 a c
In geodetic adjustments, quadratic expressions define distance constraints between points P 1 ( x 1 , y 1 ) P_1(x_1, y_1) P 1 ( x 1 , y 1 ) and P 2 ( x 2 , y 2 ) P_2(x_2, y_2) P 2 ( x 2 , y 2 ) : d 2 = ( x 2 − x 1 ) 2 + ( y 2 − y 1 ) 2 d^2 = (x_2 - x_1)^2 + (y_2 - y_1)^2 d 2 = ( x 2 − x 1 ) 2 + ( y 2 − y 1 ) 2
Logarithmic and Exponential Functions
Logarithmic and exponential functions describe atmospheric density variations, signal attenuation in electromagnetic distance measurement (EDM), and barometric leveling formulas. Atmospheric refraction corrections for zenith angles rely on exponential decay models of air density with altitude h h h : ρ ( h ) = ρ 0 e − M g h R T \rho(h) = \rho_0 e^{-\frac{M g h}{R T}} ρ ( h ) = ρ 0 e − RT M g h
where ρ 0 \rho_0 ρ 0 is sea-level air density, M M M is molar mass of dry air, g g g is gravitational acceleration, R R R is the universal gas constant, and T T T is absolute temperature in Kelvin.
3D Analytic Geometry: Lines, Planes, and Spatial Vectors
Modern geodetic systems operate in three-dimensional Cartesian space ( X , Y , Z ) (X, Y, Z) ( X , Y , Z ) , where spatial positions are defined relative to the Earth's center of mass (Geocentric Coordinate System).
Spatial Vectors and 3D Line Equations
A line passing through a point P 0 ( X 0 , Y 0 , Z 0 ) P_0(X_0, Y_0, Z_0) P 0 ( X 0 , Y 0 , Z 0 ) with direction vector v = ⟨ l , m , n ⟩ \mathbf{v} = \langle l, m, n \rangle v = ⟨ l , m , n ⟩ is expressed in parametric form as: X = X 0 + l t , Y = Y 0 + m t , Z = Z 0 + n t X = X_0 + l t, \quad Y = Y_0 + m t, \quad Z = Z_0 + n t X = X 0 + lt , Y = Y 0 + m t , Z = Z 0 + n t
where t t t is a scalar parameter. Symmetric form for the 3D line (when l , m , n ≠ 0 l, m, n \neq 0 l , m , n = 0 ) is: X − X 0 l = Y − Y 0 m = Z − Z 0 n \frac{X - X_0}{l} = \frac{Y - Y_0}{m} = \frac{Z - Z_0}{n} l X − X 0 = m Y − Y 0 = n Z − Z 0
Planes in 3D Space
The general equation of a plane in 3D space is: A X + B Y + C Z + D = 0 A X + B Y + C Z + D = 0 A X + B Y + CZ + D = 0
where N = ⟨ A , B , C ⟩ \mathbf{N} = \langle A, B, C \rangle N = ⟨ A , B , C ⟩ represents the normal vector perpendicular to the plane. The perpendicular distance d d d from an arbitrary point P 1 ( X 1 , Y 1 , Z 1 ) P_1(X_1, Y_1, Z_1) P 1 ( X 1 , Y 1 , Z 1 ) to the plane is given by: d = ∣ A X 1 + B Y 1 + C Z 1 + D ∣ A 2 + B 2 + C 2 d = \frac{|A X_1 + B Y_1 + C Z_1 + D|}{\sqrt{A^2 + B^2 + C^2}} d = A 2 + B 2 + C 2 ∣ A X 1 + B Y 1 + C Z 1 + D ∣
Worked Example: Perpendicular Distance to a Reference Datum Plane
Problem: A total station instrument offset point is located at spatial coordinates P ( 100.00 , 200.00 , 50.00 ) P(100.00, 200.00, 50.00) P ( 100.00 , 200.00 , 50.00 ) meters. A local datum reference plane is defined by the linear equation 3 X + 4 Y + 12 Z − 1300 = 0 3X + 4Y + 12Z - 1300 = 0 3 X + 4 Y + 12 Z − 1300 = 0 . Calculate the perpendicular distance from point P P P to the reference plane.
Solution:
Identify plane coefficients: A = 3 A = 3 A = 3 , B = 4 B = 4 B = 4 , C = 12 C = 12 C = 12 , D = − 1300 D = -1300 D = − 1300 .
Substitute point coordinates X 1 = 100.00 X_1 = 100.00 X 1 = 100.00 , Y 1 = 200.00 Y_1 = 200.00 Y 1 = 200.00 , Z 1 = 50.00 Z_1 = 50.00 Z 1 = 50.00 into the distance equation: d = ∣ 3 ( 100.00 ) + 4 ( 200.00 ) + 12 ( 50.00 ) − 1300 ∣ 3 2 + 4 2 + 12 2 d = \frac{|3(100.00) + 4(200.00) + 12(50.00) - 1300|}{\sqrt{3^2 + 4^2 + 12^2}} d = 3 2 + 4 2 + 1 2 2 ∣3 ( 100.00 ) + 4 ( 200.00 ) + 12 ( 50.00 ) − 1300∣
Calculate the numerator: ∣ 300 + 800 + 600 − 1300 ∣ = ∣ 1700 − 1300 ∣ = 400 |300 + 800 + 600 - 1300| = |1700 - 1300| = 400 ∣300 + 800 + 600 − 1300∣ = ∣1700 − 1300∣ = 400 .
Calculate the denominator: 9 + 16 + 144 = 169 = 13 \sqrt{9 + 16 + 144} = \sqrt{169} = 13 9 + 16 + 144 = 169 = 13 .
Compute distance: d = 400 13 ≈ 30.77 meters d = \frac{400}{13} \approx 30.77 \text{ meters} d = 13 400 ≈ 30.77 meters .
Conic Sections and Geodetic Reference Ellipsoids
Conic sections—circles, ellipses, parabolas, and hyperbolas—are generated by intersecting a plane with a double circular cone. Conic Section Types:
- 1. Circle: Plane perpendicular to cone axis (e = 0)
- 2. Ellipse: Plane angled to axis, cutting all generators (0 < e < 1)
- 3. Parabola: Plane parallel to cone generator line (e = 1)
- 4. Hyperbola: Plane parallel to cone axis, cutting both lobes (e > 1)
- The Ellipse as a Reference Geodetic Surface
An ellipse is defined by the standard cartesian equation: x 2 a 2 + y 2 b 2 = 1 \frac{x^2}{a^2} + \frac{y^2}{b^2} = 1 a 2 x 2 + b 2 y 2 = 1
where a a a is the semi-major axis and b b b is the semi-minor axis. Rotating an ellipse about its minor axis generates an oblate spheroid (reference ellipsoid), which serves as the geometric surface approximating Earth's shape.
- Key Ellipsoidal Parameters
- Flattening ( f f f ): f = a − b a f = \frac{a - b}{a} f = a a − b
First Eccentricity Squared ( e 2 e^2 e 2 ): e 2 = a 2 − b 2 a 2 = 2 f − f 2 e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2 e 2 = a 2 a 2 − b 2 = 2 f − f 2
Second Eccentricity Squared ( e ′ 2 e'^2 e ′2 ): e ′ 2 = a 2 − b 2 b 2 = e 2 1 − e 2 e'^2 = \frac{a^2 - b^2}{b^2} = \frac{e^2}{1 - e^2} e ′2 = b 2 a 2 − b 2 = 1 − e 2 e 2
Comparison of Major Reference Ellipsoids
Ellipsoid Parameter WGS 84 (Global GNSS Datum) Clarke 1866 (PRS92 Datum)
Semi-major axis ( a a a ) 6 , 378 , 137.0 m 6,378,137.0 \text{ m} 6 , 378 , 137.0 m 6 , 378 , 206.4 m 6,378,206.4 \text{ m} 6 , 378 , 206.4 m
Semi-minor axis ( b b b ) 6 , 356 , 752.3142 m 6,356,752.3142 \text{ m} 6 , 356 , 752.3142 m 6 , 356 , 583.8000 m 6,356,583.8000 \text{ m} 6 , 356 , 583.8000 m
Reciprocal Flattening ( 1 / f 1/f 1/ f ) 298.257223563 298.257223563 298.257223563 294.9786982 294.9786982 294.9786982
First Eccentricity Squared ( e 2 e^2 e 2 ) 0.00669437999014 0.00669437999014 0.00669437999014 0.00676865799729 0.00676865799729 0.00676865799729
Coordinate Transformations: 2D Conformal & 3D Helmert 7-Parameter
Coordinate transformations convert spatial position representations between different coordinate systems, map projections, or geodetic datums.
2D Similarity (Helmert 4-Parameter) Transformation
A 2D conformal transformation preserves angles and shape while applying scale adjustment s s s , rotation angle θ \theta θ , and translations ( X 0 , Y 0 ) (X_0, Y_0) ( X 0 , Y 0 ) : [ X Y ] = [ c − d d c ] [ x y ] + [ X 0 Y 0 ] \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} c & -d \\ d & c \end{bmatrix} \begin{bmatrix} x \\ y \end{bmatrix} + \begin{bmatrix} X_0 \\ Y_0 \end{bmatrix} [ X Y ] = [ c d − d c ] [ x y ] + [ X 0 Y 0 ]
where transformation coefficients c c c and d d d relate to scale and rotation as: c = s cos θ , d = s sin θ c = s \cos \theta, \quad d = s \sin \theta c = s cos θ , d = s sin θ s = c 2 + d 2 , θ = arctan ( d c ) s = \sqrt{c^2 + d^2}, \quad \theta = \arctan\left(\frac{d}{c}\right) s = c 2 + d 2 , θ = arctan ( c d )
Worked Example: 2D Helmert Transformation Parameters
Problem: Derived transformation parameters for converting local survey grid coordinates ( x , y ) (x, y) ( x , y ) to national grid coordinates ( X , Y ) (X, Y) ( X , Y ) are c = 0.99980 c = 0.99980 c = 0.99980 and d = 0.01745 d = 0.01745 d = 0.01745 . Determine the scale factor s s s and rotation angle θ \theta θ .
Solution:
Calculate scale factor s s s : s = ( 0.99980 ) 2 + ( 0.01745 ) 2 = 0.99960004 + 0.00030450 = 0.99990454 ≈ 0.99995 s = \sqrt{(0.99980)^2 + (0.01745)^2} = \sqrt{0.99960004 + 0.00030450} = \sqrt{0.99990454} \approx 0.99995 s = ( 0.99980 ) 2 + ( 0.01745 ) 2 = 0.99960004 + 0.00030450 = 0.99990454 ≈ 0.99995
Calculate rotation angle θ \theta θ : θ = arctan ( 0.01745 0.99980 ) = arctan ( 0.0174535 ) ≈ 1.00 ∘ \theta = \arctan\left(\frac{0.01745}{0.99980}\right) = \arctan(0.0174535) \approx 1.00^\circ θ = arctan ( 0.99980 0.01745 ) = arctan ( 0.0174535 ) ≈ 1.0 0 ∘
3D Helmert 7-Parameter Transformation
Converting 3D Cartesian coordinates ( x , y , z ) (x, y, z) ( x , y , z ) from one datum (e.g., WGS84) to another (e.g., PRS92) requires 3 translations ( Δ X , Δ Y , Δ Z ) (\Delta X, \Delta Y, \Delta Z) ( Δ X , Δ Y , Δ Z ) , 3 differential rotations ( R x , R y , R z ) (R_x, R_y, R_z) ( R x , R y , R z ) in radians, and 1 scale factor factor s = 1 + δ s s = 1 + \delta s s = 1 + δs : [ X Y Z ] = [ Δ X Δ Y Δ Z ] + ( 1 + δ s ) [ 1 − R z R y R z 1 − R x − R y R x 1 ] [ x y z ] \begin{bmatrix} X \\ Y \\ Z \end{bmatrix} = \begin{bmatrix} \Delta X \\ \Delta Y \\ \Delta Z \end{bmatrix} + (1 + \delta s) \begin{bmatrix} 1 & -R_z & R_y \\ R_z & 1 & -R_x \\ -R_y & R_x & 1 \end{bmatrix} \begin{bmatrix} x \\ y \\ z \end{bmatrix} X Y Z = Δ X Δ Y Δ Z + ( 1 + δs ) 1 R z − R y − R z 1 R x R y − R x 1 x y z
Solid Geometry: Polyhedra, Spheres, and Ellipsoidal Segments
Solid geometry provides volume and surface area calculations for physical terrain models, earthwork prismoids, and spherical/ellipsoidal earth segments.
Spherical Caps and Spherical Zones
Surface Area of Spherical Zone ( S S S ): S = 2 π R h S = 2 \pi R h S = 2 π R h where R R R is the sphere radius and h h h is the height of the spherical segment.
Volume of Spherical Segment ( V V V ): V = π h 2 3 ( 3 R − h ) V = \frac{\pi h^2}{3} (3R - h) V = 3 π h 2 ( 3 R − h )
Volume of Ellipsoid of Revolution
The total volume of an oblate ellipsoid of revolution with semi-axes a a a and b b b is given by: V = 4 3 π a 2 b V = \frac{4}{3} \pi a^2 b V = 3 4 π a 2 b