KR101090346B1 - A method for transformating the parcel boarderland coordinates into world geodetic system coordinates - Google Patents

Abstract

에 의해 계산으로 산출하여 각 도근점에 대응하는 세계측지계 좌표를 취득하는 단계(단계 S3) 및; 상기 단계 S3에서 각 도근점(29001∼29007)에 대응하는 세계측지계 좌표가 계산에 의해 취득되면, 지역측지계 좌표를 산출할 때 현장에서 관측한 각 도근점에서 각 필지 경계점(1~37)까지의 관측각과 거리를 경계점 관측부로부터 확인한 다음, 경계점 관측부로부터 확인된 각 도근점에서 각 필지 경계점(1~37)까지의 관측각과 거리를 이용하여 각 필지 경계점(1~37)의 세계측지계 좌표를, 상기 수학식 1 및 상기 수학식 2에 의해 계산으로 산출하여 상기 각 필지 경계점(1∼37)에 대응하는 세계측지계 좌표를 취득하는 단계(단계 S4);를 갖추어 이루어진 것을 특징으로 하는 필지 경계점 좌표를 세계측지계 좌표로 변환하기 위한 방법을 특징으로 한다.According to the present invention, a cadastral survey reference point mark is installed on the ground to determine whether there are four known points (two points of departure (98210,98211) and two points of arrival (24263,24262) of the departure side) to which the GPS surveying device is to be installed. Determining that there is a case and determining that there is no case (step S1); If it is determined in step S1 that the four known points exist, the GPS measurement apparatus is mounted on the four known points, and then GPS observation is performed on the four known points and the four points of the four points are based on the GPS observations. Acquiring world geodetic coordinates corresponding to the known point (step S2); When the global geodetic coordinates are acquired for four known points in step S2, when the existing geodetic coordinates are calculated for each approach point (29001 to 29007), the observation angle and distance observed in the field are confirmed by the Cadastral Surveying Department. The coordinates of the world geodetic system corresponding to each of the approach points 29001 to 29007 using the observation angle θ and the distance D to each of the approach points 29001 to 29007 are represented by the equations 1 Δx = Xb-Xa, Δy = Yb-Ya (where Δx is the longitudinal difference, Δy is the horizontal difference, Xa, Xb, Ya, Yb are the coordinates of the X and Y axes), and

Calculating by calculation to obtain world geodesic coordinates corresponding to each approach point (step S3); If the world geodetic coordinates corresponding to each of the approach points 29001 to 29007 are obtained by calculation in step S3, the angles of observation from each approach point observed in the field to each parcel boundary point 1 to 37 when calculating the geodetic coordinates are calculated. After confirming the distance from the boundary point observer, the coordinates of the world geodetic coordinates of each parcel boundary point (1 to 37) using the observation angle and distance from each approach point identified from the boundary point observer to each parcel boundary point (1 to 37), And calculating the geodetic coordinates corresponding to the respective parcel boundary points 1 to 37 (step S4), calculated by Equation 1 and Equation (2). A method for converting coordinates is provided.

Description

A METHOD FOR TRANSFORMATING THE PARCEL BOARDERLAND COORDINATES INTO WORLD GEODETIC SYSTEM COORDINATES}

The present invention is a parcel for converting boundary coordinates based on a geodetic system (Bessel ellipsoid) to a world geodetic system (GRS80 ellipsoid) coordinates in an area having a boundary point coordinate register (intellectual study where the boundary of land ownership is registered as a planar vertical horizontal line value). The present invention relates to a method for converting boundary point coordinates into world geodetic coordinates.

The land boundary of Korea is registered and managed in cadastral map based on geodetic geodetic system, but it is inconvenient to use new geodetic geodetic technology. Therefore, it is transformed into world geodetic coordinates based on the 3D coordinate axis of the center of mass. In the case of cadastral drawings, the physical ground boundary does not change, but if there is a boundary change in the drawing due to coordinate transformation, there is a possibility that national confusion and boundary disputes may occur.

Also, coordinate transformations suggested by public survey and technical guidelines for global geodetic transformations have different parcel boundary coordinates depending on the common point position and quantity, and when restoring the transformation results on the ground, there is no survey reference point, so it is possible to perform boundary restoration survey. There is no problem.

Here, the coordinate transformation means that the absolute position of the point does not change with respect to one point P in space, but the coordinates in one coordinate system are changed to coordinates in another coordinate system by changing the coordinate system.

In order to transform the coordinates, an existing variable is used as an independent variable, and the coordinates in the target coordinate system are created as dependent variables. Finding the coefficient of the equation is the key in the coordinate transformation. That is, if the coordinate in the existing X coordinate system of the point P is (x, y), and you want to find the coordinate (x ', y') in the X 'coordinate system of the same point P, x' = ax + by, y An equation such as' = ay-by must be established and the coefficients a and b must be obtained.

After all, coordinate transformation in two-dimensional space refers to calculating coordinates in a new (x ', y') coordinate system of the same point from coordinates in a (x, y) coordinate system of one point.

Elements of two-dimensional coordinate transformation include origin movement, rotation, and scale change. FIG. 1 is a diagram for explaining two-dimensional coordinate transformation. FIG. 1A is a reference shape, and FIG. 1B is an X-axis based on the shape shown in FIG. Fig. 1C shows a state in which rotation has occurred, and Fig. 1D shows a state in which scale change (enlargement) has occurred.

When converting to global geodetic coordinates, there is no change in position in cadastral survey control points embedded in the earth's surface or parcel boundary points established in the past. However, in order to improve the flow of international surveying technology and the accuracy of cadastral surveying, the geodetic system is applied differently to register only the numbers registered in the study.

As a result, the most important point is whether the change in positional coordinates represents the actual scope of land ownership as it was, and therefore, cadastral surveys lose credibility if the boundary of the landscape changes, even for legal and administrative reasons. . In other words, the ability of cadastral surveyors to establish boundary piles on the ground is based on the absolute credibility of the landowners, and the boundary cannot be changed for any reason.

To calculate the coordinates of the new (x ', y') coordinate system of the same point from the coordinates of the (x, y) coordinate system of one point, at least four points of common point (x, y, x ', y') The transformation coefficient of coordinates is required, the common point should not be biased in a specific part, and the deviation of the statistical confidence interval between two coordinates is excluded, and then converted into world coordinates using the distortion grid model (hereinafter, Map conversion using "conversion method using conversion software").

The coordinate conversion method using the conversion software currently used cannot change the boundary registered in the cadastral study on the ground because the coordinates of the surveying reference point, which is the basis of the survey, are changed, and the land boundary changes according to the number of common points. Can not be used for cadastral surveys.

Therefore, even if the geodetic system is applied differently and registered by changing the number registered in the cadastral study, the absolute position on the surface should be identical. In other words, even if surveying is done in a geodetic or global geodetic system, the position of the point on the ground must be the same.

The present invention has been invented in view of the above-mentioned point, and in the area having the boundary point coordinate registration unit, the world geodesic coordinates are acquired by GPS observation for the known point (starting and arriving side) of the conductor, and the approach point and boundary point of the lead Provides a method for converting parcel boundary point coordinates, which are designed to convert boundary points coordinates of each parcel based on geodetic system (Bessel ellipsoid) into world geodetic system (GRS80 ellipsoid) coordinates, in the past using a survey or distance station calculation. Has its purpose.

에 의해 계산으로 산출하여 각 도근점에 대응하는 세계측지계 좌표를 취득하는 단계(단계 S3) 및;
상기 단계 S3에서 각 도근점(29001∼29007)에 대응하는 세계측지계 좌표가 계산에 의해 취득되면, 지역측지계 좌표를 산출할 때 현장에서 관측한 각 도근점에서 각 필지 경계점(1~37)까지의 관측각과 거리를 경계점 관측부로부터 확인한 다음, 경계점 관측부로부터 확인된 각 도근점에서 각 필지 경계점(1~37)까지의 관측각과 거리를 이용하여 각 필지 경계점(1~37)의 세계측지계 좌표를, 상기 수학식 1 및 상기 수학식 2에 의해 계산으로 산출하여 상기 각 필지 경계점(1∼37)에 대응하는 세계측지계 좌표를 취득하는 단계(단계 S4);를 갖추어 이루어진 것을 특징으로 하는 필지 경계점 좌표를 세계측지계 좌표로 변환하기 위한 방법을 특징으로 한다..The present invention for achieving the above object,
The presence of four known points (two points of departure and departure (98210,98211) and two points of arrival (24263,24262) on which the GPS surveying device is to be installed) exists when the cadastral survey control point mark is installed on the ground. Determining that there is no loss and determining that it does not exist (step S1);
If it is determined in step S1 that the four known points exist, the GPS measurement apparatus is mounted on the four known points, and then GPS observation is performed on the four known points and the four points of the four points are based on the GPS observations. Acquiring world geodetic coordinates corresponding to the known point (step S2);
When the global geodetic coordinates are acquired for the four known points in step S2, when the existing geodetic coordinates are calculated for each approach point 29001 to 29007, the observation angle and distance observed in the field are confirmed by the Cadastral Surveying Department. , Using the observation angle θ and distance D to each approach point 29001 to 29007, the coordinates of the world geodetic system corresponding to each approach point 29001 to 29007,
Equation 1 Δx = Xb-Xa,
Δy = Yb-Ya
(Where Δx is the longitudinal difference, Δy is the horizontal difference, Xa, Xb, Ya, and Yb are coordinate values of the X and Y axes), and
Equation 2

Calculating by calculation to obtain world geodesic coordinates corresponding to each approach point (step S3);
If the world geodetic coordinates corresponding to each of the approach points 29001 to 29007 are obtained by calculation in step S3, the angles of observation from each approach point observed in the field to each parcel boundary point 1 to 37 when calculating the geodetic coordinates are calculated. After confirming the distance from the boundary point observer, the coordinates of the world geodetic coordinates of each parcel boundary point (1 to 37) using the observation angle and distance from each approach point identified from the boundary point observer to each parcel boundary point (1 to 37), And calculating the geodetic coordinates corresponding to the respective parcel boundary points 1 to 37 (step S4), calculated by Equation 1 and Equation (2). It features a method for converting to coordinates.

In addition, the present invention, if the GPS surveying device is impossible to install the GPS surveying device at the known point (starting, arrival side) in the step S1, the GPS surveying device is installed at the approach point capable of GPS surveying, the GPS surveying device Acquiring the global geodetic coordinates for the installed approach point (step S5); If the world geodetic coordinates are acquired for the approach point in which the GPS surveying device is installed in step S5, after calculating the azimuth and distance by inverse calculation using existing geodetic coordinate values, acquiring the world geodetic coordinates for the known point ( Step S6); characterized in that made more.

In addition, the present invention, after calculating the azimuth and distance by the inverse calculation in the step S6 to obtain the world geodetic coordinates for the known point, and proceeds to the step S3 to calculate the world geodetic coordinates for each of the approach point (29001 ~ 29007) And calculate the global geodetic coordinates of each parcel boundary point.

In addition, the present invention, the observation angle and distance from each of the approach point (29001 to 29007) to each of the parcel boundary (1 to 37) is referred to the results (boundary observation section) observed in the existing geodetic system, the world of each parcel boundary point To calculate the geodetic coordinates.

According to the present invention, when the boundary coordinates registered as geodetic system using the conversion method of the present invention are converted into world geodetic coordinates, the site survey cost can be greatly reduced, and the coordinate conversion performance is the same as the measured performance. Also, it can be used as cadastral survey results such as boundary restoration survey by managing conversion records (measuring point of a survey instrument, azimuth, distance, etc.).

Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings.

6 is a flowchart for explaining an embodiment according to the present invention. An embodiment according to the present invention will be described based on FIG. 6 with reference to FIGS. 2 to 5.

2 is a view showing a known point according to an embodiment of the present invention. Here, the known point has two points of the starting side consisting of a straight line connecting two points (98210, 98211) and two points of the arrival side consisting of a straight line connecting two points (24263, 24262).

First, it is determined whether there is a known point consisting of two points 98210 and 98211 of the starting side and two points 24263 and 24262 of the arrival side on which the GPS surveying device is to be installed (step S1). That is, if the cover of the cadastral survey control point installed on the ground is buried, it is determined that four known points exist on the ground, and if the cover of the cadastral survey control point is lost, the four known points do not exist on the ground. . In step S1, the cover of the cadastral survey control point installed on the ground surface is buried for four known points, and it is determined that there are two points (98210, 98211) of the starting side and two points (24263, 24262) of the arrival side, which are known points. If a GPS surveying device is mounted on four known points (two points of departure (98210,98211) and two points of arrival (24263,24262)), GPS observation is performed and corresponding four points of reference are based on the GPS observation. World geodetic system coordinates are acquired (step S2).

Here, GPS positioning is performed using a static or network RTK method, and acquires the geodetic coordinates of each point based on the GPS reference performance provided by the national GPS reference point service.

3 is a view for explaining the calculation of the approximate world geodetic coordinates according to an embodiment of the present invention.

In step S2, if the global geodetic coordinates are acquired for four known points (two points of departure (98210, 98211) and two points of arrival (24263, 24262)) based on the GPS survey, the results observed in the existing geodetic system That is, when the existing geodetic coordinates are calculated, the observation angles and distances recorded in the cadastral surveying department are observed in the field, and then, to each detour point (29001 to 29007) shown in FIG. 3 recorded in the cadastral surveying department. Using the observation angle θ and the distance D, the world geodesic coordinates corresponding to each approach point 29001 to 29007 are calculated by the equations 1 and 2 described later using Figs. To obtain (step S3).

Here, methods of cabinet, distance observation, azimuth calculation, calculation and allocation of angular errors, vertical and horizontal errors, and tolerances are determined by cadastral survey calculation methods determined by cadastral method.

4 is a view showing a boundary point of a parcel according to an embodiment of the present invention, Figure 5 is a view for explaining the world geodetic coordinate calculation for the parcel boundary point.

When the world geodetic coordinates corresponding to the respective approach points 29001 to 29007 are calculated and acquired in step S3, the respective parcel boundary points 1 to 37 shown in FIG. 4 are confirmed, and then, as shown in FIG. Compute the world geodetic coordinates of each parcel boundary. Here, the observation angle and distance from each approach point (29001 to 29007) to the parcel boundary points (1 to 37) are the results observed in the existing geodetic system, that is, from each approach point (29001 to 29007) to each parcel boundary point (1 to 37). When calculating the geodetic coordinates of, check the observation angle and distance observed from the site from the boundary observation and coordinate calculation unit, and then from each approach point (29001 to 29007) identified from the boundary observation and coordinate calculation unit, to each parcel boundary point (1 to 37). The geodetic coordinates corresponding to the parcel boundary points 1 to 37 are calculated by calculation using Equation 1 and Equation 2 described later using Figs. ).

If the parcel boundary point has not been observed in the existing geodetic system, the azimuth and distance are inversely calculated using the longitudinal and horizontal coordinates of the parcel boundary point registered in the boundary point coordinate register.

On the other hand, if the GPS surveying device is installed at the known point (starting and arriving side) in step 1, but the known point is lost on the surface or the GPS survey cannot be performed due to the observation obstacle, the GPS surveying is possible at the approach point where GPS surveying is possible. The device is installed and the world geodetic coordinates are acquired for the approach point in which the GPS surveying device is installed (step S5).

If the global geodetic coordinates are acquired for the approach point on which the GPS surveying device is installed in step S5, after calculating the azimuth and distance by inverse calculation using the existing geodetic coordinate values, the global geodetic coordinates for the known point are acquired (step S6). ). Subsequently, the process proceeds to step S3 to calculate the world geodetic system coordinates for the approach points 29001 to 29007 as described above, and then calculates the world geodetic system coordinates of each parcel boundary point.

On the other hand, if there is no survey result at the time of registration, such as the cadastral proximity network, known azimuth azimuth and distance calculation unit, angle measurement and distance calculation unit, cadastral survey measurement unit, boundary point observation and coordinate calculation unit shown in Figs. The azimuth and distance are obtained by inverse calculation using the coordinates (X, Y) registered in the surface register and the boundary point coordinate register.

Here, the azimuth and distance inverse calculation will be described. When there are two known points A (Xa, Ya) and B (Xb, Yb), the azimuth and distance between the two points can be obtained by the following calculation method. It is called.

Known point azimuth calculation is demonstrated using FIG.7 and FIG.8.

FIG. 7 is a diagram illustrating a vertical line difference for explaining a known point azimuth calculation according to the present invention, and since the coordinates of points A and B are known, the difference between the two points, that is, the vertical line difference Δx and the horizontal line difference Δy Can be obtained by the following equation.

Equation 1

Δx = Xb-Xa,

Δy = Yb-Ya

 In FIG. 7, cos θ, sin θ, and tan θ can be obtained by the following equation 2 around the θ angle.

Equation 2

8 is a view for explaining the azimuth and azimuth angle according to the present invention, the azimuth angle used in the cadastral survey is an angle measured in the clockwise direction with respect to North Korea, the θ angle as shown in Figure 8 I, Corresponds to the θ angle of any of the upper limits of II, III, and IV. Therefore, in the case of the upper limit of I, the angle of θ becomes the azimuth as it is, and in the case of the upper limit of II, the azimuth is less than 180 degrees by the angle of θ, so the azimuth is 180 ° -θ. Similarly, at the upper limit of III, it is 180 ° + θ because it has passed by 180 ° more than θ angle, and at the upper limit of IV, it is less than 360 ° by the angle of θ, resulting in 360 ° -θ.

In summary, as shown in Table 1 below, the angle θ is obtained by the vertical line difference Δx and the horizontal line difference Δy, and then the azimuth angle can be obtained according to the sign of the vertical line difference Δx and the horizontal line difference Δy.

                         TABLE 1

FIG. 1 is a diagram for explaining two-dimensional coordinate transformation. FIG. 1A is a reference shape, and FIG. 1B is an X-axis based on the shape shown in FIG. 1 is a state in which origin movement has occurred, FIG. 1C is a state in which rotation has occurred, and FIG. 1D is a view showing a state in which scale change (enlargement) has occurred,

2 is a view showing a known point according to an embodiment of the present invention;

3 is a view for explaining the calculation of the approximate point geodetic coordinates according to an embodiment of the present invention;

4 is a view showing the boundary point of the parcel in accordance with an embodiment of the present invention,

5 is a diagram for explaining calculation of world geodetic coordinates for parcel boundary points;

6 is a flowchart illustrating an embodiment according to the present invention;

7 is a view showing a longitudinal line difference for explaining a known point azimuth calculation according to the present invention;

8 is a view for explaining the azimuth and azimuth angle according to the present invention.

Claims (4)

The presence of four known points (two points of departure and departure (98210,98211) and two points of arrival (24263,24262) on which the GPS surveying device is to be installed) exists when the cadastral survey control point mark is installed on the ground. Determining that there is no loss and determining that it does not exist (step S1); If it is determined in step S1 that the four known points exist, the GPS measurement apparatus is mounted on the four known points, and then GPS observation is performed on the four known points and the four points of the four points are based on the GPS observations. Acquiring world geodetic coordinates corresponding to the known point (step S2); When the global geodetic coordinates are acquired for the four known points in step S2, when the existing geodetic coordinates are calculated for each approach point 29001 to 29007, the observation angle and distance observed in the field are confirmed by the Cadastral Surveying Department. , Using the observation angle θ and distance D to each approach point 29001 to 29007, the coordinates of the world geodetic system corresponding to each approach point 29001 to 29007,  Equation 1 Δx = Xb-Xa,           Δy = Yb-Ya (Where Δx is the longitudinal difference, Δy is the horizontal difference, Xa, Xb, Ya, and Yb are coordinate values of the X and Y axes), and Equation 2

Calculating by calculation to obtain world geodesic coordinates corresponding to each approach point (step S3); If the world geodetic coordinates corresponding to each of the approach points 29001 to 29007 are obtained by calculation in step S3, the angles of observation from each approach point observed in the field to each parcel boundary point 1 to 37 when calculating the geodetic coordinates are calculated. After confirming the distance from the boundary point observer, the coordinates of the world geodetic coordinates of each parcel boundary point (1 to 37) using the observation angle and distance from each approach point identified from the boundary point observer to each parcel boundary point (1 to 37), And calculating the geodetic coordinates corresponding to the respective parcel boundary points 1 to 37 (step S4), calculated by Equation 1 and Equation (2). Method for converting to coordinates. The method according to claim 1, wherein in the step S1, when the GPS surveying device cannot be installed due to the impossibility of installing the GPS surveying device at the known point (starting and arriving side), the GPS surveying device is installed at a nearby point where the GPS surveying is possible. Acquiring world geodesic coordinates for the approach point in which the device is installed (step S5); If the world geodetic coordinates are acquired for the approach point in which the GPS surveying device is installed in step S5, after calculating the azimuth and distance by inverse calculation using existing geodetic coordinate values, acquiring the world geodetic coordinates for the known point ( Step S6); further comprising parcel boundary coordinates to the world geodetic coordinates. The method according to claim 2, wherein after calculating the azimuth and distance by inverse calculation in step S6, the world geodetic coordinates for the known point are obtained, and then the process proceeds to the step S3 to determine the world geodetic coordinates for the respective approach points 29001 to 29007. And calculating the geodetic coordinates of each parcel boundary point. delete

Images