主三镜一体化离轴三反光学天线设计

任兰旭, 谭庆贵, 张缓缓, 呼新荣, 李帅

任兰旭, 谭庆贵, 张缓缓, 呼新荣, 李帅. 主三镜一体化离轴三反光学天线设计[J]. 应用光学, 2021, 42(1): 23-29. DOI: 10.5768/JAO202142.0101004
引用本文: 任兰旭, 谭庆贵, 张缓缓, 呼新荣, 李帅. 主三镜一体化离轴三反光学天线设计[J]. 应用光学, 2021, 42(1): 23-29. DOI: 10.5768/JAO202142.0101004
REN Lanxu, TAN Qinggui, ZHANG Huanhuan, HU Xinrong, LI Shuai. Design of off-axis three-mirror optical antenna based onintegration of primary-tertiary mirror[J]. Journal of Applied Optics, 2021, 42(1): 23-29. DOI: 10.5768/JAO202142.0101004
Citation: REN Lanxu, TAN Qinggui, ZHANG Huanhuan, HU Xinrong, LI Shuai. Design of off-axis three-mirror optical antenna based onintegration of primary-tertiary mirror[J]. Journal of Applied Optics, 2021, 42(1): 23-29. DOI: 10.5768/JAO202142.0101004

主三镜一体化离轴三反光学天线设计

基金项目: 装备预研航天科技联合基金(6141B06301)
详细信息
    作者简介:

    任兰旭(1987−),男,博士,工程师,主要从事空间激光通信光机技术研究。E-mail:ren_lanxu@126.com

  • 中图分类号: TN929.1

Design of off-axis three-mirror optical antenna based onintegration of primary-tertiary mirror

  • 摘要: 为提高空间激光通信离轴三反光学天线的装调效率,需要对反射镜总的装调自由度进行优化设计。基于共轴三反无焦系统的设计原理,提出一种主三镜一体化的设计方法,并推导了光学系统结构参数之间的关系,利用Zemax光学软件设计了一种结构紧凑、主三镜一体化的离轴三反光学天线。设计结果表明:全视场光学像质优于衍射极限,主镜和三镜空间位置相接近,可共母板加工和面形共基准检测,为主三镜的一体化加工提供了方法。光学天线的装调自由度减少了6个,降低了装调难度,提高了装调效率。
    Abstract: In order to improve the installation and adjustment efficiency of off-axis three-mirror optical antenna of space laser communication, it is necessary to optimize the design of total freedom of the mirror. Based on the principle of coaxial three-mirror afocal system, a method of integration of primary-tertiary mirror was proposed. The relationship between structural parameters of optical system was deduced, and an off-axis three-mirror optical antenna with compact structure and primary-tertiary mirror integration was designed by Zemax optical software. The design results show that the optical image quality of full field of view is better than diffraction limit, the primary mirror and tertiary mirror are close to each other in space, the common mother plate can be processed and surface shape can be detected by the common reference, which provides a method for integrated processing of primary mirror and three-mirror. The degree of installation and adjustment freedom of optical antenna is reduced by 6, which reduces the difficulty and improves the efficiency.
  • 核电站燃料组件受水下工作环境、装配应力和热应力等因素的影响,在工作中产生局部变形等问题,因此,核电站水下核燃料棒的高精度检测与维护是其安全运行的重要保障[1]。由于水下检测条件恶劣,水下核燃料棒高精度三维测量面临着巨大的挑战。基于线激光扫描的单目视觉测量方法采用单目摄像机在被测物表面做左右方向的扫描运动,获取线激光条纹信息,通过提取激光条纹中心,将其转换为三维点云坐标,实现被测物尺寸的高精度测量。水下检测中的激光条纹照射质量、水下检测环境和待测物表面特性都会对摄像机采集的光条图像造成影响[2],从而给中心线的提取带来一系列问题。因此,如何精确快速地提取光条中心线,是保证整个检测系统测量精度的关键一步。

    近年来,国内外学者提出了许多条纹中心提取方法[3-4],主要有几何法[5]、极值法[6]、骨架提取法[7]、高斯拟合法[8]、方向模板法[9]、灰度重心法[10]等。极值法运算速度快,但若光条出现多个极值,定位精度会明显降低;骨架提取法方法简单,但其精度受环境影响大。田庆国[11]提出基于轮廓线多边形提取方法,可快速定位光带横截面方向并进行中心提取;方向模板法稳定性较弱,中心提取效果一般,王泽浩[12]提出自适应方向模板法,提高了方向模板的提取稳定性;高世一[13]提出变边限高斯拟合提取激光条纹中心线法,可适应条纹宽度进行高斯拟合;吴家勇[14]针对灰度梯度特性,提出了梯度重心法,有效提取了光条中心;Steger[15]基于Hessian矩阵得到光条横截面上一阶导数零点位置,求得光条纹中心线的亚像素位置,其考虑了光条几何特性,定位精度高,但无法解决光条中心多极值的问题,且运算速度慢,无法满足实时测量需求;灰度重心法考虑了光条灰度特性,运算快,熊会元[16]提出应用灰度梯度计算光条法线方向,在法线方向上应用灰度重心法精确提取光条中心的方法,解决了被测物体表面曲率变化大的问题;程鑫[17]结合Hessian矩阵对光条的中心进行定位,利用高斯曲线准确拟合出激光条的中心。以上算法提高了光条中心提取的速度及鲁棒性,但对于提取水下高反光表面光条中心的问题未能很好解决。

    基于上述分析,本文提出了一种反光表面条纹自适应中心提取算法,基于光条几何特性自适应修正高反光区域轮廓,结合图像灰度卷积处理方法自适应修正高反光区域灰度分布,经灰度重心法确定中心点的亚像素坐标,实现了水下核燃料棒表面激光条纹中心的高精度提取。该方法可以有效地减少水下噪声及核燃料棒反光影响,准确、快速地提取激光条纹中心,控制三维重建精度在0.2 mm以内,满足核燃料棒变形测量精度0.3 mm的需求。

    当激光器投影到物体表面,将形成一定宽度的激光条纹,其光强从条纹中心展开,满足高斯分布。其条纹截面灰度值分布曲线的表达式为

    $$ G\left(x\right) = \dfrac{S}{\sqrt{2{\text{π}} \sigma }}·\exp\left[-\dfrac{{\left(d-\mu \right)}^{2}}{2{\sigma }^{2}}\right]$$ (1)

    式中:S为灰度值的幅值;σ为标准差;μ为光条截面坐标;d为光条中心坐标。

    因激光器自身质量、测量环境及被测量物体材质的不同,单目摄像机所采集的激光条纹灰度分布差异性较大,为此,截取核燃料棒原始图像中激光条纹部分,绘制光条左部和中部截面的灰度分布曲线,见图1所示。

    图  1  原始激光条纹及其截面灰度分布
    Figure  1.  Original laser stripes and section gray distribution

    根据图1分析可以得到原始激光条纹图像具有以下特征:1) 光条边界附近存在水下环境噪声干扰;2) 光条截面灰度极值不同,截面宽度存在差异;3) 光条中间部分存在高反光区域,轮廓改变,灰度分布不满足近似高斯分布。

    通过以上特征可以得到,受水下测量环境影响,获取的激光条纹图像存在大量噪声,影响激光条纹分离提取;受物体表面反光特性影响,光条中部截面处存在多灰度极值的情况,靠近光条中心灰度值迅速增加;受金属镜面反光特性影响,条纹中心出现垂直方向的反光噪点,通过调整图片曝光时间及拍摄角度可以减轻;通过传统提取中心算法难以获得精确激光条纹中心,中心线提取依赖光条中灰度分布符合高斯分布的截面信息。

    文中出现较多物理量,主要分为3类。第1类是像素灰度值,其中ω表示像素点经预处理的灰度值;ω′表示像素点经自适应阈值提取的灰度值;ω″表示像素点经过灰度修正的灰度值。第2类是灰度分布,其中Tn表示灰度修正卷积模板;Gij)表示光条所在矩形内的点集。第3类是光条几何特征,其中wn表示激光条纹宽度;kv表示光条轮廓的切线方向;kn表示光条的截面方向。

    实际图像中因其他因素产生的噪声,需要通过图像去噪将其消除。通过对图像中噪声研究发现,水下核燃料图像中的噪声主要分为两种:一是孤立噪点或者由几个噪点组成的假光点;二是沿核燃料棒上下边缘的线状反光噪声。通过滤波方法及卷积运算消除噪声,保证对激光条纹区域提取的精确性。

    1) 首先利用中值滤波处理图像,通过在噪点周围定义一个5×5的邻域窗口,对窗口中像素点灰度值进行升序排列,取灰度值的中值作为滤波结果,可以去除水下核燃料棒激光条纹图像中的椒盐噪点。

    2) 对原始图像进行全局阈值的二值化处理,根据原始图像选择合适的灰度阈值Ct,计算实现图像二值化:

    $$\omega (\mu ,\nu ) = \left\{ {\begin{array}{*{20}{c}} {1,\omega \left( {\mu ,\nu } \right) \geqslant {C_{{t}}}} \\ {0,\omega \left( {\mu ,\nu } \right) \leqslant {C_{{t}}}} \end{array}} \right.$$ (2)

    式中ωμν)表示二值化后的光条在坐标(μν)处的灰度值。

    在灰度值非零的点上利用3×3模板去除噪点,模板如图2所示,a1a8分别对应图像中的每个像素的位置。

    图  2  卷积模板
    Figure  2.  Convolution template

    将该卷积模板在图像上移动,根据(3)式计算判定该点是否为噪点:

    $$A = \sum\limits_{r = 0}^8 {{{\left( {1 - {a_r}} \right)}^2}} $$ (3)
    $$\omega (\mu ,\nu ) = \left\{ {\begin{array}{*{20}{c}} {1,A \leqslant 4} \\ {0,A > 4} \end{array}} \right.$$ (4)

    通过对点(μν)的A值进行判断,去除图像中的孤立噪点和沿核燃料棒上下边缘的反光噪点。将去噪图像与原始图像各像素灰度值相乘,获得具有原始灰度信息的去噪激光图像,最后使用高斯滤波平滑图像。

    3) 根据图像灰度值,实现激光条纹的自适应阈值提取。从垂直激光平面方向开始扫描,寻找包含像素点的灰度值为255的行。当第一次检测到时,记该点为(iminjmin),最后一次检测到时,记该点为(imaxjmax),框选出包含光条位置的矩形,矩阵宽度为wd。则光条所在的矩形内点集Gij)为

    $$ \begin{split} \Bigg\{ G(i,j)\Bigg|& {i_{\min }} \leqslant i \leqslant {i_{\max }},\dfrac{{{j_{\min }} + {j_{\max }} - {w_{\rm{d}}}}}{2} <j < \\ & \dfrac{{{j_{\min }} + {j_{\max }} + {w_{\rm{d}}}}}{2} \Bigg\} \end{split} $$ (5)

    根据每列像素灰度值对该光条截面的不同影响程度确定阈值,选择使用基于梯度的自适应行阈值法在矩形内计算每列的阈值。如图3所示,将激光条纹图像中任一列的灰度分布抽象为两条一次函数:

    图  3  自适应阈值原理
    Figure  3.  Principle of self-adaptive threshold

    计算两像素点间灰度梯度Vmn

    $${V_{mn}} = \left| {{\omega _{mn}} - {\omega _{(m - 1)n}}} \right|$$ (6)

    式中ωmn)是点(μν)在图像第mn列的灰度值。

    计算每行灰度梯度的平均值Bn

    $${B_n} = \dfrac{{\displaystyle\sum\limits_{m = 2}^M {{V_{mn}}} }}{{D - 1}}$$ (7)

    式中D为光条大致宽度。

    在进行测量时,D为25个像素。通过下式判断点(μν)是否属于光条有效点:

    $${g_n} = 1 - \dfrac{{{B_n} \times D}}{{{G_{\max }}}}$$ (8)
    $$ {\omega' }(\mu ,\nu )=\left\{\begin{array}{c}\omega (\mu ,\nu ),\omega (\mu ,\nu )\geqslant {g}_{n}·{G}_{\max}\\ 0,\omega (\mu ,\nu )<{g}_{n}·{G}_{\max}\end{array}\right.$$ (9)

    式中:gn表示阈值界限系数;Gmax表示每行灰度的最大值;ω'μν)表示自适应阈值提取的光条在点(μν)处的灰度值。

    该方法能使光条的提取阈值根据列像素自适应变化,快速将激光从环境背景中提取,图4(a)为经自适应分割的光条图像。通过图4(b)分析,光条附近噪声被去除,提取出的光条图像边界明显,保证了提取精度。

    图  4  自适应阈值分割及其激光截面灰度分布图
    Figure  4.  Self-adaptive threshold segmentation and gray distribution of laser section

    提取的光条已从图像背景分割,但由于反光现象,激光条纹在被测物体表面的轮廓发生改变,此时直接进行光条中心线提取无法达到精度要求。为解决这一问题,选取灰度分布满足高斯分布的光条截面端点,采用3层BP神经网络对所选坐标点进行曲线拟合,拟合出两条激光条纹轮廓,对光条进行进一步的修正提取。

    图5所示,BP神经网络设定输入层结点为X,隐藏层节点为h,输出层节点为Y,输入层和输出层节点数为1,隐藏层的节点数C为2。设定wic为输入层到隐层第c个节点之间的连接权值,ψc表示隐层第c个节点的阈值;woc为隐层第c个节点到输出层之间的连接权值,θ表示输出层节点的阈值;隐层输入变量为hicc=1,2,隐层输出变量为hocc=1,2。

    图  5  神经网络的基本原理
    Figure  5.  Basic principle of neural network

    隐层激活函数fh选用(10)式Sigmoid函数,输出层激活函数fo则选用purelin函数。

    $${f_{{h}}}(x) = \dfrac{1}{{1 + \exp \left( { - x} \right)}}$$ (10)

    隐层的第c个节点的输入为

    $${h_{ic}} = {w_{ic}}X - {\psi _c}$$ (11)

    隐层的第c个节点的输出为

    $${h_{oc}} = {f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right)$$ (12)

    输出层节点的输入为

    $${y_i} = {w_{oc}}{h_{oc}} - \theta = {w_{oc}}{f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right) - \theta $$ (13)

    输出层节点的输出为

    $${y_{{o}}} = {f_{{o}}}\left( {{y_i}} \right) = {f_{{o}}}\left[ {{w_{oc}}{f_{{h}}}\left( {{w_{ic}}X - {\psi _c}} \right) - \theta } \right]$$ (14)

    设定Yk)为第k个点的纵坐标,对图像中提取的总数为K的像素点拟合轮廓,在网络反向传播时神经网络得到的轮廓位置误差为

    $${E^{\left( t \right)}} = \dfrac{1}{2}\sum\limits_{k = 1}^K {{{\left( {{y_{{o}}}^{\left( t \right)}\left( k \right) - Y\left( k \right)} \right)}^2}} $$ (15)

    式中yotk)为神经网络权值在第t次调整后得到的第k个点的纵坐标。

    根据Levenberg-Marquardt最优化算法调整,得到最终神经网络权值。训练直到网络输出的误差减少到给定误差,或进行到预先设定的训练次数时,可以停止训练。设定目标误差为0.01,设定最大训练次数为1 000次,训练数据选取各行两端的非高反光区域像素坐标点,确定其为曲线轮廓拟合的真实值,基于条纹轮廓真实值进行预测,最终得到激光条纹的拟合曲线如图6所示。

    图  6  轮廓拟合
    Figure  6.  Contour fitting

    分别提取光条两边有效像素点进行轮廓拟合,可以获得经过轮廓修正的激光条纹图像,如图7所示。

    图  7  轮廓修正后的激光条纹
    Figure  7.  Laser stripe after contour correction

    图8所示,两组数据分别为轮廓修正前后的激光条纹各行宽度,因反光干扰,修正前条纹各行宽度变化较大。通过BP神经网络对激光条纹的轮廓进行拟合,修正表面反光造成的光条宽度变化的部分,解决了激光条纹轮廓变形问题。

    图  8  轮廓修正前后的光条宽度
    Figure  8.  Width of light bar before and after contour correction

    激光条纹反光部分轮廓已经修正,但各像素灰度值仍未改变。此时若对光条进行中心线提取,反光区域将存在多极值点的情况。本文方法通过设定自适应卷积模板修正光条灰度分布,利用条纹宽度信息生成模板对图像进行卷积运算,获得经灰度分布正常的激光条纹图像。因为线激光照射在被测物体表面后其灰度分布近似符合高斯分布,所以灰度修正系数将通过高斯函数设定。同一列中,像素修正系数由中部向列两边依次递减,激光条纹图像中任意一点p灰度修正系数表示为

    $$\begin{split}& {f}_{p}\left({w}_{n},{m}_{p},{m}_{\min},\sigma \right) = \dfrac{\delta }{\sqrt{2{\text{π}} \sigma }}\cdot \\ &\quad\quad\exp\left\{-\dfrac{{\left[{m}_{p}-\left({m}_{\min}+\dfrac{{w}_{{n}}}{2}\right)\right]}^{2}}{2{\sigma }^{2}}\right\}\end{split}$$ (16)

    式中:wnn=0,1,…,N−1)为经过轮廓修正的激光条纹图像中点p所在第n列的非零灰度值像素个数,即光条带宽;N为含有非零灰度值的总列数,表示激光条纹长度;mp为点p所在行数;mmin为点p所在的n列中第1个非零像素点的行数;б为高斯函数标准差,取值为4;δ为固定系数,取值为0.95。

    因为点p的灰度修正系数已经求得,所以激光条纹上的点p修正灰度值如(17)式所示:

    $$ {\omega'' }\left({\mu }_{p},{\nu }_{p}\right)={f}_{p}\left({w}_{n},{m}_{p},{m}_{\min},\sigma \right)·{\omega '}\left({\mu }_{p},{\nu }_{p}\right)$$ (17)

    式中:ω″(μpνp)为点(μν)经过灰度修正后的灰度值。

    已知任意点的灰度修正系数后,可设计激光条纹图像中任意列的灰度修正卷积模板。卷积模板的行数设为激光条纹宽度wn,列数设为1,此时激光条纹的每个截面都可生成与其截面宽度相适应的灰度修正卷积模板Tnn= 0,1,…,N−1)。由(16)式知fp的范围在(0~1)之内,且根据光条截面宽度wnfp在模板中间时取最高值,然后依次向两边递减。激光条纹图像中任意列n的灰度修正卷积模板为

    $${T_n} = {\left( \begin{array}{c} {f_1} \\ {f_2} \\ \vdots \\ {f_{{M_n} - 1}} \\ {f_{{M_n}}} \\ \end{array} \right)_{{M_n} \times 1}}$$ (18)

    式中Mn为激光条纹第n列含有非零灰度值的总行数。

    根据激光条纹宽度生成自适应灰度修正卷积模板,在光条原始灰度分布矩阵Iβα)上自左向右进行滑动运算,可以得到激光光条修正后的灰度分布:

    $$ {{I'}}(\beta ,\alpha ) = {\left( {\begin{array}{*{20}{c}} 0&\cdots& 0&{{f_{m1}} \cdot G(m,1)}&0&\cdots&0\\ \vdots& & & & & &\vdots \\ 0& &\vdots & \vdots & \vdots &&0\\ {{f_{11}} \cdot G(1,1)}&& & \cdots & & &{{f_{\alpha 1}} \cdot G(\alpha ,1)}\\ \vdots && \cdots &{{f_{mn}} \cdot G(m,n)}& \cdots && \vdots \\ {{f_{1{\beta _1}}} \cdot G(1,{\beta _1})}&& & \cdots & & &{{f_{\alpha {\beta _\alpha }}} \cdot G(\alpha ,{\beta _\alpha })}\\ 0&& \vdots & \vdots & \vdots&&0\\ \vdots& & & & && \vdots \\ 0&\cdots&0&{{f_{m{\beta _m}}} \cdot G(m,{\beta _m})}&0&\cdots&0 \end{array}} \right)_{\alpha \times \beta }} $$ (19)

    式中:α为光条原始图像中含有非零灰度值的总列数;βm为第m列中含有非零灰度值的总行数,且1≤mαβ为激光条纹各列中含有非零灰度值总行数的最大值。

    对经轮廓修正的水下核燃料棒扫描图像进行卷积灰度修正,其结果如图9(a)所示。

    图  9  灰度修正后的激光条纹及其截面灰度分布
    Figure  9.  Laser stripe after gray correction and section gray distribution

    激光条纹反光部分灰度值过高,灰度分布改变,因此通过自适应卷积模板处理激光条纹图像灰度,对中间高反光区域进行灰度值降低使其满足高斯分布规律,如图9(b)所示,与图4(b)对比所见,获得了良好效果。

    采用灰度重心法对完成修正的光条图像进行激光条纹亚像素中心提取。因为核燃料棒是圆柱体,所以线激光平面与核燃料棒表面相交形成的激光条纹轮廓通过单目摄像机采集后成曲线轮廓。对于非直线激光条纹的中心提取,需要在其光条截面方向进行。如图10所示,因为光条轮廓线已确定,所以可以通过计算其轮廓线上任一像素点的切线kv,确定光条的截面方向kn

    图  10  激光截面中心
    Figure  10.  Center of laser section

    取轮廓线上一点为q,通过光条轮廓线在q点的一阶导数可求得该点在光条轮廓上的切线斜率kv。此时易得出经q点的光条截面方向斜率kn

    $${k_n} = - \dfrac{1}{{{k_v}}}$$ (20)

    光条截面方向线条经过q点后向前后方向延长,直至抵达光条边界像素点,自适应光条宽度。通过判断像素点所在的行向量与列向量是否与斜率kn同时相交,确定其是否被光条截面方向线条经过。将符合条件的像素点记为一组像素点集合Xq={xh| xh=(μhνh),h=1,2,…,H},其中H为光条截面方向斜率kn所经过的像素点总数。q点所在截面的亚像素中心通过在该像素点集合Xq内使用灰度重心法求解:

    $$ {q_c}\left( {{x_q},{y_q}} \right) = \left\{ {\begin{array}{*{20}{c}} {{x_q} = \dfrac{{\displaystyle\sum\limits_{h = 1}^H {{\mu _h} \cdot \omega ''({\mu _h},{\nu _h})} }}{{\displaystyle\sum\limits_{h = 1}^H {\omega ''({\mu _h},{\nu _h})} }}}\\ {{y_q} = \dfrac{{\displaystyle\sum\limits_{h = 1}^H {{\nu _h} \cdot \omega ''({\mu _h},{\nu _h})} }}{{\displaystyle\sum\limits_{h = 1}^H {\omega ''({\mu _h},{\nu _h})} }}} \end{array}} \right. $$ (21)

    式中:xqq点所在光条截面的亚像素中心横坐标;yq为纵坐标。

    采用Visio studio 2013作为软件开发平台、结合OpenCV 3.0库及PCL 1.8.0点云库实现文中方法。如图11所示,采用的水下核燃料棒线结构光测量平台,搭载Mako G-507 PoE 工业相机,其分辨率为2 464(H)× 2 056(V)像素,像素尺寸为3.45 μm×3.45 μm;激光器波长为520 nm。通过经水下标定的水下核燃料棒测量系统对核燃料棒尺寸进行测量。

    图  11  水下核燃料棒测量系统
    Figure  11.  Underwater nuclear fuel rod detection system

    图12(a)12(c)分别选用传统灰度重心法、基于最小二乘拟合中心提取算法[18]和文中算法进行激光条纹中心提取。从图12(a)可以看出,传统灰度重心法提取的光条中心线条受反光影响,光条中心点存在前后错开的现象,线条不连续;图12(b)中基于最小二乘拟合中心提取算法对光条进行提取,光条中心线条连续性增加,但由于其提取的中心点与拟合中心点重合不足,存在线条中心不唯一的情况,容易产生噪点;图12(c)为文中算法,算法充分利用了光条的有效数据,光条中心线条连续性较好,更客观地反映了光条中心分布情况。

    图  12  不同激光条纹中心提取方法对比
    Figure  12.  Comparison of different laser stripe center extraction methods

    使用3种算法提取激光条纹中心线后生成点云进行对比。图13(a)点云为传统灰度重心法生成,图13(b)点云为基于最小二乘拟合中心提取算法生成,图13(c)点云为文中算法生成。结构光中心线提取不准确,将导致通过图像坐标点求解出的三维坐标发生较大偏差,产生点云离群现象。可以看出,通过传统灰度重心法和基于最小二乘拟合中心提取算法得到的点云含有较多的离群噪点,被测物体各部分的反光导致多数截面光条中心线条出现断点和噪点,这些会影响点云测量的准确性;通过文中算法解决因反光造成的结构光中心线提取不准的现象,可以明显看到所获得的点云大量噪点被消除,圆柱点云表面更加光滑。

    图  13  点云计算结果
    Figure  13.  Calculation results of point cloud

    通过对10组测量提取的激光条纹中心线生成的点云进行直径测量,与核燃料棒人工采用游标卡尺测量的结果进行比对。表1所示为10组拟合直径统计结果。针对水下核燃料棒检测要求,测量数据应保持高精度与稳定性,传统灰度重心法在水下测量中稳定性较好,但是其测量数据误差较大,基于最小二乘拟合中心提取算法的测量数据的最大误差和平均误差有所减少,但是其稳定性不足。文中算法相较于传统灰度重心法和基于最小二乘拟合中心提取算法平均误差明显降低,分别为97%和89%,测量稳定性分别提高9%和33%,测量效果良好,拟合圆柱的误差在0.2 mm以内。可以看出,文中提出的算法能够有效减少核燃料棒水下测量的噪声,具有较好的稳定性和精确性。

    表  1  核燃料棒测量结果
    Table  1.  Measurement results of nuclear fuel rod mm
    方法灰度重
    心法
    基于最小二乘拟合
    中心提取算法
    文中算法
    测量值9.3979.4439.528
    9.3819.4829.822
    9.7319.5449.636
    9.4899.7639.732
    9.4789.7619.548
    9.6189.6329.586
    9.5829.8879.736
    9.7219.4649.622
    9.4919.8349.612
    9.5109.5529.841
    最大误差0.28900.21700.1710
    平均误差0.13020.03380.0037
    标准差0.12140.16330.1102
    下载: 导出CSV 
    | 显示表格

    针对水下核燃料棒检测过程中激光条纹受水下环境噪声及反光影响的问题,本文根据激光条纹特征提出了结合BP神经网络轮廓拟合与卷积运算的自适应激光条纹亚像素中心提取方法。文中算法通过对核燃料棒反光部分的几何特性与灰度特征进行自适应修正,实现了水下核燃料棒表面激光条纹中心线的高精度提取。实验证明,使用该方法实现了在水下测量650 mm处燃料棒尺寸误差控制在0.3 mm以内。与其他算法进行比较,平均误差分别降低97%和89%,稳定性分别提高9%和33%,验证了该算法测量的精确性和鲁棒性。文中算法对于核燃料棒尺寸的准确测量,保证核燃料组件安全运行具有重要意义。

  • 图  1   共轴三反无焦系统结构示意图

    Figure  1.   Structure diagram of coaxial three-mirror afocal system

    图  2   次镜遮拦比${\alpha _1}$与次镜焦距${f_2}^\prime $、出瞳距$L$的关系

    Figure  2.   Relationship between ${\alpha _1}$and ${f_2}^\prime $$L$

    图  3   主三镜一体化面形检测光路示意图

    Figure  3.   Schematic of surface test optical path of primary-tertiary mirror integration

    图  4   离轴三反光学天线光路图

    Figure  4.   Optical path diagram of off-axis three-mirror optical antenna

    图  5   离轴三反光学天线像质评价图

    Figure  5.   Image quality evaluation diagram of off-axis three-mirror optical antenna

    表  1   光学系统设计参数要求

    Table  1   Parameters of optical system design

    ParameterValue
    Work spectrum/μm1.55
    FOV/mrad≥2
    Entrance pupil diameter D/mm≥150
    Compression ratio−15×
    Wave aberration RMS≤λ/15
    Optical isolation/dB≤−50
    Total length/mm≤350
    下载: 导出CSV

    表  2   优化后的光学系统结构参数

    Table  2   Structural parameters of optical system after optimization

    ElementRadius of curvature/mmThickness/mmConic type
    STO400
    Primary mirror−600 −290.9−1
    Secondary mirror−20.225 290.816−1.499
    Ternary mirror−399.738−250−3.777
    Fold mirrorInfinity341
    下载: 导出CSV

    表  3   光学天线的系统公差

    Table  3   System tolerance of optical antenna

    TolerancePrimary mirrorSecondary mirrorTernary mirrorFold mirror
    Wavefront tolerances/λ0.020.020.02
    Decentered tolerances/mm0.010.01
    Thickness tolerances/mm0.010.01
    Tilt tolerances/(″)1060
    Radius tolerance/mm0.060.010.1
    Conical coefficient tolerance0.0010.0010.001
    下载: 导出CSV

    表  4   蒙特-卡洛分析结果

    Table  4   Monte-Carlo analysis results

    Monte-Carlo probability/%0 mrad0.7 mrad1 mrad
    900.0350.0380.040
    800.0270.0300.030
    500.0150.0180.019
    下载: 导出CSV
  • [1]

    ROBERTS W T. Discovery deep space optical communication (DSOC) transceiver[J]. SPIE,2017,10096:VNSP100960V.

    [2]

    BISWAS A, KOVALIK J, WRIGHT M W. LLCD operations using the Optical Communications Telescope Laboratory (OCTL)[J]. SPIE,2014,8971:89710X.

    [3]

    YAN Peipei, DENG Xiaoguo, ZHANG Henjin. Optical design of a new type space laser communication optical system[J]. SPIE,2013,8906:890602.

    [4] 杨成龙, 颜昌翔, 杨宇飞. 星间激光通信终端光学天线的隔离度[J]. 中国光学,2017,10(4):462-468.

    YANG Chenglong, YAN Changxiang, YANG Yufei. Isolation of optical antenna of inter-satellites laser communication terminals[J]. Chinese Optics,2017,10(4):462-468.

    [5] 关姝, 王超, 佟首峰, 等. 空间激光通信离轴两镜反射光学天线自由曲面光学天线设计[J]. 红外与激光工程,2017,46(12):1222003-1-8.

    GUAN Shu, WANG Chao, TONG Shoufeng, et al. Optical antenna design of off-axis two-mirror reflective telescope with freeform surface for space laser communication[J]. Infrared and Laser Engineering,2017,46(12):1222003-1-8.

    [6] 韩培仙, 任戈, 刘永, 等. 可见/中波双波段共口径光学系统设计[J]. 应用光学,2020,41(3):435-440.

    HAN Peixian, REN Ge, LIU Yong, et al. Optical design of VIS/MWIR dual-band common-aperture system[J]. Journal of Applied Optics,2020,41(3):435-440.

    [7] 金光, 李艳杰, 钟兴, 等. 空间成像与激光通信共口径光学系统设计[J]. 光学与精密工程,2014,22(8):2067-2074.

    JIN Guan, LI Yanjie, ZHONG Xin, et al. Design of co-aperture optical system for space imaging and laser communication[J]. Optics and Precision Engineering,2014,22(8):2067-2074.

    [8]

    TUELL M T, MARTINA H M, BURGEA J H, et al. Optical testing of the LSST combined primary/tertiary mirror[J]. SPIE,2010,7739:77392V-10.

    [9] 高铎瑞, 胡辉, 汪伟, 等. 一种改进的易于加工、装调的离轴三反光学系统设计[J]. 中国激光,2015,42(6):2067-2074.

    GAO Duorui, HU Hui, WANG Wei, et al. Optical design for a kind of simply processed and installed off-axis three-reflection[J]. Chinese Journal of lasers,2015,42(6):2067-2074.

    [10] 孟庆宇, 王维, 纪振华, 等. 主三一体化离轴三反光学系统设计[J]. 红外与激光工程,2015,44(2):578-582. doi: 10.3969/j.issn.1007-2276.2015.02.033

    MENG Qingyu, WANG Wei, JI Zhenhua, et al. Design of off-axis three-mirror system based on integration of primary and tertiary mirror[J]. Infrared and Laser Engineering,2015,44(2):578-582. doi: 10.3969/j.issn.1007-2276.2015.02.033

    [11] 沙巍, 陈长征, 许艳军, 等. 离轴三反空间相机主三镜共基准一体化结构[J]. 光学与精密工程,2015,26(6):1612-1619.

    SHA Wei, CHEN Changzheng, XU Yanjun, et al. Integrated primary and tertiary mirror components from common base line of off-axis TMA space camera[J]. Optics and Precision Engineering,2015,26(6):1612-1619.

    [12] 王孝坤, 薛栋林, 张学军. 大口径非球面系统的共基准加工与检验[J]. 光学与精密工程,2018,26(4):743-748.

    WANG Xiaokun. XUE Donglin, ZHANG Xuejun. Fabrication and testing of large aspheric system based on common reference[J]. Optics and Precision Engineering,2018,26(4):743-748.

    [13] 潘君骅. 光学非球面设计、加工与检测[M]. 苏州: 苏州大学出版社, 2004: 130-135.
    [14] 孙权社, 赵发财, 陈坤峰, 等. 空间激光通信系统离轴天线设计研究[J]. 红外与激光工程,2015,44(8):2501-2505. doi: 10.3969/j.issn.1007-2276.2015.08.044

    SUN Quanshe, ZHAO Facai, CHEN Kunfeng et al. Design of off-axis optical antenna for space optical communications[J]. Infrared and Laser Engineering,2015,44(8):2501-2505. doi: 10.3969/j.issn.1007-2276.2015.08.044

    [15] 胡贞, 姜会林, 佟首峰, 等. 空问激光通信终端ATP技术与系统研究[J]. 兵工学报,2011,32(6):752-757.

    HU Zhen, JIANG Huilin, TONG Shoufeng et al. Research on ATP system technology of laser communication terminal in space[J]. Acta Armamentaii,2011,32(6):752-757.

    [16] 叶井飞. 光学自由曲面的表征方法与技术研究[D]. 南京: 南京理工大学, 2016.

    YE Jingfei. Research on the method and technique for characterizing freeform optical surface[D]. Nanjing: Nanjing University of Science & Technology, 2016.

    [17] 黎发志, 郑立功, 闫锋, 等. 自由曲面的CGH光学检测方法与实验[J]. 红外与激光工程,2012,41(4):1052-1056. doi: 10.3969/j.issn.1007-2276.2012.04.040

    LI Fazhi, ZHENG Ligong, YAN Feng, et al. Optical testing method and its experiment on freeform surface with computer-generated hologram[J]. Infrared and Laser Engineering,2012,41(4):1052-1056. doi: 10.3969/j.issn.1007-2276.2012.04.040

    [18] 赵猛, 颜昌翔, 吴从均. 激光通信地面测试终端间隔离度的仿真分析[J]. 中国光学,2020,13(3):472-481.

    ZHAO Meng, YAN Changxiang, WU Congjun. Simulation analysis of isolation between laser communication ground test equipments[J]. Chinese Optics,2020,13(3):472-481.

  • 期刊类型引用(5)

    1. 周晟刚. 通风机械仪表盘在复杂背景环境中视觉故障检测与定位研究. 计算机测量与控制. 2024(03): 106-111+145 . 百度学术
    2. 吴来新,王志. 基于机器学习的眼位测定仪视线跟踪误差自动复合补偿技术. 自动化与仪器仪表. 2024(10): 370-374 . 百度学术
    3. 洪祥,张海越,宋骐. 基于图像识别的AH36钢激光焊缝节点定位技术研究. 计算机测量与控制. 2023(11): 299-305+314 . 百度学术
    4. 付晓艳,白雪,魏博凝. 基于机器学习的物流自动分拣机器人定位精度检测算法. 自动化与仪器仪表. 2022(08): 245-248+252 . 百度学术
    5. 刘翰林,张荣福. PCB焊点及芯线端头识别. 光学仪器. 2022(05): 42-52 . 百度学术

    其他类型引用(3)

图(6)  /  表(4)
计量
  • 文章访问数:  1529
  • HTML全文浏览量:  824
  • PDF下载量:  141
  • 被引次数: 8
出版历程
  • 收稿日期:  2020-06-29
  • 修回日期:  2020-08-02
  • 录用日期:  2020-08-07
  • 网络出版日期:  2020-12-28
  • 刊出日期:  2021-01-14

目录

/

返回文章
返回