18-Geom-A2 Adjustment of Observations · December 2018
Question 7 of 7: Linearized Observation Equation for an Angle
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams — December 2018 — 04-Geom-A2, Adjustment of Observations and Data Analysis. Three-hour, CLOSED-BOOK exam; an approved Casio or Sharp calculator is permitted. Format: seven questions are given and any five (20 marks each) constitute a complete paper — all seven are solved below for completeness. Adjustments are carried out in the Canadian survey frame (heights on CGVD2013, horizontal work on NAD83(CSRS)); US-foot units are retained in Q6 wherever the printed question uses them.
Reference texts: Ghilani, Adjustment Computations: Spatial Data Analysis (6th ed., Wiley, 2017); Wolf & Ghilani, Elementary Surveying: An Introduction to Geomatics (15th ed., Pearson, 2015); Mikhail & Gracie, Analysis and Adjustment of Survey Measurements (Van Nostrand, 1981).
Method note. Every weighted least-squares result below is formed from the normal equations $N\hat{x}=A^{\mathsf T}W\ell$; the reference (unit-weight) standard deviation is $s_0=\sqrt{v^{\mathsf T}Wv/(n-u)}$ with $n$ observations and $u$ unknowns, and station precisions come from $\Sigma_{xx}=s_0^2 N^{-1}$.
Question 7: Linearized Observation Equation for an Angle (20 marks)
Given. Observed clockwise angle $\angle GAB=107^\circ29'40''$ at instrument station A, from backsight G to foresight B, with approximate coordinates $G(578.741,1103.826)$, $A(415.273,929.868)$, $B(507.934,764.652)$.
Find. The angle observation equation linearized about these coordinates — its coefficients (in $''$/m) on the six coordinate corrections and its constant term.
Figure 7.1 — Angle $\angle GAB$ at instrument A, opening clockwise from backsight G to foresight B; the observation equation is the difference of the two azimuth equations A→B and A→G.
Approach. An angle is the difference of two azimuths, $\angle GAB=\text{Az}_{AB}-\text{Az}_{AG}$; differentiate each azimuth with respect to the station coordinates (the standard azimuth partials) and combine, giving coefficients in radians per metre that are scaled by $\rho''=206{,}265''$ to arc-seconds.
Approximate azimuths and distances. $\text{Az}_{AG}=\tan^{-1}\!\big(\tfrac{\Delta X}{\Delta Y}\big)=43^\circ13'10''$ ($D_{AG}=238.711$ m); $\text{Az}_{AB}=150^\circ42'51''$ ($D_{AB}=189.426$ m).
Computed angle and misclosure.
$$\angle GAB_0=\text{Az}_{AB}-\text{Az}_{AG}=107^\circ29'41.5'',$$
$$k=\angle_{\text{obs}}-\angle_{\text{computed}}=107^\circ29'40''-107^\circ29'41.5''=\boxed{-1.5''}.$$
Azimuth partials. For a line I→J, $\partial\text{Az}/\partial X_J=\Delta Y_{IJ}/D_{IJ}^2$, $\partial\text{Az}/\partial Y_J=-\Delta X_{IJ}/D_{IJ}^2$ (and equal-and-opposite at I). Applying these to $\text{Az}_{AB}$ (foresight, $+$) and $\text{Az}_{AG}$ (backsight, $-$) and multiplying by $\rho''$ gives the coefficients (arc-seconds per metre of coordinate shift):
$$\begin{aligned}
&\text{G: } dX_G=-629.7,\ dY_G=+591.7\\
&\text{A: } dX_A=+1579.4,\ dY_A=-59.1\\
&\text{B: } dX_B=-949.7,\ dY_B=-532.6
\end{aligned}$$
(as a check, the three $dX$ coefficients sum to zero, as do the three $dY$ — the angle is invariant to a rigid shift of the whole figure).
Assemble the observation equation. Writing $v$ for the angular residual ($''$),
$$\boxed{\begin{aligned}
-629.7\,dX_G+591.7\,dY_G+1579.4\,dX_A-59.1\,dY_A\\
-949.7\,dX_B-532.6\,dY_B=-1.5''+v
\end{aligned}}$$
which is one row of the design matrix for this angle (coordinate corrections in metres, coefficients in $''$/m).