Tienstra Resection

Planimetric Positioning via the Tienstra Three-Point Resection

In engineering geodesy, establishing high-precision horizontal control is a prerequisite for total station setups and subsequent stakeout or monitoring. GNSS is often blocked in urban canyons, open pits, and forest canopy, forcing crews to rely on classical optical methods. When distances to control cannot be measured — due to hazards, obstacles, or coordinate-only networks — the surveyor must solve position from angles alone.

The three-point angular resection establishes a free station at an unknown point $P$ by measuring horizontal angles to three visible, coordinated benchmarks. The Tienstra resection calculator implements the cotangent barycentric formulation for rapid coordinate resolution from field inputs. This article covers the geometry, danger circle, field practice, computation, and verification.

Defining the Resection Problem: Resection vs. Intersection

Planimetric positioning divides into two primary configurations:

Intersection

The unknown point is an inaccessible target (peak, spire, remote monument). Instruments occupy at least two known stations; bearings from each station intersect at the target. Accurate, but requires multiple setups and access to control.

Resection

The surveyor occupies the unknown point $P$ and observes horizontal angles between three known stations $A$, $B$, and $C$. Only one setup is required; control points need not be occupied — only visible. Ideal for construction site free-stationing referencing distant stable benchmarks.

The Tienstra Three-Point Method

Classical geometric methods (Cassini, Collins) use auxiliary constructions and are error-prone in manual calculation. The Tienstra formula — popularized by J.M. Tienstra (1895–1951) at Delft — solves $P$ directly with barycentric weights from cotangents of control-triangle interior angles and the observed angles at $P$.

Once weights $f_1$, $f_2$, $f_3$ are found:

$E_P = \frac{f_1 E_A + f_2 E_B + f_3 E_C}{f_1 + f_2 + f_3}$

$N_P = \frac{f_1 N_A + f_2 N_B + f_3 N_C}{f_1 + f_2 + f_3}$

Point $P$ is the weighted balance of triangle $ABC$ — the formula handles $P$ inside or outside the reference triangle.

Control Geometry and the Danger Circle

Resection precision depends on the spatial arrangement of control and the unknown station. Poor geometry magnifies angular errors or causes failure.

Point ordering and angle closure

Sequence control points $A$, $B$, $C$ counterclockwise around triangle $ABC$ (right-handed grid: easting = $x$, northing = $y$).

At unknown station $P$, measure horizon angles clockwise around the full $360°$:

• $\alpha$ = $\angle APB$ (from $A$ to $B$)
• $\beta$ = $\angle BPC$ (from $B$ to $C$)
• $\gamma$ = $\angle CPA$ (from $C$ to $A$)

Closure requirement:

$\alpha + \beta + \gamma = 360^\circ$

The danger circle

The critical singularity is the danger circle (circumcircle of $\triangle ABC$). If $P$ lies on that circle, infinitely many positions produce the same observed angles (inscribed angle theorem).

Algebraically, barycentric denominators approach zero:

$\cot A - \cot\beta \to 0 \qquad \cot B - \cot\gamma \to 0 \qquad \cot C - \cot\alpha \to 0$

Ensuring geometric strength

• Enclosure: $P$ inside $\triangle ABC$ is most stable (positive weights, robust to noise).
• Outside triangle: avoid placing $P$ on or near the circumcircle. A quick check: if interior angle at $B$ plus observed $\gamma$ (opposite $B$) $\approx 180°$, you are near the danger circle.
• Collinearity: three collinear control points give degenerate geometry — extremely unstable.

              A -------- B
             /    *P?    \
            /   (danger    \
           C --- circle) --- 

Field Observation Practice

Face I and Face II

Eliminate collimation, tilt-axis, and eccentricity errors by observing each target in Face I (direct) and Face II (reverse) and averaging.

Horizon closure

After sighting $A \to B \to C \to A$, compare the angle sum to $360°$:

$\delta_\theta = 360^\circ - (\alpha + \beta + \gamma)$

If within instrument tolerance, distribute evenly:

$\alpha_{\text{corr}} = \alpha + \frac{\delta_\theta}{3} \qquad \text{(and similarly for } \beta, \gamma \text{)}$

Typical tolerance: $\pm c\sqrt{n}$ where $c$ is instrument accuracy (seconds) and $n = 3$ angles.

Height

Standard Tienstra resection is planimetric (2D). Record instrument and target heights for pairing with later vertical observations.

Computation Walkthrough

Step 1: Interior angles of the control triangle

Compute clockwise grid bearings between control points (use the Bearing & distance calculator for inverse checks). For CCW-ordered $A$, $B$, $C$, interior angle at $A$:

$A = \theta_{AB} - \theta_{AC}$

(normalize to $0°$–$360°$). Similarly find $B$ and $C$. Verify $A + B + C = 180°$.

Step 2: Barycentric weights

With observed angles at $P$ ($\alpha$ opposite $C$, $\beta$ opposite $A$, $\gamma$ opposite $B$):

$f_1 = \frac{1}{\cot A - \cot \beta}$

$f_2 = \frac{1}{\cot B - \cot \gamma}$

$f_3 = \frac{1}{\cot C - \cot \alpha}$

Step 3: Coordinates of $P$

$E_P = \frac{f_1 E_A + f_2 E_B + f_3 E_C}{f_1 + f_2 + f_3} \qquad N_P = \frac{f_1 N_A + f_2 N_B + f_3 N_C}{f_1 + f_2 + f_3}$

Worked Example

Given data

Observed angles at $P$ (face-averaged, horizon closed):

• $\alpha = 118.4550°$ ($\angle APB$, opposite $C$)
• $\beta = 104.3820°$ ($\angle BPC$, opposite $A$)
• $\gamma = 137.1630°$ ($\angle CPA$, opposite $B$)

Check: $118.4550 + 104.3820 + 137.1630 = 360.0000°$.

Interior angles of $\triangle ABC$

• $\theta_{AB} = 90.0000°$, $\theta_{AC} = 33.6901°$ → $A = 56.3099°$
• $\theta_{BC} = 326.3099°$, $\theta_{BA} = 270.0000°$ → $B = 56.3099°$
• $\theta_{CA} = 213.6901°$, $\theta_{CB} = 146.3099°$ → $C = 67.3802°$

Sum: $180.0000°$.

Cotangents and weights

$\cot A \approx 0.666667 \quad \cot B \approx 0.666667 \quad \cot C \approx 0.416667$

$\cot\alpha \approx -0.541951 \quad \cot\beta \approx -0.256428 \quad \cot\gamma \approx -1.078441$

$f_1 = \frac{1}{0.666667 - (-0.256428)} \approx 1.083 \qquad f_2 \approx 0.573 \qquad f_3 \approx 1.043$

$\sum f \approx 2.699$

Result

$E_P \approx 1810.962 \text{ m} \qquad N_P \approx 2579.650 \text{ m}$

Enter the same values in the Tienstra resection calculator to verify.

Checking Your Result

Back-bearing check

From computed $P$, calculate grid bearings to each control point using $\text{atan2}(\Delta E, \Delta N)$, then reconstruct $\alpha$, $\beta$, $\gamma$. Reconstructed angles should match field values within a few seconds of arc.

Alternative methods

• Collins point — auxiliary construction on the circumcircle
• Cassini — intersecting auxiliary circles through $A,B,P$ and $B,C,P$
• Least-squares resection — when distances and redundant observations are available

Limitations and Practical Warnings

Sensitivity to angular errors

Cotangents blow up near $0°$ and $180°$. Small pointing errors can shift $P$ by metres when angles are ill-conditioned.

Collinear control

Nearly collinear $A$, $B$, $C$ collapses triangle area and makes side-of-line ambiguity severe.

No redundancy

Three angles solve two unknowns exactly — no redundancy to catch misidentified targets or disturbed control. Observe a fourth point for least-squares adjustment when possible.

Next Steps

Tienstra resection establishes 2D control from angles alone — elegant and field-practical when geometry is strong and the danger circle is avoided. For rapid computation without manual cotangent steps, use the Tienstra resection calculator.