101-Article Text-715-1-10-20211115
101-Article Text-715-1-10-20211115
101-Article Text-715-1-10-20211115
ABSTRACT
This paper presents a simple way of finding the accuracy of three-point resection problem. The main idea is to
distribute known points in four-quadrants and compute the coordinates of the resected point (P), then investigate
the affection of known points positions in one quadrant or combination quadrant for the accuracy of resected point
(which is a new idea in the literature). This article aims to find out the affections of distribution control points in
different quadrants on the accuracy of the resected point, and also choosing the best positions of the three-points
in terms of their positions in one quadrant or their positions in combination quadrants. The relative positions of the
known points and the resected point are playing the main role to decide the accuracy of the resected point. The
study of the relationship between relative positions of unknown point will be explained in detail, other known-
points and the accuracy of the unknown point will be introduced, and the positions of known points (Distributed in
quadrants) of having the best accuracy of the resected point will be defined.
Keywords: Three-Point Resection Problem, Surveying, Accuracy, Three-point solution, Resected point, Four-
quadrants.
1. INTRODUCTION
The three-point resection method (also known as the “three-point problem”) is used to compute the
coordinates of an unknown point by occupying the resected point and measuring angles and distance
toward three known points, since the angles between three known points are determined at the
unknown point (ElHassan, 1986). The idea of computing the coordinates of an unknown point using
relative angular measurements to three known stations is a basic concept in surveying engineering (El
Hassan, 1989) and it is a very useful technique for fixing positions quickly where it is best needed for
setting-out purposes which has been discussed and solved analytically or graphically in the last centuries
in most previous articles and books, such as (ElHassan, 1986; Font-Llagunes & Batlle, 2009; Haralic et al.,
1991; Ligas, 2013; Problem, n.d.; Wang et al., 2018). Moreover, (El Hassan, 1989) tried to detect the
accuracy of the three-point resection method by changing the position of resected and control points,
and this paper has reached to the relative location of the resected point, indicating that the three control
points have a big affection on deciding the accuracy of the resected point.
There are various methods to solve three-point resection, yet so far all of them failed in the control
point case, and the resected point lie on the circumference of a circle. A lot of methods produce
questionable results if the known points lie on a straight line, the method of computation should be
chosen with consideration, and independent checks should be used wherever possible. (Point &
Problem, 1981)
This paper precedes the study to find out an accurate coordinates which will be obtained using three-
point resection calculations in terms of the positions of three known points in one quadrant or
combination between two or three quadrants and comparing them with the real coordinate of the
resected points, and also to find out the appearance of any relationship between the unknown point and
the three-point position. Furthermore, it should be noted that all previous articles and books concerning
to this method did not address the question “Does the position of the three points in one quadrant
affect the accuracy of the coordinates of the unknown point or not?” and the answer will be presented
in this paper.
Received Sept 15; Revised 29 Oct; Accepted 9 Nov doi 10.12962/j27745449.v2i2.101
Corresponding author: Cahyadi ([email protected])
Osman et al. (16-32)
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
2. METHODOLOGY
2.1 The resection method cases
The space resection is described as an optimization model. The function of objective model is to
reduce the number of squared residuals (Easa, 2010). The model is given by:
r11 (XA −XL ) +r12 ( YA _ YL ) +r13 (ZA −ZL)
x=(−f) ( r )
31 (XA −XL ) + r32 (YA − YL )+ r33 (ZA − ZL )
r (X −X ) + r (Y − Y )+ r (Z −Z )
y=(−f) ( r21 (XA −XL )+ r 22 ( YA −Y L) + r23( ZA−Z L) )
31 A L 32 A L 33 A L
where: (XA , YA , ZA ) are the terrestrial coordinates, (x, y) are the image coordinates,
(XL , YL , ZL ) are the coordinates of the exposure station or the camera,
(f) is the focal length of the camera, and (rij) is the rotation matrix elements, if the value of camera's
focal length and the coordinates of known ground point, then the equation (x, y) becomes a function of
the three exposure station elements (XL , YL , ZL ) and the rotation angle elements (ω, ɸ, k), will obtain
linear equations by distributing the equations in a Taylor series. It can be solved by numerical methods
or least square method as well.
Two-point resection is a method of positioning new control point (resected point) by taking
observations and measurements towards the other known control points on a network. There are two
types of resection that can be performed, including angular resection by measuring horizontal angles and
distance resection by measuring horizontal distances. Both measurements are very useful for obtaining
on-site temporary control points which known as free station points. (Topic 6 : Angle Measurement :
Intersection and Resection, n.d.)
Angular resections are used to compute the coordinate point observing the existing control points –
the advantage of this method is that the resection can be performed without occupying any control
points during measurements. The resected point can be determined in angular resection by measuring
its angles from at least three existing control points in a three-point resection, which presented
below:(Topic 6 : Angle Measurement : Intersection and Resection, n.d.)
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
There are three cases in this method as shown in Figure (1), the three cases are considered to the
position of the known and unknown points. The position of the unknown point can be located in the
middle of a polygon or in one corner of the polygon. In this paper, the position of the unknown point
located in the corner of the polygon as presented in the following figure:
A(xa, ya)
ʆ
c b
β ɸ
B(xb, yb) C(xc, yc)
α β
P(x, y)
Figure (2) the positions of unknown point P and known points (Point & Problem, 1981)
where P is the resected point (unknown point), A, B, C are the control points (known points
coordinates), α, β are the angles between the known points measured from the resected point, ɸ, β, ʆ
are the angles that will be calculated mathematically using the following formulas, and b, c are the
distances between three known points which will be calculated mathematically as well.
Resection solution
There are a lot of previous studies that produced resection three-point problem solution which have
been conducted with different ways and methods. This paper is using the solution from (Point &
Problem, 1981) and it has been found that this solution gives high accuracy compared to the other
solutions. The first step is to calculate the angles since (A, B, C) are known points, the angle (ʆ) can be
calculated by the deflection differences of the two lines.
Calculating the assisting angles ɣ and ∅. Firstly, it needs to be noted that the sum of the interior
angles is equal to 360o [the summation of interior angles of any polygon should equal (n –
2)*180o].(Point & Problem, 1981)
α + ɣ + β + ʆ + ∅ = 360 (1)
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Rearrange:
ɣ = 360 − (α + β + ʆ ) − ∅ (2)
The interior angles α, ʆ, and β are known, therefore the calculation of assisting angles ɣ and ∅ is
needed:
ɣ=R−∅ (3)
Assuming that:
b sin(α) sin(ɣ)
K1 = c ∗ sin(β) = sin(∅) (4)
Compute ∅
K1 +cos(R)
cot(∅) = (9)
sin(R)
Compute ɣ
K2 +cos(R)
cot(ɣ) = (10)
sin(R)
Figure (1), was done by solving or taking the triangle (APC) and the angle ɸ on point C that will be
calculated using equation No (9). Then, by finding the line deflection (CP) from the line deflection (CA)
(known deflection), and finally the coordinates calculation of the resected point (P) was done using the
following equations: (Point & Problem, 1981)
E = L ∗ sin(𝛼) (11)
N = L ∗ cos(𝛼) (12)
Where ΔE is the difference of Easting between resected point (P) and the control point (C), and ΔN is
the difference of Northing between resected point (P) and the control point (C).
Or it can be:
EP=EB+L*sin(𝛼) (13)
NP=NB+L*cos(𝛼) (14)
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
The triangle (APB) in Figure 1 was used to ensure that the previous computations are correct and the
angle ∅ should be computed using the equation (9), to find the coordinates point (P) using similar
equations as before, and the final equations of computing resected point will be as follows:
E = L ∗ sin(∅) (15)
N = L ∗ cos(∅) (16)
Where L is the distance between resected point P control point (B), ΔE is the difference of Easting
between resected point (P) and the control point (B), ΔN is the difference of Northing between resected
point (P) and the control point (B).
EP=EB+L*sin(∅) (17)
NP=NB+L*cos(∅) (18)
Figure (3) distributions of control points around the resected point in four-quadrants
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
The figure (3) above shows that four known points have been distributed in each quadrant except
the control points located in the vertical and horizontal axis. The known points have been divided into
eight sectors and two groups, every group contains of 4 sectors. The first group represented by the
normal quadrants and the second group represented by the combinations between two quadrants. Their
classifications are presented as follows:
Group A: Represented by the First quadrant, second quadrant, third quadrant, and fourth
quadrant.
Group Two: First, it consists of the combination between the first and second quadrants which
contained the known points named as (E, 13, 14, 15, 16, N, 1, 2, 3, 4, W). Second, the
combination between the second and third quadrants which contains of known points named as
(N, 1, 2, 3, 4, W, 5, 6, 7, 8, S). Third, the combination between the first and second quadrants
which contains of known points named as (W, 5, 6, 7, 8, S, 9, 10, 11, 12, E). Fourth, the
combination between the first and second quadrants which contains the known points named as
(S, 9, 10, 11, 12, E, 13, 14, 15, 16, N).
There are two cases in the first quadrant that gives differences greater than 1 meter in the
northern coordinate compared to the real coordinates of the known point (P), which are: The
combination (the triangle) of the (S, 11, 12) known points with the unknown point (P), and by
taking the combination of (S, 12, E) known points with the unknown point (P).
There is one case in the second quadrant with differences greater than 1 meter in the northern
coordinate and more than 2 decimeters in the eastern coordinates which is the combination (the
triangle) between the (W, 5, 6) known points with the unknown point (P).
There are three cases for the known points located in the third quadrants: based on the points
combination (2, 3, 4) with the unknown point (P), it has differences greater than 2 decimeters in
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
the northern coordinate, and the case of taking the combination of the points (3, 4, W) with the
unknown point (P) has differences greater than 4 decimeters in the northern coordinate and
finally the combination between points (N, 1, 3) with the unknown point (P) has differences
greater than 5 decimeters in the eastern coordinate.
Finally, the case in the fourth quadrant combined from the points (E, 15, 13) with the resected
point (P), have differences greater than 1 meter in the northern and eastern coordinates.
Above exception cases have not calculated on the graph accuracy of the resection points.
Computed
No. Triangle No( ) Observed Np Observed Ep computed Ep Delta N Delta E
Np
1 <>P 10,9,S 1721100.920 440332.800 1721100.921 440332.801 0.001 0.001
2 <>P 11,10,S 1721100.920 440332.800 1721100.917 440332.798 -0.003 -0.002
3 <> P 12,11,S 1721100.920 440332.800 1721102.200 440332.770 1.280 -0.030
4 <> P E,12,S 1721100.920 440332.800 1721102.193 440332.765 1.273 -0.035
5 <> P 11,10,9 1721100.920 440332.800 1721100.917 440332.799 -0.003 -0.001
6 <> P 12,11,9 1721100.920 440332.800 1721100.850 440332.810 -0.070 0.010
7 <> P E,12,9 1721100.920 440332.800 1721100.973 440332.778 0.053 -0.022
8 <> P 12,11,10 1721100.920 440332.800 1721100.892 440332.770 -0.028 -0.030
9 <> P E ,12,10 1721100.920 440332.800 1721100.933 440332.860 0.013 0.060
10 <> P E,12,11 1721100.920 440332.800 1721100.973 440332.730 0.053 -0.070
Table (1): The real and computed coordinates and the Accuracy of the resected point (P) in the first
quadrant
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (5) shows the accuracy of the resected point (P) in the first quadrant
Table (1) shows that there are two cases in the first quadrant that gives differences greater than
1 meter in the northern coordinate compared with the real coordinates of the known point (P),
which are: The combination (the triangle) of the (S, 11, 12) known points with the unknown
point (P), and by taking the combination of (S, 12, E) known points with the unknown point (P).
Those two cases have not been counted on figure (5), some of the differences in the first
quadrant have millimeters and most of them reached centimeters as presented in Table (1) and
figure (5).
2- Computing the accuracy of resected point (P) for the second quadrant:
Table (2): The real and computed coordinates and the Accuracy of the resected point (P) in the second
quadrant
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (6) shows the accuracy of the resected point (P) in the second quadrant
Table (2) shows that there is one case in the second quadrant with differences greater than 1
meter in northern coordinate and more than 2 decimeters in eastern coordinate, which is the
combination (the triangle) between (W, 5, 6) known points with the unknown point (P). This case
has not been counted on figure (6). Some of the differences in the first quadrant are in
centimeters and most of them reached millimeters as presented in Table (2) and figure (6).
3- Computing the accuracy of resected point (P) for the third quadrant:
Table (3): The real and computed coordinates and the Accuracy of the resected point (P) in the third
quadrant
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (7) shows the accuracy of the resected point (P) in the third quadrant
Table (3) shows that there are three cases for the known points located in the third quadrants:
based on the combination points (2, 3, 4) with the unknown point (P) shows differences greater
than 2 decimeters in the northern coordinate, and the case of taking the combination of the
points (3, 4, W) with the unknown point (P) shows differences greater than 4 decimeters in the
northern coordinate and finally the combination between points (N, 1, 3) with the unknown
point (P) gives differences greater than 5 decimeters in the eastern coordinate, those three cases
have not been counted on figure (7).
4- Computing the accuracy of resected point (P) for the fourth quadrant:
Computed
No. Triangle No( ) Observed Np Observed Ep computed Ep Delta N Delta E
Np
1 <> P N,16,15 1721100.920 440332.800 1721100.917 440332.803 -0.003 0.003
2 <> P 16,15,14 1721100.920 440332.800 1721100.918 440332.806 -0.002 0.006
3 <> P 15,14,13 1721100.920 440332.800 1721100.931 440332.784 0.011 -0.016
4 <> P 14,13,E 1721100.920 440332.800 1721100.917 440332.803 -0.003 0.003
5 <>P 16,15,13 1721100.920 440332.800 1721100.920 440332.801 0.000 0.001
6 <>P N,15,14 1721100.920 440332.800 1721100.910 440332.311 -0.010 -0.489
7 <>P 16,13,E 1721100.920 440332.800 1721100.886 440332.941 -0.034 0.141
8 <>P 15,13,E 1721100.920 440332.800 1721100.973 440332.732 0.053 -0.068
9 <>P N,13,E 1721100.920 440332.800 1721100.920 440332.810 0.000 0.010
10 <>p N,14,E 1721100.920 440332.800 1721100.920 440332.810 0.000 0.010
Table (4): The real and computed coordinates and the Accuracy of the resected point (P) in the fourth
quadrant
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (8) shows the accuracy of the resected point (P) in the fourth quadrant
From the Table (2) showed the case in the fourth quadrant combined from the points (E, 15, 13) with the
resected point (P), gives differences greater than 1 meter in the northern and eastern coordinates, this
case has not counted on figure (6).
B. Combination quadrants:
The Combination quadrants contains from the combination between first and second,
second and third, third and fourth, and finally fourth and first. The accuracy of all those cases are
in millimeters which means that the accuracy is higher compared to the points located in normal
quadrants as presented in the tables, graphs and explanations below:
1- Computing the accuracy of resected point (P) for the combination of the first and second
quadrants:
Table (5): The real and computed coordinates and the Accuracy of the resected point (P) in the
combination between first and second quadrants
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (9) shows the accuracy of the resected point (P) in the combination between first and second
quadrants
Table (5) and figure (9) show the case in the combination of the first and second quadrants, all of the
cases has millimeters accuracy except the triangle of (1, N,16) which has one centimeter on the delta
eastern coordinates (ΔE).
2- Computing the accuracy of resected point (P) for the combination between second and third
quadrants:
Table (6): The real and computed coordinates and the Accuracy of the resected point (P) in the
combination between second and third quadrants
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (10) shows the accuracy of the resected point (P) in the combination between second
and third quadrants
Table (6) and figure (10) showed that the case in the combination of the first and second quadrants
which is all the cases has millimeters in their accuracy, indicating that its more accurate compared to the
previous graphs.
3- Computing the accuracy of resected point (P) for the combination between third and fourth
quadrants:
Table (7): The real and computed coordinates and the Accuracy of the resected point (P) in the
combination between third and fourth quadrants
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (11) shows the accuracy of the resected point (P) in the combination between third and fourth
quadrants
Table (7) and figure (11) showed the combination case of the first and second quadrants. All of the cases
has millimeters in their accuracy, which means that its more accurate compared to the previous graphs.
4- Computing the accuracy of resected point (P) for the combination between fourth and first
quadrants:
Table (8): The real and computed coordinates and the Accuracy of the resected point (P) in the
combination between fourth and first quadrants
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
Figure (12) shows the accuracy of the resected point (P) in the combination between fourth
and first quadrants
Table (8) and figure (12) showed the case in the combination of the first and second quadrants which is
all the cases has millimeters accuracy except the triangle of (5,3,2) which it has one centimeter on the
delta eastern coordinates (ΔE) and the triangle of (7,3,2) which it has one centimeter on the delta
northern coordinates (ΔN) as well.
5. Conclusion
The procedures of computing the coordinates of the resected point (P) from several known points
(control points), and computing the differences between the real coordinates of the resected point and
the coordinates of the three-point resection solution from several quadrants (the accuracy) have been
conducted. Thus, this study concludes that there is a relationship between the positions of control points
and the accuracy of resected point (P). Concerning to the accuracy tables which have been presented
before, it shows that the positions of the known-points in combination quadrants have high accuracy
compared to their positions in only one quadrant in the three-point resection solution which has been
used in this study.
doi 10.12962/j27745449.v2i2.101
Journal of Marine-Earth Science and Technology
Vol.2 Issue 2 (2021) ISSN 2774-5449
References:
Easa, S. M. (2010). Space resection in photogrammetry using collinearity condition without linearisation.
Survey Review, 42(315), 40–49. https://doi.org/10.1179/003962609X451681
El Hassan, I. M. (1989). Analysis of the three-point resection accuracy. Australian Surveyor, 34(7), 716–
727. https://doi.org/10.1080/00050326.1989.10438605
ElHassan, I. M. (1986). An Analytical Solution of the Resection Problem. Journal of Surveying
Engineering, 112(1), 30–35. https://doi.org/10.1061/(asce)0733-9453(1986)112:1(30)
Font-Llagunes, J. M., & Batlle, J. A. (2009). New Method That Solves the Three-Point Resection
Problem Using Straight Lines Intersection. Journal of Surveying Engineering, 135(2), 39–45.
https://doi.org/10.1061/(asce)0733-9453(2009)135:2(39)
Haralic, R. M., Lee, C. nan, Ottenberg, K., & Nolle, M. (1991). Analysis and solutions of the three point
perspective pose estimation problem. July, 592–598. https://doi.org/10.1109/cvpr.1991.139759
Ligas, M. (2013). Simple Solution to the Three Point Resection Problem. Journal of Surveying
Engineering, 139(3), 120–125. https://doi.org/10.1061/(asce)su.1943-5428.0000104
Point, T., & Problem, R. (1981). Three point resection problem. Compute.
Problem, P. (n.d.). An Efficient Algebraic Solution to the. 55455, 7225–7233.
Topic 6 : Angle measurement : Intersection and resection. (n.d.).
Wang, P., Xu, G., Wang, Z., & Cheng, Y. (2018). An efficient solution to the perspective-three-point pose
problem. Computer Vision and Image Understanding, 166, 81–87.
https://doi.org/10.1016/j.cviu.2017.10.005
doi 10.12962/j27745449.v2i2.101