Study of OCT weld depth curve fitting method for laser keyhole welding
-
摘要: 为了解决激光焊接熔深在线检测及光学相干光层析成像(OCT)熔深检测信号噪声多的问题，采用了一种百分位滤波加移动平均算法的熔深曲线拟合方法。通过光学相干光层析成像系统对304不锈钢匙孔焊过程进行了熔深在线检测，在使用百分位滤波从原始散点数据中提取出百分位熔深曲线过程中，发现该方法提取的百分位熔深曲线存在尖峰干扰噪声，通过移动平均算法对百分位熔深曲线进行低通滤波处理，最终从原始散点数据中提取拟合出OCT熔深测量曲线。结果表明，将OCT熔深测量曲线与焊缝纵抛面的实际熔深曲线进行对比，该方法获取OCT熔深测量曲线的测量精度最大提升了15%。该方法能从OCT原始散点数据中有效提取拟合出更为准确的OCT熔深测量曲线。Abstract: In order to solve the online detection of laser welding depth of fusion and to solve the problem of optical coherence tomography (OCT) depth of fusion detection signal noise, a percentile filtering plus moving average algorithm was used to fit the weld depth curve. The weld depth of 304 stainless steel keyhole welding process was detected online by optical coherence photomicrographic imaging system, in the process of extracting the percentile depth curves from the original scatter data using percentile filtering; it was found that there were spikes interfering noises in the percentile depth melting curves extracted by this method, and the percentile depth curves were processed by low-pass filtering through the moving average algorithm, and the OCT weld depth curves were finally extracted and fitted from the original scatter data. The results show that, comparing the OCT weld depth curve with the actual weld depth curve on the longitudinal surface of the weld, the measurement accuracy of the OCT depth of fusion curve obtained by this method has been improved by 15% at most. The method is effective in extracting and fitting more accurate OCT weld depth curves from the original scatter data.
-
引言
食源性疾病是食品安全领域持续关注的重要问题[1]。食品生产加工设备在工作过程中会产生碎屑残留，这些残留不仅会影响后续食品生产的品质而且容易滋生食源性病菌[2]。因此对食品加工、储存器具上的食品残留物进行检测，不仅能够评判清洁程度，更可以防止病菌滋生的潜在风险。目前国内对于食品生产加工场所的卫生检测仍以目视检查为主。目视检查易受环境和检测员疲劳度的干扰，检测效果不稳定；而且对于一些透明或与背景颜色相近的残留物，目视检查难以取得较好的检测效果。荧光成像技术在食品安全检测领域具有很大的发展潜力，许多食品材料都具有荧光发射特性[3-6]。近年来，一些基于荧光成像技术的设备和系统已经被开发并应用于食品安全检查。KANG等人[7]基于食品样品的特征荧光，设计了一种快速筛查食品中多种危害因子的便携式食品安全检测仪。ZHANG[8]基于感染黄曲霉毒素的花生在365 nm紫外光照射下会发出独特的黄绿色荧光特性，设计了一款基于智能手机的黄曲霉毒素快速检测设备。在肉制品粪便污染检测应用中，SEO等人[9]利用多光谱荧光成像技术结合选定波段比(波段比为630 nm/600 nm)图像的主成分分析技术实现了鸡胴体粪便污染的在线检测；BURFOOT等人 [10]开发一种荧光成像设备用于牛羊肉胴体粪便污染检测；GORJI和SUEKER等人[11-12]将荧光成像和深度学习相结合，实现肉类胴体上粪便污染的自动识别。在食品加工场所卫生检测应用中，LEFCOURT等人[13-14]开发了一款便携式高光谱荧光成像系统，以检测加工设备表面肉眼不可见的残渣；LEE等人[15-17]开发了一款手持式荧光成像设备并应用于加工设备表面食品残留检测。但是上述设备在操作的便携性、功能的完善性方面仍有待改进，难以适应复杂环境下的食品残留物检测。
本文作者研究开发了一套手持式食品残留物荧光成像检测设备，该设备集检测-分析-显示-存储于一身，体积较小方便携带，在成本上更具优势。设备由发光二极管(light-emitting diode, LED)、散热器、超声波测距芯片、红绿蓝(red green blue, RGB)相机、触控屏、现场可编程门阵列(field-programmable gate array, FPGA)、安全数字(secure digital, SD)存储卡、大容量锂电池组成。设备充分利用FPGA实时性高、容量大以及能够进行并行流水线处理的特点[18]，在FPGA中实现基于分块阈值的大津算法[19-20](Otsu algorithm, OTSU)，提高了设备在不均匀光照下的检测效果。检测人员可以使用该设备识别食品加工设备中肉眼看不见的残留物并可以通过触控屏将当前图片保存在SD存储卡中。
1. 系统总体方案设计
本文系统装配如图 1所示，包括主控模块、相机模块、激发光源模块、超声波测距模块、触控屏、电源板。设备启动时，打开电源开关，大功率LED激发光源启动照射，检测区域的残留物荧光图像经RGB相机采集后传送到FPGA主控模块中，FPGA根据距离信息对荧光图像进行处理，标记出残留物区域并发送到触控屏中显示，此时检测员可以根据需要选择是否操作触控屏将此时的荧光图像保存至SD存储卡中。
![图 1 手持式食品残留物荧光检测系统装配图]()
图 1 手持式食品残留物荧光检测系统装配图Figure 1. Assembly drawing of a hand-held fluorescence detection system for food residues2. 系统硬件设计
2.1 FPGA主控板设计
如图 2所示，所设计的FPGA主控板集成了EP4CE22F17C8主控芯片、50 MHz高精度有源晶振时钟、256 Mbit随机动态存取内存(synchronous dynamic random access memory, SDRAM)、和16 Mbit闪存(flash eeprom memory, FLASH)。根据FPGA现场可编程的优势，所设计的FPGA主控板在FPGA芯片内部实现了串行摄像机控制总线协议(serial camera control bus, SCCB)、数字视频接口(digital video port, DVP)与RGB相机、串行外设接口(serial peripheral interface, SPI)与SD存储卡，以及集成电路总线(inter-integrated circuit, IIC)、视频图形阵列(video graphics array, VGA)与触控屏的通信与数据传输。
2.1.1 SDRAM存储电路设计
系统运行过程中需要开辟一块内存用来缓存图像处理过程中产生的数据，因此使用256 Mbit SDRAM作为图像检测过程中的缓存存储器。SDRAM存储器硬件电路连接如图 3所示，SDRAM与FPGA相连，由于相机图像输出频率为42 MHz，RGB触控屏显示频率为32 MHz，SDRAM工作频率为133 MHz，因此FPGA在缓存和读取SDRAM过程中需要分别插入两个异步先进先出队列防止数据丢失。此外由于SDRAM工作频率较高，在进行印刷电路板布线时需要进行蛇形等长处理。
2.1.2 FLASH电路设计
由于FPGA是基于随机存取存储器和查找表结构，掉电时配置信息会丢失，因此选用一片16 Mbit SPI FLASH保存配置数据。所用FLASH型号为W25Q16，该型号可以等效替代相关公司的系列配置器件。电路连接如图 4所示，程序下载时，下载数据经过DATA0口保存在FLASH芯片中，当系统再次上电时FPGA将读取FLASH中保存的用户程序和配置文件完成系统初始化配置。
2.1.3 SD卡存储电路设计
SD卡存储电路使用一张大容量SD卡作为外部存储器。SD卡存储电路如图 5所示，SD卡与FPGA采用SPI通信模式，FPGA可以根据指令将当前视频帧保存在SD卡中。
2.1.4 电源电路设计
系统主要用到12 V、5 V、3.3 V、2.5 V、1.2 V直流电源。其中大功率LED光源以及散热风扇由12 V锂电池直接供电，RGB触控屏和FPGA主板电源由12 V锂电池经电压转换芯片MC34063转换为5 V后提供，FPGA主板电源经电压转换芯片AMS1085-3.3转换为3.3 V电压给FPGA芯片供电，并分别经过电压转换芯片AMS1117-2.5和AMS1117-1.2转换为2.5 V和1.2 V电压, 作为联合测试工作组和FPGA内核供电电压。
2.2 激发光源模块
激发光源由4颗10 W LED灯珠、恒流电路、铝制散热片、光源基板组成。根据食品残留物的荧光特性，选用中心波长400 nm，波长范围395 nm~405 nm的大功率LED作为激发光源。为保证LED稳定工作，灯珠需要恒定的电流和及时散热。如图 6所示，灯珠按照两并两串方式均匀分布在光源板上，且每颗LED灯珠均搭配相应的铝制散热片，中部镂空放置相机镜头。
2.3 RGB相机模块
选用高灵敏度、低串扰、低噪声、高量子效率OV5640相机采集荧光图像，并按照RGB565格式发送给FPGA。相机镜头选用广角高清M12镜头。镜头底部具有放置滤波片的凹槽，设备可以通过更换凹槽中的滤波片实现对不同荧光特性对象的检测。
2.4 超声波测距模块
为避免误检测，设备通过超声波测距模块实现检测状态判断。如图 7所示，超声波测距模块采用高性价比超声波测距芯片CS100。
2.5 RGB触控屏模块
RGB触控屏采用800×480像素分辨率的电容触摸显示屏。该屏幕集成了电容触摸控制器，该控制器可以通过标准IIC总线接口将多点触摸信号的实时坐标传送给FPGA。
2.6 系统板级信号测试
如图 8所示，采用嵌入式逻辑分析仪SingalTap抓取系统运行过程中的时序关系。CMOS_Frame_clken为像素有效信号，CMOS_Vsync、CMOS_Href和COMS_Data分别为相机输出的场信号、行信号和像素数据，当CMOS_Href为高电平，且CMOS_Frame_clken使能时像素数据有效输出；SDRAM_Read和SDRAM_Write分别为SDRAM存储器读写控制信号，当SDRAM_Read或SDRAM_Write为高电平时，FPGA根据读写地址SDRAM_Rrd1_addr或SDRAM_Rwrd1_addr对SDRAM存储器进行读写，读写数据通过双向读写数据通道SDRAM_DQ进行传输，SDRAM_Dqm为2位数据掩码器，SDRAM_Dqm为16进制数0 x 03时, 表示当前传输数据SDRAM_DQ无效，为16进制数0 x 00时, 表示数据有效；LCD_Vs、LCD_Hs和LCD_De分别为显示器场信号、行信号和数据使能信号，当LCD_Vs、LCD_Hs和LCD_De均为高电平时，像素数据LCD_Data有效输出；SDCARD_Init_done信号为高电平表明设备检测到SD存储卡并初始化配置完成；Ultrasonic_Echo为超声波测距模块输出回响信号，其高电平脉冲宽度为超声波往返时间之和，FPGA根据回响信号高电平持续时间计算距离，当回响信号Ultrasonic_Echo产生下降沿时，FPGA内部寄存器Ultrasonic_distance更新距离信息(单位为mm)。
由硬件测试信号可知，硬件系统图像采集功能正常、SDRAM存储器读写正常、RGB显示器工作正常、SD卡初始化完成、超声波模块测距正常，系统硬件电路设计达到设计要求。
3. FPGA逻辑设计
3.1 系统流程设计
系统工作流程如图 9所示。设备启动时，大功率LED灯打开同时散热风扇开始工作，FPGA开始初始化配置并将掩模寄存器初始值清零；FPGA初始化完成后开始初始化配置相机寄存器；相机初始化完成后系统开始采集图像信息，为确保图像输出稳定系统丢弃前10帧图像；当设备距离被测物体40 cm以内时, 系统启动图像处理功能，一方面, RGB图像与掩模寄存器中保存的掩模进行融合, 并将带有残留物标记的图像传送到RGB触控屏中，另一方面，RGB图像通过分块大津算法生成新的残留物掩模并更新掩模寄存器; 而当距离超过40 cm时，RGB图像直接显示在RGB触控屏；在RGB屏中可以选择是否将融合后的荧光图像进行保存。
3.2 硬件逻辑设计
3.2.1 相机配置
当FPGA初始化完成后，FPGA需要配置相机的252个寄存器。相机配置状态机如图 10所示，状态机从idle跳转到wrreg_req状态并向IIC控制器发送写命令，IIC控制器开始向相机一个寄存器地址写入命令，写入完成后rw_done置1，状态机跳转为judge状态并判断相机寄存器配置完成个数，当计数器cnt小于252时，状态机跳转回wrreg_req状态继续写命令; 当cnt等于252时，状态机跳转到finish状态，从而完成相机的初始化配置。
相机部分配置参数如表 1所示。其中16进制数0 x 3808~0 x 380 b为图像分辨率寄存器地址范围，0 x 3503为手动曝光寄存器地址，0 x 3500~0 x 3502为曝光时长寄存器地址范围，0 x 4300为格式控制寄存器地址，0 x 501 f为图像输出格式选择寄存器地址。
表 1 相机主要配置参数Table 1. Main camera configuration parametersregister name address value TIMING DVPHO H 0x3808 0x03 TIMING DVPHO L 0x3809 0x20 TIMING DVPVO L 0x380a 0x01 TIMING DVPVO L 0x380b 0xe0 AEC PK MANUAL 0x3503 0x03 AEC PK EXPOSURE 0x3500 0x00 AEC PK EXPOSURE 0x3501 0x6e AEC PK EXPOSURE 0x3502 0x30 FORMAT CONTROL 0x4300 0x61 ISP FORMAT 0x501f 0x01 3.2.2 相机接口逻辑设计
相机接口逻辑主要完成相机RGB565格式输出数据的并行接收，其接收功能主要由10个不同位宽寄存器和2个数据选择器组成。如图 11所示，其中寄存器r_Data、r_Href、r_Vsync分别对输入数据Data、Href、Vsync、进行打一拍操作。H_count对每行像素数据进行计数，当H_count为偶数时将r_Data中8位图像数据保存在r_DataPixel高8位中；H_count为奇数时将r_Data中8位图像数据保存在r_DataPixel低8位中。在完成16位数据接收后，寄存器r_DataValid置1，表明此时r_DataPixel中16位图像数据为正确数据，可以输出到DataPixel口。为确保图像输出稳定，需要舍弃系统开始运行前10帧图像。寄存器FrameCnt对每帧图像进行计数，当FrameCnt计数到第10帧后dump_frame置1，表明此后图像有效。
3.2.3 触控屏接口逻辑设计
触控屏接口逻辑是将RGB565格式图像数据流转换成VGA模式下的并行数据并发送到触控屏中。如图 12所示，根据触控屏接收时序要求，寄存器Hcount_r和Vcount_r分别产生行同步和场同步信号；Disp_Hs、Disp_Vs和Disp_De根据Hcount_r和Vcount_r产生数据有效信号；寄存器Disp_Red、Disp_Green、Disp_Bule分别将RGB565格式视频数据转换为并行数据。寄存器Frame_begin在场同步信号Disp_Vs上升沿置1，表明开始显示一帧图像。
3.3 分块大津算法原理及FPGA实现
3.3.1 算法原理
当完成荧光图像的采集后，需要利用图像分割算法实现残留物的准确分割。由于具有计算简单、不受图像亮度和对比度影响的优点，大津算法在图像分割领域得到广泛引用。然而大津算法是一种全局阈值分割方法, 其在非均匀光照下会产生错误分割，并且当感兴趣目标灰度分布范围较大且部分目标灰度接近背景强度时，接近背景强度的部分目标会检测丢失。针对上述缺点，本文中基于分块大津算法的思想，通过对检测图像的合理分块，避免了背景子块误分割，提高了设备在不均匀光照或目标灰度差异较大情况下的检测能力。分块大津算法流程如图 13所示。首先FPGA从SDRAM缓存中读取荧光图像并经过灰度转换模块转换为灰度图像；然后FPGA按照行同步、场同步信号将一幅灰度图像划分M行N列的子块, 并运用大津算法计算每个子块分割阈值和类间平均灰度差；进一步对每个子块的类间平均灰度差进行判断，当子块类间平均灰度差较小时，该子块前景和背景灰度差异不大，判断为背景块，该子块分割阈值复位为0；当子块类间平均灰度差较大时，该子块中存在荧光残留物，保留计算的分割阈值。在系统实际使用中分块大津算法在设备距离被测表面40 cm以内启动，系统根据40 cm以内激发光源光强特点结合FPGA硬件资源占用将图像分为5行8列的子块进行处理，将判断子块是否为纯背景的类间平均灰度差设置为8。
![图 13 荧光残留物提取以及掩膜生成流程图]()
图 13 荧光残留物提取以及掩膜生成流程图Figure 13. Flow chart of fluorescence residue extraction and mask generation3.3.2 FPGA实现
分块大津算法在FPGA中实现如图 14所示。FPGA以数据流形式从SDRAM中读取的荧光图像Rgb_data缺少位置信息；为实现荧光图像准确分块，FPGA在统计模块中对经过的荧光图像数据流进行坐标统计并生成行同步信号Href、列同步信号Vsync和数据有效信号De；16位荧光图像数据经过灰度转换模块转换为8位灰度数据；系统通过调用8个大津算法模块和串行流水线设计方法实现分块大津算法计算。
分块大津算法中，使用大津算法模块计算每个子块的阈值和类间平均灰度差。如图 14中虚线框中所示，直方图统计模块具有5个输入，4个输出，其中输入Gray_data为8位图像灰度数据。输入Data_end在直方图统计完成后拉高，灰度阈值求完后拉低。输入Data_valid只在灰度直方图统计时为高，其余拉低。输入T为分割阈值, 输入Clr为清零信号，当灰度阈值求完后拉高使直方图统计模块清零复位等待下一次统计。在直方图统计模块中输出表达式[20]为：
S0=t∑i=0ni (1) S1=255∑i=t+1ni (2) G0=t∑i=0ini (3) G1=255∑i=t+1ini (4) 式中，i为灰度值，t为灰度值，ni为灰度值等于i的像素个数，S0为灰度值不大于t的像素个数，S1为灰度值大于t的像素个数，G0为灰度值不大于t的像素灰度总和，G1为灰度值大于t的像素灰度总和。直方图统计模块输出经过后续乘法器、除法器、减法模块计算后得到类间方差，在比较器模块中保存类方间差最大时的阈值T，输出此时阈值以及类间灰度差。
4. 系统测试
为了测试本文中所开发系统的适用性，本文作者选用食品加工场所中3种常见包装和加工材质(木制案板、不锈钢板、聚乙烯材质塑料板)进行测试，其中木制案板和不锈钢板不产生荧光，聚乙烯(polyethylene, PE)材质塑料菜板会产生荧光[17]。测试对象为5种不同体积分数的奶粉溶液和菠菜汁液，奶粉溶液采用全脂牛奶粉和纯净水按照不同体积分数(50%、33%、20%、10%、2%)配制；菠菜汁液采用新鲜菠菜榨汁[17]和纯净水按照不同体积分数(50%、33%、25%、20%、10%)配制。由于奶粉在516 nm处存在荧光发射峰[17]，菠菜汁在685 nm处存在荧光发射峰[17, 21]，因此针对奶粉选用520 nm窄带宽滤波片，针对菠菜汁选用670 nm窄带宽滤波片。
木制案板表面不同体积分数奶粉溶液和菠菜汁残留物检测效果如图 15所示。其中图 15a和图 15e分别是环境光下采集的奶粉残留物和菠菜汁液残留物RGB图像；图 15b和图 15f分别为激发光源下采集的奶粉残留物和菠菜汁液残留物荧光图像；图 15c和图 15g分别是分块大津算法生成的奶粉残留物和菠菜汁液残留物掩模图像；图 15d和图 15h分别是奶粉残留物和菠菜汁液残留物最终检测结果。如图 15a、图 15e所示，环境光下，当奶粉、菠菜汁液体积分数分别低于20%和10%时，肉眼几乎不可见木制案板上的残留物；而奶粉残留物、菠菜汁液与木制案板的荧光图像有着明显的区别(见图 15b和图 15f), 在图 15f中, 体积分数为50%和33%的菠菜残留物由于存在大量菠菜纤维等固形物导致其荧光强度较弱; 图 15c和图 15g中给出了分块大津算法的检测效果，可以看出，不同体积分数的残留物均被有效检出，且在图 15g中还检测出了部分在图 15e中未被标记的残留物，这些残留物是木板在简单清洗中未被发现的历史残留，这进一步体现了设备对不可见残留物的检测能力。
![图 15 木制案板表面残留物检测效果]()
图 15 木制案板表面残留物检测效果Figure 15. Detection effect of residues on wooden cutting board surface图 16中给出了材质为不锈钢板条件下，奶粉和菠菜汁残留物的检测效果。其中图 16a和图 16e分别是奶粉残留物和菠菜汁残留物RGB图像；图 16b和图 16f分别为奶粉残留物和菠菜汁残留物荧光图像；图 16c和图 16g分别是分块大津算法生成的奶粉残留物和菠菜汁残留物掩模图像；图 16d和图 16h分别是奶粉残留物和菠菜汁残留物最终检测结果。如图 16a和图 16e所示，由于不锈钢光滑镜面反射影响，体积分数为10%的菠菜汁残留物因接近无色而基本不可见；在荧光图像16b和图 16f中不同体积分数的残留物均清晰可见；分块大津算法可较好克服光照的影响，能够准确分割出污染物区域(见图 16c和图 16g)。
PE塑料板上的奶粉和菠菜汁残留物检测效果如图 17所示。其中图 17a和图 17e分别是奶粉残留物和菠菜汁残留物RGB图像；图 17b和图 17f分别为奶粉残留物和菠菜汁残留物荧光图像；图 17c和图 17g分别是分块大津算法生成的奶粉残留物和菠菜汁残留物掩模图像；图 17d和图 17h分别为奶粉残留物和菠菜汁残留物最终检测结果。从图中可以看出，在环境光下，当奶粉体积分数低于33%、菠菜汁液体积分数低于10%时，难以通过RGB图像检测出残留物；而荧光图像可以较好地实现残留物检测。但从图 17c来看，当奶粉体积分数低于20%时，残留物的分割区域存在部分缺失现象，其原因可能在于：PE塑料板会发出和奶粉波段相近的荧光，当奶粉体积分数较低时，其荧光强度与背景接近，导致分割困难。
图 18给出了大津算法与分块大津算法对典型残留物的分割结果。其中图 18a和图 18b分别是奶粉残留物在木制案板和不锈钢板上的荧光图像；图 18c和图 18d分别为菠菜汁残留物在木制案板和不锈钢板上的荧光图像；图 18e、图 18f、图 18g和图 18h分别是图 18a、图 18b、图 18c和图 18d采用分块大津算法计算得到的掩模图像；图 18i、图 18j、图 18k和图 18l分别是图 18a、图 18b、图 18c和图 18d采用大津算法计算得到的掩模图像。绿色框为成功检测出残留物；红色框为未检测或错误检测残留物。在图 18a中检测设备位于木制案板上方，案板中间部分靠近光源的反射光强度更高，左上角远离光源强度较低，这导致在使用大津算法进行分割时图 18i中部分较为明亮的背景被判定为残留物而光照强度较低的残留物则未被检出；在图 18b中, 检测设备位于不锈钢板斜上方，不锈钢板右上角光照强度更高，左下角光照强度较低，大津算法在图 18j中将靠近光源的残留物检出而远离光源且体积分数较低的残留物未检出。大津算法在检测灰度范围较大的目标时会出现漏检；在木制案板(见图 18c)和不锈钢板(见图 18d)上菠菜汁荧光强度分布范围较大, 且部分菠菜汁残留物荧光强度接近背景光强；当使用大津算法进行检测时，部分荧光强度较低的残留物漏检(见图 18k), 而尽管在图 18l中所有残留物均被检出, 但是相较于分块大津算法结果(见图 18h)，大津算法分割出的目标不够完整; 图 18e、图 18f、图 18g和图 18h为分块大津算法检测结果，不同体积分数残留物均被检出；相较于大津算法检测结果，分块大津算法提高了设备在不均匀光照下的检测效果以及对于荧光强度分布范围较大检测目标的检测能力。
![图 18 分块大津算法和大津算法检测效果对比图]()
图 18 分块大津算法和大津算法检测效果对比图Figure 18. Comparison of detection effect between block OTSU and OTSU5. 结论
基于荧光成像技术开发了一种能够协助现场检测员对食品加工场所进行视觉卫生检查的手持式荧光成像检测设备，该设备具有荧光图像采集、图像处理、图像显示、图像保存功能，能够实时地将检测到的污染区域进行标记, 并根据需要保存在SD卡中。此外，还可以通过更换不同波长的滤波片实现多种食品残留物的检测。基于改进的分块大津算法，提高了设备在不均匀光照下的检测效果。在3种不同材质表面的实测证明，设备相较于检测员目视能够更清晰地观察到食品残留物，能够显著提高卫生检测效果，可以协助检测员进行卫生安全检查。受现场检测环境的限制，所开发的设备与待检测样本为非接触测量，由于现场多种因素的不确定性(如环境光源的干扰、手持式测量装置与被测对象的距离、角度等)，导致利用荧光图像反演残留物体积分数的精度还存在较大困难。
在今后的研究中，将深入分析影响残留物体积分数反演精度的关键要素和规律，构建补偿模型，以实现残留物体积分数的高精度反演，此外，还将实现云组网功能以及将荧光图像与深度学习算法相结合实现自动检测自动分割。
-
![]()
图 1 OCT熔深检测系统原理图
Figure 1. Schematic diagram of OCT laser weld depth monitoring system
![]()
图 4 样例A的焊缝纵抛面的金相图片
Figure 4. Metallographic picture of the longitudinal throwing surface of the weld of sample A
![]()
图 6 OCT熔深原始数据与实际熔深曲线的对比
Figure 6. Comparison of OCT raw data and actual depth of fusion curves
![]()
图 7 不同百分位数的熔深曲线与实际熔深的局部对比
Figure 7. Local comparison of weld depth curves at different percentiles with actual weld depth curves
![]()
图 8 百分位数p=5%的熔深曲线与实际熔深的全局对比
Figure 8. Global comparison of weld depth extracted by percentile filter p=5% with actual weld depth
![]()
图 9 基于百分位数p=5%和移动平均算法提取的熔深曲线
Figure 9. Weld depth curve extracted based on percentile filtering p=5% and moving average algorithm
![]()
图 10 拟合曲线与实际熔深曲线的差异评估
Figure 10. Evaluation of the difference between the fitted curve and the actual melt depth curve
![]()
图 11 样例B和样例C的金相图片、百分位滤波和基于百分位滤波及移动平均算法提取的熔深曲线
Figure 11. Metallographic images of sample B and sample C, weld depth curves extracted based on percentile filtering and weld depth curves extracted based on percentile filtering and moving average algorithm
表 1 3组实验样例的平均误差
Table 1 Average error of the three experimental samples
experimental sample percentile filter/% percentile filter+MA/% accuracy improvement/% sample A 4.93 4.21 15 sample B 4.97 4.34 13 sample C 4.53 4.09 10 
下载: 导出CSV
-
[1] XIAO X, LIU X, CHENG M, et al. Towards monitoring laser welding process via a coaxial pyrometer[J]. Journal of Materials Processing Technology, 2020, 277: 116409. DOI: 10.1016/j.jmatprotec.2019.116409
[2] SANDERS P G, LEONG K H, KESKE J S, et al. Real-time monitoring of laser beam welding using infrared weld emissions[J]. Journal of Laser Applications, Laser Institute of America, 1998, 10(5): 205-211. DOI: 10.2351/1.521853
[3] WANG L, MOHAMMADPOUR M, GAO X, et al. Adjustable ring mode (ARM) laser welding of stainless steels[J]. Optics and Lasers in Engineering, 2021, 137: 106360. DOI: 10.1016/j.optlaseng.2020.106360
[4] 李竹曼, 高向东, 张南峰. 大功率碟片激光焊接状态多特征融合分析法[J]. 激光技术, 2017, 41(5): 764-768. DOI: 10.7510/jgjs.issn.1001-3806.2017.05.029 LI Zh M, GAO X D, ZHANG N F. Analysis of high-power disk laser welding status based on multi-feature fusion[J]. Laser Technology, 2017, 41(5): 764-768(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2017.05.029
[5] FAN X, GAO X, LIU G, et al. Research and prospect of welding monitoring technology based on machine vision[J]. The International Journal of Advanced Manufacturing Technology, 2021, 115(11/12): 3365-3391.
[6] ZHANG X, TANG Z, WU Y, et al. Progress in in situ X-ray imaging of welding process[J]. Review of Scientific Instruments, Melville: AIP Publishing, 2022, 93(7): 071501. DOI: 10.1063/5.0074042
[7] HUANG D, SWANSON E A, LIN C P, et al. Optical coherence tomography[J]. Science, 1991, 254(5035): 1178-1181. DOI: 10.1126/science.1957169
[8] DUPRIEZ N D, DENKL A. Advances of OCT technology for laser beam processing: Precision and quality during laser welding[J]. Laser Technik Journal, 2017, 14(4): 34-38. DOI: 10.1002/latj.201700021
[9] 刘逸飞, 苏亚, 姚晓天, 等. OCT无创血糖检测图像处理最优化方法研究[J]. 激光技术, 2023, 47(2): 178-184. DOI: 10.7510/jgjs.issn.1001-3806.2023.02.004 LIU Y F, SU Y, YAO X T, et al. An optimization method of image processing for OCT non-invasive blood glucose detection[J]. Laser Technology, 2023, 47(2): 178-184(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.02.004
[10] MIYAGI M, KAWAHITO Y, KAWAKAMI H, et al. Dynamics of solid-liquid interface and porosity formation determined through X-ray phase-contrast in laser welding of pure Al[J]. Journal of Materials Processing Technology, 2017, 250: 9-15. DOI: 10.1016/j.jmatprotec.2017.06.033
[11] LOHAUS L, BAUTZE T, DIEPOLD K. Evaluation of optical sensors for laser welding in a technical cognitive environment[C]//International Congress on Applications of Lasers & Electro-Optics. Anaheim, California, USA: Laser Institute of America, 2010: 1541-1546.
[12] BLECHER J J, GALBRAITH C M, van VLACK C, et al. Real time monitoring of laser beam welding keyhole depth by laser interferometry[J]. Science and Technology of Welding and Joining, 2014, 19(7): 560-564. DOI: 10.1179/1362171814Y.0000000225
[13] DORSCH F, HARRER T, HAUG P, et al. Process control using capillary depth measurement[C]// International Congress on Applications of Lasers & Electro-Optics. San Diego, California, USA: Laser Institute of America Press, 2016: 1505.
[14] 谢冠明, 王三宏, 张跃强, 等. 基于光学相干层析的激光焊接熔深检测方法[J]. 光学学报, 2023, 43(11): 1114002. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202311014.htm XIE G M, WANG S H, ZHANG Y Q, et al. Laser welding depth monitoring method based on optical coherence tomography[J]. Acta Optica Sinica, 2023, 43(11): 1114002(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202311014.htm
[15] 黄威威, 游德勇, 高向东, 等. 基于相关分析和神经网络的激光焊接稳态识别[J]. 激光技术, 2022, 46(3): 312-319. DOI: 10.7510/jgjs.issn.1001-3806.2022.03.004 HUANG W W, YOU D Y, GAO X D, et al. Laser welding steady status recognition method based on correlation analysis and neural network[J]. Laser Technology, 2022, 46(3): 312-319(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2022.03.004
[16] SOKOLOV M, FRANCIOSA P, SUN T, et al. Applying optical coherence tomography for weld depth monitoring in remote laser welding of automotive battery tab connectors[J]. Journal of Laser Applications, Laser Institute of America, 2021, 33(1): 012028.
[17] DUIN R, HARINGA H, ZEELEN R. Fast percentile filtering[J]. Pattern Recognition Letters, 1986, 4(4): 269-272.
[18] BOLEY M, FETZER F, WEBER R, et al. Statistical evaluation method to determine the laser welding depth by optical coherence tomography[J]. Optics and Lasers in Engineering, 2019, 119: 56-64.
[19] MITTELSTÁDT C, MATTULAT T, SEEFELD T, et al. Novel approach for weld depth determination using optical coherence tomography measurement in laser deep penetration welding of aluminum and steel[J]. Journal of Laser Applications, 2019, 31(2): 022007.
[20] SCHMOELLER M, STADTER C, LIEBL S, et al. Inline weld depth measurement for high brilliance laser beam sources using optical coherence tomography[J]. Journal of Laser Applications, Laser Institute of America, 2019, 31(2): 022409.
-
期刊类型引用(2)
1. 周涛，高承帅，王建，聂潇乾，王炳琪，吴彬彬. 一维摄像机的线结构光测头标定方法. 机电工程技术. 2024(07): 83-86 . 
百度学术
2. 李勇，周哲海，李忠祥，陈丽. 基于非线性算法的扫描振镜坐标轴标定. 激光杂志. 2024(12): 67-73 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 422
- HTML全文浏览量: 5
- PDF下载量: 7
- 被引次数: 4