CN110889221B - 基于非欧保形变换的船舶避碰模型的建立方法 - Google Patents
基于非欧保形变换的船舶避碰模型的建立方法 Download PDFInfo
- Publication number
- CN110889221B CN110889221B CN201911155257.3A CN201911155257A CN110889221B CN 110889221 B CN110889221 B CN 110889221B CN 201911155257 A CN201911155257 A CN 201911155257A CN 110889221 B CN110889221 B CN 110889221B
- Authority
- CN
- China
- Prior art keywords
- ship
- riemann
- speed
- relative
- projection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明提出一种基于非欧保形变换的船舶避碰模型的建立方法,将值班驾驶员的眼球作为真实情况下瞭望观测的研究基础,利用黎曼球面理论在本船舶上建立球面模型,将本船所在平面内的所有信息都利用空间投影的方式在球面上找到一一对应的投影点,利用此方法可将船舶的运动轨迹与球面上的圆弧相对应,更加符合人体视觉对风险的感知。在三维空间局面上对船舶的碰撞风险研究能有效反映真实情况下来船的相对速度与相对距离,有利于船舶值班驾驶员对船舶碰撞风险进行有效感知。本发明利用投影变换进行非欧空间的球面映射进行建模,提供更加符合风险感知的两船会遇相对速度和相对距离的计算方法,为值班驾驶员进行正确的避碰决策提供技术支持。
Description
技术领域
本发明属于船舶智能避碰决策的技术领域,尤其涉及一种基于非欧保形变换的船舶避碰模型的建立方法。
背景技术
随着我国航运业的迅猛发展,船舶运输贸易量显著增加。在进行国际货物运输时,船舶之间碰撞的概率将不断增大。但即使船舶智能化时代即将到来,船舶碰撞事故风险在相当长一段时间内仍无法有效降低。提出有效的船舶碰撞事故风险预测及防控机制具有紧迫性与巨大的商业价值。
现如今,船舶会遇时进行避碰决策主要依照于《1972国际海上避碰规则》(以下简称规则),《规则》中大部分的对于会遇局面的判断建立于平面的分析判断。基于平面对会遇局面的分析存在理论缺陷,主要表现为:受二维平面模型自身维度和平面的限制,对于事物的观察结果必然会出现近大远小的视觉误差。对于船舶会遇相对速度的判断于真实情况下值班驾驶员观测的来船运动情况不相吻合。值班驾驶员在观测的来船的速度应是不断变大的过程,但是平面理论分析结果是相对速度不断减小。为解决船舶会遇时平面的理论分析结果与实际的观测结果的出入,可引入视觉投影来解释相关现象。在真实的船舶航行时,值班驾驶员瞭望观测到的场景应是将船舶所在周围环境中的事物经值班驾驶员眼球投影至其视网膜上,值班驾驶员利用视网膜上外部环境投影的信息做出相应的避碰决策。
发明内容
本发明所要解决的技术问题在于针对上述存在的问题,提供一种基于非欧保形变换的船舶避碰模型的建立方法,将本船与来船的相对速度与相对距离利用黎曼球面进行视觉投影,将实时的船舶碰撞风险进行可视化呈现。
本发明解决上述技术问题所采用的技术方案是:基于非欧保形变换的船舶避碰模型的建立方法,其特征在于,包括如下步骤:
S1)数据提取:利用船舶的VDR进行VDR数据的提取;
S2)速度叠加:运用速度的合成原理,将本船的速度叠加至来船上,并运用速度的合成原理,将叠加的本船速度与他船速度进行矢量合成,本船认为处于静止状态;
S3)建立黎曼球面:以本船为空间坐标原点,建立空间直角坐标系,以本船作为黎曼球面的最低点,在本船上方建立黎曼球面;
S4)来船进行球面投影:用直线将黎曼球面的最高点与来船相连,直线于黎曼球面相交于除黎曼球面的最高点的另外一点,来船在黎曼球面进行了球面投影;
S5)相对距离与相对速度投影点:运用黎曼球面将本船平面内的一切物体投影至黎曼球面上,并在黎曼球面上可以找到对应的投影点,在黎曼球面上将来船相对于本船的相对距离与相对速度通过投影点进行反映,在黎曼球面上将来船在平面内的运动转化为圆面的圆弧运动;
S6)求解投影点的坐标:利用几何关系将黎曼球面与来船在平面运动的直线运用方程进行表述,建立黎曼球面最高点与来船位置的直线方程,利用来船的位置点与黎曼球面最高点所连的直线方程与黎曼球面的方程进行联立求解出投影点的坐标;
S7)求取相对距离与相对速度:利用投影点的坐标求取来船与本船之间的相对距离,将投影点坐标进行求导,计算出在投影点处的速度,将该速度投影至投影点与本船连线的直线上,求取出投影点相对于本船的速度,利用投影点相对于本船的距离和相对距离来反映在黎曼球面下本船观测到他船的真实运动情况;
S8)建立船舶碰撞风险模型:根据相对速度与相对距离描述船舶之间的碰撞风险,船舶碰撞的风险是与相对速度成正相关与相对距离成反比,所以建立风险评定方程以相对速度与相对距离的倒数相乘积表针风险。
按上述方案,步骤S1)中所述VDR数据包括船舶航速、航向、主机参数、AIS数据、雷达图像数据和机舱驾驶台的语音数据。
按上述方案,步骤S4)中在原点、黎曼球面最高点、来船的位置坐标和投影点建立三角形,利用几何三角相似关系可以证明出一对相似三角形关系,利用保形变换思想,来船在二维平面下的直线运动利用黎曼球面投影在球面显示出圆弧曲线运动,可将来船的运动信息通过投影点进行反映。
本发明的有益效果是:基于非欧保形变换的船舶避碰风险评估模型与现有基于欧几里得空间分析得出的船舶避碰风险评估理论模型相比具有以下优点:在非欧空间模型下,来船在黎曼球面上的投影点与本船之间的相对速度呈现先增大后减小,在碰撞点之后一段时间内继续增大,增加到最高点之后迅速的减小,matlab仿真的结果符合现实中值班驾驶员瞭望到对面来船向本船运动的真实运动情况;
在非欧空间模型下,来船在黎曼球面上的投影点与本船之间的相对距离的变化规律在碰撞时间点之前是呈减小趋势的,随之时间的不断增大,两船之间的距离不断地减小,符合现实中两船距离的实际变化,同时也符合值班驾驶员所看到的真实感受;在非欧空间的模型下,利用来船在黎曼球面上的投影点与本船之间的相对速度与相对距离拟合出的船舶碰撞风险系数随时间的不断增大,在碰撞时间点之前的一段时间内是呈现不断增大的趋势,matlab仿真的结果符合现实中驾驶员的真实感受,即当两船不断靠近时,船舶的风险是不断增大的,符合两船会遇的碰撞风险态势的变化情况。
附图说明
图1为本发明一个实施例的在平面下船舶避碰参数间的关系示意图。
图2为本发明一个实施例的在平面下两船之间相对速度、相对距离与时间关系示意图。
图3为本发明一个实施例的在基于Mobius变换的空间船舶运动情况示意图。
图4为本发明一个实施例的在黎曼球面中两船之间相对速度、相对距离与时间的关系示意图。
图5为本发明一个实施例的在黎曼球面中两船之间碰撞风险与时间的关系示意图。
具体实施方式
为更好地理解本发明,下面结合附图和实施例对本发明进一步的描述。
1.平面下船舶会遇时船舶之间相对速度与相对距离之间的关系
如图1所示,搭建出欧几里得空间下的船舶会遇模型,若某海域有两条处于运动中的船舶A、B,假设将船舶B的速度叠加至船舶A,则船舶B静止于点O,船舶A以恒定的速度V向船舶B的方向运动。设船舶B的船舶邻域是以半径为r的圆。以船舶A运动方向做一条交于圆的一条直线L。圆心O与直线L的垂直距离d(指DCPA)。t=0时两船的距离之差为S0,船舶A在运动过程中与船舶B之间的动态距离为S。
在船舶A向船舶B运动时两船之间的距离S为:
船舶A与船舶B之间的相对速度V1为:
V1=V×COSθ
通过matlab模拟船舶间距离和速度与时间的关系可得出,在船舶A向船舶B运动时其速度在两船连线方向上的分量和两船之间距离的图像(见图2)。
2.基于非欧空间下船舶会遇局面的搭建
搭建出非欧空间下的船舶会遇模型:
1)速度的叠加
运用运动速度的合成原理,将本船B的运动速度叠加至来船A上,将两速度进行矢量的合成得出来船的合速度V。
2)建立黎曼球面
运用Mobius变换和非欧几何学的原理,在本船B上方建立黎曼球面,以O点作为球的最低点和三维直角坐标系的原点,点O′作为球的最高点,球的半径暂定为5倍的船长(其长度大约为船舶邻域的半径)。由于本船相对静止,来船A将在l1上做直线运动。船舶A在l1上运动时运用黎曼球面将船舶A所在的点P与点O′相连接,此连线与球面相交的点为点Q,该点为来船A在球面上的投影点。在来船A不断在l1上运动时运动时,直线O′P会在黎曼球面上产生一段圆弧O′QQ′Q″。
3)利用保形变换求解相对速度VOQ与相对距离S
运用保形变换的思想,通过相关的几何证明易得到,在P点运动时△O′OP与△OQO′的互为相似三角形,通过变换将P点与O点连线的速度VOP进行变换到△OQO′中O点与Q点连线的速度VOQ,同时将平面上的OP之间的距离变换至OQ之间的距离。具体求解VOQ、OQ的过程如下:
现假设直线l1上任意一点P的坐标为(xo,yo,0)、点O′的坐标为(0,0,2r)。
(1)求解点Q的坐标
由P、O′的坐标可设:
球面解析式为:x2+y2+(z-r)2=r2 (Ⅰ)
将Ⅲ式代入Ⅰ式,得:(x0*t)2+(y0*t)2+(r-2r*t)2=r2
(2)求解相对速度与VOQ相对距离S
以O为原点,以平行于直线l1且过点O的直线为x轴,建立如图3所示的o-xyz空间直角坐标系。设船舶A沿直线l1向船舶B运动接近,达到最近距离后并逐渐远离。仅考虑船舶A与点P′之间距离不超过R的情况(根据参考球体M半径r的大小,确定R的大小)。船舶A的运动速度为V,运动的时间为t,过O点作l1的垂线交l1于点P′,设线段OP′的距离为d,则直线l1上任意一点P的坐标(xo,yo,0)为
由(1)得,点Q的坐标(x,y,z)为:
对点Q坐标Q(x(t),y(t),z(t))求一阶导,得到点Q处的速度:
空间下OQ之间的距离为:
运用matlab仿真出相对速度与VOQ相对距离S与时间之间的关系(见图4)。
3.基于黎曼球面的船舶碰撞风险的构建
基于以上对黎曼球面中两船相对速度与相对距离运用matlab仿真结果基本上符合驾驶员基于眼球对来船的速度与距离的认知。当前对于船舶避碰的风险评估的主要方法为相对速度与相对距离的倒数相乘积表针风险,即主要缘由为,碰撞风险是随着两船之间相对速度的增大而增加,同时随着两船之间距离的减小而增大,将碰撞风险的描述的变量划分为V即基于非欧空间的两船之间的相对速度,即基于非欧空间的两船之间距离的倒数。对于非欧空间的船舶避碰风险评估模型的具体步骤为:
碰撞风险:
对于空间情况下的风险描述,运用Matlab进行仿真得到空间下船舶风险、相对速度与时间变化如图3-图5所示。
在非欧空间的情况下,船舶的风险随时间的不断增大,在碰撞点之前的一段时间内是呈现不断增大的趋势,仿真的结果符合现实中驾驶员的的真实感受,当两船不断靠近时,船舶的风险是不断增大的,基本符合两船会遇的风险态势的变化情况,可直观并有效地对船舶会遇时碰撞风险进行预测与防控。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作任何其他形式的限制,而依据本发明的技术实质所作的任何修改或等同变化,仍属于本发明所要求保护的范围。
Claims (3)
1.基于非欧保形变换的船舶避碰模型的建立方法,其特征在于,包括如下步骤:
S1)数据提取:利用船舶的VDR进行VDR数据的提取;
S2)速度叠加:运用速度的合成原理,将本船的速度叠加至来船上,并运用速度的合成原理,将叠加的本船速度与他船速度进行矢量合成,本船认为处于静止状态;
S3)建立黎曼球面:以本船为空间坐标原点,建立空间直角坐标系,以本船作为黎曼球面的最低点,在本船上方建立黎曼球面;
S4)来船进行球面投影:用直线将黎曼球面的最高点与来船相连,直线与 黎曼球面相交于除黎曼球面的最高点的另外一点,来船在黎曼球面进行了球面投影;
S5)相对距离与相对速度投影点:运用黎曼球面将本船平面内的一切物体投影至黎曼球面上,并在黎曼球面上可以找到对应的投影点,在黎曼球面上将来船相对于本船的相对距离与相对速度通过投影点进行反映,在黎曼球面上将来船在平面内的运动转化为圆面的圆弧运动;
S6)求解投影点的坐标:利用几何关系将黎曼球面与来船在平面运动的直线运用方程进行表述,建立黎曼球面最高点与来船位置的直线方程,利用来船的位置点与黎曼球面最高点所连的直线方程与黎曼球面的方程进行联立求解出投影点的坐标;
S7)求取相对距离与相对速度:利用投影点的坐标求取来船与本船之间的相对距离,将投影点坐标进行求导,计算出在投影点处的速度,将该速度投影至投影点与本船连线的直线上,求取出投影点相对于本船的速度,利用投影点相对于本船的距离和相对距离来反映在黎曼球面下本船观测到他船的真实运动情况;
S8)建立船舶碰撞风险模型:根据相对速度与相对距离描述船舶之间的碰撞风险,船舶碰撞的风险是与相对速度成正相关与相对距离成反比,所以建立风险评定方程以相对速度与相对距离的倒数相乘积表征 风险。
2.根据权利要求1所述的基于非欧保形变换的船舶避碰模型的建立方法,其特征在于,步骤S1)中所述VDR数据包括船舶航速、航向、主机参数、AIS数据、雷达图像数据和机舱驾驶台的语音数据。
3.根据权利要求2所述的基于非欧保形变换的船舶避碰模型的建立方法,其特征在于,步骤S4)中在原点、黎曼球面最高点、来船的位置坐标和投影点建立三角形,利用几何三角相似关系可以证明出一对相似三角形关系,利用保形变换思想,来船在二维平面下的直线运动利用黎曼球面投影在球面显示出圆弧曲线运动,可将来船的运动信息通过投影点进行反映。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911155257.3A CN110889221B (zh) | 2019-11-22 | 2019-11-22 | 基于非欧保形变换的船舶避碰模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911155257.3A CN110889221B (zh) | 2019-11-22 | 2019-11-22 | 基于非欧保形变换的船舶避碰模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110889221A CN110889221A (zh) | 2020-03-17 |
CN110889221B true CN110889221B (zh) | 2022-06-17 |
Family
ID=69748428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911155257.3A Active CN110889221B (zh) | 2019-11-22 | 2019-11-22 | 基于非欧保形变换的船舶避碰模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110889221B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111861155B (zh) * | 2020-07-03 | 2022-09-20 | 南京理工大学 | 船舶碰撞风险检测方法、系统、计算机设备和存储介质 |
CN113326640B (zh) * | 2021-08-03 | 2021-11-09 | 北京海兰信数据科技股份有限公司 | 一种确定内河船舶碰撞危险度的方法及系统 |
CN114283621A (zh) * | 2021-12-24 | 2022-04-05 | 武汉理工大学 | 基于船舶交通流的内河渡线会遇风险识别及决策方法 |
CN115657665A (zh) * | 2022-09-30 | 2023-01-31 | 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) | 基于测地线的无人船航路规划方法及存储介质 |
CN115640091B (zh) * | 2022-10-28 | 2023-05-02 | 合肥通用机械研究院有限公司 | 一种基于黎曼球面的交互式多样化图案展示控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009083012A1 (en) * | 2007-12-31 | 2009-07-09 | Mohsen Negm Attia Negm | Sea bottom pilotage stars |
CN105390029A (zh) * | 2015-11-06 | 2016-03-09 | 武汉理工大学 | 基于航迹融合和航迹预测的船舶避碰辅助决策方法及系统 |
CN109887339A (zh) * | 2019-03-20 | 2019-06-14 | 集美大学 | 一种关联避让行为的船舶碰撞危险及危险度评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7860646B2 (en) * | 2007-04-16 | 2010-12-28 | The Boeing Company | Method and apparatus for routing ocean going vessels to avoid treacherous environments |
-
2019
- 2019-11-22 CN CN201911155257.3A patent/CN110889221B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009083012A1 (en) * | 2007-12-31 | 2009-07-09 | Mohsen Negm Attia Negm | Sea bottom pilotage stars |
CN105390029A (zh) * | 2015-11-06 | 2016-03-09 | 武汉理工大学 | 基于航迹融合和航迹预测的船舶避碰辅助决策方法及系统 |
CN109887339A (zh) * | 2019-03-20 | 2019-06-14 | 集美大学 | 一种关联避让行为的船舶碰撞危险及危险度评估方法 |
Non-Patent Citations (3)
Title |
---|
M.B. Mendel & P.H.A.J.M. van Gelder.Visualizing and gauging collision risk.《Risk, Reliability and Safety: Innovating Theory and Practice – Walls, Revie & Bedford (Eds)》.2016, * |
周光耀.考虑坐标系位置瞬变的船舶多自由度耦合模型及参数横摇运动分析.《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》.2018,(第06期), * |
谭忠华.船舶碰撞机理的数值仿真与试验研究.《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》.2012,(第10期),C036-11. * |
Also Published As
Publication number | Publication date |
---|---|
CN110889221A (zh) | 2020-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110889221B (zh) | 基于非欧保形变换的船舶避碰模型的建立方法 | |
Lyu et al. | COLREGS-constrained real-time path planning for autonomous ships using modified artificial potential fields | |
CN110362074B (zh) | 一种基于航迹重规划的水面无人艇动态避碰方法 | |
JP7101255B2 (ja) | 目標対象の運動の向きを予測するための方法、車両制御方法及び装置 | |
CN109753068B (zh) | 一种考虑通信情况的多usv群体协同避碰规划方法 | |
Tam et al. | Cooperative path planning algorithm for marine surface vessels | |
CN109597417B (zh) | 一种基于避碰准则的多usv群体协同避碰规划方法 | |
CN105022270B (zh) | 基于速度矢量坐标系的船自动避碰方法 | |
CN111445086B (zh) | 一种基于pidvca预测复航时间的方法 | |
Xu | Collision avoidance strategy optimization based on danger immune algorithm | |
CN111861155B (zh) | 船舶碰撞风险检测方法、系统、计算机设备和存储介质 | |
CN105807769A (zh) | 无人水下航行器ivfh避碰方法 | |
Lee et al. | A collision avoidance method for multi-ship encounter situations | |
Jincan et al. | Based on ECDIS and AIS ship collision avoidance warning system research | |
KR20210116799A (ko) | 디지털 트윈을 이용한 선박 입출항 관제 시스템 | |
Lu et al. | Preliminary study for motion pose of inshore ships based on point cloud: Estimation of ship berthing angle | |
Ni et al. | A deterministic collision avoidance decision-making system for multi-MASS encounter situation | |
Lin et al. | The fuzzy-based visual intelligent guidance system of an autonomous underwater vehicle: realization of identifying and tracking underwater target objects | |
Zhu et al. | Environment potential field modeling for ship automatic collision avoidance in restricted waters | |
CN110443368B (zh) | 一种冲突状态下船舶集群态势分析方法 | |
CN117008606A (zh) | 一种定线制水域船舶自主航行决策方法及装置 | |
Wu et al. | Deliberative collision avoidance for unmanned surface vehicle based on the directional weight | |
CN107367279B (zh) | 大型船舶高精度导航与靠泊方法 | |
CN114690772A (zh) | 用于海域非循迹航行船舶的会遇时空预测方法及系统 | |
Dai et al. | Testing method of autonomous navigation systems for ships based on virtual-reality integration scenarios |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |