CN113326640B - 一种确定内河船舶碰撞危险度的方法及系统 - Google Patents
一种确定内河船舶碰撞危险度的方法及系统 Download PDFInfo
- Publication number
- CN113326640B CN113326640B CN202110884688.4A CN202110884688A CN113326640B CN 113326640 B CN113326640 B CN 113326640B CN 202110884688 A CN202110884688 A CN 202110884688A CN 113326640 B CN113326640 B CN 113326640B
- Authority
- CN
- China
- Prior art keywords
- ship
- physical model
- target
- equal
- key feature
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000033001 locomotion Effects 0.000 claims abstract description 59
- 150000001875 compounds Chemical class 0.000 claims description 17
- QZXCCPZJCKEPSA-UHFFFAOYSA-N chlorfenac Chemical compound OC(=O)CC1=C(Cl)C=CC(Cl)=C1Cl QZXCCPZJCKEPSA-UHFFFAOYSA-N 0.000 claims description 7
- 230000009545 invasion Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- LFULEKSKNZEWOE-UHFFFAOYSA-N propanil Chemical compound CCC(=O)NC1=CC=C(Cl)C(Cl)=C1 LFULEKSKNZEWOE-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/0206—Control of position or course in two dimensions specially adapted to water vehicles
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明实施例公开了一种确定内河船舶碰撞危险度的方法,包括:基于船舶尺寸构建本船物理模型和目标船物理模型,根据本船和目标船之间相对运动态势和位置关系识别会遇场景,判断本船物理模型和目标船物理模型的多个关键特征点,以本船物理模型和目标船物理模型的多个关键特征点之间的运动趋势变化的核心参数构建相关随机变量和碰撞概率分布模型,分析船舶在不同运动趋势和不同距离下的碰撞发生概率,即碰撞危险度。本发明实施例还公开了一种确定内河船舶碰撞危险度的系统。本发明可以解决不同航行角度对碰撞危险度的影响以及过大的船舶领域范围导致的误报警问题。
Description
技术领域
本发明涉及船舶技术领域,具体而言,涉及一种确定内河船舶碰撞危险度的方法及系统。
背景技术
在智能船舶航行控制领域中,在对内河通航船舶确定碰撞危险度时,一般是先建立船舶领域模型,实时计算本船与目标船领域之间、目标船与本船领域之间的位置关系,再根据目标船对本船领域的入侵程度、本船对目标船领域的入侵程度或本船领域和目标船领域之间的相互入侵程度,来计算本船和目标船之间的碰撞危险度。其中,领域模型多基于船舶尺寸和操纵性能等经验参数建立。这种根据领域入侵计算碰撞危险度的方式,无法反应不同航行角度对碰撞危险度程度的影响,导致计算的结果不够准确。另外,对于领域模型预警范围设定通常过大,在狭窄的内河航道上会产生大量的误报警。
发明内容
本发明的目的在于提供一种确定内河船舶碰撞危险度的方法及系统,可以解决不同航行角度对碰撞危险度的影响,以及过大的船舶领域范围导致的误报警问题。
本发明实施例提供了一种确定内河船舶碰撞危险度的方法,所述方法包括:
基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,所述本船物理模型和所述目标船物理模型分别包括多个特征点;
根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点;
根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型;
根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
作为本发明进一步的改进,所述基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,包括:
以所述本船的船体中心为原点,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立所述本船的船体坐标系;
在所述本船的船体坐标系下,将覆盖所述本船的船体的最小矩形作为所述本船物理模型,该矩形的四个顶点pi作为所述本船物理模型的特征点,其中,i表示所述本船物理模型的特征点的序号,i=1,2,3,4;
在所述本船的船体坐标系下,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点tj作为所述目标船物理模型的特征点,其中,j表示所述目标船物理模型的特征点的序号,j=1,2,3,4。
作为本发明进一步的改进,所述根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点,包括:
根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点。
所述根据所述目标船的位置方向以及所述夹角所在的角度范围,识别所述本船和所述目标船的会遇场景,根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点,包括:
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为对遇且靠近场景,将所述目标船物理模型的特征点t1和t4以及所述本船物理模型的特征点p1和p4作为关键特征点;
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为尾随或被追越且靠近场景,所述目标船物理模型的特征点t2和t3以及所述本船物理模型的特征点p2和p3作为关键特征点。
作为本发明进一步的改进,所述目标船物理模型确定出两个关键特征点,所述本船物理模型确定出两个关键特征点,
所述根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型,包括:
对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
作为本发明进一步的改进,当所述两个连线在所述相对运动方向的垂直方向上的投影D11和D12正负相反时,所述目标船物理模型的边界与所述本船物理模型的边界相交,
将所述目标船物理模型的关键特征点与交点的连线在所述相对运动方向上的投影及所述连线在所述相对运动方向的垂直方向上的投影分别作为所述第一随机变量S和所述第二随机变量D,
其中,所述交点表征所述目标船物理模型的关键特征点沿所述相对运动方向延长线与所述本船物理模型的边界相交的点。
作为本发明进一步的改进,所述根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度,包括:
所述调整后的纵向船舶碰撞概率分布模型为:
当所述目标船沿所述相对运动方向向所述本船行驶时,对所述目标船物理模型的每个关键特征点,基于所述横向船舶碰撞概率分布模型和所述调整后的纵向船舶碰撞概率分布模型计算对应的碰撞概率:
将对所述目标船物理模型的两个关键特征点计算得到的两个碰撞概率中的较大值作为最终的船舶碰撞危险度P。
作为本发明进一步的改进,所述船舶碰撞危险度的取值范围为0~1。
本发明实施例还提供了一种确定内河船舶碰撞危险度的系统,所述系统包括:
物理模型构建模块,用于基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,所述本船物理模型和所述目标船物理模型分别包括多个特征点;
关键特征点识别模块,用于根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点;
概率分布模型确定模块,用于根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型;
碰撞危险度确定模块,用于根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
作为本发明进一步的改进,所述物理模型构建模块包括:以所述本船的船体中心为原点,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立所述本船的船体坐标系;
在所述本船的船体坐标系下,将覆盖所述本船的船体的最小矩形作为所述本船物理模型,该矩形的四个顶点pi作为所述本船物理模型的特征点,其中,i表示所述本船物理模型的特征点的序号,i=1,2,3,4;
在所述本船的船体坐标系下,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点tj作为所述目标船物理模型的特征点,其中,j表示所述目标船物理模型的特征点的序号,j=1,2,3,4。
根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点。
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为对遇且靠近场景,将所述目标船物理模型的特征点t1和t4以及所述本船物理模型的特征点p1和p4作为关键特征点;
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为尾随或被追越且靠近场景,所述目标船物理模型的特征点t2和t3以及所述本船物理模型的特征点p2和p3作为关键特征点。
作为本发明进一步的改进,所述目标船物理模型确定出两个关键特征点,所述本船物理模型确定出两个关键特征点, 所述概率分布模型确定模块包括:对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
作为本发明进一步的改进,当所述两个连线在所述相对运动方向的垂直方向上的投影D11和D12正负相反时,所述目标船物理模型的边界与所述本船物理模型的边界相交,
将所述目标船物理模型的关键特征点与交点的连线在所述相对运动方向上的投影及所述连线在所述相对运动方向的垂直方向上的投影分别作为所述第一随机变量S和所述第二随机变量D,
其中,所述交点表征所述目标船物理模型的关键特征点沿所述相对运动方向延长线与所述本船物理模型的边界相交的点。
作为本发明进一步的改进,所述碰撞危险度确定模块包括:所述调整后的纵向船舶碰撞概率分布模型为:
当所述目标船沿所述相对运动方向向所述本船行驶时,对所述目标船物理模型的每个关键特征点,基于所述横向船舶碰撞概率分布模型和所述调整后的纵向船舶碰撞概率分布模型计算对应的碰撞概率:
将对所述目标船物理模型的两个关键特征点计算得到的两个碰撞概率中的较大值作为最终的船舶碰撞危险度P。
作为本发明进一步的改进,所述船舶碰撞危险度的取值范围为0~1。
本发明的有益效果为:考虑了内河航行环境中船舶尺寸对碰撞危险的影响,以船舶物理模型的多个关键特征点之间的运动趋势变化的核心参数构建相关随机变量和碰撞概率分布模型,可以反应不同航行角度对碰撞危险度的影响。根据两船之间会遇时间变化动态调整碰撞概率分布模型,从概率角度分析船船在不同运动趋势和不同距离下的碰撞发生概率,可以避免误报警,提高报警准确率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一示例性实施例所述的一种确定内河船舶碰撞危险度的方法的流程示意图;
图2为本发明一示例性实施例所述的本船物理模型的示意图;
图3为本发明一示例性实施例所述的本船和目标船相对运动的示意图;
图4为本发明一示例性实施例所述的船舶碰撞概率分布模型的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明,若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后……),则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
另外,在本发明的描述中,所用术语仅用于说明目的,并非旨在限制本发明的范围。术语“包括”和/或“包含”用于指定所述元件、步骤、操作和/或组件的存在,但并不排除存在或添加一个或多个其他元件、步骤、操作和/或组件的情况。术语“第一”、“第二”等可能用于描述各种元件,不代表顺序,且不对这些元件起限定作用。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个及两个以上。这些术语仅用于区分一个元素和另一个元素。结合以下附图,这些和/或其他方面变得显而易见,并且,本领域普通技术人员更容易理解关于本发明所述实施例的说明。附图仅出于说明的目的用来描绘本发明所述实施例。本领域技术人员将很容易地从以下说明中认识到,在不背离本发明所述原理的情况下,可以采用本发明所示结构和方法的替代实施例。
本发明实施例所述的一种确定内河船舶碰撞危险度的方法,如图1所示,所述方法包括:
基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,所述本船物理模型和所述目标船物理模型分别包括多个特征点;
根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点;
根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型;
根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
在智能船舶航行控制技术研究领域,进行船舶-船舶碰撞危险度研究时,当前海上通航船舶碰撞危险度确定方法为:将本船和目标船都看作质点,根据两船运动趋势计算两者之间的最近会遇距离(DCPA)和最近会遇时间(TCPA),根据最小安全会遇距离(SDA)设置危险评判阈值,通过比较DCPA、TCPA和危险阈值的关系判定碰撞危险等级和危险度。然而内河通航环境不同于海上通航环境,其航道狭窄、弯曲、流速流态复杂,评估船船碰撞危险度时不能忽略船舶自身尺寸对危险度计算的影响。内河通航船舶碰撞危险度确定方法为:通过建立船舶领域模型,实时计算本船和目标船与对方船舶领域之间的位置关系,根据目标船对本船领域入侵程度、本船对目标船领域入侵程度、本船领域和目标船领域之间相互入侵程度(三者任选其一)计算碰撞危险度。其中,领域模型多基于船舶尺寸和操纵性能等经验参数建立。然而这种内河船舶碰撞危险度确定方法只根据领域入侵程度计算碰撞危险度,无法反应不同航行角度对碰撞危险程度的影响。
本发明所述方法,考虑了内河航行环境中船舶尺寸对碰撞危险的影响,基于船舶尺寸构建船舶物理模型,根据本船和目标船之间相对运动态势和位置关系识别会遇场景,判断本船物理模型和目标船物理模型的多个关键特征点,以本船物理模型和目标船物理模型的多个关键特征点之间的运动趋势变化的核心参数构建相关随机变量和碰撞概率分布模型,分析船舶在不同运动趋势和不同距离下的碰撞发生概率,即碰撞危险度,可以反应不同航行角度对碰撞危险度的影响,解决不同航行角度对碰撞危险度的影响。
现有技术中,对于领域模型预警范围设定通常过大,在狭窄的内河航道上会产生大量的误报警。本发明所述方法,构建船舶碰撞概率分布模型,并根据两船之间会遇时间变化动态调整碰撞概率分布模型,从概率角度分析船船在不同运动趋势和不同距离下的碰撞发生概率,可以避免误报警,提高报警准确率,解决过大的船舶领域范围导致的误报警问题。
一种可选的实施方式中,所述基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,包括:
以所述本船的船体中心为原点,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立所述本船的船体坐标系;
在所述本船的船体坐标系下,将覆盖所述本船的船体的最小矩形作为所述本船物理模型,该矩形的四个顶点pi作为所述本船物理模型的特征点,其中,i表示所述本船物理模型的特征点的序号,i=1,2,3,4;
在所述本船的船体坐标系下,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点tj作为所述目标船物理模型的特征点,其中,j表示所述目标船物理模型的特征点的序号,j=1,2,3,4。
鉴于内河航道的有限性,船身相对于内河航道并非是小到可以忽略。因此,在内河船舶的碰撞危险度评估计算中,船舶尺寸应作为重要的参数被考虑其中。以所述本船的船体中心为原点O,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立一个直角坐标系,该直角坐标系作为所述本船的船体坐标系,为了简化计算,在所述本船的船体坐标系中,以能覆盖所述本船的船体的最小矩形作为一个所述本船物理模型,矩形的四个顶点作为所述本船物理模型的平面几何特征点,也是用于后续碰撞危险度分析的特征点。例如,如图2所示,在所述本船的船体坐标系中,所述本船的船体坐标系中第二象限内的顶点p1,所述本船的船体坐标系中第一象限内的顶点p2,所述本船的船体坐标系中第三象限内的顶点p3,所述本船的船体坐标系中第四象限内的顶点p4分别作为所述本船的四个特征点,四个特征点的坐标可以根据所述本船的船长、船宽计算得到。
相应的,所述目标船物理模型同样可以用矩形表示。如图3所示,在所述本船的船体坐标系中,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点作为所述目标船物理模型的特征点。例如,如图3所示,所述目标船位于所述本船的船体坐标系中第二象限内,四个顶点t1、t2、t3、t 4均是位于第二象限内的顶点。可以理解的是,所述目标船也可以位于所述本船的船体坐标系中的其他象限内。在所述本船的船体坐标系中,根据所述本船和所述目标船的经纬度、航迹向、航速等导航数据,可推算出所述目标船物理模型的四个特征点的坐标值,图3中,为所述目标船相对于所述本船的运动速度,为所述目标船相对于所述本船的运动方向,为所述本船航速,为所述目标船航速。
一种可选的实施方式中,所述根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点,包括:
根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点。
可以理解的是,内河船舶发生碰撞时,可以视作两个矩形发生相交,可简化为所述目标船物理模型的某条矩形边界与所述本船物理模型的某条矩形边界产生相交。通过分析船船会遇场景,识别会遇时所述本船物理模型和所述目标船物理模型的相交边界。其中,两船在相交时,所述本船物理模型有一条相交边界,作为第一相交边界,所述目标船物理模型有一条相交边界,作为第二相交边界,第一相交边界上的所述本船物理模型的特征点和第二相交边界上所述目标船物理模型的特征点被识别为关键特征点,第一相交边界作为关键边界。因此,将所述目标船物理模型的关键特征点和所述本船物理模型的关键特征点之间的运动趋势变化作为碰撞危险度计算的主要依据。
一种可选的实施方式中,所述目标船的位置方向以及所述夹角的角度范围与所述会遇场景之间具有预设的对照关系,通过所述目标船的位置方向以及所述夹角的角度范围识别得到的所述本船和所述目标船的会遇场景,以及根据识别出的会遇场景得到的所述多个关键特征点如表1的对照表所示。
上述确定得到的本船的两个关键特征点,可以理解为图2中所示出的特征点p1、p2、p3、p4中得到的,相应的,上述确定得到的目标船的两个关键特征点,可以理解为图3中所示出的特征点t1、t2、t3、t 4中得到的。
还可以理解的是,上述场景中,对遇-驶离表示所述本船和所述目标船对遇且两船对遇后相互驶离,对遇-靠近表示所述本船和所述目标船对遇且两船对遇后相互靠近,船首右交叉表示所述目标船从右边方向与所述本船的船首交叉,船尾右交叉表示所述目标船从右边方向与所述本船的船尾交叉,船首左交叉表示所述目标船从左边方向与所述本船的船首交叉,船尾左交叉表示所述目标船从左边方向与所述本船的船尾交叉,船首右追越表示所述目标船从右边方向追越所述本船的船首,船尾右追越表示所述目标船从右边方向追越所述本船的船尾,船首左追越表示所述目标船从左边方向追越所述本船的船首,船尾左追越表示所述目标船从左边方向追越所述本船的船尾,尾随或被追越-驶离表示所述本船被所述目标船尾随或追越且在所述本船被所述目标船尾随或追越后两船相互驶离,尾随或被追越-靠近表示所述本船被所述目标船尾随或追越且在所述本船被所述目标船尾随或追越后两船相互靠近。
一种可选的实施方式中,所述目标船物理模型确定出两个关键特征点,所述本船物理模型确定出两个关键特征点,
所述根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型,包括:
对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
本发明以本船-目标船的相对位置作为输入,设计一个具有两个随机变量的船舶碰撞概率分布模型,如图4所示,两个随机变量分别是:所述目标船物理模型的关键特征点和所述本船物理模型的关键特征点之间的连线沿相对运动方向(所述目标船相对于所述本船的运动方向)和其垂直方向的两个分量。所述目标船物理模型的两个关键特征点中,例如图4所示的关键特征点t1和t2分别作为所述目标船物理模型的第一关键特征点和第二关键特征点,所述本船物理模型的两个关键特征点中,例如图4所示的关键特征点p1和p4分别作为所述本船物理模型的第一关键特征点和第二关键特征点。
对于所述目标船物理模型的第一关键特征点t1,t1其与所述本船物理模型的第一关键特征点p1、第二关键特征点p4之间的连线分别为t1 p1、t1 p4,两个连线t1 p1、t1 p4在相对运动速度方向上的投影分别为S11、S12,两个连线t1 p1、t1 p4在相对运动速度方向的垂直方向上的投影分别为D11、D12。从S11、S12中确定第一随机变量S,从D11、D12中确定第二随机变量D。
可以理解的是,所述投影Smn与碰撞危险度的相关关系为:Min{Smn},表征横向上的船舶碰撞危险程度,所述投影Smn越小说明航行角度导致的碰撞紧急程度越高,即横向上的实时碰撞概率越高,横向上的碰撞危险度越大。所述投影Dmn与碰撞危险度的相关关系为:Min{Dmn},表征纵向上的船舶碰撞危险程度,投影Dmn越小说明航行角导致的碰撞严重程度越高,即纵向上的实时碰撞概率越高,纵向上的碰撞危险度越大。因此,从投影中选取最小投影作为随机变量。
一种可选的实施方式中,当所述两个连线在所述相对运动方向的垂直方向上的投影D11和D12正负相反时,所述目标船物理模型的边界与所述本船物理模型的边界相交,
将所述目标船物理模型的关键特征点与交点的连线在所述相对运动方向上的投影及所述连线在所述相对运动方向的垂直方向上的投影分别作为所述第一随机变量S和所述第二随机变量D,
其中,所述交点表征所述目标船物理模型的关键特征点沿所述相对运动方向延长线与所述本船物理模型的边界相交的点。即所述目标船物理模型的关键特征点沿所述相对运动方向运动一段时间后,与所述本船物理模型的边界发生相交时的点。
一般来说,Smn和Dmn是所述目标船物理模型的关键特征点与所述本船物理模型的关键特征点的连线在相对运动方向和其垂直方向上的投影。但是,当所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点的两个连线在运动方向的垂直方向上的两个投影D11、D12正负相反时,说明该目标船物理模型的边界必然会与所述本船物理模型的边界相交,即所述目标船会与所述本船发生碰撞,此时,以该目标船物理模型的关键特征点与交点的连线在相对运动速度方向及其垂直方向上的投影分别作为第一随机变量S和第二随机变量D。
其中,根据船舶安全航行的经验要求,D的均值例如可以取一倍船长(即所述本船的船长),标准差随着所述目标船到所述本船的距离逐渐缩小也相应变小,不宜取常量。所述目标船在距离所述本船60s所对应位置附近时,标准差可以根据内河船舶航迹预测算法在弯道处的船位残差的标准差来决定,例如可以设定为27米。所述目标船在所述本船附近时,标准差可以取内河船舶GPS误差的标准差,例如可以设定为10米。
式中,TCPA为所述本船和所述目标船的最近会遇时间。
一种可选的实施方式中,所述根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度,包括:
所述调整后的纵向船舶碰撞概率分布模型为:
当所述目标船沿所述相对运动方向向所述本船行驶时,对所述目标船物理模型的每个关键特征点,基于所述横向船舶碰撞概率分布模型和所述调整后的纵向船舶碰撞概率分布模型计算对应的碰撞概率:
将对所述目标船物理模型的两个关键特征点计算得到的两个碰撞概率中的较大值作为最终的船舶碰撞危险度P。
需要指出的是,考虑到所述本船和所述目标船经过最近会遇点之后,两船之间在纵向上碰撞危险度较会遇之前下降很多,但由于两船横向距离较近,仍需保持警惕。因此,对纵向船舶碰撞概率分布模型进行调整。分别针对所述目标船物理模型的两个关键特征点,计算碰撞概率P1和P2,取两者之中的较大值作为最终碰撞概率P,碰撞危险度CRI即等于该碰撞概率P,取值范围为0~1。
本发明实施例所述的一种确定内河船舶碰撞危险度的系统,所述系统包括:
物理模型构建模块,用于基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,所述本船物理模型和所述目标船物理模型分别包括多个特征点;
关键特征点识别模块,用于根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中确定多个关键特征点;
概率分布模型确定模块,用于根据所述多个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型;
碰撞危险度确定模块,用于根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
一种可选的实施方式中,所述物理模型构建模块包括:
以所述本船的船体中心为原点,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立所述本船的船体坐标系;
在所述本船的船体坐标系下,将覆盖所述本船的船体的最小矩形作为所述本船物理模型,该矩形的四个顶点pi作为所述本船物理模型的特征点,其中,i表示所述本船物理模型的特征点的序号,i=1,2,3,4;
在所述本船的船体坐标系下,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点tj作为所述目标船物理模型的特征点,其中,j表示所述目标船物理模型的特征点的序号,j=1,2,3,4。
一种可选的实施方式中,所述关键特征点识别模块,包括:
根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点。
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为对遇且靠近场景,将所述目标船物理模型的特征点t1和t4以及所述本船物理模型的特征点p1和p4作为关键特征点;
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为尾随或被追越且靠近场景,所述目标船物理模型的特征点t2和t3以及所述本船物理模型的特征点p2和p3作为关键特征点。
一种可选的实施方式中,所述目标船物理模型确定出两个关键特征点,所述本船物理模型确定出两个关键特征点, 所述概率分布模型确定模块包括:对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
一种可选的实施方式中,当所述两个连线在所述相对运动方向的垂直方向上的投影D11和D12正负相反时,所述目标船物理模型的边界与所述本船物理模型的边界相交,
将所述目标船物理模型的关键特征点与交点的连线在所述相对运动方向上的投影及所述连线在所述相对运动方向的垂直方向上的投影分别作为所述第一随机变量S和所述第二随机变量D,
其中,所述交点表征所述目标船物理模型的关键特征点沿所述相对运动方向延长线与所述本船物理模型的边界相交的点。
一种可选的实施方式中,所述碰撞危险度确定模块包括:所述调整后的纵向船舶碰撞概率分布模型为:
当所述目标船沿所述相对运动方向向所述本船行驶时,对所述目标船物理模型的每个关键特征点,基于所述横向船舶碰撞概率分布模型和所述调整后的纵向船舶碰撞概率分布模型计算对应的碰撞概率:
将对所述目标船物理模型的两个关键特征点计算得到的两个碰撞概率中的较大值作为最终的船舶碰撞危险度P。
一种可选的实施方式中,所述船舶碰撞危险度的取值范围为0~1。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
此外,本领域普通技术人员能够理解,尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本领域技术人员应理解,尽管已经参考示例性实施例描述了本发明,但是在不脱离本发明的范围的情况下,可进行各种改变并可用等同物替换其元件。另外,在不脱离本发明的实质范围的情况下,可进行许多修改以使特定情况或材料适应本发明的教导。因此,本发明不限于所公开的特定实施例,而是本发明将包括落入所附权利要求范围内的所有实施例。
Claims (9)
1.一种确定内河船舶碰撞危险度的方法,其特征在于,所述方法包括:
基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,分别将覆盖所述本船的船体和所述目标船的船体的最小矩形作为所述本船物理模型和所述目标船物理模型,所述本船物理模型和所述目标船物理模型分别包括四个特征点;
根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中分别确定出两个关键特征点;
根据四个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型,包括:
对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
2.如权利要求1所述的方法,其中,所述基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,包括:
以所述本船的船体中心为原点,沿所述本船的船体纵向指向船艏为y轴,沿所述本船的船体横向指向右为x轴,建立所述本船的船体坐标系;
在所述本船的船体坐标系下,将覆盖所述本船的船体的最小矩形作为所述本船物理模型,该矩形的四个顶点pi作为所述本船物理模型的特征点,其中,i表示所述本船物理模型的特征点的序号,i=1,2,3,4;
在所述本船的船体坐标系下,将覆盖所述目标船的船体的最小矩形作为所述目标船物理模型,该矩形的四个顶点tj作为所述目标船物理模型的特征点,其中,j表示所述目标船物理模型的特征点的序号,j=1,2,3,4。
所述根据所述目标船的位置方向以及所述夹角所在的角度范围,识别所述本船和所述目标船的会遇场景,根据识别出的会遇场景,确定所述本船和所述目标船会遇时,所述本船物理模型和所述目标船物理模型的相交边界,将所述相交边界上的所述本船物理模型的特征点和所述目标船物理模型的特征点确定为所述多个关键特征点,包括:
当所述夹角大于或等于-180°且小于-90°或所述夹角大于90°且小于或等于180°时,所述会遇场景识别为对遇且靠近场景,将所述目标船物理模型的特征点t1和t4以及所述本船物理模型的特征点p1和p4作为关键特征点;
5.如权利要求1所述的方法,其中,当所述两个连线在所述相对运动方向的垂直方向上的投影D11和D12正负相反时,所述目标船物理模型的边界与所述本船物理模型的边界相交,
将所述目标船物理模型的关键特征点与交点的连线在所述相对运动方向上的投影及所述连线在所述相对运动方向的垂直方向上的投影分别作为所述第一随机变量S和所述第二随机变量D,
其中,所述交点表征所述目标船物理模型的关键特征点沿所述相对运动方向延长线与所述本船物理模型的边界相交的点。
7.如权利要求1所述的方法,其中,所述根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度,包括:
所述调整后的纵向船舶碰撞概率分布模型为:
当所述目标船沿所述相对运动方向向所述本船行驶时,对所述目标船物理模型的每个关键特征点,基于所述横向船舶碰撞概率分布模型和所述调整后的纵向船舶碰撞概率分布模型计算对应的碰撞概率:
将对所述目标船物理模型的两个关键特征点计算得到的两个碰撞概率中的较大值作为最终的船舶碰撞危险度P。
8.如权利要求1所述的方法,其中,所述船舶碰撞危险度的取值范围为0~1。
9.一种确定内河船舶碰撞危险度的系统,其特征在于,所述系统包括:
物理模型构建模块,用于基于本船和目标船的尺寸,分别构建本船物理模型和目标船物理模型,其中,分别将覆盖所述本船的船体和所述目标船的船体的最小矩形作为所述本船物理模型和所述目标船物理模型,所述本船物理模型和所述目标船物理模型分别包括四个特征点;
关键特征点识别模块,用于根据所述本船和所述目标船之间的相对运动态势和位置关系识别会遇场景,根据识别出的会遇场景,从所述本船物理模型和所述目标船物理模型的多个特征点中分别确定出两个关键特征点;
概率分布模型确定模块,用于根据四个关键特征点之间的相对位置确定随机变量,并基于所述随机变量构建横向船舶碰撞概率分布模型和纵向船舶碰撞概率分布模型,包括:
对所述目标船物理模型的每个关键特征点,分别确定所述目标船物理模型的关键特征点与所述本船物理模型的两个关键特征点之间的两个连线;
确定所述两个连线在相对运动速度方向上的投影Smn以及在所述相对运动方向的垂直方向上的投影Dmn,其中,m表示所述目标船物理模型的关键特征点的序号,m=1,2,n表示所述本船物理模型的关键特征点的序号,n=1,2;
从所述投影Smn中选取最小投影作为第一随机变量S,从所述投影Dmn中选取最小投影作为第二随机变量D;
基于所述第一随机变量S,构建横向船舶碰撞概率分布模型为:
基于所述第二随机变量D,构建纵向船舶碰撞概率分布模型为:
碰撞危险度确定模块,用于根据所述本船和所述目标船之间会遇时间的变化,动态调整所述纵向船舶碰撞概率分布模型,基于所述横向船舶碰撞概率分布模型和调整后的纵向船舶碰撞概率分布模型,确定船舶碰撞危险度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110884688.4A CN113326640B (zh) | 2021-08-03 | 2021-08-03 | 一种确定内河船舶碰撞危险度的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110884688.4A CN113326640B (zh) | 2021-08-03 | 2021-08-03 | 一种确定内河船舶碰撞危险度的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113326640A CN113326640A (zh) | 2021-08-31 |
CN113326640B true CN113326640B (zh) | 2021-11-09 |
Family
ID=77426890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110884688.4A Active CN113326640B (zh) | 2021-08-03 | 2021-08-03 | 一种确定内河船舶碰撞危险度的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113326640B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116564136B (zh) * | 2023-05-20 | 2024-03-19 | 广西综合交通大数据研究院 | 基于ais的船舶碰撞预测方法 |
CN117115704B (zh) * | 2023-08-03 | 2024-04-02 | 武汉理工大学 | 一种基于多传感器融合的海洋搜救系统及方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3729407B1 (en) * | 2017-12-22 | 2022-07-20 | Rolls-Royce plc | A collision avoidance method and system for marine vessels |
CN108922247B (zh) * | 2018-07-25 | 2020-11-06 | 重庆大学 | 一种基于ais的船舶-航标碰撞危险度估算方法 |
CN109887339B (zh) * | 2019-03-20 | 2020-09-15 | 集美大学 | 一种关联避让行为的船舶碰撞危险及危险度评估方法 |
CN110660272A (zh) * | 2019-09-27 | 2020-01-07 | 南京莱斯网信技术研究院有限公司 | 一种适用于内河水上船舶碰撞预警检测的方法 |
CN110889221B (zh) * | 2019-11-22 | 2022-06-17 | 武汉理工大学 | 基于非欧保形变换的船舶避碰模型的建立方法 |
CN111721321B (zh) * | 2020-05-19 | 2022-02-22 | 交通运输部东海航海保障中心上海航标处 | 一种航标智能化避碰方法及系统 |
CN111561932B (zh) * | 2020-05-27 | 2023-07-18 | 武汉理工大学 | 一种基于虚拟力的船舶导航方法 |
CN112580261B (zh) * | 2020-12-23 | 2022-10-14 | 武汉理工大学 | 一种多船避碰决策方法及装置 |
CN113112871B (zh) * | 2021-04-14 | 2022-06-24 | 上海海事大学 | 一种考虑船舶尺度的船撞桥危险度计算方法 |
-
2021
- 2021-08-03 CN CN202110884688.4A patent/CN113326640B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113326640A (zh) | 2021-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113326640B (zh) | 一种确定内河船舶碰撞危险度的方法及系统 | |
Im et al. | Potential risk ship domain as a danger criterion for real-time ship collision risk evaluation | |
CN109298708B (zh) | 一种融合雷达与光电信息的无人艇自主避障方法 | |
CN108445879A (zh) | 一种基于预测碰撞危险区域的无人艇避障方法 | |
CN113744570B (zh) | 一种桥区水域船舶防撞预警方法及装置 | |
CN105807769B (zh) | 无人水下航行器ivfh避碰方法 | |
CN108803313A (zh) | 一种基于海流预测模型的路径规划方法 | |
Blaich et al. | Fast grid based collision avoidance for vessels using A∗ search algorithm | |
CN113032896B (zh) | 一种基于船舶驾驶员偏好的避碰辅助决策方法 | |
CN114387824B (zh) | 一种符合国际海上避碰规则的避碰转向判定方法 | |
CN112285696B (zh) | 一种雷达目标跟踪方法及系统 | |
CN112561232A (zh) | 繁忙水域船舶交通冲突风险度量方法 | |
CN111984012A (zh) | 一种无人船自主避障及航线规划的方法 | |
Yan et al. | Obstacle avoidance for unmanned undersea vehicle in unknown unstructured environment | |
CN114527744A (zh) | 基于天牛须搜索优化的无人帆船路径跟踪制导方法 | |
CN114061565A (zh) | 一种无人船舶slam及其应用方法 | |
CN113851019A (zh) | 一种基于北斗定位的船舶碰撞预警方法及装置 | |
CN114384917B (zh) | 一种基于场论的船舶实时避碰方法 | |
CN117369441A (zh) | 一种考虑船舶运动学和CORLEGs的自适应智能船舶路径规划方法 | |
CN111984006B (zh) | 融合海流及尺度差异影响的无人艇多目标会遇避碰方法 | |
CN115273555B (zh) | 一种航道交汇区的船舶避碰决策方法 | |
CN115657682A (zh) | 一种关于无人艇在狭窄水道中安全航行的避碰决策方法 | |
Tengesdal et al. | Joint stochastic prediction of vessel kinematics and destination based on a maritime traffic graph | |
CN115407780A (zh) | 一种基于三阶段避障策略的无人船局部路径规划方法 | |
Lee et al. | Data association for autonomous ships based on virtual simulation environment |
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 |