CN104809350B - 海洋垃圾和生物个体的区分方法 - Google Patents
海洋垃圾和生物个体的区分方法 Download PDFInfo
- Publication number
- CN104809350B CN104809350B CN201510213209.0A CN201510213209A CN104809350B CN 104809350 B CN104809350 B CN 104809350B CN 201510213209 A CN201510213209 A CN 201510213209A CN 104809350 B CN104809350 B CN 104809350B
- Authority
- CN
- China
- Prior art keywords
- sample
- bion
- detecting object
- ocean rubbish
- study
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明提供一种海洋垃圾和生物个体的区分方法包括:建立学习样例;基于学习样例进行马尔科夫模型学习;以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断。其中,以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断具体包括:对探测物体进行初步定位;在一段时间内跟踪探测物体的位置,获得各监测时刻探测物体的移动参数,并形成移动参数变化序列;计算获得的移动参数变化序列与学习所获得马尔科夫模型的吻合程度;当计算所得吻合程度大于或等于设定阈值时,判断探测物体为海洋垃圾。
Description
技术领域
本发明涉及海洋生态环境检测技术领域,且特别涉及一种海洋垃圾和生物个体的区分方法。
背景技术
目前评估漂浮垃圾的数量、种类和分布的三种技术:来自期刊《海洋环境科学》第16卷第2期的文章《海洋环境中漂浮垃圾的监测方法》介绍了以下三种技术:
1.统计人员通过调查过往船只来估算漂浮垃圾的种类及数量。该方法根据统计人员日均抛投的垃圾数量,进行相关的数理统计后,得出某海域年均漂浮垃圾的数量。由于该方法需要收集足够的信息和专业的统计人员,由于统计人员的个体存在差异,因此统计出来的结论有较大的差异。
2.通过现场调查来确定漂浮垃圾的密度和类型。(1)使用飞机对某一海域漂浮垃圾进行调查观测,收集有关漂浮垃圾分布和数量的资料;(2)借助其它海上监测船、游船或货船共同进行,根据船只的航线进行有限观测。两种方式均要求有合适的海况(<3级)和良好的能见度,前者置信水平高而后者费用低。
3.人工海滩观测监测漂浮垃圾的状况。由于海洋中的部分漂浮垃圾受风向和海流的影响,不一定完全漂浮到岸边,该方法使总估计量偏低。
综上,漂浮垃圾的检测还缺乏经济、有效的方法;而海面下悬浮垃圾和海底垃圾的检测亟需进行相关的研究和实践。
发明内容
本发明为了克服现有的海洋垃圾检测技术存在检测精度和经济费用难平衡的问题,提供一种经济、有效的海洋垃圾和生物个体的区分方法。
为了实现上述目的,本发明提供一种海洋垃圾和生物个体的区分方法包括:
建立学习样例;
基于学习样例进行马尔科夫模型学习;
以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断,具体包括:
对探测物体进行初步定位;
在一段时间内跟踪探测物体的位置,获得各监测时刻探测物体的移动参数,并形成移动参数变化序列;
计算获得的移动参数变化序列与学习所获得马尔科夫模型的吻合程度;
当计算所得吻合程度大于或等于设定阈值时,判断探测物体为海洋垃圾。
于本发明一实施例中,建立学习样例的步骤包括:
对监测海域范围内的多个样本物体进行初步定位;
在一段时间内跟踪探测各个样本物体的位置,获得各监测时刻物体的移动参数,并形成移动参数变化序列样例;
对监测海域范围内的多个样本物体进行人工标注,识别的多个样本物体是海洋垃圾还是生物个体。
于本发明一实施例中,移动参数为移动速度和移动方向变化。
于本发明一实施例中,移动参数的获得步骤包括:
在一段时间内跟踪被初步定为后的探测物体或样本物体的位置,测量探测物体或样本物体单位时间内的移动速度;
计算各监测时刻的探测物体或样本物体的移动方向变化。
于本发明一实施例中,基于学习样例进行马尔科夫模型学习,包括计算各监测时刻探测物体移动速度的转换概率和移动方向变化的转换概率。
于本发明一实施例中,对探测物体的初步定位或样本物体的初步定位为根据探测物体或样本物体信号波,通过立体区域定位方法进行定位。
于本发明一实施例中,立体区域定位方法为基于“椭球腔”模型的立体区域定位方法。
于本发明一实施例中,基于“椭球腔”模型的立体区域定位方法包括:
超声波传感器任意一节点发射超声波信号,其它节点接收直接发射的超声波以及探测物体或样本物体反射回来的超声波;
分别包含有探测物体或样本物体距发射点最近点和最远点且绕发射点与任意接收点连线旋转形成的两个椭球面,交集形成“椭球腔”;
改变接收点取得更多的“椭球腔”,“椭球腔”之间进行交集得出质心;
所得的质心形成质心群,求出质心群的质心,所述质心为探测物体或样本物体的位置。
综上所述,本发明提供的海洋垃圾和生物个体的区分方法与现有技术相比,具有以下优点:
由于海洋中的生物具有自主移动性,因此它们的状态变化是没有规律可循的。然而,海洋中垃圾是不具自主移动性的,因此它们的状态变化是具有一定规律可循的。本发明通过在特定海域的水文特征下,通过不断测量采样,记录海洋中各个样本物体在不同时刻下的移动参数,形成移动参数变化序列样例,多个样本的移动参数变化序列样例形成了学习样例,进而基于学习样例进行马尔科夫模型学习。在后续检测中,通过测量各个时刻探测物体的移动参数,并形成移动参数变化序列,通过计算形成的移动参数变化序列于马尔科夫模型的吻合程度来判断探测物体为海洋垃圾或是生物个体。
本发明提供的海洋垃圾和生物个体的区分方法是一种基于马尔科夫模型的区分方法,具有很高的区分精度。同时,区分方法可借助于计算机来完成,不仅具有很高的区分效率,同时具有极低的区分成本,有效解决现有海洋垃圾检测方法检测精度和经费难平衡的问题。
为让本发明的上述和其它目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合附图,作详细说明如下。
附图说明
图1所示为本发明一实施例提供的海洋垃圾和生物个体区分方法的流程图。
图2所示为图1中建立学习样例的流程图。
图3所示为图1中以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断的流程图。
具体实施方式
如图1所示,本实施例提供的海洋垃圾和生物个体的区分方法包括:
步骤S10、建立学习样例,包括:
步骤S11、对监测海域范围内的多个样本物体进行初步定位。
于本实施例中,采用基于“椭球腔”模型的立体区域定位方法对样本物体进行初步定位。然而,本发明对此不作任何限定。于其它实施例中,可采用无线传感器网络节点定位技术等其它定位方式进行定位。
“椭球腔”模型的定位原理为:
以3×3阵列为例:以海岸水平线为Z轴正方向,竖直向
下为Y轴正方向垂直于发射阵列的法向量方向为X轴正方向。以I点发射
为例,首先取I(x1,y1,z1,),J(x2,y2,z2),K(x3,y3,z3),B(x4,y4,z4)为参考点,O(x,y,z)为待
测物体距发射点I(x1,y1,z1,)最近点,P(x′,y′,z′)为发射点I(x1,y1,z1,)最远可达点,J点
接收到I点发射的信号时刻和O点反射回来的信号时刻分别为t1和t1′,K点接收到I点发射的
信号时刻和O点反射回来的信号时刻分别为t2和t2′,B点接收到I点发射的信号时刻和O点反
射回来的信号时刻分别为t3和t3′;则:
dioj=v·(t1′-t1)+dij,
diok=v·(t2′-t2)+dik (1-1)
diob=v·(t3′-t3)+dib
则以I,J为焦点,O为椭球上任意一点可得椭球TQIOJ方程为:
同理椭球TQIOK和椭球TQIOB方程分别为:
联立方程(1-2)、(1-3)、(1-4)可得O(x,y,z)坐标。
同理,对于点I(x1,y1,z1,),J(x2,y2,z2),K(x3,y3,z3),B(x4,y4,z4)为参考点,最远可达点P(x′,y′,z′)为椭球上一点,J点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为t4和t4′,K点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为t5和t5′,B点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为t6和t6′,则可得反射路径长度分别为:
dipj=v·(t4′-t4)+dij
dipk=v·(t5′-t5)+dik (1-5)
dipb=v·(t6′-t6)+dib
基于三点的椭球TQIPJ、TQIPK和TQIPB方程分别为:
联立方程(1-6)、(1-7)、(1-8)可得P(x′,y′,z′)坐标。
任意焦点相同,以O,P表面一点所形成的两个椭球表面作为内外椭球面,即形成一个“椭球腔”,如上述方程即可得到三个“椭球腔”分别为:(TQIOJ,TQIPJ)、 (TQIOK,TQIPK)和(TQIOB,TQIPB),则可得三个“椭球腔”的交集,并求出质心,记为Q(X1,Y1,Z1)。
Q(X1,Y1,Z1)即为样本物体的位置。
步骤S12、在一段时间内跟踪探测各个样本物体的位置,获得各监测时刻物体的移动参数,并形成移动参数变化序列样例。于本实施例中,移动参数为样本物体的移动速度和移动方向变化。得到该移动参数的步骤为:
设某一样本物体在不同时刻的三维坐标为Q(x,y,z)。如T0时刻该样本物体的坐标是Q(x0,y0,z0),Ti时刻的坐标是Q(xi,yi,zi)。在一段时间内跟踪被初步定为后的样本物体的位置,测量样本物体单位时间内的移动速度。具体而言,在Ti-1到Ti时间物体的移动距离为:
Ti-1到Ti时间物体O的平均移动速度:
Vi-1,i=Di-1,i/(Ti-Ti-1) (2)
已知样本物体在不同时间段的平均速度v0,v1,v2,……。已知T0,T1和T2时刻样本物体的坐标,则可以用向量来表示样本物体的移动方向:
两个向量的夹角为:
则单位时间内样本物体的移动方向的变化情况可以用以下公式表达:
其中t0=(T2-T0)/2 (5)
步骤S13、对监测海域范围内的多个样本物体进行人工标注,识别的多个样本物体是海洋垃圾还是生物个体。每个样本物体形成移动参数变化序列样例<(V1,V2,…,Vn),(M1,M2,…,Mn)>,建立学习样例。其中Vi为第i时刻样本物体的移动速度,Mi为第i时刻样本物体的移动方向变化。
步骤S20、基于学习样例进行马尔科夫模型学习。具体而言:
首先,设Vmax为最大的速度,Vmin为最小的速度,Vavg为平均速度。计状态α1为高速区α2为中速区α3为低速区根据每个
海洋垃圾的记录,计算各状态转换的概率矩阵其中aij为状态αi转移为
αj的概率,如a12为状态α1转移到状态α2的概率。本实施例划分三个速度区域。然而,本发明对
此不作限定。
接着,设Mmax为最大的方向变化,Mmin为最小的方向变化,Mavg为平均方向变化。计状
态β1为快速变化区β2为中速变化区 β3为
低速变化区根据每个海洋垃圾的记录,计算各状态转换的概率矩阵其中bij为状态βi转移为βj的概率。本实施例划分三个方向变化区域。然
而,本发明对此不作限定。
最后,计算样本物体速度序列吻合程度阈值PS1和方向变化序列吻合程度阈值PS2。
PS1的计算方式如下:
设Vi∈αk,Vj∈αl,则PV(Vi,Vj)=akl。 (6)
速度序列吻合程度阈值
其中,k,l=1,2,3;PV(Vi,Vj)为Vi所属状态αk转移到Vj所属状态αl的概率。
PS2的计算方式如下:
设Mi∈βk,Mj∈βl,则PM(Mi,Mj)=bkl。 (8)
方向变化序列吻合程度阈值
其中,k,l=1,2,3;PM(Mi,Mj)为Mi所属状态βk转移到Mj所属状态βl的概率。
步骤S30、以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断,具体包括:
步骤S31、对探测物体进行初步定位。该步骤与步骤S11相同,采用基于“椭球腔”模型的立体区域定位方法对探测物体进行定位,获得QT(X1,Y1,Z1)。
步骤S32、在一段时间内跟踪探测物体的位置,获得各监测时刻探测物体的移动参数,并形成移动参数变化序列。该步骤与步骤S12相同,经公式(1)-(2)-(3)-(4)-(5)计算后形成的移动参数变化序列为<(VT1,VT2,…,VTn),(MT1,MT2,…,MTn)〉。其中VTi为探测物体T在第i时刻的移动速度,MTi为探测物体T在第i时刻移动方向变化。
步骤S33、计算获得的移动参数变化序列<(VT1,VT2,…,VTn),(MT1,MT2,…,MTn)〉与学习所获得马尔科夫模型的吻合程度,包括速度序列吻合程度P速度吻合和P方向变化吻合。其中P速度吻合是将序列(VT1,VT2,…,VTn)代入公式(6)和(7)后获得;P 方向变化吻合是将序列(MT1,MT2,…,MTn)代入公式(8)和(9)后获得。
步骤S34、当计算所得吻合程度大于或等于设定阈值时,判断探测物体为海洋垃圾。具体而言,当P速度吻合>=PS1且P方向变化吻合>=PS2,则该物体被判断为是海洋垃圾。
综上所述,本发明提供的海洋垃圾和生物个体的区分方法与现有技术相比,具有以下优点:
由于海洋中的生物具有自主移动性,因此它们的状态变化是没有规律可循的。然而,海洋中垃圾是不具自主移动性的,因此它们的状态变化是具有一定 规律可循的。本发明通过在特定海域的水文特征下,通过不断测量采样,记录海洋中各个样本物体在不同时刻下的移动参数,形成移动参数变化序列样例,多个样本的移动参数变化序列样例形成了学习样例,进而基于学习样例进行马尔科夫模型学习。在后续检测中,通过测量各个时刻探测物体的移动参数,并形成移动参数变化序列,通过计算形成的移动参数变化序列于马尔科夫模型的吻合程度来判断探测物体为海洋垃圾或是生物个体。
本发明提供的海洋垃圾和生物个体的区分方法是一种基于马尔科夫模型的区分方法,具有很高的区分精度。同时,区分方法可借助于计算机来完成,不仅具有很高的区分效率,同时具有极低的区分成本,有效解决现有海洋垃圾检测方法检测精度和经费难平衡的问题。
虽然本发明已由较佳实施例揭露如上,然而并非用以限定本发明,任何熟知此技艺者,在不脱离本发明的精神和范围内,可作些许的更动与润饰,因此本发明的保护范围当视权利要求书所要求保护的范围为准。
Claims (7)
1.一种海洋垃圾和生物个体的区分方法,其特征在于,包括:
建立学习样例;
基于学习样例进行马尔科夫模型学习;
以学习获得的马尔科夫模型对探测物体进行海洋垃圾和生物个体的分类判断,具体包括:
对探测物体进行初步定位;
在一段时间内跟踪探测物体的位置,获得各监测时刻探测物体的移动参数,并形成移动参数变化序列;
计算获得的移动参数变化序列与学习所获得马尔科夫模型的吻合程度;
当计算所得吻合程度大于或等于设定阈值时,判断探测物体为海洋垃圾;
所述移动参数为移动速度和移动方向变化。
2.根据权利要求1所述的海洋垃圾和生物个体的区分方法,其特征在于,所述建立学习样例的步骤包括:
对监测海域范围内的多个样本物体进行初步定位;
在一段时间内跟踪探测各个样本物体的位置,获得各监测时刻物体的移动参数,并形成移动参数变化序列样例;
对监测海域范围内的多个样本物体进行人工标注,识别所述的多个样本物体是海洋垃圾还是生物个体。
3.根据权利要求1或2所述的海洋垃圾和生物个体的区分方法,其特征在于,所述移动参数的获得步骤包括:
在一段时间内跟踪被初步定为后的探测物体或样本物体的位置,测量探测物体或样本物体单位时间内的移动速度;
计算各监测时刻的探测物体或样本物体的移动方向变化。
4.根据权利要求1或2所述的海洋垃圾和生物个体的区分方法,其特征在于,所述基于学习样例进行马尔科夫模型学习,包括计算各监测时刻探测物体移动速度的转换概率和移动方向变化的转换概率。
5.根据权利要求1或2所述的海洋垃圾和生物个体的区分方法,其特征在于,对探测物体的初步定位或样本物体的初步定位为根据探测物体或样本物体信号波,通过立体区域定位方法进行定位。
6.根据权利要求5所述的海洋垃圾和生物个体的区分方法,其特征在于,所述立体区域定位方法为基于“椭球腔”模型的立体区域定位方法。
7.根据权利要求6所述的海洋垃圾和生物个体的区分方法,其特征在于,所述基于“椭球腔”模型的立体区域定位方法包括:
超声波传感器任意一节点发射超声波信号,其它节点接收直接发射的超声波以及探测物体或样本物体反射回来的超声波;
分别包含有探测物体或样本物体距发射点最近点和最远点且绕发射点与任意接收点连线旋转形成的两个椭球面,交集形成“椭球腔”;
改变接收点取得更多的“椭球腔”,“椭球腔”之间进行交集得出质心;
所得的质心形成质心群,求出质心群的质心,所述质心为探测物体或样本物体的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510213209.0A CN104809350B (zh) | 2015-04-29 | 2015-04-29 | 海洋垃圾和生物个体的区分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510213209.0A CN104809350B (zh) | 2015-04-29 | 2015-04-29 | 海洋垃圾和生物个体的区分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104809350A CN104809350A (zh) | 2015-07-29 |
CN104809350B true CN104809350B (zh) | 2018-11-09 |
Family
ID=53694166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510213209.0A Expired - Fee Related CN104809350B (zh) | 2015-04-29 | 2015-04-29 | 海洋垃圾和生物个体的区分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104809350B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115560736B (zh) * | 2022-09-22 | 2024-04-05 | 山东省地质测绘院 | 一种海洋工程测绘用的水下测绘装置及方法 |
CN115372975B (zh) * | 2022-10-25 | 2023-01-06 | 青州鑫聚隆装备制造有限公司 | 一种基于河道保洁船的超声定位*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508287A (zh) * | 2011-10-17 | 2012-06-20 | 大连海事大学 | 一种水下物体探测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502915B (zh) * | 2014-12-17 | 2017-04-26 | 河海大学常州校区 | 基于主动式目标探测原理的超短基线水下定位方法及*** |
-
2015
- 2015-04-29 CN CN201510213209.0A patent/CN104809350B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508287A (zh) * | 2011-10-17 | 2012-06-20 | 大连海事大学 | 一种水下物体探测装置 |
Non-Patent Citations (4)
Title |
---|
A multi-resolution hidden markov model for optimal detection, tracking, separation, and classification of marine mammal vocalizations;Brian F. Harrison e.t.;《OCEANS 2008》;20080918;第2008年卷;全文 * |
基于时空联合的水下运动目标检测方法;彭复员 等;《华中科技大学学报》;20071231;第35卷(第12期);摘要 * |
基于马尔可夫过程的水下运动目标启发式搜索;吴芳 等;《电子与信息学报》;20100531;第32卷(第5期);全文 * |
基于高斯-马尔可夫模型的海洋环境辨识方法;杜金燕 等;《计算机仿真》;20091231;第26卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104809350A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901806B (zh) | 一种智能船舶靠泊辅助***及方法 | |
CN109212540A (zh) | 基于激光雷达***的测距方法、装置及可读存储介质 | |
CN108008099A (zh) | 一种污染源定位方法 | |
CN103969639B (zh) | 五波束探鱼仪的信号处理***及其信号处理方法 | |
CN107942336B (zh) | 适用于复杂水环境的鱼类超声波标记精密定位方法及*** | |
CN103885030A (zh) | 一种无线传感器网络移动节点的定位方法 | |
CN109901205A (zh) | 一种水下机器人多传感器融合与运动轨迹预测方法 | |
Sutarto | Application of large scale particle image velocimetry (LSPIV) to identify flow pattern in a channel | |
Vortmeyer-Kley et al. | Detecting and tracking eddies in oceanic flow fields: a Lagrangian descriptor based on the modulus of vorticity | |
CN104457703A (zh) | 基于多视觉小通道气液两相流参数测量***及方法 | |
CN104809350B (zh) | 海洋垃圾和生物个体的区分方法 | |
CN105387842A (zh) | 基于感知驱动的自航式海底地形地貌测绘***及测绘方法 | |
CN110207688A (zh) | 一种基于特征矢量的磁信标快速定向方法与*** | |
KR102082263B1 (ko) | 수중위치 추정 시스템 및 방법 | |
CN110058245A (zh) | 基于云模型的低频主动拖线阵声呐浅海探测效能评估方法 | |
CN110187302A (zh) | 一种基于单信标的水下蛙人自导航方法 | |
Pidcock et al. | The spatial variability of vertical velocity in an Iceland basin eddy dipole | |
CN116595870A (zh) | 一种基于动平台的尾迹识别与流场探测方法、*** | |
Viikmäe et al. | Quantification and characterization of mesoscale eddies with different automatic identification algorithms | |
CN109212598A (zh) | 基于直达波反演的三维空间二次定位方法 | |
O'Connor et al. | A multi-modal event detection system for river and coastal marine monitoring applications | |
CN109490868A (zh) | 一种基于分布式垂直线列阵的海上目标运动分析方法 | |
KR102561270B1 (ko) | 구조 구난을 위한 해상 표류물의 위치 추정 방법 및 장치 | |
Wang et al. | Reconstruction of 3-D surface waves generated by moving submerged sphere based on stereo imaging principle | |
Wang et al. | Multibeam Line Finding Problem Based on Geometric Analysis and Genetic Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20191104 Address after: Room A1244, Building 2889 Baoan Road, Jiading District, Shanghai, 201801 Patentee after: Shanghai Guangchen Information Technology Co.,Ltd. Address before: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park is 18 street. Patentee before: Zhejiang Gongshang University |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181109 |
|
CF01 | Termination of patent right due to non-payment of annual fee |