CN113504306A - 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 - Google Patents
一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 Download PDFInfo
- Publication number
- CN113504306A CN113504306A CN202110788227.7A CN202110788227A CN113504306A CN 113504306 A CN113504306 A CN 113504306A CN 202110788227 A CN202110788227 A CN 202110788227A CN 113504306 A CN113504306 A CN 113504306A
- Authority
- CN
- China
- Prior art keywords
- matrix
- low
- rank matrix
- phased array
- recovery
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明属于无损检测处理技术领域,具体涉及一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法。通过差分等声程采样的方式对目标超声相控阵信号进行采样,得到关于目标相控阵的结构性缺失的回波信号矩阵MD;然后对结构性缺失的回波信号矩阵MD进行添零补全,得到低秩矩阵M;通过低秩矩阵恢复ADMN算法对低秩矩阵进行恢复。在本发明中,通过差分等声程采样的方式有效减少超声相控阵回波矩阵的采样数据量,从而提升了回波信号的采集速度。在进行缺陷检测时,也可以将缺陷的回波信息保留,可以进行缺陷报警。本发明在提升检测速率的同时,还保证了成像效果。
Description
技术领域
本发明属于无损检测处理技术领域,具体涉及一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法。
背景技术
无损检测方法在工业质量检测中得到了很好的应用。超声相控阵检测技术是近年来发展起来和广泛应用的一项新兴无损检测技术,是常规超声检测技术的延申,与常规超声的单晶探头不同,超声相控阵探头是一组阵列分布的压电晶片,通过不同的延迟法则激发各个压电晶片,使得声场声束进行偏转聚焦。
现有的基于超声相控阵进行钢轨缺陷检测方法在采集数据时,通常采用全采样的方式获取采集数据。全采样的方式是对一次波束总数据的每一行回波进行相同起点和相同声程的区域进行采样,数据采集量大,导致采集速度低,相比于后数据处理的过程,其已经是局限超声相控阵检测技术检测速率的主要原因。近年来,随着高铁建设的加速,我国的高铁保有量巨大,这就对集成于探伤小车的超声相控阵检测技术的检测速率提出了新的要求。为符合发展需求,目前已经有很多学者提出通过稀疏发射矩阵的方式来减少采集数据量,从而加快检测速率。但是在实际应用中发现,这种方式仅适用于全聚焦成像方法,且该方式会损失本身的成像效果,影响检测结果准确性。
发明内容
本发明的目的在于:提供一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,以解决目前利用超声相控阵进行钢轨缺陷检测时,因检测速率提升后造成的成像效果差,影响检测结果准确性的问题。
为解决上述问题,本发明采取如下技术方案:
一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,包括以下步骤:
步骤1、构建低秩矩阵模型
1.1、对于已知的目标超声相控阵信号采用差分等声程采样,得到结构性缺失的回波信号矩阵MD;
1.2、对结构性缺失的回波信号矩阵MD进行添零补全,即对未采样部分赋值为零,通过赋值补充实现添零补全,得到低秩矩阵M;构建低秩矩阵X且X=M,并根据低秩矩阵X构建差分等声程采样回波信号恢复模型:
s.t.Xij=Mij,(i,j)∈Ω
其中,Ω为低秩矩阵M中超声相控阵回波信号的已知位置,||X||*为矩阵X的核范数(奇异值之和),s.t.Xij表示矩阵X第i行第j列的值,Mij表示矩阵M第i行第j列的值,(i,j)∈Ω表示只对已知位置进行相等赋值
步骤2、利用低秩矩阵恢复算法对低秩矩阵进行恢复,具体过程如下:
2.1、构造增广拉格朗日函数,该函数表示为:
2.2、将步骤1.2构建的差分等声程采样超声回波信号恢复模型输入增广拉格朗日函数中,然后利用交替方向算法ADMM对增广拉格朗日函数中的变量X、E、P进行优化更新,对低秩矩阵M进行恢复。
步骤3、运用超声相控阵扇扫成像对步骤2.2恢复后的低秩矩阵M进行成像;通过对图像分析比对,实现对钢轨内部缺陷情况的分析、判断,完成钢轨内部缺陷检测。
进一步的,所述步骤2.2详细过程为:
根据恢复的精度需求设定初始参数、算法迭代停止阈值∈=1e-5、交替方向乘子法惩罚参数初始值μ0=1e-4、交替方向乘子法惩罚参数最大值ρmax=1e10、交替方向乘子法惩罚参数增长率η=1.1,
利用交替方向算法ADMM循环更新增广拉格朗日函数中的变量X、E、P,直至更新的变量满足设定的算法迭代停止阈值∈,即||X(k+1)-X(k)||∞<∈,||E(k+1)-E(k)||∞<∈时,则结束循环以得到恢复矩阵。其中,X(k+1)和X(k)代表相邻状态的矩阵X,E(k+1)和E(k)代表相邻状态的矩阵E,|| ||∞代表矩阵的无穷范数。
本发明提供的一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,通过差分等声程采样的方式对目标超声相控阵信号进行采样,得到关于目标相控阵的结构性缺失的回波信号矩阵MD;然后对结构性缺失的回波信号矩阵MD进行添零补全,得到低秩矩阵M;通过低秩矩阵恢复ADMN算法对低秩矩阵进行恢复。在本发明中,所述差分等声程采样是指对一次波束总数据的每一行回波,进行交替行起点相同和声程相同的区域进行的采样,但差分等声程采样的相邻行回波起点差值与声程长度相加应与全采样的声程长度相等。差分等声程采样方式相比于全采样减少了起点差值部分的数据,因此能够有效减少超声相控阵回波矩阵的采样数据量,从而提升了回波信号的采集速度。在进行缺陷检测时,也可以将缺陷的回波信息保留,可以进行缺陷报警。
与现有的超声相控阵缺陷检测技术相比,本发明在提升检测速率的同时,还保证了成像效果,使检测的准确性更高。
附图说明
图1是的钢轨缺陷检测方法流程图;
图2(a)是现有技术的超声相控阵回波信号采样的示意图;
图2(b)是超声相控阵回波信号差分等声程采样的示意图;
图2(c)是实施例恢复结构图;
图3(a)是全采样时的成像图;
图3(b)是差分值x=25mm时差分等声程采样低秩矩阵恢复后的成像图;
图3(c)是差分值x=33mm时差分等声程采样低秩矩阵恢复后的成像图;
图3(d)是差分值x=37.5mm时差分等声程采样低秩矩阵恢复后的成像图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
如图1所示,本发明提供的一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,包括以下步骤:
步骤1、构建低秩矩阵模型
1.1、对于已知的目标超声相控阵信号采用差分等声程采样,得到结构性缺失的回波信号矩阵MD。
为了更直观的对比全采样与差分等声程采样的区别,在本实施例中,图2(a)是现有技术的超声相控阵回波信号全采样示意图,全采样是对一次波束总数据的每一行回波进行相同起点和相同声程的区域进行采样。如图2(a)所示,全采样得到的矩阵MF中,超声相控阵的扇扫所有角度,采样的起始点和长度都相同。图2(b)是超声相控阵回波信号差分等声程采样的示意图,差分等声程采样是对一次波束总数据的每一行回波,进行交替行起点相同和回波相同声程的区域进行采样,但差分声程采样的相邻行回波起点差值与声程长度相加应与全采样的声程长度相等;如图2(b)所示,差分等声程采样得到的结构性缺失的回波信号矩阵MD中,一个波束内的相邻回波信号采样起始点相差xmm,但采样声程长度相同。对比图2(a)和图2(b)可以看出,差分等声程采样方式相较于全采样减少了起点差值部分的数据,数据量更少,因此能够减少采样时间。
1.2、对结构性缺失的回波信号矩阵MD进行添零补全,即对未采样部分赋值为零,通过赋值补充实现添零补全,得到低秩矩阵M;构建低秩矩阵X且X=M,并根据低秩矩阵X构建差分等声程采样回波信号恢复模型:
s.t.Xij=Mij,(i,j)∈Ω
其中,Ω为低秩矩阵M中,相超声控阵回波信号的已知位置,||X||*为矩阵X的核范数(奇异值之和),s.t.Xij表示矩阵X第i行第j列的值,Mij表示矩阵M第i行第j列的值,(i,j)∈Ω表示只对已知位置进行相等赋值。
步骤2、利用低秩矩阵恢复算法对低秩矩阵进行恢复,具体过程如下:
2.1、构造增广拉格朗日函数,该函数表示为:
2.2、将步骤1.2构建的差分等声程采样超声回波信号恢复模型输入增广拉格朗日函数中,然后利用交替方向算法ADMM对增广拉格朗日函数中的变量X、E、P进行优化更新,对低秩矩阵M进行恢复。具体过程可以参阅下表:
由上表可以清楚的看出本实施例对低秩矩阵M的恢复过程,该过程具体阐述如下;
2.2.1、根据恢复的精度需求设定初始参数、算法迭代停止阈值∈、交替方向乘子法惩罚参数初始值μ0、交替方向乘子法惩罚参数最大值ρmax、交替方向乘子法惩罚参数增长率η;
2.2.2、利用交替方向算法ADMM循环更新增广拉格朗日函数中的变量X、E、P;变量X更新方式为:
变量E更新方式为:
变量P更新方式为:
Pt+1=Pt+μ(Xt+1-M-Et+1)
2.2.3、直至更新的变量满足设定的算法迭代停止阈值∈,即||X(k+1)-X(k)||∞<∈,||E(k+1)-E(k)||∞<∈时,则结束循环得到恢复矩阵。恢复完成的矩阵如图2(c)所示,图2(c)的上半部分为本实施例差分等声程采样的回波信号,2(c)的下半部分为恢复完成的矩阵,该矩阵的虚线部分为原始信号,点线部分为矩阵恢复的未采样回波信号。
步骤3、运用超声相控阵扇扫成像对恢复后的低秩矩阵M进行成像;通过对图像分析,并将其与标准的钢轨图像比对,实现对钢轨内部缺陷情况的分析、判断,完成钢轨内部缺陷检测。需要注意的是,此处的标准钢轨图像即为未损坏过的图像,这类图像再行业中是可以通过常规手段获得的。
为验证本实施例方案的可行性,本实施例针对同一超声相控阵回波信号(10mm至110mm)分别采用全采样成像和差分采样成像:
图3(a)为全采样图,图中给出的是基于超声相控阵低秩矩阵恢复的钢轨缺陷检测系统对B型便携式相控阵试块超声回波声程10mm至110mm全采样的扇扫成像图;图3(b)是差分值x=25mm时差分等声程采样低秩矩阵恢复后的成像图;图3(c)是差分值x=33mm时差分等声程采样低秩矩阵恢复后的成像图;图3(d)是差分值x=37.5mm时差分等声程采样低秩矩阵恢复后的成像图。图3(b)、图3(c)、图3(d)均能进行有效的结构低秩矩阵恢复,其各自的PSNR值为44.48dB、44.01dB、43.99dB,数值非常接近且恢复效果好。由此可见,本发明所采用的算法对低秩矩阵恢复过程中,当恢复部分存在缺陷时,恢复后得出的理论值能表现出缺陷的特性。本发明的方案是可以用于缺陷检测。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,其特征在于:包括以下步骤:
步骤1、构建低秩矩阵模型
1.1、对于已知的目标超声相控阵信号采用差分等声程采样,得到结构性缺失的回波信号矩阵MD;
1.2、对结构性缺失的回波信号矩阵MD进行添零补全,即对未采样部分赋值为零,通过赋值补充实现添零补全,得到低秩矩阵M;构建低秩矩阵X且X=M,并根据低秩矩阵X构建差分等声程采样回波信号恢复模型:
s.t.Xij=Mij,(i,j)∈Ω
其中,Ω为低秩矩阵M中相超声控阵回波信号的已知位置,||X||*为矩阵X的核范数(奇异值之和),s.t.Xij表示矩阵X第i行第j列的值,Mij表示矩阵M第i行第j列的值,(i,j)∈Ω表示只对已知位置进行相等赋值;
步骤2、利用低秩矩阵恢复算法对低秩矩阵进行恢复,具体过程如下:
2.1、构造增广拉格朗日函数,该函数表示为:
2.2、将步骤1.2构建的差分等声程采样超声回波信号恢复模型输入增广拉格朗日函数中,然后利用交替方向算法ADMM对增广拉格朗日函数中的变量X、E、P进行优化更新,对低秩矩阵M进行恢复。
步骤3、运用超声相控阵扇扫成像对恢复后的低秩矩阵M进行成像;通过对图像分析比对,实现对钢轨内部缺陷情况的分析、判断,完成钢轨内部缺陷检测。
2.根据权利要求1所述的一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法,其特征在于:所述步骤2.2详细过程为:
根据恢复的精度需求设定初始参数、算法迭代停止阈值∈=1e-5、交替方向乘子法惩罚参数初始值μ0=1e-4、交替方向乘子法惩罚参数最大值ρmax=1e10、交替方向乘子法惩罚参数增长率η=1.1,
利用交替方向算法ADMM循环更新增广拉格朗日函数中的变量X、E、P,直至更新的变量满足设定的算法迭代停止阈值∈,即X(k+1)-X(k)||∞<∈,E(k+1)-E(k)||∞<∈时,则结束循环以得到恢复矩阵;其中,X(k+1)和X(k)代表相邻状态的矩阵X,E(k+1)和E(k)代表相邻状态的矩阵E,|| ||∞代表矩阵的无穷范数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110788227.7A CN113504306B (zh) | 2021-07-13 | 2021-07-13 | 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110788227.7A CN113504306B (zh) | 2021-07-13 | 2021-07-13 | 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113504306A true CN113504306A (zh) | 2021-10-15 |
CN113504306B CN113504306B (zh) | 2023-03-28 |
Family
ID=78012845
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110788227.7A Active CN113504306B (zh) | 2021-07-13 | 2021-07-13 | 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113504306B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111768023A (zh) * | 2020-05-11 | 2020-10-13 | 国网冀北电力有限公司电力科学研究院 | 一种基于智慧城市电能表数据的概率峰值负荷估计方法 |
CN115239632A (zh) * | 2022-06-23 | 2022-10-25 | 西南交通大学 | 一种融合巡检图像和超声波图像的钢轨表面伤损检测方法 |
CN115656343A (zh) * | 2022-12-07 | 2023-01-31 | 汕头市超声检测科技有限公司 | 一种基于串列式矩阵扫查的钢轨焊缝缺陷定位方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123702A (zh) * | 2014-07-08 | 2014-10-29 | 浙江传媒学院 | 一种基于交替方向法的高效图像修复方法 |
-
2021
- 2021-07-13 CN CN202110788227.7A patent/CN113504306B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123702A (zh) * | 2014-07-08 | 2014-10-29 | 浙江传媒学院 | 一种基于交替方向法的高效图像修复方法 |
Non-Patent Citations (2)
Title |
---|
仲小伟: "低秩矩阵近似理论与应用", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 * |
王淑琴: "低秩矩阵逼近模型的改进及应用研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111768023A (zh) * | 2020-05-11 | 2020-10-13 | 国网冀北电力有限公司电力科学研究院 | 一种基于智慧城市电能表数据的概率峰值负荷估计方法 |
CN111768023B (zh) * | 2020-05-11 | 2024-04-09 | 国网冀北电力有限公司电力科学研究院 | 一种基于智慧城市电能表数据的概率峰值负荷估计方法 |
CN115239632A (zh) * | 2022-06-23 | 2022-10-25 | 西南交通大学 | 一种融合巡检图像和超声波图像的钢轨表面伤损检测方法 |
CN115239632B (zh) * | 2022-06-23 | 2024-04-09 | 西南交通大学 | 一种融合巡检图像和超声波图像的钢轨表面伤损检测方法 |
CN115656343A (zh) * | 2022-12-07 | 2023-01-31 | 汕头市超声检测科技有限公司 | 一种基于串列式矩阵扫查的钢轨焊缝缺陷定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113504306B (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113504306B (zh) | 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法 | |
EP3382386B1 (en) | Defect detection using ultrasound scan data | |
US12007362B2 (en) | Method for reconstructing crack profiles based on composite-mode total focusing method | |
GB2610449A (en) | Efficient high-resolution non-destructive detecting method based on convolutional neural network | |
CN111175381B (zh) | 基于全矩阵数据的复合构件界面快速成像定量检测方法 | |
JP2008139304A (ja) | 任意の表面輪郭を有する部材の超音波浸漬検査 | |
CN111122700A (zh) | 一种提高激光超声saft缺陷定位速度的方法 | |
CN111239246B (zh) | 一种分步筛选有效信号的曲面结构缺陷全聚焦成像方法 | |
CN106802323A (zh) | 一种基于全矩阵数据的超声全聚焦成像系统 | |
CN112508901A (zh) | 一种水下结构病害识别方法、系统、装置及存储介质 | |
CN111257426A (zh) | 火箭燃料贮箱焊缝的多模式全聚焦检测方法、系统及介质 | |
CN113075297A (zh) | 钛合金相控阵线阵超声检测声场模型构建方法 | |
CN111398426A (zh) | 一种全聚焦相控阵三维超声场的测量及补偿校准方法 | |
CN112684005A (zh) | 基于二维矩阵换能器的全聚焦检测方法 | |
CN115223044A (zh) | 基于深度学习的端到端三维探地雷达目标识别方法及系统 | |
CN115856087B (zh) | 基于纵波一发一收超声相控阵探头的全聚焦成像方法 | |
CN113777166A (zh) | 一种基于超声平面波成像和时间反转算子相结合的高分辨缺陷无损检测方法 | |
CN112581456B (zh) | 一种适用于超声成像测井的图像裂缝实时提取方法 | |
CN113552571B (zh) | 基于psm算法的水下激光致声saft成像方法 | |
CN113887454A (zh) | 基于卷积神经网络点源识别的非接触激光超声检测方法 | |
CN115587291A (zh) | 一种基于裂纹超声散射矩阵的去噪表征方法及系统 | |
CN108693252B (zh) | 三角矩阵聚焦成像的工件探伤检测方法 | |
CA3152727A1 (en) | Method and system for generating merged b-scan | |
CN117420209B (zh) | 基于深度学习的全聚焦相控阵超声快速高分辨率成像方法 | |
CN111047547A (zh) | 一种基于多视图tfm的联合缺陷定量方法 |
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 |