CN112418521A - 一种短时海洋鱼群和鱼量的预测方法 - Google Patents
一种短时海洋鱼群和鱼量的预测方法 Download PDFInfo
- Publication number
- CN112418521A CN112418521A CN202011316548.9A CN202011316548A CN112418521A CN 112418521 A CN112418521 A CN 112418521A CN 202011316548 A CN202011316548 A CN 202011316548A CN 112418521 A CN112418521 A CN 112418521A
- Authority
- CN
- China
- Prior art keywords
- fish
- data
- fishing
- points
- grid
- 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
- 241000251468 Actinopterygii Species 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 32
- 150000003839 salts Chemical class 0.000 claims abstract description 20
- 238000012549 training Methods 0.000 claims description 29
- 208000037170 Delayed Emergence from Anesthesia Diseases 0.000 claims description 18
- 238000012360 testing method Methods 0.000 claims description 16
- 238000011176 pooling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- AGVAZMGAQJOSFJ-WZHZPDAFSA-M cobalt(2+);[(2r,3s,4r,5s)-5-(5,6-dimethylbenzimidazol-1-yl)-4-hydroxy-2-(hydroxymethyl)oxolan-3-yl] [(2r)-1-[3-[(1r,2r,3r,4z,7s,9z,12s,13s,14z,17s,18s,19r)-2,13,18-tris(2-amino-2-oxoethyl)-7,12,17-tris(3-amino-3-oxopropyl)-3,5,8,8,13,15,18,19-octamethyl-2 Chemical compound [Co+2].N#[C-].[N-]([C@@H]1[C@H](CC(N)=O)[C@@]2(C)CCC(=O)NC[C@@H](C)OP(O)(=O)O[C@H]3[C@H]([C@H](O[C@@H]3CO)N3C4=CC(C)=C(C)C=C4N=C3)O)\C2=C(C)/C([C@H](C\2(C)C)CCC(N)=O)=N/C/2=C\C([C@H]([C@@]/2(CC(N)=O)C)CCC(N)=O)=N\C\2=C(C)/C2=N[C@]1(C)[C@@](C)(CC(N)=O)[C@@H]2CCC(N)=O AGVAZMGAQJOSFJ-WZHZPDAFSA-M 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Images
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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- 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
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Evolutionary Computation (AREA)
- Human Resources & Organizations (AREA)
- Artificial Intelligence (AREA)
- Strategic Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Mining & Mineral Resources (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Computing Systems (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Agronomy & Crop Science (AREA)
- Marine Sciences & Fisheries (AREA)
- Software Systems (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
Abstract
本发明公开了一种短时海洋鱼群和鱼量的预测方法,通过挖掘渔船轨迹数据,从而获取渔船进行捕鱼活动的轨迹点,利用3D‑DBSCAN聚类算法,将时间因素考虑在内,从而避免了一些错误的聚类,可以在不划分渔船单次出海轨迹的情况下完成捕鱼点的提取,获得的捕鱼点更加准确。将海域进行栅格区域划分,将取得的捕鱼轨迹点在该海域进行投影,同时以渔船在每一栅格区域内的作业总时间来表示该区域鱼量的多少。本发明选取的残差网络这一特殊的卷积网络,通过加入海洋温、盐数据,使得预测的准确度进一步提高,最终得到短时海洋鱼群和鱼量的预测模型,利用短时海洋鱼群和鱼量的预测模型就可以对短时间内的栅格区域内是否有鱼群及鱼群数量进行预测。
Description
技术领域
本发明涉及海洋渔业领域,具体涉及一种短时海洋鱼群和鱼量的预测方法。
背景技术
海洋渔业资源在人类生存和发展中起着不可替代的作用,是现代社会主要经济活动的一部分。但随着捕捞区域的扩大和超出渔业资源再生速度的高强度捕捞,近海渔业资源环境遭到严重破坏,海洋渔业资源岌岌可危。为了维持渔业资源和环境的可持续发展,其中最重要的一点就是实现渔业资源的合理配置,保障适度捕捞。而未来渔业资源的合理预测则对于进行以上工作具有重要的意义。
预测未来的鱼群的位置与鱼量,需要有关的历史数据,但是由于海洋中的传感器的数量相对较少,无法获得海洋中准确的鱼群位置与数量,且每艘渔船的捕捞量对于各公司也都是保密的,这给鱼群预测带来了极大的挑战。因此,以往的研究仅从长远的角度探索海洋渔业资源的变化趋势与鱼群的迁徙状况,但该视角的时间跨度很长,往往以年为单位,对于现时的指导意义并不大。此外,有研究者通过渔船的行驶轨迹来预测未来捕鱼区域。从公开的AIS渔船轨迹数据集着手,需要从轨迹中尽可能准确地提取出渔船作业区域。目前,对于渔船作业状态(停泊、行驶、捕鱼)的识别多采用基于规则的方法(如基于速度阈值、角度阈值的方法)和基于密度的聚类方法(如DBSCAN算法)。由于渔船在驶离港口时载有燃油、食物等,而在返回港口时又满载鱼类,这就导致渔船的行驶速度较为缓慢,基于速度的方法很难过滤掉近港区域的轨迹点。渔船在进行捕鱼时,渔船的行驶方向往往会不断变化,基于角度阈值的方法可以获取捕鱼区域,但该方法需设置多个参数,如角度阈值,最小持续时间(过滤掉角度变化大,但并没有持续很长时间的轨迹点),最大容忍时间(两个大于角度阈值的轨迹点的时间差若小于该容忍时间,则将这两点之间的所有轨迹点置为作业状态点),参数的设置会对结果产生重大的影响。传统的DBSCAN聚类算法在单次出海轨迹聚类中取得了较好的效果,但该方法需要首先划分渔船的单次出海轨迹。由于传统的DBSCAN算法没有将时间考虑在内,若是选取的时间段内包含多次渔船出海轨迹,那么就会导致港口附近轨迹密集,从而形成错误的聚类。
发明内容
针对现有技术存在的海洋短时鱼群和鱼量无法有效预测的问题,本发明提供了一种短时海洋鱼群和鱼量的预测方法。
本发明采用以下的技术方案:
一种短时海洋鱼群和鱼量的预测方法,包括以下步骤:
步骤1:选取海洋中要预测的区域,获得该区域的所有渔船轨迹数据并进行预处理;
步骤2:利用3D-DBSCAN聚类算法对预处理后的渔船轨迹数据进行聚类获得捕鱼点,将捕鱼点投影到栅格矩阵进行数据的格式转换;
步骤3:搭建残差网络模型,将转换后的栅格矩阵数据分为训练集和测试集,将训练集送入残差网络模型进行模型的训练;
步骤4:测试集送入训练好的残差网络模型进行验证,最终得到短时海洋鱼群和鱼量的预测模型。
优选地,步骤1中的所有渔船轨迹数据为所有渔船完整的轨迹数据。
优选地,3D-DBSCAN聚类算法对预处理后的渔船轨迹数据进行聚类获得捕鱼点的具体过程为:
步骤2.1:设定3D-DBSCAN聚类算法用参数:eps、minPts和t,其中,eps为扫描半径,minPts为最小包含点数,t为扫描最小时间间隔;
步骤2.2:对一条渔船在单位时间长度内的渔船轨迹数据进行聚类,任选一个未被访问的点为出发点,找出与出发点距离在小于等于eps之内的所有附近点,如果附近点的数量≥minPts,且附近点与出发点的时间间隔小于t,则出发点与其附近点形成一个簇,并且出发点被标记为已访问,然后递归,以相同的方法处理该簇内所有未被标记为已访问的点,从而对簇进行扩展;如果簇充分地被扩展,即簇内的所有点被标记为已访问,然后用同样的算法去处理簇外未被访问的点,最终形成多个有用簇;
步骤2.3:对所有渔船在单位时间长度内的渔船轨迹数据利用步骤2.2进行聚类,得到所有渔船轨迹数据聚类后的有用簇,有用簇内的点即为捕鱼点。
优选地,将预测的区域划分为(i×j)个栅格,将获得的捕鱼点在(i×j)个栅格中进行投影,求得在第T天在第(i×j)栅格区域内渔船捕鱼的总时间间接代表栅格区域内的鱼量多少;进而得到在第T天所有栅格区域的鱼量多少情况XT;
进而得到总天数内的栅格区域的鱼量多少情况(X1,X2,…,XT)。
优选地,步骤3具备包括以下步骤:
步骤3.1:将总天数内的栅格区域的鱼量多少情况(X1,X2,…,XT)进行划分,每连续k天为一组,用每连续k天为一组的数据输入残差网络模型预测出接下来一天的数据,这样就能得到总的数据集为[T-k,k,i,j],其中选取80%作为训练集,20%作为测试集,将数据归一化到[-1,1]区间;
其中,*为卷积操作,W(1)和b(1)表示在第一层卷积中的权重和偏置系数;
外部因素获取:
由于只能获得每个栅格区域的四个角坐标点处的温、盐数据,所以将每个栅格的四个角坐标点处的温、盐数据经过一层均值池化层,得到四个角坐标中温、盐数据的平均值,用平均值来表示对应栅格区域的温、盐数据值;
E表示外部因素的输入,即栅格区域的温、盐数据值;
则经过均值池化层与卷积层后变为:
E(1)=f卷积(f均值池化(E))
步骤3.4:由于残差学习单元允许原始输入的信息直接传输到后层的执行机构,因此可以搭建深层网络;在这项工作中,每个残差块由卷积层和线性整流单元组成,第l个残差块的输出为:
Yl=f残差(Yl-1;θl)+Yl-1,l=1,…,L
其中f残差为残差函数,表示为两个线性整流单元和两个卷积层的叠加,θl表示第l个残差块的所有可训练参数,L为残差块的个数。Yl-1与Yl分别表示第l个残差块的输入与输出;
采用均方误差MSE来衡量预测值与真实值之间的差距;
步骤3.5:利用训练集和外部因素数据一直训练残差网络模型,直到残差网络模型中的参数达到最优,均方误差MSE得到的预测值与真实值之间的差距最小。
优选地,步骤4具体包括:
将测试集送入训练好的模型进行测试,测试完成后就能得到短时海洋鱼群和鱼量的预测模型:
其中,模型中的参数均已通过训练求得最优;
通过短时海洋鱼群和鱼量的预测模型就能在输入要预测天之前的连续k天的栅格区域的鱼量多少情况下,得到要预测天的栅格区域中的鱼量分布情况。
本发明具有的有益效果是:
本发明提供的短时海洋鱼群和鱼量的预测方法,通过挖掘渔船轨迹数据,从而获取渔船进行捕鱼活动的轨迹点,这里利用了全新的3D-DBSCAN聚类算法,将时间因素考虑在内,从而避免了一些错误的聚类,可以在不划分渔船单次出海轨迹的情况下完成捕鱼点的提取,获得的捕鱼点更加准确。原始的轨迹序列数据不能直观的表现出鱼群的时空分布状况,且不容易被神经网络提取时空依赖性。因此将海域进行栅格区域划分,将取得的捕鱼轨迹点在该海域进行投影,同时以渔船在每一栅格区域内的作业总时间来表示该区域鱼量的多少。这样原始的序列数据就转换为了栅格矩阵形式,转换后的数据更容易被神经网络所接受,并且能更加直观地表现鱼群的时空分布状况。在捕获空间依赖性上,卷积网络具有明显的优越性,同时为了搭建深层的网络,本发明选取的残差网络这一特殊的卷积网络,模型的结构简单,参数较少,训练、预测时间短且准确率高。同时,海洋温、盐会对鱼群的时空分布产生影响,本发明通过加入海洋温、盐数据,使得预测的准确度进一步提高,最终得到短时海洋鱼群和鱼量的预测模型,利用短时海洋鱼群和鱼量的预测模型就可以对短时间内的栅格区域内是否有鱼群及鱼群数量进行预测。
附图说明
图1为实施例1中一艘渔船一次出海的轨迹聚类图。
图2为残差网络模型的结构图。
具体实施方式
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
实施例1
结合图1和图2,一种短时海洋鱼群和鱼量的预测方法,包括以下步骤:
步骤1:选取海洋中要预测的区域,获得该区域的所有渔船轨迹数据并进行预处理。
上述所有渔船轨迹数据为所有渔船完整的轨迹数据。
在该实施例中,该预测的区域范围为:经度:121°E~128°E,纬度:28°N~34°N。
选取的是该预测的区域内所有的渔船1年时间内完整的轨迹数据。
步骤2:利用3D-DBSCAN聚类算法对预处理后的渔船轨迹数据进行聚类获得捕鱼点,将捕鱼点投影到栅格矩阵进行数据的格式转换。
3D-DBSCAN聚类算法的概念是在传统的DBCSAN算法中加入了时间参数,原因在于:
1、在一年时间范围内,一艘渔船往往有多次出海轨迹,这就会使港口附近轨迹重叠,导致轨迹点密度增加,从而形成错误的聚类。
2、一艘渔船的一次出海轨迹会存在轨迹重叠,导致出现错误的聚类。
下面用一个例子进行说明,如图1为一艘渔船一次出海的轨迹,利用传统的DBSCAN算法将eps设为圆的半径,minPts设置为3,可以得到两个停留点的聚类。但是会发现,聚类1是一个错误的聚类,原因在于仅仅考虑了点之间的平面距离,而没有考虑时间。聚类1中两个实心点与两个空心点在时间上并不是联系的。
而3D-DBSCAN聚类算法将时间考虑在内,给定扫描最小时间间隔t,若两点间的时间间隔大于t,且不能够通过它点连续可达,那么这两点将不能成簇。即聚类1中的实心点和空心点将不会被聚为一类,只有聚类2是正确的。
由此,3D-DBSCAN聚类算法就可以将一条渔船在单位时间长度内的渔船轨迹数据进行很好的聚类,去掉了港口出多次出海轨迹带来的错误聚类,也去掉了单次出海轨迹重叠带来的错误聚类,能在不划分渔船单次出海轨迹的情况下完成捕鱼点的提取。
具体包括:
步骤2.1:设定3D-DBSCAN聚类算法用参数:eps、minPts和t,其中,eps为扫描半径,minPts为最小包含点数,t为扫描最小时间间隔;
在本实施例中,由于渔船在正常行驶的情况下,速度大于10节/小时,即18.52千米/小时。轨迹数据的采样点时间间隔为1分钟,两个采样点之间的距离间隔大于308米。因此将这一距离值作为扫描半径eps。航行状态下,以308为扫描半径的区域圆中将最多包含3个点,而在渔船处于捕鱼状态时,轨迹点数会大于3,因此将3作为临界值,通过逐次增加点数与可视化结果分析得到最佳minPts等于7。对于时间间隔t的设置,由于渔船出海捕鱼往往以天为单位,同时渔船回港,卸载、整休再到下次出海也会间隔几天,因此将时间间隔设置为一天,即24小时。
步骤2.2:对一条渔船在1年时间长度内的渔船轨迹数据进行聚类,任选一个未被访问的点为出发点,找出与出发点距离在小于等于eps之内的所有附近点,如果附近点的数量≥minPts,且附近点与出发点的时间间隔小于t,则出发点与其附近点形成一个簇,并且出发点被标记为已访问,然后递归,以相同的方法处理该簇内所有未被标记为已访问的点,从而对簇进行扩展;如果簇充分地被扩展,即簇内的所有点被标记为已访问,然后用同样的算法去处理簇外未被访问的点,最终形成多个有用簇;
步骤2.3:对所有渔船在单位时间长度内的渔船轨迹数据利用步骤2.2进行聚类,得到所有渔船轨迹数据聚类后的有用簇,有用簇内的点即为捕鱼点。
捕鱼点获取后,为了描述鱼群位置与数量的时空分布,同时能更好的利用神将网络提取时空特征,本实施例中将预测的区域海域划分为(24×28)个栅格,每个栅格的精度为0.25°×0.25°。
进而得到总天数内的栅格区域的鱼量多少情况(X1,X2,…,XT)。
在本实施例中捕鱼的总天数为200,即T的最大取值为200。
步骤3:搭建残差网络模型,将转换后的栅格矩阵数据分为训练集和测试集,将训练集送入残差网络模型进行模型的训练。
步骤3.1:将总天数内的栅格区域的鱼量多少情况(X1,X2,…,XT)进行划分,每连续k天为一组,用每连续k天为一组的数据输入残差网络模型预测出接下来一天的数据,这样就能得到总的数据集为[T-k,k,i,j],其中选取80%作为训练集,20%作为测试集,为了加快训练,将数据归一化到[-1,1]区间。
在本实施例中,k=3,即以连续3天的鱼量多少情况数据为一组数据。例如,用第1到第3天的鱼量多少情况数据预测第4天的鱼量多少情况数据,用第2天到第4天的鱼量多少情况数据预测第5天的鱼量多少情况数据,依次类推,直到用第197天至199天的鱼量多少情况数据预测第200天的鱼量多少情况数据,总计197组数据。
上述每组数据有(3×24×28)的张量。
步骤3.2:设定残差网络模型的输入,残差网络模型的输入包括训练集和外部因素数据,外部因素数据包括海洋的温、盐数据。已有的研究表明,海洋温度、盐度会对鱼群的分布产生影响,因此本发明考虑了海洋温度和盐度两种外部因素。
其中,*为卷积操作,W(1)和b(1)表示在第一层卷积中的权重和偏置系数。
外部因素获取:
由于只能获得每个栅格区域的四个角坐标点处的温、盐数据,所以将每个栅格的四个角坐标点处的温、盐数据经过一层均值池化层,得到四个角坐标中温、盐数据的平均值,用平均值来表示对应栅格区域的温、盐数据值。
E表示外部因素的输入,即栅格区域的温、盐数据值;
则经过均值池化层与卷积层后变为:
E(1)=f卷积(f均值池化(E))
其中+为张量的拼接操作,即将(3×24×28)的轨迹数据与(2×24×28)的外部因素数据拼接为(5×24×28)的张量。
步骤3.4:由于残差学习单元允许原始输入的信息直接传输到后层的执行机构,因此可以搭建深层网络;在这项工作中,每个残差块由卷积层和线性整流单元(ReLU)组成,第l个残差块的输出为:
Yl=f残差(Yl-1;θl)+Yl-1,l=1,…,L
其中f残差为残差函数,表示为两个线性整流单元(ReLU)和两个卷积层的叠加,θl表示第l个残差块的所有可训练参数,L为残差块的个数。Yl-1与Yl分别表示第l个残差块的输入与输出。
因为需要的输出为(24×28)的张量,即接下来一天内的栅格区域的鱼量多少情况,Yl再经过只有一个卷积核的卷积层,最后再通过Tanh激活层将结果映射到[-1,1],这样就得到了连续k天后接下来一天的鱼量预测的结果
采用均方误差MSE来衡量预测值与真实值之间的差距;
步骤3.5:利用训练集和外部因素数据一直训练残差网络模型,直到残差网络模型中的参数达到最优,均方误差MSE得到的预测值与真实值之间的差距最小。
步骤4:测试集送入训练好的残差网络模型进行验证,最终得到短时海洋鱼群和鱼量的预测模型。
测试完成后就能得到短时海洋鱼群和鱼量的预测模型:
其中,模型中的参数均已通过训练求得最优;
通过短时海洋鱼群和鱼量的预测模型就能在输入要预测天之前的连续k天的栅格区域的鱼量多少情况下,得到要预测天的栅格区域中的鱼量分布情况。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (6)
1.一种短时海洋鱼群和鱼量的预测方法,其特征在于,包括以下步骤:
步骤1:选取海洋中要预测的区域,获得该区域的所有渔船轨迹数据并进行预处理;
步骤2:利用3D-DBSCAN聚类算法对预处理后的渔船轨迹数据进行聚类获得捕鱼点,将捕鱼点投影到栅格矩阵进行数据的格式转换;
步骤3:搭建残差网络模型,将转换后的栅格矩阵数据分为训练集和测试集,将训练集送入残差网络模型进行模型的训练;
步骤4:测试集送入训练好的残差网络模型进行验证,最终得到短时海洋鱼群和鱼量的预测模型。
2.根据权利要求1所述的一种短时海洋鱼群和鱼量的预测方法,其特征在于,步骤1中的所有渔船轨迹数据为所有渔船完整的轨迹数据。
3.根据权利要求1所述的一种短时海洋鱼群和鱼量的预测方法,其特征在于,3D-DBSCAN聚类算法对预处理后的渔船轨迹数据进行聚类获得捕鱼点的具体过程为:
步骤2.1:设定3D-DBSCAN聚类算法用参数:eps、minPts和t,其中,eps为扫描半径,minPts为最小包含点数,t为扫描最小时间间隔;
步骤2.2:对一条渔船在单位时间长度内的渔船轨迹数据进行聚类,任选一个未被访问的点为出发点,找出与出发点距离在小于等于eps之内的所有附近点,如果附近点的数量≥minPts,且附近点与出发点的时间间隔小于t,则出发点与其附近点形成一个簇,并且出发点被标记为已访问,然后递归,以相同的方法处理该簇内所有未被标记为已访问的点,从而对簇进行扩展;如果簇充分地被扩展,即簇内的所有点被标记为已访问,然后用同样的算法去处理簇外未被访问的点,最终形成多个有用簇;
步骤2.3:对所有渔船在单位时间长度内的渔船轨迹数据利用步骤2.2进行聚类,得到所有渔船轨迹数据聚类后的有用簇,有用簇内的点即为捕鱼点。
5.根据权利要求4所述的一种短时海洋鱼群和鱼量的预测方法,其特征在于,步骤3具备包括以下步骤:
步骤3.1:将总天数内的栅格区域的鱼量多少情况(X1,X2,…,XT)进行划分,每连续k天为一组,用每连续k天为一组的数据输入残差网络模型预测出接下来一天的数据,这样就能得到总的数据集为[T-k,k,i,j],其中选取80%作为训练集,20%作为测试集,将数据归一化到[-1,1]区间;
其中,*为卷积操作,W(1)和b(1)表示在第一层卷积中的权重和偏置系数;
外部因素获取:
由于只能获得每个栅格区域的四个角坐标点处的温、盐数据,所以将每个栅格的四个角坐标点处的温、盐数据经过一层均值池化层,得到四个角坐标中温、盐数据的平均值,用平均值来表示对应栅格区域的温、盐数据值;
E表示外部因素的输入,即栅格区域的温、盐数据值;
则经过均值池化层与卷积层后变为:
E(1)=f卷积(f均值池化(E))
步骤3.4:由于残差学习单元允许原始输入的信息直接传输到后层的执行机构,因此可以搭建深层网络;在这项工作中,每个残差块由卷积层和线性整流单元组成,第l个残差块的输出为:
Yl=f残差(Yl-1;θl)+Yl-1,l=1,…,L
其中f残差为残差函数,表示为两个线性整流单元和两个卷积层的叠加,θl表示第l个残差块的所有可训练参数,L为残差块的个数。Yl-1与Yl分别表示第l个残差块的输入与输出;
采用均方误差MSE来衡量预测值与真实值之间的差距;
步骤3.5:利用训练集和外部因素数据一直训练残差网络模型,直到残差网络模型中的参数达到最优,均方误差MSE得到的预测值与真实值之间的差距最小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011316548.9A CN112418521B (zh) | 2020-11-23 | 2020-11-23 | 一种短时海洋鱼群和鱼量的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011316548.9A CN112418521B (zh) | 2020-11-23 | 2020-11-23 | 一种短时海洋鱼群和鱼量的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112418521A true CN112418521A (zh) | 2021-02-26 |
CN112418521B CN112418521B (zh) | 2023-02-24 |
Family
ID=74777169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011316548.9A Active CN112418521B (zh) | 2020-11-23 | 2020-11-23 | 一种短时海洋鱼群和鱼量的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112418521B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113598131A (zh) * | 2021-07-24 | 2021-11-05 | 南京添益越科技有限公司 | 一种基于大数据区块链的内湖养殖系统 |
CN114004428A (zh) * | 2022-01-04 | 2022-02-01 | 中科曙光南京研究院有限公司 | 一种基于蜂窝网格集成机制与时空Transformer的渔情预测方法及系统 |
CN114065520A (zh) * | 2021-11-17 | 2022-02-18 | 西安理工大学 | 一种鱼类洄游通道确定方法及系统 |
CN114429309A (zh) * | 2022-02-08 | 2022-05-03 | 闫雪 | 基于人工智能的区域数据检测平台 |
CN114444819A (zh) * | 2022-04-11 | 2022-05-06 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种渔业资源预测方法、装置、存储介质以及电子设备 |
CN115860284A (zh) * | 2023-03-01 | 2023-03-28 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 渔业资源密度识别方法、装置、存储介质以及电子设备 |
CN116579508A (zh) * | 2023-07-13 | 2023-08-11 | 海煜(福州)生物科技有限公司 | 一种鱼类预测方法、装置、设备及存储介质 |
CN117576839A (zh) * | 2024-01-17 | 2024-02-20 | 深圳市雷铭科技发展有限公司 | 针对无电地区的渔场状态监测方法及相关装置 |
CN117576839B (zh) * | 2024-01-17 | 2024-05-10 | 深圳市雷铭科技发展有限公司 | 针对无电地区的渔场状态监测方法及相关装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108844A (zh) * | 2017-12-25 | 2018-06-01 | 儒安科技有限公司 | 一种城市人流量预测方法及系统 |
CN109448442A (zh) * | 2018-10-12 | 2019-03-08 | 天津大学 | 一种海域内多船碰撞危险度的高效监测方法 |
CN110033051A (zh) * | 2019-04-18 | 2019-07-19 | 杭州电子科技大学 | 一种基于多步聚类的拖网渔船行为判别方法 |
CN110310474A (zh) * | 2018-05-14 | 2019-10-08 | 桂林远望智能通信科技有限公司 | 一种基于时空残差网络的车流量预测方法及装置 |
CN111144666A (zh) * | 2020-01-02 | 2020-05-12 | 吉林大学 | 一种基于深度时空残差网络的海洋温跃层预测方法 |
CN111583052A (zh) * | 2020-05-21 | 2020-08-25 | 北京天合睿创科技有限公司 | 渔船轨迹追踪及渔情分析系统 |
CN111680734A (zh) * | 2020-06-02 | 2020-09-18 | 河海大学常州校区 | 基于改进型elm的海洋鱼群迁移路线预测方法及系统 |
CN111858565A (zh) * | 2020-06-10 | 2020-10-30 | 北京航空航天大学 | 一种基于线段聚类的多维度高精度航迹智能预测方法 |
EP3739295A1 (en) * | 2019-05-13 | 2020-11-18 | S.A.T.E. - Systems and Advanced Technologies Engineering S.R.L. | Method for determining an optimal naval navigation routes from historical gnss data of naval trajectories |
-
2020
- 2020-11-23 CN CN202011316548.9A patent/CN112418521B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108844A (zh) * | 2017-12-25 | 2018-06-01 | 儒安科技有限公司 | 一种城市人流量预测方法及系统 |
CN110310474A (zh) * | 2018-05-14 | 2019-10-08 | 桂林远望智能通信科技有限公司 | 一种基于时空残差网络的车流量预测方法及装置 |
CN109448442A (zh) * | 2018-10-12 | 2019-03-08 | 天津大学 | 一种海域内多船碰撞危险度的高效监测方法 |
CN110033051A (zh) * | 2019-04-18 | 2019-07-19 | 杭州电子科技大学 | 一种基于多步聚类的拖网渔船行为判别方法 |
EP3739295A1 (en) * | 2019-05-13 | 2020-11-18 | S.A.T.E. - Systems and Advanced Technologies Engineering S.R.L. | Method for determining an optimal naval navigation routes from historical gnss data of naval trajectories |
CN111144666A (zh) * | 2020-01-02 | 2020-05-12 | 吉林大学 | 一种基于深度时空残差网络的海洋温跃层预测方法 |
CN111583052A (zh) * | 2020-05-21 | 2020-08-25 | 北京天合睿创科技有限公司 | 渔船轨迹追踪及渔情分析系统 |
CN111680734A (zh) * | 2020-06-02 | 2020-09-18 | 河海大学常州校区 | 基于改进型elm的海洋鱼群迁移路线预测方法及系统 |
CN111858565A (zh) * | 2020-06-10 | 2020-10-30 | 北京航空航天大学 | 一种基于线段聚类的多维度高精度航迹智能预测方法 |
Non-Patent Citations (2)
Title |
---|
初晓等: "基于残差网络的海洋温跃层分析方法", 《吉林大学学报(理学版)》 * |
袁红春等: "基于空间自回归和空间聚类的渔情预报模型", 《海洋科学》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113598131A (zh) * | 2021-07-24 | 2021-11-05 | 南京添益越科技有限公司 | 一种基于大数据区块链的内湖养殖系统 |
CN114065520A (zh) * | 2021-11-17 | 2022-02-18 | 西安理工大学 | 一种鱼类洄游通道确定方法及系统 |
CN114065520B (zh) * | 2021-11-17 | 2022-07-05 | 西安理工大学 | 一种鱼类洄游通道确定方法及系统 |
CN114004428A (zh) * | 2022-01-04 | 2022-02-01 | 中科曙光南京研究院有限公司 | 一种基于蜂窝网格集成机制与时空Transformer的渔情预测方法及系统 |
CN114429309A (zh) * | 2022-02-08 | 2022-05-03 | 闫雪 | 基于人工智能的区域数据检测平台 |
CN114429309B (zh) * | 2022-02-08 | 2022-11-08 | 深圳市润信数据技术有限公司 | 基于人工智能的区域数据检测平台 |
CN114444819A (zh) * | 2022-04-11 | 2022-05-06 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种渔业资源预测方法、装置、存储介质以及电子设备 |
CN115860284A (zh) * | 2023-03-01 | 2023-03-28 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 渔业资源密度识别方法、装置、存储介质以及电子设备 |
CN116579508A (zh) * | 2023-07-13 | 2023-08-11 | 海煜(福州)生物科技有限公司 | 一种鱼类预测方法、装置、设备及存储介质 |
CN116579508B (zh) * | 2023-07-13 | 2023-10-13 | 海煜(福州)生物科技有限公司 | 一种鱼类预测方法、装置、设备及存储介质 |
CN117576839A (zh) * | 2024-01-17 | 2024-02-20 | 深圳市雷铭科技发展有限公司 | 针对无电地区的渔场状态监测方法及相关装置 |
CN117576839B (zh) * | 2024-01-17 | 2024-05-10 | 深圳市雷铭科技发展有限公司 | 针对无电地区的渔场状态监测方法及相关装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112418521B (zh) | 2023-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112418521B (zh) | 一种短时海洋鱼群和鱼量的预测方法 | |
CN109214107A (zh) | 一种船舶航行行为在线预测方法 | |
CN109409499A (zh) | 一种基于深度学习和卡尔曼滤波修正的轨迹恢复方法 | |
US20150094904A1 (en) | System and method for underwater vehicle simulation | |
Li et al. | Long-term vessel motion predication by modeling trajectory patterns with AIS data | |
Marie et al. | 1 Multi-objective optimization of motor vessel route | |
CN114842180B (zh) | 一种点云补全方法、装置、设备及介质 | |
CN114936332A (zh) | 一种风电场水域船舶交通流分析的方法、装置和存储介质 | |
Molchanov et al. | Active drifters: Towards a practical multi-robot system for ocean monitoring | |
CN113393027B (zh) | 一种基于深度学习的航标漂移智能预测方法 | |
CN114139642A (zh) | 联合船舶行程语义对象与图论理论的船舶航线提取方法 | |
CN112556698B (zh) | 船舶行为变化及异常识别方法、设备及介质 | |
CN117494871A (zh) | 一种考虑船舶交互影响的船舶轨迹预测方法 | |
Ware et al. | A system for cleaning high volume bathymetry | |
CN111612217A (zh) | 一种用于海上风电场巡检的水面无人艇改进ga-sa路径规划方法 | |
CN114565176B (zh) | 一种长期船舶轨迹预测方法 | |
CN116306790A (zh) | 一种基于cnn-gru和注意力机制的近海船舶轨迹实时预测方法、系统、设备及介质 | |
CN115859116A (zh) | 一种基于径向基函数回归插值法的海洋环境场重建方法 | |
Eljabu et al. | Destination port detection for vessels: An analytic tool for optimizing port authorities resources | |
KR20230001867A (ko) | 구조 구난을 위한 해상 표류물의 위치 추정 방법 및 장치 | |
Teng et al. | HIT: Hybrid interpolation for trawlers based on deep learning | |
CN113392884A (zh) | 基于lstm网络和注意力机制的坐标融合方法 | |
CN117076893B (zh) | 一种基于长短期记忆神经网络的声速分布预报方法 | |
Li et al. | DCSP-Yolov5: Improved Yolov5 Based on Dilated Convolution for Object Detection of Forward-Looking Sonar Images | |
CN113434617B (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 |