留言板

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

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

漫反射成像法的激光参数测量系统设计

程乙轮 谭逢富 何枫 侯再红 秦来安 王浩 黄志刚 吴德成

程乙轮, 谭逢富, 何枫, 侯再红, 秦来安, 王浩, 黄志刚, 吴德成. 漫反射成像法的激光参数测量系统设计[J]. 红外与激光工程, 2022, 51(9): 20210921. doi: 10.3788/IRLA20210921
引用本文: 程乙轮, 谭逢富, 何枫, 侯再红, 秦来安, 王浩, 黄志刚, 吴德成. 漫反射成像法的激光参数测量系统设计[J]. 红外与激光工程, 2022, 51(9): 20210921. doi: 10.3788/IRLA20210921
Cheng Yilun, Tan Fengfu, He Feng, Hou Zaihong, Qin Laian, Wang Hao, Huang Zhigang, Wu Decheng. Design of laser parameter measurement system based on diffuse reflection imaging[J]. Infrared and Laser Engineering, 2022, 51(9): 20210921. doi: 10.3788/IRLA20210921
Citation: Cheng Yilun, Tan Fengfu, He Feng, Hou Zaihong, Qin Laian, Wang Hao, Huang Zhigang, Wu Decheng. Design of laser parameter measurement system based on diffuse reflection imaging[J]. Infrared and Laser Engineering, 2022, 51(9): 20210921. doi: 10.3788/IRLA20210921

漫反射成像法的激光参数测量系统设计

doi: 10.3788/IRLA20210921
基金项目: 国家自然科学基金(41875033)
详细信息
    作者简介:

    程乙轮,男,博士生,主要从事激光参数测量方面的研究

    侯再红,男,研究员,博士,主要从事激光大气传输、激光大气参数测量等方面的研究

  • 中图分类号: TN247

Design of laser parameter measurement system based on diffuse reflection imaging

Funds: National Natural Science Foundation of China (41875033)
  • 摘要: 精确地测量激光在大气传输后的光斑参数,是研究激光大气传播效应和分析激光发射系统性能的关键技术手段。测量激光远场参数的方法主要包括阵列探测法和相机成像法,目前在激光大气传输效应的测量评估中大都采用阵列探测法。由于探测器阵列靶受物理空间和研发成本等因素的限制不能均匀且高分辨率紧密排布,将造成采样光斑的失真,难以精确地测量远场光斑参数。针对此问题,利用相机分辨率高的特点,设计了一套基于漫反射屏成像法的激光参数测量系统。该系统最小测量分辨力小于0.39 mm,质心位置平均偏差为0.05 mm,测量光斑到靶功率不确定度优于10%。该系统能有效地测量激光发射系统的跟瞄精度和到靶功率,为分析激光大气传输效应和分析激光发射系统性能提供有效手段。
  • 图  1  测量系统结构图

    Figure  1.  Measurement system structure diagram

    图  2  漫反射靶面示意图

    Figure  2.  Schematic diagram of diffuse reflection target surface

    图  3  畸变校正前后光斑图像。 (a) 畸变校正前光斑图像;(b) 畸变校正后光斑图像

    Figure  3.  Spot image before and after distortion correction. (a) Spot image before distortion correction; (b) Spot image after distortion correction

    图  4  系统跟瞄精度测试示意图

    Figure  4.  Schematic diagram of system tracking accuracy test

    图  5  不同功率下光斑质心位置曲线

    Figure  5.  Curve of spot centroid position under different powers

    图  6  ADU值与衰减系统透过率关系

    Figure  6.  Relationship between ADU and attenuation system transmittance

    表  1  漫反射靶面半球反射比数据

    Table  1.   Diffuse target surface hemispherical reflectance data

    Wavelength/nmDirectional hemispherical reflectance
    10300.918
    10400.918
    10500.919
    10600.919
    10640.920
    10700.920
    下载: 导出CSV

    表  2  畸变校正后特征点位置坐标

    Table  2.   Position coordinates of feature points after distortion correction

    Actual coordinate in Y direction/pixelActual coordinate in X direction/pixelIdeal
    coordinate in Y direction/pixel
    Ideal
    coordinate in X direction/pixel

    Y direction deviation/pixel

    X direction deviation/pixel
    132.92167.39132.54167.390.380
    133.28852.32135.93854.68−3.65−2.36
    818.89167.82820.93169.16−2.04−1.34
    817.74852.38820.93854.68−3.19−2.3
    下载: 导出CSV

    表  3  光斑质心位置坐标数据

    Table  3.   Coordinate data of spot centroid position

    Laser output power/WFirst location
    coordinates
    (X-direction)/mm
    Second location coordinates
    (X-direction)/mm
    Deviation/mmThird location coordinates
    (X-direction)/mm
    Deviation/mm
    ≈3000−19.03−29.2410.21−39.2810.04
    ≈1500−18.93−29.0410.10−39.0510.01
    ≈1000−18.92−29.0110.09−39.0310.02
    ≈560−18.84−28.9010.06−38.899.99
    下载: 导出CSV

    表  4  激光功率测量结果

    Table  4.   Results of laser power measurement

    Laser output
    power/W
    Attenuation system
    transmittance
    Measured laser
    power/W
    Uncertainty
    of measurement
    18701.66E-0518700
    13402.72E-0514004.29%
    12204.13411E-0511208.93%
    6209.26837E-055846.20%
    下载: 导出CSV
  • [1] Chen Y Q, Wang J H. Laser Principle [M]. Hangzhou: Zhejiang University Electronic Press, 2002: 1-39. (in Chinese)
    [2] Zuo T C. Manufacturing Laser Beam Quality Transmission Quality and Focus Quality [M]. Beijing: Science Press, 2008: 1-6. (in Chinese)
    [3] Li Bo, Wang Tingfeng, Wang Dinan, et al. Simulation of laser beam propagation through turbulence [J]. Chinese Optics, 2012, 5(3): 289-295. (in Chinese)
    [4] Zhu Wenyue, Qian Xianmei, Rao Ruizhong, et al. Evaluation technology of high energy laser atmospheric propagation performance [J]. Infrared and Laser Engineering, 2019, 48(12): 1203002. (in Chinese)
    [5] Li Gaoping, Chen Chao, Li Dong, et al. Study on parameters measurement technology of high energy and high powerlaser [J]. Journal of Applied Optics, 2020, 41(4): 645-650. doi:  10.5768/JAO202041.0409001
    [6] Feng Guobin, Yang Pengling, Wang Qunshu, et al. Measyreing technology for far-field beam profile of high power laser [J]. High Power Laser and Particle Beams, 2013, 25(7): 1615-1619. (in Chinese) doi:  10.3788/HPLPB20132507.1615
    [7] Wei Jifeng, Hu Xiaoyang, Zhang Kai, et al. Technologies and development trends ofdirectly measuring high energy laser energy [J]. Infrared and Laser Engineering, 2017, 46(7): 0706004. (in Chinese)
    [8] Yang Penglin, Feng Guobin, Wang Zhenbao, et al. Detector array for measuring far-field power density distribution of mid-infrared laser [J]. Chinese Journal of Lasers, 2010, 37(2): 521-525. (in Chinese) doi:  10.3788/CJL20103702.0521
    [9] Tan Fengfu, Huang Zhigang, Zhang Silong, et al. Detector array target calibration system based on point-by-point scanning [J]. Infrared and Laser Engineering, 2020, 49(2): 0213003. (in Chinese)
    [10] Wang Zhenbao, Feng Guobin, Yang Pengling, et al. Detector array for measuring spatio-temporal distribution of near infrared laser power density [J]. Infrared and Laser Engineering, 2011, 40(5): 935-938. (in Chinese)
    [11] Guan Youguang, Fu Shuzhen, Gao Xueyan, et al. A novel high energy laser beam detector for measuring energy distribution [J]. High Power Laser and Particle Beams, 2000, 12(2): 175-177. (in Chinese)
    [12] Hu Linting, Lu Xiankui, Jin Junkun, et al. Laser facula measurement with CCD [J]. Laser Technology, 2001(2): 154-157. (in Chinese)
    [13] ESumachev K, Sharkov V V, Savikin A P, et al. Measurement of 1.94 μm YAP: Tm3+ laser beam quality using a CCD camera [J]. Journal of Optical Technology, 2020, 87(5): 262-265. doi:  10.1364/JOT.87.000262
    [14] Meng Lingqiang, Ji Kunhao, Kong Qingqing, et al. Suppressing the influence of charge-coupled device veritical blooming on the measurement of laser beam quality factor (M2) of a near-infrared laser [J]. Applied Optics, 2018, 57(2): 130-137. doi:  10.1364/AO.57.000130
    [15] Pang Miao, Yuan Xuewen, Gao Xueyan, et al. Spot distortion calibration in measurement of laser intensity distribution based on imaging by scattering [J]. Acta Optica Sinica, 2010, 30(2): 451-455. (in Chinese) doi:  10.3788/AOS20103002.0451
    [16] Fan Yong, Zhang Jiacheng, Chen Niannian, et al. Image geometric distortion correction method [J]. Computer Engineering and Applications, 2009, 45(29): 184-197. (in Chinese)
  • [1] 纪小辉.  一种激光搜跟器的光学设计 . 红外与激光工程, 2020, 49(11): 20200088-1-20200088-7. doi: 10.3788/IRLA20200088
    [2] 魏合理, 戴聪明, 唐朝礼, 武鹏飞, 黄宏华, 李学彬, 朱文越, 饶瑞中, 王英俭.  大气参数模式及其在辐射大气传输计算中的应用 . 红外与激光工程, 2019, 48(12): 1203001-1203001(8). doi: 10.3788/IRLA201948.1203001
    [3] 朱文越, 钱仙妹, 饶瑞中, 王辉华.  高能激光大气传输性能评估技术 . 红外与激光工程, 2019, 48(12): 1203002-1203002(12). doi: 10.3788/IRLA201948.1203002
    [4] 林正国, 金星, 常浩.  脉冲激光大光斑辐照空间碎片冲量耦合特性研究 . 红外与激光工程, 2018, 47(12): 1243001-1243001(6). doi: 10.3788/IRLA201847.1243001
    [5] 马鹏阁, 陈恩庆, 庞栋栋, 羊毅.  基于快速贝叶斯匹配追踪的激光大气信道估计 . 红外与激光工程, 2017, 46(9): 922002-0922002(7). doi: 10.3788/IRLA201746.0922002
    [6] 国爱燕, 高文军, 周傲松, 程竟爽, 何善宝.  星间光通信无信标捕跟瞄技术 . 红外与激光工程, 2017, 46(10): 1022002-1022002(7). doi: 10.3788/IRLA201786.1022002
    [7] 吴佳彬, 陈云善, 高世杰, 吴志勇.  高精度近红外光斑位置检测模型研究 . 红外与激光工程, 2016, 45(7): 717001-0717001(7). doi: 10.3788/IRLA201645.0717001
    [8] 李艳晓, 张磊, 张红刚, 胡磊力.  基于双CCD 探测的外场高精度激光光斑测试技术 . 红外与激光工程, 2015, 44(1): 59-64.
    [9] 艾勇, 段梦云, 徐洁洁, 单欣, 陈晶, 熊准, 姜茹.  LC-SLM激光大气传输湍流模拟及通信实验分析 . 红外与激光工程, 2015, 44(10): 3103-3109.
    [10] 赵婧鑫, 周富强.  小尺寸光斑中心的高精度定位算法 . 红外与激光工程, 2014, 43(8): 2690-2693.
    [11] 王锋, 罗建军, 李明, 邓玥.  激光光斑整形算法 . 红外与激光工程, 2014, 43(12): 3937-3940.
    [12] 马超, 廖鹤, 周文龙, 李文峰.  卫星激光大气遥感数据反演探讨 . 红外与激光工程, 2014, 43(S1): 13-17.
    [13] 魏合理, 戴聪明.  辐射特性测量大气传输修正研究:大气辐射传输模式和关键大气参数分析 . 红外与激光工程, 2014, 43(3): 884-890.
    [14] 丛日进, 汪井源, 徐智勇, 王荣, 王喆.  长波长红外光大气信道传输 . 红外与激光工程, 2014, 43(3): 927-932.
    [15] 王飞, 徐作冬, 戢运峰, 姜畅.  采用扫描式漫反射成像法的激光强度分布测量装置 . 红外与激光工程, 2014, 43(7): 2219-2222.
    [16] 张勇, 曹喜滨, 吴龙, 赵远.  线性增益调制激光雷达系统参数对测距精度的影响 . 红外与激光工程, 2013, 42(4): 925-929.
    [17] 戴聪明, 魏合理, 陈秀红.  通用大气辐射传输软件(CART)分子吸收和热辐射计算精度验证 . 红外与激光工程, 2013, 42(1): 174-180.
    [18] 戴聪明, 魏合理, 陈秀红.  通用大气辐射传输软件(CART)大气散射辐射计算精度验证 . 红外与激光工程, 2013, 42(6): 1575-1581.
    [19] 姚佰栋, 时晶晶, 侯再红, 谭逢富, 吴毅.  基于超磁致伸缩材料的高速倾斜镜的设计与应用 . 红外与激光工程, 2013, 42(1): 96-101.
    [20] 王彦斌, 刘永志, 黄成功, 李华, 姚梅, 亓凤杰, 王敏, 邹前进.  地面对卫星瞄准精度受大气湍流的影响分析 . 红外与激光工程, 2013, 42(10): 2687-2691.
  • 加载中
图(6) / 表(4)
计量
  • 文章访问数:  133
  • HTML全文浏览量:  35
  • PDF下载量:  53
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-02-20
  • 修回日期:  2022-03-25
  • 刊出日期:  2022-09-28

漫反射成像法的激光参数测量系统设计

doi: 10.3788/IRLA20210921
    作者简介:

    程乙轮,男,博士生,主要从事激光参数测量方面的研究

    侯再红,男,研究员,博士,主要从事激光大气传输、激光大气参数测量等方面的研究

基金项目:  国家自然科学基金(41875033)
  • 中图分类号: TN247

摘要: 精确地测量激光在大气传输后的光斑参数,是研究激光大气传播效应和分析激光发射系统性能的关键技术手段。测量激光远场参数的方法主要包括阵列探测法和相机成像法,目前在激光大气传输效应的测量评估中大都采用阵列探测法。由于探测器阵列靶受物理空间和研发成本等因素的限制不能均匀且高分辨率紧密排布,将造成采样光斑的失真,难以精确地测量远场光斑参数。针对此问题,利用相机分辨率高的特点,设计了一套基于漫反射屏成像法的激光参数测量系统。该系统最小测量分辨力小于0.39 mm,质心位置平均偏差为0.05 mm,测量光斑到靶功率不确定度优于10%。该系统能有效地测量激光发射系统的跟瞄精度和到靶功率,为分析激光大气传输效应和分析激光发射系统性能提供有效手段。

English Abstract

    • 由于激光具有良好的方向性好,突出的单色性、高亮度等特点,激光技术被广泛应用到工业制造领域、国防军事领域以及生物医学领域[1-2]。准确的测量激光光斑远场分布情况是评价激光发射系统性能和激光大气传输特性最直接有效的手段[3-4]

      目前对于激光远场光斑参数测量的方法有许多,主要包括:烧蚀法、感光法、扫描法、阵列探测法、基于相机的成像法[5-7]。烧蚀法作为一种粗略的光斑尺寸测量方法,不能获得光斑的质心抖动、光斑能量等光斑参数,也不能获得激光在照射过程中光斑的实时变化;感光法同样也不能定量的获得各项测试参数;扫描法虽然能定量地获得远场光斑参数,但是其自身结构复杂需要精密的机械结构设计和复杂的逻辑系统设计。

      阵列探测法是将对特定波长敏感的探测器按照一定的规律进行排布形成阵列靶板,入射激光通过衰减取样系统后照射在探测器上,探测器通过光电效应产生电信号,再对电信号进行数字转换处理从而获得激光强度分布的测量数据[8-11]。但是由于探测器受空间物理限制和研发成本的限制不能均匀且紧密地排布,这将造成采样光斑的失真导致远场激光参数测量的精度不够等问题。

      基于相机成像法的激光参数测量系统根据光斑的取样方式又可以分为直接成像法和间接成像法。直接测量法是使入射激光直接进入相机镜头,相机实现光电信号的转换后传输至数据处理终端得到数字图像,最后通过一些数据处理方法得到光斑参数[12-13]。间接测量法是将入射激光束照射至标准的漫反射靶板上,相机在固定位置拍摄反射光后得到光斑图像数据,最终经过数据处理得到远场光斑的各项参数。

      因此,阵列探测法和基于相机成像的方法仍然是激光远场参数测量领域的研究的一个热点。由于相机具有阵列探测法所不能达到的高分辨力,在一些测量场景的应用中具有更高的精度优势。但是采用漫反射成像法时由于相机与靶面夹角的影响需要对图像进行校正[14-16],且不能保证试验中夹角大小相同,导致校正精度的降低。文中针对漫反射成像法难以系统集成和系统可能需要对采集图像需要多次校正的问题,设计了一套基于漫反射成像的激光参数测量系统。

    • 漫反射成像法由漫反射屏、相机、衰减片组合以及数据处理终端等组成。在测量过程中激光照射至标准的漫反射屏上,反射光激发相机图像传感器的光电二极管并释放电荷产生电信号,经过内部的A/D模块将电信号转化为数字信号,从而实现激光强度到灰度图像的映射,在这一过程中图像的灰度值和入射激光的光强成正相关。

      当入射光束垂直于漫反射屏时,相机与漫反射屏的法线方向会存在一定的夹角。若相机与漫反射屏法线方向夹角较小时,由相机与漫反射板间角度带来的误差对测量系统的误差影响可忽略不计。但是在激光远场光斑测量过程中,由于光束受到湍流起伏的影响光束发生扩展、漂移,到靶光斑尺寸变大,这将可能导致入射光束照射到相机,对相机造成损害。因此相机与漫反射屏的夹角需预留足够的余量(一般认为大于15°),此时测到的光斑图像会有明显的几何畸变,需要对其进行校正处理。此外,当入射激光束的能量较高时,需要在相机镜头前添加衰减组合,否则采集得到的光斑图像将饱和,无法显示图像的细节并对计算参数带来很大的误差,更严重的将会导致相机光敏面损坏。

    • 测量系统结构如图1所示,测量系统硬件主要由漫反射屏、图像采集测量系统、信标系统、安装调整支架等组成。为了使系统在降低质量的同时使各个分系统间能保持刚性连接,整体框架选用碳纤维材料和铝合金材料组合而成。漫反射屏由一块260 mm×260 mm的正方形铝板制作而成,并在表面采用喷砂镀金的方式来使靶面具有良好的均匀性和一定的抗强光能力。三根碳纤维管构成一个三角形结构,铝合金结构与碳纤维之间采用压接的形式,使铝合金固定件与碳纤维管间无相对运动,使得连接支架具有足够的刚性,进而使相机与靶面保持固定角度,并且不会产生相对运动,使测量系统具有较高的测量精度。

      图  1  测量系统结构图

      Figure 1.  Measurement system structure diagram

      该系统的相机选用一款分辨率为2 048×1088 pixel的近红外扩展型硅探测器相机,像元尺寸为5.5 μm,搭配焦距为12.5 mm的近红外镜头。根据相机的成像原理,相机与靶面的距离和采集图像分辨力的关系可表示为:

      $$ M = \frac{{5.5 \times {{10}^{ - 3}}}}{{12.5}} \times L $$ (1)

      式中:L为相机与靶面中心位置的距离(单位为mm);M为相机采集图像的分辨力。为了确保图像分辨力优于0.5 mm/pixel,则L最大为1136.4 mm。另外,在相机与靶面一定距离的情况下,为使整个漫反射靶面在相机视场中采集图像的分辨率应尽可能的大,但是随着图像分辨率的增加也会造成采样帧频的降低。

      在激光参数测量系统中,图像采集系统中单帧图像的采样时间应与一个湍流冻结时间接近(一般认为一个湍流冻结时间为10~20 ms)。所以系统中相机距靶面距离设计为880 mm,图像窗口分辨率为645×645 pixel,此时系统采样帧频大于52 Hz。

      (1)漫反射靶面及信标系统

      该系统采用漫反射靶板作为远场光斑的接收设备,如图2所示。靶面尺寸设计为260 mm×260 mm的正方形,靶面采用铝板喷砂镀金的加工方式。由于该系统在测量时与激光发射系统是一种合作信标的工作模式,相机与靶面的位置需要进行校正,所以在靶面的四个边角位置安装了四个信标光源,单灯功率大于1.5 W,发散角大于30°,信标灯的亮度可根据天空背景亮度进行调节。由于相机量子效率在650 nm波段较高且能与被测1064 nm波段激光在波长上间隔较远,通过长波通滤光片即可对其进行过滤,不会对测量系统的精度造成不利影响。

      图  2  漫反射靶面示意图

      Figure 2.  Schematic diagram of diffuse reflection target surface

      靶面的半球方向反射比(Directional HemisphericalReflectance, DHR)在1064 nm的中心波段接近92%,具有较高的反射率,说明靶面可以承受较高功率密度的激光辐照,具有较好的抗激光损失能力。在20°的范围内,双向反射分布函数(Bidirectional Reflectance Distribution Function, BRDF)因子小于0.5%,从表1可以看出,面板在喷砂镀金处理后能满足系统设计需求。

      表 1  漫反射靶面半球反射比数据

      Table 1.  Diffuse target surface hemispherical reflectance data

      Wavelength/nmDirectional hemispherical reflectance
      10300.918
      10400.918
      10500.919
      10600.919
      10640.920
      10700.920

      (2)衰减组合

      由于系统使用的是可见光相机,为了避免可见光波段的杂散光对系统测量精度造成影响,在相机镜头前加入长波通滤光片进行杂散光的滤除。在入射光束的中心波段为1050 nm,半带宽为20 nm时,透过率为85%左右,动态误差为1.46%。具有较高的均匀性,能满足系统设计需求。

    • 测量系统软件部分由数据采集模块和数据处理分析模块组成。数据采集模块实现采集参数的设置、图像的畸变校正以及储存等功能;数据处理分析部分实现图像数据的处理分析功能,得到激光束远场光斑图像的跟瞄精度、光束质量β因子、到靶功率、光斑直径、环围功率比等一系列光斑参数。

      (1) 数据采集模块

      在数据采集模块设计中,最重要的是图像的畸变校正处理,由于相机与靶面存在一定的固定角度偏差,相机在对靶面成像时会形成光斑图像的几何畸变和能量分布的畸变,所以在图像采集前要进行图像的几何畸变校正。图像几何校正的方法首先是根据图像畸变前后的图像特征点位置坐标信息,建立畸变教正模型,其次是利用已知条件确定校正模型参数,最后根据畸变前后的映射关系重新对图像进行赋值。假设畸变图像为f(x,y),畸变校正处理后的图像为F(X,Y),那么畸变前后图像坐标点(x,y)和(X,Y)的关系可表示为:

      $$ \left\{\begin{array}{l} X=j(x, y) \\ Y=i(x, y) \end{array}\right. $$ (2)

      式中:j(x,y)和i(x,y)分别表示图像上x方向坐标和y方向坐标产生几何畸变时的畸变函数。这种畸变可分为简单的线性畸变和非线性畸变,两种情况下像素的变换可由下式给出:

      $$ \left\{ {\begin{array}{*{20}{l}} {j(x,y) = {{{k}}_1}X + {k_2}X + {k_3}X} \\ {i(x,y) = {k_4}Y + {k_5}Y + k{}_6Y} \end{array}} \right. $$ (3)
      $$ \left\{ {\begin{array}{*{20}{l}} {j(x,y) = {k_1}{\text{ + }}{{{k}}_2}X + {k_3}Y + k{}_4{X^2} + {k_5}XY + {k_6}{Y^2}} \\ {i(x,y) = {k_7}{\text{ + }}{{{k}}_8}X + {k_9}Y + {k_{10}}{X^2} + {k_{11}}XY + k{}_{12}{Y^2}} \end{array}} \right. $$ (4)

      在图像畸变校正约束点的选取中,通常选用一个四边形的顶点作为约束对应点,在这种情况下畸变模型可以用一对双线性等式来表示,其表达式为:

      $$ \left\{\begin{array}{l} j(x, y)=k_{1} X+k_{2} Y+k_{9} X Y+k_{4} \\ i(x, y)=k_{5} X+k_{6} Y+k_{7} X Y+k_{8} \end{array}\right. $$ (5)

      在该系统的校正中可选取靶面四个信标灯为特征点,点亮信标灯后采集靶面图像并对其进行二值化处理利用opencv中的霍夫变换检测圆形并获取其圆心位置坐标。由于信标灯的位置在靶面上构成标准的正四方形,所以选取四点中的任意一点为起始点,任意相邻两点间的距离为边长做正四边形即可得到图像上特征点的理想位置坐标将四个点的位置坐标代入公式(5)中即可求取畸变校正函数,最后根据映射关系对其重新赋值即可完成图像的校正,图像畸变校正前后的图像如图3所示。

      图  3  畸变校正前后光斑图像。 (a) 畸变校正前光斑图像;(b) 畸变校正后光斑图像

      Figure 3.  Spot image before and after distortion correction. (a) Spot image before distortion correction; (b) Spot image after distortion correction

      图3中蓝色方框为畸变校正前获取特征点的所围四边形,红色方框为特征点的理想位置坐标。校正后特征点的理想位置坐标和实际位置坐标如表2所示,最大偏差仅为3.65 pixel,从而可以看出这种畸变校正方法是积极有效的。

      表 2  畸变校正后特征点位置坐标

      Table 2.  Position coordinates of feature points after distortion correction

      Actual coordinate in Y direction/pixelActual coordinate in X direction/pixelIdeal
      coordinate in Y direction/pixel
      Ideal
      coordinate in X direction/pixel

      Y direction deviation/pixel

      X direction deviation/pixel
      132.92167.39132.54167.390.380
      133.28852.32135.93854.68−3.65−2.36
      818.89167.82820.93169.16−2.04−1.34
      817.74852.38820.93854.68−3.19−2.3

      (2)数据处理模块

      测量系统在获取远场激光强度时空分布后即可计算激光到靶功率、光斑质心位置、质心位置抖动、光斑环围直径、激光衍射极限倍数等参数。不同的参数表征着不同的物理意义,例如,激光到靶功率表示的是激光大气传输透过率特性;光斑的质心位置和质心抖动则表示的激光发射系统的跟瞄性能;光斑环围直径、衍射极限倍数等参数则是综合的效应,受激光大气传输效应和激光发射系统性能的影响。该系统主要用于测试激光发射系统的动态跟瞄性能,所以分析光斑的质心位置和质心抖动是最直接有效的方法。

      光斑的质心位置可以通过单帧图像获得其瞬时的质心位置坐标(单帧光斑质心位置坐标)和整个出光时间内长曝光图像的质心位置坐标(积分光斑质心位置坐标),通常情况下计算光斑的质心位置坐标选用的都是重心法来计算,将光斑的重心等效为光斑的质心。单帧光斑的质心位置坐标和积分光斑的质心位置坐标可由公式(6)和公式(7)求得:

      $$ \left\{ {\begin{array}{*{20}{c}} {x(f) = \dfrac{1}{{p(f)}}\displaystyle\sum\nolimits_{i = 1}^m {\displaystyle\sum\nolimits_{j = 1}^n {{x_{ij}}{A_{ij}}(f)} } } \\ {y(f) = \dfrac{1}{{p(f)}}\displaystyle\sum\nolimits_{i = 1}^m {\displaystyle\sum\nolimits_{j = 1}^n {{y_{ij}}{A_{ij}}(f)} } } \end{array}} \right. $$ (6)
      $$ \left\{ {\begin{array}{*{20}{c}} {{{\bar x = }}\dfrac{1}{p}\displaystyle\sum\nolimits_{i = 1}^m {\displaystyle\sum\nolimits_{j = 1}^n {{x_{ij}}{A_{ij}}} } } \\ {{{\bar y = }}\dfrac{1}{p}\displaystyle\sum\nolimits_{i = 1}^m {\displaystyle\sum\nolimits_{j = 1}^n {{y_{ij}}{A_{ij}}} } } \end{array}} \right. $$ (7)

      激光发射系统的跟踪精度定义为激光远场光斑瞬时质心位置围绕其平均质心位置变化的均方根值,即激光到目标上实际的跟踪误差,可用公式(8)表示:

      $$ \left\{ {\begin{array}{*{20}{c}} {{X_{{\text{rms}}}}{\text{ = }}\sqrt {\dfrac{{\displaystyle\sum\limits_{{{f = 1}}}^F {{{\left( {x(f) - \bar x} \right)}^2}} }}{F}} } \\ {{Y_{{\text{rms}}}}{\text{ = }}\sqrt {\dfrac{{\displaystyle\sum\limits_{{{f = 1}}}^F {{{\left( {y(f) - \bar y} \right)}^2}} }}{F}} } \end{array}} \right. $$ (8)

      激光发射系统的瞄准精度可定义为激光远场光斑平均质心位置坐标与预先设定的瞄准点之间的偏差,即目标上的实际光束瞄准偏差,可表示为:

      $$ \left\{ {\begin{array}{*{20}{c}} {\Delta {\theta _{{t{p}},x}} = \arctan \left(\dfrac{{{x_c} - {x_0}}}{L}\right)} \\ {\Delta {\theta _{{{p}},y}} = \arctan \left(\dfrac{{{y_c} - {y_0}}}{L}\right)} \end{array}} \right. $$ (9)

      式中:${{{x}}_{{c}}}$为积分光斑质心位置在靶面上的X坐标;${{{y}}_{{c}}}$为积分光斑质心位置在靶面上的Y坐标;${{{x}}_{\text{0}}}$为系统设定瞄准点在靶面上的X坐标;${{{y}}_{\text{0}}}$为系统设定瞄准点在靶面上的Y坐标。

    • 对于该系统相机与靶面的距离为880 mm,相机镜头的焦距为12 mm,相机像元尺寸为5.5 μm,因此系统的分辨力为:

      $$ \frac{{5.5\; {\text{μm}}}}{{12.5\;{\text{mm}}}} \times 880\;{\text{mm}} \approx 0.39\;{\text{mm}} $$ (10)

      为了验证系统的测量精度,设计了系统跟瞄精度的测量精度测试实验。实验示意图如图4所示,激光器输出激光束经扩束系统扩束后照射在靶面上,数据处理系统对数据进行分析处理。

      图  4  系统跟瞄精度测试示意图

      Figure 4.  Schematic diagram of system tracking accuracy test

      在测试系统准备完毕的情况下通过上位机控制软件使平移台以1 cm的步长在水平方向移动,垂直方向保持静止,在不同功率下多次实验得到光斑的质心位置曲线如图5所示。实验中平移台水平移动三次,因此图5中质心偏移曲线出现三个明显的“台阶”,四次实验得到光斑质心位置在移动过程中的数据如表3所示。

      图  5  不同功率下光斑质心位置曲线

      Figure 5.  Curve of spot centroid position under different powers

      表 3  光斑质心位置坐标数据

      Table 3.  Coordinate data of spot centroid position

      Laser output power/WFirst location
      coordinates
      (X-direction)/mm
      Second location coordinates
      (X-direction)/mm
      Deviation/mmThird location coordinates
      (X-direction)/mm
      Deviation/mm
      ≈3000−19.03−29.2410.21−39.2810.04
      ≈1500−18.93−29.0410.10−39.0510.01
      ≈1000−18.92−29.0110.09−39.0310.02
      ≈560−18.84−28.9010.06−38.899.99

      综上所述,平移台在以10 mm的步长移动时系统测量得到光斑位置移动的平均偏差为0.05 mm,所以系统具有较高的测量精度。此外该系统在一定程度上兼具了测试到靶功率的能力,首先需要利用功率计对其功率进行标定,在功率标定后得到ADU (Analog Digital Unit) 值与能量的转换系数从而实现能量的测试。

      在相同的曝光时间下(9000 μs)改变系统的衰减倍率采集激光器以40%的输出功率出光得到的测量数据如图6所示,由于衰减系统采用的是中性衰减片,系统测量得到ADU值与衰减系统的透过率线性相关性达到0.99。

      图  6  ADU值与衰减系统透过率关系

      Figure 6.  Relationship between ADU and attenuation system transmittance

      图6可以看出,系统测量得到ADU值与衰减系统的透过率成正相关,而ADU值与曝光时间的大小也成正相关因此可以得到如下的能量计算公式:

      $$ {{P}} ={\rm{ ADU}} \times \left( {\frac{{{\text{Expose Tim}}{{\text{e}}_{\text{0}}}}}{{{\text{Expose Tim}}{{\text{e}}_{{N}}}}}} \right)\left( {\frac{{{\text{Transmittanc}}{{\text{e}}_{\text{0}}}}}{{{\text{Transmittanc}}{{\text{e}}_{{N}}}}}} \right) \times K $$ (11)

      式中:ADU为实际测试时得到ADU值;Transmittance0和Expose Time0为系统进行能量标定时衰减系统的透过率和系统曝光时间;k为标定功率下激光器实际出光功率与ADU值的比值。实际使用过程中通过改变曝光时间去调整系统的动态范围时会造成系统采样频率的变换,所以在固定的曝光时间下改变衰减系统的透过率更为合适,那么激光到靶功率可由公式(12)求得:

      $$ {{P}} = {\rm{ADU}} \times \left( {\frac{{{\text{Transmittanc}}{{\text{e}}_{\text{0}}}}}{{{\text{Transmittanc}}{{\text{e}}_{{N}}}}}} \right) \times K $$ (12)

      以激光器70%的出光功率、曝光时间为20000 μs、衰减透过率为1.66E-5的情况为标定条件得到的出光功率与ADU值的系数k=6.30E-5。改变激光器的出光功率,并通过改变衰减系统的透过率测量得到的激光功率如表4所示,表中激光输出功率为校准后的功率计测量得到,且其测量不确定度为2.1%。

      表 4  激光功率测量结果

      Table 4.  Results of laser power measurement

      Laser output
      power/W
      Attenuation system
      transmittance
      Measured laser
      power/W
      Uncertainty
      of measurement
      18701.66E-0518700
      13402.72E-0514004.29%
      12204.13411E-0511208.93%
      6209.26837E-055846.20%
    • 针对漫反射成像法的激光参数测量系统,设计了以铝板为基底的喷砂镀金强光接收面板、相机以及衰减系统的激光参数测量系统。各分系统以碳纤维管为骨架进行连接,具有较高的稳定性和可靠性,避免了利用相机在图像采集过程中需要进行校正的问题。在跟瞄精度的测试实验中,系统测量分辨率为0.39 mm,质心位置平均偏差为0.05 mm;根据衰减系统的衰减特性拟合得到激光到靶能量的计算公式,不同功率下激光到靶能量测量不确定度优于10%。该系统的成功研制为准确测量激光发射系统的跟踪精度和瞄准精度提供了一种新的方法,对评估激光发射系统的性能和研究激光大气传输特性具有重要的意义和工程应用价值。

参考文献 (16)

目录

    /

    返回文章
    返回