CN102938029A - Method for generating extrapolation boundary of continental shelf 200 nautical miles away - Google Patents

Abstract

The invention discloses a method for generating an extrapolation boundary of a continental shelf which is 200 nautical miles away and in particular relates to a method for generating a boundary of a continental shelf in the mode that based on a known fixed point, a certain distance is extrapolated outwards in a specific direction. The method comprises the following steps of: generating extrapolation arcs; intersecting the extrapolation arcs; cutting the extrapolation arcs; deleting internal arcs; and deleting surplus arcs, combining the extrapolation arcs and the like so as to realize the automatic generation of a line of FOS (factor of safety)+60M, a line of 350M and a line of 2500m+100M. The method has the advantages of accuracy in surveying and mapping data, convenience in operation, small influence of natural weather and the like in the process of coastal surveying and mapping and maritime delimitation in China.

Description

A kind of 200 nautical miles with outer shelf extrapolation boundary generation method

Technical field

The present invention relates to a kind of 200 nautical miles with outer shelf boundary generation method, specifically refer to a kind of based on known fixed point, to specific direction extrapolation certain distance, thereby generate the necessary technical method of delimitation of continental shelf beyond 200 nautical miles, FOS+60M line, 350M line and 2500m+100M line etc., relate to computer graphics, submarine science and description field, ocean.

Background technology

" Seal Law Convention " come into force according to the nineteen eighty-two signature, " continental platform of coastal state comprises beyond its territorial waters and naturally extending according to the whole of its territory, land; expand to sea bed and the subsoil of the sea-bed area of limit, continent outer rim; if from baseline from which the breadth of the territorial sea is measured, measure the limit, continent outer rim apart from less than 200 nautical miles, expand to the distance of 200 nautical miles "; If coastal state advocates to measure over 200 nautical miles from baseline from which the breadth of the territorial sea is measured, Outer Limit of Continental Shelfs beyond should delimiting 200 nautical miles according to " pact " the 76 related request with the 4th of annex 2, and submit the description case to Commission on the Limits of the Continental Shelf.

FOS+60M line, 350M line and 2500m+100M line (M is abbreviation in the sea, is the international practice in the establishment of delimitation of continental shelf case beyond 200 nautical miles) are the important foundations of Outer Limit of Continental Shelf beyond determining 200 nautical miles.Particularly; the FOS+60M line is continent toe point (being called for short FOS) the 60 nautical miles boundaries that generate of offward extrapolating; the 350M line is from the territorial waters basic point 350 nautical miles formed boundaries of offward extrapolating, and the 2500m+100M line is the 2500m sea-bottom contour 100 nautical miles boundaries that generate of offward extrapolating.The FOS+60M line is called for short the formula line, and 350M line and 2500m+100M line are called for short again restraining line.Article three, the common trait of boundary is the certain distance of extrapolating from given fixed point to assigned direction, so be collectively referred to as the extrapolation boundary.

Open source information analysis from retrieval, current domestic without the proven technique method for 200 nautical miles generations with description boundaries such as outer shelf FOS+60M line, 350M line and 2500m+100M lines, although the minority document is arranged to be introduced delimitation of continental shelf beyond 200 nautical miles, but, without holding within full and accurate boundary generation method, art methods is not enough to support the generation of the description boundaries such as FOS+60M line, 350M line and 2500m+100M line.

Summary of the invention

The present invention is directed to deficiency of the prior art, propose a kind of set point that utilizes for the center of circle, is that radius is looked unfamiliar into series extrapolation arc at the earth according to certain distance, by intersecting and union operation of extrapolation arc, thus the technical method of generation delimit.Specifically from 60 nautical miles of given continent toe point extrapolations generate the FOS+60M lines, from 350 nautical miles of given territorial waters basic point extrapolations generate 350 nautical miles restraining lines, from the technical methods of 100 nautical miles generation 2500m+100M lines of given 2500m sea-bottom contour point extrapolation.The boundary points such as the necessary FOS+60M line of delimitation of continental shelf, 350M line and 2500m+100M line beyond this invention can be applicable to automatically generate 200 nautical miles, can be applicable to the generation of other kind extrapolation lines of description field, ocean, also can be applicable to the generation of the field extrapolation envelopes such as computer graphics.

The present invention is achieved by following technical proposals:

A kind of 200 nautical miles with outer shelf extrapolation boundary generation method, it is characterized in that, comprise the following steps:

(1) generate the extrapolation arc.

Given initial data set G={g _i(x _i, y _i), each data point g _i(x _i, y _i) all need by the given terrestrial coordinate (x of the mode of longitude and latitude _i, y _i).With (x _i, y _i) put as the center of circle, with identical spherical distance r, generate sphere extrapolation arc a _i.Traversal initial data set G={g _i(x _i, y _i), form data set Arc ₀={ a _i, each extrapolation arc a _iComprise central coordinate of circle (x _i, y _i) and extrapolated radius value r.

(2) the extrapolation arc intersects.

By dual searching loop data set Arc ₀={ a _i, obtain the mutual intersection point set A rc of extrapolation arc ₁={ a1 _i.Ergodic data collection Arc ₀={ a _i, take out successively extrapolation arc a _i.With extrapolation arc a _iHeaded by arc, enter second and recirculate, successively with data set Arc ₀={ a _iIn except a _iAll the other outer extrapolation arcs intersect judgement, and record intersection point forms data set Arc ₁={ a1 _i, each data point a1 _iComprise central coordinate of circle, extrapolated radius and intersection point set C _i={ c _j.Arc ₁={ a1 _iData point a1 _iWith inherited data set Arc ₀={ a _iMiddle data point a _iAll properties.

Two extrapolation arc whether crossing bases for estimation are: the two distance of center circle is from d≤2r, and r is extrapolated radius, and when d<2r, the two has two intersection points, and when d=2r, the two has an intersection point; D > during 2r without intersection point.

Indefinite because of the intersection point sum of each extrapolation arc and other extrapolation arcs, therefore, need the data set C of design stores intersection point _i={ c _j.The all corresponding intersection point data set C of each extrapolation arc _i={ c _j, intersection point c _jShould comprise terrestrial coordinate value (x _j, y _j) and with the azimuth value α of extrapolation arc circle center line connecting _j.Intersection point c _jAccording to azimuth value α _jSize order deposit into data set C _i={ c _j.Intersection point data set C _i={ c _jChain enters data acquisition Arc ₁={ a1 _iThe new data acquisition Arc of rear formation ₂={ a2 _i.With Arc ₁={ a1 _iCompare, except intersection point C _i={ c _jOutside, Arc ₂={ a2 _iInherited Arc ₁={ a1 _iAll properties.

(3) cutting extrapolation arc.

Adopt dual circulation, cutting extrapolation arc forms little segmental arc.Ergodic data collection Arc ₂={ a2 _i, take out successively data point a2 _i, then travel through corresponding with it intersection point data acquisition C _i={ c _j.According to the intersection point order, take out successively two adjacent intersection points and form subarc section a3 _i, and form new data set Arc ₃={ a3 _i.Each data point a3 _iStarting point, terminal, intermediate point and the center of circle that all should comprise the subarc section, and the azimuth angle alpha of starting point and circle center line connecting ₁, terminal and circle center line connecting azimuth angle alpha ₂, and record subarc section a3 _iFrom original arc.So-called intermediate point, refer to and be positioned at subarc section a3 _iStarting point and the point between terminal, for judgement and the mutual relationship of original extrapolation arc.

(4) delete inner arc.

Adopt dual circulation to delete inner subarc section, namely by certain original arc, be completely contained in inner subarc section.Ergodic data set A rc ₀={ a _i, take out successively original arc data point a pointed _i.Each data point a _iEqual ergodic data collection Arc ₃={ a3 _iOnce, judge successively data point a3 _iWith a _iRelation, work as a3 _iSubarc section pointed is by original arc a _iInstitute comprises and not from original arc a _iThe time, this subarc section is inner arc, need be from data acquisition Arc ₃={ a3 _iThe middle deletion.Thereby at data set Arc ₃={ a3 _iForm new data set Arc on basis ₄={ a4 _i.A4 _iInherited Arc ₃={ a3 _iIn a3 _iAll properties.

Subarc section a3 _iWith data acquisition Arc ₀={ a _iMiddle a _iThe relation of inclusion basis for estimation is: group segmental arc a3 _iIntermediate point coordinate and a _iWhen central coordinate of circle is less than extrapolated radius r apart from d, subarc section a3 _iBe contained in a fully _i.

(5) Delete superfluous arc.

The extrapolation direction of FOS+60M line, 350M line and 2500m+100M line is offward, and is not the envelope arc of sealing, and in this and computer graphics, buffer zone has essential difference.Generally speaking, the extrapolation direction D of FOS+60M line, 350M line and 2500m+100M line is set as due east, due south, positive west or positive north.Extrapolation direction D is outside input variable.

Ergodic data collection Arc ₄={ a4 _i, by extrapolation direction D and subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂Relation judgement, α namely ₁Or α ₂Point to when consistent and retain this extrapolation segmental arc with extrapolation direction D, delete the segmental arc that has nothing to do, thereby form new data acquisition Arc ₅={ a5 _i, a5 _iInherited data set Arc ₄={ a4 _iMiddle a4 _iAll properties.Subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂By step (3), calculate.

Subarc section a4 _iAzimuth angle alpha ₁Or α ₂With the extrapolation direction, D concerns that basis for estimation is as follows:

(a) when extrapolation direction D is due east, and α ₁∈ | 0,180| or α ₂∈ | during 0,180|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(b) when extrapolation direction D is due south, and α ₁∈ | 90,270| or α ₂∈ | during 90,270|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(c) as extrapolation direction D when just western, and α ₁∈ | 180,360| or α ₂∈ | during 180,360|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(d) when extrapolation direction D is positive north, and α ₁∈ | 270,360| ∪ | 0,90| or α ₂∈ | 270,360| ∪ | during 0,90|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain.

(6) the extrapolation arc merges.

Adopt dual searching loop data acquisition Arc ₅={ a5 _i, take out at first successively subarc section a5 _i, ergodic data set A rc again ₅={ a5 _i, take out except a5 successively _iAll the other subarc sections, judge whether the two exists neighbouring relations.Whether the beginning or end coordinate that basis for estimation is the two is identical.Starting point and terminal point coordinate calculate by step (3).

By this step, set up all contacting directly of neighbouring relations segmental arc that have, can review all relevant segmental arcs by any segmental arc.Ergodic data set A rc again ₅={ a5 _i, by the mutual relationship between segmental arc, sort in order and form a complete data set Arc ₆={ a6 _i.Extrapolation arc Arc is linked in sequence ₆={ a6 _iForm final boundary L.

(7) the FOS+60M line generates

FOS={fos _iCount according to collection for the continent toe, with FOS={fos _iBe initial data set G _fos, take 60 nautical miles as extrapolated radius r _fos, adopt step (1)～(6), according to given extrapolation direction D, obtain the FOS+60M line;

(8) the 350M line generates

BASE={b _iBe territorial waters basic point data set, with BASE={b _iBe initial data set G _b, take 350 nautical miles as extrapolated radius r _b, adopt step (1)～(6), according to given extrapolation direction D, obtain the 350M line;

(9) the 2500m+100M line generates

TOPO={t _iBe 2500m sea-bottom contour data set, with TOPO={t _iBe initial data set G _t, take 100 nautical miles as extrapolated radius r _t, adopt step (1)～(6), according to given extrapolation direction D, obtain the 2500m+100M line.

Regulation according to the 76, the United Nations's " Seal Law Convention " and delimitation of continental shelf council technical criteria, FOS+60M line, 350M line and 2500m+100M line etc. are all the ocean important boundaries of demarcating, and common feature is to form half envelope arc from set point according to certain distance and direction extrapolation.The present invention extrapolates that arc, extrapolation arc intersect, the cutting of extrapolation arc, deletes the automatic generation that the steps such as inner arc, Delete superfluous arc and the merging of extrapolation arc realize FOS+60M line, 350M line and 2500m+100M line by generation.

Beneficial effect: the present invention, in the mapping and ocean description process of China coast, has that surveying and mapping data is accurate, easy to operate and affected by natural weather more first-class.

The accompanying drawing explanation

Fig. 1 algorithm overall procedure of the present invention schematic diagram

Step (1) input variable schematic diagram in Fig. 2 example 2 of the present invention

In Fig. 3 example 2 of the present invention, step (2) generates extrapolation arc schematic diagram

In Fig. 4 example 2 of the present invention, step (3) extrapolation arc intersects schematic diagram

Step (4) cutting extrapolation arc schematic diagram in Fig. 5 example 2 of the present invention

In Fig. 6 example 2 of the present invention, step (5) is deleted the inner arc schematic diagram

In Fig. 7 example 2 of the present invention, step (6) Delete superfluous arc and step (7) merge extrapolation arc schematic diagram

Embodiment

Below in conjunction with enforcement of the present invention, illustrate:

Embodiment 1

Flow process as shown in Figure 1:

(1) input variable.Input initializaing variable by outside: extrapolation arc starting point G, extrapolated radius r and extrapolation direction D.Namely the extrapolate central coordinate of circle point of arc of extrapolation arc starting point, provided by the user, and the series fixed point forms primary data set G={g _i(x _i, y _i).Extrapolated radius is by extrapolation arc type decided, and FOS+60M line, 350M line and 2500m+100M line extrapolated radius r are respectively 60 nautical miles, 350 nautical miles and 100 nautical miles.Extrapolation direction D is the propagation direction of extrapolation arc, and FOS+60M line, 350M line and 2500m+100M line be offward extrapolation all.

(2) generate the extrapolation arc.

According to initial fixed point set G={g _i(x _i, y _i), with (x _i, y _i) put as the center of circle, with identical spherical distance r, generate sphere extrapolation arc a _i.Ergodic data set G={g _i(x _i, y _i), form data set Arc ₀={ a _i, each extrapolation arc a _iComprise central coordinate of circle (x _i, y _i) and extrapolated radius value r.

(3) the extrapolation arc intersects.

By dual searching loop data set Arc ₀={ a _i, obtain the mutual intersection point set A rc of extrapolation arc ₁={ a1 _i.Ergodic data collection Arc ₀={ a _i, take out successively extrapolation arc a _i.With extrapolation arc a _iHeaded by arc, enter second and recirculate, successively with data set Arc ₀={ a _iIn except a _iAll the other outer extrapolation arcs intersect judgement, and record intersection point forms data set Arc ₁={ a1 _i, each data point a1 _iComprise central coordinate of circle (x1 _i, y1 _i), extrapolated radius r and intersection point set C _i={ c _j.Arc ₁={ a1 _iData point a1 _iWith inherited data set Arc ₀={ a _iMiddle data point a _iAll properties.

Two extrapolation arc whether crossing bases for estimation are: the two distance of center circle is from d≤2r, and r is extrapolated radius, and when d<2r, the two has two intersection points, and when d=2r, the two has an intersection point; D > during 2r without intersection point.

Indefinite because of the intersection point sum of each extrapolation arc and other extrapolation arcs, therefore, need the data set C of design stores intersection point _i={ c _j.The all corresponding intersection point data set C of each extrapolation arc _i={ c _j, intersection point c _jShould comprise terrestrial coordinate value (x _j, y _j) and with the azimuth value α of extrapolation arc circle center line connecting _j.Intersection point c _jAccording to azimuth value α _jSize order deposit into data set C _i={ c _j.Intersection point data set C _i={ c _jChain enters data acquisition Arc ₁={ a1 _iThe new data acquisition Arc of rear formation ₂={ a2 _i.With Arc ₁={ a1 _iCompare, except intersection point C _i={ c _jOutside, Arc ₂={ a2 _iInherited Arc ₁={ a1 _iAll properties.

(4) cutting extrapolation arc.

Adopt dual circulation, cutting extrapolation arc forms little segmental arc.Ergodic data collection Arc ₂={ a2 _i, take out successively data point a2 _i, then travel through corresponding with it intersection point data acquisition C _i={ c _j.According to the intersection point order, take out successively two adjacent intersection points and form subarc section a3 _i, and form new data set AArc ₃={ a3 _i.Each data point a3 _iStarting point, terminal, intermediate point and the center of circle that all should comprise the subarc section, and the azimuth angle alpha of starting point and circle center line connecting ₁, terminal and circle center line connecting azimuth angle alpha ₂, and record subarc section a3 _iFrom original arc.So-called intermediate point, refer to and be positioned at subarc section a3 _iStarting point and the point between terminal, for judgement and the mutual relationship of original extrapolation arc.

(5) delete inner arc.

Adopt dual circulation to delete inner subarc section, namely by certain original arc, be completely contained in inner subarc section.Ergodic data set A rc ₀={ a _i, take out successively original arc data point a pointed _i.Each data point a _iEqual ergodic data collection Arc ₃={ a3 _iOnce, judge successively data point a3 _iWith a _iRelation, work as a3 _iSubarc section pointed is by original arc a _iInstitute comprises and not from original arc a _iThe time, this subarc section is inner arc, need be from data acquisition Arc ₃={ a3 _iThe middle deletion.Thereby at data set Arc ₃={ a3 _iForm new data set Arc on basis ₄={ a4 _i.A4 _iInherited Arc ₃={ a3 _iIn a3 _iAll properties.

Subarc section a3 _iWith data acquisition Arc ₀={ a _iMiddle a _iThe relation of inclusion basis for estimation is: group segmental arc a3 _iIntermediate point coordinate and a _iWhen central coordinate of circle is less than extrapolated radius r apart from d, subarc section a3 _iBe contained in a fully _i.

(6) Delete superfluous arc.

The extrapolation direction of FOS+60M line, 350M line and 2500m+100M line is offward, and is not the envelope arc of sealing, and in this and computer graphics, buffer zone has essential difference.Generally speaking, the extrapolation direction D of FOS+60M line, 350M line and 2500m+100M line is set as due east, due south, positive west or positive north.Extrapolation direction D is outside input variable.

Ergodic data collection Arc ₄={ a4 _i, by extrapolation direction D and subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂Relation judgement, α namely ₁Or α ₂Point to when consistent and retain this extrapolation segmental arc with extrapolation direction D, delete the segmental arc that has nothing to do, thereby form new data acquisition Arc ₅={ a5 _i, a5 _iInherited data set Arc ₄={ a4 _iMiddle a4 _iAll properties.Subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂By step (3), calculate.

Subarc section a4 _iAzimuth angle alpha ₁Or α ₂With the extrapolation direction, D concerns that basis for estimation is as follows:

(a) when extrapolation direction D is due east, and α ₁∈ | 0,180| or α ₂∈ | during 0,180|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(b) when extrapolation direction D is due south, and α ₁∈ | 90,270| or α ₂∈ | during 90,270|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(c) as extrapolation direction D when just western, and α ₁∈ | 180,360| or α ₂∈ | during 180,360|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(d) when extrapolation direction D is positive north, and α ₁∈ | 270,360| ∪ 0,90| or α ₂∈ 270,360| ∪ 0, and during 90|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain.

(7) the extrapolation arc merges.

Adopt dual searching loop data acquisition Arc ₅={ a5 _i, take out at first successively arc a5 _i, ergodic data set A rc again ₅={ a5 _i, take out except a5 successively _iOuter remainder data point, judge whether the two exists neighbouring relations.Whether the beginning or end coordinate that basis for estimation is the two is identical.Starting point and terminal point coordinate calculate by step (3).

By this step, set up all contacting directly of neighbouring relations segmental arc that have, can review all relevant segmental arcs by any segmental arc.Ergodic data set A rc again ₅={ a5 _i, by the mutual relationship between segmental arc, sort in order and form a complete data set Arc ₆={ a6 _i.Extrapolation arc Arc is linked in sequence ₆={ a6 _iForm final boundary L.

(8) the FOS+60M line generates.

The continent toe is counted according to collection FOS={fos _iProvided belong to system outside input variable the same as extrapolation direction D by the user.With FOS={fos _iBe initial data set G _fos, take 60 nautical miles as extrapolated radius r _fos, adopting step (2)～(7), the extrapolation direction D according to given, can obtain the FOS+60M line.

(9) the 350M line generates.

Territorial waters basic point data set BASE={b _iProvided belong to system outside input variable the same as extrapolation direction D by the user.With BASE={b _iBe initial data set G _b, take 350 nautical miles as extrapolated radius r _b, adopting step (2)～(7), the extrapolation direction D according to given, can obtain the 350M line.

(10) the 2500m+100M line generates.

2500m sea-bottom contour data set TOPO={t _iProvided belong to system outside input variable the same as extrapolation direction D by the user.With TOPO={t _iBe initial data set G _t, take 100 nautical miles as extrapolated radius r _t, adopting step (2)～(7), the extrapolation direction D according to given, can obtain the 2500m+100M line.

Embodiment 2

Fig. 2～Fig. 7 has showed from given continent toe point (G ₁～G ₄) generate the process of FOS+60M line:

(1) input variable.Input initializaing variable by outside: extrapolation arc starting point G, extrapolated radius r and extrapolation direction D.Primary data set G={g _i(x _i, y _i), in the sea, the extrapolation direction is due east to extrapolated radius r=60.

As shown in Figure 2.

(2) generate the extrapolation arc.

According to primary data set G={g _i(x _i, y _i), with g _i(x _i, y _i) put as the center of circle, take 60 nautical miles as radius r, generate sphere extrapolation arc a _i.Traversal initial data set G={g _i(x _i, y _i), form data set Arc ₀={ a _i, each extrapolation arc a _iComprise central coordinate of circle (x _i, y _i) and extrapolated radius value r.

As shown in Figure 3.

(3) the extrapolation arc intersects.

By dual searching loop data set Arc ₀={ a _i, obtain the mutual intersection point set A rc of extrapolation arc ₁={ a1 _i.Ergodic data collection Arc ₀={ a _i, take out successively extrapolation arc a _i.With extrapolation arc a _iHeaded by arc, enter second and recirculate, successively with data set Arc ₀={ a _iIn except a _iAll the other outer extrapolation arcs intersect judgement, and record intersection point forms new data set Arc ₁={ a1 _i, each data point a _iComprise central coordinate of circle, extrapolated radius and serial intersecting point coordinate value.

Each extrapolation arc a _iAn all corresponding intersection point data link table C _i={ c _j, intersection point c _jShould comprise the terrestrial coordinate value and with extrapolation arc a _iThe azimuth value α of circle center line connecting _j, intersection point c _iAccording to azimuth value α _jSize order deposit into data acquisition C _i={ c _j.Form new data acquisition Arc after inserting chained list ₂={ a _i.

As shown in Figure 4.

(4) cutting extrapolation arc.

Adopt dual circulation, cutting extrapolation arc forms little segmental arc.Traversal extrapolation arc data set Arc ₂={ a2 _i, take out successively data point a2 _i.Each extrapolation arc a2 _iEqual corresponding intersection point data acquisition C of traversal _i={ c _j.According to the intersection point order, take out successively adjacent intersection point and form the subarc section, and form new extrapolation segmental arc data acquisition Arc ₃={ a3 _i.Each subarc section a3 _iAll should comprise arc starting point coordinate, terminal point coordinate, intermediate point coordinate, central coordinate of circle, and the azimuth angle alpha of starting point and circle center line connecting ₁, the azimuth angle alpha of terminal and circle center line connecting ₂, and record subarc a3 _iFrom original arc a _i.As shown in Figure 5.

(5) delete inner arc.

Adopt dual circulation to delete the residue segmental arc, namely be comprised in inner segmental arc.Travel through original segmental arc data acquisition Arc ₀={ a _i, take out successively extrapolation arc a _i.Each extrapolation arc a _iAll travel through subarc segment data set Arc ₃={ a3 _iOnce, judge successively subarc section a3 _iWith each original arc a _iRelation, the group segmental arc is involved and not from original arc a _iThe time, subarc section a3 _iFor inner arc, need be from data acquisition Arc ₃={ a3 _iMiddle this subarc section of deleting.Thereby at subarc set A rc ₃={ a3 _iForm new data acquisition Arc on basis ₄={ a4 _i.

As shown in Figure 6.

(6) Delete superfluous arc.

Calculated the starting point azimuth angle alpha of each subarc by step (3) ₁With home orientation angle α ₂.By extrapolation direction D and azimuth angle alpha ₁And α ₂Relation judgement, α namely ₁Or α ₂Point to when consistent and retain this extrapolation segmental arc with extrapolation direction D, delete the segmental arc that has nothing to do, thereby form new data acquisition Arc ₅={ a5 _i.

(7) the extrapolation arc merges.

Adopt dual searching loop data acquisition Arc ₅={ a5 _i, take out at first successively arc a5 _i, ergodic data set A rc again ₅={ a5 _i, take out except a5 successively _iThe remainder data point, judge whether the two exists neighbouring relations.Whether the beginning or end coordinate that basis for estimation is the two is identical.Starting point and terminal point coordinate calculate by step (3).By this step, set up all contacting directly of neighbouring relations segmental arc that have, can review all relevant segmental arcs by any segmental arc.Ergodic data set A rc again ₅={ a5 _i, by the mutual relationship between segmental arc, sort in order and form a complete data set Arc ₆={ a6 _i.Extrapolation arc Arc is linked in sequence ₆={ a6 _iForm final boundary L, namely FOS+60M line.As shown in Figure 7.

Claims (1)

1. one kind 200 nautical miles with outer shelf extrapolation boundary generation method, it is characterized in that, comprises the following steps:

(1) generate the extrapolation arc

Given initial data set G={g _i(x _i, y _i) }: each data point g _i(x _i, y _i) all need by the given terrestrial coordinate (x of the mode of longitude and latitude _i, y _i); With (x _i, y _i) put as the center of circle, with identical spherical distance r, generate sphere extrapolation arc a _iTraversal initial data set G={g _i(x _i, y _i), form data set Arc ₀={ a _i, each extrapolation arc a _iComprise central coordinate of circle (x _i, y _i) and extrapolated radius value r;

(2) the extrapolation arc intersects

By dual searching loop data set Arc ₀={ a _i, obtain the mutual intersection point set A rc of extrapolation arc ₁={ a1 _i, ergodic data collection Arc ₀={ a _i, take out successively extrapolation arc a _i, with extrapolation arc a _iHeaded by arc, enter second and recirculate, successively with data set Arc ₀={ a _iIn except a _iAll the other outer extrapolation arcs intersect judgement, and record intersection point forms data set Arc ₁={ a1 _i, each data point a1 _iComprise central coordinate of circle, extrapolated radius and intersection point set C _i={ c _j; Arc ₁={ a1 _iData point a1 _iWith inherited data set Arc ₀={ a _iMiddle data point a _iAll properties;

Two extrapolation arc whether crossing bases for estimation are: the two distance of center circle is from d≤2r, and r is extrapolated radius, and when d<2r, the two has two intersection points, and when d=2r, the two has an intersection point; D > during 2r without intersection point;

Because of each extrapolation arc indefinite with the intersection point sum of other arcs of extrapolating, therefore, the data set C of design stores intersection point _i={ c _j.The all corresponding intersection point data set C of each extrapolation arc _i={ c _j, intersection point c _jComprise terrestrial coordinate value (x _j, y _j) and with the azimuth value α of extrapolation arc circle center line connecting _jIntersection point c _jAccording to azimuth value α _jSize order deposit into data set C _i={ c _j; Intersection point data set C _i={ c _jChain enters data acquisition Arc ₁={ a1 _iThe new data acquisition Arc of rear formation ₂={ a2 _i, with Arc ₁={ a1 _iCompare, except intersection point C _i={ c _jOutside, Arc ₂={ a2 _iInherited Arc ₁={ a1 _iAll properties;

(3) cutting extrapolation arc

Adopt dual circulation, cutting extrapolation arc forms little segmental arc; Ergodic data collection Arc ₂={ a2 _i, take out successively data point a2 _i, then travel through corresponding with it intersection point data acquisition C _i={ c _j, according to the intersection point order, take out successively two adjacent intersection points and form subarc section a3 _i, and form new data set Arc ₃={ a3 _i, each data point a3 _iStarting point, terminal, intermediate point and the center of circle that all should comprise the subarc section, and the azimuth angle alpha of starting point and circle center line connecting ₁, terminal and circle center line connecting azimuth angle alpha ₂, and record subarc section a3 _iFrom original arc; So-called intermediate point, refer to and be positioned at subarc section a3 _iStarting point and the point between terminal, for judgement and the mutual relationship of original extrapolation arc;

(4) delete inner arc

Adopt dual circulation to delete inner subarc section, namely by certain original arc, be completely contained in inner subarc section; Ergodic data set A rc ₀={ a _i, take out successively original arc data point a pointed _i, each data point a _iEqual ergodic data collection Arc ₃={ a3 _iOnce, judge successively data point a3 _iWith a _iRelation, work as a3 _iSubarc section pointed is by original arc a _iInstitute comprises and not from original arc a _iThe time, this subarc section is inner arc, from data acquisition Arc ₃={ a3 _iThe middle deletion, thereby at data set Arc ₃={ a3 _iForm new data set Arc on basis ₄={ a4 _i, a4 _iInherited Arc ₃={ a3 _iIn a3 _iAll properties;

Subarc section a3 _iWith data acquisition Arc ₀={ a _iMiddle a _iThe relation of inclusion basis for estimation is: group segmental arc a3 _iIntermediate point coordinate and a _iWhen central coordinate of circle is less than extrapolated radius r apart from d, subarc section a3 _iBe contained in a fully _i

(5) Delete superfluous arc

The extrapolation direction of FOS+60M line, 350M line and 2500m+100M line is offward, and is all nonocclusive envelope arc; The extrapolation direction D of FOS+60M line, 350M line and 2500m+100M line is set as due east, due south, positive west or positive north; Extrapolation direction D is outside input variable;

Ergodic data collection Arc ₄={ a4 _i, by extrapolation direction D and subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂Relation judgement, α namely ₁Or α ₂Point to when consistent and retain this extrapolation segmental arc with extrapolation direction D, delete the segmental arc that has nothing to do, thereby form new data acquisition Arc ₅={ a5 _i, a5 _iInherited data set Arc ₄={ a4 _iMiddle a4 _iAll properties; Subarc section a4 _iThe starting point azimuth angle alpha ₁With home orientation angle α ₂By step (3), calculate;

Subarc section a4 _iAzimuth angle alpha ₁Or α ₂With the extrapolation direction, D concerns that basis for estimation is as follows:

(a) when extrapolation direction D is due east, and α ₁∈ | 0,180| or α ₂∈ | during 0,180|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(b) when extrapolation direction D is due south, and α ₁∈ | 90,270| or α ₂∈ | during 90,270|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(c) as extrapolation direction D when just western, and α ₁∈ | 180,360| or α ₂∈ | during 180,360|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(d) when extrapolation direction D is positive north, and α ₁∈ | 270,360| ∪ | 0,90| or α ₂∈ | 270,360| ∪ | during 0,90|, subarc section a4 _iWith the extrapolation direction, D is consistent, should retain;

(6) the extrapolation arc merges

Adopt dual searching loop data acquisition Arc ₅={ a5 _i, take out at first successively subarc section a5 _i, ergodic data set A rc again ₅={ a5 _i, take out except a5 successively _iAll the other subarc sections, judge whether the two exists neighbouring relations; Whether the beginning or end coordinate that basis for estimation is the two is identical, and starting point and terminal point coordinate calculate by step (3);

By this step, set up all contacting directly of neighbouring relations segmental arc that have, can review all relevant segmental arcs, ergodic data set A rc again by any segmental arc ₅={ a5 _i, by the mutual relationship between segmental arc, sort in order and form a complete data set Arc ₆={ a6 _i, extrapolation arc Arc is linked in sequence ₆={ a6 _iForm final boundary L;

(7) the FOS+60M line generates

FOS={fos _iCount according to collection for the continent toe, with FOS={fos _iBe initial data set G _fos, take 60 nautical miles as extrapolated radius r _fos, adopt step (1)～(6), according to given extrapolation direction D, obtain the FOS+60M line;

(8) the 350M line generates

BASE={b _iBe territorial waters basic point data set, with BASE={b _iBe initial data set G _b, take 350 nautical miles as extrapolated radius r _b, adopt step (1)～(6), according to given extrapolation direction D, obtain the 350M line;

(9) the 2500m+100M line generates

TOPO={t _iBe 2500m sea-bottom contour data set, with TOPO={t _iBe initial data set G _t, take 100 nautical miles as extrapolated radius r _t, adopt step (1)～(6), according to given extrapolation direction D, obtain the 2500m+100M line.

Images