CN107330803A - 一种基于多功能雷达船只违法捕捞的行为判别方法 - Google Patents
一种基于多功能雷达船只违法捕捞的行为判别方法 Download PDFInfo
- Publication number
- CN107330803A CN107330803A CN201710483906.7A CN201710483906A CN107330803A CN 107330803 A CN107330803 A CN 107330803A CN 201710483906 A CN201710483906 A CN 201710483906A CN 107330803 A CN107330803 A CN 107330803A
- Authority
- CN
- China
- Prior art keywords
- mrow
- ship
- illegal
- alarm
- time
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 241000251468 Actinopterygii Species 0.000 claims abstract description 14
- 230000033001 locomotion Effects 0.000 claims abstract description 10
- 230000006399 behavior Effects 0.000 claims description 51
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000003542 behavioural effect Effects 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims description 3
- 230000001960 triggered effect Effects 0.000 claims description 2
- 239000003643 water by type Substances 0.000 abstract description 6
- 241000209094 Oryza Species 0.000 description 3
- 235000007164 Oryza sativa Nutrition 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 235000009566 rice Nutrition 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000009372 pisciculture Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/886—Radar or analogous systems specially adapted for specific applications for alarm systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Primary Health Care (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Electromagnetism (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Traffic Control Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于多功能雷达船只违法捕捞的行为判别方法，首先预设违法捕捞的告警规则和预设监管区域集合S；其次，通过实时处理雷达扫描回波数据，获取预设监管区域范围内船只的位置、大小及运动状态等信息；再次，根据预设违法捕捞告警规则，判别目标船只是否违法，并输出告警级别；最后，根据告警级别，触发不同类型告警。本发明基于实际执法经验，结合雷达探测技术，通过违法捕捞判别规则实现违法船只的自动识别，有效提升基于多功能雷达***的水域监管能力，及时、有效的掌握水上船只动态，为执法人员和监管部门提供及时、有效的船只违法信息，提高执法监管效率，降低执法成本，有效打击水上违法捕捞行为。
Description
技术领域
本发明涉及海洋渔业违法捕捞监管技术领域，特别是一种基于多功能雷达船只违法捕捞的行为判别方法。
背景技术
中国是一个海洋大国，拥有300多万平方公里海洋领土，蕴藏着丰富的渔业资源，为人类提供了大量优质食物，同时带来了巨大的经济效益；然而近年来，随着违法捕捞日益猖獗，违法捕捞手段日益现代化，造成渔业资源结构失衡，破坏了水域生态环境，挑战我国渔业资源的可持续发展，严重影响了我国生态文明建设的进程。
目前传统渔政执法，依托群众举报，结合执法船巡逻，实现水域的监控和管理，这样的执法方式，因无法全面实时掌握水域态势，缺乏有效的信息互通，造成执法效率低，执法成本高和执法被动；伴随着现代传感器手段的发展，特别是雷达、红外光电和AIS等设备，可实现水域全天候、全时空监控，实现执法由被动向主动的转变变为可能，但仅依托雷达、红外光电和AIS等设备，只能获取水面船只实时情况，无法有效判断船只是否违法，需要对这些设备采集的数据进行进一步的分析和处理。
本发明正是基于实际执法经验，结合雷达探测技术，通过对违法捕捞的行为实现违法船只的自动识别，有效提升基于多功能雷达***的水域监管能力，实现全天候、全时空水面违法船只的监控，解决了违法捕捞监管困难的问题，为执法人员、监管部门提供及时、有效的船只违法信息，提高执法监管效率，降低执法成本，有效打击水上违法捕捞行为。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足，而提供一种基于多功能雷达船只违法捕捞的行为判别方法，该基于多功能雷达船只违法捕捞的行为判别方法，基于实际执法经验，并结合多功能雷达探测技术，自动判别目标是否存在违法捕捞的行为，继而发出相关告警信息，为执法人员、管理部门提供决策依据，从而解决渔政监管困难的问题。
为解决上述技术问题，本发明采用的技术方案是：
一种基于多功能雷达船只违法捕捞的行为判别方法，包括如下步骤。
步骤1，预设违法捕捞的告警规则：违法捕捞的告警规则包括四种类型，分别为禁渔告警，禁入告警，超时告警，超速告警。
每种告警规则的预设方法如下：
——禁渔告警：当目标船只在监管区域中的滞留时间超过预设阈值时，再结合目标船只的行为，进行智能判断是否违法。
目标船只的行为包括单船捕鱼行为和双船捕鱼行为。
是否违法的智能判断方法，包括如下两个步骤：
步骤11，单船捕鱼禁渔告警：采用滑窗算法一，也即先设定一个滑动时间间隔一和一个时间窗口一，时间每向前滑动一个时间间隔一，统计一次时间窗口一内目标船只的运动轨迹，并计算曲线曲率，当曲线曲率超过预设阈值时，则判断目标船只的行为是单船捕鱼行为，且目标船只违法。
步骤12，双船捕鱼禁渔告警：采用滑窗算法二，也即先设定一个滑动时间间隔二和一个时间窗口二，时间每向前滑动一个时间间隔二，在时间窗口二内：首先，提取与目标船只间距小于预设距离阈值的船只集合Tdis；然后，计算目标船只与船只集合Tdis内所有船只的航迹相关系数、航速相关系数，若航迹相关系数与航速相关系数均超过预设阈值，则判定为两艘目标船只的行为是双船捕鱼行为；若双船捕鱼行为持续时间超过预设时间阈值时，则定义为违法船只，并设定为一级告警。
——禁入告警：在监管区域中只要有目标船只进入，即定义为违法船只，并设定为一级告警。
——超时告警：目标船只在监管区域中的停留时间超过预设的时间阈值时，则定义为违法船只。
——限速告警，当目标船只的当前航速低于监管区域设定的最低限速阈值或者高于最高限速阈值时，判断为违法船只，并设定为三级告警。
步骤2，预设监管区域集合S：监管区域为多个，每个监管区域均为多边形区域，每个监管区域均包括不在同一条直线上的至少三个顶点，每个顶点坐标均用经纬度表示；每个监管区域均进行步骤1中告警规则类型的预设，所有监管区域及监管区域预设的告警规则类型，构成预设的监管区域集合S。
步骤3，获取监管区域船只信息：自动雷达标绘仪ARPA自动接收雷达回波数据并计算获取所有目标船只的信息，形成目标船只集合T；目标集合T中每个目标船只的信息均包含位置信息和运动状态信息，其中，位置信息以经纬度表示，运动状态信息包括航速和航向。
步骤4，判别目标船只是否违法：目标船只是否违法的具体判别方法如下：
步骤41，获取步骤1预设的监管区域集合S。
步骤42，选取监管区域集合S中的某一个监管区域S(i)，同时获取步骤2中的目标船只集合T；遍历目标船只集合T，当目标船只信息符合监管区域S(i)预设的告警规则时，则提供报警。
步骤43，重复步骤42，遍历监管区域集合S，完成对所有违法船只的判断。
步骤5，触发告警：输入目标船只告警信息，触发告警装置。
步骤11，单船捕鱼禁渔告警时，曲线曲率的计算公式如下：
式中：angRate为曲线曲率；n为时间窗口一内航迹的点数；i为第i个航迹点；angi为第i个航迹点对应的航向，单位为弧度；angi+1为第i+1个航迹点对应的航向，单位为弧度。
步骤11，单船捕鱼禁渔告警时，曲线曲率的预设阈值包括第一级预设阈值、第二级预设阈值和第三级预设阈值，且第三级预设阈值＞第二级预设阈值＞第一级预设阈值；曲线曲率超过第一级预设阈值时，设定为三级告警；曲线曲率超过第二级预设阈值时，设定为二级告警；曲线曲率超过第三级预设阈值时，设定为一级告警。
第一级预设阈值取值范围为20％～30％；第二级预设阈值取值范围为30％～40％；第三级预设阈值取值范围为40％～50％。
步骤1中禁渔告警时，滞留时间的预设阈值取值范围为10～120分钟。
步骤11，双船捕鱼禁渔告警时，采用Pearson系数法计算集合Tdis内所有船只的航迹相关系数、航速相关系数；
航迹相关系数计算方法为：航迹的位置以经纬度表示，假设两艘目标船只分别为A船和B船，则A船与B船航迹的Pearson系数，计算公式如下：
式中：r1为A船与B船航迹的Pearson系数；n为时间窗口二内航迹的点数；i为第i个航迹点；m取值为2，当j＝1时，Amn为时间窗口二内A船航迹的第n个航迹点的经度，Bmn为时间窗口二内B船航迹的第n个航迹点的经度；当j＝2时，Amn为时间窗口二内A船航迹的第n个航迹点的纬度，Bmn为时间窗口二内B船航迹的第n个航迹点的纬度；A为时间窗口二内A船航迹所有经度与纬度的均值；B为时间窗口二内B船航迹所有经度与纬度的均值。
航速相关系数计算公式如下：
式中：r为速度的Pearson系数；x、y分别为两艘目标船只的速度，单位为m/s；n为时间窗口二内航迹的点数；i为第i个航迹点。
步骤1中，禁渔告警且目标船只的行为是双船捕鱼行为时，两艘船只的距离阈值取值范围为50～1500米，双船捕鱼行为持续时间的预设阈值取值范围为8～30分钟；两艘船只航迹相关系数的预设阈值取值范围为0.5～0.9，两艘船只航速相关系数的预设阈值取值范围为0.5～0.9。
步骤1中，超时告警时，停留时间的预设时间阈值包括第一级时间阈值、第二级时间阈值和第三级时间阈值；且第一级时间阈值＜第二级时间阈值＜第三级时间阈值；当目标船只在监管区停留时间超过第一级时间阈值时，则定义为违法船只，并设定为三级告警；目标船只在监管区停留时间超过第二级时间阈值时，升级为二级告警；目标船只在监管区停留时间超过第三级时间阈值时，升级为一级告警。
步骤1中，超时告警时，第一级时间阈值取值范围为0～60分钟，第二级时间阈值取值范围为60～120分钟，第三级时间阈值取值范围为120～240分钟。
步骤1中，限速告警时，监管区域设定的最低限速阈值取值范围为0.1～5节；最高限速阈值取值范围为5～15节。
本发明采用上述方法后，能同时判断多个监管区域中，是否存在违法船只，并产生告警信息。首先预设违法捕捞的告警规则和预设监管区域集合S；其次，通过实时处理雷达扫描回波数据，获取预设监管区域范围内船只的位置、大小及运动状态等信息；再次，根据预设违法捕捞告警规则，判别目标船只是否违法，并输出告警级别；最后，根据告警级别，触发不同类型告警，从而有效提高执法监管效率，降低执法成本，打击水上违法捕捞行为。
附图说明
图1是本发明一种基于多功能雷达船只违法捕捞的行为判别方法的流程示意图。
具体实施方式
下面结合附图和具体较佳实施方式对本发明作进一步详细的说明。
如图1所示，一种基于多功能雷达船只违法捕捞的行为判别方法，包括如下步骤。
步骤1，预设违法捕捞的告警规则，也即输入违法捕捞告警规则。
违法捕捞的告警规则包括四种类型，分别为禁渔告警，禁入告警，超时告警，超速告警。
第一种：禁渔告警。
先将目标船只在监管区域中的滞留时间与预设阈值进行比较，本发明中滞留时间的预设阈值的优选取值范围为10～120分钟，进一步优选为10分钟。
目标船只的行为包括单船捕鱼行为和双船捕鱼行为。
当目标船只在监管区域中的滞留时间超过预设阈值，结合目标船只的行为，进行智能判断是否违法。
是否违法的智能判断方法，包括如下两个步骤。
步骤11，单船捕鱼禁渔告警
采用滑窗算法一，也即先设定一个滑动时间间隔一和一个时间窗口一，时间每向前滑动一个时间间隔一，统计一次时间窗口一内目标船只的运动轨迹，并计算曲线曲率，当曲线曲率超过预设阈值时，则判断目标船只的行为是单船捕鱼行为，且目标船只违法。
曲线曲率的预设阈值包括第一级预设阈值、第二级预设阈值和第三级预设阈值，且第三级预设阈值＞第二级预设阈值＞第一级预设阈值。其中，第一级预设阈值优选取值范围为20％～30％，进一步优选为25％；第二级预设阈值优选取值范围为30％～40％，进一步优选为35％；第三级预设阈值优选取值范围为40％～50％，进一步优选为45％。
曲线曲率超过第一级预设阈值时，设定为三级告警；曲线曲率超过第二级预设阈值时，设定为二级告警；曲线曲率超过第三级预设阈值时，设定为一级告警。
根据两船经纬度计算两船距离公式如下：
式中：LAT、LON分别为纬度和经度，单位为弧度；a为地球长半径(WGS-84坐标系下，取值为6378137.0米)；b为地球短半径，单位为米；e为地球偏心率；X和Y分别为目标在空间直角坐标系中的横纵坐标，单位为米；d为两点(X0，Y0)、(X1，Y1)之间的距离，单位为米。
其中，曲线曲率的计算公式如下：
式中：angRate为曲线曲率；n为时间窗口一内航迹的点数；i为第i个航迹点；angi为第i个航迹点对应的航向，单位为弧度；angi+1为第i+1个航迹点对应的航向，单位为弧度。
步骤12，双船捕鱼禁渔告警
采用滑窗算法二，也即先设定一个滑动时间间隔二和一个时间窗口二，其中，时间间隔二可以与时间间隔一相等，也可以不等，如时间间隔一为2分钟，而时间间隔二为5分钟等。
时间每向前滑动一个时间间隔二，在时间窗口二内，采用相关性法，进行一次所有目标船只距离相关性的判断。
相关性法的具体判断过程为：假设当前船只为目标船只，以目标船只为基准，对一个固定的间隔距离范围内(如100米的范围内)。
首先，提取与目标船只间距小于预设距离阈值的船只集合Tdis；两艘船只的上述预设距离阈值优选取值范围为50～1500米，进一步优选为100米。
然后，计算目标船只与船只集合Tdis内所有船只的航迹相关系数、航速相关系数，本发明中，采用Pearson系数法计算集合Tdis内所有船只的航迹相关系数、航速相关系数。作为替换，也可采用现有技术中已知的其他方法计算航迹相关系数、航速相关系数，均在本发明的保护范围之内。
若航迹相关系数与航速相关系数均超过预设阈值，认为两艘目标船只的航迹和速度基本一致，则判定为两艘目标船只的行为是双船捕鱼行为。
若双船捕鱼行为持续时间超过预设时间阈值时，则定义为违法船只，并设定为一级告警。
双船捕鱼行为持续时间的预设时间阈值优选取值范围为8～30分钟，进一步优选为15分钟。
两艘船只航迹相关系数的阈值优选取值范围为0.5～0.9，进一步优选为0.65。
两艘船只航速相关系数的阈值优选取值范围为0.5～0.9，进一步优选为0.7。
航迹相关系数计算方法方法如下。
航迹的位置以经纬度表示，假设两艘目标船只分别为A船和B船，则A船与B船航迹的Pearson系数，计算公式如下：
式中：r1为A船与B船航迹的Pearson系数；n为时间窗口二内航迹的点数(也即：在时间窗口二内，A船与B船同时被跟踪到的点迹集合)；i为第i个航迹点；m取值为2，当j＝1时，Amn为时间窗口二内A船航迹的第n个航迹点的经度，Bmn为时间窗口二内B船航迹的第n个航迹点的经度；当j＝2时，Amn为时间窗口二内A船航迹的第n个航迹点的纬度，Bmn为时间窗口二内B船航迹的第n个航迹点的纬度；为时间窗口二内A船航迹所有经度与纬度的均值；为时间窗口二内B船航迹所有经度与纬度的均值。
航速相关系数的计算公式如下：
式中：r为速度的Pearson系数；x、y分别为两艘目标船只的速度，单位为m/s；n为时间窗口二内航迹的点数；i为第i个航迹点。
第二种：禁入告警。
在监管区域中只要有目标船只进入，即定义为违法船只，并设定为一级告警。
判断一个点是否位于多边形区域内的判断方法为：顺序遍历多边形所有的线段，取线段的两个端点分别为(x1，y1)(x2，y2)，ARPA目标的位置为(x0，y0)，线段所在直线与y＝y0的交点为：
若x1<po int x≤x2且y1<y0≤y2，则该顶点与线段相交。统计交点的个数，若交点个数为偶数，则点在多边形外；若交点个数为奇数，则点在多边形内。
第三种：超时告警。
目标船只在监管区域中的停留时间超过预设的时间阈值时，则定义为违法船只。
超时告警时，停留时间的预设时间阈值包括第一级时间阈值、第二级时间阈值和第三级时间阈值；且第一级时间阈值＜第二级时间阈值＜第三级时间阈值。本发明中，第一级时间阈值优选取值范围为0～60分钟，进一步优选为30分钟；第二级时间阈值优选取值范围为60～120分钟，进一步优选为90分钟；第三级时间阈值优选取值范围为120～240分钟，进一步优选为180分钟。
当目标船只在监管区停留时间超过第一级时间阈值时，则定义为违法船只，并设定为三级告警；目标船只在监管区停留时间超过第二级时间阈值时，升级为二级告警；目标船只在监管区停留时间超过第三级时间阈值时，升级为一级告警。
第四种：限速告警。
当目标船只的当前航速低于监管区域设定的最低限速阈值或者高于最高限速阈值时，判断为违法船只，并设定为三级告警。
限速告警时，监管区域设定的最低限速阈值优选取值范围为0.1～5节，进一步优选为3节；最高限速阈值优选取值范围为5～15节，进一步优选为7节。
步骤2，预设监管区域集合S。
步骤21，预设监管区域范围，也即输入监管区域范围。
监管区域为多个，多个监管区域相互之间可以重叠，互不影响。
每个监管区域均为多边形区域，每个监管区域均包括不在同一条直线上的至少三个顶点，每个顶点坐标均用经纬度表示。
步骤22，预设告警规则类型：每个监管区域均进行步骤1中告警规则类型的预设。
步骤23，形成预设监管区域集合S：步骤21预设的所有监管区域及步骤22中监管区域预设的告警规则类型，构成预设的监管区域集合S。
步骤3，获取监管区域船只信息。
自动雷达标绘仪ARPA自动接收雷达回波数据并计算获取所有目标船只的信息，形成目标船只集合T；目标集合T中每个目标船只的信息均包含位置信息和运动状态信息，其中，位置信息以经纬度表示，运动状态信息包括航速和航向。
步骤4，判别目标船只是否违法：目标船只是否违法的具体判别方法如下。
步骤41，获取步骤1预设的监管区域集合S。
步骤42，选取监管区域集合S中的某一个监管区域S(i)，同时获取步骤2中的目标船只集合T；遍历目标船只集合T，当目标船只信息符合监管区域S(i)预设的告警规则时，则提供报警。
步骤43，重复步骤42，遍历监管区域集合S，完成对所有违法船只的判断。
步骤5，触发告警：输入步骤4产生的目标船只告警信息，包括目标船只的位置、航速、航向和告警类型，触发告警装置，并显示船只位置、航速和航向等信息，为监管人员提供告警信息。
综上所述，本发明基于实际执法经验，结合雷达探测技术，通过违法捕捞判别规则实现违法船只的自动识别，有效提升基于多功能雷达***的水域监管能力，及时、有效的掌握水上船只动态，为执法人员和监管部门提供及时、有效的船只违法信息，提高执法监管效率，降低执法成本，有效打击水上违法捕捞行为。
以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。
Claims (10)
1.一种基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：包括如下步骤：
步骤1，预设违法捕捞的告警规则：违法捕捞的告警规则包括四种类型，分别为禁渔告警，禁入告警，超时告警，超速告警；
每种告警规则的预设方法如下：
——禁渔告警：当目标船只在监管区域中的滞留时间超过预设阈值时，再结合目标船只的行为，进行智能判断是否违法；目标船只的行为包括单船捕鱼行为和双船捕鱼行为；
是否违法的智能判断方法，包括如下两个步骤：
步骤11，单船捕鱼禁渔告警：采用滑窗算法一，也即先设定一个滑动时间间隔一和一个时间窗口一，时间每向前滑动一个时间间隔一，统计一次时间窗口一内目标船只的运动轨迹，并计算曲线曲率，当曲线曲率超过预设阈值时，则判断目标船只的行为是单船捕鱼行为，且目标船只违法；
步骤12，双船捕鱼禁渔告警：采用滑窗算法二，也即先设定一个滑动时间间隔二和一个时间窗口二，时间每向前滑动一个时间间隔二，在时间窗口二内：首先，提取与目标船只间距小于预设距离阈值的船只集合Tdis；然后，计算目标船只与船只集合Tdis内所有船只的航迹相关系数、航速相关系数，若航迹相关系数与航速相关系数均超过预设阈值，则判定为两艘目标船只的行为是双船捕鱼行为；若双船捕鱼行为持续时间超过预设时间阈值时，则定义为违法船只，并设定为一级告警；
——禁入告警：在监管区域中只要有目标船只进入，即定义为违法船只，并设定为一级告警；
——超时告警：目标船只在监管区域中的停留时间超过预设的时间阈值时，则定义为违法船只；
——限速告警，当目标船只的当前航速低于监管区域设定的最低限速阈值或者高于最高限速阈值时，判断为违法船只，并设定为三级告警；
步骤2，预设监管区域集合S：监管区域为多个，每个监管区域均为多边形区域，每个监管区域均包括不在同一条直线上的至少三个顶点，每个顶点坐标均用经纬度表示；每个监管区域均进行步骤1中告警规则类型的预设，所有监管区域及监管区域预设的告警规则类型，构成预设的监管区域集合S；
步骤3，获取监管区域船只信息：自动雷达标绘仪ARPA自动接收雷达回波数据并计算获取所有目标船只的信息，形成目标船只集合T；目标集合T中每个目标船只的信息均包含位置信息和运动状态信息，其中，位置信息以经纬度表示，运动状态信息包括航速和航向；
步骤4，判别目标船只是否违法：目标船只是否违法的具体判别方法如下：
步骤41，获取步骤1预设的监管区域集合S；
步骤42，选取监管区域集合S中的某一个监管区域S(i)，同时获取步骤2中的目标船只集合T；遍历目标船只集合T，当目标船只信息符合监管区域S(i)预设的告警规则时，则提供报警；
步骤43，重复步骤42，遍历监管区域集合S，完成对所有违法船只的判断；
步骤5，触发告警：输入目标船只告警信息，触发告警装置。
2.根据权利要求1所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤11，单船捕鱼禁渔告警时，曲线曲率的计算公式如下：
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>R</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>ang</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>ang</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中：angRate为曲线曲率；n为时间窗口一内航迹的点数；i为第i个航迹点；angi为第i个航迹点对应的航向，单位为弧度；angi+1为第i+1个航迹点对应的航向，单位为弧度。
3.根据权利要求1或2所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤11，单船捕鱼禁渔告警时，曲线曲率的预设阈值包括第一级预设阈值、第二级预设阈值和第三级预设阈值，且第三级预设阈值＞第二级预设阈值＞第一级预设阈值；曲线曲率超过第一级预设阈值时，设定为三级告警；曲线曲率超过第二级预设阈值时，设定为二级告警；曲线曲率超过第三级预设阈值时，设定为一级告警。
4.根据权利要求3所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：第一级预设阈值取值范围为20％～30％；第二级预设阈值取值范围为30％～40％；第三级预设阈值取值范围为40％～50％。
5.根据权利要求1所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中禁渔告警时，滞留时间的预设阈值取值范围为10～120分钟。
6.根据权利要求1所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中，禁渔告警且目标船只的行为是双船捕鱼行为时，采用Pearson系数法计算集合Tdis内所有船只的航迹相关系数、航速相关系数；
航迹相关系数计算方法为：航迹的位置以经纬度表示，假设两艘目标船只分别为A船和B船，则A船与B船航迹的Pearson系数，计算公式如下：
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>&times;</mo>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</msqrt>
</mfrac>
</mrow>
式中：r1为A船与B船航迹的Pearson系数；n为时间窗口二内航迹的点数；i为第i个航迹点；m取值为2，当j＝1时，Amn为时间窗口二内A船航迹的第n个航迹点的经度，Bmn为时间窗口二内B船航迹的第n个航迹点的经度；当j＝2时，Amn为时间窗口二内A船航迹的第n个航迹点的纬度，Bmn为时间窗口二内B船航迹的第n个航迹点的纬度；为时间窗口二内A船航迹所有经度与纬度的均值；为时间窗口二内B船航迹所有经度与纬度的均值；
航速相关系数计算公式如下：
<mrow>
<mi>r</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>n</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msqrt>
<mrow>
<mi>n</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<mi>n</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>y</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
式中：r为速度的Pearson系数；x、y分别为两艘目标船只的速度，单位为m/s；n为时间窗口二内航迹的点数；i为第i个航迹点。
7.根据权利要求6所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中，禁渔告警且目标船只的行为是双船捕鱼行为时，两艘船只的距离阈值取值范围为50～1500米，双船捕鱼行为持续时间的预设阈值取值范围为8～30分钟；两艘船只航迹相关系数的预设阈值取值范围为0.5～0.9，两艘船只航速相关系数的预设阈值取值范围为0.5～0.9。
8.根据权利要求1所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中，超时告警时，停留时间的预设时间阈值包括第一级时间阈值、第二级时间阈值和第三级时间阈值；且第一级时间阈值＜第二级时间阈值＜第三级时间阈值；当目标船只在监管区停留时间超过第一级时间阈值时，则定义为违法船只，并设定为三级告警；目标船只在监管区停留时间超过第二级时间阈值时，升级为二级告警；目标船只在监管区停留时间超过第三级时间阈值时，升级为一级告警。
9.根据权利要求8所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中，超时告警时，第一级时间阈值取值范围为0～60分钟，第二级时间阈值取值范围为60～120分钟，第三级时间阈值取值范围为120～240分钟。
10.根据权利要求1所述的基于多功能雷达船只违法捕捞的行为判别方法，其特征在于：步骤1中，限速告警时，监管区域设定的最低限速阈值取值范围为0.1～5节；最高限速阈值取值范围为5～15节。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710483906.7A CN107330803B (zh) | 2017-06-23 | 2017-06-23 | 一种基于多功能雷达船只违法捕捞的行为判别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710483906.7A CN107330803B (zh) | 2017-06-23 | 2017-06-23 | 一种基于多功能雷达船只违法捕捞的行为判别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107330803A true CN107330803A (zh) | 2017-11-07 |
| CN107330803B CN107330803B (zh) | 2020-05-05 |
Family
ID=60194074
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710483906.7A Active CN107330803B (zh) | 2017-06-23 | 2017-06-23 | 一种基于多功能雷达船只违法捕捞的行为判别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107330803B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108776937A (zh) * | 2018-04-16 | 2018-11-09 | 上海海事大学 | 一种基于双目云台摄像机的船只违法捕捞行为判别方法及*** |
| CN110288855A (zh) * | 2019-06-05 | 2019-09-27 | 上海鹰觉科技有限公司 | 针对海上船只行为的分类方法及其分类*** |
| CN110390288A (zh) * | 2019-04-26 | 2019-10-29 | 上海鹰觉科技有限公司 | 基于计算机视觉的目标智能搜索、定位和取证***及方法 |
| CN110736971A (zh) * | 2019-11-05 | 2020-01-31 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种岸基雷达海杂波测量区域非合作目标实时识别方法 |
| CN111028545A (zh) * | 2019-12-18 | 2020-04-17 | 广州大学华软软件学院 | 一种湖泊船只管理与监控*** |
| CN112435449A (zh) * | 2020-11-24 | 2021-03-02 | 中国水产科学研究院南海水产研究所 | 一种全数字智能化渔业违规船舶识别排查***及方法 |
| CN114217307A (zh) * | 2021-11-29 | 2022-03-22 | 福建新继船舶服务有限公司 | 一种基于雷达和ais数据的船舶异常行为判定的方法 |
| CN115423225A (zh) * | 2022-11-07 | 2022-12-02 | 中国水产科学研究院南海水产研究所 | 一种基于大数据的渔港作业管理方法及*** |
| CN116743970A (zh) * | 2023-08-14 | 2023-09-12 | 安徽塔联智能科技有限责任公司 | 一种具有视频ai预警分析的智能管理平台 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040046689A1 (en) * | 2000-08-12 | 2004-03-11 | Rees Huw D | Signal processing |
| CN102854503A (zh) * | 2011-06-27 | 2013-01-02 | 镇江光宁航海电子科技有限公司 | 一种应用于边防的基于数字化航海雷达技术海域监控*** |
| CN103812942A (zh) * | 2014-02-21 | 2014-05-21 | 北京国交信通科技发展公司 | 渔政监控*** |
| CN105894860A (zh) * | 2016-05-05 | 2016-08-24 | 国家***第三海洋研究所 | 一种海上设施警戒*** |
| CN105905259A (zh) * | 2016-04-18 | 2016-08-31 | 太仓弘杉环保科技有限公司 | 一种多功能雷达桅杆信号灯及其工作方法 |
-
2017
- 2017-06-23 CN CN201710483906.7A patent/CN107330803B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040046689A1 (en) * | 2000-08-12 | 2004-03-11 | Rees Huw D | Signal processing |
| CN102854503A (zh) * | 2011-06-27 | 2013-01-02 | 镇江光宁航海电子科技有限公司 | 一种应用于边防的基于数字化航海雷达技术海域监控*** |
| CN103812942A (zh) * | 2014-02-21 | 2014-05-21 | 北京国交信通科技发展公司 | 渔政监控*** |
| CN105905259A (zh) * | 2016-04-18 | 2016-08-31 | 太仓弘杉环保科技有限公司 | 一种多功能雷达桅杆信号灯及其工作方法 |
| CN105894860A (zh) * | 2016-05-05 | 2016-08-24 | 国家***第三海洋研究所 | 一种海上设施警戒*** |
Non-Patent Citations (1)
| Title |
|---|
| 刘娟: ""空管***中雷达告警技术的研究与应用"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108776937A (zh) * | 2018-04-16 | 2018-11-09 | 上海海事大学 | 一种基于双目云台摄像机的船只违法捕捞行为判别方法及*** |
| CN108776937B (zh) * | 2018-04-16 | 2022-03-01 | 上海海事大学 | 一种基于双目云台摄像机的船只违法捕捞行为判别方法及*** |
| CN110390288A (zh) * | 2019-04-26 | 2019-10-29 | 上海鹰觉科技有限公司 | 基于计算机视觉的目标智能搜索、定位和取证***及方法 |
| CN110288855A (zh) * | 2019-06-05 | 2019-09-27 | 上海鹰觉科技有限公司 | 针对海上船只行为的分类方法及其分类*** |
| CN110736971A (zh) * | 2019-11-05 | 2020-01-31 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种岸基雷达海杂波测量区域非合作目标实时识别方法 |
| CN111028545A (zh) * | 2019-12-18 | 2020-04-17 | 广州大学华软软件学院 | 一种湖泊船只管理与监控*** |
| CN112435449A (zh) * | 2020-11-24 | 2021-03-02 | 中国水产科学研究院南海水产研究所 | 一种全数字智能化渔业违规船舶识别排查***及方法 |
| CN114217307A (zh) * | 2021-11-29 | 2022-03-22 | 福建新继船舶服务有限公司 | 一种基于雷达和ais数据的船舶异常行为判定的方法 |
| CN115423225A (zh) * | 2022-11-07 | 2022-12-02 | 中国水产科学研究院南海水产研究所 | 一种基于大数据的渔港作业管理方法及*** |
| CN116743970A (zh) * | 2023-08-14 | 2023-09-12 | 安徽塔联智能科技有限责任公司 | 一种具有视频ai预警分析的智能管理平台 |
| CN116743970B (zh) * | 2023-08-14 | 2023-11-21 | 安徽塔联智能科技有限责任公司 | 一种具有视频ai预警分析的智能管理平台 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107330803B (zh) | 2020-05-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107330803A (zh) | 一种基于多功能雷达船只违法捕捞的行为判别方法 | |
| CN112182133B (zh) | 一种基于ais数据的船舶徘徊检测方法 | |
| CN107622231B (zh) | 一种水天一体智能漂浮物收集***及其收集方法 | |
| Welcker et al. | Displacement of seabirds by an offshore wind farm in the North Sea | |
| Zhang et al. | Enhance the AIS data availability by screening and interpolation | |
| CN115603466B (zh) | 一种基于人工智能视觉识别的船舶岸电*** | |
| CN204667645U (zh) | 智能船桥碰撞监测预警*** | |
| CN113362654B (zh) | 一种海缆防锚损的船舶监控预警*** | |
| CN104537891B (zh) | 一种船舶轨迹实时预测方法 | |
| Wang et al. | Use of AIS data for performance evaluation of ship traffic with speed control | |
| CN108776937B (zh) | 一种基于双目云台摄像机的船只违法捕捞行为判别方法及*** | |
| CN107909601A (zh) | 一种适用于航标的船舶防撞预警视频检测***及检测方法 | |
| CN110208803A (zh) | 桥梁处船舶防撞雷达监测识别方法 | |
| CN108663686A (zh) | 一种基于激光雷达的游泳池溺水监控装置及方法 | |
| CN104485023B (zh) | 船舶冲突解脱的规划方法 | |
| CN110163519A (zh) | 面向基地攻防任务的uuv红蓝方威胁评估方法 | |
| CN107609099A (zh) | 一种船舶违章取证与信息查询*** | |
| Masroeri et al. | IUU fishing and transhipment identification with the miss of AIS data using Neural Networks | |
| CN203025963U (zh) | 监测内河航道船舶交通量的*** | |
| CN113542697B (zh) | 基于深度学习的采砂船自适应监控***的自适应监控方法 | |
| CN104849705B (zh) | 局部均匀杂波协方差矩阵自适应估计方法 | |
| Masroeri et al. | Analysis of fuzzy logic systems types 1 and 2 in identifying of IUU fishing and transshipment: A case study in Indonesia’s vulnerable waters | |
| Li et al. | [Retracted] The Key Technologies of Marine Multiobjective Ship Monitoring and Tracking Based on Computer Vision | |
| CN206773731U (zh) | 基于图像处理的识别*** | |
| CN104730509A (zh) | 一种基于知识辅助置换检验的雷达检测方法 |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: No.32, Changqing street, Jiangning District, Nanjing City, Jiangsu Province 211153 Patentee after: China Shipbuilding Pengli (Nanjing) Atmospheric and Ocean Information System Co.,Ltd. Address before: No.32, Changqing street, Jiangning District, Nanjing City, Jiangsu Province 211153 Patentee before: CSIC PRIDE (NANJING) ATMOSPHERE MARINE INFORMATION SYSTEM Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |