NivaarExam PrepOfficial exam papers ↗

18-Geom-A3 Geodesy and Positioning · December 2017

Question 6 of 7: Network Post-analysis — Confidence Ellipses and Ellipsoids

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format: National Exams, December 2017 — 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); the per-part marking scheme printed on page 4 of the paper is reproduced in each answer. Most answers are required in essay format. 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); Natural Resources Canada geodetic references for NAD83(CSRS), CGVD2013 and the CGG2013 geoid model. Canadian datums/regulators throughout (NRCan, Canadian Geodetic Survey).

Question 6: Network Post-analysis — Confidence Ellipses and Ellipsoids (20 marks — a 8, b 8, c 4)

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) the meaning of out-of-context vs in-context (simultaneous) regions; (c) how to scale the standard ellipse to 95%.

a (semi-major)b (semi-minor)θStandard error ellipse — 39.4% confidence in 2-D
Standard (k = 1) error ellipse from the station covariance matrix: semi-major axis a and semi-minor axis b along the eigenvector directions, orientation θ from the eigenvectors of the 2×2 covariance sub-matrix. It encloses about 39.4% probability in 2-D.

(a) Standard error ellipse / ellipsoid. After adjustment, each station has a \(2\times2\) (or \(3\times3\)) covariance sub-matrix. The standard error ellipse is the region whose semi-axes are the square roots of the eigenvalues of that covariance matrix, 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) Out-of-context vs in-context (simultaneous). An out-of-context (or one-at-a-time) ellipse describes the uncertainty of a single point considered in isolation, ignoring that many points were estimated together. An in-context (or simultaneous) ellipse accounts for the fact that we are making a confidence statement about all the network’s points (or parameters) at once: to keep the joint confidence at the stated level, the region must be enlarged. The in-context region is therefore always larger than the out-of-context one for the same confidence, because guarding against error in any of many parameters simultaneously demands a bigger multiplier: with the variance factor known, the out-of-context factor uses \(\chi^2\) with 2 (or 3) degrees of freedom, whereas the in-context factor uses \(\chi^2\) with degrees of freedom equal to the total number of coordinates \(u\) in the network (for an estimated variance factor the \(\chi^2\) value is replaced by the corresponding \(F\)-based value).

(c) From standard to 95% ellipse. The shape and orientation stay the same; only the size scales. Multiply both semi-axes of the standard ellipse by the factor \(c\), where \(c^2\) is the 95% quantile of the \(\chi^2\) distribution with 2 degrees of freedom:

$$c = \sqrt{\chi^2_{2,\,0.95}} = \sqrt{5.991} = \boxed{2.448}.$$

So the 95% error ellipse has semi-axes \(a_{95}=2.448\,a\) and \(b_{95}=2.448\,b\). For example, a standard ellipse with \(a = 12.0\ \text{mm}\), \(b = 7.0\ \text{mm}\) becomes \(a_{95}=29.4\ \text{mm}\), \(b_{95}=17.1\ \text{mm}\) at 95% confidence. (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 factors with a known a-priori variance factor; with the a-posteriori variance factor \(\hat\sigma_0^2\) (\(\nu\) degrees of freedom) the factor becomes \(\sqrt{2F_{2,\nu,0.95}}\), and the simultaneous (in-context) 95% region uses \(\sqrt{\chi^2_{u,0.95}}\) over all \(u\) network coordinates — e.g. 5.61 for a 10-station 2-D network (\(u=20\)).

RegionEnclosed 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\)