CN113607068B - 一种光声测量信号特征的识别模型建立方法及提取方法 - Google Patents
一种光声测量信号特征的识别模型建立方法及提取方法 Download PDFInfo
- Publication number
- CN113607068B CN113607068B CN202110814894.8A CN202110814894A CN113607068B CN 113607068 B CN113607068 B CN 113607068B CN 202110814894 A CN202110814894 A CN 202110814894A CN 113607068 B CN113607068 B CN 113607068B
- Authority
- CN
- China
- Prior art keywords
- photoacoustic measurement
- measurement signal
- dimensional time
- time sequence
- signals
- 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.)
- Active
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 190
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000012549 training Methods 0.000 claims abstract description 45
- 238000003062 neural network model Methods 0.000 claims abstract description 37
- 239000007787 solid Substances 0.000 claims abstract description 34
- 238000000605 extraction Methods 0.000 claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 238000012795 verification Methods 0.000 claims abstract description 10
- 238000011176 pooling Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 20
- 238000002372 labelling Methods 0.000 claims description 14
- 239000012528 membrane Substances 0.000 claims description 12
- 230000009467 reduction Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 12
- 239000000523 sample Substances 0.000 description 40
- 230000000694 effects Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000004913 activation Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000010606 normalization Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000010895 photoacoustic effect Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 238000013136 deep learning model Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000011343 solid material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
- G01B11/0616—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material of coating
-
- 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/067—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using optical means
- G06N3/0675—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using optical means using electro-optical, acousto-optical or opto-electronic means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Neurology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种光声测量信号特征的识别模型建立方法及提取方法,属于信号处理领域,包括:获得测量自不同种类的固体样品的多组一维时序光声测量信号,并标注出每组信号中的各特征信号所在的时域区间和类型,得到多个样本;将U‑Net模型中所有的卷积层、最大池化层和上采样均修改为1D形式,以建立一维U‑Net神经网络模型,用于预测一维时序光声测量信号中各类特征信号的时域分布概率曲线;将所有的样本划分为训练集、验证集和测试集,对一维U‑Net神经网络模型进行训练、验证和测试后,得到光声测量信号特征的识别模型。本发明能够有效解决现有的一维时序光声测量信号特征识别与提取方法的鲁棒性差、准确度低、无法区分不同膜层回声信号的技术问题。
Description
技术领域
本发明属于信号处理领域,更具体地,涉及一种光声测量信号特征的识别模型建立方法及提取方法。
背景技术
光声效应(Photoacoustic Effect)是指当媒质受到周期性强度调制的光源照射时,内部温度改变引起区域结构和体积变化从而产生声信号的现象。基于光声效应的测量技术结合了光学测量的高分辨率和声学测量的高穿透性优势,已被广泛应用于生物医学、军工航天、半导体产业等多个领域,而该技术实现良好应用的关键在于对测量信号进行准确的特征识别与提取。
面向非透明固体材料的光声测量技术是实现快速无损缺陷探测、膜厚测量、物性表征的有效途径,然而目前关于固体光声测量中一维时序信号处理方法的研究较少,而在生物医学领域中的二维光声图像的处理算法相对成熟,主要包括光声图像的重构和反演、光声图像的分类分割算法等,如发明专利CN201911341072.1中公开的一种采用ResNet对光声图像进行超分辨率重构的技术,以此来提高图像所能提供的生物生理信息。在论文 U-Net:Convolutional Networks for Biomedical Image Segmentation中,Long 等人提出一种专门针对小数据集医学图像的语义分割模型,能够提供对图像像素级别的分类分割,从而对图像中的关键特征成分进行标注和区分。
基于固体光声效应的无损探伤、膜厚测量等技术通常是利用脉冲激光在样品内部近表面处激发出声脉冲,并通过超声换能器、光学探测等方法对此声脉冲在样品中的传播过程进行监测,获得一维时序光声测量信号。声脉冲产生的瞬间会形成零点信号,其传播过程中会在样品界面处发生部分反射从而回到样品表面,形成回声信号,通过对时序测量信号中的零点信号和回声信号进行特征识别与提取,能够反演计算出样品的膜厚等信息。
现有的方法往往利用求导或极值函数抓取时序测量信号中的零点信号和回声信号,但对信号的信噪比要求较高,存在鲁棒性差和准确率低的问题,此外,这类方法无法对多层膜样品中不同膜层反射的回声信号进行区分,难以实现对多层膜的膜厚测量。
发明内容
针对现有技术的缺陷和改进需求,本发明提供了一种光声测量信号特征的识别模型建立方法及提取方法,旨在解决固体光声测量技术中现有一维时序信号特征识别与提取方法的鲁棒性差、准确度低,以及无法区分多膜层样品中不同膜层回声信号的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种光声测量信号特征的识别模型建立方法,包括:
分别以不同种类的固体样品为测量对象,获得各测量对象的至少一组一维时序光声测量信号,并标注出每一组一维时序光声测量信号中的每一个特征信号所在的时域区间和对应的类型;不同种类的固体样品的膜层数量均为N,N为正整数,特征信号包括:零点信号和各膜层返回的回声信号;
将U-Net神经网络模型中所有的卷积层、最大池化层和上采样均修改为1D形式,以建立一维U-Net神经网络模型,用于预测目标一维时序光声测量信号中各类特征信号的时域分布概率曲线;目标一维时序光声测量信号由具有N个膜层的固体样品在光声测量中产生;
将每一组一维时序光声测量信号及其对应的标注结果作为一个样本,将所有的样本划分为训练集、验证集和测试集,利用训练集、验证集和测试集分别对一维U-Net神经网络模型进行训练、验证和测试后,得到光声测量信号特征的识别模型。
本发明在用于图像处理的U-Net神经网络模型的基础上,建立了可用于处理一维时序光声测量信号的一维U-Net神经网络模型,并构建了相应的数据集对该模型进行训练、验证和测试,使该模型可以预测一维时序光声测量信号中各类特征信号的时域分布概率曲线;一方面,一维U-Net神经网络模型为深度学习模型,其在训练过程中,能够学习到数据内在的规律,对于复杂噪声环境下获取的一维时序光声测量信号也能准确识别出其中的特征信号,具有较高的准确率和鲁棒性;另一方面,由于构建训练数据时,对一维时序光声测量信号中不同膜层的回声信号都进行了分类标注,因此U-Net神经网络模型经过训练之后,能够分别预测出每一个膜层返回的回声信号的时域分布概率,基于该预测结果,可以准确识别出一维时序光声测量信号中的每一个膜层返回的回声信号所在的时域区间,从而完成多层膜样品中不同膜层回声信号的分类识别和提取。
进一步地,一维U-Net神经网络模型中,下采样和上采样部分中卷积层的卷积核均大于或等于3×1,由此能够获得最好的预测效果。
进一步地,建立一维U-Net神经网络模型时,还包括:在U-Net神经网络模型的下采样部分和/或上采样部分中,加入至少一个dropout层;通过在上采样部分和/或下采样部分中加入dropout层,能够有效防止过拟合。
进一步地,在对一维时序光声测量信号进行标注之前,还包括:对一维时序光声测量信号进行数据增强;
数据增强包括:添加脉冲噪声,和/或添加低频扰动;
低频扰动的频率小于一维时序光声测量信号中特征信号的频率最小值,且低频扰动的幅值小于一维时序光声测量信号中特征信号的幅值最小值。
实际的光声测量环境十分复杂,往往存在噪声和干扰,而在构建训练数据时,很难覆盖到全部的场景,而且目前在光声测量领域,可用的训练数据较少;本发明对一维时序光声测量信号添加脉冲噪声和添加低频扰动,一方面,能够使训练数据尽可能覆盖到全部的场景,提高训练所得的模型的泛化能力,进一步增强后续特征信号识别与提取的鲁棒性和准确性;另一方面可以扩充训练数据量,提高模型的训练效果。
进一步地,数据增强还包括:平移变换,和/或裁剪变换。
本发明通过对模型进行平移变换和裁剪变换,可以进一步扩充训练数据量,提高模型的训练效果。
进一步地,在对一维时序光声测量信号进行标注之前,还包括:对一维时序光声测量信号进行预处理;预处理包括以下操作中的至少一项:
在对一维时序光声测量信号进行数据增强之前,对一维时序光声测量信号进行数字降噪;通过数字降噪处理,可以有效避免噪声对模型训练效果的干扰,进一步降低对信号的信噪比要求,提高模型的鲁棒性和准确率;
在对一维时序光声测量信号进行数据增强之后,对一维时序光声测量信号进行标准化;通过标准化处理,可以加速模型的收敛。
进一步地,数字降噪为小波去噪;通过小波去噪,能够有效滤除一维时序光声测量信号中的噪声。
按照本发明的另一个方面,提供了一种光声测量信号特征的提取方法,包括:
对于具有N个膜层的待测固体样品,获得其一维时序光声测量信号,将一维时序光声测量信号输入至光声测量信号特征的识别模型,以预测一维时序光声测量信号中各类特征信号的时域分布概率曲线;光声测量信号特征的识别模型由本发明提供的光声测量信号特征的识别模型建立方法得到;
对于每一类特征信号,提取出其时域分布概率曲线中大于相应概率阈值的曲线部分,并将所提出的曲线部分所对应的时域区间识别为该类特征信号所在的时域区间,从一维时序光声测量信号中提取出该时域区间内的信号,完成对该类特征信号的提取;
其中,N为正整数。
进一步地,本发明提供的光声测量信号特征的提取方法,一维时序光声测量信号输入至光声测量信号特征的识别模型之前,还包括:对一维时序光声测量信号进行数字降噪和标准化。
按照本发明的又一个方面,提供了一种计算机可读存储介质,包括存储的计算机程序;计算机程序被处理器执行时,控制计算机可读存储介质所在设备执行本发明提供的光声测量信号特征的识别模型建立方法,和/或本发明提供的光声测量信号特征的提取方法。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
(1)本发明在用于图像处理的U-Net神经网络模型的基础上,建立了可用于处理一维时序光声测量信号的一维U-Net神经网络模型,并构建了相应的数据集对该模型进行训练、验证和测试,使该模型可以预测一维时序光声测量信号中各类特征信号的时域分布概率曲线,对于复杂噪声环境下获取的一维时序光声测量信号也能准确识别和提取出其中的特征信号,具有较高的准确率和鲁棒性,并且能够分别预测出每一个膜层返回的回声信号的时域分布概率,基于该预测结果,可以准确识别出一维时序光声测量信号中的每一个膜层返回的回声信号所在的时域区间,从而完成多层膜样品中不同膜层回声信号的分类识别和提取。总体而言,本发明能够有效解决固体光声测量技术中现有一维时序信号特征识别与提取方法的鲁棒性差、准确度低,以及无法区分多膜层样品中不同膜层回声信号的技术问题。
(2)本发明在构建训练数据时,会通过对一维时序光声测量信号添加脉冲噪声和添加低频扰动来进行数据增强,能够使训练数据尽可能覆盖到全部的场景,提高训练所得的模型的泛化能力,进一步增强后续特征信号识别与提取的鲁棒性和准确性,并且能够扩充训练数据量,提高模型的训练效果。
(3)目前光声测量还处在起步阶段,可用的训练数据较少,本发明所建立的一维U-Net神经网络模型,是在经典U-Net神经网络模型的基础上改进得到的,由于U-Net神经网络模型即使在小数据集的情况下,依然具有较好的预测效果,因此,本发明所建立的模型在训练数据较少的情况下也具有较好的特征识别效果;本发明在构建数据时,进行了数据增强操作,能够有效解决训练数据过少的问题,进一步提高训练后模型的预测效果。
附图说明
图1为本发明实施例提供的光声测量信号特征的识别模型建立方法流程图;
图2为本发明实施例提供的一维时序光声测量信号及利用Labelme进行分类标记后的标注结果示意图;
图3为本发明实施例提供的一维U-Net神经网络模型示意图;
图4为本发明实施例提供的一维时序光声测量信号的各类特征信号以及背景信号的时域分布概率曲线示意图;
图5为本发明实施例提供的一维时序光声测量信号的特征提取结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明中,本发明及附图中的术语“第一”、“第二”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。
为了解决固体光声测量技术中现有一维时序信号特征识别与提取方法的鲁棒性差、准确度低,以及无法区分多膜层样品中不同膜层回声信号的技术问题,本发明提供了一种光声测量信号特征的识别模型建立方法及提取方法,其整体思路在于:建立可用于处理一维时序光声测量信号的深度学习模型,并构建相应的训练数据对模型进行训练、验证和测试,以通过深度学习的方式,预测光声测量的一维时序光声测量信号中每一类特征信号在时域上的分布概率曲线,使得模型对于复杂环境下获取的一维时序光声测量信号也能准确完成特征的识别,并且可以区分每一个膜层返回的回声信号,从而完成一维时序光声测量信号中特征信号的识别和提取。
本发明所提供的光声测量信号特征的识别模型建立方法及提取方法,可适用于具有任意数量膜层的固体样品,为便于描述,若无特殊说明,在以下实施例中,均以双膜层固体样品作为测量对象。
以下为实施例。
实施例1:
一种光声测量信号特征的识别模型建立方法,如图1所示,包括以下步骤(S1)~(S7):
(S1)分别以不同种类的固体样品为测量对象,获得各测量对象的至少一组一维时序光声测量信号,共得到n组一维时序光声测量信号:
本实施例,不同种类的固体样品的膜层数量均为N=2,固体样品具体是金属双层膜样品;
利用现有的光声测量系统,即可获取到测量对象的一维时序光声测量信号;可选地,本实施例中,所选用的光声测量系统包括光电探测器和锁相放大器;光声测量的具体过程为:将一束脉冲激光照射到固体样品表面测量点的位置,用于在样品内部近表面处激发出声脉冲,此脉冲激光称作激发光;同时用另一束激光照射到样品表面相同位置处,用于对样品表面变化进行实时监测,此激光称作探测光,经样品表面反射后的探测光由光电探测器接收并转化为电信号,再经过锁相放大器完成降噪及信号放大输出为一维时序光声测量信号;应当说明的是,此处关于光声测量的描述,仅为本发明可选的实施方式,不应理解为对本发明的唯一限定,在本发明其他的一些实施例中,也可以使用其他的测量系统或方法获取测量对象的一维时序光声测量信号;
对于同一固体样品,通过更换测量点,即可获得该固体样品的多组一维时序光声测量信号;
本实施例中,固体样品的种类不同,具体是指固体样品的膜层材料、膜层厚度等不同;
为了为后续训练神经网络提供足够的训练数据,保证模型训练效果,本实施例中,所获取的一维时序光声测量信号的组数为n≥30。
(S2)对每一组一维时序光声测量信号进行数字降噪处理,保证合理的信噪比:
通过数字降噪处理,可以有效避免噪声对模型训练效果的干扰,使一维时序测量数据达到较高的信噪比,从而一维时序光声测量信号中的特征信号更加明显,从而有利于信号特征识别的准确性,提高模型的鲁棒性和准确率;
为了达到较好的降噪效果,作为一种优选的实施方式,本实施例中,所选用的数字降噪为小波去噪,具体地,选择sym4小波基进行6阶小波分解,重构规则采用了Bayes估计阈值,阈值函数选用中值阈值函数;
应当说明的是,此处的小波去噪仅为本发明优选的数字降噪方式,不应理解为对本发明的唯一限定,在本发明其他的一些实施例中,也可以使用经验模态分解(empiricalmode decomposition,EMD)等其他方式进行数字降噪处理;
(S3)对降噪后的一维时序光声测量信号进行数据增强及标准化:
从测量的时间成本考虑,S1步骤中获取的一维时序光声测量信号的组数往往是有限的,且很难覆盖掉全部的场景,因此需要通过数据增强的方式来扩充测量信号的样本规模,进而提高后续训练出的神经网络模型的泛化能力;
可选地,本实施例中,对于每一组一维时序光声测量信号都会进行四种数据增强操作,具体为平移变换、剪裁变换、添加脉冲噪声以及添加低频扰动;
通过数据增强,训练数据从n组一维时序光声测量信号扩展到了n×4 组一维时序光声测量信号;
考虑到信号为一维信号,按照xi_stan=(xi-μ)/σ对数据增强后的一维时序光声测量信号进行标准化处理,其中xi代表单组测量信号中的第i个数据,μ代表该组一维时序光声测量信号中所有数据的均值,σ代表该组一维时序光声测量信号中所有数据的方差,而xi_stan代表标准化后该组一维时序光声测量信号中的第i个数据;
可选地,最后将经过标准化的一维时序光声测量信号保存为.mat格式的signal数据文件,总共获得n×4个signal数据文件。
(S4)标注出每一组一维时序光声测量信号中的每一个特征信号所在的时域区间和对应的类型;特征信号包括:零点信号和各膜层返回的回声信号:
声脉冲在样品内部近表面处产生的瞬间,会引起样品表面光学性质的剧烈变化,从而使探测光的光强发生明显变化,反映到一维时序光声测量信号中表现为零点信号;声脉冲传播的过程中会在样品膜层界面处发生部分反射从而回到样品表面,引起样品表面光学性质的变化以及探测光光强的波动,反映到一维时序光声测量信号中表现为回声信号;一维时序光声测量信号中的零点信号和回声信号共同构成了一维时序光声测量信号的特征信号;本实施例中,一维时序光声测量信号的特征信号具体有三种,分别是零点信号、自第一个膜层与第二个膜层交界处反射的回声信号以及自第二个膜层与基底交界处反射的回声信号;
对于已知的固体样品,其中每一类特征信号的形貌已知,且相同膜层返回的回声信号的间距相同,基于此,在保证所标注的时域区间内包含特征信号的完整形貌特征,且每个特征信号的时域区间不与其他特征信号的时域区间重叠的情况下,即可完成对一维时序光声测量信号中的每一个特征信号所在的时域区间和对应的类型的标注;
在实际应用中,上述标注工作可直接利用深度学习图像标注工具 Labelme软件完成,将用于保存一维时序光声测量信号的signal数据文件转换成图像文件后,在Labelme中利用多边形工具对图像中上述三种特征信号所在区域进行四边形框选,图2所示为一组一维时序光声测量信号,及利用Labelme进行分类标记后的标注结果示意图,其中,零点信号的类型标记为zero,第一个膜层与第二个膜层交界处反射的回声信号以及自第二个膜层与基底交界处反射的回声信号的类型分别标记为echo1和echo2;
本实施例中,标注完成后,会将所得到的标注结果存储为.json文件,共得到与n×4个signal数据文件一一对应的n×4个.json文件;
上述.json文件中包含了标记框各个顶点的横纵坐标信息,但实际上只需关注其中的横坐标信息,也即特征信号的时域信息,并将其解析为数据长度与对应的signal数据文件一致的一维数组label,从而获得n×4个与 signal数据文件一一对应的label数据文件。实际上,一维时序光声测量信号中除前述三种特征信号以外的部分皆为背景信号(background),因此标记框以外的横坐标区域就对应了背景信号的时域信息,而这一信息也会被包含到label数据文件中。
(S5)将每一组一维时序光声测量信号及其对应的标注结果作为一个样本,将所有的样本划分为训练集、验证集和测试集:
将对应的一组signal文件和label文件中的数据提取出来,便得到了一个样本;
可选地,本实施例中,训练集、验证集和测试集的数据比例为7:2:1,相应地,n×4个signal数据文件和n×4个label数据文件分别按照此比例随机抽样存放于训练集、验证集和测试集中。
(S6)建立一维U-Net神经网络模型,用于预测目标一维时序光声测量信号中各类特征信号的时域分布概率曲线:
本实施例中,目标一维时序光声测量信号由具有N(N=2)个膜层的固体金属样品在光声测量中产生;
本实施例中,一维U-Net神经网络模型是在用于处理二维图像的经典 U-Net神经网络模型的基础上改进得到的,其建立方式包括:
将U-Net神经网络模型中所有的卷积层、最大池化层和上采样均修改为1D形式;
本实施例中,下采样和上采样部分中卷积层的卷积核大小为3×1,并且一维U-Net神经网络模型的建立还包括以下操作:
在U-Net神经网络模型的下采样和/或上采样部分中,加入至少一个 dropout层;通过在下采样和/或上采样部分中加入dropout层,能够有效防止过拟合;
容易理解的是,为了支持对不同数据长度的一维时序光声测量信号的特征识别与提取,在进行下采样之前,需要将一维时序光声测量信号调整至统一的长度(resize);在实际应用中,可以通过数据预处理的方式完成数据长度的统一,也可以直接在一维U-Net神经网络模型的输入层中完成数据长度的统一;可选地,本实施例中,由U-Net神经网络模型的输入层完成数据长度的统一,且具体使用线性插值的方式调节数据的长度;应当说明的是,在调节数据长度时,也可以采用多项式插值、二次插值、Cubic 插值等其他的插值方式;
最终,本实施例所建立的一维U-Net神经网络模型如图3所示,首先,输入信号在输入层进行resize,具体地,是通过线性插值的方法将数据大小统一为2360*1,然后进入下采样过程。下采样过程中采用1D卷积配合最大池化实现特征提取,其中,设置卷积层的卷积核大小为3×1,步长为1, padding方式为same,激活函数为ReLU;设置池化层大小为2×1,步长为 2;经过3次特征提取完成下采样的信号会进入上采样过程;上采样过程中采用upsampling1D方法配合concatenate方法以及1D卷积实现数据扩展,具体地,信号经upsampling处理后需先进行卷积核大小为2×1、步长为1、 padding方式为same、激活函数为ReLU的卷积,再与下采样过程中对称位置的信号直接进行拼接操作(concatenate),然后再进行卷积核大小为3×1、步长为1、padding方式为same、激活函数为ReLU的卷积;经过3次数据扩展完成上采样的信号会最后经过一次卷积核大小为1×1、步长为1、 padding方式为same、激活函数为Softmax的卷积,最终输出信号特征识别及分类的结果——数据长度为2360的四种信号(前述三种特征信号以及背景信号)的预测概率曲线;除此之外,本施例还在下采样和上采样过程中适当加入了dropout,以防止过拟合。
(S7)利用训练集、验证集和测试集分别对一维U-Net神经网络模型进行训练、验证和测试后,得到光声测量信号特征的识别模型。
上述神经网络模型构建完成后,将训练集中的signal数据文件中的数据作为该模型的输入信号,同时将模型输出结果与对应的label数据文件中的数据进行对比、计算损失,并据此对模型权重等学习参数进行反向传播的参数更新,以实现对模型的训练;
同时,利用验证集对训练出的模型进行验证,并根据训练及验证结果调整卷积核大小等模型超参数、特征提取次数和dropout位置等模型结构特征以及epoch和learningrate等训练参数,以优化模型的训练效果;
最后,利用测试集对训练并验证后的模型进行测试,判断该模型信号特征识别及分类的准确率是否达到标准,如未达标则需继续调整模型超参数、结构以及训练参数,再重新进行训练、验证和测试,直至模型的测试结果达标,从而得到光声测量信号特征的识别模型;该模型可用于预测具有2个膜层的固体样品的一维时序光声测量信号中各类特征信号的时域分布的概率曲线;
上述训练、验证及测试过程中,采用交叉熵损失作为损失函数,采用 accuracy作为模型优劣的评价指标;
可选地,本实例中,在得到光声测量信号特征的识别模型后,将该模型保存为.h5文件。
应当说明的是,当固体样品的膜层数量N发生变化时,其一维时序光声测量信号中的特征信号种类也会发生变化,具体为N+1种,即零点信号和N个膜层中每个膜层返回的回声信号,此时在构建训练数据时,进行标注即可;最终所得到的模型即可用于预测具有N个膜层的固体样品的一维时序光声测量信号中各类特征信号的时域分布的概率曲线。
实施例2:
一种光声测量信号特征的提取方法,包括:
对于具有2个膜层的待测固体金属样品,获得其在光声测量中产生的一维时序光声测量信号,对其进行数字降噪和标准化;数字降噪和标准化的方法可参考上述实施例1中的描述;
将数字降噪和标准化之后的一维时序光声测量信号输入至光声测量信号特征的识别模型,以获得一维时序光声测量信号中各类特征信号的时域分布概率曲线;光声测量信号特征的识别模型由上述实施例1提供的光声测量信号特征的识别模型建立方法得到,实际应用中,调用相应的.h5文件,即可获得模型;本实施例中,模型预测得到的各类特征信号以及背景信号的时域分布概率曲线如图4所示;
对于每一类特征信号,设定相应的概率阈值,若某一信号类别的预测概率大于该概率阈值,便可将一维时序光声测量信号中对应时域区间内的部分归于这一信号类别;最后再将属于同一类别的信号从一维时序光声测量信号中分离出来,便完成了光声测量信号的特征提取;
可选地,本实施例中,对于每一类特征信号设定的概率阈值相同,且均为0.8,相应地,所提取出的各类特征信号如图5所示。
提取出一维时序光声测量信号中的各类特征信号后,可进一步分析得到待测固体样品的膜层厚度等信息。
实施例3:
一种计算机可读存储介质,包括存储的计算机程序;计算机程序被处理器执行时,控制计算机可读存储介质所在设备执行上述实施例1提供的光声测量信号特征的识别模型建立方法,和/或上述实施例2提供的光声测量信号特征的提取方法。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种光声测量信号特征的识别模型建立方法,其特征在于,包括:
分别以不同种类的固体样品为测量对象,获得各测量对象的至少一组一维时序光声测量信号,并标注出每一组一维时序光声测量信号中的每一个特征信号所在的时域区间和对应的类型;所述不同种类的固体样品的膜层数量均为N,N为正整数;不同种类的固体样品的膜层参数不同,膜层参数包括膜层材料和膜层厚度;所述特征信号包括:零点信号和各膜层返回的回声信号;
将U-Net神经网络模型中所有的卷积层、最大池化层和上采样均修改为1D形式,以建立一维U-Net神经网络模型,用于预测目标一维时序光声测量信号中各类特征信号的时域分布概率曲线;所述目标一维时序光声测量信号由具有N个膜层的固体样品在光声测量中产生;
将每一组一维时序光声测量信号及其对应的标注结果作为一个样本,将所有的样本划分为训练集、验证集和测试集,利用所述训练集、所述验证集和所述测试集分别对所述一维U-Net神经网络模型进行训练、验证和测试后,得到所述光声测量信号特征的识别模型。
2.如权利要求1所述的光声测量信号特征的识别模型建立方法,其特征在于,所述一维U-Net神经网络模型中,下采样和上采样部分中卷积层的卷积核均大于或等于3×1。
3.如权利要求1或2所述的光声测量信号特征的识别模型建立方法,其特征在于,建立所述一维U-Net神经网络模型时,还包括:在U-Net神经网络模型的下采样部分和/或上采样部分中,加入至少一个dropout层。
4.如权利要求1所述的光声测量信号特征的识别模型建立方法,其特征在于,在对一维时序光声测量信号进行标注之前,还包括:对一维时序光声测量信号进行数据增强;
所述数据增强包括:添加脉冲噪声,和/或添加低频扰动;
所述低频扰动的频率小于一维时序光声测量信号中特征信号的频率最小值,且所述低频扰动的幅值小于一维时序光声测量信号中特征信号的幅值最小值。
5.如权利要求4所述的光声测量信号特征的识别模型建立方法,其特征在于,所述数据增强还包括:平移变换,和/或裁剪变换。
6.如权利要求4或5所述的光声测量信号特征的识别模型建立方法,其特征在于,在对一维时序光声测量信号进行标注之前,还包括:对一维时序光声测量信号进行预处理;所述预处理包括以下操作中的至少一项:
在对一维时序光声测量信号进行数据增强之前,对一维时序光声测量信号进行数字降噪;
在对一维时序光声测量信号进行数据增强之后,对一维时序光声测量信号进行标准化。
7.如权利要求6所述的光声测量信号特征的识别模型建立方法,其特征在于,所述数字降噪为小波去噪。
8.一种光声测量信号特征的提取方法,其特征在于,包括:
对于具有N个膜层的待测固体样品,获得其一维时序光声测量信号,将所述一维时序光声测量信号输入至光声测量信号特征的识别模型,以预测所述一维时序光声测量信号中各类特征信号的时域分布概率曲线;所述光声测量信号特征的识别模型由权利要求1~7任一项所述的光声测量信号特征的识别模型建立方法得到;
对于每一类特征信号,提取出其时域分布概率曲线中大于相应概率阈值的曲线部分,并将所提出的曲线部分所对应的时域区间识别为该类特征信号所在的时域区间,从所述一维时序光声测量信号中提取出该时域区间内的信号,完成对该类特征信号的提取;
其中,N为正整数。
9.如权利要求8所述的光声测量信号特征的提取方法,其特征在于,所述一维时序光声测量信号输入至光声测量信号特征的识别模型之前,还包括:对所述一维时序光声测量信号进行数字降噪和标准化。
10.一种计算机可读存储介质,其特征在于,包括存储的计算机程序;所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行权利要求1~7任一项所述的光声测量信号特征的识别模型建立方法,和/或权利要求8或9所述的光声测量信号特征的提取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110814894.8A CN113607068B (zh) | 2021-07-19 | 2021-07-19 | 一种光声测量信号特征的识别模型建立方法及提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110814894.8A CN113607068B (zh) | 2021-07-19 | 2021-07-19 | 一种光声测量信号特征的识别模型建立方法及提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113607068A CN113607068A (zh) | 2021-11-05 |
CN113607068B true CN113607068B (zh) | 2022-08-05 |
Family
ID=78337937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110814894.8A Active CN113607068B (zh) | 2021-07-19 | 2021-07-19 | 一种光声测量信号特征的识别模型建立方法及提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113607068B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114548191B (zh) * | 2022-04-27 | 2022-08-05 | 之江实验室 | 光声成像环形稀疏阵列信号预测方法与装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900633A (en) * | 1997-12-15 | 1999-05-04 | On-Line Technologies, Inc | Spectrometric method for analysis of film thickness and composition on a patterned sample |
WO2003006918A2 (en) * | 2001-07-13 | 2003-01-23 | Rudolph Technologies Inc. | Method and apparatus for increasing signal to noise ratio in a photoacoustic film thickness measurement system |
JP2005338063A (ja) * | 2004-04-28 | 2005-12-08 | Japan Science & Technology Agency | 試料の物理的な特性を測定するための装置 |
CN110764064A (zh) * | 2019-11-08 | 2020-02-07 | 哈尔滨工业大学 | 基于深度卷积神经网络集成的雷达干扰信号识别方法 |
CN111242276A (zh) * | 2019-12-27 | 2020-06-05 | 国网山西省电力公司大同供电公司 | 一种用于负荷电流信号识别的一维卷积神经网络构建方法 |
CN112364779A (zh) * | 2020-11-12 | 2021-02-12 | 中国电子科技集团公司第五十四研究所 | 信号处理与深-浅网络多模型融合的水声目标识别方法 |
CN112507881A (zh) * | 2020-12-09 | 2021-03-16 | 山西三友和智慧信息技术股份有限公司 | 一种基于时间卷积神经网络的sEMG信号分类方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6504618B2 (en) * | 2001-03-21 | 2003-01-07 | Rudolph Technologies, Inc. | Method and apparatus for decreasing thermal loading and roughness sensitivity in a photoacoustic film thickness measurement system |
US20060256916A1 (en) * | 2005-05-13 | 2006-11-16 | Rudolph Technologies, Inc. | Combined ultra-fast x-ray and optical system for thin film measurements |
US11832969B2 (en) * | 2016-12-22 | 2023-12-05 | The Johns Hopkins University | Machine learning approach to beamforming |
-
2021
- 2021-07-19 CN CN202110814894.8A patent/CN113607068B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900633A (en) * | 1997-12-15 | 1999-05-04 | On-Line Technologies, Inc | Spectrometric method for analysis of film thickness and composition on a patterned sample |
WO2003006918A2 (en) * | 2001-07-13 | 2003-01-23 | Rudolph Technologies Inc. | Method and apparatus for increasing signal to noise ratio in a photoacoustic film thickness measurement system |
JP2005338063A (ja) * | 2004-04-28 | 2005-12-08 | Japan Science & Technology Agency | 試料の物理的な特性を測定するための装置 |
CN110764064A (zh) * | 2019-11-08 | 2020-02-07 | 哈尔滨工业大学 | 基于深度卷积神经网络集成的雷达干扰信号识别方法 |
CN111242276A (zh) * | 2019-12-27 | 2020-06-05 | 国网山西省电力公司大同供电公司 | 一种用于负荷电流信号识别的一维卷积神经网络构建方法 |
CN112364779A (zh) * | 2020-11-12 | 2021-02-12 | 中国电子科技集团公司第五十四研究所 | 信号处理与深-浅网络多模型融合的水声目标识别方法 |
CN112507881A (zh) * | 2020-12-09 | 2021-03-16 | 山西三友和智慧信息技术股份有限公司 | 一种基于时间卷积神经网络的sEMG信号分类方法及系统 |
Non-Patent Citations (3)
Title |
---|
"U-Net: Convolutional Networks for Biomedical";N. Navab et al.;《Springer International Publishing Switzerland》;20151231;第234–241页 * |
基于U-Net神经网络的肥厚型心肌病与高血压性左心室肥厚磁共振图像定量分析与鉴别;焦梓灵等;《磁共振成像》;20200920;第11卷(第09期);第741-746页 * |
激光超声缺陷统计特征神经网络识别技术研究;郭华玲等;《应用激光》;20171215;第37卷(第06期);第888-892页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113607068A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11449757B2 (en) | Neural network system for non-destructive optical coherence tomography | |
US10176588B2 (en) | System and method for specular surface inspection | |
Kershenbaum et al. | An image processing based paradigm for the extraction of tonal sounds in cetacean communications | |
Liu et al. | Pulsed eddy current data analysis for the characterization of the second-layer discontinuities | |
CN111783616B (zh) | 一种基于数据驱动自学习的无损检测方法 | |
Karthikeyan et al. | Explainable AI-infused ultrasonic inspection for internal defect detection | |
Wang et al. | The aircraft skin crack inspection based on different-source sensors and support vector machines | |
CN113607068B (zh) | 一种光声测量信号特征的识别模型建立方法及提取方法 | |
CN117783051B (zh) | 一种基于多传感器数据融合的甲烷气体泄漏检测方法 | |
CN116524313A (zh) | 基于深度学习和多模态图像的缺陷检测方法及相关装置 | |
CN113962968B (zh) | 面向复杂电磁环境下多源混合干扰雷达图像目标检测系统 | |
CN106682604B (zh) | 一种基于深度学习的模糊图像检测方法 | |
CN114660180A (zh) | 一种基于声发射与1D CNNs的中小桥梁轻量化健康监测方法与系统 | |
CN118151119A (zh) | 一种面向检索任务的毫米波雷达开集步态识别方法 | |
CN112818762B (zh) | 一种大尺寸复合材料及其夹层结构快速无损检测方法 | |
Cantero-Chinchilla et al. | A data-driven approach to suppress artefacts using PCA and autoencoders | |
US20190139214A1 (en) | Interferometric domain neural network system for optical coherence tomography | |
CN112071423A (zh) | 一种基于机器学习的免疫层析浓度检测方法及系统 | |
Jayasudha et al. | Weld defect segmentation and feature extraction from the acquired phased array scan images | |
CN114201993B (zh) | 一种检测超声缺陷的三分支注意力特征融合方法与系统 | |
CN114881938A (zh) | 基于小波分析和神经网络的晶粒度检测方法及系统 | |
Nebaba et al. | Patterns detection in saed images of transmission electron microscopy | |
Hader et al. | On noise-sensitive automatic tuning of gate-defined sensor dots | |
Momtaz Dargahi et al. | NDE data correlation using encoder–decoder networks with wavelet scalogram images | |
Zhang et al. | Nondestructive Testing of Steel Wire Rope Based on Gagnetic Signal and Infrared Information |
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 |