CN110135559A - 一种基于船位数据获取张网网位和捕捞努力量的方法 - Google Patents
一种基于船位数据获取张网网位和捕捞努力量的方法 Download PDFInfo
- Publication number
- CN110135559A CN110135559A CN201910349833.1A CN201910349833A CN110135559A CN 110135559 A CN110135559 A CN 110135559A CN 201910349833 A CN201910349833 A CN 201910349833A CN 110135559 A CN110135559 A CN 110135559A
- Authority
- CN
- China
- Prior art keywords
- net
- fishing
- time
- accommodation point
- flight number
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000004308 accommodation Effects 0.000 claims abstract description 103
- 238000013528 artificial neural network Methods 0.000 claims abstract description 52
- 238000012549 training Methods 0.000 claims abstract description 21
- 238000012163 sequencing technique Methods 0.000 claims abstract description 12
- 239000000284 extract Substances 0.000 claims abstract description 8
- 238000012216 screening Methods 0.000 claims abstract description 7
- 210000002569 neuron Anatomy 0.000 claims description 41
- 238000003306 harvesting Methods 0.000 claims description 10
- FWQHNLCNFPYBCA-UHFFFAOYSA-N fluoran Chemical compound C12=CC=CC=C2OC2=CC=CC=C2C11OC(=O)C2=CC=CC=C21 FWQHNLCNFPYBCA-UHFFFAOYSA-N 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 239000004744 fabric Substances 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 description 21
- 239000013535 sea water Substances 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 238000003062 neural network model Methods 0.000 description 5
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 4
- 239000002775 capsule Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 210000005036 nerve Anatomy 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/06—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
- G06N3/061—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using biological neurons, e.g. biological neurons connected to an integrated circuit
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Business, Economics & Management (AREA)
- Neurology (AREA)
- Algebra (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Databases & Information Systems (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
Abstract
本发明涉及一种基于船位数据获取张网网位和捕捞努力量的方法，包括以下步骤：提取每个船位点的速度、偏角、时间顺序上相邻船位点的距离的特征作为训练集，输入到BP神经网络中进行训练；将陌生渔船的完整航次的每个船位点数据采用上述相同的方式进行预处理，输入到已训练好的BP神经网络中，通过BP神经网络判断每个船位点可能的作业状态；从经过筛选后的行驶状态和收网具状态的船位点中，确定张网渔船放网的开始时间、结束时间和网位坐标，计算每幅张网网具的放网时长和捕捞努力量，统计航次中全部网次的捕捞努力量获得该航次的累计捕捞努力量。本发明能够适用于从北斗船位数据中提取张网渔船每航次作业网位坐标和捕捞努力量。
Description
技术领域
本发明涉及渔船船位数据应用技术领域，特别是涉及一种基于船位数据获取张网网位和捕捞努力量的方法。
背景技术
渔船监测系统(VMS，Vessel Monitoring System)能够提供渔船时间、位置和动态信息，目前我国渔船安装船舶自动识别设备(AIS)近5万艘、北斗卫星船位监控近5万艘，CDMA公众移动通讯设备11万艘，初步实现对海洋渔船船位的实时联络及跟踪监控。国际上研究现有渔船监测位置分析渔船捕捞状态，其方法可以概括为三种：(1)通过分析船速的变化判断捕捞状态；(2)通过分析轨迹变化判断捕捞状态；(3)通过分析船速、航向等特征数据组成向量判断捕捞状态。由于北斗卫星传送的经纬度船位数据的时间分辨率为3分钟记录一次，空间分辨率约为10米，时空精度高，实时性强。已有的方法不适用于高时间分辨率的北斗卫星的船位数据分析，目前没有专门针对北斗卫星的渔船作业网次位置的判断方法。需要找到针对特定网具的计算方法，用来判断基于高时间分辨率的北斗船位数据的渔船作业网位并计算捕捞努力量。
发明内容
本发明所要解决的技术问题是提供一种基于船位数据获取张网网位和捕捞努力量的方法，能够适用于从北斗船位数据中提取张网渔船每航次作业网位坐标和捕捞努力量。
本发明解决其技术问题所采用的技术方案是：提供一种基于船位数据获取张网网位和捕捞努力量的方法，包括以下步骤：
(1)将已标记作业状态的完整的作业航次的船位点数据进行预处理，提取每个船位点的速度、偏角、时间顺序上相邻船位点的距离的特征作为训练集，输入到BP神经网络中进行训练；
(2)将陌生渔船的完整航次的每个船位点数据采用与步骤(1)相同的方式进行预处理，输入到已训练好的BP神经网络中，通过BP神经网络判断每个船位点可能的作业状态；
(3)从经过筛选后的行驶状态和收网具状态的船位点中，确定张网渔船放网的开始时间、结束时间和网位坐标，计算每幅张网网具的放网时长和捕捞努力量，统计航次中全部网次的捕捞努力量获得该航次的累计捕捞努力量。
所述步骤(1)中已标记的作业状态包括行驶状态、布网状态、收渔获状态、收网具状态和抛锚状态。
所述步骤(1)中的预处理具体为：统计并计算每个船位点的前后速度变化和速度变化率、行驶偏转角和向量偏转角变化、时间顺序上相邻船位点之间的距离和距离变化率。
所述步骤(1)中的BP神经网络包括输入层、隐含层和输出层，其中，输入层的神经元至隐含层的神经元的计算方式为隐含层的神经元至输出层的神经元的计算方式为φ1和φ2为神经网络激活函数，均为tanh函数，Xi为输入层的神经元，hj为隐含层的神经元，Wij为输入层第i个神经元与隐含层第j个神经元之间的权值，Gjk为隐含层第j个神经元与输出层第k个神经元之间的权值，l为输入层神经元的个数，m为隐含层神经元的个数，bj与均为偏差。
所述步骤(2)和步骤(3)中包括以下步骤：设置时间顺序上相邻船位点之间的距离和向量偏转角的阈值，在BP神经网络已判断作业状态的船位点中提取行驶状态和收网具状态的船位点，并根据设置的阈值过滤神经网络判断错误的船位点；所述向量偏转角是在时间顺序上相邻的三个船位点Ni-1、Ni和Ni+1形成的向量和之间的夹角。
所述步骤(3)中以航次中最后一个连续行驶状态的船位点记录时间作为放网开始时间，以每网次第一个收网具状态的船位点坐标作为该网次的网位坐标，并以该船位点记录时间作为该网次的放网结束时间。
所述步骤(3)中以各网次的放网结束时间减放网开始时间，获得各网次的放网时长，并于该网具网口迎流面积相乘获得各网次的捕捞努力量，累计各网次的捕捞努力量获得张网渔船该航次的捕捞努力量。
有益效果
由于采用了上述的技术方案，本发明与现有技术相比，具有以下的优点和积极效果：本发明通过BP神经网络自主学习训练人工标记的部分航次不同作业状态的船位点特征，用于判断陌生航次的船位点的作业状态，通过阈值筛选的方式过滤神经网络判断状态错误的船位点，在通过筛选的船位点中确定每个航次的网位坐标并计算该航次的捕捞努力量。
附图说明
图1A是张网船155某日作业中各船位点速度和速度变化率的特征变化图；
图1B是张网船155某日作业中各船位点航向偏转角的特征变化图；
图1C是张网船155某日作业中各船位点向量偏转角的特征变化图；
图1D是张网船155某日作业中各船位点距离等特征变化图；
图2是船位数据预处理方法示意图；
图3是BP神经网络拓扑结构示意图；
图4是隐含层神经元个数与均方误差关系图。
具体实施方式
下面结合具体实施例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明的实施方式涉及一种基于船位数据获取张网网位和捕捞努力量的方法，包括以下步骤：将已标记作业状态的完整的作业航次的船位点数据进行预处理，提取每个船位点的速度、偏角、时间顺序上相邻船位点的距离的特征作为训练集，输入到BP神经网络中进行训练；将陌生渔船的完整航次的每个船位点数据采用上述相同的方式进行预处理，输入到已训练好的BP神经网络中，通过BP神经网络判断每个船位点可能的作业状态；从经过筛选后的行驶状态和收网具状态的船位点中，确定张网渔船放网的开始时间、结束时间和网位坐标，计算每幅张网网具的放网时长和捕捞努力量，统计航次中全部网次的捕捞努力量获得该航次的累计捕捞努力量。
下面以具体的实施案例进一步说明本发明。
每航次作业时，张网渔船航行至目标渔场后随即开始布网工作，以渔船抵达目标海域结束长时间行驶状态的最后一个船位点记录时间作为放网开始时间；渔船收网具作业时，首先捞起固定于水底的网锚，该阶段以近乎垂直的方式打捞网锚，以渔船各网次收网具状态的第一个船位点坐标作为该网次的网位坐标，以该船位点的记录时间，作为该网次的放网截至时间。各网次的放网截止时间与放网开始时间相减，计算出各网次的放网时长，各网次放网时长与该网具的网口迎流面积相乘获得该网次捕捞努力量，统计该航次各网次的捕捞努力量，即为该航次累计捕捞努力量。
浙江省某张网渔船(以下称“张网船155”)，2016年12月～2017年3月有5个连续航次。每个航次的作业过程可划分为行驶、布网、收渔获、收网具、抛锚等5个状态。从 5个航次中，随机选出4个航次的船位数据，根据人工经验标记每个船位点的作业状态，将以选出的航次数据进行预处理，以放大每个船位点的特征，将预处理的航次数据作为训练集数据输入到BP神经网络中，对照以标记的状态训练神经网络模型。将剩余的一个航次数据进行与训练集数据相同处理方法的数据预处理，并输入到已训练好的神经网络中进行判断，通过阈值筛选作业开始的船位点和各网次的网位点坐标并记录时间，计算各网次的捕捞时长和捕捞努力量，并统计该航次的捕捞努力量。
1.张网作业特点分析
张网渔船每个航次的捕捞努力量与该航次作业的网次数，每副网具网口迎流面积和每副网具的放网时间直接相关。张网渔船作业方式较为特殊，每艘渔船通常搭载5～10副网具，其作业原理是利用浮、沉子作为网具网口的垂直扩展装置，在海水潮流的作用下伸展全网，同时利用海水潮流迫使捕捞对象进入网中。所以张网渔船每个航次作业均在海洋大潮汐期间进行(农历每月初一至初五或十六至二十)。因此张网渔船捕捞努力量也与海水潮流的流速有间接关系，但海水潮流方向和速度大小具有不确定性，指标难以量化，不纳入捕捞努力量考虑范围。张网渔船每个航次作业开始时，从渔港航行至目标渔场；到达渔场后，随即进行布网工作，先后释放网具的网锚、网身完成布网工作；通常每12小时对网身末端网囊进行一次起网，将渔获打捞至船舱后把网囊重新沉入水中，渔船抛锚等待渔获入网；多日捕捞工作结束后回收网具，先后收回网锚、网口和网身；回收网具工作完成后，随即返回渔港，整个航次作业结束，共有行驶、布网、收渔获、收网具、抛锚五种作业状态。张网网具的特点是网锚沉于水底静止不动，网口和网身随海水潮流飘荡，依靠海水潮流压力，迫使鱼类进入网中。因此，以每副张网的网锚坐标作为每副网具的网位坐标较为合理。
张网渔船抵达渔场后随即开始布网工作，以渔船结束长时间行驶状态后的第一个船位点作为放网开始时间，并记录该点坐标。张网渔船收网具工作时，先打捞网锚，此时渔船基本以垂直姿态进行打捞，以收网具状态的第一个船位点坐标来表示网锚位置，且该时刻渔船已停止该网具的捕捞作业，可作为该网具的放网截至时间。考虑到渔船作业的随机性 (同一航次可能先后在多地开展捕捞作业，导致收网作业次数大于渔船搭载网具数量)，但由于同一航次异地作业海域相距较远(通常20nmi以上)，在渔船航行到下一个目标渔场时，同样可以记录到一个新的捕捞作业开始时间和坐标。计算每副网的放网时长时，选择距离该网位坐标最近的作业开始船位点，将网位坐标记录的时间与作业开始时间相减，即可获得每副网的放网时长，每副网具的放网时长与该网具的网口迎流面积的乘积即为该网次的捕捞努力量，统计整个每个网次的捕捞努力量，计算得到该航次的累计捕捞努力量。
2.张网渔船作业状态划分
张网渔船根据作业环节的不同，可划分为5种状态，随机分析张网船155于2017年某日的船位数据为例：往返于渔场与渔港之间和捕捞作业时巡航于各网具之间的行驶状态(如图1A-1D中13:34～19:53时段)；抵达目标渔场后进行的布网状态(如图1A-1D中 19:54～23:34时段)；打捞网囊的进行捕捞工作的收渔获状态(如图1A-1D中5:16～6:32时段)；每航次结束前的收网具状态(如图1A-1D中6:34～13:32时段)；等待渔获入网的抛锚状态(如图1A-1D中0:00～5:15时段)。其中，图1A为渔船各作业状态速度与速度变化率的变化曲线，图1B为各渔船状态航向偏转角变化曲线，图1C为渔船各状态的向量偏转角变化曲线，图1D为渔船各状态的在时间顺序上相邻船位点距离和距离变化率的变化曲线。
航速变化范围为0～6m/s；航向随机性大，变化范围为0～360°，渔船长时间静置时航向为0°，航向偏转角是除首个船位点外，任意船位点航向与前一个船位点航向之差的绝对值，变化范围为0°～360°；向量偏转角为除首尾两个船位点外的任意船位点与前后两个船位点形成向量夹角，变化范围为-180°～180°；相邻距离为除首个船位点外，从任意船位点到其时间顺序上相邻的前一个船位点之间平均每分钟移动距离多为0～400米。不同作业状态的航速、航向偏转角、向量偏转角和相邻距离变化有不同特征，相同作业状态具有相似性。
渔船航向偏转角以正东方向为0°，以顺时针方向旋转为正方向；向量偏转角以向逆时针偏转为正，向顺时针偏转为负。图1A-1D表明了渔船作业时不同状态的特征差异。当渔船处于行驶状态时，其航速多为2.5～6m/s，航速变化率多为-0.04～0.04；航向偏转角多为 -10°～10°；向量偏转角多为0°～10°；平均每分钟移动距离多为200～360m，平均每分钟移动距离变化率多为0～2。当渔船处于布网状态时，其航速多为0～2m/s，航速随时间的变化率多为-0.96～0.14；航向偏转角为0°～342°；向量偏转角为-179°～177°；平均每分钟移动距离多为0.5～50m，平均每分钟移动距离变化率多为0～15，少数极端情况下平均每分钟移动距离变化率可高达45～65。当渔船处于收渔获状态时，其航速为多0～2.5m/s，航速随时间变化率多为-0.8～0.7；航向偏转角多为0°～346°；向量偏转角多为-177°～175°；平均每分钟移动距离多为5～50m，平均每分钟移动距离变化率多为0.3～5。当渔船处于收网具状态时，其航速多为0～3m/s，速度变化率多为-0.84～0.37，航向偏转角多为0°～300°；向量偏转角多为-90°～90°；平均每分钟移动距离多为5～120m，平均每分钟移动距离变化率多为0.01～6。当渔船处于抛锚状态时，其航速多为0m/s，航速变化率多为-0.38～0.5；航向偏转角多为0°；向量偏转角多为-180°～180°；平均每分钟移动距离多为0～10m，平均每分钟移动距离变化率多为0～3。由此可根据渔船多种特征来判断渔船作业状态。
3.张网渔船船位数据预处理
张网船在2016年12月～2017年3月连续的5个航次，共28770条VMS数据，把这些数据作为分析和处理的数据集，人工给每个船位点添加状态标记，用于训练和测试神经网络。北斗VMS数据主要包含渔船经纬度坐标、时间、航速、航向等信息，数据输入神经网络训练和测试前，需要进行预处理，以放大每个船位点的数据特征，图2为部分数据预处理方法示意图。
图2中，Ni-2、Ni-1、Ni、Ni+1、Ni+2点位是北斗VMS数据中按时间顺序排列的连续的5个船位坐标点。以Ni点为例，需要分别计算Ni与Ni-2、Ni-1、Ni+1、Ni+2四点之间的欧拉距离(式1)和Ni Ni+1距离相对Ni Ni-1距离变化率(式2)；Ni-1点到Ni点、Ni点到 Ni+1点的航向偏转角(式3)和与 与 与之间的向量夹角(式4)，若相邻向量呈逆时针方向偏转，则记为向量正偏转角，若相邻向量呈顺时针方向偏转，则记为向量负偏转角；Ni点航速相对Ni-1点航速变化率、Ni点航速相对Ni+1点航速变化率(式5)。并于Ni-1点、Ni点、Ni+1点处的航速共同构成Ni点的15个特征属性。
(xn,yn)与(xm,ym)分别为相邻两船位点Nn和Nm的经纬度坐标，其中，x代表经度，y代表纬度，R为地球半径，以6371公里代入计算，距离单位为米。
Lnn+1为船位点Nn与船位点Nn+1之间的距离，Lnn-1为船位点Nn与船位点Nn-1之间的距离。
hn为船位点Nn处的航向，为相邻船位点Nn和Nn-1航向差的绝对值。
若相对顺时针偏转小于180°，则θ取正值；若相对逆时针偏转小于180°，则θ取负值。
其中，Vn+1与Vn分别为船位点Nn与Nn+1处的航速。
4.神经网络模型构建
BP神经网络目前已被已广泛应用于模式识别和信号处理等领域。本实施例中采用BP 神经网络学习已标记作业状态的各船位点特征，自主判断陌生航次的船位点作业状态，其拓扑结构如图3所示。
如图3所示，Xi代表输入层神经元，其数据为训练集或测试集的输入数据，hj代表隐含层神经元，Yk代表输出层神经元，Wij代表输入层第i个神经元于隐含层第j个神经元之间的权值，Gjk代表隐含层第j个神经元于输出层第k个神经元之间的权值。输入层神经元到隐含层神经元的计算表达式为式(6)，隐含层神经元到输出层神经元的计算表达式为式(7)。
其中，bj与均为偏差，φ1和φ2均为神经网络的激活函数，本实施例中，根据输入层和输出层的数据特点，φ1和φ2均选用tanh函数，其表达式为式(8)。
本实施例中，将每条VMS数据中经纬度坐标、航速、航向等数据进行预处理，使每个船位点有15个属性，将每个船位点的15个属性数据，按北斗VMS记录船位点的时间顺序，依次输入到神经网络输入层中。输出层为5个神经元，在训练时分别对应由人工标记的5种作业状态，用布尔类型标记，即目标状态为1，其余为0，其分类结果为行驶状态[1,0,0,0,0]；布网状态[0,1,0,0,0]；收渔获状态[0,0,1,0,0]；收网具状态[0,0, 0,1,0]；抛锚状态[0,0,0,0,1]。BP神经网络不断训练学习中，输出结果与期望之间的误差不断改变，经过若干次迭代训练，每次迭代都根据神经网络设置的学习率对各神经元之间的权值进行一定程度的调整，最后输出误差达到预设范围或迭代次数达到预设值即停止训练。
5.神经网络模型训练
本实施例使用MATLAB R2014a作为实验平台，将预处理后的15维数据内部进行归一化处理，使每个船位点的特征数据的数值均控制在[-1,1]，调用MATLAB中神经网络工具箱建立网络网络模型，使用共轭梯度下降法进行训练，设置学习率为0.02，最大均方误差0.001。隐含层神经元个数以经验公式其中，h为BP神经网路中隐含层神经元个数，m为输入层神经元个数，n为输出层神经元个数，a为1～10之间整数为调节常数。在本实施例中，输入层神经元数为15，输出层神经元数位5，隐含层神经元个数最佳取值范围为5～15个。任意选择张网船155船5个航次中若干个航次的VMS数据，并对数据进行预处理，输入到已搭建的BP神经网络中训练10次，保存每次训练均方误差，对10 次均方误差取平均，选择平均均方误差最小的隐含层神经元个数作为本研究最终的BP神经网络模型。图4表示了BP神经网络输出结果的均方误差与隐含层神经元个数之间的关系直方图。当选择13个隐含层神经元时，10次训练的平均均方误差最小，因此确定隐含层神经元数量为13。
6.网位坐标和计算放网时长获取
张网渔船每航次作业过程中，先快速行驶至目标渔场，通常行驶时间为10～36小时。以渔船抵达渔场后的第一个船位点的记录时间作为放网开始时间，以渔船收网具状态的第一个船位点坐标作为该网次的网位坐标，并以该船位点的记录时间作为该网次的放网截至时间。将放网截至时间与开始时间相减，可计算每副网的放网时长时。每个网次的放网时长与该网次网具的网口迎流面积的向乘，获得该网次在该航次中的捕捞努力量，统计每航次中的所有网次的捕捞努力量，可得到该航次的累计捕捞努力量。
本实施例以张网船155于2016年12月～2017年3月的5个航次共28770条VMS数据作为分析数据集，记2016年12月9日～2016年12月19日为第一航次；记2016年12月 28日～2017年1月3日为第二航次；记2017年1月9日～2017年1月17日为第三航次；记2017年2月11日～2017年2月17日为第四航次；记2017年2月24日～2017年3月4 日为第五航次。使用人工标记每个船位点的状态，通过若干航次的VMS数据训练BP神经网络，使BP神经网络具有较高的船位状态判断能力，抽取其中的行驶状态和收网具状态，采用阈值降噪的方法排除由于神经网络错误标记而产生的干扰，从而确定张网渔船每个航次作业的网位坐标并计算每副网具的放网时长，通过已记录的渔船搭载网具网口迎流面积大小，计算出该渔船每个网次的捕捞努力量，进而统计每个航次的渔船累计捕捞努力量。
本实施例采用的阈值降噪方法具体为：设置时间顺序上相邻船位点之间的距离和向量偏转角的阈值，在BP神经网络已判断作业状态的船位点中提取行驶状态和收网具状态的船位点，并根据设置的阈值过滤神经网络判断错误的船位点。其中，所述向量偏转角是在时间顺序上相邻的三个船位点Ni-1、Ni和Ni+1形成的向量和之间的夹角。
将本实施例中采用不同训练集训练并搭建BP神经网络进行实验，并将神经网络提取出的测试集中每个航次作业的网位坐标、放网时长和累计捕捞努力量展示如下。表1～表4 分别为张网船155第二航次、第三航次、第四航次、第五航次作业中的网位坐标和每副网的放网时长的统计。张网船155每个航次携带10副网具，每副网具的网口迎流面积为200m2，结合各航次中每副网具的放网时长，计算该网次的捕捞努力量和该航次累计捕捞努力量(m2h)。
在4个航次的船位数据测试中，第二航次丢失一个网次数据，第三航次丢失两个网次数据，这是因为神经网络在这两个航次中的判断有一定失误，阈值降噪时将判断失误的船位点一并过滤，导致计算出的结果比实际情况少个别网次数据；第四航次和第五航次中，网次数大于该船实际搭载的网具数，是因为这两个航次都曾先后异地作业，导致网次数量较大。
本实施例中，主要用到BP神经网络自主判断的行驶和收网具两种状态，表5为结合BP神经网络在第二、第三、第四、第五航次的状态判断中，航行状态和收网具状态的精确度混淆矩阵表，其中矩阵横向表示BP神经网络判断各个船位点的分类和正确分类所占比例；矩阵纵向表示各个状态的船位点的真实数量和BP神经网络分类结果与精确度。BP神经网络总体分类准确度由矩阵对角线上分类正确的船位点总数与测试集总体样本数之间的比值确定。BP神经网络在这4个航次中的判断准确度分别为89.83％、92.06％、91.88％和89.39％。表6为统计第二、第三、第四、第五个航次的累计捕捞努力量。
表1第二航次作业网位、放网时长和捕捞努力量
表2第三航次作业网位、放网时长和捕捞努力量
表3第四航次作业网位、放网时长和捕捞努力量
表4第五航次作业网位、放网时长和捕捞努力量
表5神经网络判断第二、第三、第四、第五航次各状态的混淆矩阵
表6第二、第三、第四、第五航次的捕捞努力量统计
| 航次 | 网次数 | 累计捕捞时长(时) | 累计捕捞努力量(m<sup>2</sup>h) |
| 第二航次 | 9 | 1126.89 | 225378 |
| 第三航次 | 8 | 799.51 | 159902 |
| 第四航次 | 10 | 931.13 | 186226 |
| 第五航次 | 20 | 1592.74 | 318548 |
7.张网网位和捕捞努力量获取方法分析
本实施例使用BP神经网络模型，对人工判断标记的船位点状态进行学习，训练神经网络模型，使网络可以依据陌生船位点的特征判断其状态，并提取神经网络判断的行驶状态和收网具状态的船位点，通过阈值处理和筛选每个航次作业的网为坐标，并计算每个网次的放网时长，进而统计渔船各航次的累计捕捞努力量。在本实施例分析的4个航次中，作业时间最长的航次是2017年月24日～2017年3月4日，该航次作业时长200.57小时，累计放网时长1592.74小时，捕捞努力量318548，日均放网时长176.97小时，因为该航次作业中曾进行异地捕捞，导致作业网次较多，使该航次成为本研究分析的4个航次中捕捞努力量最大。2016年12月28日～2017年1月3日航次中，累计放网时长799.51小时，累计捕捞努力量115734，日均放网时长88.83小时，是4个航次中捕捞努力量最小的一个航次，由于该航次出海时间正处于东海寒潮期，海上风浪较大，渔船寻找合适的布网位置用了较长时间，但由于该船的续航能力和生活维持能力有限，不得不压缩该航次的捕捞时间，导致累计放网时长较短，进而影响捕捞努力量。
不难发现，本发明利用北斗VMS数据对传统作业方式中的张网渔船作业行为进行判断，并利用判断的结果提取渔船每航次的网位坐标并计算每副网的放网时长，从而量化统计渔船每个航次的捕捞努力量。
Claims (7)
1.一种基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，包括以下步骤：
(1)将已标记作业状态的完整的作业航次的船位点数据进行预处理，提取每个船位点的速度、偏角、时间顺序上相邻船位点的距离的特征作为训练集，输入到BP神经网络中进行训练；
(2)将陌生渔船的完整航次的每个船位点数据采用与步骤(1)相同的方式进行预处理，输入到已训练好的BP神经网络中，通过BP神经网络判断每个船位点可能的作业状态；
(3)从经过筛选后的行驶状态和收网具状态的船位点中，确定张网渔船放网的开始时间、结束时间和网位坐标，计算每幅张网网具的放网时长和捕捞努力量，统计航次中全部网次的捕捞努力量获得该航次的累计捕捞努力量。
2.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(1)中已标记的作业状态包括行驶状态、布网状态、收渔获状态、收网具状态和抛锚状态。
3.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(1)中的预处理具体为：统计并计算每个船位点的前后速度变化和速度变化率、行驶偏转角和向量偏转角变化、时间顺序上相邻船位点之间的距离和距离变化率。
4.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(1)中的BP神经网络包括输入层、隐含层和输出层，其中，输入层的神经元至隐含层的神经元的计算方式为隐含层的神经元至输出层的神经元的计算方式为φ1和φ2为神经网络激活函数，均为tanh函数，Xi为输入层的神经元，hj为隐含层的神经元，Wij为输入层第i个神经元与隐含层第j个神经元之间的权值，Gjk为隐含层第j个神经元与输出层第k个神经元之间的权值，l为输入层神经元的个数，m为隐含层神经元的个数，bj与均为偏差。
5.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(2)和步骤(3)中包括以下步骤：设置时间顺序上相邻船位点之间的距离和向量偏转角的阈值，在BP神经网络已判断作业状态的船位点中提取行驶状态和收网具状态的船位点，并根据设置的阈值过滤神经网络判断错误的船位点；所述向量偏转角是在时间顺序上相邻的三个船位点Ni-1、Ni和Ni+1形成的向量和之间的夹角。
6.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(3)中以航次中最后一个连续行驶状态的船位点记录时间作为放网开始时间，以每网次第一个收网具状态的船位点坐标作为该网次的网位坐标，并以该船位点记录时间作为该网次的放网结束时间。
7.根据权利要求1所述的基于船位数据获取张网网位和捕捞努力量的方法，其特征在于，所述步骤(3)中以各网次的放网结束时间减放网开始时间，获得各网次的放网时长，并于该网具网口迎流面积相乘获得各网次的捕捞努力量，累计各网次的捕捞努力量获得张网渔船该航次的捕捞努力量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910349833.1A CN110135559A (zh) | 2019-04-28 | 2019-04-28 | 一种基于船位数据获取张网网位和捕捞努力量的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910349833.1A CN110135559A (zh) | 2019-04-28 | 2019-04-28 | 一种基于船位数据获取张网网位和捕捞努力量的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110135559A true CN110135559A (zh) | 2019-08-16 |
Family
ID=67575603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910349833.1A Pending CN110135559A (zh) | 2019-04-28 | 2019-04-28 | 一种基于船位数据获取张网网位和捕捞努力量的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110135559A (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110796288A (zh) * | 2019-09-29 | 2020-02-14 | 宁波海上鲜信息技术有限公司 | 一种信息推送方法、装置及存储介质 |
| CN111652753A (zh) * | 2020-04-16 | 2020-09-11 | 中国水产科学研究院东海水产研究所 | 一种磷虾渔业中适宜连续捕捞渔场位置的多参数识别方法 |
| CN112434465A (zh) * | 2020-11-19 | 2021-03-02 | 江苏省海洋水产研究所 | 基于船位数据提取毛虾网有效布网网长的方法 |
| CN113360544A (zh) * | 2021-06-04 | 2021-09-07 | 中国海洋大学 | 基于捕捞时序关系的短时捕捞努力量分布预测方法 |
| CN113627359A (zh) * | 2021-08-13 | 2021-11-09 | 上海磐启微电子有限公司 | 一种基于双向长短时神经网络的渔船行为智能监控算法 |
| CN113887562A (zh) * | 2021-09-06 | 2022-01-04 | 江苏省海洋水产研究所 | 一种基于船位数据提取毛虾网船捕捞行为特征的方法 |
| CN113987838A (zh) * | 2021-11-19 | 2022-01-28 | 上海海洋大学 | 一种基于船位数据提取灯光敷网网次数的方法 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090105952A1 (en) * | 2007-10-19 | 2009-04-23 | Grace Ted V | Watercraft automation and aquatic effort data utilization |
| CN103678902A (zh) * | 2013-12-05 | 2014-03-26 | 中国水产科学研究院东海水产研究所 | 一种北斗船位数据提取拖网捕捞状态的方法 |
| CN103699781A (zh) * | 2013-12-05 | 2014-04-02 | 中国水产科学研究院东海水产研究所 | 一种北斗船位追溯拖网捕捞的方法 |
| CN104463433A (zh) * | 2014-11-19 | 2015-03-25 | 中国水产科学研究院东海水产研究所 | 一种基于渔船船位数据提取刺网捕捞状态的方法 |
| CN104484826A (zh) * | 2014-11-14 | 2015-04-01 | 中国水产科学研究院东海水产研究所 | 基于北斗船位数据提取拖网捕捞努力量时空特征的方法 |
| EP2910120A1 (en) * | 2014-02-25 | 2015-08-26 | Simrad Spain, S.L. | Method and system for controlling the fishing effort of a craft in a trawling fleet |
| CN105279702A (zh) * | 2015-09-23 | 2016-01-27 | 中国水产科学研究院东海水产研究所 | 一种基于北斗卫星数据拖网捕捞累计值获取方法 |
| CN106295833A (zh) * | 2015-05-18 | 2017-01-04 | 上海海洋大学 | 一种太平洋褶柔鱼资源补充量预测方法及其应用 |
| CN108170737A (zh) * | 2017-12-15 | 2018-06-15 | 华侨大学 | 一种渔业产值查询分析方法及系统 |
| CN108427824A (zh) * | 2018-01-29 | 2018-08-21 | 中国水产科学研究院东海水产研究所 | 一种利用船位数据提取刺网网长与角度的方法 |
| CN109188464A (zh) * | 2018-07-20 | 2019-01-11 | 中国水产科学研究院东海水产研究所 | 一种基于北斗船位数据的渔船在港和出海时间的分析方法 |
-
2019
- 2019-04-28 CN CN201910349833.1A patent/CN110135559A/zh active Pending
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090105952A1 (en) * | 2007-10-19 | 2009-04-23 | Grace Ted V | Watercraft automation and aquatic effort data utilization |
| CN103678902A (zh) * | 2013-12-05 | 2014-03-26 | 中国水产科学研究院东海水产研究所 | 一种北斗船位数据提取拖网捕捞状态的方法 |
| CN103699781A (zh) * | 2013-12-05 | 2014-04-02 | 中国水产科学研究院东海水产研究所 | 一种北斗船位追溯拖网捕捞的方法 |
| EP2910120A1 (en) * | 2014-02-25 | 2015-08-26 | Simrad Spain, S.L. | Method and system for controlling the fishing effort of a craft in a trawling fleet |
| CN104484826A (zh) * | 2014-11-14 | 2015-04-01 | 中国水产科学研究院东海水产研究所 | 基于北斗船位数据提取拖网捕捞努力量时空特征的方法 |
| CN104463433A (zh) * | 2014-11-19 | 2015-03-25 | 中国水产科学研究院东海水产研究所 | 一种基于渔船船位数据提取刺网捕捞状态的方法 |
| CN106295833A (zh) * | 2015-05-18 | 2017-01-04 | 上海海洋大学 | 一种太平洋褶柔鱼资源补充量预测方法及其应用 |
| CN105279702A (zh) * | 2015-09-23 | 2016-01-27 | 中国水产科学研究院东海水产研究所 | 一种基于北斗卫星数据拖网捕捞累计值获取方法 |
| CN108170737A (zh) * | 2017-12-15 | 2018-06-15 | 华侨大学 | 一种渔业产值查询分析方法及系统 |
| CN108427824A (zh) * | 2018-01-29 | 2018-08-21 | 中国水产科学研究院东海水产研究所 | 一种利用船位数据提取刺网网长与角度的方法 |
| CN109188464A (zh) * | 2018-07-20 | 2019-01-11 | 中国水产科学研究院东海水产研究所 | 一种基于北斗船位数据的渔船在港和出海时间的分析方法 |
Non-Patent Citations (8)
| Title |
|---|
| MILLS, C. M. ET AL: "Estimating high resolution trawl fishing effort from satellite-based vessel monitoring system data", 《ICES JOURNAL OF MARINE SCIENCE》 * |
| 张胜茂等: "北斗船位数据提取拖网捕捞努力量算法研究", 《水产学报》 * |
| 张胜茂等: "基于船位监控系统的拖网捕捞努力量提取方法研究", 《海洋科学》 * |
| 郑国富: "底拖网作业捕捞努力量标准化方法研究", 《福建水产》 * |
| 郑巧玲等: "基于神经网络和VMS的渔船捕捞类型辨别", 《南方水产科学》 * |
| 陈伟南等: "基于BP神经网络的船舶主机能效状态评估", 《中国船舰研究》 * |
| 黄冬梅等: "《海洋科学前沿系列丛书 海洋信息技术与应用》", 31 July 2016 * |
| 齐广瑞: "海州湾近岸张网的渔获结构及对渔业资源的影响", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110796288A (zh) * | 2019-09-29 | 2020-02-14 | 宁波海上鲜信息技术有限公司 | 一种信息推送方法、装置及存储介质 |
| CN110796288B (zh) * | 2019-09-29 | 2022-05-03 | 宁波海上鲜信息技术有限公司 | 一种信息推送方法、装置及存储介质 |
| CN111652753A (zh) * | 2020-04-16 | 2020-09-11 | 中国水产科学研究院东海水产研究所 | 一种磷虾渔业中适宜连续捕捞渔场位置的多参数识别方法 |
| CN112434465A (zh) * | 2020-11-19 | 2021-03-02 | 江苏省海洋水产研究所 | 基于船位数据提取毛虾网有效布网网长的方法 |
| CN112434465B (zh) * | 2020-11-19 | 2021-11-09 | 江苏省海洋水产研究所 | 基于船位数据提取毛虾网有效布网网长的方法 |
| CN113360544A (zh) * | 2021-06-04 | 2021-09-07 | 中国海洋大学 | 基于捕捞时序关系的短时捕捞努力量分布预测方法 |
| CN113627359A (zh) * | 2021-08-13 | 2021-11-09 | 上海磐启微电子有限公司 | 一种基于双向长短时神经网络的渔船行为智能监控算法 |
| CN113887562A (zh) * | 2021-09-06 | 2022-01-04 | 江苏省海洋水产研究所 | 一种基于船位数据提取毛虾网船捕捞行为特征的方法 |
| CN113987838A (zh) * | 2021-11-19 | 2022-01-28 | 上海海洋大学 | 一种基于船位数据提取灯光敷网网次数的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110135559A (zh) | 一种基于船位数据获取张网网位和捕捞努力量的方法 | |
| Gunderson | Surveys of fisheries resources | |
| CN110706827B (zh) | 基于船舶ais大数据的通航水域水流信息的提取方法及系统 | |
| Zimmerman et al. | Patterns of distribution of sound‐scattering zooplankton in warm‐and cold‐core eddies in the Gulf of Mexico, from a narrowband acoustic Doppler current profiler survey | |
| Satoh | Horizontal and vertical distribution of larvae of Pacific bluefin tuna Thunnus orientalis in patches entrained in mesoscale eddies | |
| Prozorkevich et al. | Survey report from the joint Norwegian/Russian ecosystem survey in the Barents Sea and adjacent waters, August-October 2017 | |
| KR20220025537A (ko) | 항적 데이타를 이용한 어선 조업 판별 방법 | |
| CN106097320A (zh) | 水下海参图像自动分割方法及装置 | |
| Stokesbury et al. | Impact of limited short-term sea scallop fishery on epibenthic community of Georges Bank closed areas | |
| Doray et al. | ICES survey protocols–Manual for acoustic surveys coordinated under ICES working group on acoustic and egg surveys for small pelagic fish (WGACEGG) | |
| Liao et al. | Research on intelligent damage detection of far-sea cage based on machine vision and deep learning | |
| Lough | Larval fish trophodynamic studies on Georges Bank: sampling strategy and initial results | |
| van der Meeren et al. | Survey report from the joint Norwegian/Russian ecosystem survey in the Barents Sea and adjacent waters, august-November 2020 | |
| CN112434465A (zh) | 基于船位数据提取毛虾网有效布网网长的方法 | |
| Butler et al. | The use of multibeam sonar mapping techniques to refine population estimates of the endangered white abalone (Haliotis sorenseni) | |
| Feng et al. | The study of identification of fishing vessel behavior based on VMS data | |
| CN112766308A (zh) | 一种远洋船舶分类方法 | |
| Kristmundsson et al. | Fish monitoring in aquaculture using multibeam echosounders and machine learning. | |
| Hassel et al. | Reaction of sandeel to seismic shooting: A field experiment and fishery statistics study | |
| Asher | A deeper look at Hawaiian coral reef fish assemblages: a comparison of survey approaches and assessments of shallow to mesophotic communities | |
| Yeh et al. | Influence of velocity and types of beam trawl towing on deep-sea demersal fish and decapod crustacean samples | |
| CN113887562B (zh) | 一种基于船位数据提取毛虾网船捕捞行为特征的方法 | |
| Xue et al. | Busy fishing area recognition based on improved K-means with random walking centroid | |
| Zhang et al. | Fishing vessel type recognition based on ship position data | |
| Zeinali et al. | CPUE, CPUA and distribution patterns of four demersal fishes in coastal waters of the northern Persian Gulf, Iran |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190816 |