CN113820016B - 一种相位调制热波信号全变分去噪方法 - Google Patents
一种相位调制热波信号全变分去噪方法 Download PDFInfo
- Publication number
- CN113820016B CN113820016B CN202110953030.4A CN202110953030A CN113820016B CN 113820016 B CN113820016 B CN 113820016B CN 202110953030 A CN202110953030 A CN 202110953030A CN 113820016 B CN113820016 B CN 113820016B
- Authority
- CN
- China
- Prior art keywords
- thermal
- noise
- thermal wave
- wave signal
- phase modulation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 230000005284 excitation Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 239000011158 industrial composite Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000001931 thermography Methods 0.000 description 7
- 230000006835 compression Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011208 reinforced composite material Substances 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 229920000049 Carbon (fiber) Polymers 0.000 description 2
- 241001494479 Pecora Species 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000004917 carbon fiber Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004918 carbon fiber reinforced polymer Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000003685 thermal hair damage Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0003—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiant heat transfer of samples, e.g. emittance meter
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明提供了一种相位调制热波信号全变分去噪方法,该方法利用稀疏优化原理对相位调制热波信号进行降噪。将红外热像仪测得的相位调制热波信号建模为可稀疏或稀疏导数表示的热波分量和噪声信号,并可根据1‑范数定义转化为成本函数不可微的无约束最小化问题,并利用最大‑最小迭代优化算法进行求解,以从含噪声热波信号中重构可稀疏或稀疏导数表示的热波分量,进而实现对相位调制热波信号进行降噪的目的。此外,由于所引入的最大‑最小迭代算法利用了带状系统快速求解的优势,可显著提高计算效率。该发明不仅适用于单一的相位调制信号,如巴克码波形,还适用于更为复杂的同时进行频率和相位调制的信号,因此具有广泛应用范围。
Description
技术领域
本发明涉及多物理场光热无损检测技术领域,尤其涉及一种相位调制热波信号全变分去噪方法。
背景技术
脉冲压缩热成像技术近年来由于具有高信噪比、大动态探测深度,已被用于碳纤维增强复合材料、牙齿等诸多领域的无损检测,尤其是该技术即使在仅使用低功率外部激励源情形下,也能显著提高信噪比和增大热成像的探测范围/深度分辨率,避免了对待测样品表面产生热损伤。距离/深度分辨率是制约脉冲压缩热成像技术快速发展的主要瓶颈,除了对待测样品施加合适激励波形外,还可采用先进后处理去噪算法来提高信噪比,提升脉冲压缩热成像的距离分辨率。
巴克码相位调制信号是目前脉冲压缩热成像技术中较为常用的激励波形之一,虽然具有良好的抗噪性能,然而,由于载波频率较低且单一,不利于实现待测样品的一次性动态探测,为了弥补单一巴克码相位调制信号的不足,可同时进行频率和相位调制的正交相位调制信号最近也被提出。然而,干扰噪声的存在往往也不利于利用相位调制类信号实现高分辨率脉冲压缩热成像。
发明内容
本发明目的在于针对现有技术的缺陷,提供一种相位调制热波信号全变分去噪方法,该方法利用稀疏优化原理对相位调制热波信号进行降噪。
为解决上述技术问题,本发明提供技术方案如下:
一种相位调制热波信号全变分去噪方法,其特征在于:包括以下步骤:
步骤10:以低载波频率的巴克码相位调制或正交相位编码线性调频信号作为激励波形s(t)施加到被测样品表面,利用红外热像仪获取样品表面的含噪声热回波信号T(t),所述含噪声热回波信号T(t)通常包括与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号x(t)和噪声w,其中t表示时间变量;
步骤20:根据所施加的所述激励波形s(t),通过网格搜索法确定正则化参数λ和最大迭代次数Nmax,同时定义大小为N×(N-1)的一阶差分矩阵D;
步骤30:通过最优化方程从所述含噪声的热回波信号T(t)来重构所述低载波频率热波信号x(t),所述最优化方程为:
式中,‖·‖1和‖·‖2分别表示1-范数和2-范数,正则化参数λ大于0。
步骤40:通过采用最大-最小迭代优化算法求解式(1)并获得不含噪声的x(t),从而实现所述含噪声热回波信号T(t)的降噪目的。
进一步的,所述含噪声热回波信号T(t)为:
T(t)=x(t)+w
式中,x(t)表示为与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号、w表示噪声,其中t表示时间变量。
进一步的,所述步骤40中的最大-最小迭代优化算法为:
进一步的,所述步骤10中的激励波形s(t)为正交相位编码线性调频激励波形,其表达式为:
式中,foPCLFM(t)=f1+(f2-f1)/T·t表示瞬时频率,C表示13位二进制相位编码,C=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1],每节长度为j/fc,fc为载波频率,n为单个子脉冲的数目,P表示码元系数,
进一步的,所述被测样品为工业复合材料或生物组织。
与现有技术相比,本发明的有益效果是:1、将红外热像仪测得的相位调制热回波信号建模为可稀疏或稀疏导数表示的热波分量和噪声信号,并可根据1-范数定义转化为成本函数不可微的无约束最小化问题,并利用最大-最小迭代算法进行求解,以从含噪声热波信号中重构可稀疏或稀疏导数表示的热波分量,进而实现对相位调制热波信号进行降噪的目的。2、由于所引入的最大-最小迭代优化算法利用了带状系统快速求解的优势,可显著提高计算效率。3、该发明不仅适用于单一的相位调制信号,如巴克码波形,还适用于更为复杂的同时进行频率和相位调制的信号,因此具有广泛应用范围。4、本发明对碳纤维增强复合材料等工业复合材料、牙齿和羊骨等生物组织的高分辨率脉冲压缩热成像具有重要价值。
附图说明
图1为本发明流程示意图;
图2为本发明实施例一的巴克码相位调制激励波形及不含噪声的热回波信号波形图;
图3为本发明实施利一的含噪声的巴克码相位调制热回波信号及全变分去噪重建后的热回波信号波形图;
图4为本发明实施例二的正交相位编码线性调频激励波形及不含噪声热回波信号波形图;
图5为本发明实施例二的含噪声的正交相位编码线性调频热回波信号及全变分去噪重建后的热回波信号波形图。
具体实施方式
为了加深本发明的理解,下面我们将结合附图对本发明作进一步详述,该实施例仅用于解释本发明,并不构成对本发明保护范围的限定。
如图1所示,一种相位调制热波信号全变分去噪方法,其特征在于:包括以下步骤:
步骤10:以低载波频率的巴克码相位调制或正交相位编码线性调频信号作为激励波形s(t)施加到被测样品表面,利用红外热像仪获取样品表面的含噪声热回波信号T(t),含噪声热回波信号T(t)通常包括与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号x(t)和噪声w,其中t表示时间变量,含噪声热回波信号T(t)的表达式为:
T(t)=x(t)+w;
步骤20:根据所施加的激励波形s(t),通过网格搜索法确定正则化参数λ和最大迭代次数Nmax,同时定义大小为N×(N-1)的一阶差分矩阵D;
步骤30:通过最优化方程从含噪声的热回波信号T(t)来重构所述低载波频率热波信号x(t),所述最优化方程为:
式中,‖·‖1和‖·‖2分别表示1-范数和2-范数,正则化参数λ大于0。
步骤40:通过采用最大-最小迭代优化算法求解式(1)并获得不含噪声的x(t),从而实现所述含噪声热回波信号T(t)的降噪目的。
优选的,步骤40中的最大-最小迭代优化算法为:
下面结合具体的实施例进一步说明本发明的方法:
以下实施例以单层厚度为0.2mm、总厚度为3.2mm的碳纤维增强聚合板为被测样品,其热扩散率为α=6.16×10-7m2/s,热传导率为k=0.46W/mK。
实施例一:
步骤10:向被测样品表面施加低载波频率巴克码相位调制激励波形。
该激励波形如图2中的实线所示,所获得的理想状态下的不含噪声的热回波信号波形如图2中虚线所示,含有噪声的热回波信号如图3中灰色线所示。
步骤20:根据所施加的相位调制激励波形,通过网格搜索法确定正则化参数λ为20、最大迭代次数Nmax为300,同时定义大小为N×(N-1)的一阶差分矩阵D。
步骤30:通过以下最优化方程
从含噪声的热回波信号T(t)中重构与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号x(t)。
步骤40:通过采用最大-最小迭代优化算法求解上述最优化方程并获得x(t),如图3中的黑色实线所示,从而实现相位调制热波信号降噪目的。
由表1中的对比可知,在不同的噪声水平下,采用全变分去噪方法处理后的巴克码相位调制热回波波形的均方根误差值均低于传统低通滤波方法,从图2和图3中可以看到经过全变分去噪后的波形几乎能够完美地重构理想状态下的巴克码相位调制热回波信号,实现很好的去噪效果。
实施例二:
步骤10:向被测样品表面施加正交相位编码线性调频激励波形,该波形的起始频率f1=0.1×(1-0.33)Hz,终止频率f2=0.1×(1+0.33)Hz,其相位调制部分是采用具有180°相移的13位巴克码,激励持续时间T=130s,设计出可同时进行频率和相位调制的正交相位编码线性调频激励波形,其表达式为:
式中,fOPCLFM(t)=f1+(f2-f1)/T·t表示瞬时频率,C表示13位二进制相位编码,C=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1],每节长度为j/fc,fc为载波频率,n为单个子脉冲的数目,P表示码元系数,
该激励波形如图4中实线所示,所获得的理想状态下不含噪声的热回波信号波形如图4中虚线所示,含有噪声的热回波信号如图5中灰色线所示。
步骤20:根据所施加的相位调制激励波形,通过网格搜索法确定正则化参数λ为20和最大迭代次数Nmax为300,同时定义大小为N×(N-1)的一阶差分矩阵D。
步骤30:通过以下最优化方程
从含噪声的热回波信号T(t)中重构与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号x(t)。
步骤40:通过采用最大-最小迭代优化算法求解上述最优化方程并获得x(t),如图5的黑色实线所示,从而实现正交相位编码线性调频热波信号降噪目的。
由表1可知,相比于传统低通滤波方法,全变分去噪方法几乎仍能完美地重构正交相位编码线性调频热波信号,实现很好的去噪效果。由表1中的对比可知,在不同的噪声水平下,采用全变分去噪方法处理后的正交相位编码线性调频热回波波形的均方根误差值均低于传统低通滤波方法,从图4和图5中可以看到经过全变分去噪后的波形几乎能够完美地重构理想状态下的正交相位编码线性调频热回波信号,实现很好的去噪效果。
表1不同噪声水平下全变分去噪与传统低通滤波后的RMSE值对比
上述具体实施方式,仅为说明本发明的技术构思和结构特征,目的在于让熟悉此项技术的相关人士能够据以实施,但以上内容并不限制本发明的保护范围,凡是依据本发明的精神实质所作的任何等效变化或修饰,均应落入本发明的保护范围之内。
Claims (5)
1.一种相位调制热波信号全变分去噪方法,其特征在于:包括以下步骤:
步骤10:以低载波频率的巴克码相位调制或正交相位编码线性调频信号作为激励波形s(t)施加到被测样品表面,利用红外热像仪获取样品表面的含噪声热回波信号T(t),所述含噪声热回波信号T(t)包括与激励波形相似但已失真,且可稀疏或稀疏导数表示的低载波频率热波信号x(t)和噪声w,其中t表示时间变量;
步骤20:根据所施加的所述激励波形s(t),通过网格搜索法确定正则化参数λ和最大迭代次数Nmax,同时定义大小为N×(N-1)的一阶差分矩阵D;
步骤30:通过最优化方程从所述含噪声热回波信号T(t)来重构所述低载波频率热波信号x(t),所述最优化方程为:
式中,‖·‖1和‖·‖2分别表示1-范数和2-范数,正则化参数λ大于0;
步骤40:通过采用最大-最小迭代优化算法求解式(1)并获得不含噪声的x(t),从而实现所述含噪声热回波信号T(t)的降噪目的。
2.根据权利要求1所述一种相位调制热波信号全变分去噪方法,其特征在于:所述含噪声热回波信号T(t)为:T(t)=x(t)+w。
5.根据权利要求4所述一种相位调制热波信号全变分去噪方法,其特征在于:所述被测样品为工业复合材料或生物组织。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110953030.4A CN113820016B (zh) | 2021-08-19 | 2021-08-19 | 一种相位调制热波信号全变分去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110953030.4A CN113820016B (zh) | 2021-08-19 | 2021-08-19 | 一种相位调制热波信号全变分去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113820016A CN113820016A (zh) | 2021-12-21 |
CN113820016B true CN113820016B (zh) | 2022-11-11 |
Family
ID=78913285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110953030.4A Active CN113820016B (zh) | 2021-08-19 | 2021-08-19 | 一种相位调制热波信号全变分去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113820016B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114858285A (zh) * | 2022-04-02 | 2022-08-05 | 电子科技大学 | 一种线性调频红外无损检测系统 |
CN115508407A (zh) * | 2022-08-29 | 2022-12-23 | 电子科技大学 | 一种巴克码调制激光激励光声压电热扩散率测量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5745837A (en) * | 1995-08-25 | 1998-04-28 | Terayon Corporation | Apparatus and method for digital data transmission over a CATV system using an ATM transport protocol and SCDMA |
WO2009052663A1 (fr) * | 2007-10-23 | 2009-04-30 | Jianhua Luo | Procédé de débruitage de signal par reconstruction de signal remplaçant des données de spectre |
CN102033081A (zh) * | 2010-10-15 | 2011-04-27 | 哈尔滨工业大学 | 基于图像序列处理的红外锁相热波无损检测方法 |
CN102818786A (zh) * | 2012-08-17 | 2012-12-12 | 中国科学院上海光学精密机械研究所 | 正弦相位调制的并行复频域光学相干层析成像系统和方法 |
CN103345621A (zh) * | 2013-07-09 | 2013-10-09 | 东南大学 | 一种基于稀疏浓度指数的人脸分类方法 |
CN103854263A (zh) * | 2014-03-27 | 2014-06-11 | 湖州师范学院 | 一种基于正交约束投影解混模型的消噪解混算法 |
CN104036456A (zh) * | 2013-03-08 | 2014-09-10 | 株式会社日立制作所 | 图像去噪装置和图像去噪方法 |
-
2021
- 2021-08-19 CN CN202110953030.4A patent/CN113820016B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5745837A (en) * | 1995-08-25 | 1998-04-28 | Terayon Corporation | Apparatus and method for digital data transmission over a CATV system using an ATM transport protocol and SCDMA |
WO2009052663A1 (fr) * | 2007-10-23 | 2009-04-30 | Jianhua Luo | Procédé de débruitage de signal par reconstruction de signal remplaçant des données de spectre |
CN102033081A (zh) * | 2010-10-15 | 2011-04-27 | 哈尔滨工业大学 | 基于图像序列处理的红外锁相热波无损检测方法 |
CN102818786A (zh) * | 2012-08-17 | 2012-12-12 | 中国科学院上海光学精密机械研究所 | 正弦相位调制的并行复频域光学相干层析成像系统和方法 |
CN104036456A (zh) * | 2013-03-08 | 2014-09-10 | 株式会社日立制作所 | 图像去噪装置和图像去噪方法 |
CN103345621A (zh) * | 2013-07-09 | 2013-10-09 | 东南大学 | 一种基于稀疏浓度指数的人脸分类方法 |
CN103854263A (zh) * | 2014-03-27 | 2014-06-11 | 湖州师范学院 | 一种基于正交约束投影解混模型的消噪解混算法 |
Also Published As
Publication number | Publication date |
---|---|
CN113820016A (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113820016B (zh) | 一种相位调制热波信号全变分去噪方法 | |
CN204495996U (zh) | 宽带声纳接收器 | |
CN105281779B (zh) | 多观测值向量稀疏度自适应压缩采样匹配追踪方法 | |
Bi et al. | Multifrequency compressed sensing for 2-D near-field synthetic aperture radar image reconstruction | |
CN103512960B (zh) | 一种超声阵列成像方法 | |
Duan et al. | Relevant modes selection method based on Spearman correlation coefficient for laser signal denoising using empirical mode decomposition | |
CN106772270A (zh) | 一种雷达回波信号的采样方法及重构方法 | |
CN111650654A (zh) | 联合emd与wt算法的地面磁共振信号尖峰噪声剔除方法 | |
CN111934690A (zh) | 基于自适应压缩感知的核信号重构新方法 | |
CN112731306A (zh) | 基于CS和简化FrFT的UWB-LFM信号参数估计方法 | |
CN109586728B (zh) | 基于稀疏贝叶斯的调制宽带转换器框架下信号盲重构方法 | |
Zhang et al. | A novel ECG signal denoising method based on Hilbert-Huang Transform | |
CN110032968A (zh) | 基于双树复小波和自适应半软阈值法的去噪方法 | |
CN113589253A (zh) | 一种基于伪时域的小波变换算法对微弱回波信号进行检测的方法 | |
Peng et al. | MR image reconstruction with convolutional characteristic constraint (CoCCo) | |
Wang et al. | Orthogonal sparse dictionary based on Chirp echo for ultrasound imaging | |
CN113820310B (zh) | 一种频率调制热波信号稀疏辅助去噪方法 | |
CN116383605A (zh) | 基于小波分解的车辆震动信号去噪方法 | |
CN112965066B (zh) | 一种单比特带通采样合成孔径雷达成像方法 | |
CN113820360B (zh) | 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 | |
CN104883192A (zh) | 基于混合压缩采样的非合作超宽带信号的重构方法 | |
Chen et al. | Study of threshold setting for rapid detection of multicomponent LFM signals based on the fourth-order origin moment of fractional spectrum | |
Faisal et al. | Suppression of false-terms in wigner-ville distribution using time and frequency windowing | |
CN108535542B (zh) | 一种寻峰鉴相方法 | |
Meng et al. | A method based on random demodulator and Waveform matching dictionary to estimate LFM signal parameter |
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 |