CN113159400A - 基于原位遥感信息的远洋渔类生境分布预测方法及其应用 - Google Patents
基于原位遥感信息的远洋渔类生境分布预测方法及其应用 Download PDFInfo
- Publication number
- CN113159400A CN113159400A CN202110359602.6A CN202110359602A CN113159400A CN 113159400 A CN113159400 A CN 113159400A CN 202110359602 A CN202110359602 A CN 202110359602A CN 113159400 A CN113159400 A CN 113159400A
- Authority
- CN
- China
- Prior art keywords
- cpue
- remote sensing
- sensing information
- prediction model
- time period
- 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
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000009826 distribution Methods 0.000 title claims abstract description 51
- 241000251468 Actinopterygii Species 0.000 claims abstract description 78
- 238000012360 testing method Methods 0.000 claims abstract description 57
- 238000012549 training Methods 0.000 claims abstract description 52
- 238000013528 artificial neural network Methods 0.000 claims abstract description 31
- 230000007613 environmental effect Effects 0.000 claims abstract description 12
- 238000004519 manufacturing process Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims description 10
- 238000012952 Resampling Methods 0.000 claims description 6
- 230000015654 memory Effects 0.000 claims description 6
- 229930002868 chlorophyll a Natural products 0.000 claims description 4
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 4
- 238000013480 data collection Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims 1
- 238000000605 extraction Methods 0.000 abstract 1
- 241000238366 Cephalopoda Species 0.000 description 8
- 230000006870 function Effects 0.000 description 8
- 238000004088 simulation Methods 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 238000003062 neural network model Methods 0.000 description 4
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000208175 Daucus Species 0.000 description 1
- 241001022083 Flemingia Species 0.000 description 1
- 241001290310 Sciaenops Species 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001617 migratory effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012795 verification Methods 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
- 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/044—Recurrent networks, e.g. Hopfield 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/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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Agronomy & Crop Science (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Farming Of Fish And Shellfish (AREA)
Abstract
本发明公开了一种基于原位遥感信息的远洋渔类生境分布预测方法及其应用,方法:处理时间段A内海域B鱼类C的生产统计数据及时间段A内海域B的包括遥感反射比和亮温在内的原位遥感信息得到时间段A内海域B鱼类C的时空同步的CPUE与原位遥感信息的数据;针对任意1~m种原位遥感信息利用BP神经网络建立多个CPUE预测模型,使用所得数据训练并测试建立的各CPUE预测模型得到最终预测模型;将待预测时间段内海域B鱼类C的原位遥感信息输入最终预测模型完成预测。本发明的预测方法,首次提出直接以原位遥感信息作为模型因子考察CPUE,其与传统方法相比,既避免因海洋环境因子提取损失必要的海洋信息,又提高了预测精度。
Description
技术领域
本发明属于远洋渔类生境时空分布的预测技术领域,涉及一种基于原位遥感信息的远洋渔类生境分布预测的方法及其应用。
背景技术
柔鱼是一种大洋性头足类种类,广泛分布于北太平洋海域,主要分布范围为20°-50°N。柔鱼具有季节性洄游特征,包含2个繁殖群体:冬春生群和秋生群体,其中冬春生群体的西部种群是我国鱿钓渔船主要捕捞对象。在产卵季节,柔鱼向南进行繁殖洄游且主要在30°N以南水域进行产卵;而在夏秋季,柔鱼主要向北洄游在亚北极边界附近和过渡区域进行索饵。柔鱼仅有1年生命周期,是短生命周期种类,其资源丰度和空间分布与海洋环境具有紧密关联。为了提高我国柔鱼捕捞效率和产量,深入了解海洋环境和柔鱼资源丰度之间的关系具有重要意义。
国内外有众多学者基于卫星遥感的海洋环境因子产品,开展了柔鱼资源丰度预测模型的研究。这些模型主要考虑的海洋环境因子主要包括叶绿素a (Chl-a)、海表温度(SST)、海表面高度、海表风速等,但是在这些海洋环境信息的获取过程中,需要对卫星遥感手段获取的原始信息进行校正和预处理,这一过程会不可避免地去除部分海洋环境信息,使得基于海洋环境因子产品建立的柔鱼资源丰度预测模型难以更加准确的体现实际的柔鱼资源分布特征,进而导致了对柔鱼资源丰度的预报精度难以得到进一步的提升。
因此,开发一种能够更加准确的体现实际的柔鱼资源分布特征进而提高预测精度的远洋渔类(适用于柔鱼)资源丰度预测方法极具现实意义。
发明内容
本发明的目的在于克服现有的预测方法难以更加准确的体现实际的柔鱼资源分布特征,预报精度较差的缺陷,提供一种能够更加准确的体现实际的柔鱼资源分布特征进而提高预测精度的远洋渔类(适用于柔鱼)资源丰度预测方法,具体是利用较为原始的遥感观测数据完成对远洋渔类生境分布的预测。
为实现上述目的,本发明提供如下技术方案:
基于原位遥感信息的远洋渔类生境分布预测方法,首先,处理时间段A 内海域B鱼类C的生产统计数据及时间段A内海域B的包括遥感反射比(Rrs) 和亮温(BT)(其可用于反演Chl-a和SST)在内的原位遥感信息得到时间段 A内海域B鱼类C的时空同步的CPUE与原位遥感信息的数据,按照1、2、 3……m的次序依次对所有类别的原位遥感信息进行编号,将这些原位遥感信息依次记为x1、x2、x3……xm;然后,针对x1、x2、x3……xm中的任意1~m种原位遥感信息利用基于误差反向传播BP神经网络建立至多2m-1个CPUE预测模型,使用所得数据训练并测试建立的各CPUE预测模型得到最终预测模型;最后,将待预测时间段内海域B鱼类C的原位遥感信息输入最终预测模型即可完成对待预测时间段内海域B鱼类C的CPUE的预测。
本发明的基于原位遥感信息的远洋渔类生境分布预测方法,首次提出直接以原位遥感信息作为模型因子考察CPUE,与使用海洋环境因子的传统方法相比,一方面可以避免因为海洋环境因子提取过程中损失必要的海洋信息,另一方面引入原位遥感信息后预测资源丰度的精度可得到有效提高。本发明通过获取待预测时间段内海域B鱼类C的原位遥感信息即可完成对远洋渔业资源丰度,便捷性好,此外,本发明的CPUE预测模型是基于BP神经网络建立的,BP神经网络无论在网络理论还是在性能方面已比较成熟,其具有很强的非线性映射能力和柔性的网络结构,网络的中间层数、各层的神经元个数可根据具体情况任意设定,适用性好,能够一定程度上降低实际操作及应用时的难度。
作为优选的技术方案:
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,步骤如下:
(1)对时间段A内海域B鱼类C的生产统计数据、时间段A内海域B 的原位遥感信息的数据以及时间段A内海域B的海洋环境因子的数据(原位遥感信息及海洋环境因子的相关数据是通过MODIS和AMSR-E/2卫星数据获取的)进行处理,获取时间段A内海域B鱼类C的时空同步的CPUE、海洋环境因子与原位遥感信息的数据,所述原位遥感信息包括遥感反射比(Rrs) 和亮温(BT),按照1、2、3……m的次序依次对所有类别的原位遥感信息及海洋环境因子进行编号,将这些原位遥感信息及海洋环境因子依次记为x1、 x2、x3……xm;
(2)将步骤(1)获取的时间段A内海域B鱼类C的时空同步的CPUE、海洋环境因子与原位遥感信息的数据分为训练集和测试集;
(3)针对x1、x2、x3……xn中的任意1~n种原位遥感信息及海洋环境因子利用基于误差反向传播BP神经网络建立至多2n-1个CPUE预测模型;
(4)使用训练集训练各CPUE预测模型并使用测试集对各CPUE预测模型进行测试,选取精度最高的模型作为最终预测模型;
(5)将待预测时间段内海域B鱼类C的最终预测模型所需的原位遥感信息输入最终预测模型即可完成对待预测时间段内海域B鱼类C的CPUE的预测。针对不同海域、不同时间的海洋环境的差异性,可以指定分辨率单独建立每个海域单独时间段的预测模型进而进行对应海域鱼类的资源丰度的预测。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,其步骤(4) 中一CPUE预测模型的训练及测试具体步骤为:
(a)从训练集中选取训练样本,判断训练集中训练样本是否已经全部被选取过;
(b)训练集中训练样本未被全部被选取过,则采用步骤(a)选取的训练样本对CPUE预测模型I进行训练,所述CPUE预测模型I的输入为海洋环境因子和原位遥感信息中的任意项,输出为CPUE,反之,则以前次训练得到的CPUE预测模型I作为CPUE预测模型I的优选模型;
(c)从测试集中选取测试样本,将测试样本的对应CPUE预测模型I的输入的相关数据输入步骤(b)训练得到的CPUE预测模型I的优选模型得到该测试样本的测试CPUE;
(c)测试样本的测试CPUE与测试样本的真实CPUE的偏差小于阈值,则CPUE预测模型I的优选模型即为CPUE预测模型I的最终模型,反之,则调整BP神经网络的参数并返回步骤(a)。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,步骤(4) 中所述最终预测模型的获取过程如下:从测试集中选取测试样本,将测试样本的对应各CPUE预测模型的输入的相关数据分别输入各CPUE预测模型的最终模型得到各测试样本的测试CPUE;计算各测试样本的测试CPUE与真实 CPUE的偏差;选取最小偏差值对应的CPUE预测模型的最终模型作为最终预测模型。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,所述海洋环境因子包括海表温度(SST)和叶绿素a浓度(Chl-a)。本发明的海洋环境因子并不仅限于此,此处列举部分常规的海洋环境因子为例而已,包括太平洋年代际震荡指数、海水质量变化和地转流等在内的海洋环境因子也可作为海洋环境因子,针对不同海域不同的鱼类,海洋环境因子的选择并不相同,本领域技术人员根据实际需求调整所选取的海洋环境因子。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,步骤(1) 中的对时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据以及时间段A内海域B的海洋环境因子的数据是指在时间域上按照指定的时间分辨率进行平均处理,在空间域上按照指定的空间分辨率进行重采样。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,所述CPUE 预测模型包括输入层、隐含层和输出层;输入层为海洋环境因子和原位遥感信息中的任意项,输出层为CPUE。
BP神经网络是一种按误差逆传播算法训练的多层前馈网络,是目前应用最广泛的神经网络模型之一。它属于前向网络,可包含多个隐含层,具备处理线性不可分问题的能力。BP神经网络包括信号的前向传播和误差的反向传播两个过程,具备以下特点:
(1)BP神经网络实现了一个从输入到输出的映射功能,一个三层的神经网络能够以任意精度逼近任何非线性连续函数,使得其特别适合于求解内部机制复杂的问题,具有较强的非线性映射能力。
(2)BP神经网络在训练时,能够通过学习自动提取输出、输出数据间的“合理规则”,并自适应的将学习内容记忆于网络的权值中,具有高度自学习和自适应的能力。
(3)BP神经网络在设计模式分类器时,既考虑了网络对所需分类对象的正确分类,还关注了网络在经过训练后对未见过的模式或有噪声污染的模式进行正确的分类,具有将学习成果应用于新知识的能力。
(4)BP神经网络在其局部的或者部分的神经元受到破坏后对全局的训练结果不会造成很大的影响,具有一定的容错能力。
在本专利中,BP神经网络的输入层为上述海洋环境因子或原位遥感信息,经神经网络隐含层传递函数处理后,由输出层输出CPUE;由于输出层输出的CPUE值与实际的CPUE值存在偏差,神经网络通过评估偏差大小,如果偏差过大,则将偏差反向传播至输入层,并通过迭代更新并调整输入层与隐含层之间的连接权值,直到输出层输出的CPUE值与实际的CPUE值之间的偏差小于指定的阈值或达到设置的迭代次数,即可认为BP神经网络的训练学习达到了所需精度要求,训练得到的神经网络模型即可以用于CPUE的预报。
本专利进行BP神经网络的实例构建中,选用的训练函数为 Levenberg-Marquardt(LM)反向传播算法。该算法类似牛顿法,是为了在修正速率是避免计算Hessian矩阵而设计的。当误差性能函数e具有平方和误差的形式时,Hessian矩阵可近似表示为:
H=JTJ (1)
梯度可以表示为
g=JTe (2)
式中J为包含误差性能函数对网络权值一阶导数的雅可比矩阵。LM算法根据下式修正网络权值:
ω(n+1)=ω(n)-[JTJ+μI]-1JTe (3)
当μ=0时,LM算法退化为牛顿法;当μ很大时,上式相当于步长较小的梯度下降法。柔鱼雅可比矩阵的计算比Hessian矩阵易于计算,因此速度非常快。
本专利实施的BP神经网络的构建及实例分析均是在Matlab 2018b的环境下进行并完成的。所设置的BP神经网络的其他主要参数包括:学习速率0.01,训练精度0.001,最大训练次数为10000次,隐含层层数选用程序默认的10 层,隐含层和输出层的传递函数均为S型函数logsig。
如上所述的基于原位遥感信息的远洋渔类生境分布预测方法,步骤(d) 中的调整BP神经网络的参数是指调整输入层、隐含层和输出层之间的连接权值。在符合本领域常识的基础上,上述各优选条件,可任意组合,即得本发明各较佳实例。
本发明还提供基于原位遥感信息的远洋渔类生境分布预测方法的应用,具体为一种电子设备,包括一个或多个处理器、一个或多个存储器、一个或多个程序及数据搜集装置;
所述数据搜集装置用于获取时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据及待预测时间段内海域B鱼类C的原位遥感信息的数据,所述一个或多个程序被存储在所述存储器中,当所述一个或多个程序被所述处理器执行时,使得所述电子设备执行如上所述的基于原位遥感信息的远洋渔类生境分布预测方法。
有益效果:
(1)本发明的基于原位遥感信息的远洋渔类生境分布预测方法,首次提出直接以原位遥感信息作为模型因子考察CPUE,与使用海洋环境因子的传统方法相比,一方面可以避免因为海洋环境因子提取过程中损失必要的海洋信息,另一方面引入原位遥感信息后预测资源丰度的精度可得到有效提高;
(2)本发明的电子设备,结构简单,成本低廉,实现了对远洋渔类生境分布的快速精确预测。
附图说明
图1为北纬44°、东经156°处2004-2018年柔鱼月平均CPUE实测值的空间分布;
图2~7分别为模型I、II、III、IV、V及VI对北纬44°、东经156°处柔鱼月平均CPUE分布的模拟(2004~2013年)和预测(2014~2018年)结果;
图8为2004年9月西北太平洋柔鱼平均CPUE实测值的空间分布;
图9~14分别为模型I、II、III、IV、V及VI对2004年9月西北太平洋柔鱼平均CPUE值的模拟结果;
图15为2017年10月西北太平洋柔鱼平均CPUE实测值的空间分布;
图16~21分别为模型I、II、III、IV、V及VI对2017年10月西北太平洋柔鱼平均CPUE值的预测结果;
图22~27分别为模型I、II、III、IV、V及VI模拟2004-2013年西北太平洋柔鱼月平均CPUE的RMSE分布;
图28~33分别为模型I、II、III、IV、V及VI预测2014-2018年西北太平洋柔鱼月平均CPUE的RMSE分布。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
一种基于原位遥感信息的远洋渔类生境分布预测方法,其步骤如下:
(1)对2004~2018年5~12月,坐标范围为35°~50°N、150°~175°E 内柔鱼(即对应前文中时间段A内海域B鱼类C)的生产统计数据、2004~2018 年5~12月,坐标范围为35°~50°N、150°~175°E内的原位遥感信息(遥感反射比Rrs和亮温BT)的数据以及2004~2018年5~12月,坐标范围为35°~50°N、150°~175°E内的海洋环境因子(海表温度SST和叶绿素a浓度 Chl-a)的数据进行处理(即在时间域上按照指定的时间分辨率(月)进行平均处理,在空间域上按照1°×1°的空间分辨率进行重采样),获取2004~2018 年5~12月,坐标范围为35°~50°N、150°~175°E内柔鱼的时空同步的 CPUE、海洋环境因子与原位遥感信息的数据,按照1、2、3……m的次序依次对所有类别的原位遥感信息及海洋环境因子进行编号,将这些原位遥感信息及海洋环境因子依次记为x1、x2、x3……xm。
数据的获取及处理具体如下:
本专利中的柔鱼生产统计数据来自上海海洋大学中国远洋渔业数据中心。选取2004~2018年的西北太平洋海域的柔鱼历史渔业捕捞数据,统计获取逐月平均、空间分辨率为1°×1°的柔鱼CPUE值。对不同卫星平台的海洋环境因子的数据来源及处理方法如下:
本专利所使用的Chl-a和Rrs数据均来自于中分辨率成像光谱仪MODIS 传感器(http://oceandata.sci.gsfc.nasa.gov/)。MODIS是由美国宇航局研制的大型空间传感器,搭载于Terra和Aqua上。MODIS数据有36个波段,每日或每两日可获得一次全球观测数据。为了保证研究区域内数据的完整性,我们选取了Terra和Aqua卫星平台上的月平均SST和Rrs数据产品,空间分辨率为9km。利用MODIS生成Chl-a产品时,主要利用了465、555和645nm的 Rrs数据,实际操作中,我们也仅分析了这三个波段的Rrs数据。将Terra和 Aqua得到的数据产品进行平均后,然后对平均后的SST和Rrs进行重采样处理,获得了2004-2018年5~12月的空间分辨率为1°×1°的月平均Chl-a和Rrs。
SST和BT数据来自于高级微波成像扫描辐射计(AMSR-E)及其后继卫星AMSR-2(https://sharaku.eorc.jaxa.jp/AMSR/index.html)。AMSR-E是搭载在Aqua卫星上的被动微波辐射计,有6.9、10.7、18.7,23.8、36.5和89.0GHz 共6个频率、总计12个通道的数据,空间分辨率为25km,数据有效时间段为2002年6月至2011年9月。AMSR2是继AMSR-E之后又一新型的传感器,它搭载在日本的GCOMW1(Global Change Observation Mission 1st-Water)“SHIZUKU”卫星上,该传感器包括6.9、7.3、10.65、18.7、23.8、36.5GHz 和89GHz共7个频率,每个频率有水平(H)和垂直(V)两种极化方式,空间分辨率为10km,数据有效时间段为2012年7月至今。本发明直接使用了月平均的SST和BT数据,然后对SST和BT在空间域上进行了重采样处理,获得2004-2018年5~12月的空间分辨率为1°×1°的月平均SST和BT。需要进一步说明的是,本专利仅选择处理和分析了6.9H、6.9V、10.7H、10.7V、 23.8V、36.5H和36.5V四个频率共7个通道的BT数据。
(2)将步骤(1)获取的2004~2018年5~12月,坐标范围为35°~50° N、150°~175°E内柔鱼的时空同步的CPUE、海洋环境因子与原位遥感信息的数据中分别选用2004-2013年5~12月和2014-2017年5~12月的相关数据,分别记为T1和T2,而后分别将T1和T2分为训练集和测试集;
(3)针对x1、x2、x3……xn中的任意1~n种原位遥感信息及海洋环境因子利用BP神经网络建立6个CPUE预测模型,这六个CPUE预测模型的输入层不同,其输出层均为CPUE,CPUE预测模型包括输入层、隐含层和输出层, CPUE预测模型编号与输入层具体对应关系如下表所示:
(4)使用训练集训练各CPUE预测模型并使用测试集对各CPUE预测模型进行测试,选取精度最高的模型作为最终预测模型;
其中一CPUE预测模型的训练及测试具体步骤为:
(a)从训练集中选取训练样本,判断训练集中训练样本是否已经全部被选取过;
(b)训练集中训练样本未被全部被选取过,则采用步骤(a)选取的训练样本对CPUE预测模型I进行训练,所述CPUE预测模型I的输入为海洋环境因子和原位遥感信息中的任意项,输出为CPUE,反之,则以前次训练得到的CPUE预测模型I作为CPUE预测模型I的优选模型;
(c)从测试集中选取测试样本,将测试样本的对应CPUE预测模型I的输入的相关数据输入步骤(b)训练得到的CPUE预测模型I的优选模型得到该测试样本的测试CPUE;
(c)测试样本的测试CPUE与测试样本的真实CPUE的偏差小于阈值,则CPUE预测模型I的优选模型即为CPUE预测模型I的最终模型,反之,则调整BP神经网络的参数(调整输入层、隐含层和输出层之间的连接权值)并返回步骤(a)。
最终预测模型的获取过程具体如下:从测试集中选取测试样本,将测试样本的对应各CPUE预测模型的输入的相关数据分别输入各CPUE预测模型的最终模型得到各测试样本的测试CPUE;计算各测试样本的测试CPUE与真实CPUE的偏差;选取最小偏差值对应的CPUE预测模型的最终模型作为最终预测模型。
本专利实施的BP神经网络的构建及实例分析均是在Matlab 2018b的环境下进行并完成的。所设置的BP神经网络的其他主要参数包括:学习速率0.01,训练精度0.001,最大训练次数为10000次,隐含层层数选用程序默认的10 层,隐含层和输出层的传递函数均为S型函数logsig。
对以上模型I~模型VI进行优选:
图1是北纬44°、东经156°处2004~2018年柔鱼月平均CPUE实测值的空间分布,图2~7分别给出了模型I~VI对北纬44°、东经156°处的柔鱼月平均CPUE进行模拟(2004~2013年)和预测(2014~2018年)的结果。将图2~7的结果与图1的实际结果进行比较分析后发现,基于T1时间段的Chl-a 对北纬44°、东经156°处模拟的柔鱼月平均CPUE的标准差(STD)和均方根误差(RMSE)分别为0.75和0.78t/d,而基于Chl-a的原位遥感信息Rrs 进行该处柔鱼月平均CPUE模拟的RMSE和STD分别为0.54和0.66t/d,分别减少了27%和15%。基于T2时间段的Chl-a对北纬44°、东经156°处的柔鱼月平均CPUE预测的RMSE和STD分别为2.23和2.27t/d,对应的基于原位遥感信息Rrs的预测的RMSE和STD分别降低到2.11和2.21t/d,分别减少了6%和2%。
基于T1时间段的SST对北纬44°、东经156°处的柔鱼月平均CPUE 进行模拟的RMSE和STD分别为0.68和0.76t/d,预测的RMSE和STD分别为1.88和2.0t/d。而利用其原位遥感信息BT进行模拟的RMSE和STD分别为0.54和0.68t/d,预测的RMSE和STD分别为1.35和1.32t/d。因此,相比于SST,其原位遥感信息BT在该网格处模拟的柔鱼CPUE的RMSE和STD 分别降低了21%和11%,预测结果的RMSE和STD分别降低了28%和34%。这说明原位遥感信息对柔鱼CPUE的模拟和预测的精度比海洋环境因子产品更高。
为了验证多种海洋环境因子或原位遥感信息组合时对柔鱼CPUE的模拟和预测情况,我们还比较分析了Chl-a和SST联合(模型V),以及Rrs和BT 联合(模型VI)对柔鱼CPUE进行模拟和预测的结果。结果发现,多种海洋环境因子组合和多种原位遥感信息组合模拟北纬44°、东经156°处的柔鱼 CPUE的结果较为一致(RMSE=0.55t/d;STD=0.68t/d),但多种原位遥感信息组合比多种海洋环境因子组合预测的柔鱼CPUE的RMSE和STD提高了 29%和32%。
图8为2004年9月西北太平洋柔鱼平均CPUE实测值的空间分布,图9~14 分别给出了模型I~VI模拟的柔鱼月平均CPUE在2004年9月的分布情况。结果显示,基于Rrs建模模拟的柔鱼CPUE的RMSE和STD比基于Chl-a建模时在2004年9月改善了54%和51%,基于BT建模模拟的柔鱼CPUE的 RMSE和STD比基于SST建模时在2004年9月改善了29%和28%,而基于Rrs和BT建模模拟的柔鱼CPUE的RMSE和STD也比基于Chl-a和SST建模时在2004年9月改善了9%和6%。
为了进一步验证预测结果的空间分布的精度,图15为2017年10月西北太平洋柔鱼平均CPUE实测值的空间分布,图9~14分别给出了模型I~VI预测的柔鱼月平均CPUE在2017年10月的分布情况。结果发现,基于Rrs建模预测的柔鱼CPUE的RMSE和STD比基于Chl-a建模时在2004年9月改善了65%和58%,基于BT建模预测的柔鱼CPUE的RMSE和STD比基于SST 建模时在2004年9月改善了18%和30%,而基于Rrs和BT建模预测的柔鱼 CPUE的RMSE和STD也比基于Chl-a和SST建模时在2004年9月改善了 5%和8%。
在每个1°×1°的网格点上,也可以利用每个月的CPUE模拟和预测结果,分析不同方案建立的BP神经网络模型模拟和预测的柔鱼CPUE的精度分布情况。图22~27给出了模型I~VI模拟的柔鱼CPUE的RMSE分布情况。结果发现,采用海洋环境因子及其组合,以及采用原位遥感信息及其组合均能以较高的精度(平均RMSE<0.5t/d)模拟的西北太平洋地区的柔鱼CPUE。
图28~33还给出了模型I~VI预测的柔鱼CPUE的RMSE的空间分布。结果发现,Rrs预测的柔鱼CPUE比Chl-a预测结果的平均RMSE减少了4%, BT预测的柔鱼CPUE比SST预测结果的平均RMSE减少了9%,而Rrs和 BT预测的柔鱼CPUE比Chl-a和SST预测结果的平均RMSE减少了14%。即证实了将原位遥感信息引入BP神经网络模型,比使用传统的海洋环境因子在预测西北太平洋柔鱼CPUE时的精度得到了有效提高。
综上所述,最终预测模型为模型VI。
(5)将待预测时间段内坐标范围为35°~50°N、150°~175°E内柔鱼的最终预测模型(模型VI)所需的原位遥感信息(Rrs和BT)输入最终预测模型(模型VI)即可完成对待预测时间段内坐标范围为35°~50°N、150°~175°E内柔鱼的CPUE的预测。
经验证,本发明的基于原位遥感信息的远洋渔类生境分布预测方法,首次提出直接以原位遥感信息作为模型因子考察CPUE,与使用海洋环境因子的传统方法相比,一方面可以避免因为海洋环境因子提取过程中损失必要的海洋信息,另一方面引入原位遥感信息后预测资源丰度的精度可得到有效提高,极具应用前景。
实施例2
一种电子设备,包括一个或多个处理器、一个或多个存储器、一个或多个程序及数据搜集装置;
数据搜集装置用于获取时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据及待预测时间段内海域B鱼类C的原位遥感信息的数据,一个或多个程序被存储在存储器中,当一个或多个程序被处理器执行时,使得电子设备执行如实施例1所述的基于原位遥感信息的远洋渔类生境分布预测方法。
虽然以上描述了本发明的具体实施方式,但是本领域的技术人员应该理解,这些仅是举例说明,在不违背本发明的原理和实质的前提下,可以对这些实施方式做出多种变更或修改。
Claims (9)
1.基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,首先,处理时间段A内海域B鱼类C的生产统计数据及时间段A内海域B的包括遥感反射比和亮温在内的原位遥感信息得到时间段A内海域B鱼类C的时空同步的CPUE与原位遥感信息的数据,按照1、2、3……m的次序依次对所有类别的原位遥感信息进行编号,将这些原位遥感信息依次记为x1、x2、x3……xm;然后,针对x1、x2、x3……xm中的任意1~m种原位遥感信息利用BP神经网络建立至多2m-1个CPUE预测模型,使用所得数据训练并测试建立的各CPUE预测模型得到最终预测模型;最后,将待预测时间段内海域B鱼类C的原位遥感信息输入最终预测模型即可完成对待预测时间段内海域B鱼类C的CPUE的预测。
2.根据权利要求1所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,步骤如下:
(1)对时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据以及时间段A内海域B的海洋环境因子的数据进行处理,获取时间段A内海域B鱼类C的时空同步的CPUE、海洋环境因子与原位遥感信息的数据,所述原位遥感信息包括遥感反射比和亮温,按照1、2、3……m的次序依次对所有类别的原位遥感信息及海洋环境因子进行编号,将这些原位遥感信息及海洋环境因子依次记为x1、x2、x3……xm;
(2)将步骤(1)获取的时间段A内海域B鱼类C的时空同步的CPUE、海洋环境因子与原位遥感信息的数据分为训练集和测试集;
(3)针对x1、x2、x3……xn中的任意1~n种原位遥感信息及海洋环境因子利用BP神经网络建立至多2n-1个CPUE预测模型;
(4)使用训练集训练各CPUE预测模型并使用测试集对各CPUE预测模型进行测试,选取精度最高的模型作为最终预测模型;
(5)将待预测时间段内海域B鱼类C的最终预测模型所需的原位遥感信息输入最终预测模型即可完成对待预测时间段内海域B鱼类C的CPUE的预测。
3.根据权利要求2所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,步骤(4)中一CPUE预测模型的训练及测试具体步骤为:
(a)从训练集中选取训练样本,判断训练集中训练样本是否已经全部被选取过;
(b)训练集中训练样本未被全部被选取过,则采用步骤(a)选取的训练样本对CPUE预测模型I进行训练,所述CPUE预测模型I的输入为海洋环境因子和原位遥感信息中的任意项,输出为CPUE,反之,则以前次训练得到的CPUE预测模型I作为CPUE预测模型I的优选模型;
(c)从测试集中选取测试样本,将测试样本的对应CPUE预测模型I的输入的相关数据输入步骤(b)训练得到的CPUE预测模型I的优选模型得到该测试样本的测试CPUE;
(c)测试样本的测试CPUE与测试样本的真实CPUE的偏差小于阈值,则CPUE预测模型I的优选模型即为CPUE预测模型I的最终模型,反之,则调整BP神经网络的参数并返回步骤(a)。
4.根据权利要求3所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,步骤(4)中所述最终预测模型的获取过程如下:从测试集中选取测试样本,将测试样本的对应各CPUE预测模型的输入的相关数据分别输入各CPUE预测模型的最终模型得到各测试样本的测试CPUE;计算各测试样本的测试CPUE与真实CPUE的偏差;选取最小偏差值对应的CPUE预测模型的最终模型作为最终预测模型。
5.根据权利要求2所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,所述海洋环境因子包括海表温度和叶绿素a浓度。
6.根据权利要求2所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,步骤(1)中的对时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据以及时间段A内海域B的海洋环境因子的数据是指在时间域上按照指定的时间分辨率进行平均处理,在空间域上按照指定的空间分辨率进行重采样。
7.根据权利要求3所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,所述CPUE预测模型包括输入层、隐含层和输出层;输入层为海洋环境因子和原位遥感信息中的任意项,输出层为CPUE。
8.根据权利要求7所述的基于原位遥感信息的远洋渔类生境分布预测方法,其特征在于,步骤(d)中的调整BP神经网络的参数是指调整输入层、隐含层和输出层之间的连接权值。
9.一种电子设备,其特征在于,包括一个或多个处理器、一个或多个存储器、一个或多个程序及数据搜集装置;
所述数据搜集装置用于获取时间段A内海域B鱼类C的生产统计数据、时间段A内海域B的原位遥感信息的数据及待预测时间段内海域B鱼类C的原位遥感信息的数据,所述一个或多个程序被存储在所述存储器中,当所述一个或多个程序被所述处理器执行时,使得所述电子设备执行如权利要求1~8任一项所述的基于原位遥感信息的远洋渔类生境分布预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110359602.6A CN113159400A (zh) | 2021-04-02 | 2021-04-02 | 基于原位遥感信息的远洋渔类生境分布预测方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110359602.6A CN113159400A (zh) | 2021-04-02 | 2021-04-02 | 基于原位遥感信息的远洋渔类生境分布预测方法及其应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113159400A true CN113159400A (zh) | 2021-07-23 |
Family
ID=76886284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110359602.6A Pending CN113159400A (zh) | 2021-04-02 | 2021-04-02 | 基于原位遥感信息的远洋渔类生境分布预测方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113159400A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630238A (zh) * | 2013-11-22 | 2014-03-12 | 中国科学院南京地理与湖泊研究所 | 一种浅水湖泊水面成像光谱采集系统及同步自动监测方法 |
CN106157162A (zh) * | 2016-07-22 | 2016-11-23 | 上海海洋大学 | 一种北太平洋鱿鱼中心渔场预测方法 |
CN108665104A (zh) * | 2018-05-14 | 2018-10-16 | 华际卫星通信有限公司 | 一种基于lstm的渔场预报方法 |
CN111652443A (zh) * | 2020-06-05 | 2020-09-11 | 上海海洋大学 | 综合多源卫星遥感预测远洋渔业资源丰度的方法及其应用 |
-
2021
- 2021-04-02 CN CN202110359602.6A patent/CN113159400A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630238A (zh) * | 2013-11-22 | 2014-03-12 | 中国科学院南京地理与湖泊研究所 | 一种浅水湖泊水面成像光谱采集系统及同步自动监测方法 |
CN106157162A (zh) * | 2016-07-22 | 2016-11-23 | 上海海洋大学 | 一种北太平洋鱿鱼中心渔场预测方法 |
CN108665104A (zh) * | 2018-05-14 | 2018-10-16 | 华际卫星通信有限公司 | 一种基于lstm的渔场预报方法 |
CN111652443A (zh) * | 2020-06-05 | 2020-09-11 | 上海海洋大学 | 综合多源卫星遥感预测远洋渔业资源丰度的方法及其应用 |
Non-Patent Citations (3)
Title |
---|
唐浩;许柳雄;陈新军;朱国平;周成;王学昉;: "基于GAM模型研究时空及环境因子对中西太平洋鲣鱼渔场的影响", 海洋环境科学, no. 04, 20 August 2013 (2013-08-20) * |
徐升;张鹰;: "长江口水域多光谱遥感水深反演模型研究", 地理与地理信息科学, no. 03, 30 May 2006 (2006-05-30) * |
阎雨, 陈圣波, 田静, 谭笑平: "卫星遥感估产技术的发展与展望", 吉林农业大学学报, no. 02, 30 April 2004 (2004-04-30) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lehodey et al. | Modelling climate‐related variability of tuna populations from a coupled ocean–biogeochemical‐populations dynamics model | |
JP2019511783A (ja) | 灰色システムに基づく遠洋イカ類資源の存在度の予測方法 | |
Gregg et al. | Skill assessment in ocean biological data assimilation | |
Liu et al. | San Francisco Bay nutrients and plankton dynamics as simulated by a coupled hydrodynamic-ecosystem model | |
Kearney et al. | Coupling planktonic ecosystem and fisheries food web models for a pelagic ecosystem: Description and validation for the subarctic Pacific | |
Bricker et al. | Integration of ecosystem-based models into an existing interactive web-based tool for improved aquaculture decision-making | |
Johnson et al. | agTrend: A Bayesian approach for estimating trends of aggregated abundance | |
Filgueira et al. | Ecosystem modelling for ecosystem-based management of bivalve aquaculture sites in data‑poor environments | |
JP7022456B2 (ja) | 掃過面積に基づいて推定された漁業資源量の補正方法及びその応用 | |
Melià et al. | Density and temperature-dependence of vital rates in the Manila clam Tapes philippinarum: a stochastic demographic model | |
Luque et al. | Direct ageing of Thunnus thynnus from the eastern Atlantic Ocean and western Mediterranean Sea using dorsal fin spines | |
Swannack et al. | A robust, spatially explicit model for identifying oyster restoration sites: case studies on the Atlantic and Gulf coasts | |
CN111784034B (zh) | 一种影响智利海域美洲赤鱿渔场的关键环境因子的筛选及探测方法 | |
Jansen | First‐year survival of North East Atlantic mackerel (Scomber scombrus) from 1998 to 2012 appears to be driven by availability of Calanus, a preferred copepod prey | |
Leahy et al. | River flows affect the growth of a tropical finfish in the wet-dry rivers of northern Australia, with implications for water resource development | |
CN107169679B (zh) | 一种利用食物网稳定性评价湿地生态补水实施效果的方法 | |
CN115758074A (zh) | 基于多源数据的高空间分辨率海水二氧化碳分压重构方法 | |
Cwynar et al. | Depth distribution of chironomids and an evaluation of site-specific and regional lake-depth inference models: a good model gone bad? | |
CN113159400A (zh) | 基于原位遥感信息的远洋渔类生境分布预测方法及其应用 | |
Kruse et al. | Exploratory simulation of English sole recruitment mechanisms | |
Giacoletti et al. | Functional trait-based layers-an aquaculture siting tool for the Mediterranean Sea | |
Duarte et al. | Space‐for‐time is not necessarily a substitution when monitoring the distribution of pelagic fishes in the San Francisco Bay‐Delta | |
Richard et al. | Relating annual increments of the endangered Blanding's turtle plastron growth to climate | |
Cerco | A multi-module approach to calculation of oyster (Crassostrea virginica) environmental benefits | |
CN114330917A (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 |