Processing math: 100%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高动态范围条纹结构光在机检测技术及应用进展

刘泽隆 李茂月 卢新元 张明垒

刘泽隆, 李茂月, 卢新元, 张明垒. 高动态范围条纹结构光在机检测技术及应用进展[J]. 中国光学(中英文), 2024, 17(1): 1-18. doi: 10.37188/CO.2023-0068
引用本文: 刘泽隆, 李茂月, 卢新元, 张明垒. 高动态范围条纹结构光在机检测技术及应用进展[J]. 中国光学(中英文), 2024, 17(1): 1-18. doi: 10.37188/CO.2023-0068
LIU Ze-long, LI Mao-yue, LU Xin-yuan, ZHANG Ming-lei. On-machine detection technology and application progress of high dynamic range fringe structured light[J]. Chinese Optics, 2024, 17(1): 1-18. doi: 10.37188/CO.2023-0068
Citation: LIU Ze-long, LI Mao-yue, LU Xin-yuan, ZHANG Ming-lei. On-machine detection technology and application progress of high dynamic range fringe structured light[J]. Chinese Optics, 2024, 17(1): 1-18. doi: 10.37188/CO.2023-0068

高动态范围条纹结构光在机检测技术及应用进展

基金项目: 国家自然科学基金资助项目(No. 51975169);黑龙江省自然科学基金资助项目(No. LH2022E085)
详细信息
    作者简介:

    刘泽隆(1996—),男,黑龙江大庆人,博士研究生,2022 年于哈尔滨理工大学获得硕士学位,主要从事图像处理和机器视觉方面的研究。E-mail:LZL_LOUIS1231@163.com

    李茂月(1981—),男,山东青岛人,博士,教授,博士生导师,2004 年于南京林业大学获得学士学位,2007 年于长安大学获得硕士学位,2012 年于哈尔滨工业大学获得博士学位,主要从事智能加工与光学检测技术方面的研究。E-mail:lmy0500@163.com

  • 中图分类号: TH741

On-machine detection technology and application progress of high dynamic range fringe structured light

Funds: Supported by National Natural Science Foundation of China (No. 51975169); Natural Science Foundation of Heilongjiang Province(No. LH2022E085)
More Information
  • 摘要:

    条纹结构光技术是近年来发展迅速的非接触式测量方法,为机械加工在机检测提供了新的解决方案。由于加工环境光线复杂且金属零件本身具有高反光的特性,造成结构光在机检测的精度降低。将高动态范围(High Dynamic Range,HDR)技术应用于结构光检测中,可抑制高反光的影响,实现金属零件在复杂场景的测量。本文首先介绍了结构光测量原理,总结出HDR结构光在机检测面临的难点;其次,对HDR结构光技术进行了全面综述,以机械加工在机检测为背景,对基于硬件设备的HDR技术和基于条纹算法的HDR技术分别进行了归纳分析;然后,根据在机检测的条件需求,对各类技术进行总结,并比较不同方法的优缺点和在机检测的适用性;最后,结合近年来先进制造技术和精密测量的研究热点,对潜在应用进行分析,提出技术展望。

     

  • 航空发动机涡轮叶片、核电汽轮机大叶片、大口径光学镜面等典型的复杂曲面零件,几何精度和物理性能要求高,在机械加工过程中极易产生变形[1]。因此,对零件进行在机不拆卸检测,高效地检测加工质量,确保其加工精度满足要求,成为近年来先进制造领域的研究热点。传统的测量方式是将零件从机床上卸载,置于三坐标测量机(Coordinate Measuring Machine, CMM)上,通过接触式测量原理进行检测。这种方法需要将工件多次拆卸重新装夹,容易带来定位误差和装夹变形,且测量探头易对工件表面造成损伤。随着先进制造技术的推进,出现了可对零件进行非接触式的在机检测,可有效避免上述缺点,在机检测技术即在工件装夹定位状态不变的情况下,对工件进行原位检测。以光学检测为代表的非接触式测量可有效避免反复拆卸带来的影响,同时测量速度快,具备一定的实时性,测量结果可用于指导加工参数的修正[2]

    光学测量设备主要有激光干涉仪、结构光测量仪、显微共焦仪、白光干涉仪等。其中,结构光技术具有测量原理简单、硬件成本低、测量精度高、测量速度快等优点,在工业和医学等领域受到广泛应用。另外,结构光技术算法易于集成,相关硬件可以通过与机械臂配合进行测量视点和路径规划[3],受待测件尺寸的影响较小。

    基于结构光技术的测量效果受到相机分辨率、投影分辨率、标定精度等因素影响。其测量精度很大程度上取决于被测物体表面的光学反射特性。由于加工过程中环境光线复杂,且金属表面具有较强的镜面反射,使用传统结构光方法测量时,工业相机拍摄的条纹图像易出现过饱和,导致条纹数据信息丢失,造成重建的点云出现大面积的波纹和孔洞,无法准确得到待测物的三维信息。在商用设备市场中,著名测量系统供应商德国GOM和瑞典Hexagon通过在待测件表面喷涂一种抑制强反射的涂层,来改变材料的反射特性,使其表面由镜面反射变为漫反射,从而削弱反光,但是涂层的厚度和均匀程度主要取决于操作者的经验,极容易带来测量误差,降低测量精度。国内外学者为了解决这一问题,进行了多方面的研究。研究发现相机所拍摄条纹图像呈现过饱和现象,是因为待测件的反射强度变化程度高于相机的0~255灰度强度范围。因此用于高反光表面测量的方法被称为高动态范围(High Dynamic Range, HDR)技术。

    在结构光检测方面,现有的综述大多介绍了不同的检测原理和方法[4-5],针对HDR结构光测量技术,现有文献虽然进行了技术总结,讨论了不同技术的优缺点[6-7],但未与在机加工检测结合,且随着智能技术的成熟,近两年国内外学者将深度学习等智能算法应用于结构光检测中,解决了实时在机检测、动态物体检测和HDR检测等问题。本文以机械加工在机检测为背景,对近年来用于测量高反光表面的HDR技术进行了全面综述。本文将现有技术方案归纳为基于硬件设备的HDR技术和基于条纹算法的HDR技术两大类,针对不同方法的特点,结合对加工场景在机检测的需求,总结了不同方法的技术特点,并且对HDR结构光在机检测进行了技术展望。

    下文安排如下:第二部分介绍HDR条纹测量原理和在机检测关键问题;第三、四部分分别介绍基于硬件设备的HDR技术和基于条纹算法的HDR技术;第五部分结合加工在机检测方面面临的问题,对现有方法进行分析比较,总结特点;第六部分介绍HDR结构光在机检测技术的潜在应用和技术展望;第七部分为结论。

    与传统接触式测量技术相比,三维视觉测量技术具有速度快、成本低的特点,为我国先进制造、精密工程发展战略提供新的技术支撑。如表1所示,三维视觉测量技术根据测量过程中是否投射光源,可分为被动视觉测量和主动视觉测量。其中,面结构光式测量方法在所有结构光方法中测量效率最高,广泛应用于物体的三维测量领域。

    表  1  三维视觉测量技术分类
    Table  1.  Classification of three-dimensional visual measurement technology
    视觉测
    量分类
    是否投
    射光源
    具体分类 特点
    被动视
    觉测量
    单目视觉测量 基于图像聚焦程度完成三维重建,多用于显微视觉测量中。
    双目视觉测量 根据三角测量原理实现三维重建,应用于双目立体摄像头。
    多目视觉测量 增加辅助相机,通过光束平差提高测量精度。
    主动视
    觉测量
    点扫描式 激光器投射光点,根据光标中心坐标和标定数据进行重建,测量效率低。
    线扫描式 激光器投射光条代替光点,提高效率,广泛应用于激光扫描仪中。
    面扫描式 通过投影仪投射二维结构光,单次投射覆盖区域大,测量效率最高。
    下载: 导出CSV 
    | 显示表格

    结构光测量系统主要由工业相机、投影仪和计算机组成,其原理示意图如图1(彩图见期刊电子版)所示。

    图  1  结构光测量原理示意图[8]
    Figure  1.  Schematic diagram of monocular structured light measurement principle[8]

    条纹轮廓术(Fringe Projection Profilometry, FPP)是常用的面结构光测量方法。其通过计算机获得编码条纹,经投影仪投射到待测物表面,解析变形条纹即可得到相位信息。在经典的FPP测量系统中,常采用相移法求解相位,条纹图案为正弦光栅,条纹强度Ii(x,y)可以表示为:

    Ii(x,y)=a(x,y)+b(x,y)cos[φ(x,y)+2πiN](i=0,1,N1),
    (1)

    式中,(x,y)是图像像素坐标,a(x,y)是条纹图像的背景光强;b(x,y)是受物体表面反射率影响的调制光强;i表示相移步数;N表示相移总步数;φ(x,y)为待求相位主值,其表达式为:

    ϕ(x,y)=arctanNi=1In(x,y)sin2πiNNi=1In(x,y)cos2πiN.
    (2)

    得到相位主值后,通过多频外差法获得连续相位,将连续相位与标定数据相结合,即可得到待测物的三维信息。

    在信息互联的智能制造背景下,利用AGV小车和六自由度机器人搭载结构光测量设备,对加工过程中的零件进行在机检测(如图2所示)。该方法可避免反复拆卸产生定位误差和装夹变形,同时可提高零件整体的加工效率,有助于推进智能化加工(如图3所示)。

    图  2  结构光在机检测示意图
    Figure  2.  Schematic diagram of structured light on-machine detection
    图  3  结构光在机检测系统架构
    Figure  3.  Architecture of on-machine detection system of structured light

    在加工环境中,传统结构光技术易受到复杂环境光线和金属零件高反光的影响。但将HDR技术应用于结构光在机检测中,可有效抑制相关影响,提升测量精度。结构光系统中,相机成像模型如公式(3)所示[8]

    I=αt(βIα1+βIp+Iα2)+μ,
    (3)

    式中:I是图像像素强度;α是相机灵敏度系数;t是相机的曝光时间;β是待测物的反射系数;I p是投影仪的光强;$ {I^{{\alpha _1}}} $是被测物体反射的环境光;$ {I^{{\alpha _2}}} $是直接进入相机的环境光;µ是相机的噪声误差。

    根据该模型可知,通过改变相机的曝光量、条纹的光强和物体的反射特性,可以改进拍摄条纹图像的成像效果。基于此原理,国内外学者通过改变相机曝光度、调整投影图案强度,加置偏振滤光片,利用相位测量偏折术、光度立体技术,更改编码解码方案和深度学习智能算法等多种方法,对HDR测量技术进行研究。在使用HDR结构光技术进行在机检测时,需重点考虑如下问题:

    (1)零件加工时,整体处于暗房环境,但受到环境光和金属工件镜面反射影响,光线条件较为复杂,在抑制高反光的同时,不能忽略暗区域的影响。

    (2)测量系统的硬件设备不能过于复杂,要易于配合机床、机械臂、AGV小车,同时零件检测的路径要易于规划,使检测方案更加集成化、智能化。

    (3)测量方案简单的同时,对测量速度和效率要求较高,要具有实时性,可及时反馈测量结果,指导加工参数的补偿修正。

    下面结合在机检测需求,详细综述各类方法的研究进展。

    基于硬件设备的HDR技术是通过对结构光硬件设备进行改进,从而降低高反光对其测量结果的影响,主要可分为相机曝光法、偏振滤光片法、相位测量偏折术以及光度立体技术。

    相机的曝光是成像质量的重要影响因素,若曝光过度,则图像整体过饱和;若曝光不足,则图像整体过暗,两者都会使图像丢失细节,无法准确获得条纹信息。Zhang等人[9]提出了一种快速自动曝光的方法。通过捕获一次曝光的条纹图像,来确定最佳曝光时间。所提方法虽然可以快速找到最佳曝光时间,但是适用范围受限,对于整体偏亮或偏暗的图像有较好的效果,对于加工过程中明暗混合的复杂场景,难以确定一个统一的最佳曝光时间,影响在机检测质量。

    多重曝光[10-11]技术是指将不同曝光时间下采集的图像,融合成一幅图像以避免图像饱和,从而达到较高的信噪比(Signal Noise Ratio, SNR)。Zhang等人[12]提出将多重曝光技术应用于三维测量,通过调整相机镜头光圈获得一组不同曝光程度的条纹图像。对陶瓷花瓶检测时,将高曝光的饱和像素用低曝光下相应的像素替换。所提方法不仅适用于相移法,对于莫尔条纹、激光干涉法同样适用。

    在传统多重曝光技术中,对于曝光度的选择主要依靠经验。若曝光范围选取不好,直接影响最终测量结果。为了避免依据经验选择曝光度的不可靠性,国内外学者通过相机响应函数曲线选择适合的曝光时间。Song等人[13]提出利用相机响应函数来压缩动态范围,使用中等曝光数据作为参考进行图像融合。对不锈钢冲压件通过6次曝光进行检测实验,平均偏差(Mean Absolute Error, MAE)为0.06 mm。Feng等人[14]基于灰度直方图分布将测量出的表面反射率细分为几个组。根据相机响应函数自适应地预测每个组的最佳曝光时间。Cui等人[15]利用不同曝光时间的序列图像,标定相机响应曲线函数。基于最高灰度值和最低灰度值双阈值原理,提出了曝光时间的自适应选择算法。实验表明,所提方法自动计算了4次曝光时间,避免了基于经验调整曝光时间的缺点,提高了曝光的自动化水平和测量精度。

    考虑获取相机响应曲线并记录曝光时间较为繁琐,且相机响应曲线受环境光影响较大,存在误差。Rao[16]通过分析每个像素的条纹调制强度,设定调制阈值来自动计算多个曝光时间。实验表明,该方法最多在5次曝光下,即可完成HDR三维测量。Wu等人[17]提出了一种曝光融合方法,在水平和竖直两个方向投射二值格雷码,算法流程如图4所示。通过拉普拉斯金字塔分解图像的加权平均值,以高斯金字塔分解图像权重,得到融合金字塔,从而获得最终图像。实验对圆柱形电池通过8次曝光进行检测,MAE值为0.0880 mm,测量时间为1436 s,所提方法虽然无需拟合相机响应曲线,但是需要多组曝光数据,较为耗时。

    图  4  曝光融合算法流程[17]
    Figure  4.  The process of exposure fusion algorithm[17]

    基于相机曝光技术进行结构光在机检测时,无需额外添加硬件,只需对图像进行融合处理。其核心在于曝光参数的选择,单次最佳曝光难以适用于复杂场景,而利用多重曝光的方式,需要合理选择曝光时间实现对待测工件的多次测量,影响了检测效率,且曝光时间不能随意设定,应保证与投影仪的刷新率精确同步。

    自然光经物体反射会变为部分偏振光,利用偏振成像获得反射光的偏振状态成为近年来计算机视觉的热点方向。将偏振片置于工业相机镜头前进行在机检测,可有效地达到过滤高光的目的。

    使用单组偏振片会将图像整体的亮度降低,影响图像的SNR,国内外学者采用多组偏振片和图像融合的方式,来实现对高光的抑制并提高条纹图像的SNR。Salahieh等人[18]通过选择不同的偏振测量值或偏振角和曝光时间的正确组合,消除饱和或低对比度条纹区域,以保持条纹良好的可见性。平茜茜等人[19]在两个CCD相机前加装偏振片,将一个固定成0°,另一个调成4个角度。通过4幅图像的偏振信息得到梯度信息,对其积分后得到用于重构的深度信息。根据杂散光和漫反射光具有不同的光谱特性和偏振特性,郝婧蕾[20]对待测目标进行高光消除预处理,从而提高三维重构精度。Wang等人[21]基于条纹图像生成辅助图像,确定最佳曝光时间和偏振方向,在适当的曝光时间内提取4种偏振方向生成复合相移图像,消除饱和以及低对比度区域。Zhu等人[22]建立了相机在偏振系统下的强度响应函数,避免了复杂的偏振双向反射分布模型,可直接计算偏振片之间的角度,通过图像融合算法生成最佳条纹。该方法显著消除了高光的影响,并减少了黑白条纹之间的模糊过渡区域。

    上述方式与多重曝光类似,需要相机进行多次拍摄,影响在机检测的效率。因此,国内外学者利用不同偏振方向的多组相机和投影仪,同时获取不同偏振角度的条纹图像。Maeda等人[23]将偏振相机前的偏振片以0°、45°、90°和135°放置,通过1/4波形板和空间光调制器与投影仪组合产生空间变换的偏振图案。偏振相机可同时获取4个偏振方向的相移图像,从而得到三维信息,MAE结果为0.04 mm。Xiang等人[24]设计了一种由两个相机和两个投影仪构成的偏振测量系统,如图5所示。在测量过程中,利用入射光偏振的特性,使用两个同轴排列的偏振相机同时捕获被测表面上的变形条纹。采用两步相移算法计算相位,得到三维形状。实验均方根偏差(Root Mean Square Error, RMSE)为0.17 mm。

    图  5  Xiang设计的偏振测量系统[24]
    Figure  5.  Polarization measurement system designed by Xiang[24]

    基于偏振滤光片的HDR技术易降低图像的整体亮度。这是因为其在削弱镜面反射的同时也削弱了漫反射,导致图像SNR下降。由于加工环境光线较为复杂,此类方法不适用于加工现场的在机检测。虽可通过选择多组偏振片进行数据融合的方式得到较好的条纹图像,但是降低了整体的测量效率。另外,对于偏振片角度的选择,需依靠人为经验并仔细调整光路,整体结构较为复杂。

    工业生产中的镜面和类镜面零件具有极强的镜面反射,传统结构光技术中,工业相机无法直接清晰捕获投影仪投射的条纹结构光。相位测量偏折术(Phase Measuring Deflectometry, PMD)是基于条纹反射原理,将传统的结构光测量系统中的投影仪用一个LCD显示屏代替,使计算机生成的条纹图显示在LCD显示屏上。相机通过被测物的镜面反射,捕捉所显示条纹图案的反射图像,再通过对捕获的图像求解相位信息,来重建被测物形状[25]

    传统的PMD技术中,LCD显示屏必须与被测物平行放置,占用了很大的测量空间,不适合紧凑的在机测量。Zhuang等人[26]提出了一种紧凑的在机PMD技术。该方法将液晶屏垂直于被测镜面放置,用于实时测量表面形状畸变(如图6所示)。该方法显著扩展了传统PMD方法的应用范围,适用于封闭式激光设备和组装式激光系统中的实时测量。Gao等人[27]利用平板分束器使系统更加紧凑。与传统配置相比显著减小了系统体积,测量精度在0.001 mm之内。传统的PMD系统中,通常采用平面LCD显示屏,这限制了被测表面的可测量曲率。针对这一问题,Han等人[28]将弯曲的液晶显示器用于PMD系统中,并利用标定参数和种子点将两个摄像头视角下重建的三维模型进行点云配准。通过仿真分析,验证了系统的曲率测量范围大于平面LCD显示屏的曲率测量范围。

    图  6  封闭激光腔内COPMD测量系统[26]
    Figure  6.  COPMD measurement system in an enclosed laser cavity[26]

    在大口径光学镜面的加工检测领域,美国亚利桑那大学Su等人[29]在2010年基于条纹反射原理,首次提出了软件可调式光学检测系统(Software Configurable Optical Test System, SCOTS)。在对130 mm离轴抛物面进行测量时,得到了RMSE为1 μm的测量结果。中国科学院长春光学精密机械与物理研究所的邵山川团队[30]利用条纹反射法对超精密金刚石车削反射镜进行在位面形检测。对口径为100 mm凹球面金属反射镜进行检测,结果表明:面形检测精度(峰谷值)优于1 μm;与干涉仪检测结果相比,RMSE测量结果为0.089 μm,其检测系统装置如图7所示。在大口径非球面镜面测量方面,袁婷[31]等人用条纹反射法对口径为821.65 mm抛物面反射镜进行了检测,实验装置如图8所示。检测结果与干涉测量结果偏差仅为0.183 μm。这证明条纹反射法应用于大口径非球面检测的可行性。Navarro等人[32]利用SCOTS对直径为4.2 m的Daniel K.Inouye太阳望远镜进行检测,其RMSE结果小于20 nm。

    图  7  条纹反射在位面形检测系统中的装置[30]
    Figure  7.  Fringe reflection setup in on-machine surface measurement system[30]
    图  8  大口径抛物面反射镜实验装置图[31]
    Figure  8.  Experimental device diagram of large diameter parabolic reflector[31]

    由于PMD技术测量原理特殊,一般用于测量连续的镜面工件,如手机屏幕、天文望远镜镜面等时,测量精度较高,可达到亚微米级别,且不需太多的前期准备和后续处理工作。通过合理设计硬件布置方案,即可实现镜面工件的在机检测。但对于具有复杂特征的金属零件的加工过程,零件表面反射特性表现为镜面反射与漫反射相结合,此时,无法反射出变形的条纹,在这种情况下,此类方法不适用。

    光度立体(Photometric Stereo,PS)技术于上世纪80年代由Woodham提出[33],其原理是保持观察方向不变,在不同方向照明下,将多个图像的深度信息与法向信息相结合,重建物体表面的三维形貌。

    光度立体技术的重点在于对物体表面法向量的计算。Lu等人[34]使用卷积神经网络(Convolutional Neural Network, CNN)预测物体的初始法线,在多光谱光度立体框架中进行迭代优化。近年来,国内外学者将双向反射分布函数(Bidirectional Reflectance Distribution Function, BRDF)应用于光度立体技术中,进行法向量求解,用于检测和分离高反射成分[35-36]

    国内外学者还对照明系统进行了重点研究,通过合理设计光源数量和角度,得到准确的法线信息。Pei等人[37]提出的测量系统由投影仪、工业相机和29个发光二极管构成。利用FPP测量得到待测物点云,通过近场光度立体系统得到像素的法向信息,将两者融合,进而得到完整的点云。测量结果与CMM测量数据相比,RMSE仅为0.0069 mm,峰谷值偏差为0.0094 mm。Meng等人[38]将照明系统放置在30个不同方向上,所得均方根误差为0.0053 mm。然而,由于系统测量角度约为±7°,限制了可以重建的曲面范围。

    光度立体技术具有多个照明光路,能够实现相互补盲,精度可达到0.01 mm之内,有助于实现HDR测量。但是,测量结果依赖于反射模型,对于加工过程中的零件,由于其表面的反射情况复杂及形状未知,对反射模型进行标定不具有普适性,容易带来测量误差,且测量系统较为复杂,测量场景具有较大的限制性,难以与加工机床和机械臂配合实现在机测量。

    基于硬件设备的HDR技术在硬件上对传统FPP技术进行了改进,使相机可以拍摄到较为清晰的条纹图像,从而降低了反光对测量精度的影响。表2分析了上述技术的测量精度以及有无额外硬件方面的对比情况。对于测量精度,本文仅对各参考文献中的精度进行了总结。由于不同实验时,所测工件尺寸、表面反射属性和测量环境不同,因此测量精度的比较不具有普遍性。

    表  2  基于硬件设备的HDR技术对比
    Table  2.  Comparison of HDR technologies based on hardware devices
    HDR技术 额外硬件 测量精度 参考文献
    相机曝光 MAE<0.1 mm [13, 17]
    偏振滤光片 偏振片 MAE<0.1 mm [23-24]
    相位偏折术 LCD显示屏 MAE<0.001 mm [27, 30-31]
    光度立体法 多个光源 MAE<0.01 mm [37-38]
    下载: 导出CSV 
    | 显示表格

    基于条纹算法的HDR结构光在机检测技术不依赖于硬件设备,主要是在算法层面改进条纹成像质量,从而提高测量精度。目前,国内外学者通过改变条纹图案亮度、颜色、频率以及更改条纹编码和解码等方式展开了深入研究。

    调整条纹图案强度是根据物体表面反射强度的不同,逐像素确定投影条纹的灰度,从而使整体图像具有较好的SNR且未达到饱和状态,也称为自适应条纹法。

    Waddington等人[39]首次提出自适应调整投影条纹图最大输入灰度的方法,并合成了不同强度的图像,以避免图像饱和。然而,机械加工环境大多属于暗房环境,仅改变条纹的最大灰度值虽可削弱金属工件反光的效果,但同时也降低了暗区亮度,导致在机检测精度降低。该团队还根据物体的局部反射率,将适当强度的条纹投影到物体的相应区域[40],但这需要进行复杂的预校准。Jiang等人[41]开发了一种自动选择参数(曝光时间、条纹颜色、条纹调制强度)的方法。他们通过调整不同相机曝光时间和不同强度的条纹,提高最终条纹图像的SNR,算法处理参数及效果如图9所示。

    图  9  Jiang所提方法的测量结果[41]
    Figure  9.  The measurement results by Jiang’s method[41]

    实验测量结果的MAE为0.003 mm,并且对环境光不敏感。但是此种方法需要6组条纹才能合成出一组SNR较好的条纹,基于三频四步相移的方法,双目相机需要拍摄144(2×6×12=144)张图像,需30 s完成,测量效率低且参数选择的程序较为复杂。Wang等人[42]在测量前先投影48幅图像,以建立投影与捕获图像之间的像素映射关系,然后投影了两幅图像灰度分别为0和100的图像,来估计反射率和背景强度。该方法的测量实验结果显示RMSE为0.109 mm。这种方法的测量效率依赖于硬件设备,为了满足在机检测的实时性需求。需要高速相机以及高刷新率的投影仪才能在短时间内完成映射关系。Sun等人[43]通过对投影仪像素建立投影强度模型,在高曝光和低曝光下分别拍摄两幅图像,以计算投影强度,再将亚像素坐标映射到投影仪图像上,并对其进行滤波,以实现精确的坐标对应。同时,提出了点云完整性(Point Cloud Integrality,PCI)作为评价量化指标。在采用三频四步相移方法的情况下,通过38幅图像获得了被测物的完整三维信息,高反射区域PCI指标达到94.1%。

    除了通过投影一系列条纹计算最佳强度外,国内外学者还提出利用多次迭代和多项式拟合等方法获得最优的投影强度。Babaie等人[44]使用多次迭代来计算适当的投影强度,并在小范围内测量有光泽的物体,但测量精度取决于迭代次数。Chen等人[45]通过拟合多项式函数将捕获的强度转换为投影强度,确定每个饱和像素簇对应的最佳投影强度。当簇中饱和像素对应的曲面区域反射率变化较大时,该函数将不准确。冯维等人[46]首先估计物体表面反射率、环境光和表面结构相互作用产生的反射光后,建立相机-投影仪间的强度关系式,以求解像素点的最佳投影灰度值。相较于传统方法,该方法的平均误差减少了61.9%,标准偏差减少了67.7%。李乾等人[47]依据待测元件表面的反射率分布情况,利用相机响应曲线逐像素计算元件上每一点对应的投影灰度值范围和最低投影灰度。实验RMSE值为0.103 mm,可完成金属材质高光工件的三维形貌测量。

    调整投影图案强度的优点在于无需额外添置硬件,算法容易集成,可以逐像素地进行灰度调整,避免了图像SNR整体降低的问题,因此测量精度较高,但由于需要投影多组条纹来确定最佳强度和映射关系,影响了测量效率,实时性较差。随着硬件设备的提升,可通过采用高刷新率的投影仪和相机改善这一问题,以满足在机检测的要求。

    通过颜色特征的镜面反射分离方法是由Shafer[48]提出的。该方法基于双色反射模型,对光线和颜色进行了讨论,认为像素点的颜色由物体表面的颜色和光源的颜色决定,由此分离出像素点的镜面反射部分。

    彩色图像通过拜耳滤波器可分离成R、G、B三个通道,不同通道的亮度衰减程度不同。对于单色照明,以蓝光为例,如图10(彩图见期刊电子版)所示,其量子效率存在以下关系eb>eg>er,高衰减红色通道可以接受非常高的亮度而不饱和,而低衰减蓝色通道则对低亮度更敏感。因此,国内外学者研究了彩色条纹和彩色相机,用于HDR测量。

    图  10  彩色图像不同通道亮度衰减原理图[51]
    Figure  10.  Schematic diagram of brightness attenuation in different channels of color image[51]

    Wang等人[49]在保持相机曝光度不变的前提下,通过投射7组不同颜色的条纹,从中选择灰度范围分布较大的4组条纹(蓝色、黄色、青色和白色)进行合成。该方法的缺点是对于投影仪的刷新率有较高的要求。Chua等人[50]使用RGB通道控制投影强度,来隔离不同颜色通道,将通道分为最高强度通道(蓝色通道)、中等强度通道(绿色通道)和最低强度通道(红色通道)。在每个通道下都可独立计算3D数据,达到增加测量动态范围的目的。

    使用不同颜色条纹进行多次投影时,需要高刷新率的投影仪才能保证测量效率,对硬件要求较高。因此,国内外学者还提出利用彩色相机拍摄单色条纹,再通过分离不同颜色通道实现HDR测量。Yin等人[51]从条纹图像中分离出R、G、G和B通道的4个单色子图像,校准R&G、G&B通道之间的衰减比后,用4个子图像合成一幅HDR的图像,从而避免了不良曝光的影响。实验分析相位平均误差(Phase Average Error, PAE)在0.03 rad内。Zheng等人[52]将绿色二值条纹投射到被测物上,使用彩色相机拍摄两张不同强度的条纹图像,进而分离出6张单通道图像,合成出1幅HDR图像后进行解码。Liu等人[53]将蓝色条纹投影到待测物上,并利用彩色相机采集图像。利用R、G、B三通道对光线的响应不同这一特性,将1幅条纹图像和1幅常规图像分成6幅具有不同亮度级别的单色图像。从每组图像中选择最亮但不饱和的对应像素生成最终的HDR图案,对其进行傅立叶解码,可以对具有光泽表面的动态物体进行测量,实验对撕A4纸的过程进行测量,测量效果如图11(彩图见期刊电子版)所示。

    图  11  Liu所提方法的动态物体测量结果[53]
    Figure  11.  The measurement results of dynamic objects by Liu’s method[53]

    基于颜色信息的HDR技术需要对不同颜色通道进行条纹编码,后续还需逐个通道对相位信息进行求解,再最终融合,因此算法编程复杂程度较高。但编程后算法运行较快,可在一定程度上抑制反光的影响,满足在机检测的需求。当被测的加工件具有复杂的纹理特征和多种颜色时,此类方法会受到限制,测量精度有所下降。

    除了对条纹图案亮度进行调整之外,国内外学者还通过对投影图案的编码和解码进行调整,从而防止高反射引起的相位丢失。

    Zhang等人[8]首次提出了混合质量的概念,以更全面地评估相位质量。在混合质量这一概念的指导下,对所有初始相位进行加权融合,以获得更精确的相位作为最终相位。Chen等人[54]提出了一种相位求解方法,将同一频率下存在饱和现象的帧去除,至少保留3张有效条纹,用余下条纹进行相位求解,其PAE结果可在0.01 rad内。Jiang等人[55]利用反向条纹和常规条纹的组合代替饱和条纹进行相位恢复,以减少相位误差,相位RMSE结果从0.29 rad减少到0.02 rad。然而,在估算饱和像素时,缺乏通用的公式。在Jiang方法基础上,Wang等人[56]将反转条纹和常规条纹的互补技术与广义相移算法相结合,投影规则条纹和反向条纹两组互补相移条纹。从两组条纹图中选择相同相机像素处的所有非饱和强度值,并使用广义相移算法检索相位,可以同时减少因饱和像素和伽马效应引起的误差。

    使用正弦条纹投射时,与条纹图像在饱和区域的灰度值达到255时,出现截断现象,不满足正弦规律,使该区域的相位信息丢失,引起相位求解误差。对于这一情况,在文献[57]中提到,当绝对相位范围在[0,2π)时,相位误差σΦ可以表示为:

    σ2Φ=2σ2Nf2B2,
    (4)

    式中,σ是高斯分布噪声的方差;N是相移步长;f是条纹频率;B是条纹调制。

    根据式(4)可知,在条纹调制固定的情况下,可以通过增加相移步长和条纹频率的方式降低相位误差。基于此原理,Chen等人[58]通过模拟分析得出以下结论:当相移步数N与条纹频率P满足整数倍关系时,饱和情况下的条纹也可恢复精确的相位。He等人[59]设计了一种棋盘式的高频条纹图案,在解码过程中,通过交叉位置而不是单个图像强度来提取编码信息。所提方法具有较好的鲁棒性,其测量精度在0.04 mm内。对于显微结构光三维测量方法,由于其景深较小,密集的高频条纹容易出现散焦,从而影响测量结果。针对这一问题,Hu等人[60]提出一种多频相移方案,采用相位融合的方式,用密度较低的条纹图像计算绝对相位,代替饱和区域中不可用的相位,从而降低相位误差的影响。所提方法可有效用于手表零件、封装的芯片等微小反光物体的测量,测量结果如图12(彩图见期刊电子版)所示。

    图  12  Hu所提方法对微小反光物体测量结果[60]
    Figure  12.  The measurement results of tiny reflective objects by Hu’s method[60]

    除此之外,国内外学者还曾通过投射格雷码图案,解决因高反光引起的解码误差。格雷码采用的是二值编码方式,具有较好的保边性,能很好地测量具有非均匀反射率和高亮区域表面的三维形貌[61]。Song等人[62]在传统格雷码的基础上,额外投影了一组黑白翻转的条纹图像。他们利用改进的边缘检测算法实现了亚像素定位,解决了因高反光导致条纹周期性模糊的问题。Lu等人[63]提出一种错位格雷码的编码策略。通过获取边界相互交错的编码图案,并通过降低单周期解码次序的子域解相阶数,获取正确的条纹阶次顺序。该方法有效降低了不均匀反射所导致的格雷码边界相位的阶跃误差。

    对于条纹图像编码和解码的方式,虽然无需复杂的算法,但是对于加工场景的复杂光线环境适用性较差。格雷码作为一种二值化编码方式所含信息量较少,且边缘部分易产生阶跃误差,不适合高精度的测量。此外,利用高频率正弦条纹和增加相移的方法,会使整体测量时间增加,不利于在机检测。

    近些年,随着人工智能的快速发展,国内外学者将深度学习等方法应用于结构光测量领域[64]。在图像预处理[65-66]、彩色条纹颜色通道分离[67-68]、高速动态测量[69]等方面进行了大量研究。深度学习的应用,将基于物理模型的传统思维转为基于数据模型驱动的智能化思维,为HDR结构光在机测量提供了新的解决方案。

    Liu等人[70]提出了一种基于支持向量机算法的智能条纹投影技术。通过训练不同表面的曝光时间和饱和像素数,仅将一个灰度模式投影到测试表面上,即可预测合适的曝光时间。所提方法用于电子束增材制造(Electron Beam Additive Manufacturing, EBAM)的在机检测中(如图13所示),该方法可以有效测量高反光曲面的金属粉末熔化状态和表面缺陷,并将测量结果反馈到制造过程中,从而提高了加工质量。

    图  13  EBAM机床结构光在机检测设备[70]
    Figure  13.  On-machine detection equipment with structured light in EBAM machine tool[70]

    彭广泽等人[71]基于CNN对条纹高光图像进行修复。将两个曝光条件下的条纹图像进行融合,确定迭代修复算法的初值,然后通过所提算法对局部高光区域进行快速修复。Yang等人[72]设计了一种基于条纹调制的“检测-修复”网络结构(如图14所示)。该结构由低调制区检测模块和条纹增强模块组成。通过低调制区检测模块准确分割出低调制区域,然后引入到条纹增强模块中,再将条纹图的高阶特征和低阶特征进行融合,从而可以预测低调制区的条纹分布。两个模块的结合可以恢复条纹图的饱和区,增强暗区的条纹信息。对量块的测量结果表明,采用该结构可以使RMSE结果由0.55 mm降至0.06 mm。

    图  14  高动态范围条纹图像改善网络[72]
    Figure  14.  High dynamic range fringe pattern improvement network[72]

    上述文献是通过深度学习的方法对高动态条纹图像进行修复,以消除高反射和低亮度的影响,从而根据修复的条纹图像求解准确的相位,得到三维模型。除此类方法之外,国内外学者还提出通过直接预测相位的方式实现HDR测量。Qiao等人[73]将深度学习应用于PMD技术中,从单个条纹图中提取镜面的相位信息,该方法用到两个神经网络(CNN1和CNN2)。CNN1用于预测输入条纹图像的背景强度,再将原条纹和背景强度输入到CNN2中,预测相位主值的分子和分母。Zhang等人[74]利用深度学习强大的计算能力,基于CNN预测相位主值,消除HDR引起的相位误差。训练的数据集为2880个条纹图像。所提方法不仅适用于高反射率的被测物,而且对低反射率物体也适用。RMSE测量结果为0.0578 mm,与三步相移法相比降低了53.5%,同时满足高速动态测量要求。

    通过深度学习方法解决HDR测量,可以提高测量效率,并在一定程度上提高测量精度,可用于加工时的在机检测。但是,相关算法的复杂程度较高,建立数据集和神经网络模型需要大量时间。而且,前期准备工作量较大,且测量结果过于依赖训练模型的种类和数量;训练模型少会导致测量精度降低;当训练模型足够多时,又会增加训练时间以及其他成本。

    根据算法思路、算法复杂程度和测量精度,对基于条纹算法的HDR技术对比列于表3。由于不同学者采用的精度指标有所差异,根据不同文献从平均偏差、相位偏差、均方根偏差方面对精度范围进行了总结。

    表  3  基于条纹算法的HDR技术对比
    Table  3.  Comparison of HDR technologies based on fringe algorithm
    HDR技术 算法思路 算法复杂程度 测量精度 参考文献
    调整条纹强度 逐像素改变条纹灰度值 MAE<0.03 mm [39, 44-45]
    颜色信息 颜色通道分离 PAE<0.03 rad [48, 49]
    图案编码、解码 增加相移步数、多频条纹 PAE<0.02 rad [52, 55, 57]
    智能算法 神经网络处理高光图像、预测相位 RMSE<0.06 mm [70, 72]
    下载: 导出CSV 
    | 显示表格

    上述已对HDR结构光测量技术从硬件设备和条纹算法两方面进行了分类和综述。表4结合加工零件在机检测的关键问题,如光线条件适应性、系统硬件设备和检测效率,对技术特点进行总结。

    表  4  各类HDR测量技术总结
    Table  4.  Summary of various HDR measurement technologies
    HDR技术 优点 缺点 光线条件
    适应性
    系统硬件设备 检测效率 加工在机检测
    相机曝光 无需添置额外硬件、无后续其他处理。 选择曝光时间具有一定盲目性,需多次测量合成最优数据。 单次曝光适应性较差,多重曝光适应性好。 简单 不适用
    偏振滤光片 额外硬件较为简单、无其他复杂算法。 单偏振通道易降低整体图像的SNR,使用多个偏振通道时,需多次调整偏振片角度合成最优数据。 单通道适应性较差,多通道适应性好。 单通道简单,多通道较复杂。 不适用
    相位偏折术 适用于类镜面物体的测量,测量精度高,无其他复杂算法 空间摆放位置受限制,不适用于金属等反光件。 简单 适用(镜面、类镜面工件)
    光度立体法 利用多照明系统实现视角补盲 建立的反射模型不具有普适性。 复杂 不适用
    调整条纹强度 逐像素调整图像亮度,条纹图像具有较高的SNR 对于场景和反射区域的标定需投射多组条纹确定映射关系,算法的效率需依靠投影仪的帧率决定。 简单
    (配合高速投影)
    适用(配合高速投影)
    颜色信息 算法简单,无其他复杂算法 对于带有颜色和纹理特征的被测物,测量精度会受到影响。 简单 适用
    图案编码、解码 算法简单 增加条纹频率和相移步数,影响了测量的效率,且测量精度较低。 简单 不适用
    智能算法 测量效率高,可以实现动态测量 算法复杂,成本较高,需要高度定制的训练样本。 简单 适用
    下载: 导出CSV 
    | 显示表格

    将HDR结构光测量技术应用于在机检测中,与传统的FPP技术相比,该技术更为灵活,可实现对复杂场景的检测。应用于机械加工、航空航天、汽车装备等领域时,可提高作业效率,推动工业智能化,具有较高的应用价值。

    (1)智能制造加工-检测一体化

    机械制造业是工业发展的基础,将结构光设备与机床配合,应用于工件的切削、抛光、轮廓识别、装配等工艺中,实现加工-检测一体化,避免了零件反复拆卸产生变形和定位误差,从而提高了加工效率和精度。HDR技术可有效抑制反光带来的不利影响,有助于判断刀具磨损、零件过切、欠切和表面质量,实现超精密加工和无人化智能工厂。

    (2)极端尺寸工件的精密测量

    对极大和极小零件进行精密测量是仪器测量领域的难点。对于飞机机翼、汽车车身、汽轮机叶片等极大尺寸的待测件,无需将零件拆卸,可通过结构光设备与机械臂配合,进行路径规划后逐视点检测,所测数据再通过点云拼接得到完整三维模型;而对于芯片、光学镜片等毫米级尺寸的待测件,可将FPP技术与显微镜结合,实现小视场高精度的测量[75]。HDR技术可有效提高测量精度,为零件缺陷检测及维护、增材制造提供准确的三维数据。

    (3)医工结合领域的精密测量

    FPP技术广泛应用于生物医学领域,为人脸识别、医学诊断提供了新的技术支持。相关学者采用双目结构光自适应光栅算法对术中肝脏表面进行三维重建,并对肿瘤进行精确定位[76]。使用HDR条纹投影技术对口腔、骨骼等部位进行三维重建,可直接应用于临床诊断和手术中,无需指定暗房场地进行扫描,提高了诊断效率和建模精度,为后续治疗提供了新的途径。

    (4)大口径镜面、类镜面物体精密测量

    天文望远镜镜片、反射镜镜片等大口径镜面、类镜面物体在军用、民用、天文等领域起到重要作用,其加工主要包括铣磨成型、研磨和抛光三个阶段,每道工序完成后,都需对表面进行精密测量,判断是否满足加工要求。研磨阶段后常采用轮廓检测法进行面型测量,抛光后采用光学干涉法进行测量。对于从研磨到抛光的过渡阶段,轮廓检测仪受到精度限制,且测量时间较长;而光学干涉仪动态范围小,难以覆盖全口径检测,两种方法都具有局限性。此时,考虑到条纹反射法具有成本低、精度高、测量速度快的特点,可有效衔接轮廓检测和光学干涉检测的面形误差测量范围,实现对大口径镜面、类镜面物体的全口径在机检测。

    (1)复杂属性及复杂光照环境的适应性

    复杂属性表面和复杂光照环境对HDR结构光技术的影响是具有重要研究意义和挑战性的课题。在未来技术研究中,针对被测物尺寸(极大、极小)、表面属性(如纯镜面反射、镜面反射与漫反射结合、大曲率镜面以及透明属性等)和测量环境(日光环境、暗房环境、环境光影响等)进行参数选择,使所研究技术具有更高的适用性是重要突破方向。

    (2)使用机械臂进行视点、路径规划

    结构光设备与机械臂配合是实现在机检测的重要条件,对待测件进行合理的视点、路径规划,可提高检测效率。在已知CAD模型的条件下,需要综合考虑视点数量和测量精度,要在保证测量精度的前提下保证测量效率。在未知CAD模型的条件下,需要准确预判下一视点范围,保证视点选择的合理性。同时,如何在机械臂运动时避免机床颤振、机械臂抖动对测量精度的影响,也是近些年的研究热点。

    (3)深度学习智能算法的应用

    深度学习技术的成熟也为工业领域提供了新的解决方案。在HDR结构光技术中,将深度学习应用于位置标定、图像处理和相位预测,对于数据集的采集、训练机制和优化等方面都是重点研究对象。目前,此类方法处于研究初期,虽然给结构光测量技术带来了实质性的改进,但同时也需评估在特定任务中使用大量数据和计算资源的经济性与可靠性。

    (4)动态测量与在机检测结合

    结合结构光系统的测量原理,目前的在机加工测量更多是针对某道工序结束后,暂停加工,在短时间内进行形貌等质量检测。使用现有的商用硬件设备进行在机检测时,常采用相移法模式,且对于具有高反光特性的金属零件需在表面喷涂层,影响了加工效率、自动化和智能化程度,实时性较差,不利于及时反馈测量结果以及对加工参数进行补偿修正。将动态测量应用于在机检测中,可边加工边检测,克服了现有力传感器等的安装尺寸的局限,有助于进一步提高效率并及时修正加工参数。使用傅立叶轮廓术(Fourier Transform Profilometry, FTP),可通过单张条纹图像求解得到三维模型,但是测量精度较低。在未来的研究中,可将高速相机、高速投影仪的硬件设备与深度学习、HDR技术、FTP技术等算法相结合。在完成动态测量的同时,进一步提高测量精度。

    近年来,FPP测量技术广泛应用于航空发动机叶片、汽轮机叶片等复杂曲面零件的在机检测中,是以满足国家需求为导向的先进制造技术。针对金属零件加工时,所处环境光线复杂且自身具有高反光的问题,采用HDR技术与FPP在机测量技术结合,可有效抑制相关影响。

    本文首先介绍了条纹结构光的测量原理,总结出HDR结构光在机检测面临的重点问题。其次,回顾了国内外学者近几年对HDR结构光技术的研究进展,将其综述为基于硬件设备的HDR技术和基于条纹算法的HDR技术两大类。然后,结合在机检测的条件需求,对各类技术进行总结,比较不同方法的优缺点和在机检测的适用性,为不同技术的应用提供一定的指导。最后,结合近年来工业制造和精密测量的研究热点,进行潜在应用分析,并提出技术展望。HDR结构光在机检测技术是一种极具研究价值的测量技术,对智能制造、精密测量的发展起到推动作用。在未来,进一步提高HDR结构光在机检测技术的测量效率和精度,具有更高的研究价值和发展潜力。

  • 图 1  结构光测量原理示意图[8]

    Figure 1.  Schematic diagram of monocular structured light measurement principle[8]

    图 2  结构光在机检测示意图

    Figure 2.  Schematic diagram of structured light on-machine detection

    图 3  结构光在机检测系统架构

    Figure 3.  Architecture of on-machine detection system of structured light

    图 4  曝光融合算法流程[17]

    Figure 4.  The process of exposure fusion algorithm[17]

    图 5  Xiang设计的偏振测量系统[24]

    Figure 5.  Polarization measurement system designed by Xiang[24]

    图 6  封闭激光腔内COPMD测量系统[26]

    Figure 6.  COPMD measurement system in an enclosed laser cavity[26]

    图 7  条纹反射在位面形检测系统中的装置[30]

    Figure 7.  Fringe reflection setup in on-machine surface measurement system[30]

    图 8  大口径抛物面反射镜实验装置图[31]

    Figure 8.  Experimental device diagram of large diameter parabolic reflector[31]

    图 9  Jiang所提方法的测量结果[41]

    Figure 9.  The measurement results by Jiang’s method[41]

    图 10  彩色图像不同通道亮度衰减原理图[51]

    Figure 10.  Schematic diagram of brightness attenuation in different channels of color image[51]

    图 11  Liu所提方法的动态物体测量结果[53]

    Figure 11.  The measurement results of dynamic objects by Liu’s method[53]

    图 12  Hu所提方法对微小反光物体测量结果[60]

    Figure 12.  The measurement results of tiny reflective objects by Hu’s method[60]

    图 13  EBAM机床结构光在机检测设备[70]

    Figure 13.  On-machine detection equipment with structured light in EBAM machine tool[70]

    图 14  高动态范围条纹图像改善网络[72]

    Figure 14.  High dynamic range fringe pattern improvement network[72]

    表  1  三维视觉测量技术分类

    Table  1.   Classification of three-dimensional visual measurement technology

    视觉测
    量分类
    是否投
    射光源
    具体分类 特点
    被动视
    觉测量
    单目视觉测量 基于图像聚焦程度完成三维重建,多用于显微视觉测量中。
    双目视觉测量 根据三角测量原理实现三维重建,应用于双目立体摄像头。
    多目视觉测量 增加辅助相机,通过光束平差提高测量精度。
    主动视
    觉测量
    点扫描式 激光器投射光点,根据光标中心坐标和标定数据进行重建,测量效率低。
    线扫描式 激光器投射光条代替光点,提高效率,广泛应用于激光扫描仪中。
    面扫描式 通过投影仪投射二维结构光,单次投射覆盖区域大,测量效率最高。
    下载: 导出CSV

    表  2  基于硬件设备的HDR技术对比

    Table  2.   Comparison of HDR technologies based on hardware devices

    HDR技术 额外硬件 测量精度 参考文献
    相机曝光 MAE<0.1 mm [13, 17]
    偏振滤光片 偏振片 MAE<0.1 mm [23-24]
    相位偏折术 LCD显示屏 MAE<0.001 mm [27, 30-31]
    光度立体法 多个光源 MAE<0.01 mm [37-38]
    下载: 导出CSV

    表  3  基于条纹算法的HDR技术对比

    Table  3.   Comparison of HDR technologies based on fringe algorithm

    HDR技术 算法思路 算法复杂程度 测量精度 参考文献
    调整条纹强度 逐像素改变条纹灰度值 MAE<0.03 mm [39, 44-45]
    颜色信息 颜色通道分离 PAE<0.03 rad [48, 49]
    图案编码、解码 增加相移步数、多频条纹 PAE<0.02 rad [52, 55, 57]
    智能算法 神经网络处理高光图像、预测相位 RMSE<0.06 mm [70, 72]
    下载: 导出CSV

    表  4  各类HDR测量技术总结

    Table  4.   Summary of various HDR measurement technologies

    HDR技术 优点 缺点 光线条件
    适应性
    系统硬件设备 检测效率 加工在机检测
    相机曝光 无需添置额外硬件、无后续其他处理。 选择曝光时间具有一定盲目性,需多次测量合成最优数据。 单次曝光适应性较差,多重曝光适应性好。 简单 不适用
    偏振滤光片 额外硬件较为简单、无其他复杂算法。 单偏振通道易降低整体图像的SNR,使用多个偏振通道时,需多次调整偏振片角度合成最优数据。 单通道适应性较差,多通道适应性好。 单通道简单,多通道较复杂。 不适用
    相位偏折术 适用于类镜面物体的测量,测量精度高,无其他复杂算法 空间摆放位置受限制,不适用于金属等反光件。 简单 适用(镜面、类镜面工件)
    光度立体法 利用多照明系统实现视角补盲 建立的反射模型不具有普适性。 复杂 不适用
    调整条纹强度 逐像素调整图像亮度,条纹图像具有较高的SNR 对于场景和反射区域的标定需投射多组条纹确定映射关系,算法的效率需依靠投影仪的帧率决定。 简单
    (配合高速投影)
    适用(配合高速投影)
    颜色信息 算法简单,无其他复杂算法 对于带有颜色和纹理特征的被测物,测量精度会受到影响。 简单 适用
    图案编码、解码 算法简单 增加条纹频率和相移步数,影响了测量的效率,且测量精度较低。 简单 不适用
    智能算法 测量效率高,可以实现动态测量 算法复杂,成本较高,需要高度定制的训练样本。 简单 适用
    下载: 导出CSV
  • [1] 李文龙, 李中伟, 毛金城. iPoint3D曲面检测软件开发与工程应用综述[J]. 机械工程学报,2020,56(7):127-150. doi: 10.3901/JME.2020.07.127

    LI W L, LI ZH W, MAO J CH. The development and application review of iPoint3D software for surface inspection[J]. Journal of Mechanical Engineering, 2020, 56(7): 127-150. (in Chinese) doi: 10.3901/JME.2020.07.127
    [2] 李茂月, 刘泽隆, 赵伟翔, 等. 面结构光在机检测的叶片反光抑制技术[J]. 中国光学,2022,15(3):464-475. doi: 10.37188/CO.2021-0194

    LI M Y, LIU Z L, ZHAO W X, et al. Blade reflection suppression technology based on surface structured light on-machine detection[J]. Chinese Optics, 2022, 15(3): 464-475. (in Chinese) doi: 10.37188/CO.2021-0194
    [3] 任明阳, 王立忠, 赵建博, 等. 复杂曲面零件面结构光扫描视点规划[J]. 中国光学(中英文),2023,16(1):113-126. doi: 10.37188/CO.2022-0026

    REN M Y, WANG L ZH, ZHAO J B, et al. Viewpoint planning of surface structured light scanning for complex surface parts[J]. Chinese Optics, 2023, 16(1): 113-126. (in Chinese) doi: 10.37188/CO.2022-0026
    [4] 张宗华, 刘巍, 刘国栋, 等. 三维视觉测量技术及应用进展[J]. 中国图象图形学报,2021,26(6):1483-1502.

    ZHANG Z H, LIU W, LIU G D, et al. Overview of the development and application of 3D vision measurement technology[J]. Journal of Image and Graphics, 2021, 26(6): 1483-1502. (in Chinese)
    [5] XU J, ZHANG S. Status, challenges, and future perspectives of fringe projection profilometry[J]. Optics and Lasers in Engineering, 2020, 135: 106193. doi: 10.1016/j.optlaseng.2020.106193
    [6] 张宗华, 于瑾, 高楠, 等. 高反光表面三维形貌测量技术[J]. 红外与激光工程,2020,49(3):0303006. doi: 10.3788/IRLA202049.0303006

    ZHANG Z H, YU J, GAO N, et al. Three-dimensional shape measurement techniques of shiny surfaces[J]. Infrared and Laser Engineering, 2020, 49(3): 0303006. (in Chinese) doi: 10.3788/IRLA202049.0303006
    [7] FENG S J, ZHANG L, ZUO CH, et al. High dynamic range 3D measurements with fringe projection profilometry: a review[J]. Measurement Science and Technology, 2018, 29(12): 122001. doi: 10.1088/1361-6501/aae4fb
    [8] ZHANG P, ZHONG K, LI ZH W, et al. Hybrid-quality-guided phase fusion model for high dynamic range 3D surface measurement by structured light technology[J]. Optics Express, 2022, 30(9): 14600-14614. doi: 10.1364/OE.457305
    [9] Zhang S. Rapid and automatic optimal exposure control for digital fringe projection technique[J]. Optics and Lasers in Engineering, 2020, 128: 106029. doi: 10.1016/j.optlaseng.2020.106029
    [10] 马泽龙, 高慧斌, 余毅, 等. 采用图像直方图特征函数的高速相机自动曝光方法[J]. 光学 精密工程,2017,25(4):1026-1035. doi: 10.3788/OPE.20172504.1026

    MA Z L, GAO H B, YU Y, et al. Auto exposure control for high frame rate camera using image histogram feature function[J]. Optics and Precision Engineering, 2017, 25(4): 1026-1035. (in Chinese) doi: 10.3788/OPE.20172504.1026
    [11] 雷经发, 陆宗胜, 李永玲, 等. 基于投影栅相位法和多曝光图像融合技术的强反射表面轮廓检测[J]. 光学 精密工程,2022,30(18):2195-2204. doi: 10.37188/OPE.20223018.2195

    LEI J F, LU Z SH, LI Y L, et al. High reflection surface topography measurement based on fringe projection phase method and multi-exposure image fusion technology[J]. Optics and Precision Engineering, 2022, 30(18): 2195-2204. (in Chinese) doi: 10.37188/OPE.20223018.2195
    [12] ZHANG S, YAU S T. High dynamic range scanning technique[J]. Optical Engineering, 2009, 48(3): 033604. doi: 10.1117/1.3099720
    [13] SONG ZH, JIANG H L, LIN H B, et al. A high dynamic range structured light means for the 3D measurement of specular surface[J]. Optics and Lasers in Engineering, 2017, 95: 8-16. doi: 10.1016/j.optlaseng.2017.03.008
    [14] FENG SH J, ZHANG Y ZH, CHEN Q, et al. General solution for high dynamic range three-dimensional shape measurement using the fringe projection technique[J]. Optics and Lasers in Engineering, 2014, 59: 56-71. doi: 10.1016/j.optlaseng.2014.03.003
    [15] CUI H H, LI ZH J, TIAN W, et al. Multiple-exposure adaptive selection algorithm for high dynamic range 3D fringe projection measurement[J]. Proceedings of SPIE, 2019, 11053: 110530M.
    [16] RAO L, DA F P. High dynamic range 3D shape determination based on automatic exposure selection[J]. Journal of Visual Communication and Image Representation, 2018, 50: 217-226. doi: 10.1016/j.jvcir.2017.12.003
    [17] WU K, TAN J, XIA H L, et al. An exposure fusion-based structured light approach for the 3D measurement of a specular surface[J]. IEEE Sensors Journal, 2021, 21(5): 6314-6324. doi: 10.1109/JSEN.2020.3027317
    [18] SALAHIEH B, CHEN ZH Y, RODRIGUEZ J J, et al. Multi-polarization fringe projection imaging for high dynamic range objects[J]. Optics Express, 2014, 22(8): 10064-10071. doi: 10.1364/OE.22.010064
    [19] 平茜茜, 刘勇, 董欣明, 等. 基于偏振双目视觉的无纹理高反光目标三维重构[J]. 红外与毫米波学报,2017,36(4):432-438.

    PING X X, LIU Y, DONG X M, et al. 3-D reconstruction of textureless and high-reflective target by polarization and binocular stereo vision[J]. Journal of Infrared and Millimeter Waves, 2017, 36(4): 432-438. (in Chinese)
    [20] 郝婧蕾, 赵永强, 赵海盟, 等. 偏振多光谱机器视觉的高反光无纹理目标三维重构方法[J]. 测绘学报,2018,47(6):816-824.

    HAO J L, ZHAO Y Q, ZHAO H M, et al. 3D Reconstruction of High-reflective and Textureless Targets Based on Multispectral Polarization and Machine Vision[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(6): 816-824. (in Chinese)
    [21] WANG Y H, ZHANG Q, HU Y, et al. Rapid 3D measurement of high dynamic range surface based on multi-polarization fringe projection[J]. Optical Engineering, 2021, 60(8): 084107.
    [22] ZHU ZH M, ZHU W T, ZHOU F Q, et al. Three-dimensional measurement of fringe projection based on the camera response function of the polarization system[J]. Optical Engineering, 2021, 60(5): 055105.
    [23] MAEDA Y, SHIBATA S, HAGEN N, et al. Single shot 3D profilometry by polarization pattern projection[J]. Applied Optics, 2020, 59(6): 1654-1659. doi: 10.1364/AO.382690
    [24] XIANG G, ZHU H J, GUO H W. Spatial phase-shifting profilometry by use of polarization for measuring 3D shapes of metal objects[J]. Optics Express, 2021, 29(13): 20981-20994. doi: 10.1364/OE.427407
    [25] 王月敏, 张宗华, 高楠. 基于全场条纹反射的镜面物体三维面形测量综述[J]. 光学 精密工程,2018,26(5):1014-1027. doi: 10.3788/OPE.20182605.1014

    WANG Y M, ZHANG Z H, GAO N. Review on three-dimensional surface measurements of specular objects based on full-field fringe reflection[J]. Optics and Precision Engineering, 2018, 26(5): 1014-1027. (in Chinese) doi: 10.3788/OPE.20182605.1014
    [26] ZHUANG Y CH, ZHENG Y M, LIN SH B, et al. Surface shape distortion online measurement method for compact laser cavities based on phase measuring Deflectometry[J]. Photonics, 2022, 9(3): 151. doi: 10.3390/photonics9030151
    [27] GAO F, XU Y J, JIANG X Q. Near optical coaxial phase measuring deflectometry for measuring structured specular surfaces[J]. Optics Express, 2022, 30(10): 17554-17566. doi: 10.1364/OE.457198
    [28] HAN H, WU SH Q, SONG ZH. Curved LCD based deflectometry method for specular surface measurement[J]. Optics and Lasers in Engineering, 2022, 151: 106909. doi: 10.1016/j.optlaseng.2021.106909
    [29] SU P, PARKS R E, WANG L R, et al. Software configurable optical test system: a computerized reverse Hartmann test[J]. Applied Optics, 2010, 49(23): 4404-4412. doi: 10.1364/AO.49.004404
    [30] 邵山川, 陶小平, 王孝坤. 基于条纹反射的超精密车削反射镜的在位面形检测[J]. 激光与光电子学进展,2018,55(7):071203.

    SHAO SH CH, TAO X P, WANG X K. On-machine surface shape measurement of reflective mirrors by ultra-precision turning based on fringe reflection[J]. Laser & Optoelectronics Progress, 2018, 55(7): 071203. (in Chinese)
    [31] 袁婷. 基于条纹反射法的大口径非球面反射镜面形检测技术研究[D]. 长春: 中国科学院研究生院(长春光学精密机械与物理研究所), 2016.

    YUAN T. Study on fringe-reflection optical surface shape measurement technology for large aspheric mirror[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, CAS, 2016. (in Chinese)
    [32] OH C J, LOWMAN A E, SMITH G A, et al. Fabrication and testing of 4.2m off-axis aspheric primary mirror of Daniel K. Inouye Solar Telescope[J]. Proceedings of SPIE, 2016, 9912: 99120O.
    [33] WOODHAM R J. Photometric method for determining surface orientation from multiple images[J]. Optical Engineering, 1980, 19(1): 191139.
    [34] LU L, QI L, LUO Y S, et al. Three-dimensional reconstruction from single image base on combination of CNN and multi-spectral photometric stereo[J]. Sensors, 2018, 18(3): 764. doi: 10.3390/s18030764
    [35] 张颖, 李金龙, 黄趾维, 等. 基于BRDF模型的金属表面反射特性及相变特性研究[J]. 光电技术应用,2017,32(3):32-35. doi: 10.3969/j.issn.1673-1255.2017.03.008

    ZHANG Y, LI J L, HUANG ZH W, et al. Research on reflection and phase shift characters of metal surface based on BRDF model[J]. Electro-Optic Technology Application, 2017, 32(3): 32-35. (in Chinese) doi: 10.3969/j.issn.1673-1255.2017.03.008
    [36] 王金海, 李华, 魏力. 基于C-T模型的光学元件加工表面的光学特性研究[J]. 光学技术,2021,47(2):172-177.

    WANG J H, LI H, WEI L. Study on optical properties of machining surface of optical element based on C-T model[J]. Optical Technique, 2021, 47(2): 172-177.
    [37] PEI X H, REN M J, WANG X, et al. Profile measurement of non-Lambertian surfaces by integrating fringe projection profilometry with near-field photometric stereo[J]. Measurement, 2022, 187: 110277. doi: 10.1016/j.measurement.2021.110277
    [38] MENG L F, LU L Y, BEDARD N, et al. Single-shot specular surface reconstruction with gonio-plenoptic imaging[C]. IEEE International Conference on Computer Vision, IEEE, 2015: 3433-3441.
    [39] WADDINGTON C, KOFMAN J. Analysis of measurement sensitivity to illuminance and fringe-pattern gray levels for fringe-pattern projection adaptive to ambient lighting[J]. Optics and Lasers in Engineering, 2010, 48(2): 251-256. doi: 10.1016/j.optlaseng.2009.07.001
    [40] LI D, KOFMAN J. Adaptive fringe-pattern projection for image saturation avoidance in 3D surface-shape measurement[J]. Optics Express, 2014, 22(8): 9887-9901. doi: 10.1364/OE.22.009887
    [41] JIANG H ZH, ZHAO H J, LI X D. High dynamic range fringe acquisition: a novel 3-D scanning technique for high-reflective surfaces[J]. Optics and Lasers in Engineering, 2012, 50(10): 1484-1493. doi: 10.1016/j.optlaseng.2011.11.021
    [42] WANG J H, YANG Y X, ZHOU Y G. 3-D shape reconstruction of non-uniform reflectance surface based on pixel intensity, pixel color and camera exposure time adaptive adjustment[J]. Scientific Reports, 2021, 11(1): 4700. doi: 10.1038/s41598-021-83779-9
    [43] SUN J H, ZHANG Q Y. A 3D shape measurement method for high-reflective surface based on accurate adaptive fringe projection[J]. Optics and Lasers in Engineering, 2022, 153: 106994. doi: 10.1016/j.optlaseng.2022.106994
    [44] BABAIE G, ABOLBASHARI M, FARAHI F. Dynamics range enhancement in digital fringe projection technique[J]. Precision Engineering, 2015, 39: 243-251. doi: 10.1016/j.precisioneng.2014.06.007
    [45] CHEN CH, GAO N, WANG X J, et al. Adaptive pixel-to-pixel projection intensity adjustment for measuring a shiny surface using orthogonal color fringe pattern projection[J]. Measurement Science and Technology, 2018, 29(5): 055203. doi: 10.1088/1361-6501/aab07a
    [46] 冯维, 徐仕楠, 王恒辉, 等. 逐像素调制的高反光表面三维测量方法[J]. 中国光学,2022,15(3):488-497. doi: 10.37188/CO.2021-0220

    FENG W, XU SH N, WANG H H, et al. Three-dimensional measurement method of highly reflective surface based on per-pixel modulation[J]. Chinese Optics, 2022, 15(3): 488-497. (in Chinese) doi: 10.37188/CO.2021-0220
    [47] 李乾, 薛俊鹏, 张启灿, 等. 利用相机响应曲线实现高反光元件三维面形测量[J]. 光学学报,2022,42(7):0712001. doi: 10.3788/AOS202242.0712001

    LI Q, XUE J P, ZHANG Q C, et al. Three dimensional shape measurement of high reflective elements using camera response curve[J]. Acta Optica Sinica, 2022, 42(7): 0712001. (in Chinese) doi: 10.3788/AOS202242.0712001
    [48] SHAFER S A. Using color to separate reflection components[J]. Color Research & Application, 1985, 10(4): 210-218.
    [49] WANG J H, YANG Y X. High-speed three-dimensional measurement technique for object surface with a large range of reflectivity variations[J]. Applied Optics, 2018, 57(30): 9172-9182. doi: 10.1364/AO.57.009172
    [50] CHUA S Y, LIM C C, ENG S K, et al. Improved high dynamic range for 3D shape measurement based on saturation of the coloured fringe[J]. Pertanika Journal of Science & Technology, 2021, 29(2): 759-770.
    [51] YIN Y K, CAI Z W, JIANG H, et al. High dynamic range imaging for fringe projection profilometry with single-shot raw data of the color camera[J]. Optics and Lasers in Engineering, 2017, 89: 138-144. doi: 10.1016/j.optlaseng.2016.08.019
    [52] ZHENG Y, WANG Y J, SURESH V, et al. Real-time high-dynamic-range fringe acquisition for 3D shape measurement with a RGB camera[J]. Measurement Science and Technology, 2019, 30(7): 075202. doi: 10.1088/1361-6501/ab0ced
    [53] LIU Y ZH, FU Y J, ZHUAN Y H, et al. High dynamic range real-time 3D measurement based on Fourier transform profilometry[J]. Optics & Laser Technology, 2021, 138: 106833.
    [54] CHEN Y M, HE Y M, HU E Y. Phase deviation analysis and phase retrieval for partial intensity saturation in phase-shifting projected fringe profilometry[J]. Optics Communications, 2008, 281(11): 3087-3090. doi: 10.1016/j.optcom.2008.01.070
    [55] JIANG C F, BELL T, ZHANG S. High dynamic range real-time 3D shape measurement[J]. Optics Express, 2016, 24(7): 7337-7346. doi: 10.1364/OE.24.007337
    [56] WANG M M, DU G L, ZHOU C L, et al. Enhanced high dynamic range 3D shape measurement based on generalized phase-shifting algorithm[J]. Optics Communications, 2017, 385: 43-53. doi: 10.1016/j.optcom.2016.10.023
    [57] ZUO CH, HUANG L, ZHANG M L, et al. Temporal phase unwrapping algorithms for fringe projection profilometry: A comparative review[J]. Optics and Lasers in Engineering, 2016, 85: 84-103. doi: 10.1016/j.optlaseng.2016.04.022
    [58] CHEN B, ZHANG S. High-quality 3D shape measurement using saturated fringe patterns[J]. Optics and Lasers in Engineering, 2016, 87: 83-89. doi: 10.1016/j.optlaseng.2016.04.012
    [59] HE Z X, LI P L, ZHAO X Y, et al. Chessboard-like high-frequency patterns for 3D measurement of reflective surface[J]. IEEE Transactions on Instrumentation and Measurement, 2021, 70: 5009712.
    [60] HU Y, CHEN Q, LIANG Y CH, et al. Microscopic 3D measurement of shiny surfaces based on a multi-frequency phase-shifting scheme[J]. Optics and Lasers in Engineering, 2019, 122: 1-7. doi: 10.1016/j.optlaseng.2019.05.019
    [61] 张启灿, 吴周杰. 基于格雷码图案投影的结构光三维成像技术[J]. 红外与激光工程,2020,49(3):0303004. doi: 10.3788/IRLA202049.0303004

    ZHANG Q C, WU ZH J. Three-dimensional imaging technique based on Gray-coded structured illumination[J]. Infrared and Laser Engineering, 2020, 49(3): 0303004. (in Chinese) doi: 10.3788/IRLA202049.0303004
    [62] SONG ZH, CHUNG R, ZHANG X T. An accurate and robust strip-edge-based structured light means for shiny surface micromeasurement in 3-D[J]. IEEE Transactions on Industrial Electronics, 2013, 60(3): 1023-1032. doi: 10.1109/TIE.2012.2188875
    [63] LU L L, WU ZH J, ZHANG Q C, et al. High-efficiency dynamic three-dimensional shape measurement based on misaligned Gray-code light[J]. Optics and Lasers in Engineering, 2022, 150: 106873. doi: 10.1016/j.optlaseng.2021.106873
    [64] ZUO CH, QIAN J M, FENG SH J, et al. Deep learning in optical metrology: a review[J]. Light:Science & Applications, 2022, 11(1): 39.
    [65] REYES-FIGUEROA A, FLORES V H, RIVERA M. Deep neural network for fringe pattern filtering and normalization[J]. Applied Optics, 2021, 60(7): 2022-2036. doi: 10.1364/AO.413404
    [66] JEON W, JEONG W, SON K, et al. Speckle noise reduction for digital holographic images using multi-scale convolutional neural networks[J]. Optics Letters, 2018, 43(17): 4240-4243. doi: 10.1364/OL.43.004240
    [67] ZHANG B W, LIN SH N, LIN J Y, et al. Single-shot high-precision 3D reconstruction with color fringe projection profilometry based BP neural network[J]. Optics Communications, 2022, 517: 128323. doi: 10.1016/j.optcom.2022.128323
    [68] NGUYEN H, WANG ZH Y. Accurate 3D shape reconstruction from single structured-light image via fringe-to-fringe network[J]. Photonics, 2021, 8(11): 459. doi: 10.3390/photonics8110459
    [69] FENG SH J, ZUO CH, YIN W, et al. Micro deep learning profilometry for high-speed 3D surface imaging[J]. Optics and Lasers in Engineering, 2019, 121: 416-427. doi: 10.1016/j.optlaseng.2019.04.020
    [70] LIU Y, BLUNT L, GAO F, et al. High-dynamic-range 3D measurement for E-beam fusion additive manufacturing based on SVM intelligent fringe projection system[J]. Surface Topography:Metrology and Properties, 2021, 9(3): 034002. doi: 10.1088/2051-672X/ac0c62
    [71] 彭广泽, 陈文静. 基于卷积神经网络去噪正则化的条纹图修复[J]. 光学学报,2020,40(18):1810002. doi: 10.3788/AOS202040.1810002

    PENG G Z, CHEN W J. Fringe pattern inpainting based on convolutional neural network denoising regularization[J]. Acta Optica Sinica, 2020, 40(18): 1810002. (in Chinese) doi: 10.3788/AOS202040.1810002
    [72] YANG G W, YANG M, ZHOU N, et al. High dynamic range fringe pattern acquisition based on deep neural network[J]. Optics Communications, 2022, 512: 127765. doi: 10.1016/j.optcom.2021.127765
    [73] QIAO G, HUANG Y Y, SONG Y P, et al. A single-shot phase retrieval method for phase measuring deflectometry based on deep learning[J]. Optics Communications, 2020, 476: 126303. doi: 10.1016/j.optcom.2020.126303
    [74] ZHANG L, CHEN Q, ZUO CH, et al. High-speed high dynamic range 3D shape measurement based on deep learning[J]. Optics and Lasers in Engineering, 2020, 134: 106245. doi: 10.1016/j.optlaseng.2020.106245
    [75] HU Y, CHEN Q, TAO T Y, et al. Absolute three-dimensional micro surface profile measurement based on a Greenough-type stereomicroscope[J]. Measurement Science and Technology, 2017, 28(4): 045004. doi: 10.1088/1361-6501/aa5a2d
    [76] 陈龙, 王文聪, 张峰峰, 等. 基于双目结构光的术中肝脏表面局部亮度饱和分区投影[J]. 光学 精密工程,2021,29(11):2590-2602. doi: 10.37188/OPE.20212911.2590

    CHEN L, WANG W C, ZHANG F F, et al. Zonal projection based on binocular structured light for localized luminance saturation of intraoperative liver surface[J]. Optics and Precision Engineering, 2021, 29(11): 2590-2602. (in Chinese) doi: 10.37188/OPE.20212911.2590
  • 期刊类型引用(4)

    1. 汪运,郭建英,梁浚哲,朱峰,陈光希,任茂栋,梁晋. 面向高反光表面的结构光面形测量方法. 中国光学(中英文). 2025(01): 42-52 . 百度学术
    2. 张康洋,倪梓浩,董博,白玉磊. 基于贝叶斯神经网络的相位梯度计算方法. 中国光学(中英文). 2024(04): 842-851 . 百度学术
    3. 李娜. 高动态范围(HDR)技术在广播电视系统中的应用. 家庭影院技术. 2024(12): 36-38 . 百度学术
    4. 卫自强,丁苏东,胡岩,牟双,王克惠,桂坤,冯世杰,左超,陈钱. 基于过曝连通域投影强度自适应分布的三维成像(特邀). 激光与光电子学进展. 2024(20): 18-27 . 百度学术

    其他类型引用(1)

  • 加载中
图(14) / 表(4)
计量
  • 文章访问数:  970
  • HTML全文浏览量:  274
  • PDF下载量:  369
  • 被引次数: 5
出版历程
  • 收稿日期:  2023-04-16
  • 修回日期:  2023-05-15
  • 网络出版日期:  2023-09-21

目录

/

返回文章
返回