3.1 Advanced Algebra, Analytic Geometry, and Solid Geometry in Geodetic Calculations

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 ($ax^2 + bx + 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 = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}$

In geodetic adjustments, quadratic expressions define distance constraints between points $P_1(x_1, y_1)$ and $P_2(x_2, y_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$:

$\rho(h) = \rho_0 e^{-\frac{M g h}{R T}}$

where $\rho_0$ is sea-level air density, $M$ is molar mass of dry air, $g$ is gravitational acceleration, $R$ is the universal gas constant, and $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)$, 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)$ with direction vector $\mathbf{v} = \langle l, m, n \rangle$ is expressed in parametric form as:

$X = X_0 + l t, \quad Y = Y_0 + m t, \quad Z = Z_0 + n t$

where $t$ is a scalar parameter. Symmetric form for the 3D line (when $l, m, n \neq 0$) is:

$\frac{X - X_0}{l} = \frac{Y - Y_0}{m} = \frac{Z - Z_0}{n}$

Planes in 3D Space

The general equation of a plane in 3D space is:

$A X + B Y + C Z + D = 0$

where $\mathbf{N} = \langle A, B, C \rangle$ represents the normal vector perpendicular to the plane. The perpendicular distance $d$ from an arbitrary point $P_1(X_1, Y_1, Z_1)$ to the plane is given by:

$d = \frac{|A X_1 + B Y_1 + C Z_1 + D|}{\sqrt{A^2 + B^2 + C^2}}$

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)$ meters. A local datum reference plane is defined by the linear equation $3X + 4Y + 12Z - 1300 = 0$. Calculate the perpendicular distance from point $P$ to the reference plane.

Solution:

1. Identify plane coefficients: $A = 3$, $B = 4$, $C = 12$, $D = -1300$.
2. Substitute point coordinates $X_1 = 100.00$, $Y_1 = 200.00$, $Z_1 = 50.00$ into the distance equation:

$d = \frac{|3(100.00) + 4(200.00) + 12(50.00) - 1300|}{\sqrt{3^2 + 4^2 + 12^2}}$

1. Calculate the numerator: $|300 + 800 + 600 - 1300| = |1700 - 1300| = 400$.
2. Calculate the denominator: $\sqrt{9 + 16 + 144} = \sqrt{169} = 13$.
3. Compute distance: $d = \frac{400}{13} \approx 30.77 \text{ 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:

$\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1$

where $a$ is the semi-major axis and $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 = \frac{a - b}{a}$
• First Eccentricity Squared ($e^2$): $e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2$
• Second Eccentricity Squared ($e'^2$): $e'^2 = \frac{a^2 - b^2}{b^2} = \frac{e^2}{1 - e^2}$

Comparison of Major Reference Ellipsoids

---

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$, rotation angle $\theta$, and translations $(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}$

where transformation coefficients $c$ and $d$ relate to scale and rotation as:

$c = s \cos \theta, \quad d = s \sin \theta$ $s = \sqrt{c^2 + d^2}, \quad \theta = \arctan\left(\frac{d}{c}\right)$

Worked Example: 2D Helmert Transformation Parameters

Problem: Derived transformation parameters for converting local survey grid coordinates $(x, y)$ to national grid coordinates $(X, Y)$ are $c = 0.99980$ and $d = 0.01745$. Determine the scale factor $s$ and rotation angle $\theta$.

Solution:

1. Calculate scale factor $s$: $s = \sqrt{(0.99980)^2 + (0.01745)^2} = \sqrt{0.99960004 + 0.00030450} = \sqrt{0.99990454} \approx 0.99995$
2. Calculate rotation angle $\theta$: $\theta = \arctan\left(\frac{0.01745}{0.99980}\right) = \arctan(0.0174535) \approx 1.00^\circ$

3D Helmert 7-Parameter Transformation

Converting 3D Cartesian coordinates $(x, y, z)$ from one datum (e.g., WGS84) to another (e.g., PRS92) requires 3 translations $(\Delta X, \Delta Y, \Delta Z)$, 3 differential rotations $(R_x, R_y, R_z)$ in radians, and 1 scale factor factor $s = 1 + \delta s$:

$\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}$

---

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 = 2 \pi R h$ where $R$ is the sphere radius and $h$ is the height of the spherical segment.
• Volume of Spherical Segment ($V$): $V = \frac{\pi h^2}{3} (3R - h)$

Volume of Ellipsoid of Revolution

The total volume of an oblate ellipsoid of revolution with semi-axes $a$ and $b$ is given by:

$V = \frac{4}{3} \pi a^2 b$