18-Geom-A3 Geodesy and Positioning · December 2018
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Paper format: National Exams, December 2018 — 3 hours, closed book (approved Casio/Sharp calculators only). SEVEN numbered questions; six constitute a complete paper and each is of equal value (20 marks). Most answers are required in essay format; clarity and organization are explicitly marked. All seven questions are solved below for completeness.
Reference texts: Vaníček & Krakiwsky, Geodesy: The Concepts (2nd ed., North-Holland); Hofmann-Wellenhof, Lichtenegger & Wasle, GNSS — Global Navigation Satellite Systems (Springer, 2008); Torge & Müller, Geodesy (4th ed., de Gruyter); Heiskanen & Moritz, Physical Geodesy (Freeman); Snyder, Map Projections — A Working Manual (USGS PP 1395); Ghilani & Wolf, Elementary Surveying (15th ed.); Natural Resources Canada geodetic references for NAD83(CSRS), CGVD2013, the CGG2013 geoid model and the CACS/CBN networks. Canadian datums/regulators throughout (NRCan, Canadian Geodetic Survey).
Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.
Given. A least-squares network adjustment yielding a covariance matrix for each station, from which error ellipses (2-D) or ellipsoids (3-D) are formed.
Find. (a) the definition and confidence level of the standard error ellipse/ellipsoid; (b) how to scale the standard ellipse to 95%.
(a) Standard error ellipse / ellipsoid. After adjustment, each station has a \(2\times2\) (or \(3\times3\)) covariance sub-matrix \(\Sigma\). The standard error ellipse is the region whose semi-axes are the square roots of the eigenvalues of \(\Sigma\), oriented along its eigenvectors; the semi-major axis \(a\) lies in the direction of maximum positional uncertainty and the semi-minor axis \(b\) in the direction of minimum uncertainty. It is the “\(k=1\)” (one-sigma) ellipse. The standard error ellipsoid is the 3-D analogue, with three semi-axes from the eigenvalues of the \(3\times3\) covariance matrix. The confidence level they define is not the familiar 1-D 68% value: in 2-D the standard ellipse encloses only about 39.4% probability (\(1-e^{-1/2}\)), and in 3-D the standard ellipsoid encloses about 19.9%. These low percentages are exactly why a scaled (e.g. 95%) region is normally reported.
(b) From standard to 95% ellipse. The shape and orientation stay the same; only the size scales.
Approach. Multiply both semi-axes by \(c=\sqrt{\chi^2_{2,0.95}}\), the square root of the 95% quantile of the chi-square distribution with 2 degrees of freedom.
For a 3-D ellipsoid the corresponding factor is \(\sqrt{\chi^2_{3,0.95}}=\sqrt{7.815}=2.796\). These are out-of-context (one-point-at-a-time) factors for a known a-priori variance factor. If the a-posteriori variance factor \(\hat\sigma_0^2\) (with \(df\) degrees of freedom) is used to scale the covariance matrix, the factor becomes \(\sqrt{2F_{2,\,df,\,0.95}}\) instead. A simultaneous (in-context) region that must hold for all network points at once uses the larger factor \(\sqrt{\chi^2_{u,\,0.95}}\), with \(u\) the total number of network coordinates.
| Region | Enclosed probability / scale factor |
|---|---|
| Standard error ellipse (2-D, \(k=1\)) | 39.4% (\(1-e^{-1/2}\)) |
| Standard error ellipsoid (3-D, \(k=1\)) | 19.9% |
| Standard → 95% ellipse (2-D) | × \(\sqrt{\chi^2_{2,0.95}} = 2.448\) |
| Standard → 95% ellipsoid (3-D) | × \(\sqrt{\chi^2_{3,0.95}} = 2.796\) |