一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置及方法

1.本发明属于光热科学与探测技术、信号分析与特征提取、图像处理与识别技术领域，特别是涉及一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置及方法。该成像装置与方法适用于航空航天、微电子及微纳结构等材料缺陷/损伤的精准无损检测与评价领域。


背景技术：

2.红外热波成像检测技术作为一种新兴无损检测技术，具有非接触、直观、探测面积大等优势，适用于薄金属构件、复合材料构件以及高分子材料的无损检测与质量评价等方面。红外热波成像检测技术采用外部激励源对被检测试件进行主动热激励加载，热流由被检测试件表面向内部渗透，当被检测试件内部存在缺陷时，缺陷与被检测试件材料在其热物理参数上存在明显区别，进而加速或阻碍热流扩散过程。此时缺陷的阻碍作用通过热传导扩散到被检测试件表面，利用红外热像仪可以识别出被检测试件的红外辐射差异，进而实现缺陷的有效识别。热流的传播完全依靠被检测试件内部的温度梯度，由高温区域到低温区域扩散，存在对缺陷的放大检测效果，热流的横向热扩散必然导致缺陷被检测尺寸大于实际尺寸。虽然，热波成像检测方法利用热流的非均匀扩散检测手段可以实现对小尺度缺陷的放大检测，但是缺陷尺寸检测偏差大，因此，基于红外热波成像检测技术的缺陷精准检测面临较大挑战，如果在缺陷检测过程中对热流横向热扩散作用进行有效抑制成为实现缺陷精准检测的关键技术瓶颈。
3.在已有研究中，贵州大学魏嘉呈等提出了一种相似光流法(likeness optical flow,lof)的方法，该方法可以在二维对应的热特征图像中来估计热流梯度扩散场，并利用该方法减少横向热扩散对检测精度的影响。通过仿真与试验证明了lof方法对红外图像序列热流横向热扩撒效应起到了较好的抑制作用，同时该方法对于特定特征提取方法具有较好的特征增强效果，但是，在进行lof迭代运算过程中，初始图像与终止图像选取并未制定严格准则，因此，在未知缺陷尺度时，获得的结果准确度无实在意义。哈尔滨工业大学王飞等对lof方法进行优化，基于梯度变化的horn-schunk算法，在光流基本约束方程的基础上附加了全局平滑假设，实现了方程的完备性，通过特征提取算法提取缺陷特征，并经理论计算获取缺陷深度信息，将深度信息反馈到特征图像迭代模型中，最终实现缺陷对横向热扩散的有效抑制。但是由于horn-schunk算法应用的前提为图像缓慢运动(即热流扩散慢)，因此，缺陷等效直径精度依然与实际缺陷直径存在较大偏差。哈尔滨商业大学卜迟武等(专利号：zl201810891371.1)提出一种热障涂层结构脱粘缺陷脉冲飞秒激光激励红外热波检测装置与方法，在该专利中为实现热障涂层结构脱粘缺陷的精准测量，提出采用飞秒激光对试件进行主动热激励加载，缩短热流横向热扩散时间，实现缺陷尺寸的精准测量，但该方法也带来了缺陷检测信噪比低的问题，对于较小的脱粘缺陷无法识别。首都师范大学与北京维泰凯信新技术有限公司的冯立春等(专利号：zl200510077750.x)提出一种红外热波检测
层析图像的重建方法，在该专利中仅根据热流在深度方向扩散特性理论公式实现缺陷的三维重构，并未对热流的平面扩散放大作用进行考虑，因此，实现的层析成像为仅考虑深度方向的伪层析成像，获取的缺陷等效直径与实际直径偏差较大。综上，目前针对红外热波成像检测方法的研究主要集中在如何对缺陷进行三维定位，由于热流梯度扩散的特殊性，对于缺陷尺寸的精准量化表征存在较大挑战，进而导致在该方面的研究总体较少。
4.为了充分考虑热流梯度扩散特性进而对缺陷尺寸进行准确量化表征，本发明涉及一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置与方法，该发明涉及热波干涉与衍射特性、信号分析与特征提取、光热特征融合、特征图像变化追踪及缺陷识别量化与表征等多方面内容，能够有效抑制热流横向热扩散导致的缺陷检测精度差的问题，实现针对金属材料以及复合材料缺陷平面尺寸的有效识别。


技术实现要素：

5.本发明目的是为了解决现有技术中的问题，提出了一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置及方法。本发明适用于航空航天、微电子及微纳结构等领域复合材料、金属材料等无损检测方面。
6.本发明是通过以下技术方案实现的，本发明提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置，所述成像装置包括红外热像仪、第一bnc控制线、第二bnc控制线、第一以太网线、传导光纤、风冷制冷器、半导体制冷器、激光器、控制线、激光器电源、第三bnc控制线、函数发生器、第一usb通讯线、第四bnc控制线、锁相放大器、第五bnc控制线、数据采集卡、第二以太网线、光学相机、计算机、第二usb通讯线、运动控制器、运动控制线、二维移动台、柱透镜、第一偏振片、工程漫射体、套筒、准直镜、夹具、试件和第二偏振片；
7.所述的计算机设有四个信号端，计算机的其中一个信号输入端通过第一以太网线与红外热像仪的信号输出端连接，计算机的第二个信号输入端通过第二以太网线与光学相机的信号输出端连接，计算机的第三个信号输出端通过第一usb通讯线与函数发生器输入端相连接，计算机的第四个信号输出端通过第二usb通讯线与运动控制器输入端连接，红外热像仪的其中一个信号输入端通过第一bnc控制线与函数发生器输出端连接，红外热像仪的第二信号输出端通过第二bnc控制线与数据采集卡输入端连接，风冷制冷器、半导体制冷器以及激光器通过螺栓机械连接，激光器输出端通过传导光纤与准直镜相连接，柱透镜、第一偏振片、工程漫射体、套筒和准直镜通过螺纹连接在套筒上，激光器电源的信号输出端通过控制线与激光器输入端相连接，激光器电源的信号输入端通过第三bnc控制线与函数发生器信号输出端相连接，函数发生器信号输出端通过第四bnc控制线与锁相放大器信号输入端相连接，数据采集卡输出端通过第五bnc控制线与锁相放大器信号输入端相连接，运动控制器输出端通过运动控制线与二维移动台信号输入端相连接，夹具通过螺纹连接固定于二维移动台之上，夹具用于装夹试件，第二偏振片通过螺纹连接与红外热像仪相连接。
8.本发明提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置的成像方法，所述成像方法具体包括：
9.步骤一：明确待检测试件，将试件放置于二维移动台上；
10.步骤二：开启基于热流扩散跟踪的缺陷精准检测光热融合成像装置；
11.步骤三：开启激光器电源若干分钟后，激光器、风冷制冷器以及半导体制冷器已将
激光器温度控制于预设温度，此时打开激光器；
12.步骤四：开启红外热像仪和光学相机并对试验样件进行实时成像，并采用激光器照射试件，计算机控制运动控制器，进而调整二维移动台以确保激光束照射位置全部处于红外热像仪和光学相机视野；
13.步骤五：计算机控制函数发生器产生脉冲触发信号，该脉冲触发信号一路控制激光器产生固定功率的常脉冲激光，另一路控制红外热像仪按照固定频率采集试件热辐射图像，与此同时，计算机控制二维移动台及光学相机工作，其中二维移动台移动方向垂直于线激光方向，红外热像仪、激光器以及光学相机位置是相对静止的，试件随着二维移动台的某一方向匀速运动，线扫描激光作用于试场的位置固定，扫描完成后进行图像序列重构；
14.步骤六：重构后的图像序列等同于受到脉冲热源激励后的热图像响应，将重构的图像利用离散余弦变化获得其特征响应图像序列；
15.步骤七：经过离散余弦变换可以得到不同频率分量的特征图像序列，此时选取试件缺陷位置与无缺陷位置相位特征差异与不同频率对应关系曲线；
16.步骤八：基于两张特征图像，采用卢卡斯卡纳德方法获取图像沿x,y方向的运动速度v
x
,vy；
17.步骤九：获取最终的特征图像；
18.步骤十：将获取的特征图像与光学相机获取的图像进行配准融合，采用拟合方法获取原红外特征图像像素间隙之间的特征值；
19.步骤十一：基于获取的光热融合成像检测结果进行二值化阈值分割与连通域标记，实现缺陷的精准检测；
20.步骤十二：试验结束后，间隔5分钟后，关闭激光器电源、激光器、函数发生器、数据采集卡、锁相放大器、光学相机以及红外热像仪设备。
21.进一步地，在步骤二中，对计算机、激光器电源、数据采集卡、函数发生器和锁相放大器进行开启。
22.进一步地，在步骤三中，开启激光器电源2分钟，激光器温度控制于15℃。
23.进一步地，在步骤四中，激光器前端安置的偏振片与红外热像仪前的偏振片偏振方向一致，以此实现信号的有效滤波，激光器前端安置的工程漫射体将高斯光源转为均匀光，同时柱透镜将均匀光源转化为线激光光源。
24.进一步地，在步骤五中，取视场中线激光热源与其运动方向垂直的δpixel_x列像素，定义为用于重构的采集像素数，在扫描过程中δpixel_x列像素会跟随热源扫过被检测物体的表面并按照一定频率采集记录，要想使δpixel_x列像素能够记录到扫过的全部范围的热波信号，则需要满足下列关系式：
[0025][0026]
式中，v为二维移动台移动速度，pl为单个像素对应的实际物理位置，f为红外热像仪采集频率；
[0027]
当线激光完全扫描完被检测试件后，利用如下公式进行图像序列重构，
[0028][0029]
式中，t
seqi
为红外热像仪采集到的第i帧图像，t
seq11
为重构后的第一张图像。
[0030]
进一步地，在步骤六中，所述重构图像后的离散余弦变换具体为：
[0031][0032][0033]
式中，cu[t
seq
(tn)]为重构图像后的离散余弦变换，a(u)为正交变换系数，n为重构图像帧数。
[0034]
进一步地，在步骤七中，所述相位特征差异定义为，
[0035]
diph＝|ph
nondefect
{cu[t
seq
(tn)]}-ph
defect
{cu[t
seq
(tn)]}|
ꢀꢀ
(4)
[0036]
式中，diph为相位特征差异，ph
nondefect
为无缺陷位置相位特征均值，ph
defect
为缺陷位置相位特征均值；选取diph峰值位置与1/2位置对应频率的相位特征图像分别为p
peak
与p
peak/2
。
[0037]
进一步地，在步骤八中，获得两个方向的变化速率：
[0038][0039][0040]
进一步地，所述最终获取的特征图像为：
[0041]
ph＝ph
peak/2-ph
peak
·vx-ph
peak
·
vyꢀꢀ
(6)。
[0042]
本发明的有益效果为：
[0043]
(1)本发明公开了基于热流扩散跟踪的缺陷精准检测光热融合成像方法，提出一种基于热波时频特征的离散余弦-卢卡斯卡纳德方法，该方法通过对热波特征的时频域变化过程的准确跟踪，并采用逆向重构的方法构建缺陷的平面尺寸。最终可以实现缺陷尺寸检测误差＜2.5％。
[0044]
(2)本发明提出了一种基于红外热像仪与单目光学相机多源多传感数据融合成像检测方法，将基于红外热波成像获取的特征图像与光学图像进行高精度融合，进而可以提高红外热波成像检测技术的缺陷检测分辨率，最终可以实现缺陷定位精度误差＜0.1％。
[0045]
(3)本发明公开了一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置，采用锁相放大器实时监测热流调制信号与红外热像仪采集信号的实时同步，将同步延迟信息反馈给函数发生器实现信号的精准控制。
附图说明
[0046]
图1是线激光扫描示意图；
[0047]
图2是基于热流扩散跟踪的缺陷精准检测光热融合成像装置；
[0048]
图3是检测结果示意图。
[0049]
图中：1-红外热像仪、2-第一bnc控制线、3-第二bnc控制线、4-第一以太网线、5-传导光纤、6-风冷制冷器、7-半导体制冷器、8-激光器、9-控制线、10-激光器电源、11-第三bnc控制线、12-函数发生器、13-第一usb通讯线、14-第四bnc控制线、15-锁相放大器、16-第五bnc控制线、17-数据采集卡、18-第二以太网线、19-光学相机、20-计算机、21-第二usb通讯线、22-运动控制器、23-运动控制线、24-二维移动台、25-柱透镜、26-第一偏振片、27-工程漫射体、28-套筒、29-准直镜、30-夹具、31-试件、32-第二偏振片。
具体实施方式
[0050]
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0051]
本发明的目的是为了提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像方法，以解决目前红外热波成像检测过程中热流横向热扩散作用造成的缺陷检测精度低的问题，本发明提出基于热波时频特征的离散余弦-卢卡斯卡纳德方法，该方法通过对热波特征的时频域变化过程的准确跟踪，并采用逆向重构的方法构建缺陷的平面平面尺寸。最终可以实现复合材料及金属材料等缺陷尺寸检测误差＜2.5％。
[0052]
本发明提出了基于红外热像仪与单目光学相机多源多传感融合成像检测方法，以解决目前单纯依靠红外热像仪进行检测过程中缺陷位置定位精度较低的问题。目前红外热像仪的像素普遍较低(一般常用的红外热像仪为640
×
512像素)，而光学相机的分辨率一般较高，将基于红外热波成像获取的特征图像与光学图像进行高精度融合，进而可以提高红外热波成像检测技术的缺陷检测分辨率，最终可以实现复合材料及金属材料等缺陷定位精度误差＜0.1％。
[0053]
结合图1-图3，本发明提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置，所述成像装置包括红外热像仪、第一bnc控制线、第二bnc控制线、第一以太网线、传导光纤、风冷制冷器、半导体制冷器、激光器、控制线、激光器电源、第三bnc控制线、函数发生器、第一usb通讯线、第四bnc控制线、锁相放大器、第五bnc控制线、数据采集卡、第二以太网线、光学相机、计算机、第二usb通讯线、运动控制器、运动控制线、二维移动台、柱透镜、第一偏振片、工程漫射体、套筒、准直镜、夹具、试件和第二偏振片；
[0054]
所述的计算机设有四个信号端，计算机的其中一个信号输入端通过第一以太网线与红外热像仪的信号输出端连接，计算机的第二个信号输入端通过第二以太网线与光学相机的信号输出端连接，计算机的第三个信号输出端通过第一usb通讯线与函数发生器输入端相连接，计算机的第四个信号输出端通过第二usb通讯线与运动控制器输入端连接，红外热像仪的其中一个信号输入端通过第一bnc控制线与函数发生器输出端连接，红外热像仪的第二信号输出端通过第二bnc控制线与数据采集卡输入端连接，风冷制冷器、半导体制冷器以及激光器通过螺栓机械连接，激光器输出端通过传导光纤与准直镜相连接，柱透镜、第一偏振片、工程漫射体、套筒和准直镜通过螺纹连接在套筒上，激光器电源的信号输出端通
过控制线与激光器输入端相连接，激光器电源的信号输入端通过第三bnc控制线与函数发生器信号输出端相连接，函数发生器信号输出端通过第四bnc控制线与锁相放大器信号输入端相连接，数据采集卡输出端通过第五bnc控制线与锁相放大器信号输入端相连接，运动控制器输出端通过运动控制线与二维移动台信号输入端相连接，夹具通过螺纹连接固定于二维移动台之上，夹具用于装夹试件，第二偏振片通过螺纹连接与红外热像仪相连接。
[0055]
本发明提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置的成像方法，所述成像方法具体包括：
[0056]
步骤一：明确待检测试件，将试件放置于二维移动台上；
[0057]
步骤二：开启基于热流扩散跟踪的缺陷精准检测光热融合成像装置；在步骤二中，对计算机、激光器电源、数据采集卡、函数发生器和锁相放大器进行开启。
[0058]
步骤三：开启激光器电源若干分钟后，激光器、风冷制冷器以及半导体制冷器已将激光器温度控制于预设温度，此时打开激光器；在步骤三中，开启激光器电源2分钟，激光器温度控制于15℃。
[0059]
步骤四：开启红外热像仪和光学相机并对试验样件进行实时成像，并采用激光器照射试件，计算机控制运动控制器，进而调整二维移动台以确保激光束照射位置全部处于红外热像仪和光学相机视野；在步骤四中，激光器前端安置的偏振片与红外热像仪前的偏振片偏振方向一致，以此实现信号的有效滤波，激光器前端安置的工程漫射体将高斯光源转为均匀光，同时柱透镜将均匀光源转化为线激光光源。
[0060]
步骤五：计算机控制函数发生器产生脉冲触发信号，该脉冲触发信号一路控制激光器产生固定功率的常脉冲激光，另一路控制红外热像仪按照固定频率采集试件热辐射图像，与此同时，计算机控制二维移动台及光学相机工作，其中二维移动台移动方向垂直于线激光方向，红外热像仪、激光器以及光学相机位置是相对静止的，试件随着二维移动台的某一方向匀速运动，线扫描激光作用于试场的位置固定，如图1所示，扫描完成后进行图像序列重构；
[0061]
在步骤五中，取视场中线激光热源与其运动方向垂直的δpixel_x列像素，定义为用于重构的采集像素数，在扫描过程中δpixel_x列像素会跟随热源扫过被检测物体的表面并按照一定频率采集记录，要想使δpixel_x列像素能够记录到扫过的全部范围的热波信号，则需要满足下列关系式：
[0062][0063]
式中，v为二维移动台移动(扫描)速度，pl为单个像素对应的实际物理位置，f为红外热像仪采集频率；
[0064]
当线激光完全扫描完被检测试件后，利用如下公式进行图像序列重构，
[0065][0066]
式中，t
seqi
为红外热像仪采集到的第i帧图像，t
seq11
为重构后的第一张图像。
[0067]
步骤六：重构后的图像序列等同于受到脉冲热源激励后的热图像响应，将重构的
图像利用离散余弦变化获得其特征响应图像序列；
[0068]
在步骤六中，所述重构图像后的离散余弦变换具体为：
[0069][0070][0071]
式中，cu[t
seq
(tn)]为重构图像后的离散余弦变换，a(u)为正交变换系数，n为重构图像帧数。
[0072]
步骤七：经过离散余弦变换可以得到不同频率分量的特征图像序列，此时选取试件缺陷位置与无缺陷位置相位特征差异与不同频率对应关系曲线；
[0073]
在步骤七中，所述相位特征差异定义为，
[0074]
diph＝|ph
nondefect
{cu[t
seq
(tn)]}-ph
defect
{cu[t
seq
(tn)]}|
ꢀꢀ
(4)
[0075]
式中，diph为相位特征差异，ph
nondefect
为无缺陷位置相位特征均值，ph
defect
为缺陷位置相位特征均值；选取diph峰值位置与1/2位置对应频率的相位特征图像分别为p
peak
与p
peak/2
。
[0076]
步骤八：基于两张特征图像，采用卢卡斯卡纳德方法获取图像沿x,y方向的运动速度v
x
,vy；
[0077]
在步骤八中，获得两个方向的变化速率：
[0078][0079][0080]
步骤九：获取最终的特征图像；
[0081]
所述最终获取的特征图像为：
[0082]
ph＝ph
peak/2-ph
peak
·vx-ph
peak
·
vyꢀꢀ
(6)。
[0083]
步骤十：将获取的特征图像与光学相机获取的图像进行配准融合，采用拟合方法获取原红外特征图像像素间隙之间的特征值，实现像素扩容与图像分辨率的提高；
[0084]
步骤十一：基于获取的光热融合成像检测结果进行二值化阈值分割与连通域标记，实现缺陷的精准检测；
[0085]
步骤十二：试验结束后，间隔5分钟后，关闭激光器电源、激光器、函数发生器、数据采集卡、锁相放大器、光学相机以及红外热像仪设备。
[0086]
实施例
[0087]
本发明所述的计算机20设有四个信号端，计算机20的其中一个信号输入端通过第一以太网线4与红外热像仪1的信号输出端连接，计算机20的第二个信号输入端通过第二以太网线18与光学相机19的信号输出端连接，计算机20的第三个信号输出端通过第一usb通讯线13与函数发生器12输入端相连接，计算机20的第四个信号输出端通过第二usb通讯线21与运动控制器22输入端连接，红外热像仪1的其中一个信号输入端通过第一bnc控制线2与函数发生器12输出端连接，红外热像仪1的第二信号输出端通过第二bnc控制线3与数据
采集卡17输入端连接，风冷制冷器6、半导体制冷器7以及激光器8通过螺栓机械连接，激光器8输出端通过传导光纤5与准直镜29相连接，柱透镜25、第一偏振片26、工程漫射体27、套筒28、准直镜29通过螺纹连接在套筒28上，激光器电源10的信号输出端通过控制线9与激光器8输入端相连接，激光器电源10的信号输入端通过第三bnc控制线11与函数发生器12信号输出端相连接，函数发生器12信号输出端通过第四bnc控制线14与锁相放大器15信号输入端相连接，数据采集卡17输出端通过第五bnc控制线16与锁相放大器15信号输入端相连接，运动控制器22输出端通过运动控制线23与二维移动台24信号输入端相连接，夹具30通过螺纹连接固定于二维移动台24之上，夹具30用于装夹试件31，第二偏振片32通过螺纹连接与红外热像仪1相连接。
[0088]
根据图2所述的基于热流扩散跟踪的缺陷精准检测光热融合成像装置，本实施方式中红外热像仪1的型号为flir sc6520，响应波长为3.6～5.2μm，像素尺寸为320
×
256，最大帧频率为100hz，光学相机像素为1920
×
1080像素，采用的试件31为预制平底孔模拟脱粘缺陷的碳纤维复合材料。
[0089]
本发明提出一种基于热流扩散跟踪的缺陷精准检测光热融合成像装置的成像方法，所述成像方法具体包括：
[0090]
步骤一：明确待检测试件31，将试件31放置于二维移动台上；
[0091]
步骤二：开启基于热流扩散跟踪的缺陷精准检测光热融合成像装置，此步骤包括计算机20、激光器电源10、数据采集卡17、函数发生器12、锁相放大器15等设备的开启；
[0092]
步骤三：开启激光器电源10约2分钟后，激光器风冷制冷器6以及半导体制冷器7已将激光器8温度控制于15℃，此时打开激光器8；
[0093]
步骤四：开启红外热像仪1和光学相机19并对试验样件进行实时成像，并采用激光器8照射试件，计算机20控制运动控制器22，进而调整二维移动台24以确保激光束照射位置全部处于红外热像仪1和光学相机19视野，其中激光器前端安置的偏振片26与红外热像仪1前的偏振片32偏振方向一致，以此实现信号的有效滤波，激光器前端安置的工程漫射体27将高斯光源转为均匀光，同时柱透镜25将均匀光源转化为线激光光源。
[0094]
步骤五：计算机20控制函数发生器12产生脉冲触发信号，该脉冲触发信号一路控制激光器8产生固定功率的常脉冲激光，另一路控制红外热像仪1按照固定频率采集试件31热辐射图像，与此同时，计算机20控制控制二维移动台24及光学相机19工作，其中二维移动台24移动方向垂直与线激光方向，红外热像仪1、激光器8以及光学相机19位置是相对静止的，试件31随着二维移动台24的某一方向匀速运动，线扫描激光作用于试场的位置固定。
[0095]
取视场中线激光热源与其运动方向垂直的δpixel_x列像素，定义为用于重构的采集像素数，在扫描过程中δpixel_x列像素会跟随热源扫过被检测物体的表面并按照一定频率采集记录，要想使δpixel_x列像素能够记录到扫过的全部范围的热波信号，则需要满足下列关系式：
[0096][0097]
式中，v为二维移动台移动(扫描)速度，pl为单个像素对应的实际物理位置，f为红外热像仪采集频率。
[0098]
当线激光完全扫描完被检测试件后，利用如下公式进行图像序列重构，
[0099][0100]
式中，t
seqi
为红外热像仪1采集到的第i帧图像，t
seq11
为重构后的第一张图像。
[0101]
步骤六：重构后的图像序列等同于受到脉冲热源激励后的热图像响应，将重构的图像利用离散余弦变化获得其特征响应图像序列，
[0102][0103][0104]
式中，cu[t
seq
(tn)]为重构图像后的离散余弦变换，a(u)为正交变换系数，n为重构图像帧数。
[0105]
步骤七：经过离散余弦变换可以得到不同频率分量的特征图像序列，此时选取试件缺陷位置与无缺陷位置相位特征差异与不同频率对应关系曲线，其中相位特征差异定义为，
[0106]
diph＝|ph
nondefect
{cu[t
seq
(tn)]}-ph
defect
{cu[t
seq
(tn)]}|
ꢀꢀ
(10)
[0107]
式中，diph为相位特征差异，ph
nondefect
为无缺陷位置相位特征均值，ph
defect
为缺陷位置相位特征均值。选取diph峰值位置与1/2位置对应频率的相位特征图像分别为p
peak
与p
peak/2
。
[0108]
步骤八：基于两张特征图像，采用卢卡斯卡纳德方法获取图像沿x,y方向的运动速度v
x
,vy,其中，利用该方法可以描述两张特征图像的变换历程，进而获得两个方向的变化速率，
[0109][0110][0111]
步骤九：最终获取的特征图像为，
[0112]
ph＝ph
peak/2-ph
peak
·vx-ph
peak
·
vyꢀꢀ
(12)
[0113]
步骤十：将获取的特征图像与光学相机19获取的图像进行配准融合，采用拟合方法获取原红外特征图像像素间隙之间的特征值，实现像素扩容与图像分辨率的提高。
[0114]
步骤十一：基于获取的光热融合成像检测结果进行二值化阈值分割与连通域标记，实现缺陷的精准检测。
[0115]
步骤十二：试验结束后，间隔5分钟后，关闭激光器电源10、激光器8、函数发生器12、数据采集卡17、锁相放大器15、光学相机19以及红外热像仪1等设备。
[0116]
对检测试件的检测结果如图3所示，其中针对三个模拟平底孔缺陷直径检测误差小于2％。
[0117]
以上对本发明所提出的一种基于热流扩散跟踪的缺陷精准检测光热融合成像装
置及方法进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。