高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

一种激光三角测量的标定方法及误差分析

王振宇 李烨 郁丰 史海青

引用本文:
Citation:

一种激光三角测量的标定方法及误差分析

    作者简介: 王振宇(1989-), 男, 硕士研究生, 主要研究方向为控制工程、导航制导与控制.
    通讯作者: 郁丰, yufeng@nua.edu.cn
  • 基金项目:

    国家自然科学基金资助项目 61203197

  • 中图分类号: V448.25+1

Calibration method and error analysis of laser triangulation measurement

    Corresponding author: YU Feng, yufeng@nua.edu.cn ;
  • CLC number: V448.25+1

  • 摘要: 为了解决激光三角测量法在标定时,不能直接准确地测量激光器端口与标定物的距离的问题,采用了一种基于距离差进行标定的方法。该方法利用标定物移动的距离作为标定输入,改进了测量系统的标定方法;采用高斯-牛顿迭代法计算测量系统的参量,分析了测量系统在标定时可能产生的激光像点提取误差和非垂直误差。结果表明,测量距离的精度可达到4.000mm。该方法能够准确标定激光三角测量系统的参量。
  • Figure 1.  Vertical laser triangulation measurement

    Figure 2.  Relationship between gray value and cross section width

    Figure 3.  The diagram of the moving direction of calibration object

    Figure 4.  Experimental device of laser triangulation distance measurement

    Figure 5.  Moving process of laser line on camera imaging plane

    Figure 6.  Relationship between calibration distance error and α

    Figure 7.  Laser triangulation measurement with non-vertical error

    Table 1.  Recalculate of calibration distance

    number of moves distance between laser beam port and calibration object/mm actual moving distance of calibration object /mm
    1 430.993 0.000
    2 477.480 46.487
    3 526.487 49.007
    4 578.472 51.985
    5 628.450 49.978
    6 679.815 51.365
    7 732.246 52.431
    8 782.891 50.645
    9 831.363 48.472
    10 878.855 47.492
    下载: 导出CSV

    Table 2.  Two recalculation of calibration distance

    number of moves check distance with no improvement method/mm the actual moving distance without improved method/mm check distance with improvement method/mm two times calculation calibration distance/mm
    1 370.653 0.000 430.993 60.340
    2 414.078 43.425 477.480 63.402
    3 460.716 46.638 526.487 65.771
    4 511.181 50.465 578.472 67.291
    5 560.690 49.509 628.450 67.760
    6 612.618 51.928 679.815 67.197
    7 666.752 54.134 732.246 65.494
    8 720.159 53.470 782.891 62.732
    9 772.333 52.174 831.363 59.030
    10 824.489 52.156 878.855 54.366
    下载: 导出CSV

    Table 3.  Comparison of actual distance and measurement distance for α=5°

    number of moves actual distance/mm measuring distance/mm
    1 450.0 451.723
    2 500.0 501.914
    3 550.0 552.105
    4 600.0 602.296
    5 650.0 652.487
    6 700.0 702.678
    7 750.0 752.869
    8 800.0 803.059
    9 850.0 853.250
    10 900.0 903.441
    下载: 导出CSV
  • [1]

    LI D M. Dynamics and control of space robot in self-assembling on orbit[D]. Harbin: Harbin Institute of Technology, 2012: 1-8(in Chinese).
    [2]

    LI Y, CHANG D P. The development of orbital servicing technology in space[J]. Ordnance Industry Automation, 2012, 31(5):79-82(in Chinese).
    [3]

    ZENG C. Research on motion and mission punning methods for on-orbit servicing space manipulator[D]. Dalian: Dalian University of Technology, 2013: 2-18 (in Chinese).
    [4]

    FENG C Z. Target depth detection based on the linear laser triangulation [D]. Qingdao: Ocean University of China, 2012: 10-16(in Chinese).
    [5]

    ZHU W J, JIAO K H, XU C G, et al. Extraction of laser stripe center line based on genetic algorithm and NURBS interpolation[J]. Journal of Beijing Institute of Technology, 2008, 17(2):143-147(in Chinese). 
    [6]

    ZHANG X J, LI Z K, WANG X Z, et al. Research of 3-D point cloud data registration algorithms based on feature points and improved ICP [J]. Transducer and Microsystem Technologies, 2012, 31(9):116-118(in Chinese).
    [7]

    ZHANG G J. Vision measurement [M]. Beijing: Beijing Science Press, 2008:175-201(in Chinese).
    [8]

    ZHANG Z Y. Flexible camera calibration by viewing a plane from unknown orientations. international conference on computer vision[C]//IEEE International Conference on Computer Vision, 1999. New York, USA: IEEE, 1999: 666-673.
    [9]

    SCARAMUZZA D, MARTINELLI A, SIEGWART R. A flexible technique for accurate omnidirectional camera calibration and structure from motion [C]//IEEE International Conference on Computer Vision Systems. New York, USA: IEEE, 2006: 45-52.
    [10]

    ZHANG Z. A flexible new technique for camera calibration[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2000, 22(11):1330-1334. 
    [11]

    ZHANG Y, WANG Z. A flexible camera calibration method for computer visual 3-D reconstruction system[C]//International Confe-rence on Signal Processing. New York, USA: IEEE, 2006: 1580-1583.
    [12]

    REN W M, SUN P M, WAGN Y L, et al. A calibration method for laser displacement system based on triangulation[J]. Optical Technology, 1997, 23(3):10-12(in Chinese). 
    [13]

    DENG C M. Real time filtering algorithm for laser center based on image processing [J].Journal of Huazhong University of Science & Technology, 1999, 27(8):10-12 (in Chinese). 
    [14]

    PU J Y. Design of laser triangulation measurement system based on PSD [D].Wuhan: Huazhong University of Science and Technology, 2008: 6-42(in Chinese).
    [15]

    FENG J C. The design of a real-time and high-precision laser triangulation system [D]. Nanjing: Nanjing University of Science & Technology, 2008: 16-22(in Chinese).
  • [1] 周莉萍赵斌李柱 . 无衍射光束在激光三角测量系统中的应用研究. 激光技术, 1998, 22(1): 22-25.
    [2] 亢甲杰张福民曲兴华 . 激光雷达坐标测量系统的测角误差分析. 激光技术, 2016, 40(6): 834-839. doi: 10.7510/jgjs.issn.1001-3806.2016.06.013
    [3] 梁琨马泳程飞王宏远 . 基于布里渊激光雷达测量大气声速误差分析. 激光技术, 2009, 33(1): 94-96.
    [4] 陈曼龙 . 机器视觉螺纹测量的误差分析. 激光技术, 2014, 38(1): 109-113. doi: 10.7510/jgjs.issn.1001-3806.2014.01.024
    [5] 黄进袁晓东王成程贾怀庭秦兴武邓武徐冰唐晓东 . 利用测量显微镜开展激光束角漂实验研究. 激光技术, 2005, 29(3): 328-331.
    [6] 龙龙李宗峰 . 基于激光位移传感器的3维位置测量算法研究. 激光技术, 2017, 41(4): 531-536. doi: 10.7510/jgjs.issn.1001-3806.2017.04.015
    [7] 罗铂靓王成程杜惊雷马驰 . 用二元位相光栅进行激光束叠加的研究. 激光技术, 2006, 30(4): 381-384,389.
    [8] 王红敏罗自赢李蕾彭俊玲 . 平行度误差激光准直法测量技术研究. 激光技术, 2022, 46(5): 674-679. doi: 10.7510/jgjs.issn.1001-3806.2022.05.016
    [9] 刘辉张春波卢进南张锦茹 . 基于激光三角法的煤矸石装载体积测量方法. 激光技术, 2021, 45(5): 630-635. doi: 10.7510/jgjs.issn.1001-3806.2021.05.016
    [10] 王正家盛文婷解家月何涛 . 基于线性激光三角法的圆柱对象定位测量研究. 激光技术, 2019, 43(4): 464-470. doi: 10.7510/jgjs.issn.1001-3806.2019.04.006
    [11] 李秋实罗洪胡永明倪明 . 三角波形激光脉冲在液体中致声的特性分析. 激光技术, 2010, 34(3): 409-412. doi: 10.3969/j.issn.1001-3806.2010.03.035
    [12] 林宏董天临马泳梁速 . 基于边缘探测技术的布里渊散射测量误差分析. 激光技术, 2008, 32(1): 30-32,36.
    [13] 高亚丽王存山张梅熊党生 . 镁合金激光表面熔凝技术分析. 激光技术, 2009, 33(4): 362-365. doi: 10.3969/j.issn.1001-3806.2009.04.008
    [14] 林静焕戴勇林旭焕许敏界 . 激光扫描球形偏心引起误差的理论分析. 激光技术, 2014, 38(6): 813-816. doi: 10.7510/jgjs.issn.1001-3806.2014.06.019
    [15] 陈玉宝王箫鹏步志超王一萌郭泽勇熊峰王宣 . 超大城市试验气溶胶激光雷达标定及结果分析. 激光技术, 2022, 46(4): 435-443. doi: 10.7510/jgjs.issn.1001-3806.2022.04.001
    [16] 王诗昊向阳申屠功伟 . 外鼓式热敏CTP制版机光学系统设计. 激光技术, 2023, 47(1): 98-102. doi: 10.7510/jgjs.issn.1001-3806.2023.01.015
    [17] 王绍龙王阳恩陈奇陈善俊 . 激光诱导击穿光谱技术定量分析原油金属元素. 激光技术, 2015, 39(1): 104-108. doi: 10.7510/jgjs.issn.1001-3806.2015.01.021
    [18] 樊玉杰周建忠黄舒卫登辉王敏 . MEMS金属微构件的激光微喷丸强化技术分析与展望. 激光技术, 2010, 34(5): 665-669. doi: 10.3969/j.issn.1001-3806.2010.O5.024
    [19] 田国周欧群飞钟鸣叶大华吕百达 . 2kJ高能量钕玻璃激光器的热管理技术分析. 激光技术, 2007, 31(3): 253-256.
    [20] 马泳林宏冀航董天临 . 基于边缘探测技术的激光单稳频指标分析. 激光技术, 2007, 31(5): 469-472.
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  3142
  • HTML全文浏览量:  1818
  • PDF下载量:  272
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-07-27
  • 录用日期:  2016-07-28
  • 刊出日期:  2017-07-25

一种激光三角测量的标定方法及误差分析

    通讯作者: 郁丰, yufeng@nua.edu.cn
    作者简介: 王振宇(1989-), 男, 硕士研究生, 主要研究方向为控制工程、导航制导与控制
  • 1. 南京航空航天大学 航天学院, 南京 210016
  • 2. 航天恒星科技有限公司 导航事业部, 北京 100086
  • 3. 中国人民解放军 91362部队, 舟山 316200
基金项目:  国家自然科学基金资助项目 61203197

摘要: 为了解决激光三角测量法在标定时,不能直接准确地测量激光器端口与标定物的距离的问题,采用了一种基于距离差进行标定的方法。该方法利用标定物移动的距离作为标定输入,改进了测量系统的标定方法;采用高斯-牛顿迭代法计算测量系统的参量,分析了测量系统在标定时可能产生的激光像点提取误差和非垂直误差。结果表明,测量距离的精度可达到4.000mm。该方法能够准确标定激光三角测量系统的参量。

English Abstract

    • 大型航天器结构复杂、功能众多,其在轨检修、维护一直是航天界的焦点问题,是进一步提高大型航天器的可靠性与安全性的关键[1]。“哥伦比亚”号航天飞机在归途中解体,事后分析表明,失事原因之一是覆盖在航天飞机外层的一片隔热瓦脱离引起飞机外表融化解体。因此,对于大型航天器表面结构的在轨检测是保障航天器在轨正常运行的有效途径之一[2]。但是有航天员参与的在轨检测面临诸多困难;空间机械臂曾被应用于航天飞机的隔热瓦损伤在轨检测,具备精确操作能力和视觉识别能力,而空间机械臂也存在作用范围有限,不能全方位地对航天器进行观察和监视的不足[3]。因此,空间自由飞行机器人在轨服务中具有明显优势。美国国家航空航天局(National Aeronautics and Space Administration,NASA)约翰逊空间中心的工程师设计开发了一颗可以在航天器外自由飞行的微小航天器(miniature autonomous extravehicular robotic camera, Mini AERCam),旨在帮助宇航员及地面成员在任务执行期间观察航天器外部情况。该Mini AERCam采用的是惯性/全球定位系统(global positioning system,GPS)导航系统进行导航,但是精度有限,并且不适合在没有全球卫星导航系统(global navigation satellite system,GNSS)信号的高轨使用。

      区别于AERCam使用差分GPS,本文中提出利用线激光三角测量装置提取目标航天器表面结构点云信息[4-5],并与已知的外形结构数据匹配来获取两者之间的相对位姿关系[6]。其中线激光三角测量装置中的结构参量需要精确标定,才能实现精确可靠的测量。

      对于激光三角测量系统,常见的标定方法有拉丝法、锯齿靶法、微分法、已知直线与图像点匹配法、基于3维靶标法和本文中采用的基于自由移动标定物的激光测量标定方法[7-11]。这些方法在标定过程中都需要测量激光器端口与标定物间的距离作为系统的标定输入,但是激光器端口的位置无法准确确定,所以这个距离的测量值存在较大的误差[12]。为了提高测量系统的精度,本文中提出了基于距离差的标定方法,建立了基于距离差的标定模型,回避了绝对距离的测量,利用容易获取的标定物的移动距离作为标定输入,标定过程实施简单、精度较高。

      本文中的激光三角测量系统,测量距离在1m左右,测距精度为4.000mm。激光三角测量法作为一种传统的激光测距方法,有着结构简单、测量频率高、数据可实时更新、体积小、质量轻的优点[13-15]。在实验中使用的基于距离差的标定方法有效提高了激光三角测量的精度,测距范围及实时性均满足本文中的要求。

    • 在标定物表面上,用一束激光以一定的角度照射,激光在标定物表面发生反射,在另一个角度用成像系统对激光反射光进行汇聚成像。由于入射和反射光构成一个三角形,如果标定物沿着激光束的方向发生位移时,则相机成像平面上的像点也会发生相应位移,通过物像位移之间的关系可以计算出被测物面各点的坐标。

      按入射激光光束与被测物面的角度不同,可以把激光三角测量法分为斜射式和直射式两种。斜入射和直入射都可以实现对标定物进行非接触测量。一方面,直入射较斜入射式体积更小,光强集中,光斑变形小;另一方面,直入射式与斜入射式的测量模型在本质上是一致的,所以采用直入射法进行标定方法的研究。

    • 垂直激光三角测量原理如图 1所示。图中,s为激光器光束端面A与基准面R的垂直距离,d为端面A与相机焦点O′的距离,d0OO′的距离,f是相机焦距,可以使用基于2-D平面靶标的摄像机标定方法标定得到,Δ为被测物体距离基准面的距离,θ是激光束与OO′的夹角,δ是在成像平面上成像点偏移的距离,可以求得公式:

      Figure 1.  Vertical laser triangulation measurement

      $\mathit{\Delta = }{\mathit{d}_{\rm{0}}}/\left( {{\mathit{\delta }^{ - 1}}f\sin \mathit{\theta + }{\rm{cos}}\mathit{\theta }} \right) $

      (1)

      式中,虽然OO′的长度d0、镜头的焦距f、激光束与镜头的夹角θ均为定值,在测量时是固定不变的,但是在测距系统装配时很难获得其准确值。因此,为了让激光三角测量系统精确地测量位移,在激光三角测量系统装配完毕后需要对这些参量进行精密标定。

    • 图 1所示,由于激光器发出激光束端面A无法准确判定,所以在标定时端面A与被测物体的距离无法直接测得。本文中提出了一种根据距离差进行标定的新方法,根据(1)式得到如下公式:

      $\begin{array}{l} {\mathit{\Delta }_{i + 1}} - {\mathit{\Delta }_i} = {\mathit{\delta }_{i + 1}}/\left[ {{d_0}^{ - 1}\left( {f\sin \mathit{\theta + }{\mathit{\delta }_{i + 1}}\cos \mathit{\theta }} \right)} \right] - \\ \;\;\;\;\;\;\;\;\;\;\;\;{\mathit{\delta }_i}/\left[ {{d_0}^{ - 1}\left( {f\sin \mathit{\theta + }{\mathit{\delta }_i}\cos \mathit{\theta }} \right)} \right] \end{array} $

      (2)

      式中,i和(i+1)为两次标定物移动的次数,δi, δi+1为像点在成像平面上的横坐标,可以通过对所拍图像的处理得到,Δi+1Δi为任意两次标定物移动的距离,可以直接通过测量得到。令a=f sinθ/d0b=cosθ/d0y=Δi+1Δixi+1=1/δi+1xi=1/δi,得到公式:

      $y = {\left( {a{x_{i + 1}} + b} \right)^{ - 1}} - {\left( {a{x_i} + b} \right)^{ - 1}} $

      (3)

      由于Δi+1Δiδiδi+1的值之间是非线性关系,本文中利用非线性回归分析的方法来求ab的值,并由公式θ=arctan(af-1b-1)和d0=b-1cosθ计算得到θd0的值,从而完成激光三角测量系统的准确标定。

      本文中采用高斯-牛顿迭代法来计算ab的估计值。首先,选取a0b0作为ab的初始估计值,将(3)式在点(a0, b0)处展开成泰勒级数,并忽略2阶以上的高阶项,得到如下方程:

      $y - \mathit{\hat y = m}{\rm{d}}\mathit{a + n}{\rm{d}}\mathit{b} $

      (4)

      式中,da, db为参量a, b的微分量,$\mathit{\hat y}$为函数y泰勒展开式的第1项,mn为函数y泰勒展开式1阶项的系数;其中:

      $\left\{ \begin{array}{l} m = \frac{{ - {x_{i + 1}}}}{{{{\left( {a{x_{i + 1}} + b} \right)}^2}}} + \frac{{{x_i}}}{{{{\left( {a{x_i} + b} \right)}^2}}}\\ n = \frac{{ - 1}}{{{{\left( {a{x_{i + 1}} + b} \right)}^2}}} + \frac{1}{{{{\left( {a{x_i} + b} \right)}^2}}}\\ \mathit{\hat y = }\frac{1}{{{a_0}{x_{i + 1}} + {b_0}}} - \frac{1}{{{a_0}{x_i} + {b_0}}} \end{array} \right. $

      (5)

      根据(4)式给出标定公式的矩阵描述:

      $\mathit{\boldsymbol{Z = HX + \varepsilon }} $

      (6)

      ${\rm{其中, }}\left\{ \begin{array}{l} \mathit{\boldsymbol{Z = }}{\left[ {{y_1} - {{\mathit{\hat y}}_1}{\rm{\;\;\;\;}}{\mathit{y}_2} - {{\mathit{\hat y}}_2}{\rm{\;\;\;\;}}{\mathit{y}_3} - {{\mathit{\hat y}}_3}{\rm{\;\;\;\;}} \cdots } \right]^{\rm{T}}}\\ \mathit{\boldsymbol{H = }}{\left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}}&{{m_3}}& \cdots \\ {{n_1}}&{{n_2}}&{{n_2}}& \cdots \end{array}} \right]^{\rm{T}}}\\ \mathit{\boldsymbol{X = }}{\left[ {{\rm{d}}\mathit{a}{\rm{\;\;\;\;d}}\mathit{b}} \right]^{\rm{T}}} \end{array} \right. $

      (7)

      式中,ε为误差项矩阵,(6)式可使用加权最小二乘法进行计算,当计算出X后,执行ai+1=ai+dabi+1=bi+db,可迭代固定次数或者改正量小于某一阈值, 即可迭代计算,直到第n次计算得到比较理想的参量估计值anbn,并以此次的计算值作为ab的标定值。

    • 本文中的测量系统使用的是线激光光源,线激光光源照射到被测物面的光带是有一定宽度的,因而它在相机上的成像并不是一个像素的宽度,激光束在成像平面上的成像宽度是占用几个像素甚至几十个像素,如图 2所示。所以在图像处理的过程中就涉及到激光线中心点提取的问题。

      Figure 2.  Relationship between gray value and cross section width

      一般认为,激光光带的灰度分布是成高斯对称分布的。但由于受到外界环境和传感器内部噪声影响,接收的光斑波形发生改变。根据这个特性,本文中采用重心法求激光带中心点坐标[2]。首先,确定一个光带边界的灰度阈值T,然后在图像的每一行上判断出灰度值大于T的像素,利用重心公式计算光心的位置。假设在图像的一行上,灰度值大于阈值T的像素点坐标为uj,其灰度值为gj,其中,j=1, 2, 3, …, M(M为大于阈值T的像素的个数)。利用重心公式可以得到光心的位置u为:

      $\mathit{u = }\sum\limits_{j = 0}^M {{g_j}{u_j}} /\sum\limits_{j = 0}^M {{g_\mathit{j}}} $

      (8)

      利用这种方法得到的光带中心是亚像素级的,该算法在图像背景较暗的情况下具有很好的效果。

    • 测量系统标定过程中,标定物移动的方向与激光束的方向不会完全一致,如图 3所示。当标定物移动方向与激光束方向夹角为α时,移动距离BC′=BC,有效移动距离为BD,其长度小于BC′,使得标定物的实际移动距离Δ2′-Δ1小于标定计算中使用的移动距离Δ2Δ1,从而产生误差。本文中将标定物移动方向与激光束方向之间的误差角α导致的标定误差称为非垂直误差。

      Figure 3.  The diagram of the moving direction of calibration object

      令标定物移动的距离为l,则标定物移动方向与激光束方向夹角α和标定距离误差e的关系如下所示:

      $\mathit{e = l - l} \times {\rm{cos}}\mathit{\alpha } $

      (9)

      对非垂直误差的分析在第2节中详细分析说明。

    • 首先搭建激光三角测量系统的实验装置,本文中研究的测量系统硬件系统包括:线激光器、红外滤镜、标定物和相机。激光器和相机固定在工作台上,如图 4所示。

      Figure 4.  Experimental device of laser triangulation distance measurement

      本系统需要采集多组相对深度数据,取一个标定模板放置于线激光器正前方,使标定物与激光束的方向垂直,沿激光束方向多次移动标定物,移动距离分别为30.000mm, 50.000mm, 100.000mm, 150.000mm, 200.000mm, 250.000mm, 300.000mm, 350.000mm, 400.000mm, 450.000mm,共移动10次。随着标定物移动距离的增加,激光束在相机成像平面上的成像位置在不断移动,如图 5所示。

      Figure 5.  Moving process of laser line on camera imaging plane

    • 相机的焦距经标定得到f=15.332mm,利用第1.2节中的原理处理所拍图像,得到a=0.00520, b=0.00159,由此计算可得测量系统的参量θ=12.041°, d=128.316mm,激光器与基准面距离s=601.560mm。

      同样采用距离差的方法来验证标定的准确性,重新在实验中每隔50.000mm移动一次标定物并采集一组图像。通过分析图像得到的参量δ和标定实验中计算得到的测量系统的参量,利用(1)式计算得到激光束端口与标定物距离,从而检验这些参量是否准确,见表 1

      Table 1.  Recalculate of calibration distance

      number of moves distance between laser beam port and calibration object/mm actual moving distance of calibration object /mm
      1 430.993 0.000
      2 477.480 46.487
      3 526.487 49.007
      4 578.472 51.985
      5 628.450 49.978
      6 679.815 51.365
      7 732.246 52.431
      8 782.891 50.645
      9 831.363 48.472
      10 878.855 47.492

      表 1中的数据可以看出,在400.000mm~900.000mm的范围内测距,实验中由非垂直误差以及实验仪器本身的精度等一系列因素引起的误差优于4.000mm,满足要求。上述实验数据表明,本文中提出的实验标定方法的测量精度能够达到毫米级。

      如果采用未改进的测量系统的标定方法,直接测量激光器端口与标定物的距离,假定激光器端口与标定物的距离测量时误差为比实际距离减少5.000mm,利用前面计算得到的激光器与基准面距离为已知值,通过第1.2节中的原理处理所拍图像, 并采用最小二乘法对(1)式进行计算的得到测量系统的参量θ=10.725°,d=101.075mm,激光器与基准面距离s=533.648mm。

      采用前面改进方法中验证标定的准确性所采集的数据进行复核计算。利用(1)式计算得到激光束端口与标定物距离,与改进方法中复核计算的标定距离作比较,从而检验这些参量是否准确,如表 2所示。

      Table 2.  Two recalculation of calibration distance

      number of moves check distance with no improvement method/mm the actual moving distance without improved method/mm check distance with improvement method/mm two times calculation calibration distance/mm
      1 370.653 0.000 430.993 60.340
      2 414.078 43.425 477.480 63.402
      3 460.716 46.638 526.487 65.771
      4 511.181 50.465 578.472 67.291
      5 560.690 49.509 628.450 67.760
      6 612.618 51.928 679.815 67.197
      7 666.752 54.134 732.246 65.494
      8 720.159 53.470 782.891 62.732
      9 772.333 52.174 831.363 59.030
      10 824.489 52.156 878.855 54.366

      从表中的数据可以看出,由于测量激光器端口与标定物的距离减少了5.000mm,最终标定出来的测量系统的测量误差很大, 其测量出的移动距离精度也略低于改进方法的精度; 同时, 与改进后的测量方法测量标定物的实际距离相差了50.000mm以上,而改进后的测量系统测量标定物的距离经测量验证与实际距离误差很小,所以未改进的方法测量出的距离误差比较大。上述实验数据表明,本文中提出的实验标定方法的测量精度比与未改进的测量方法准确了很多。

    • 第1.3节中已经分析了当标定物移动方向与激光束方向存在误差角α时,标称移动距离比实际移动距离要大。

      假设每次标定物移动的距离l=50.000mm,随着α的增加,标定距离误差如图 6所示。

      Figure 6.  Relationship between calibration distance error and α

      图 6可以看出,标定距离误差与α大小呈非线性变化。同时开展仿真研究,探索非垂直误差对最终测量结果的影响规律。令相机参量为f=15.000mm,三角测量系统参量θ=12.000°, d=120.000mm,激光器与基准面距离设定为s=564.560mm。假设α=5.000°时,初始标定物距离为450.000mm时,每隔50.000mm采集一次图像,即实际标定的距离为450.000mm, 499.809mm, 549.619mm, 599.429mm, 649.239mm, 699.049mm, 748.859mm, 798.669mm, 848.479mm, 898.289mm,经标定后得到的三角测量系统的参量为θ=12.0001°,d=120.46mm,激光器与基准面距离s=566.718mm。

      图 7为含有非垂直误差的激光三角测量系统。利用含有非垂直误差的标定参量进行测量计算,结果与实际距离的比较如表 3所示。

      Figure 7.  Laser triangulation measurement with non-vertical error

      Table 3.  Comparison of actual distance and measurement distance for α=5°

      number of moves actual distance/mm measuring distance/mm
      1 450.0 451.723
      2 500.0 501.914
      3 550.0 552.105
      4 600.0 602.296
      5 650.0 652.487
      6 700.0 702.678
      7 750.0 752.869
      8 800.0 803.059
      9 850.0 853.250
      10 900.0 903.441

      由此可见,当激光三角测量系统存在非垂直误差时,每次测量距离的间隔变为50.191mm,随着测量距离的增加,测量距离与实际的误差也在增加。这是由于在标定的时候,存在非垂直误差,使得计算中将每次标定物的实际移动距离49.809mm作为50.000mm进行计算,从而使得到的激光三角测量系统测量的距离比实际距离要大。在测量时,标定物移动距离50.000mm,测量得到的距离间隔为50.191mm,即测量得到的距离间隔比实际距离间隔要大。

    • 本文中提出了利用距离差对激光三角测量系统进行标定的方法,解决了三角测量法在标定时,无法直接准确地测量激光器与标定物距离的问题。同时,作者通过仿真对激光像点的提取误差和非垂直误差进行了分析。实验结果表明,该方法结构简单并能够准确地对激光三角测量系统进行标定,系统精度可以达到4.000mm,为基于惯性/激光三角测量的航天器相对导航的研究提供了基础。

参考文献 (15)

目录

    /

    返回文章
    返回