Guideline For Adjustment and Evaluation of Survey Control - v2.2
Guideline For Adjustment and Evaluation of Survey Control - v2.2
Special Publication 1
Version 2.2
Document History
DATE VERS ISSUE AMENDMENTS AUTHOR(S)
7/12/2020 2 2 Update to copyright, ICSM Geodesy Working
update of term AGRS, Group
revision of references,
Section 2.1 and Section 6.1
update datum and geoid
model references
24/09/2014 2 1 Values in Table 10 updated, ICSM Permanent
copyright statement Committee on Geodesy
updated
24/10/2013 2 0 Document available ICSM Permanent
Committee on Geodesy
With the exception of the Commonwealth Coat of Arms, and where otherwise noted, this
product is provided under a Creative Commons Attribution 4.0 International Licence.
http://creativecommons.org/licenses/by/4.0/legalcode
Table of contents
List of figures
Figure 1: Combined GNSS and terrestrial control survey ........................................................ 9
Guideline for the Adjustment and Evaluation of Survey Control – SP1 iii
Version 2.2
Intergovernmental Committee on Surveying and Mapping
List of tables
Table 1: Station information .................................................................................................... 9
Table 2: Terrestrial measurement information ..................................................................... 10
Table 3: GNSS measurement information ............................................................................. 11
Table 4: AUSGeoid2020 interpolated values ......................................................................... 11
Table 5: Minimally constrained adjustment results .............................................................. 12
Table 6: Baseline components ............................................................................................... 12
Table 7: 95% SU values (minimally constrained GNSS-only adjustment).............................. 13
Table 8: Minimally constrained combined adjustment results ............................................. 13
Table 9: 95% SU values (minimally constrained combined adjustment) .............................. 14
Table 10: AGRS constraints (1 sigma) .................................................................................... 14
Table 11: Constrained adjustment results............................................................................. 15
Table 12: 95% PU values (constrained adjustment) .............................................................. 15
Term/Acronym Definition
Coverage factor (k) A statistical scalar used to modify a standard uncertainty to reach a
greater (or expanded) level of confidence. The value assigned to k
is chosen according to the desired level of confidence.
Critical value Used in the context of testing reliability, it is the value at the
extremity of a confidence interval or range. Estimated results
exceeding this value are deemed to have failed the particular test
of reliability.
Sigma zero The magnitude of the sum of the squares of the weighted
measurement corrections arising from a least squares adjustment.
1.1 Introduction
The availability of accurate and reliable information relating to the position and uncertainty
of Australia’s survey control marks is critical to the integrity of the Australian Geospatial
Reference System (AGRS). The purpose of this Guideline is to promote the adoption of
uniform least squares adjustment and evaluation procedures to achieve the highest level of
rigour and integrity in Australia’s survey control mark network.
The generation of position and uncertainty for Australia’s survey control marks relies upon
a large number of measurements to these survey control marks using a variety of
instruments, measurement techniques and processing software. Whilst a surveyor will
always endeavour to obtain the most accurate and precise measurements to survey
control marks, the true value of a survey control mark’s position can never be measured
nor derived with absolute certainty due to the inescapable presence of measurement
error. In order to determine the single, most reliable position from a range of
measurements whilst at the same time detecting and removing unacceptable
measurement errors, the technique of least squares estimation should be used.
For geodetic control surveys, the least squares technique provides a reliable means for
estimating the most probable value of survey control mark coordinates and uncertainties
from a redundant set of measurements by satisfying predefined functional and stochastic
models. The technique of least squares also provides a rigorous means for testing the
estimated values and the measurements used to derive them.
This Guideline outlines ICSM’s recommended procedures for using least squares to
estimate and evaluate the position and uncertainty of survey control marks, to evaluate
measurement uncertainty, and to evaluate network reliability. Examples are provided at
the end of this Guideline.
The following documents may have relevance to the application of this Guideline.
International Guidelines
SP1 Standard
ICSM (2020), Standard for the Australian Survey Control Network – Special Publication 1,
Version 2.2, Intergovernmental Committee on Surveying and Mapping, Canberra,
Australia.
SP1 Guidelines
ICSM (2020), Guideline for Control Surveys by Differential Levelling, Version 2.2,
Intergovernmental Committee on Surveying and Mapping, Canberra, Australia.
ICSM (2020), Guideline for Control Surveys by GNSS, Version 2.2, Intergovernmental
Committee on Surveying and Mapping, Canberra, Australia.
ICSM (2020), Guideline for Conventional Traverse Surveys, Version 2.2, Intergovernmental
Committee on Surveying and Mapping, Canberra, Australia.
2 Connection to datum
Survey control marks established for the AGRS shall be coordinated relative to the datums
set out in Section 2 of the Standard.
For Datum Control Surveys, propagation of datum and uncertainty to new and existing
survey control marks can be achieved by including the survey within AGRS adjustments
(State, Territory and Australian Government).
For General Purpose Control Surveys, the jurisdictional agency responsible for maintaining
the survey control mark network should make available, where possible, the relevant
survey control information to facilitate the rigorous propagation of datum and uncertainty.
Such information should include, at the least, datum aligned coordinates and heights,
uncertainties and metadata for the survey control marks to which connection is being
made.
Similarly, the jurisdictional agency responsible for maintaining the Australian Height Datum
(AHD) should make available the relevant survey control information to facilitate the
rigorous propagation of AHD. Such information should include AHD heights, quality
estimates and metadata for the survey control marks to which a connection to AHD is
being made.
The preferred means for propagating datum and uncertainty throughout an adjustment is
to provide coordinates and a variance-covariance matrix which rigorously expresses the
datum and quality (absolute and relative) of the existing survey control marks. However, in
practice a full variance-covariance matrix is often unavailable. In such instances, a priori
statistical information may be derived using empirical means or estimated values.
3 Expression of uncertainty
Section 4.2.1 of the Standard states that Survey Uncertainty (SU), Relative Uncertainty (RU)
and Positional Uncertainty (PU) shall be expressed in terms of the 95% confidence level
using any one of the following as appropriate:
standard deviations;
To express standard deviations for one dimensional components at the 95% confidence
level, the uncertainty value is simply computed by scaling the estimated (1 sigma) standard
deviation by coverage factor k = 1.960.
To express the (two dimensional) standard error ellipse at the 95% confidence level, the
axes of the 95% error ellipse are obtained by scaling the (1 sigma) axes by coverage factor k
= 2.448.
Similarly, the axes of the (three dimensional) 95% error ellipsoid are obtained by scaling
the (1 sigma) axes by coverage factor k = 2.796.
For the horizontal circular confidence region, the 95% uncertainty value is calculated from
the standard (1 sigma) error ellipse and is expressed as a single quantity, being the radius
of the circular confidence region. The radius (r) of the circular confidence region is
computed by:
𝑟 = 𝑎×𝐾
𝐾 = 𝑞0 + 𝑞1 𝐶 + 𝑞2 𝐶 2 + 𝑞3 𝐶 3
𝐶 = 𝑏⁄𝑎
Where:
a = semi-major axis of the standard error ellipse
b = semi-minor axis of the standard error ellipse
q0 = 1.960790
q1 = 0.004071
q2 = 0.114276
q3 = 0.371625
Values for a and b shall be derived from the full a-posteriori variance matrix obtained from
least squares adjustment.
4.1 Purpose
In the context of control surveys, the primary purpose of least squares adjustment is to
estimate survey control mark coordinates from a set of measurements. In addition to this,
least squares provides for the following additional outcomes:
the necessary statistics to evaluate the quality of the measurements and the
estimated coordinates.
(a) All survey control measurements should be corrected for all known calibration
corrections and systematic error sources, and be accompanied by reliable
values of uncertainty (or weights). To test the control survey for errors,
redundant measurements sufficient to identify errors shall be used. The larger
the degrees of freedom (DoF), the greater confidence can be gained from a
survey.
(b) A minimally constrained adjustment should be tested using the local test (see
Section 5.1) and global test (see Section 5.2).
(c) If required, estimated SU and associated RU values (or other reliable statistical
methods) should be examined to assess whether the survey has achieved any
predefined uncertainty or quality thresholds.
(d) When attempting to propagate datum and uncertainty, a fully constrained and
appropriately weighted adjustment should be undertaken.
(e) The fully constrained adjustment should be tested using the local test (see
Section 5.1), to verify that the imposed constraints do not result in
measurement failure(s).
(f) The fully constrained adjustment should be tested using the global test (see
Section 5.2). If this adjustment test fails, the quality of survey measurements
(g) If required, estimated PU and associated RU values (or other reliable statistical
methods) should be examined to assess whether the survey control network
has achieved any predefined uncertainty or quality thresholds.
Whilst ICSM regards the circular form of PU as an acceptable means for simplifying the
expression of a-posteriori uncertainty at the 95% confidence level, the circular form of
uncertainty is not acceptable for use as a-priori statistical information to be used in
constraining an adjustment in procedure (d).
Local tests should be conducted using the Normal distribution at the 95% confidence level.
A correction which exceeds the upper 95% confidence limit indicates a failure and the need
to re-evaluate the assumed uncertainty of the measurements and/or imposed constraints.
If the adjustment result is equal to the DoF, the sigma-zero (or variance factor) will be
unity, indicating that the system of survey measurements, uncertainties and constraints is
statistically reliable. Values larger than one indicate that one or more of the a-priori
measurement uncertainties are over-optimistic, or a larger-than-expected correction has
resulted. Values less than one suggest that the measurements were better than assumed
by the combined set of measurement uncertainties. Values which exceed the upper
confidence limit indicate a failure and the need to re-evaluate the uncertainty of the
measurements (via the local test) and/or the uncertainty of the imposed constraints.
Whilst it is not always possible, nor essential, to achieve sigma-zero values equal to unity, a
pass for both local and global tests is a minimum requirement for demonstrating survey
control network reliability.
ICSM recommends the following procedure to evaluate the coordinates and uncertainties
of marks in a survey control network using least squares:
(a) Confirm that the adopted procedures for each measurement technique are
commensurate with the desired quality outcome. Refer to the relevant survey
Guidelines for guidance on selecting survey equipment and procedures
appropriate for desired quality outcomes.
(b) Using the global test (see Section 5.2), confirm the sigma-zero value obtained
from the minimally constrained adjustment is within the upper and lower
limits of the 95% confidence interval.
(c) Confirm all measurements pass the local test at the 95% confidence level (see
Section 5.1) in the minimally constrained adjustment.
(d) When propagating datum and uncertainty from the AGRS or including a survey
within the AGRS, undertake a fully constrained adjustment using appropriate
constraints, and:
(i) Confirm the sigma-zero value obtained from the global test is within the
upper and lower limits of the 95% confidence interval
(ii) Using the local test, confirm the measurements still satisfy the expected
measurement uncertainties at the 95% confidence level
(e) As required, evaluate the standard deviations, error ellipses, PU and RU at the
95% confidence level for all survey control marks.
All GNSS baseline measurements and uncertainties are based upon a two-hour observation
period and have been extracted directly from baseline processing software. Each horizontal
angle has been reduced from a set of three angle measurements.
Table 1: Station information
Measurement Mark 1 Mark 2 Mark 3 Value Std. dev. Inst. ht. Targ. ht.
Level diff. 1 21 22 0.506 0.010 m
Level diff. 2 22 23 -0.045 0.011 m
Level diff. 3 23 25 -1.320 0.011 m
Level diff. 4 24 21 0.112 0.011 m
Level diff. 5 25 24 0.754 0.015 m
Level diff. 6 22 21 -0.506 0.010 m
Level diff. 7 23 22 0.042 0.011 m
Level diff. 8 25 23 1.318 0.011 m
Level diff. 9 21 24 -0.112 0.011 m
Level diff. 10 24 25 -0.740 0.015 m
Slope dist. 1 21 22 306.790 0.010 m 1.650 1.651
Slope dist. 2 22 23 410.031 0.010 m 1.650 1.651
Slope dist. 3 24 21 356.003 0.010 m 1.650 1.651
Slope dist. 4 25 23 371.568 0.010 m 1.650 1.651
Slope dist. 5 25 24 715.488 0.010 m 1.650 1.651
Vt. angle 1 21 22 0 05 35.651 2.0” 1.650 1.651
Vt. angle 2 22 23 -0 00 25.657 2.0” 1.650 1.651
Vt. angle 3 23 25 -0 12 20.092 2.0” 1.650 1.651
Vt. angle 4 24 21 0 01 09.908 2.0” 1.650 1.651
Vt. angle 5 25 24 0 03 20.601 2.0” 1.650 1.651
Hz. angle 1 24 21 25 91 18 43.522 1.0”
Hz. angle 2 25 24 23 88 55 20.712 1.0”
Hz. angle 3 23 25 22 89 47 58.951 1.0”
Hz. angle 4 22 23 21 180 04 31.199 1.0”
Hz. angle 5 21 22 24 89 53 25.614 1.0”
To cater for the influence of gravity on the angle and spirit levelling measurements,
ellipsoid-geoid separations (N) and deflections of the vertical in the prime meridian (ξ) and
prime vertical (η) have been interpolated from AUSGeoid2020. The respective
AUSGeoid2020 values for all marks are shown in Table 4.
Table 4: AUSGeoid2020 interpolated values
Station N ξ η
21 4.512 -2.934 -2.535
22 4.515 -2.950 -2.541
23 4.518 -2.974 -2.576
24 4.507 -2.931 -2.513
25 4.512 -2.963 -2.544
26 4.506 -2.943 -2.508
Unknown parameters 9
Measurements 18
Degrees of freedom 9
Sigma zero 1.380
Upper limit 2.114
Lower limit 0.300
Using the global test, the network adjustment is deemed a success since the sigma-zero
value (1.380) is less than the upper limit (2.114). Due to the high precision of the GNSS
baseline measurements, the SU for all coordinates is less than 1 millimetre.
Using the local test, each baseline is examined to assess whether the correction exceeds
the 95% critical value. A procedure for this test is to examine whether the normalised
residual (calculated by dividing the measurement correction by the standard deviation of
the correction) exceeds the critical value of the unit Normal distribution at 95%, which is
1.96. From this test, it is discovered that the X and Y components of baseline 1 fail (Table 6)
since the normalised residuals exceed 1.96.
Table 6: Baseline components
In this example, it is assumed that all GNSS observations were taken under the same set of
conditions and that these failures are most likely the consequence of over-optimistic
measurement precisions, rather than baseline 1 being a gross error. In this case, it is
appropriate to modify (i.e. scale) the a-priori uncertainties.
As a first attempt to rectify these failures, the sigma zero value (1.380) is used to scale all
GNSS measurement uncertainties. This is based on the assumption that all GNSS baselines
were derived under the same conditions and there are no gross errors. Following scaling of
all measurements, the minimally constrained adjustment produces a new sigma zero value
of unity. Repeating the adjustment yields a pass in the global test. However, the Y
component of baseline 1 still fails the local test.
From this analysis, a better way to deal with this problem would be to rescale the variance
matrix for baseline 1 by partial matrix scalars of 1.0, 1.0 and 5.0, and to rescale all other
baseline variance matrices back to their original values. Re-running the minimally
constrained adjustment leads to all GNSS measurements passing the local test, and a sigma
zero of 1.139 – which is a much more stable outcome and confirms that the precision for
baseline 1 was indeed over optimistic.
95% SU and circular radius values for all marks arising from the minimally constrained
GNSS-only adjustment are shown in Table 7.
Table 7: 95% SU values (minimally constrained GNSS-only adjustment)
Having confirmed all GNSS measurements pass the local test, the remaining terrestrial
measurements are introduced in order to verify the set of network measurements as a
whole. The ellipsoid-geoid separations (N) and deflections of the vertical in the prime
meridian (ξ) and prime vertical (η) shown in Table 4 are also introduced.
Unknown parameters 15
Measurements 43
Degrees of freedom 28
Sigma zero 0.778
Upper limit 1.588
Lower limit 0.547
This adjustment yields a pass in the global test and a pass in the local test for all GNSS and
terrestrial measurements. As indicated by the lower than expected sigma zero value, many
of the terrestrial measurements were better than anticipated. Whilst much time could be
spent modifying the a-priori standard deviations to achieve a more stable result, the
outcome is within the lower and upper limits and thereby deemed a success.
95% SU and circular radius values for all marks arising from the minimally constrained
combined adjustment are shown in Table 9.
Having confirmed all measurements pass the local test, constraints are imposed on the
adjustment to propagate datum and uncertainty from the AGRS. The adopted AGRS
constraints are survey control mark 23 and GNSS CORS site (mark 26). The respective
GDA2020 coordinates, AHD height and uncertainties for these marks are shown in Table
10. All uncertainties are at the 1 sigma (68.3% confidence level).
The published GDA2020 coordinates and standard deviations for survey control mark 23
have been estimated from a jurisdictional (state-wide) adjustment. The standard
deviations in seconds translate to approximately 0.025 m in the north-south and east-west
axes. Since a rigorously derived ellipsoid height value for mark 23 is unavailable, the AHD
height is introduced as an approximation. A 0.2 m standard deviation has been adopted for
AHD height for mark 23.
Table 10: AGRS constraints (1 sigma)
The GDA2020 (Cartesian) coordinates and uncertainties for the GNSS CORS site have been
derived from routine analysis of Australia’s GNSS CORS sites. Hence, a rigorous, fully
populated variance matrix is available. As shown by the lower VCV values in Table 10, the
standard deviations are in the order of 1 mm. Since this is a rather optimistic estimate, in
this example the VCV has been scaled by 7.5 to yield uncertainties at the 3 mm level.
For the constrained adjustment, all marks are held free and the constraints are introduced
as measurements with their assumed uncertainties. The introduction of these constraints
as measurements translates to a total measurement count of 49. Setting mark 22 as free
leads to 18 unknowns, resulting in a DoF of 31. The constrained adjustment at the 95%
confidence level yields the results in Table 11.
Table 11: Constrained adjustment results
Unknown parameters 18
Measurements 49
Degrees of freedom 31
Sigma zero 0.735
Upper limit 1.556
Lower limit 0.566
This adjustment yields a pass in the global test and a pass in the local test for all
measurements and constraints. As inferred from the lower than expected sigma zero value,
the system of measurements were better than indicated by the prescribed uncertainties,
and none of the constraints was shown to bias the adjustment in a significant way or to
cause any of the measurements to fail.
95% PU and circular radius values for all marks arising from the constrained adjustment are
shown in Table 12. From the foregoing analysis, the adjustment has proven to be a success.
The benefits of including the survey measurements within the jurisdiction adjustment
cannot be over emphasised. For instance, consider AGRS survey control mark 23. Prior to
this adjustment, the best published GDA2020 coordinates for mark 23 were stated as being
precise to 0.025 m, 68.3% confidence level (or 0.049 m, 95%). However, according to Table
12 the estimated uncertainty for mark 23 based on the new system of measurements is
0.003 m, 68.3% (or 0.007 m 95%). Unless this survey is integrated within the AGRS,
subsequent users of mark 23 will not be able to benefit from this increased precision.