专利名称：一种200海里以外大陆架外推界限生成方法
技术领域：
本发明涉及一种200海里以外大陆架界限生成方法，具体是指一种基于已知定点，向特定方向外推一定距离，从而生成F0S+60M线、350M线和2500m + 100M线等200海里以外大陆架划界所必须的技术方法，涉及到计算机图形学、海底科学和海洋划界领域。
背景技术：
根据1982年签署生效的《海洋法公约》，“沿海国的大陆架包括其领海以外依其陆地领土的全部自然延伸，扩展到大陆边外缘的海底区域的海床和底土，如果从测算领海宽度的基线量起到大陆边的外缘的距离不到二百海里，则扩展到二百海里的距离”；沿海国如果主张从测算领海宽度的基线量起超过200海里，则应按照《公约》第七十六条和附件二第四条的相关要求划定200海里以外大陆架外部界限，并向大陆架界限委员会提交划界案。F0S+60M线、350M线和2500m+100M线(M是海里的简称，是200海里以外大陆架划界案编制中的国际习惯)是确定200海里以外大陆架外部界限的重要基础。具体而言，F0S+60M线是大陆坡脚点(简称F0S)向海方向外推60海里所生成的界限，350M线是从领海基点向海方向外推350海里所形成的界限，2500m+100M线是2500m等深线向海方向外推100海里所生成的界限。F0S+60M线简称公式线，350M线和2500m + 100M线又简称限制线。三条界限的共同特征是从所给定的定点向指定方向外推一定距离，所以被统称为外推界限。从检索的公开资料分析，目前国内无成熟的技术方法用于200海里以外大陆架F0S+60M线、350M线和2500m+100M线等划界界限的生成，尽管有少数文献对200海里以外大陆架划界进行了介绍，但无翔实的界限生成方法之内容，现有技术方法不足以支持F0S+60M线、350M线和2500m+100M线等划界界限的生成。
发明内容
本发明针对现有技术中的不足，提出一种利用给定点为圆心，按照一定距离为半径在地球面生成系列外推弧，通过外推弧的相交和合并运算，从而生成指定界限的技术方法。具体而言是从给定大陆坡脚点外推60海里生成F0S+60M线、从给定领海基点外推350海里生成350海里限制线、从给定2500m等深线点外推100海里生成2500m+100M线的技术方法。该发明可应用于自动生成200海里以外大陆架划界所必须的F0S+60M线、350M线和2500m+100M线等界限点，可应用于海洋划界领域其他种类外推线的生成，也可应用于计算机图形学等领域外推包络线的生成。本发明是通过下述技术方案得以实现的一种200海里以外大陆架外推界限生成方法，其特征在于，包括下列步骤(I)生成外推弧。给定初始数据集G= {gi (Xi, Yi)},每个数据点gi (Xi, Yi)均需按经纬度的方式给定大地坐标(Xi，yi)。以(Xi，yi)点为圆心，以相同的球面距离r，生成球面外推弧％。遍历初始数据集G=IgiUi, Yi) I，形成数据集Arctl=IaJ ,每个外推弧Si包括圆心坐标(Xi, Yi)和外推半径值r。(2)外推弧相交。通过双重循环遍历数据集Arctl=IaJ ,获取外推弧相互间的交点集合Arc1=IalJ。遍历数据集Arctl= {aj，依次取出外推弧％。以外推弧％为首弧，进入第二重循环，依次和数据集Arctl=IaJ中除Si外的其余外推弧进行相交判断,并记录交点形成数据集Arc1= {alj ,每个数据点ah包括圆心坐标，外推半径和交点集合Ci=IcJ。Arc1= {alj数据点Bli与继承了数据集Arc。= {aj中数据点Bi的所有属性。二个外推弧是否相交的判断依据为二者圆心距离d < 2r，r为外推半径，当d〈2r时，二者有两个交点，当d=2r时，二者有一个交点；d>2r时无交点。因每个外推弧与其他外推弧的交点总数不定，因此，需要设计存储交点的数据集Ci= {Cj}。每个外推弧均对应一个交点数据集Ci= {Cj},交点Cj应包括大地坐标值(Xj, Yj)及其与外推弧圆心连线的方位角值α”交点h按照方位角值的大小顺序存放入数据集Ci= {c」}。交点数据集Ci= {c」}链入数据集合Arc1=IalJ后形成新的数据集合Arc2= {a2j。与Arc1=IalJ相比，除交点Ci= {c」}之外，Arc2= {a2j继承了 Arc1=IalJ所有属性。(3)切割外推弧。采用双重循环，切割外推弧形成小的弧段。遍历数据集Arc2= {a2j，依次取出数据点a2i，然后遍历与之对应的交点数据集合(；={(^}。按照交点顺序，依次取出相邻的两个交点形成子弧段a3i;并形成新的数据集Arc3= {a3j。每个数据点a3i均应包括子弧段的起点、终点、中间点和圆心，以及起点与圆心连线的方位角a i、终点与圆心连线的方位角α2，且记录子弧段a3i来自的原始弧。所谓中间点，是指位于子弧段a3i的起点和终点之间的点，用于判断与原始外推弧的相互关系。(4)删除内部弧。采用双重循环删除内部子弧段，也就是被某个原始弧完全包含在内部的子弧段。遍历数据集合Arctl= {aj，依次取出原始弧所指向的数据点％。每个数据点％均遍历数据集Arc3= {a3j 一次,依次判断数据点a3i与Si的关系，当a3i所指向的子弧段被原始弧Si所包含且不来自原始弧％时，该子弧段为内部弧，需从数据集合ArC3={a3j中删除。从而在数据集Arc3={a3i}基础上形成新的数据集Arc4= {a4j。a4jS承了 Arc3={a3j中的a3i的所有属性。子弧段a3i与数据集合Arc。= {aj中％包含关系判断依据为当子弧段a3i的中间点坐标与％圆心坐标距离d小于外推半径r时，子弧段a3i完全包含于％。·(5)删除多余弧。F0S+60M线、350M线和2500m+100M线的外推方向均为向海方向，且都不是封闭的包络弧，这与计算机图形学中缓冲区有本质的区别。一般而言，F0S+60M线、350M线和2500m+100M线的外推方向D设定为正东、正南、正西或正北。外推方向D是外部输入变量。遍历数据集Arc4= {a4j ,通过外推方向D与子弧段a4i起点方位角a j和终点方位角α2的关系判断，也就是\或Ci2与外推方向D指向一致时保留该外推弧段，删除无关弧段，从而形成新的数据集合Arc5={a5i}，a5i继承了数据集Arc4= {a4j中a4i的所有属性。子弧段a4i起点方位角α !和终点方位角α 2已通过步骤(3)计算。
子弧段a4i方位角a i或α 2与外推方向D关系判断依据如下(a)当外推方向D为正东时，且ai e O, 180或者a 2 e O, 180时，子弧段
与外推方向D —致，应当保留；(b)当外推方向D为正南时，且ai e I 90，270 I或者a2 e 90, 270时，子弧段
al与外推方向D —致，应当保留；(C)当外推方向D为正西时，且CIi e I 180，360 I或者a2 e 180, 360时，子弧
段a4i与外推方向D—致，应当保留；(d)当外推方向D为正北时，且Ci1 e 1270，3601 U 0,90或者 a2e I 270，360 I U 0，90时，子弧段与外推方向D—致，应当保留。(6)外推弧合并。采用双重循环遍历数据集合Arc5= {a5j，首先依次取出子弧段a5i，再次遍历数据集合Arc^iaSj，依次取出除a5i的其余子弧段，判断二者是否存在相邻关系。判断依据为二者的起点或终点坐标是否相同。起点和终点坐标已通过步骤(3)计算。通过该步骤，已建立了所有有相邻关系弧段的直接联系，通过任意弧段可以追溯所有相关弧段。再次遍历数据集合Arc5= {a5j，通过弧段间的相互关系，按顺序进行排序形成一个完整的数据集Arc6= {a6j。顺序连接外推弧Arc6= {a6j形成最终的界限L。(7)F0S+60M 线生成FOS= {fosj为大陆坡脚点数据集，以FOS= {f0Si}为初始数据集Gf()S，以60海里为外推半径rf()S，采用步骤(I) (6)，按照给定的外推方向D，获得F0S+60M线；(8) 35OM 线生成BASE= IbJ为领海基点数据集，以BASE= {bj为初始数据集Gb，以350海里为外推半径rb，采用步骤(I) (6)，按照给定的外推方向D，获得350M线；(9)2500m+100M 线生成TOPO= {tj为2500m等深线数据集，以TOPO= {tj为初始数据集Gt，以100海里为外推半径rt，采用步骤(I) (6)，按照给定的外推方向D，获得2500m+100M线。根据联合国《海洋法公约》第七十六条及大陆架划界委员会技术准则之规定，F0S+60M线、350M线和2500m+100M线等均是海洋划界重要界限，共同特点是从给定点按照一定距离和方向外推形成半包络弧。本发明通过生成外推弧、外推弧相交、外推弧切割、删除内部弧、删除多余弧和外推弧合并等步骤实现F0S+60M线、350M线和2500m+100M线的自动生成。有益效果本发明在我国沿海的测绘和海洋划界过程中，具有测绘数据准确、操作方便、以及受自然天气影响较上等。
图I本发明的算法总体流程示意2本发明实例2中步骤(I)输入变量示意3本发明实例2中步骤(2)生成外推弧示意4本发明实例2中步骤(3)外推弧相交示意5本发明实例2中步骤(4)切割外推弧示意图
图6本发明实例2中步骤(5)删除内部弧示意7本发明实例2中步骤(6)删除多余弧和步骤(7)合并外推弧示意图
具体实施例方式下面结合本发明的实施作具体说明
实施例I如图I所示的流程(I)输入变量。由外部输入初始变量外推弧起始点G、外推半径r和外推方向D。外推弧起点也就是外推弧的圆心坐标点，由用户给出，系列定点形成初始数据集合G=Igi(XilYi)Io外推半径由外推弧类型决定，F0S+60M线、350M线和2500m+100M线外推半径r分别为60海里、350海里和100海里。外推方向D是外推弧的扩展方向，F0S+60M线、350M线和2500m+100M线均向海方向外推。(2)生成外推弧。根据初始的定点集合G={gi(Xi，yi)}，以(Xi，yi)点为圆心，以相同的球面距离r，生成球面外推弧Si。遍历数据集合G= {gi (Xi, Yi)},形成数据集Arctl= {aj ,每个外推弧Si包括圆心坐标(X^yi)和外推半径值r。(3)外推弧相交。通过双重循环遍历数据集Arcci=IaJ ,获取外推弧相互间的交点集合Arc1=IalJ。遍历数据集Arctl= {aj，依次取出外推弧％。以外推弧％为首弧，进入第二重循环，依次和数据集Arctl=IaJ中除Si外的其余外推弧进行相交判断,并记录交点形成数据集Arc1= {alj ,每个数据点ali包括圆心坐标(Xli, yl),外推半径r和交点集合Ci=IcjI。Arc1= {alj数据点ah与继承了数据集Arctl= {aj中数据点％的所有属性。二个外推弧是否相交的判断依据为二者圆心距离d < 2r，r为外推半径，当d〈2r时，二者有两个交点，当d=2r时，二者有一个交点；d>2r时无交点。因每个外推弧与其他外推弧的交点总数不定，因此，需要设计存储交点的数据集Ci= {cj}。每个外推弧均对应一个交点数据集Ci= {C」},交点Cj应包括大地坐标值(χ」，Yj)及其与外推弧圆心连线的方位角值α”交点h按照方位角值的大小顺序存放入数据集Ci= {c」}。交点数据集Ci= {c」}链入数据集合Arc1=IalJ后形成新的数据集合Arc2= {a2j。与Arc1=IalJ相比，除交点Ci= {c」}之外，Arc2= {a2j继承了 Arc1=IalJ所有属性。(4)切割外推弧。采用双重循环，切割外推弧形成小的弧段。遍历数据集Arc2= {a2j，依次取出数据点a2i，然后遍历与之对应的交点数据集合(；={(^}。按照交点顺序，依次取出相邻的两个交点形成子弧段a3i;并形成新的数据集AArc3= {a3j。每个数据点a3i均应包括子弧段的起点、终点、中间点和圆心，以及起点与圆心连线的方位角a i、终点与圆心连线的方位角α2，且记录子弧段&31来自的原始弧。所谓中间点，是指位于子弧段a3i的起点和终点之间的点，用于判断与原始外推弧的相互关系。(5)删除内部弧。采用双重循环删除内部子弧段，也就是被某个原始弧完全包含在内部的子弧段。遍历数据集合Arctl= {aj，依次取出原始弧所指向的数据点％。每个数据点％均遍历数据集Arc3= {a3j 一次,依次判断数据点a3i与Si的关系，当a3i所指向的子弧段被原始弧Si所包含且不来自原始弧％时，该子弧段为内部弧，需从数据集合ArC3={a3j中删除。从而在数据集Arc3={a3i}基础上形成新的数据集Arc4= {a4j。a4jS承了 Arc3={a3j中的a3i的所有属性。子弧段a3i与数据集合Arctl= {aj中％包含关系判断依据为当子弧段a3i中间点坐标与％圆心坐标距离d小于外推半径r时，子弧段a3i完全包含于％。(6)删除多余弧。 F0S+60M线、350M线和2500m+100M线的外推方向均为向海方向，且都不是封闭的包络弧，这与计算机图形学中缓冲区有本质的区别。一般而言，F0S+60M线、350M线和2500m+100M线的外推方向D设定为正东、正南、正西或正北。外推方向D是外部输入变量。遍历数据集Arc4= {a4j ,通过外推方向D与子弧段a4i起点方位角a j和终点方位 角α2的关系判断，也就是\或Ci2与外推方向D指向一致时保留该外推弧段，删除无关弧段，从而形成新的数据集合Arc5={a5i}，a5i继承了数据集Arc4= {a4j中a4i的所有属性。子弧段a4i起点方位角α !和终点方位角α 2已通过步骤(3)计算。子弧段a4i方位角a i或α 2与外推方向D关系判断依据如下(a)当外推方向D为正东时，且α e O, 180或者a2 e O, 180时，子弧段与外推方向D —致，应当保留；(b)当外推方向D为正南时，且ai e I 90，270 I或者a2 e 90, 270时，子弧段
al与外推方向D —致，应当保留；(C)当外推方向D为正西时，且CIi e I 180，360 I或者a2 e 180, 360时，子弧
段a4i与外推方向D—致，应当保留；(d)当外推方向D为正北时，且a i e 1270，3601 U O, 901或者α2 e 270，360 I U 0,90时，子弧段与外推方向D—致，应当保留。(7)外推弧合并。采用双重循环遍历数据集合Arc5= {a5j，首先依次取出弧a5i，再次遍历数据集合Arc5HaSJ，依次取出除&51外的其余数据点，判断二者是否存在相邻关系。判断依据为二者的起点或终点坐标是否相同。起点和终点坐标已通过步骤(3)计算。通过该步骤，已建立了所有有相邻关系弧段的直接联系，通过任意弧段可以追溯所有相关弧段。再次遍历数据集合Arc5= {a5j，通过弧段间的相互关系，按顺序进行排序形成一个完整的数据集Arc6= {a6j。顺序连接外推弧Arc6= {a6j形成最终的界限L。(8)F0S+60M 线生成。大陆坡脚点数据集FOS= {f0Si}由用户给出，和外推方向D —样属于系统外部输入变量。以F0S={fOSi}为初始数据集Gfos，以60海里为外推半径rf()S，采用步骤(2) (7)，按照给定的外推方向D，可以获得F0S+60M线。(9) 35OM 线生成。领海基点数据集BASE= IbJ由用户给出，和外推方向D —样属于系统外部输入变量。以BASE=IbJ为初始数据集Gb，以350海里为外推半径rb，采用步骤(2) (7)，按照给定的外推方向D，可以获得350M线。(10)2500m+100M 线生成。
2500m等深线数据集TOPO=ItJ由用户给出，和外推方向D—样属于系统外部输入变量。以TOPO=ItJ为初始数据集Gt，以100海里为外推半径!·,，采用步骤(2) (7)，按照给定的外推方向D，可以获得2500m+100M线。实施例2图2 图7展示 了从给定大陆坡脚点(G1 G4)生成F0S+60M线的过程(I)输入变量。由外部输入初始变量外推弧起始点G、外推半径r和外推方向D。初始数据集合G= {gi (Xi, Yi)}，外推半径r=60海里，外推方向为正东。如图2所示。(2)生成外推弧。根据初始数据集合G= {gi (Xi, yi)}，以gi (Xi, Yi)点为圆心，以60海里为半径r，生成球面外推弧Si。遍历初始数据集G= {gi (Xi, Yi)},形成数据集Arctl= {aj ,每个外推弧Si包括圆心坐标(X^yi)和外推半径值r。如图3所示。(3)外推弧相交。通过双重循环遍历数据集Arcci=IaJ ,获取外推弧相互间的交点集合Arc1=IalJ。遍历数据集Arctl= {aj，依次取出外推弧％。以外推弧％为首弧，进入第二重循环，依次和数据集Arctl= {aj中除％外的其余外推弧进行相交判断，并记录交点形成新的数据集Arc1=IalJ ,每个数据点Si包括圆心坐标,外推半径和系列交点坐标值。每个外推弧ai均对应一个交点数据链表Ci= {cj},交点Cj应包括大地坐标值及其与外推弧％圆心连线的方位角值a j,交点Ci按照方位角值a j的大小顺序存放入数据集合Ci=Icj。插入链表后形成新的数据集合Arc2=IaJ。如图4所示。(4)切割外推弧。采用双重循环，切割外推弧形成小的弧段。遍历外推弧数据集Arc2= {a2j，依次取出数据点a2it)每个外推弧a2i均遍历之对应的交点数据集合Ci=Icj。按照交点顺序，依次取出相邻交点形成子弧段，并形成新的外推弧段数据集合Arc3= {a3j。每个子弧段a3i均应包括弧起点坐标、终点坐标、中间点坐标、圆心坐标，以及起点与圆心连线的方位角α 1；终点与圆心连线的方位角α 2，且记录子弧a3i来自的原始弧a”如图5所示。(5)删除内部弧。采用双重循环删除剩余弧段，也就是被包含在内部的弧段。遍历原始弧段数据集合Arctl=IaJ ,依次取出外推弧a”每个外推弧Si均遍历子弧段数据集合Arc3= {a3j —次，依次判断子弧段a3i与每个原始弧％的关系，当子弧段被包含且不来自原始弧％时，子弧段a3iS内部弧，需从数据集合ArC3={a3j中删除该子弧段。从而在子弧集合Arc3= {a3j基础上形成新的数据集合Arc4= {a4j。如图6所示。(6)删除多余弧。通过步骤(3)已计算了每个子弧的起点方位角Ci1和终点方位角α2。通过外推方向D与方位角\和α2的关系判断，也就是\或α2与外推方向D指向一致时保留该外推弧段，删除无关弧段，从而形成新的数据集合Arc5= {a5j。
(7)外推弧合并。
采用双重循环遍历数据集合Arc5= {a5j，首先依次取出弧a5i，再次遍历数据集合Arc5= {a5j，依次取出除a5i的其余数据点，判断二者是否存在相邻关系。判断依据为二者的起点或终点坐标是否相同。起点和终点坐标已通过步骤(3)计算。通过该步骤，已建立了所有有相邻关系弧段的直接联系，通过任意弧段可以追溯所有相关弧段。再次遍历数据集合Arc5= {a5j，通过弧段间的相互关系，按顺序进行排序形成一个完整的数据集Arc6= {a6i}。顺序连接外推弧Arc6= {a6j形成最终的界限L，也就是F0S+60M线。如图7所示。
权利要求
1.一种200海里以外大陆架外推界限生成方法，其特征在于，包括下列步骤 (1)生成外推弧 给定初始数据集G=IgiUi, yj}:每个数据点gi (Xi, y)均需按经纬度的方式给定大地坐标(Xpyi);以(Xpyi)点为圆心，以相同的球面距离r，生成球面外推弧％ ;遍历初始数据集G=IgiUi, Yi) I，形成数据集Arctl=IaJ ,每个外推弧Si包括圆心坐标(Xi, Yi)和外推半径值r ； (2)外推弧相交 通过双重循环遍历数据集Arctl=IaJ ,获取外推弧相互间的交点集合Arc1=IalJ ,遍历数据集Arctl= {aj，依次取出外推弧ai;以外推弧％为首弧，进入第二重循环，依次和数据集Arc0= {aj中除Si外的其余外推弧进行相交判断,并记录交点形成数据集Arc1=IalJ ,每个数据点ah包括圆心坐标，外推半径和交点集合Ci=IcjI ;ArCl={alJ数据点Bli与继承了数据集ArcfiaJ中数据点％的所有属性； 二个外推弧是否相交的判断依据为二者圆心距离d < 2r，r为外推半径，当d〈2r时，二者有两个交点，当d=2r时,二者有一个交点；d>2r时无交点； 因每个外推弧与其他外推弧的交点总数不定，因此，设计存储交点的数据集Ci= Icj。每个外推弧均对应一个交点数据集Ci=Ic^,交点包括大地坐标值(\，yp及其与外推弧圆心连线的方位角值a j ;交点按照方位角值a j的大小顺序存放入数据集Ci=Icj ;交点数据集Ci= {C」}链入数据集合Arc1=IalJ后形成新的数据集合Arc2= {a2j ,与Arc1= {alj相t匕，除交点Ci=IcjI之外，Arc2= {a2j继承了 Arc1=IalJ所有属性； (3)切割外推弧 采用双重循环，切割外推弧形成小的弧段；遍历数据集Arc2= {a2j，依次取出数据点a2i，然后遍历与之对应的交点数据集合Ci=Ic^,按照交点顺序，依次取出相邻的两个交点形成子弧段a3i;并形成新的数据集Arc3= {a3j，每个数据点a3i均应包括子弧段的起点、终点、中间点和圆心，以及起点与圆心连线的方位角a i、终点与圆心连线的方位角α 2，且记录子弧段a3i来自的原始弧；所谓中间点，是指位于子弧段a3i的起点和终点之间的点，用于判断与原始外推弧的相互关系； (4)删除内部弧 采用双重循环删除内部子弧段，也就是被某个原始弧完全包含在内部的子弧段；遍历数据集合Arctl= {aJ，依次取出原始弧所指向的数据点ai;每个数据点％均遍历数据集Arc3= {a3j 一次,依次判断数据点a3i与Si的关系，当a3i所指向的子弧段被原始弧Si所包含且不来自原始弧％时，该子弧段为内部弧，从数据集合Arc3= {a3j中删除，从而在数据集Arc3= {a3j基础上形成新的数据集Arc4= {a4j ,继承了 Arc3= {a3j中的a3i的所有属性； 子弧段a3i与数据集合Arctl= {aj中％包含关系判断依据为当子弧段a3i的中间点坐标与％圆心坐标距离d小于外推半径r时，子弧段a3i完全包含于％ ； (5)删除多余弧 F0S+60M线、350M线和2500m+100M线的外推方向均为向海方向，且都为非封闭的包络弧；F0S+60M线、350M线和2500m+100M线的外推方向D设定为正东、正南、正西或正北；外推方向D是外部输入变量； 遍历数据集Arc4= {a4j ,通过外推方向D与子弧段a4i起点方位角Ct1和终点方位角α 2的关系判断，也就是a i或α 2与外推方向D指向一致时保留该外推弧段，删除无关弧段，从而形成新的数据集合Arc5=MJ，a5i继承了数据集Arc4= {a4j中a4i的所有属性；子弧段a4i起点方位角α 和终点方位角α 2已通过步骤(3)计算； 子弧段al方位角a i或α 2与外推方向D关系判断依据如下 (a)当外推方向D为正东时，且Ci1e O, 180或者a 2 e O, 180时，子弧段8七与外推方向D—致，应当保留； (b)当外推方向D为正南时，且CI1e I 90，270 I或者a 2 e 90, 270时，子弧段a4i与外推方向D —致，应当保留； (C)当外推方向D为正西时，且ai e I 180，360 I或者a 2 e 180, 360时，子弧段与外推方向D —致，应当保留； (d)当外推方向D为正北时，且ai e 1270，3601 U 0,90或者a2e I 270，360 I U 0，90时，子弧段与外推方向D—致，应当保留； (6)外推弧合并 采用双重循环遍历数据集合Arc5= {a5j，首先依次取出子弧段a5i，再次遍历数据集合Arc5= {a5j，依次取出除a5i的其余子弧段，判断二者是否存在相邻关系；判断依据为二者的起点或终点坐标是否相同，起点和终点坐标已通过步骤(3)计算； 通过该步骤，已建立了所有有相邻关系弧段的直接联系，通过任意弧段可以追溯所有相关弧段，再次遍历数据集合Arc5= {a5j，通过弧段间的相互关系，按顺序进行排序形成一个完整的数据集Arc6= {a6j，顺序连接外推弧Arc6= {a6j形成最终的界限L ； (7)F0S+60M线生成 FOS= {fosj为大陆坡脚点数据集，以FOS= {fosj为初始数据集Gf()S，以60海里为外推半径rf()S，采用步骤(I) (6)，按照给定的外推方向D，获得F0S+60M线； (8)350M线生成 BASE=IbJ为领海基点数据集，以BASE=IbJ为初始数据集Gb，以350海里为外推半径rb，采用步骤(I) (6)，按照给定的外推方向D，获得350M线； (9)2500m+100M线生成 TOPO= {tj为2500m等深线数据集，以TOPO=ItJ为初始数据集Gt，以100海里为外推半径rt，采用步骤(I) (6)，按照给定的外推方向D，获得2500m+100M线。
全文摘要
本发明公开了大陆架界限生成方法，具体是指一种基于已知定点，向特定方向外推一定距离的大陆架界限生成方法。本发明是通过若干步骤包括生成外推弧、外推弧相交、切割外推弧、删除内部弧、删除多余弧、外推弧合并等步骤实现FOS+60M线、350M线和2500m+100M线的自动生成。本发明的优点是在我国沿海的测绘和海洋划界过程中，具有测绘数据准确、操作方便、以及受自然天气影响较上等。
文档编号G06F19/00GK102938029SQ20121053672
公开日2013年2月20日 申请日期2012年12月12日 优先权日2012年12月12日
发明者吴自银, 李家彪, 方银霞, 黎明碧, 李守军, 尚继宏 申请人:国家海洋局第二海洋研究所