CN106376511A - 一种基于被动声学信息的淡水鱼存活率预测装置及方法 - Google Patents
一种基于被动声学信息的淡水鱼存活率预测装置及方法 Download PDFInfo
- Publication number
- CN106376511A CN106376511A CN201610801989.5A CN201610801989A CN106376511A CN 106376511 A CN106376511 A CN 106376511A CN 201610801989 A CN201610801989 A CN 201610801989A CN 106376511 A CN106376511 A CN 106376511A
- Authority
- CN
- China
- Prior art keywords
- fish
- survival rate
- band energy
- energy
- sample
- 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 220
- 238000000034 method Methods 0.000 title claims abstract description 101
- 230000004083 survival effect Effects 0.000 title claims abstract description 64
- 239000013505 freshwater Substances 0.000 title claims abstract description 42
- 230000005236 sound signal Effects 0.000 claims description 56
- 241000252229 Carassius auratus Species 0.000 claims description 37
- 238000000354 decomposition reaction Methods 0.000 claims description 35
- 241001519451 Abramis brama Species 0.000 claims description 34
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 22
- 239000000284 extract Substances 0.000 claims description 17
- 230000009467 reduction Effects 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 13
- 238000012417 linear regression Methods 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 230000006978 adaptation Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000011835 investigation Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 10
- 238000001514 detection method Methods 0.000 abstract description 9
- 230000008569 process Effects 0.000 abstract description 3
- 238000009395 breeding Methods 0.000 abstract 2
- 230000001488 breeding effect Effects 0.000 abstract 2
- 229920000742 Cotton Polymers 0.000 abstract 1
- 230000001413 cellular effect Effects 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000252335 Acipenser Species 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 210000004712 air sac Anatomy 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012258 culturing Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000013305 food Nutrition 0.000 description 2
- 230000009931 harmful effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 230000037361 pathway Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 241001609213 Carassius carassius Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000796 flavoring agent Substances 0.000 description 1
- 235000019634 flavors Nutrition 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 235000013372 meat Nutrition 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 235000008935 nutritious Nutrition 0.000 description 1
- 230000000050 nutritive effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K63/00—Receptacles for live fish, e.g. aquaria; Terraria
- A01K63/003—Aquaria; Terraria
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K61/00—Culture of aquatic animals
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K63/00—Receptacles for live fish, e.g. aquaria; Terraria
- A01K63/003—Aquaria; Terraria
- A01K63/006—Accessories for aquaria or terraria
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/80—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in fisheries management
- Y02A40/81—Aquaculture, e.g. of fish
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Environmental Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Biodiversity & Conservation Biology (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Farming Of Fish And Shellfish (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明属于农产品品质分析技术领域,具体涉及一种基于被动声学信息的淡水鱼存活率预测装置及方法。提取原始鱼声信号的装置包括鱼箱,所述鱼箱内各设有水听器,所述水听器分别与声学记录仪连接,所述鱼箱外设有隔音棉。淡水鱼存活率预测方法,包括采集原始鱼声信号,对鱼声信号进行消噪处理,提取鱼声信号的特征参数:短时平均能量、短时平均过零率、频段能量,构建特征向量,划分样本集,挑选鱼声信号特征频段,建立淡水鱼存活率预测模型。本发明建立的淡水鱼存活率预测模型能实现淡水鱼存活率的在线检测,可用于淡水鱼养殖与活鱼运输过程中鱼的存活率检测,对提高淡水鱼养殖和活鱼运输过程中鱼的存活率具有重要意义。
Description
技术领域
本发明属于农产品品质分析技术领域,具体涉及一种基于被动声学信息的淡水鱼存活率预测装置及方法。
背景技术
我国是传统的渔业大国,水产品总产量和出口量均居世界首位,其中淡水鱼的养殖产量占内陆水产品总养殖产量的较大比重,是我国主要的水产养殖种类。根据中国渔业统计年鉴公布的数据可知,2009~2013年之间,我国的淡水鱼年产量平均增长率为5.5%,到2013年,我国淡水鱼年产量就已达到2635.08万吨,且会持续增长(梁焕秋等,2014;刘佳,2014)。鱼类肉味鲜美,营养丰富,含有人体所必需的动物蛋白、钙、多种维生素等营养物质,是日常生活中不可或缺的食材。而冷冻鱼和鲜活鱼无论是在营养价值还是在销售价格上都有较大的差别,普通消费者更多的倾向于购买鲜活鱼类。目前,国内和国际的淡水鱼市场需求一直在持续增长,其中90%以上淡水鱼以鲜活形式在市面上销售。
相关研究表明,鱼类在水下有着各种各样的发声方式,如骨骼发出的声音、鳃盖的闭合声、游动时激起的水花和旋窝、鱼在沉浮时鱼鳔受到撞击震动以及鱼鳔换气时发出的声音、鱼体与鱼体之间发生碰撞所发出的声音等,声音信号丰富。而之所以我们听不到鱼所发出的任何声音,是因为水的密度是相当于空气密度的7500倍,鱼在水里面发出的声音几乎都在水的介质中逐渐消失。相关研究表明,鱼类的声音是用来实现种间或者种内的信息传递,包括生殖时的集群声、躲避敌害发生的声音、寻找食物的试探声以及识别同类的呼叫声等,其特征与鱼的种类、生理状态等生命活动紧密关联,有着特定的生物学意义。国内外对声学信息在渔业中的应用研究主要集中在鱼类发声特性及发声机理、鱼群发声信号处理、鱼类个体的水声目标特性和基于水声信号的鱼类资源评估方法等方面。在针对鱼类声音信号的研究中,大都围绕在通过分析水下的鱼声信号判断鱼群的地理分布等;个别针对鱼声信号本身的研究也以海鱼和船舶辐射噪声为主,提取信号的频域特征,目的是进行海鱼种类识别以及对鱼声信号与船舶声音信号进行识别;近年来,鱼探仪已成为渔业资源调查与评估的一个重要工具,广泛地应用于监测鱼群行为、鉴别鱼类性别、评估其它水生生物量、探测水体底部类型、监测水质和水生生态系统。
鱼声检测方式主要包括主动声呐检测和被动声呐检测,现有研究中采用主动声呐方式较多,将被动水声信息检测技术应用于淡水鱼存活率检测的应用目前尚未发现。被动声学探测具有成本适中、灵敏度高、发展较成熟的优点,因此,基于被动声学信息的淡水鱼状态检测,进而实现对淡水鱼存活率预测更具有应用潜力,也是该领域的发展趋势。
发明内容
本发明的目的采用语音分析技术,提取了鱼声信号特征值,采用多元线性回归方法建立了淡水鱼存活率的预测模型,实现了淡水鱼存活率的快速检测。
本发明技术方案:
一种基于被动声学信息的淡水鱼存活率预测装置,提取原始鱼声信号的装置包括第一鱼箱,第二鱼箱,所述第一鱼箱内设有第一水听器,所述第二鱼箱内设有第二水听器,所述第一水听器、第二水听器分别与声学记录仪连接。
优选地,所述第一鱼箱、第一水听器、第二水听器设有若干个。
优选地,所述第一水听器、第二水听器分别位于第一鱼箱、第二鱼箱水面以下,所述第一鱼箱、第二鱼箱外设有隔音棉。
优选地,所述隔音棉为橡胶海绵。
一种基于声呐声学信息的淡水鱼存活率预测方法,所述方法具体包括以下步骤:
1)消噪处理:提取原始鱼声信号,并对所得的鱼声信号进行消噪处理;
2)提取特征参数:根据所述步骤1)消噪处理后的鱼声信号,提取鱼声短时平均能量、鱼声短时平均过零率,再对消噪处理后的鱼声信号进行分解,采用四层、五层、六层和七层小波包分解方法对鱼声信号进行了频段划分,提取各频段能量;
3)构建特征向量:根据所述步骤2)提取的短时平均能量、短时平均过零率和各频段能量构建特征向量;
4)样本集划分:采集不同条数不同时间的鱼声信号,构成一个鱼声样本集,将样本集划分为训练集和验证集;
5)特征频段挑选、特征向量降维:对经过步骤4)划分样本集的鱼声音信号进行Z-score标准化预处理,并采用竞争自适应重加权采样法、多元线性回归对特征频段再次进行优选,剔除不显著的特征频段,得到鱼声信号特征频段能量,对步骤3)特征向量降维,得到降维后的鱼声信号特征向量;
6)模型建立:对所述步骤5)降维后的鱼声信号特征向量,采用多元线性回归法和偏最小二乘法分别建立鱼存活率预测模型;
7)预测存活率:对未知淡水鱼样本鱼声信号进行检测,并将鱼声样本的特征向量经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数;
完成鱼存活率的预测。
优选地,所述步骤4)采用Rank-SPXY法进行样本集的划分,该方法由两部分组成,首先是Rank部分,即将样本按因变量活鱼条数的升序排序,然后将样本等分为m份;其次是SPXY法部分,即在等分m份的每个区间内采用SPXY法选取出训练集,其余的样本自动归为验证集,m分别取值为5和10。
进一步优选地,所述步骤6)鲫鱼存活率预测模型,选取短时平均能量、短时平均过零率以及基于6层小波包分解的频段能量为鲫鱼声音信号样本的最优特征提取方法,Rank-SPXY法样本集划分方法中m=10,建立鲫鱼存活率MLR预测模型,回归方程如下:
y=40.471-3.095x2+1.710x4-1.981x5-1.768x6+2.349x7-10.883x11+7.301x12-1.306x16-2.187x18+19.417x25+9.734x28+70.133x35-25.264x43-79.860x47-40.098x50+26.155x55-24.005x61+31.320x62,
y为存活率,xi为鲫鱼声音信号经Z-score标准化预处理后特征值,鱼声短时平均能量x1、鱼声短时平均过零率x2,xi中i≥3时为基于6层小波包分解的频段能量,8~16Hz频段能量x4,16~24Hz频段能量x5,24~32Hz频段能量x6,32~40Hz频段能量x7,64~72Hz频段能量x11,72~80Hz频段能量x12,104~112Hz频段能量x16,120~128Hz频段能量x18,176~184Hz频段能量x25,200~208Hz频段能量x28,256~264Hz频段能量x35,320~328Hz频段能量x43,352~360Hz频段能量x47,376~384Hz频段能量x50,416~424Hz频段能量x55,464~472Hz频段能量x61,472~480Hz频段能量x62。
更进一步优选地,鲫鱼声音信号的特征频段为256~264Hz、352~360Hz、376~384Hz。
进一步优选地,所述步骤6)鳊鱼存活率预测模型,选取短时平均能量、短时平均过零率以及基于6层小波包分解的频段能量为鳊鱼声音信号样本的最优特征提取方法,建立鳊鱼存活率MLR预测模型,Rank-SPXY法样本集划分方法中m=5,回归方程如下:
y=4.384+1.415x5+2.681x9+8.356x14+2.694x18-3.290x24,
y为存活率,xi为鳊鱼声音信号经Z-score标准化预处理后特征值,鱼声短时平均能量x1、鱼声短时平均过零率x2,xi中i≥3时为基于6层小波包分解的频段能量,16~24Hz频段能量x5,48~56Hz频段能量x9,88~96Hz频段能量x14,120~128Hz频段能量x18,168~178Hz频段能量x24。
更进一步优选地,鳊鱼声音信号的特征频段为48~56Hz、88~96Hz、120~128Hz、168~178Hz。
所述方法在淡水鱼养殖、活鱼运输以及渔业资料调查中的应用。
本发明提供的一种基于被动声学信息的淡水鱼存活率预测装置和方法,有益效果如下:
1、本发明确定了鱼声信号的特征参数和特征频段,建立的淡水鱼存活率预测模型的方法,实现了活鱼的存活状况的快速、无损检测。
2、应用本发明建立的淡水鱼存活率模型,可以实现淡水鱼养殖和活鱼运输过程中存活率的在线检测和预警,即当存活率降低的时候,发出预警。本发明有利于准确判断养殖及运输过程中鱼的存活状况,即时采取相应的对策,避免对鱼类产生不良影响,提高鱼的存活率,从而降低养殖或运输成本。
附图说明
图1:本发明鱼声信号采集系统结构图;
图2:本发明的六层小波包分解结构方法;
其中1为第一鱼箱,2为第二鱼箱,3声学记录仪,4为第一水听器,5为第二水听器,6为隔音棉。
具体实施方式
基于小波包分解的频段能量提取
用小波包对淡水鱼声音信号进行分解时要选择适合的分解尺度和小波包基函数,分解尺度的确定与鱼声信号的主要频段以及采样频率有关。用不同的小波包基函数对同一个鱼声信号进行分解将会得到不同的结果,所以在选用小波包对信号进行分解时,要根据不同信号的特征以及反复对比分析来选择一个最好的小波包基函数,选择的小波包基函数需要满足以下几点要求:在时域和频域都具有一定的局部化分析能力;在时域具有紧支撑性,在频域具有快速衰减性;至少具有一阶消失矩;具有良好的分解与重构性。满足上述要求的常用小波有SymletsA(symN)小波、Coiflet(coifN)小波、Daubechies(dbN)小波等。鱼声信号分析处理实际应用中一般选择dbN小波,它是Daubechies从双尺度方程系数{hk}中定义出来的离散正交小波,是离散小波变换的良好工具。N指小波的阶数。
本发明采集鱼声信号时的采样频率是4000Hz,而常见淡水鱼声音信号的主要频率成分是500Hz以下的低频部分,因此,本发明只针对0-500Hz以内的鱼声信号进行分析研究。采用小波包分解对鱼声信号进行特征参数提取,提取信号各频段内能量作为分类识别的特征参数。其步骤如下:
第一步:选取合适的分解尺度和小波基函数,对鱼声信号S进行分解。在对海洋鱼类声音信号进行小波包分解时,分解尺度一般选取为三层或四层,而由于常见淡水鱼声音信号比较微弱,且频率较低,因此,本发明分别选取四层、五层、六层作为分解尺度,选取db1小波作为小波包基函数,以六层小波包分解为例,其分解结构如图1所示。
图2中,每个节点都代表了一定的特征,比如,节点(0,0)代表原始鱼声信号S;节点(1,0)代表第一层小波包分解第0个节点的系数;节点(1,1)代表第一层小波包分解第1个节点的系数;依次类推。
第二步:利用分解的节点系数对分解尺度上的信号进行单尺度重构,得到各频段内的小波包重构信号。以Sij表示节点(i,j)的小波包重构信号,则原始鱼声信号S可以表示为:
S=S60+S61+S62+S63+…+S662+S663
四层、五层、和七层小波包分解方法与此方法相同。
实施例1
一种基于被动声学信息的淡水鱼存活率预测装置,提取原始鱼声信号的装置包括第一鱼箱1,第二鱼箱2,所述第一鱼箱1内设有第一水听器4,所述第二鱼箱2内设有第二水听器5,所述第一水听器4、第二水听器5分别与声学记录仪3连接。第一鱼箱用来采集淡水鱼声音信号,第二鱼箱用来采集背景环境噪音,用作后期数据处理时的背景环境噪音消噪。
所述第一鱼箱1、第一水听器4、第二水听器5设有若干个。当用于运输或养殖时,放置鱼的第一鱼箱有多个,因此每个放置鱼的鱼箱都需要采集鱼声信号。
在大型养殖或运输过程中,鱼箱很大,放置一个水听器不能采集到完整的鱼声信息,因此需要设置若干个。
所述第一水听器4、第二水听器5分别位于第一鱼箱1、第二鱼箱2水面以下,所述第一鱼箱1、第二鱼箱2外设有隔音棉6。为了保证鱼声信号采集的完整性,减少外界的干扰,因此水听器设置于水面以下。以最大程度减少来自外界的噪声和地面振动对试验带来的干扰,在鱼箱外表面和底部设有隔音棉。
所述隔音棉6为橡胶海绵。橡胶海绵制品具有优异的缓冲以及隔音性能,将鱼箱外表面均匀的覆盖,并在鱼箱底部衬垫两层橡胶海绵,以最大程度减少来自外界的噪声和地面振动对测试带来的干扰。
实施例2鲫鱼存活率预测模型
1)鱼声信号采集及消噪处理:分别在图1中第一鱼箱和第二鱼箱中注入500L的水,水温度为10~15℃,溶氧量为7-8mg/L,pH为7.2-7.5,将第一水听器和第二水听器布置在水面以下20cm处。将鲫鱼放入1号鱼箱中,鲫鱼的活鱼条数范围为1-50条,最小鱼水比为1:999,最大鱼水比为1:19,静置5min,待鱼在水中比较稳定时,设置声学记录仪进行鱼声信号采集。
设置参数如下:采集时长:1min;采样频率:4000Hz;采集通道:双通道;采集次数:3次。共采集鲫鱼声音信号样本1363个。
2)提取特征参数:根据所述步骤1)消噪处理后的鱼声信号,提取鱼声短时平均能量、鱼声短时平均过零率,再对消噪处理后的鱼声信号进行分解,采用四层、五层、六层和七层小波包分解方法对鱼声信号进行了频段划分,提取各频段能量;
3)构建特征向量:根据所述步骤2)提取的短时平均能量、短时平均过零率和各频段能量构建特征向量;特征向量维数如表1所示。
表1不同分解尺度的特征向量维数
4)样本集划分:对步骤1)采集的鲫鱼声音信号样本划分为训练集和验证集,样本数为1363个。
分别采用了SPXY法和Rank-SPXY法将鲫鱼声音信号样本集按照4:1的比例划分为训练集和验证集,并进行了分析比较,其中Rank-SPXY法的m分别取值为5和10。样本集划分结果如表2所示。Rank-SPXY法,该方法由两部分组成,首先是“Rank”部分,即将样本按因变量(活鱼条数)的升序排序,然后将样本等分为m份;其次是“SPXY法”部分,即在等分的每个区间内采用SPXY法选取出训练集,其余的样本自动归为验证集。其中m也是重要参数,当m=1时,即为SPXY法;当m较大时,得到的训练集活鱼条数更加均匀,但特征值的代表性有所下降。
表2不同样本集划分方法划分结果
由表2可知,采用SPXY法划分的样本集中,29-50条的鲫鱼声音样本全部划分为了训练集,这造成了样本集划分不均匀,而Rank-SPXY法所划分的验证集的数据范围包含在训练集的数据范围内,并且验证集平均值小于训练集平均值,因此选择Rank-SPXY法划分样本集更加合理。在采用Rank-SPXY法对样本集进行划分时,通过比较m=5和m=10时验证集的标准差可知,采用m=10对鲫鱼声音信号样本集的划分更加均匀,但m值较大也会影响特征值的代表性。因此,本实施例同时采用m两种不同取值对鲫鱼声音信号样本集进行划分,并建立预测模型,比较哪种取值更佳。
5)特征值挑选、特征向量降维:对经过步骤4)划分样本集的鲫鱼声音信号进行Z-score(标准分数)标准化预处理,并采用运用竞争自适应重加权采样(CARS)法对鲫鱼声音信号样本集进行特征值优选,并采用10折交叉验证选择模型交叉验证均方差(RMSECV)值最小的特征值变量子集;
运用多元线性回归(MLR)对特征频率再次进行优选,MLR建模后,剔除不显著的特征频率段,得到鱼声信号特征频段能量,对步骤3)特征向量降维,得到降维后的鱼声信号特征向量,优选结果如表3所示。
表3鲫鱼声音信号特征值优选结果
6)模型建立:对所述步骤5)降维后的鱼声信号特征向量,采用多元线性回归(MLR)法和偏最小二乘(PLSR)法分别建立鲫鱼存活率预测模型。预测模型的相关系数如表4所示。
表4预测模型相关系数
综合上述的采用不同特征提取、不同样本集划分方法以及不同建模方法得到的建模结果可知,“7层+Rank-SPXY(m=10)”得到的样本集建立的MLR预测模型的相关系数最高,但计算量是“6层+Rank-SPXY(m=10)”时的二倍,因此选取短时平均能量和短时平均过零率以及基于6层小波包分解的频段能量为鲫鱼声音信号样本的最优特征提取方法;Rank-SPXY(m=10)法为鲫鱼声音信号样本的最优样本集划分方法。对鲫鱼声音信号进行特征提取后,采用Rank-SPXY(m=10)法划分样本集,然后进行Z-score(标准分数)标准化预处理,并采用CARS法对样本集进行特征值优选,最后建立鲫鱼存活率MLR预测模型,回归方程如下:
y=40.471-3.095x2+1.710x4-1.981x5-1.768x6+2.349x7-10.883x11+7.301x12-1.306x16-2.187x18+19.417x25+9.734x28+70.133x35-25.264x43-79.860x47-40.098x50+26.155x55-24.005x61+31.320x62
方程的参数及其显著性见表5,其中,回归常数项b=40.071,xi为鲫鱼声音信号经Z-score标准化预处理后特征值,x1是指鱼声短时平均能量、x2是指鱼声短时平均过零率,xi中i≥3时为基于6层小波包分解的频段能量,具体含义见表6,ai为各特征值的回归系数。
表5回归方程的参数及其显著性
表6鱼声信号短时平均能量、短时平均过零率及鱼六层分解频段能量特征值
鲫鱼存活率预测模型的相关系数R值为0.835,定标标准差RMSECV值为10.096,说明模型具有较好的稳定性和预测性。由表5可知,在x35、x47、x50处,回归系数绝对值最大,其t值相对较大,P值均为0.000,说明这些特征值对预测模型的影响比较显著,其分别代表了鲫鱼声音信号的特征频段为256~264Hz、352~360Hz、376~384Hz。
7)预测存活率:对未知淡水鱼样本鱼声信号进行检测,并将鱼声样本的特征向量经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数;
运用所述步骤6)建立的回归方程,对验证集中的274个鲫鱼声音样本进行预测。将274个鲫鱼声音样本的18个特征值经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数。部分鲫鱼验证集的实际条数以及预测条数见表7。鲫鱼存活率预测模型验证的相关系数R为0.816,校验标准差RMSEP值为8.015,相对分析误差RPD值为1.79,说明该预测模型较可靠。
表7鲫鱼存活率预测结果
由表7可知,靠近两端(1条和50条)的样本预测条数误差较大,在12条左右;而中间的样本预测条数误差较小,在3条左右,模型精度有待进一步提高。
实施例3鳊鱼存活率预测模型
1)鱼声信号采集及消噪处理:分别在图1中第一鱼箱和第二鱼箱中注入500L的水,水温度为10~15℃,溶氧量为7-8mg/L,pH为7.2-7.5,将第一水听器和第二水听器布置在水面以下20cm处。将鲫鱼放入1号鱼箱中,鳊鱼的活鱼条数范围为1-30条,最小鱼水比为1:666,最大鱼水比为1:21,静置5min,待鱼在水中比较稳定时,设置声学记录仪进行鱼声信号采集。
设置参数如下:采集时长:1min;采样频率:4000Hz;采集通道:双通道;采集次数:3次。鳊鱼声音信号样本294个。
2)提取特征参数:同实施例2;
3)构建特征向量:同实施例2;
4)样本集划分:对步骤1)采集的鲫鱼声音信号样本划分为训练集和验证集,样本数为294个。
采用了SPXY法和Rank-SPXY法将鳊鱼声音信号样本集按照4:1的比例划分为训练集和验证集,并进行了分析比较,其中Rank-SPXY法的m分别取值为5和10。样本集划分结果如表8所示。
表8不同样本集划分方法划分结果
由表8可知,采用SPXY法划分的样本集中,训练集和验证集的平均值和标准差均相差较大,而采用Rank-SPXY法划分样本集的平均值和标准差均相差不大,说明采用SPXY法划分样本集造成了划分不均匀,而采用Rank-SPXY法对样本集进行划分更加均匀。因此选择Rank-SPXY法划分鳊鱼声音信号样本集更加合理。在采用Rank-SPXY法对样本集进行划分时,通过比较m=5和m=10时验证集的标准差可知,采用m=10对样本集的划分更加均匀,但m值较大也会影响特征值的代表性。因此本实施例同时采用m两种不同取值对鳊鱼声音信号样本集进行划分。
5)特征值挑选、特征向量降维:对经过步骤4)划分样本集的鲫鱼声音信号进行Z-score(标准分数)标准化预处理,并采用运用竞争自适应重加权采样(CARS)法对对鳊鱼声音信号样本集进行特征值优选,并采用10折交叉验证选择模型交叉验证均方差(RMSECV)值最小的特征值变量子集。
运用MLR对特征频率再次进行优选,MLR建模后,剔除不显著的特征频率段,得到鱼声信号特征频段能量,对步骤3)特征向量降维,得到降维后的鱼声信号特征向量,优选结果如表9所示。
表9鳊鱼声音信号特征值优选结果
6)模型建立:对所述步骤5)降维后的鱼声信号特征向量,采用多元线性回归(MLR)法和偏最小二乘(PLSR)法分别建立鳊鱼存活率预测模型。预测模型的相关系数如表10所示。
表10预测模型相关系数
综合上述的采用不同特征提取、不同样本集划分方法以及不同建模方法得到的建模结果可知,“7层+Rank-SPXY(m=5)”得到的样本集建立MLR预测模型的相关系数最高,但计算量是“6层+Rank-SPXY(m=5)”时的二倍,因此本实施例选取短时平均能量和短时平均过零率以及基于6层小波包分解的频段能量为鳊鱼声音信号样本的最优特征提取方法;Rank-SPXY(m=5)法为鳊鱼声音信号样本的最优样本集划分方法。对鳊鱼声音信号进行特征提取后,采用Rank-SPXY(m=5)法划分样本集,然后进行Z-score标准化预处理,并采用CARS法对样本集进行特征值优选,最后建立鳊鱼存活率MLR预测模型,回归方程如下:
y=4.384+1.415x5+2.681x9+8.356x14+2.694x18-3.290x24
方程的参数及其显著性见表11,其中,回归常数项b=4.384,xi为鳊鱼声音信号经预处理后的特征值,ai为各特征值的回归系数。
表11鳊鱼存活率预测模型的参数及其显著性
鳊鱼存活率模型的相关系数R值为0.894,定标标准差RMSECV值为3.83,说明模型具有较好的稳定性和预测性。由表4-10可知,在x9、x14、x18、x24处,回归系数绝对值最大,其t值相对较大,P值最大为0.034,小于0.05,说明这些特征值对预测模型的影响比较显著,其分别代表了鳊鱼声音信号的特征频段为48~56Hz、88~96Hz、120~128Hz、168~178Hz。
7)预测存活率:对未知淡水鱼样本鱼声信号进行检测,并将鱼声样本的特征向量经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数;
运用步骤6)所建立的回归方程,对验证集中60个鳊鱼声音样本进行预测。将60个鳊鱼声音样本的5个特征值经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数。部分鳊鱼验证集的实际条数以及预测条数见表12。鳊鱼存活率预测模型验证的相关系数R为0.865,校验标准差RMSEP值为4.54,相对分析误差RPD值为2.01,说明该预测模型非常可靠。
表12鳊鱼存活率预测结果
由表12可知,偏差最小为0条,最大为7条,这种误差是由于噪音干扰造成,有待进一步提高信号的信噪比,进而提高模型的预测精度。
本实施例建立了淡水鱼鳊鱼存活率预测模型,运用模型对验证集样本进行了预测,并研究了不同分解尺度的频段分解以及不同样本集划分方法对存活率预测模型性能的影响,结果表明:采用“短时平均能量+短时平均过零率+6层小波包分解各频段能量”提取特征结合Rank-SPXY(m=5)样本划分方法建立的鳊鱼存活率预测模型预测性能最优(R=0.894,RPD=2.01)。
上述的实施例仅为本发明的优选技术方案,而不应视为对于本发明的限制,本发明中的实施例及实施例中的特征在不冲突的情况下,可以相互任意组合。本发明的保护范围应以权利要求记载的技术方案,包括权利要求记载的技术方案中技术特征的等同替换方案为保护范围。即在此范围内的等同替换改进,也在本发明的保护范围之内。
Claims (10)
1.一种基于被动声学信息的淡水鱼存活率预测装置,提取原始鱼声信号的装置包括第一鱼箱(1),第二鱼箱(2),其特征在于,所述第一鱼箱(1)内设有第一水听器(4),所述第二鱼箱(2)内设有第二水听器(5),所述第一水听器(4)、第二水听器(5)分别与声学记录仪(3)连接。
2.根据权利要求1所述基于被动声学信息的淡水鱼存活率预测装置,其特征在于:所述第一鱼箱(1)、第一水听器(4)、第二水听器(5)设有若干个。
3.根据权利要求1所述基于被动声学信息的淡水鱼存活率预测装置,其特征在于:所述第一水听器(4)、第二水听器(5)分别位于第一鱼箱(1)、第二鱼箱(2)水面以下,所述第一鱼箱(1)、第二鱼箱(2)外设有隔音棉(6)。
4.一种基于被动声学信息的淡水鱼存活率预测方法,其特征在于:所述方法具体包括以下步骤:
1)消噪处理:提取原始鱼声信号,并对所得的鱼声信号进行消噪处理;
2)提取特征参数:根据所述步骤1)消噪处理后的鱼声信号,提取鱼声短时平均能量、鱼声短时平均过零率,再对消噪处理后的鱼声信号进行分解,采用四层、五层、六层和七层小波包分解方法对鱼声信号进行了频段划分,提取各频段能量;
3)构建特征向量:根据所述步骤2)提取的短时平均能量、短时平均过零率和各频段能量构建特征向量;
4)样本集划分:采集不同条数不同时间的鱼声信号,构成一个鱼声样本集,将样本集划分为训练集和验证集;
5)特征频段挑选、特征向量降维:对经过步骤4)划分样本集的鱼声音信号进行Z-score标准化预处理,并采用竞争自适应重加权采样法、多元线性回归对特征频段再次进行优选,剔除不显著的特征频段,得到鱼声信号特征频段能量,对步骤3)特征向量降维,得到降维后的鱼声信号特征向量;
6)模型建立:对所述步骤5)降维后的鱼声信号特征向量,采用多元线性回归法和偏最小二乘法分别建立鱼存活率预测模型;
7)预测存活率:对未知淡水鱼样本鱼声信号进行检测,并将鱼声样本的特征向量经Z-score标准化预处理后带入到回归方程中,计算出活鱼条数;
完成鱼存活率的预测。
5.根据权利要求4所述的方法,其特征在于:所述步骤4)采用Rank-SPXY法进行样本集的划分,该方法由两部分组成,首先是Rank部分,即将样本按因变量活鱼条数的升序排序,然后将样本等分为m份;其次是SPXY法部分,即在等分m份的每个区间内采用SPXY法选取出训练集,其余的样本自动归为验证集,m分别取值为5和10。
6.根据权利要求5所述的方法,其特征在于:所述步骤6)鲫鱼存活率预测模型,选取短时平均能量、短时平均过零率以及基于6层小波包分解的频段能量为鲫鱼声音信号样本的最优特征提取方法,Rank-SPXY法样本集划分方法中m=10,建立鲫鱼存活率MLR预测模型,回归方程如下:
y=40.471-3.095x2+1.710x4-1.981x5-1.768x6+2.349x7-10.883x11+7.301x12-1.306x16-2.187x18+19.417x25+9.734x28+70.133x35-25.264x43-79.860x47-40.098x50+26.155x55-24.005x61+31.320x62,
y为存活率,xi为鲫鱼声音信号经Z-score标准化预处理后特征值,鱼声短时平均能量x1、鱼声短时平均过零率x2,xi中i≥3时为基于6层小波包分解的频段能量,8~16Hz频段能量x4,16~24Hz频段能量x5,24~32Hz频段能量x6,32~40Hz频段能量x7,64~72Hz频段能量x11,72~80Hz频段能量x12,104~112Hz频段能量x16,120~128Hz频段能量x18,176~184Hz频段能量x25,200~208Hz频段能量x28,256~264Hz频段能量x35,320~328Hz频段能量x43,352~360Hz频段能量x47,376~384Hz频段能量x50,416~424Hz频段能量x55,464~472Hz频段能量x61,472~480Hz频段能量x62。
7.根据权利要求6所述的方法,其特征在于:鲫鱼声音信号的特征频段为256~264Hz、352~360Hz、376~384Hz。
8.根据权利要求5所述的方法,其特征在于:所述步骤6)鳊鱼存活率预测模型,选取短时平均能量、短时平均过零率以及基于6层小波包分解的频段能量为鳊鱼声音信号样本的最优特征提取方法,Rank-SPXY法样本集划分方法中m=5,建立鳊鱼存活率MLR预测模型,回归方程如下:
y=4.384+1.415x5+2.681x9+8.356x14+2.694x18-3.290x24,
y为存活率,xi为鳊鱼声音信号经Z-score标准化预处理后特征值,鱼声短时平均能量x1、鱼声短时平均过零率x2,xi中i≥3时为基于6层小波包分解的频段能量,16~24Hz频段能量x5,48~56Hz频段能量x9,88~96Hz频段能量x14,120~128Hz频段能量x18,168~178Hz频段能量x24。
9.根据权利要求8所述的方法,其特征在于:鳊鱼声音信号的特征频段为48~56Hz、88~96Hz、120~128Hz、168~178Hz。
10.权利要求1-37任意一项所述的装置或权利要求2-9任意一项所述方法在淡水鱼养殖、活鱼运输以及渔业资料调查中的应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610801989.5A CN106376511B (zh) | 2016-09-05 | 2016-09-05 | 一种基于被动声学信息的淡水鱼存活率预测装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610801989.5A CN106376511B (zh) | 2016-09-05 | 2016-09-05 | 一种基于被动声学信息的淡水鱼存活率预测装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106376511A true CN106376511A (zh) | 2017-02-08 |
CN106376511B CN106376511B (zh) | 2019-04-02 |
Family
ID=57939060
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610801989.5A Active CN106376511B (zh) | 2016-09-05 | 2016-09-05 | 一种基于被动声学信息的淡水鱼存活率预测装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106376511B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107942027A (zh) * | 2018-01-08 | 2018-04-20 | 江苏省农业科学院 | 一种蟹类水下生活环境适宜性检测方法及装置 |
CN109566478A (zh) * | 2018-10-18 | 2019-04-05 | 厦门大学 | 水产养殖生物体健康状态的在线监测方法 |
CN109934370A (zh) * | 2017-12-15 | 2019-06-25 | 中国农业大学 | 一种鱼类存活率预测方法及装置 |
CN111721754A (zh) * | 2020-06-24 | 2020-09-29 | 天津大学 | 基于激光诱导击穿光谱的液态奶中矿物质元素检测方法 |
-
2016
- 2016-09-05 CN CN201610801989.5A patent/CN106376511B/zh active Active
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109934370A (zh) * | 2017-12-15 | 2019-06-25 | 中国农业大学 | 一种鱼类存活率预测方法及装置 |
CN109934370B (zh) * | 2017-12-15 | 2024-03-05 | 中国农业大学 | 一种鱼类存活率预测方法及装置 |
CN107942027A (zh) * | 2018-01-08 | 2018-04-20 | 江苏省农业科学院 | 一种蟹类水下生活环境适宜性检测方法及装置 |
CN107942027B (zh) * | 2018-01-08 | 2023-03-28 | 江苏省农业科学院 | 一种蟹类水下生活环境适宜性检测方法及装置 |
CN109566478A (zh) * | 2018-10-18 | 2019-04-05 | 厦门大学 | 水产养殖生物体健康状态的在线监测方法 |
CN111721754A (zh) * | 2020-06-24 | 2020-09-29 | 天津大学 | 基于激光诱导击穿光谱的液态奶中矿物质元素检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106376511B (zh) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106376511B (zh) | 一种基于被动声学信息的淡水鱼存活率预测装置及方法 | |
Jech et al. | A multifrequency method to classify and evaluate fisheries acoustics data | |
CN106417143B (zh) | 一种基于被动声学信息的淡水鱼品种识别装置和方法 | |
Korneliussen et al. | An operational system for processing and visualizing multi-frequency acoustic data | |
Kloser et al. | Species identification in deep water using multiple acoustic frequencies | |
Martin et al. | Acoustic classification of zooplankton | |
Smith et al. | The use of passive acoustics to measure feed consumption by Penaeus monodon (giant tiger prawn) in cultured systems | |
Herman et al. | Detection and abundance estimation of euphausiids using an optical plankton counter | |
CN112578725A (zh) | 一种基于物联网的渔业水上养殖设施安全预警系统及方法 | |
CN110243936B (zh) | 一种原位无损伤测定沉水植物生物量的方法 | |
Tracey et al. | Fish species composition on seamounts and adjacent slope in New Zealand waters | |
CN109479750A (zh) | 一种基于声学信息的梅山猪发情监测方法 | |
Fang et al. | The source parameters of echolocation clicks from captive and free-ranging Yangtze finless porpoises (Neophocaena asiaeorientalis asiaeorientalis) | |
CN115578678A (zh) | 一种鱼类摄食强度分类方法及系统 | |
Overholtz et al. | Empirical comparisons of survey designs in acoustic surveys of Gulf of Maine-Georges Bank Atlantic herring | |
Hassan et al. | Fish as a source of acoustic signal measurement in an aquaculture tank: Acoustic sensor based time frequency analysis | |
Barange et al. | Determination of composition and vertical structure of fish communities using in situ measurements of acoustic target strength | |
CN116541672A (zh) | 基于声学信息和水面波动信息的鱼类养殖投喂量决策方法 | |
ARNAYA et al. | Studies on Acoustic Target Strength of Squid: Ⅲ. Measurement of the mean target strength of small live squid | |
Baidai et al. | Recent advances on the use of supervised learning algorithms for detecting tuna aggregations under fads from echosounder buoys data | |
CN115690777A (zh) | 基于成熟度识别装置的香菇菌棒成熟度识别方法及系统 | |
Stanton et al. | A field examination of acoustical scattering from marine organisms at 70 kHz | |
Li et al. | SIMULATION OF FISH MIGRATION AT DIFFERENT WATER DEPTHS BASED ON BACKPROPAGATION NEURAL NETWORK. | |
CN112883773A (zh) | 一种基于声学影像数据评估的物种判别方法 | |
CN105486758B (zh) | 一种霉变玉米种子检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |