Skip to main content

Derive three point analytical resection

Short answer

From known coordinates of A,B,C compute AB=c, BC=a and ∠ABC=B. With measured α=∠APB, β=∠BPC, let x=∠PAB, y=∠PCB. In quadrilateral sum, x+y = 360°-(α+β+B)=S. Using sine rules in triangles APB and BPC, equating PB gives sin x/sin y = (a sinα)/(c sinβ)=k. Apply identity to get tan((x-y)/2) = ((k-1)/(k+1)) tan(S/2), whence x and y. Then compute AP or BP length and bearing to obtain P coordinates.

Derive the expression for the three point resection problem.

In the three-point resection problem the coordinates of three stations A, B, C are known and the angles α=∠APB\alpha = \angle APB and β=∠BPC\beta = \angle BPC are measured at the unknown station P.

         B
        /|\
     c / | \ a
      /  |  \
     A   |   C
      \  |  /
       \ | /
        \|/ 
         P      angles at P: alpha (APB), beta (BPC)

Derivation

From the coordinates, find the sides c=ABc = AB, a=BCa = BC and the angle B=∠ABCB = \angle ABC (from the bearings of BA and BC). Let x=∠PABx = \angle PAB and y=∠PCBy = \angle PCB be the unknown angles. In the quadrilateral PABC the sum of angles is 360∘360^\circ:

x+y+α+β+B=360∘⇒x+y=S=360∘−(α+β+B)x + y + \alpha + \beta + B = 360^\circ \quad\Rightarrow\quad x + y = S = 360^\circ - (\alpha + \beta + B)

In triangle APB, by the sine rule, PB=csin⁡xsin⁡αPB = \dfrac{c\sin x}{\sin\alpha}. In triangle BPC, PB=asin⁡ysin⁡βPB = \dfrac{a\sin y}{\sin\beta}. Equating,

sin⁡xsin⁡y=asin⁡αcsin⁡β=k\frac{\sin x}{\sin y} = \frac{a\sin\alpha}{c\sin\beta} = k

Using sin⁡x−sin⁡ysin⁡x+sin⁡y=tan⁡x−y2tan⁡x+y2\dfrac{\sin x - \sin y}{\sin x + \sin y} = \dfrac{\tan\frac{x-y}{2}}{\tan\frac{x+y}{2}}:

tan⁡x−y2=k−1k+1tan⁡S2\tan\frac{x-y}{2} = \frac{k-1}{k+1}\tan\frac{S}{2}

So x−y2\dfrac{x-y}{2} is found, and with x+y=Sx + y = S we get xx and yy separately.

Coordinates of P

  1. In triangle APB: ∠ABP=180∘−α−x\angle ABP = 180^\circ - \alpha - x, so AP=csin⁡∠ABPsin⁡αAP = \dfrac{c\sin\angle ABP}{\sin\alpha} (or BP=csin⁡xsin⁡αBP = \dfrac{c\sin x}{\sin\alpha}).
  2. Bearing of APAP = bearing of ABAB ± x\pm\, x (sign depends on the side on which P lies).
  3. EP=EA+APsin⁡(bearing AP)E_P = E_A + AP\sin(\text{bearing }AP),   NP=NA+APcos⁡(bearing AP)\;N_P = N_A + AP\cos(\text{bearing }AP).
  4. Check by computing P from C: CP=asin⁡∠CBPsin⁡βCP = \dfrac{a\sin\angle CBP}{\sin\beta} with the bearing of CPCP = bearing of CB∓yCB \mp y. Both must give the same coordinates.

Note: If P lies on the circle through A, B and C (the "danger circle"), then α+β+B=180∘\alpha + \beta + B = 180^\circ, tan⁡S2\tan\frac{S}{2} is infinite, and the problem has no unique solution. The point P should be chosen away from this circle.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Surveying II chapter 5 questions →