MEDIAN LINE METHOD IN MARITIME DELIMITATION BASED ON EARTH ELLIPSOID

Abstract

The present disclosure discloses a median line method in maritime delimitation based on an earth ellipsoid. The method is implemented by calculating a coordinate point of a median line of a maritime space, including: determining two delimitation base points A and B on a coastline of a certain country and a delimitation base point C on a coastline of the other country; calculating to obtain an equidistant point O with the same geodetic distance from the two points A and B; solving the geodetic distance from the point O to the point C according to the geodetic coordinates of the equidistant point O; determining whether the difference between respective geodetic distances of OC and OA is less than a first target error; if the difference is less than the first target error, then taking O as a point on a solved median line; and otherwise, adjusting the geodetic distance, and recalculating the coordinates of the point O.

Claims

1. A median line method in maritime delimitation based on an earth ellipsoid, comprising: determining respective geodetic coordinates of two delimitation base points A and B and a geodetic distance between the two points on a coastline of a certain country; and determining geodetic coordinates of a delimitation base point C on a coastline of the other country; calculating to obtain geodetic coordinates of an equidistant point O with the same geodetic distance from the two points A and B; solving the geodetic distance from the point O to the point C in accordance with a formula of an inverse solution of a geodetic problem according to the geodetic coordinates of the equidistant point O; determining whether the difference between respective geodetic distances of OC and OA is less than a first target error; if so, taking the point O as an equidistant point of the three delimitation base points A, B, and C; and if not, recalculating the geodetic coordinates of the point O after adjusting the geodetic distance of OA.

2. The median line method in maritime delimitation based on an earth ellipsoid according to claim 1, wherein the geodetic distance of OA is adjusted through formula 4:
Δs=s+λ/2  4 in the formula: s is the geodetic distance of OA, and λ is the difference between respective geodetic distances of OC and OA.

3. The median line method in maritime delimitation based on an earth ellipsoid according to claim 1, comprising: scanning whether there are other delimitation base points in an area by taking the point O as a center and taking the geodetic distance of the two points A and B as a radius, wherein if there are no other delimitation base point, then point O is taken as a solved inflection point of a boundary; and if there are other delimitation base points, then the point O is not a desired boundary point, the geodetic coordinates of points A, B, and C are redetermined.

4. The median line method in maritime delimitation based on an earth ellipsoid according to claim 1, wherein the calculating to obtain geodetic coordinates of an equidistant point O with the same geodetic distance from the two points A and B comprises that: respective geodetic coordinates of the two points A and B and the geodetic distance between the two points on the earth ellipsoid are known; point O is set as an equidistant point with the same geodetic distance from the two points A and B, and the distance is known; geodetic coordinates of three points A, B, and point O′ equidistant from the two points A and B are calculated in accordance with a formula of a direct solution of a geodetic problem according to the geodetic coordinates of the point A, an approximate value of a geodetic azimuth of AO, and the geodetic distance of AO; respective geodetic distance of O′A and O′B are calculated in accordance with a formula of an inverse solution of a geodetic problem according to the geodetic coordinates of three points A, B, and C; whether the difference between the respective geodetic distances of the O′A and O′B is less than a second target error is determined; if so, the geodetic coordinates of the point are taken as the geodetic coordinates of the point O; and if not, the geodetic coordinates of the point O′ are recalculated after adjusting an approximate value of the geodetic azimuth of AO.

5. The median line method in maritime delimitation based on an earth ellipsoid according to claim 4, wherein the geodetic coordinates (B.sub.O′, L.sub.O′) of the approximate point and the reverse geodetic azimuth of A.sub.O′A of O′A are calculated through formula 1: $\begin{array}{cc}\{\begin{array}{c}{B}_{O^{\prime }}=f_B\left(B_A,L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\\ L_{O^{\prime }}=f_L\left(B_A,L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\\ A_{O^{\prime }A}=f_A\left(B_{\mathrm{PA}},L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\end{array}& 1\end{array}$ in the formula, B.sub.A and L.sub.A are respectively the geodetic coordinates of point A, S.sub.AO is the geodetic distance between the point A and the point O, and A.sub.AO is an approximate value of the geodetic azimuth of AO.

6. The median line method in maritime delimitation based on an earth ellipsoid according to claim 5, wherein the approximate value A.sub.O′A of the geodetic azimuth of AO is calculated by the method as follows: approximating a triangle ABO as a plane triangle; setting ∠BAO=α, obtaining cos α=s/2 according to a trigonometric cosine theorem a.sup.2 b.sup.2 c.sup.2−2bc*cos A, and then obtaining α=arccos (s/2r), wherein the approximate value of the geodetic azimuth of AO is A.sub.AO, when A.sub.AB−A.sub.AC>0 then A.sub.Ao=A.sub.AB−α; when A.sub.AB−A.sub.AC<0, A.sub.Ao=A.sub.AB+α; A.sub.AB is a geodetic azimuth of a geodetic line AB; A.sub.AC is an azimuth of a geodetic line AC; s is the geodetic distance between the two points A and B; and r is the geodetic distance from point O to point A or point B.

7. The median line method in maritime delimitation based on an earth ellipsoid according to claim 5, wherein respective geodetic distances S.sub.O′A and S.sub.O′Q of O′A and O′B are calculated through formula 2 and formula 3: $\begin{array}{cc}\{\begin{array}{c}{S}_{O^{\prime }A}=f_s^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_A\right)\\ A_{O^{\prime }A}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_A\right)\\ A_{{\mathrm{AO}}^{\prime }}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_{\mathrm{PA}}\right)\end{array}& 2\end{array}$ $\begin{array}{cc}\{\begin{array}{c}{S}_{O^{\prime }B}=f_s^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_B\right)\\ A_{O^{\prime }B}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_{\mathrm{PB}}\right)\\ A_{{\mathrm{BO}}^{\prime }}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_B\right)\end{array}& 3\end{array}$ in the formula: A.sub.O′A is a geodetic azimuth of O′A, A.sub.AO′ is a reverse geodetic azimuth of O′A (B.sub.B,L.sub.B) is geodetic coordinates of the point B, A.sub.O′B is a geodetic azimuth of O′B, and A.sub.BO′ is a reverse geodetic azimuth of O′B.

8. The median line method in maritime delimitation based on an earth ellipsoid according to claim 7, wherein the approximate value of the geodetic azimuth of AO is adjusted to obtain ΔA.sub.AO, and the geodetic coordinates of the point O′ is recalculated by substituting ΔA.sub.AO into formula 1.

9. The median line method in maritime delimitation based on an earth ellipsoid according to claim 8, wherein ΔA.sub.AO=A.sub.AO+180*δ/π*r, δ is the difference between S.sub.O′A and S.sub.O′Q, and r is the geodetic distance from the point O to the point A or the point B.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0037] To describe the technical solutions in the embodiments of the present disclosure or in the prior art more clearly, the following briefly describes the drawings required for describing the embodiments. Apparently, the drawings in the following description are merely some embodiments of the present disclosure, and those of ordinary skill in the art may still derive other drawings from these drawings without creative efforts.

[0038] FIG. 1 is a schematic diagram of a construction process of a median line method in maritime delimitation based on an earth ellipsoid provided by the embodiment of the present disclosure.

[0039] FIG. 2 is a schematic diagram of a calculation process of a method for solving an equidistant point of two points on an earth ellipsoid provided by the embodiment of the present disclosure.

[0040] FIG. 3 is a schematic diagram of a coordinate position of a median line of a maritime area obtained by two methods in the embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0041] Technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in the embodiments of the present disclosure. Apparently, the described embodiments are merely part rather than all of the embodiments of the present disclosure. On the basis of the embodiments of the present disclosure, all other embodiments obtained by those of ordinary skill in the art without creative effort fall within the scope of protection of the present disclosure.

[0042] First, the calculation of the distance and the azimuth on an earth ellipsoid is described:

[0043] The calculation of the distance and the azimuth performed on the earth ellipsoid needs to use an earth ellipsoid geodetic problem calculation method, including a direct solution of a geodetic problem and an inverse solution of a geodetic problem. Assuming that a point on the earth ellipsoid has a geodetic longitude L and a geodetic latitude B, a geodetic length between two points is S, and a positive geodetic azimuth and a reverse geodetic azimuth are A and A. Since there are many formulas that can realize the calculation of geodetic problems on the earth ellipsoid. The embodiment of the application only describes by taking Vincenty formula as an example. It can be understood that any other formula that can realize the calculation of a geodetic problem on the earth ellipsoid is also applicable to the method provided by the present application.

[0044] The direct solution of the geodetic problem is a process of calculating geodetic coordinates (B.sub.Q,L.sub.Q) of the other point and the reverse geodetic azimuth A.sub.QP according to the known geodetic coordinates (B.sub.P,L.sub.P) of a certain point P on a geodetic line, the geodetic length S.sub.PQ from the point P to the other point Q on the geodetic line, and a geodetic azimuth A.sub.pQ. That is, an equation set is solved:

[00003]$\{\begin{array}{c}{B}_Q=f_B\left(B_P,L_P,S_{\mathrm{PQ}},A_{\mathrm{PQ}}\right)\\ L_Q=f_L\left(B_P,L_P,S_{\mathrm{PQ}},A_{\mathrm{PQ}}\right)\\ A_{\mathrm{QP}}=f_A\left(B_P,L_P,S_{\mathrm{PQ}},A_{\mathrm{PQ}}\right)\end{array}$

[0045] The inverse solution of the geodetic problem is a process of calculating the geodetic length S.sub.PQ between the two points and the positive geodetic azimuth A.sub.PQ and the reverse geodetic azimuth A.sub.QP according to the known geodetic coordinates (B.sub.P, L.sub.P) and (B.sub.Q, L.sub.Q) of two different points P and Q on the geodetic line. That is, an equation set is solved:

[00004]$\{\begin{array}{c}{S}_{\mathrm{PQ}}=f_S^{\prime }\left(B_P,L_P,B_Q,L_Q\right)\\ A_{\mathrm{PQ}}=f_A^{\prime }\left(B_P,L_P,B_Q,L_Q\right)\\ A_{\mathrm{QP}}=f_{A^{\prime }}^{\prime }\left(B_P,L_P,B_Q,L_Q\right)\end{array}$

EMBODIMENT

[0046] The embodiment of the present disclosure provides a median line method in maritime delimitation based on an earth ellipsoid. The method may include that:

[0047] Respective geodetic coordinates of two delimitation base points A and B and a geodetic distance between the two points on a coastline of a certain country are determined. Geodetic coordinates of a delimitation base point C on a coastline of the other country is determined.

[0048] Geodetic coordinates of an equidistant point O with the same geodetic distance from the two points A and B are calculated by a method for solving an equidistant point of two points on the earth ellipsoid provided by Embodiment 1.

[0049] The geodetic distance from point O to point C is solved in accordance with a formula of an inverse solution of a geodetic problem according to the geodetic coordinates of the equidistant point O.

[0050] Whether the difference between respective geodetic distances of the OC and OA is less than a first target error is determined.

[0051] If so, the point O is taken as an equidistant point of the three delimitation base points A, B, and C.

[0052] If not, the geodetic coordinates of the point O is recalculated after adjusting the geodetic distance of OA.

[0053] In the abovementioned mode, the equidistant point of the three points is calculated by finding other combinations of delimitation base points on the coastlines of both countries, until all combinations are used, and finally, various points are connected in sequence to form a median line of a maritime space.

[0054] Specifically, the geodetic distance of OA is adjusted through formula 4:

Δs=s+λ/2  4

[0055] In the formula, s is the geodesic distance of OA, and λ is the difference between respective geodesic distances of OC and OA.

[0056] Further, whether there are other delimitation base points in an area is scanned by taking the point O as a center and taking the geodetic distance of the two points A and B as a radius. If there are no other delimitation base points, then point O is a desired inflection point of a boundary; and if there are other delimitation base points, then the point O is not a desired boundary point, and the geodetic coordinates of points A, B, and C are redetermined.

[0057] The embodiment of the application provides a new method for generating a median line of a maritime area based on an earth ellipsoid, that is, an equidistant point of three points is directly calculated on the earth ellipsoid, which completely gets rid of the influence of map projection. The core of the method is to convert a three-point equidistant problem of an earth ellipsoid into solving an equidistant point of two points first, then solving an equidistant point of a third point, and finally performing condition determination. The point that meets an error requirement is an equidistant point of the three points.

[0058] Respective geodetic coordinates of two points P and Q and the geodetic distance between the two points may be calculated by a plurality of methods in practice. For example, in an implementation mode, the embodiment of the present application may provide the calculation of geodetic coordinates of an equidistant point O with the same geodetic distance from the two points A and B, which includes that:

[0059] respective geodetic coordinates of the two points A and B and the geodetic distance between the two points on the earth ellipsoid are known;

[0060] the point O is set as an equidistant point with the same geodetic distance from the two points A and B, and the distance is known;

[0061] geodetic coordinates of an approximate point O′ equidistant from the two points A and B are calculated in accordance with a formula of a direct solution of a geodetic problem according to the geodetic coordinates of point A, an approximate value of a geodetic azimuth of AO, and the geodetic distance of AO.

[0062] Specifically, the geodetic coordinates (B.sub.O′, L.sub.O′) of the approximate point O′ and the reverse geodetic azimuth A.sub.O′A of AO′ are calculated through formula 1:

[00005]$\begin{array}{cc}\{\begin{array}{c}{B}_{O^{\prime }}=f_B\left(B_A,L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\\ L_{O^{\prime }}=f_L\left(B_A,L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\\ A_{O^{\prime }A}=f_A\left(B_{\mathrm{PA}},L_A,S_{\mathrm{AO}},A_{\mathrm{AO}}\right)\end{array}& 1\end{array}$

[0063] In the formula, B.sub.A, L.sub.A are respectively the geodetic coordinates of point A, S.sub.AO is the geodesic distance between the point A and the point O, and A.sub.AO is an approximate value of the geodetic azimuth of AO.

[0064] The approximate value A.sub.AO of the geodetic azimuth of AO is calculated by the method as follows.

[0065] A triangle ABO is as approximated as a plane triangle:

[0066] ∠BAO=α is set, cos α=s/2r is obtained according to a trigonometric cosine theorem a.sup.2=b.sup.2+c.sup.2−2bc*cos A, and then α=arccos (s/2r) is obtained. The approximate value of the geodetic azimuth of AO is A.sub.AO, when A.sub.AB−A.sub.AC>0, A.sub.Ao=A.sub.AB−α; and when A.sub.AB−A.sub.AC<0, A.sub.Ao=A.sub.AB+α. A.sub.AB is a geodetic azimuth of a geodesic line AB; A.sub.AC is an azimuth of a geodesic line AC; s is the geodesic distance of the two points A and B; and r is the geodesic distance from point O to point A or point B.

[0067] Respective geodetic distances S.sub.O′A and S.sub.O′Q of the O′A and O′B are calculated through formula 2 and formula 3:

[00006]$\begin{array}{cc}\{\begin{array}{c}{S}_{O^{\prime }A}=f_s^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_A\right)\\ A_{O^{\prime }A}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_A\right)\\ A_{{\mathrm{AO}}^{\prime }}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_A,L_{\mathrm{PA}}\right)\end{array}& 2\end{array}$$\begin{array}{cc}\{\begin{array}{c}{S}_{O^{\prime }B}=f_s^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_B\right)\\ A_{O^{\prime }B}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_{\mathrm{PB}}\right)\\ A_{{\mathrm{BO}}^{\prime }}=f_A^{\prime }\left(B_{O^{\prime }},L_{O^{\prime }},B_B,L_B\right)\end{array}& 3\end{array}$

[0068] In the formula: A.sub.O′A is a geodetic azimuth of O′A, A.sub.AO is a reverse geodetic azimuth of O′A, (B.sub.B,L.sub.B) is geodetic coordinates of point B, A.sub.O′B is a geodetic azimuth of O′B, and A.sub.BO′ is a reverse geodetic azimuth of O′B.

[0069] Respective geodetic distances of O′A and O′B are calculated in accordance with a formula of an inverse solution of a geodetic problem according to the geodetic coordinates of three points A, B, and O′.

[0070] Whether the difference between the geodetic distances of the O′A and O′B is less than a second target error is determined. The second target error may be determined according to the actually required calculation accuracy. For example, the distance between two selected two points is 200 nautical miles, and then the first target error may be set as 0.1 meter.

[0071] If so, the geodetic coordinates of the point O′ are taken as the geodetic coordinates of the point O.

[0072] If not, the geodetic coordinates of the point O′ are recalculated after adjusting the approximate value of the geodetic azimuth of AO.

[0073] The approximate value of the geodetic azimuth of AO is adjusted to obtain ΔA.sub.40, and the geodetic coordinates of the point O′ is recalculated by substituting ΔA.sub.AO into formula 1. ΔA.sub.AO=A.sub.AO+180*δ/π*r, S is the difference between S.sub.O′A and S.sub.O′Q, and r is the geodetic distance from the point O to the point A or the point B.

[0074] The median line method in maritime delimitation based on an earth ellipsoid of the present disclosure will be described in detail below through a specific example.

[0075] The method for solving a boundary point of a maritime area provided by the present application is described in detail below through a specific example.

[0076] Referring to FIG. 1, two delimitation base points A((B.sub.A, L.sub.A) and B ((B.sub.B,L.sub.B) on a coast of a certain country, and one delimitation base point C(B.sub.C, L.sub.C) on a coast of the other country are known, the geodesic distance of the two points A and B is s, and a point P with equal geodesic distance from point A, point B, and point C is solved.

[0077] A calculation process is as follows.

[0078] Step one: an equidistant point between point A and point B is solved. In order to facilitate calculation and understanding, an equilateral triangle □ABO is constructed with by taking two points A and B on the coast of a country as vertices and the length s of AB as the length of a side. According to a method provided by the embodiment to solve the equidistant point of two points on an earth ellipsoid, an equidistant point of points A and B are iteratively calculated until the point that meets a specified error is the solved point O.

[0079] Step two: the distance from a third point is solved. After the accurate position of point O is solved, the distance d from point O to point C of the other country is solved in accordance with a formula of an inverse solution of a geodetic problem, and λ=d−s is calculated. If 2 is less than the specified error, the point O is an equidistant point of the three points. On the contrary, s is adjusted, s=s+λ/2, and point O is recalculated by substituting s into step one, and iterative calculation is performed until point O that satisfies an error condition is solved.

[0080] Step three: condition determination. Assuming that the finally solved equidistant point of the three points is P, then whether there are other delimitation base points in an area is scanned by taking point P as a center and s as a radius. If there are no other delimitation base points, then point P is a desired inflection point of a boundary. On the contrary, point P is not a solved boundary point, and three points must be reselected to calculate according to the abovementioned steps.

[0081] Method for calculating an equidistant point of two points on an earth ellipsoid

[0082] Referring to FIG. 2, assuming that two points A(B.sub.A, L.sub.A) and B((B.sub.B,L.sub.B) are on the earth ellipsoid, the geodetic distance of the two points A and B is s, the geodesic azimuth from the point A to the point B is A.sub.AB, and the point O is an equidistant point with the geodetic distance r from point A and point B, the geodetic coordinates of point O are solved. The triangle ABO is approximated as a plane triangle, and the geodetic coordinates of point O are calculated through iterative approximation. A calculation process is as follows:

[0083] (1) ∠BAO is calculated. ∠BAO=α is set, cos α=s/2r is obtained according to trigonometric cosine theorem a.sup.2=b.sup.2+c.sup.2−2bc*cos A, and then α=arccos (s/2r).

[0084] (2) An approximate point O′ with equal distance from the point A and the point B is calculated. The geodetic coordinates (B.sub.O′,L.sub.O′) of the approximate point O′ of point O is calculated in accordance with a formula of a direct solution of a geodetic problem according to the geodetic coordinates (B.sub.A, L.sub.A) of the point A. The approximate value of the geodetic azimuth of AO is A.sub.AO, when A.sub.AB−A.sub.AC>0, A.sub.AO=A.sub.AB−α; when A.sub.AB−A.sub.AC<0, A.sub.Ao=A.sub.AB+α; A.sub.AB is a geodetic azimuth of a geodetic line AB; A.sub.AC is an azimuth of a geodetic line AC; s is the geodetic distance of the two points A and B; and r is the geodetic distance from the point O to point A or point B.

[0085] (3) An actual point O with equal distance from point A and point B is calculated. The geodetic distance r1 and r2 of O′A and O′B are calculated in accordance with a formula of an inverse solution of a geodetic problem according to the geodetic coordinates of three points A(B.sub.A, L.sub.A), B((B.sub.B,L.sub.B), and O′(B.sub.O′,L.sub.O′), where δ=r.sub.1−r.sub.2. If δ is less than a specific error, the point O is an equidistant point from point A and point B. Otherwise, A.sub.AO=A.sub.AO+180*δ/π*r, then A.sub.AO is substituted in step (2) and step (3) to recalculate the point O′ until δ meets a specified error, and the solved point O is an equidistant point from the point A and the point B.

[0086] The method for solving an equidistant point from two points on an earth ellipsoid provided by the present application is more concise through an iterative approximation calculation algorithm, and can calculate a result that meets the requirement of accuracy without complex high-power operation. Meanwhile, distance calculation is performed on the earth ellipsoid, which gets rid of the influence of map projection and does not need to consider the selection of various map projection methods.

[0087] The calculation accuracy of the method for solving a coordinate point of a median line of a maritime area based on an earth ellipsoid provided by the application is verified.

[0088] A maritime area with east-west opposite coasts of both parties is selected as a virtual delimitation maritime space. The west side of the maritime area is party A, and the east side of the maritime area is party B. A median line of an overlapped maritime area of both parties is reckoned by taking this as an example. Nine delimitation base points are selected along the coast of party A, and eleven delimitation base points are selected along the coast of party B (see Table 1). The coordinate unit is decimal degree, the calculation result retains 6 decimal places, “-” means that the point is located in south latitude or west longitude, and the distance calculation error is 10.sup.−4 m.

TABLE-US-00001 TABLE 1 WGS-84 coordinate system of starting coordinates of delimitation between both parties Coordinates of party A Coordinates of party B Number Longitude (°) Latitude (°) Longitude (°) Latitude (°) 1 119.018231 0.945006 120.081850 0.752270 2 118.804079 0.795100 120.003328 0.688024 3 118.025996 0.766546 119.782038 0.202614 4 117.647661 −0.432702 119.617855 0.002739 5 117.633385 −0.803898 119.289489 −1.325000 6 116.612596 −2.210159 118.789803 −2.645601 7 116.41986 −3.423684 118.761249 −2.845475 8 116.298507 −3.966201 118.832633 −3.445100 9 116.06294 −4.087554 118.932570 −3.637836 10 — — 119.068200 −3.873403 11 — — 118.989677 −4.201768

[0089] The delimitation base points of both parties are taking as starting data, and coordinates of a median line of a maritime area generated by the method and the Geocap MLB 2.3.4 for ArcGIS10.2, which was developed by Geocap Company in Norway are compared (see Table 2).

TABLE-US-00002 TABLE 2 Delimitation result comparison “Three-point method” of Method provided by Geocap software the present application Number Longitude (°) Latitude (°) Longitude (°) Latitude (°) 1 119.545193 0.950153 119.620540 1.242777 2 119.465975 0.642444 119.465975 0.642444 3 119.281997 0.450620 119.281997 0.450620 4 119.225873 0.414424 119.225873 0.414424 5 118.630625 −0.205166 118.630625 −0.205166 6 118.688214 −0.469236 118.688214 −0.469236 7 118.588502 −0.655162 118.588502 −0.655162 8 118.267305 −1.689246 118.267305 −1.689246 9 117.786731 −1.994681 117.786731 −1.994681 10 117.669983 −2.585315 117.669983 −2.585315 11 117.553674 −2.983007 117.553674 −2.983007 12 117.627550 −3.285265 117.627550 −3.285265 13 117.621451 −3.978886 117.621451 −3.978886 14 117.643663 −4.087804 117.643663 −4.087804 15 117.633767 −4.201768 117.435776 −6.473682

[0090] It can be seen from the analysis that the coordinates of an inflection point of a median line of a maritime area calculated by the method and the Geocap software are completely consistent except for the coordinates of a starting point and an ending point of the median line of the maritime area (see FIG. 3). “Note: FIG. 3 is a schematic diagram, is a virtual delimitation scenario, is only intended to show the median line delimitation method provided by Embodiment 2 of the present application, and does not prejudice to any claims of the parties in the maritime space.” There is a difference between coordinates of the starting point and the ending point of the median line in the two methods due to the adoption of different delimitation methods: the present method calculates by a boundary extension method, while Geocap directly takes a midpoint of the starting base points of both parties as the coordinates of the starting point of the median line, and takes the midpoint of the ending base points as the coordinates of the ending point of the median line.

[0091] It is also to be noted that relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation herein, and do not necessarily require or imply the existence of any such actual relationship or order between these entities or operations. Moreover, the terms “include”, “contain” or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, article or device including a series of elements not only includes those elements, but also includes those elements that are not explicitly listed, or includes elements inherent to such a process, method, article or device. In the absence of more restrictions, elements defined by the phrase “include a/an . . . ” do not exclude the existence of additional identical elements in the process, method, commodity, or device that includes the elements.

[0092] The above is only a preferred embodiment of the present disclosure and is not intended to limit the scope of the present disclosure. Any modifications, equivalent replacements, improvements and the like made within the spirit and principle of the present disclosure shall fall within the scope of protection of the present disclosure.