一种激光摄影测量系统及相机标定方法
一种激光摄影测量系统及相机标定方法
【技术领域】
[0001] 本发明设及一种激光摄影测量系统及相机标定方法，属于测量装置标定方法技术 领域。
【背景技术】
[0002] 随着计算机技术的发展，摄影测量与Ξ维重构技术在机器人视觉、工业制造、工程 建筑、文物保护和医学图像分析等领域得到广泛应用。相比于利用Ξ坐标测量仪等接触式 测量仪器对实物外表面进行测量，激光扫描等非接触式测量利用非接触的光电方法对曲面 的Ξ维形貌进行快速测量，具有非接触、速度快、精度高、应用广泛等优点，在工业工程等领 域有着一定的应用前景。在测量过程中，基于视觉传感器(如摄像机)获取的图像，建立二维 像点与Ξ维物点的对应关系，利用图像中的二维数据反算出对应的Ξ维物体的形貌和相关 位置等信息。
[0003] 要想获得高精度的测量结果，首先需要设计的测量系统具有较高的精度;其次建 立二维像点与Ξ维物点的对应关系，也即摄像机标定技术，又包括摄像机建模与摄像机标 定方法。相机标定主要是标定相机固有的内部几何参数和相机相对物方坐标之间的位置关 系等外部参数，其标定结果的好坏直接决定测量得到的Ξ维坐标点的精确程度，尤其是对 度量领域，因此需要研究相机的标定方法和研制相应的装置。
[0004] 目前相机标定技术根据解算相机内部几何参数和外部位置关系参数时所使用的 方法，将其分为线性方法、非线性方法和自标定方法。线性方法如直接线性变换方法(DLT变 换），根据环境点与相机图像对应点的关系，通过线性方法求得线性模型参数，进而建立线 性模型来描述相机的成像过程，该方法估计相机参数过程效率较高，但是没有考虑成像过 程中的非线性崎变问题，因此精度较低。非线性方法对提高摄像机标定精度提供了一条很 好的途径，非常适合于多参数的非线性崎变模型的标定，在能提供较好初值的情况下，能更 好地模拟和补偿各种像差，可W较快地收敛，得到高精度优化结果，但是它对标定时摄像机 模型参数的初始值要求较高，计算量大，计算速度慢。自标定方法在场景未知和相机运动参 数未知的情况下也可W标定相机的内外参数，在标定过程中仅依靠多幅图像对应点之间的 关系进行，其灵活性高、应用广泛，但是标定的过程主要是基于绝对二次曲线和绝对二次曲 面方法，此外还需要求解多元非线性方程，因此只适用于一些对精度要求不高的场合。
[0005] 造成目前所采用的标定方法存在上述问题的原因主要是其对应的测量系统较为 复杂。不同的测量系统需要采用不同的标定方法来标定，因此为了将标定方法简单化，尽可 能避免复杂标定方法造成的标定误差，首先需要使设计的测量系统简单，然后在针对设计 的测量系统设计标定方法。
[0006] 此外，针对大型的被测物，实际测量时往往是利用多台相机同时测量，而多台相机 之间没有相对位置关系，因此需要单独对每台相机进行标定，运样就增加了标定工作量，同 时还会增加标定误差，造成测量精度降低。所W在设计测量系统及标定方法时还应考虑同 时标定两台及W上的相机。
【发明内容】
[0007] 针对现有技术的不足，本发明提供一种激光摄影测量系统。
[0008] 本发明还提供上述一种激光摄影测量系统中相机的标定方法。
[0009] 本发明的技术方案如下：
[0010] -种激光摄影测量系统，包括测量装置，所述测量装置包括底座、支架及计算机； 所述支架设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与 线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相机，所述计算 机分别与相机、激光测距仪、线激光发射器电连接。
[0011] 优选的，所述支架的底部设置有水平分度盘。此设计的好处在于，当旋转支架时， 通过水平分度盘能够清楚直观地观察到支架旋转的角度。
[0012] 优选的，所述激光摄影测量系统还包括标定装置，所述标定装置包括物方标识点 和Ξ个物方板，Ξ个物方板的外形均为正方形，Ξ个物方板前后依次同轴平行设置且Ξ个 物方板的边长依次增大，在每个物方板上设置八个所述的物方标识点，八个物方标识点按 如下方式分布:在每条边的中点设置一个物方标识点，在相邻两条边的交点设置一个物方 标识点。
[001引优选的，所述立个物方板的尺寸分别为：100mm X 100mm X 10mm，200mm X 200mm X 10mm,300mmX300mmX10mm0
[0014] 优选的，所述物方标识点采用直径5mm的圆形薄片。
[0015] 优选的，所述标定装置还包括导轨，Ξ个物方板通过滑动支架设置在导轨上。此设 计的好处在于，将Ξ个物方板通过底部的滑动支架安装在导轨上，方便调节Ξ个物方板之 间的间距。
[0016] -种激光摄影测量系统的相机标定方法，包括W下步骤，
[0017] (1)将测量装置放置在标定装置的正前方，建立物方坐标系，并记录各物方标识点 在物方坐标系下的物方坐标(Xci，Yci，Zci)(i = l，2·-24):
[0018] (2)通过计算机控制相机对Ξ个物方板拍照，得到物方标识点的图像数据，通过计 算机内部执行程序识别出图像数据中物方标识点在相机的像坐标系下的实际像坐标(Χ01, y〇i)(i = l ,2---24)；
[0019] (3)对每台相机的9个相机参数赋初值，9个相机参数包括相机的线坐标(Xs，Ys， Zs)、角坐标(知W，K)、像主点相对于影像中屯、的位置X0、yo和焦距f，各参数初始化赋值如下：
[0022] 利用相机角坐标林，W，K)的初始值确定旋转矩阵R:
[0023]
[0024] (4)根据公式(2)计算各物方标识点(沿，打，21)。= 1，2-，24)在像坐标系下的理论 像坐标(xi，yi)(i = l，2···24)，
[0027] 然后根据旋转矩阵R、理论像坐标^1，7〇(1 = 1，2-，24)和实际像坐标^(^，7()〇(1 =1，2-'24)按照最小二乘法计算得到误差方程式系数矩阵心。=1，2-'24)，并利用物方标 识点的实际像坐标^〇1，7〇0。= 1，2，，，24)减去理论像坐标值^1，7〇(1 = 1,2，，，24)构成误 差矩阵k(i = l，2...24):
[002引 [lx, ]_y]T= [x0i-Xi,y0i_yi]T;
[0029] (5)根据各物方标识点的误差方程式系数矩阵Αια = 1，2···24)和误差矩阵^(i = 1，2···24)分别构造两大矩阵A和レA=[Al A2……A24]t，L=[Li L2……L24]t;
[0030] 然后构造法方程，X=(aTa)-1aTl=[Ms ΔΚ,么爲匈> Aw Δκ·Δ/Δ?0 Δ細 （雖
[0031] (6)求解法方程(3)便可从公式右端得到各相机参数的修正量，判断Ξ个角坐标（ 9，w，K)对应的修正量(Δ;0，Aw，Δκ)是否都小于规定的阔值Δ，若满足条件则计算结束;若 不满足条件，通过解算出来的9个相机参数的修正量修改相机参数值，然后重复步骤(3)~ (6);若迭代解算次数超过20次仍然不满足条件，则说明硬件摄影测量系统有问题。
[0032] 优选的，步骤（1)中，在物方板的中屯、建立物方坐标系ΧΥΖ，其中水平向右为X轴正 方向，竖直向上为Υ轴正方向，垂直于物方板为Ζ轴。
[0033] 本发明的有益效果在于：
[0034] (1)相比于目前已有的摄影测量系统的标定装置，本发明的标定装置由物方板和 物方标识点组成，结构简单，成本低、安装方便、操作简单、精度高。
[0035] (2)传统方法中将定位定姿数据转换到相机中屯、，如果存在两个或两个W上的相 机时，需要建立多个相机坐标系，给计算增加了难度，且不便于数据的归一化处理;本发明 中将坐标系建立在标定装置上，将两台相机 的位姿坐标都统一到物方坐标系下，便于计算， 减少了编程工作量，不再需要对两台相机解算的数据进行坐标变换处理。
[0036] (3)根据本发明相机标定方法编写的程序简单、易于操作、求解的精度较高，在图 像数据中识别物方标识点的像坐标时，快速、精确;迭代过程简单，整个标定过程不设及自 标定方法中由于二次曲线和二次曲面造成的多元非线性方程组。
【附图说明】
[0037] 图1为本发明中测量装置的结构示意图；
[0038] 图2为本发明中测量装置和标定装置的位置关系示意图；
[0039] 图3为本发明中测量装置的激光中屯、轴线垂直于物方板的关系图。
[0040] 其中：1、底座;2、水平分度盘;3、支架;4、相机;5、激光测距仪;6、线激光发射器;7、 相机;8、数据线;9、计算机；10、Ξ角架；11、测量装置；12、点激光；13、第一物方板；14、第二 物方板;15、第Ξ物方板;16、物方标识点；17、导轨;18、滑动支架。
【具体实施方式】
[0041] 下面通过实施例并结合附图对本发明做进一步说明，但不限于此。
[0042] 相机的标定，主要是标定相机在物方坐标系下的外方位元素和相机的内方位元 素。其中每台相机都有3个内方位元素和6个外方位元素，内方位元素包括像主点相对于影 像中屯、（即照片中屯、点）的位置xo、y〇W及镜头中屯、到影像面的焦距f (也称主距）。外方位元 素包括3个用于描述镜头中屯、相对于物方坐标系位置的线元素 Xs，Ys，Zs和3个用于描述影像 面在摄影瞬间空中姿态的角元素护、ω、k。激光测距仪点激光中屯、在物方坐标系中的坐标为 (0,0,Zcj)。
[0043] 实施例1:
[0044] 如图1和图2所示，一种激光摄影测量系统，该激光摄影测量系统包括测量装置和 标定装置，其中，测量装置包括底座1、水平分度盘2、支架3、两台相机4和7、线激光发射器6、 激光测距仪5、数据线8、计算机9。其中，相机4和相机7对称安装在激光测距仪5两侧，Ξ者都 安装在支架3上，计算机9通过数据线8与支架3内的控制设备连接。测量系统还包括Ξ脚架 10，底座1通过Ξ脚架10上的螺纹锁紧件与Ξ脚架10连接固定。
[0045] 相机4、相机7、线激光发射器6和激光测距仪5均由计算机9通过数据线8来传输控 制命令，并将获得的图像数据传递到计算机9中，方便数据的存储、计算和分析。
[0046] 计算机9装有相机4、相机7、线激光发射器6和激光测距仪5的驱动及控制程序，负 责控制发射激光、相机转动和拍摄照片等。
[0047] 计算机9还装有实现相机标定方法的执行程序，通过读取图像数据解算得到相机 位置和姿态坐标(位置坐标也即相机的线元素、姿态坐标也即相机的角元素，简称位姿坐 标)。
[0048] 如图2所示，标定装置是一个物方模型，该物方模型包括Ξ个用于布置物方标识点 16的物方板，Ξ个物方板包括第一物方板13、第二物方板14和第Ξ个物方板15,第一物方板 13和第二物方板14底部设置有滑动支架18,通过滑动支架1則尋第一物方板13和第二物方板 14放置导轨17上，第Ξ物方板15与导轨17的一端固定在一起(第Ξ物方板也可通过滑动支 架设置在导轨上）。第一物方板13、第二物方板14和第Ξ物方板15上均布置有8个物方标识 点16，均匀布置在每个物方板四条边上。
[0049] 其中，第一物方板13的尺寸(长宽厚)为lOOmmXlOOmmXlOmm，第二物方板14的尺 寸为200mmX200mmX10mm，第Ξ物方板15的尺寸为300mmX300mmX10mm，Ξ个物方板前后 依次平行放置，物方板的中屯、W及测量装置11的激光中屯、均在同一条线上，因此测量装置 11正对着物方模型获取物方标识点16的图像时，在一张照片上便能够同时显示所有物方标 识点的图像。
[0050] 物方标识点16是直径5mm的圆形薄片，为了与背景色有明显的区别，将物方标识点 16表面涂成黑色。
[0051] 第一物方板13和第二物方板14通过与其相连接的滑动支架18可在导轨17上滑动， 从而可W调节Ξ个物方板之间的距离。
[0052] 所述的物方模型上，在第Ξ物方板15的中屯、建立物方坐标系，水平向右为X轴正方 向，竖直向上为Y轴正方向，垂直于第Ξ物方板15为Z轴，方向遵循右手定则。
[0053]本实施例中的激光摄影测量系统结构简单、方便操作、精度高，是进行相机标定方 法的基础，同时高精度的激光摄影测量系统也是后续进行激光摄影测量而获得高精度测量 结果的保证。
[0化4] 实施例2:
[0055] 在实施例1的基础上，本实施例详细说明利用实施例1所述的激光摄影测量系统进 行相机标定方法的过程。
[0056] 本实施例提供一种激光摄影测量系统的相机标定方法，首先得到在物方坐标系下 各物方标识点16的坐标，然后调整测量装置11，使两台相机正对物方模型拍照，得到24个物 方标识点16的图像数据。通过计算机9上实现相机标定方法的执行程序读出各物方标识点 16的像坐标，利用物方标识点16的物方坐标与像坐标的对应关系，由执行程序计算得到两 台相机和激光测距仪点激光中屯、在物方坐标系中的位姿坐标及内方位元素参数。具体标定 方法如下：
[0057] 步骤1:确定各物方标识点16在物方坐标系下的坐标值(Xci，Yci，Zci) (i = 1，2… 24);
[0058] 步骤2:将测量装置11通过螺纹锁紧件固定在Ξ脚架10上，然后将计算机9通过数 据线8与测量装置相连，启动计算机9上控制相机与激光测距仪点激光中屯、的驱动及控制程 序，通过控制程序调整测量系统位姿，使激光测距仪轴线水平；
[0059] 步骤3:调整物方模型的位置，使其正对着测量装置11，第一物方板13在最前面；
[0060] 步骤4:如图2、图3所示，调整测量装置11的位姿，使点激光12与物方模型的第一物 方板13垂直。具体步骤如下：
[0061] (1)通过计算机9上的控制程序，打开测量装置11的激光测距仪5，激光测距仪5发 射的点激光12，通过肉眼观察使点激光12近似垂直于物方模型的第一物方板13，点激光12 打在第一物方板13的化位置，如图3所示；
[00创 （2似0102线为基准位置，通过旋转支架3，使0102绕点激光中屯、01点逆时针转过目1 到OiA的位置，通过激光测距得到化A的长度La;
[0063] (3) W化化线为基准位置，通过旋转支架3,使化化绕测量装置中屯、化点顺时针转过 θι至化iB的位置，通过激光测距得到化B的长度Lb;
[0064] (4)比较La与Lb的大小，如果La〉Lb，转动水平支架3使〇2右移向B靠近;如果La<Lb，转 动水平支架使化左移向A靠近；
[0065] (5)如果La = Lb，则点激光12在水平方向垂直于第一物方板13;否则重复步骤(2)~ (4)，最终使点激光12与物方模型的第一物方板13在水平方向上垂直；
[0066] (6)保证点激光12在竖直方向上与第一物方板13垂直的方法同上述步骤，只需要 将方向改为在竖直方向上微调点激光12的位置。
[0067] 步骤5:通过计算机9上的控制程序，打开相机4和相机7，对物方模型拍照，通过程 序分别获取各个物方标识点16在像坐标系下的实际像坐标(xoi，yoi) (i = 1，2· · · 24);
[0068] 步骤6:给每台相机的9个参数赋初值，9个参数包括相机的线坐标(Xs，Ys，Zs)、角坐 标(r，w，K)、像主点相对于影像中屯、（即照片中屯、点）的位置x〇、y〇和焦距f，各参数初始化赋 值如下：
[0069] 口二w='A=.v(尸.!',)=0，30</.:? 120;
[0070] Xs = ( Σ Xci) /24，Ys = ( Σ Yci) /24，Zs = ( Σ Zci) /24;
[0071] 步骤7:根据公式(1)，利用相机角坐标元素值界、*、1〇确定旋转矩阵R中各元素的值：
[0072]
[0073] 步骤8:利用基于空间后方交会原理编写的程序计算得到理想情况下物方标识点 在相机的像坐标系下的理论像坐标，根据公式(2)计算各物方标识点(Χι，Υι，Ζι)在像坐标系 下的理论像坐标(xi，yi)(i = l，2···24)，
[0076] 步骤9:然后根据旋转矩阵R、各物方标识点的理论像坐标(Χι，γι)α = 1，2···24^Ρ 实际像坐标(x()l，y()l)(i = l，2…24)按照最小二乘法计算误差方程式系数矩阵Al(i = l，2… 24);
[0077] 步骤10:利用物方标识点像坐标的实际值(xoi，y日i)减去理论像坐标值(xi，y i)构成 误差矩阵以。=1，2-'24):
[007引 [lx, ；Ly]T= [x0i-Xi，y0i-yi]T;
[0079] 步骤11:根据各物方标识点系数矩阵Αια = 1，2···24)和误差矩阵^α = 1，2···24) 分别构造成两个大矩阵4和1^:4=[41 Α2 ...... A24]T，L=[Li L2 ...... L24]T;
[0080] 步骤 12:构造法方程，X=(aTa)-iATl=[A^ ΔΚ Δ& Aw ΔκΔ/細ο Δ础，求解法方 程便可从公式右端得到相机各参数的修正量；
[0081 ] 步骤13:判断Ξ个角坐标作，W，Κ)对应的修正量(Δ口，Δ W，Δ Κ)是否都小于规定的 阔值A，如果满足条件则计算结束;如果不满足条件，通过解算出来的9个相机参数的修正 量修改相机参数值，然后重复步骤7~12;如果迭代解算次数超过20次仍然不满足条件，说 明硬件摄影测量系统有问题。
[0082] 步骤14:两台相机都按照上述步骤进行解算，不同之处在于物方标识点的像坐标 分别取自两台相机各自拍到的图像数据，最终解算得到两台相机的9个参数。由于激光测距 仪点激光中屯、位于两台相机的对称中屯、，因此其线坐标为两台相机线坐标的平均值。
[0083] 实际在拍摄照片进行相机参数解算时，为了避免随机误差，每个相机拍摄3~5张 照片，并解算每张照片的数据，通过分析数据，剔除偏差较大的数据，对剩下的几组数据取 平均值得到相机的9个参数。
[0084] 本发明在实施过程中借鉴了现有技术，限于篇幅，未对现有技术部分进行详细描 述;凡是本发明未提及的技术部分，均可W采用现有技术实现。
[0085] W上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技 术人员来说，本发明可W有各种更改和变化。凡在本发明的精神和原则之内，所做的任何修 改，等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种激光摄影测量系统，包括测量装置，其特征在于，所述测量装置包括底座、支架 及计算机;所述支架设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激 光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相 机，所述计算机分别与相机、激光测距仪、线激光发射器电连接。2. 如权利要求1所述的激光摄影测量系统，其特征在于，所述支架的底部设置有水平分 度盘。3. 如权利要求1所述的激光摄影测量系统，其特征在于，所述激光摄影测量系统还包括 标定装置，所述标定装置包括物方标识点和三个物方板，三个物方板的外形均为正方形，三 个物方板前后依次同轴平行设置且三个物方板的边长依次增大，在每个物方板上设置八个 所述的物方标识点，八个物方标识点按如下方式分布:在每条边的中点设置一个物方标识 点，在相邻两条边的交点设置一个物方标识点。4. 如权利要求3所述的激光摄影测量系统，其特征在于，所述三个物方板的尺寸分别 为：100mm X 100mm X 10mm，200mm X 200mm X 10mm，300mm X 300mm X 1 Omnin5. 如权利要求3所述的激光摄影测量系统，其特征在于，所述物方标识点采用直径5mm 的圆形薄片。6. 如权利要求3所述的激光摄影测量系统，其特征在于，所述标定装置还包括导轨，三 个物方板通过滑动支架设置在导轨上。7. -种如权利要求1-6任一项所述的激光摄影测量系统的相机标定方法，包括以下步 骤， (1) 将测量装置放置在标定装置的正前方，建立物方坐标系，并记录各物方标识点在物 方坐标系下的物方坐标(Xc^Ycd,Zed)(i = l，2···24); (2) 通过计算机控制相机对三个物方板拍照，得到物方标识点的图像数据，通过计算机 内部执行程序识别出图像数据中物方标识点在相机的像坐标系下的实际像坐标(XQ 1J01) (i = l,2···24); (3) 对每台相机的9个相机参数赋初值，9个相机参数包括相机的线坐标(XS，YS，ZS)、角 坐标0，W，4、:像主点相对于影像中心的位置XtKyo和焦距f，各参数初始化赋值如下：利用相机角坐标?，W，A·)的初始值确定旋转矩阵R:(4) 根据公式(2)计算各物方标识点(X1, Y130(1 = 1,2-44)在像坐标系下的理论像坐 标(xi，yi)(i = l，2···24)，然后根据旋转矩阵R、理论像坐标(Xi，yi) (i = I，2…24)和实际像坐标(xoi，y〇i) (i = I， 2···24)按照最小二乘法计算得到误差方程式系数矩阵仏（1 = 1，2 - 24)，并利用物方标识点 的实际像坐标(xoi，y〇i) (i = 1，2…24)减去理论像坐标值(Xi，yi) (i = 1，2··· 24)构成误差矩 阵 Li(i = l，2."24): Li=[lx，ly] =[XOi_Xi，y〇i_yi]; (5) 根据各物方标识点的误差方程式系数矩阵A1Q = I, 2···24)和误差矩阵L1Q = I, 2··· 24)分别构造两大矩阵A和L = A=U1 A2……L2……L24]T; 然后构造法方程，(6) 求解法方程（3)便可从公式右端得到各相机参数的修正量，判断三个角坐标(… w.，功对应的修正量(?, Aw，Aa:).是否都小于规定的阈值Δ，若满足条件则计算结束;若 不满足条件，通过解算出来的9个相机参数的修正量修改相机参数值，然后重复步骤(3)~ (6);若迭代解算次数超过20次仍然不满足条件，则说明硬件摄影测量系统有问题。8.如权利要求7所述的相机标定方法，其特征在于，步骤（1)中，在物方板的中心建立物 方坐标系XYZ，其中水平向右为X轴正方向，竖直向上为Y轴正方向，垂直于物方板为Z轴。
【专利摘要】本发明涉及一种激光摄影测量系统及相机标定方法，包括测量装置和标定装置，测量装置包括底座、支架及计算机；所述支架转动设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相机，所述计算机分别与相机、激光测距仪、线激光发射器电连接；标定装置包括三个物方板和均布设置的物方标识点。通过测量装置和标定装置，利物方板和物方标识点，通过相机拍照以及后续编写的相应程序，可以有效实现测量系统相机的精确标定，且该相机标定方法步骤简单，易于操作，为后续的激光摄影测量提供了可靠的基础。
【IPC分类】G01C25/00, G01C11/02
【公开号】CN105486289
【申请号】CN201610065915
【发明人】李学艺, 吕永刚, 曾庆良, 张鑫, 钟佩思, 王宁宁, 赵丹丹, 张庆雪
【申请人】山东科技大学
【公开日】2016年4月13日
【申请日】2016年1月31日
【技术领域】
[0001] 本发明设及一种激光摄影测量系统及相机标定方法，属于测量装置标定方法技术 领域。
【背景技术】
[0002] 随着计算机技术的发展，摄影测量与Ξ维重构技术在机器人视觉、工业制造、工程 建筑、文物保护和医学图像分析等领域得到广泛应用。相比于利用Ξ坐标测量仪等接触式 测量仪器对实物外表面进行测量，激光扫描等非接触式测量利用非接触的光电方法对曲面 的Ξ维形貌进行快速测量，具有非接触、速度快、精度高、应用广泛等优点，在工业工程等领 域有着一定的应用前景。在测量过程中，基于视觉传感器(如摄像机)获取的图像，建立二维 像点与Ξ维物点的对应关系，利用图像中的二维数据反算出对应的Ξ维物体的形貌和相关 位置等信息。
[0003] 要想获得高精度的测量结果，首先需要设计的测量系统具有较高的精度;其次建 立二维像点与Ξ维物点的对应关系，也即摄像机标定技术，又包括摄像机建模与摄像机标 定方法。相机标定主要是标定相机固有的内部几何参数和相机相对物方坐标之间的位置关 系等外部参数，其标定结果的好坏直接决定测量得到的Ξ维坐标点的精确程度，尤其是对 度量领域，因此需要研究相机的标定方法和研制相应的装置。
[0004] 目前相机标定技术根据解算相机内部几何参数和外部位置关系参数时所使用的 方法，将其分为线性方法、非线性方法和自标定方法。线性方法如直接线性变换方法(DLT变 换），根据环境点与相机图像对应点的关系，通过线性方法求得线性模型参数，进而建立线 性模型来描述相机的成像过程，该方法估计相机参数过程效率较高，但是没有考虑成像过 程中的非线性崎变问题，因此精度较低。非线性方法对提高摄像机标定精度提供了一条很 好的途径，非常适合于多参数的非线性崎变模型的标定，在能提供较好初值的情况下，能更 好地模拟和补偿各种像差，可W较快地收敛，得到高精度优化结果，但是它对标定时摄像机 模型参数的初始值要求较高，计算量大，计算速度慢。自标定方法在场景未知和相机运动参 数未知的情况下也可W标定相机的内外参数，在标定过程中仅依靠多幅图像对应点之间的 关系进行，其灵活性高、应用广泛，但是标定的过程主要是基于绝对二次曲线和绝对二次曲 面方法，此外还需要求解多元非线性方程，因此只适用于一些对精度要求不高的场合。
[0005] 造成目前所采用的标定方法存在上述问题的原因主要是其对应的测量系统较为 复杂。不同的测量系统需要采用不同的标定方法来标定，因此为了将标定方法简单化，尽可 能避免复杂标定方法造成的标定误差，首先需要使设计的测量系统简单，然后在针对设计 的测量系统设计标定方法。
[0006] 此外，针对大型的被测物，实际测量时往往是利用多台相机同时测量，而多台相机 之间没有相对位置关系，因此需要单独对每台相机进行标定，运样就增加了标定工作量，同 时还会增加标定误差，造成测量精度降低。所W在设计测量系统及标定方法时还应考虑同 时标定两台及W上的相机。
【发明内容】
[0007] 针对现有技术的不足，本发明提供一种激光摄影测量系统。
[0008] 本发明还提供上述一种激光摄影测量系统中相机的标定方法。
[0009] 本发明的技术方案如下：
[0010] -种激光摄影测量系统，包括测量装置，所述测量装置包括底座、支架及计算机； 所述支架设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与 线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相机，所述计算 机分别与相机、激光测距仪、线激光发射器电连接。
[0011] 优选的，所述支架的底部设置有水平分度盘。此设计的好处在于，当旋转支架时， 通过水平分度盘能够清楚直观地观察到支架旋转的角度。
[0012] 优选的，所述激光摄影测量系统还包括标定装置，所述标定装置包括物方标识点 和Ξ个物方板，Ξ个物方板的外形均为正方形，Ξ个物方板前后依次同轴平行设置且Ξ个 物方板的边长依次增大，在每个物方板上设置八个所述的物方标识点，八个物方标识点按 如下方式分布:在每条边的中点设置一个物方标识点，在相邻两条边的交点设置一个物方 标识点。
[001引优选的，所述立个物方板的尺寸分别为：100mm X 100mm X 10mm，200mm X 200mm X 10mm,300mmX300mmX10mm0
[0014] 优选的，所述物方标识点采用直径5mm的圆形薄片。
[0015] 优选的，所述标定装置还包括导轨，Ξ个物方板通过滑动支架设置在导轨上。此设 计的好处在于，将Ξ个物方板通过底部的滑动支架安装在导轨上，方便调节Ξ个物方板之 间的间距。
[0016] -种激光摄影测量系统的相机标定方法，包括W下步骤，
[0017] (1)将测量装置放置在标定装置的正前方，建立物方坐标系，并记录各物方标识点 在物方坐标系下的物方坐标(Xci，Yci，Zci)(i = l，2·-24):
[0018] (2)通过计算机控制相机对Ξ个物方板拍照，得到物方标识点的图像数据，通过计 算机内部执行程序识别出图像数据中物方标识点在相机的像坐标系下的实际像坐标(Χ01, y〇i)(i = l ,2---24)；
[0019] (3)对每台相机的9个相机参数赋初值，9个相机参数包括相机的线坐标(Xs，Ys， Zs)、角坐标(知W，K)、像主点相对于影像中屯、的位置X0、yo和焦距f，各参数初始化赋值如下：
[0022] 利用相机角坐标林，W，K)的初始值确定旋转矩阵R:
[0023]
[0024] (4)根据公式(2)计算各物方标识点(沿，打，21)。= 1，2-，24)在像坐标系下的理论 像坐标(xi，yi)(i = l，2···24)，
[0027] 然后根据旋转矩阵R、理论像坐标^1，7〇(1 = 1，2-，24)和实际像坐标^(^，7()〇(1 =1，2-'24)按照最小二乘法计算得到误差方程式系数矩阵心。=1，2-'24)，并利用物方标 识点的实际像坐标^〇1，7〇0。= 1，2，，，24)减去理论像坐标值^1，7〇(1 = 1,2，，，24)构成误 差矩阵k(i = l，2...24):
[002引 [lx, ]_y]T= [x0i-Xi,y0i_yi]T;
[0029] (5)根据各物方标识点的误差方程式系数矩阵Αια = 1，2···24)和误差矩阵^(i = 1，2···24)分别构造两大矩阵A和レA=[Al A2……A24]t，L=[Li L2……L24]t;
[0030] 然后构造法方程，X=(aTa)-1aTl=[Ms ΔΚ,么爲匈> Aw Δκ·Δ/Δ?0 Δ細 （雖
[0031] (6)求解法方程(3)便可从公式右端得到各相机参数的修正量，判断Ξ个角坐标（ 9，w，K)对应的修正量(Δ;0，Aw，Δκ)是否都小于规定的阔值Δ，若满足条件则计算结束;若 不满足条件，通过解算出来的9个相机参数的修正量修改相机参数值，然后重复步骤(3)~ (6);若迭代解算次数超过20次仍然不满足条件，则说明硬件摄影测量系统有问题。
[0032] 优选的，步骤（1)中，在物方板的中屯、建立物方坐标系ΧΥΖ，其中水平向右为X轴正 方向，竖直向上为Υ轴正方向，垂直于物方板为Ζ轴。
[0033] 本发明的有益效果在于：
[0034] (1)相比于目前已有的摄影测量系统的标定装置，本发明的标定装置由物方板和 物方标识点组成，结构简单，成本低、安装方便、操作简单、精度高。
[0035] (2)传统方法中将定位定姿数据转换到相机中屯、，如果存在两个或两个W上的相 机时，需要建立多个相机坐标系，给计算增加了难度，且不便于数据的归一化处理;本发明 中将坐标系建立在标定装置上，将两台相机 的位姿坐标都统一到物方坐标系下，便于计算， 减少了编程工作量，不再需要对两台相机解算的数据进行坐标变换处理。
[0036] (3)根据本发明相机标定方法编写的程序简单、易于操作、求解的精度较高，在图 像数据中识别物方标识点的像坐标时，快速、精确;迭代过程简单，整个标定过程不设及自 标定方法中由于二次曲线和二次曲面造成的多元非线性方程组。
【附图说明】
[0037] 图1为本发明中测量装置的结构示意图；
[0038] 图2为本发明中测量装置和标定装置的位置关系示意图；
[0039] 图3为本发明中测量装置的激光中屯、轴线垂直于物方板的关系图。
[0040] 其中：1、底座;2、水平分度盘;3、支架;4、相机;5、激光测距仪;6、线激光发射器;7、 相机;8、数据线;9、计算机；10、Ξ角架；11、测量装置；12、点激光；13、第一物方板；14、第二 物方板;15、第Ξ物方板;16、物方标识点；17、导轨;18、滑动支架。
【具体实施方式】
[0041] 下面通过实施例并结合附图对本发明做进一步说明，但不限于此。
[0042] 相机的标定，主要是标定相机在物方坐标系下的外方位元素和相机的内方位元 素。其中每台相机都有3个内方位元素和6个外方位元素，内方位元素包括像主点相对于影 像中屯、（即照片中屯、点）的位置xo、y〇W及镜头中屯、到影像面的焦距f (也称主距）。外方位元 素包括3个用于描述镜头中屯、相对于物方坐标系位置的线元素 Xs，Ys，Zs和3个用于描述影像 面在摄影瞬间空中姿态的角元素护、ω、k。激光测距仪点激光中屯、在物方坐标系中的坐标为 (0,0,Zcj)。
[0043] 实施例1:
[0044] 如图1和图2所示，一种激光摄影测量系统，该激光摄影测量系统包括测量装置和 标定装置，其中，测量装置包括底座1、水平分度盘2、支架3、两台相机4和7、线激光发射器6、 激光测距仪5、数据线8、计算机9。其中，相机4和相机7对称安装在激光测距仪5两侧，Ξ者都 安装在支架3上，计算机9通过数据线8与支架3内的控制设备连接。测量系统还包括Ξ脚架 10，底座1通过Ξ脚架10上的螺纹锁紧件与Ξ脚架10连接固定。
[0045] 相机4、相机7、线激光发射器6和激光测距仪5均由计算机9通过数据线8来传输控 制命令，并将获得的图像数据传递到计算机9中，方便数据的存储、计算和分析。
[0046] 计算机9装有相机4、相机7、线激光发射器6和激光测距仪5的驱动及控制程序，负 责控制发射激光、相机转动和拍摄照片等。
[0047] 计算机9还装有实现相机标定方法的执行程序，通过读取图像数据解算得到相机 位置和姿态坐标(位置坐标也即相机的线元素、姿态坐标也即相机的角元素，简称位姿坐 标)。
[0048] 如图2所示，标定装置是一个物方模型，该物方模型包括Ξ个用于布置物方标识点 16的物方板，Ξ个物方板包括第一物方板13、第二物方板14和第Ξ个物方板15,第一物方板 13和第二物方板14底部设置有滑动支架18,通过滑动支架1則尋第一物方板13和第二物方板 14放置导轨17上，第Ξ物方板15与导轨17的一端固定在一起(第Ξ物方板也可通过滑动支 架设置在导轨上）。第一物方板13、第二物方板14和第Ξ物方板15上均布置有8个物方标识 点16，均匀布置在每个物方板四条边上。
[0049] 其中，第一物方板13的尺寸(长宽厚)为lOOmmXlOOmmXlOmm，第二物方板14的尺 寸为200mmX200mmX10mm，第Ξ物方板15的尺寸为300mmX300mmX10mm，Ξ个物方板前后 依次平行放置，物方板的中屯、W及测量装置11的激光中屯、均在同一条线上，因此测量装置 11正对着物方模型获取物方标识点16的图像时，在一张照片上便能够同时显示所有物方标 识点的图像。
[0050] 物方标识点16是直径5mm的圆形薄片，为了与背景色有明显的区别，将物方标识点 16表面涂成黑色。
[0051] 第一物方板13和第二物方板14通过与其相连接的滑动支架18可在导轨17上滑动， 从而可W调节Ξ个物方板之间的距离。
[0052] 所述的物方模型上，在第Ξ物方板15的中屯、建立物方坐标系，水平向右为X轴正方 向，竖直向上为Y轴正方向，垂直于第Ξ物方板15为Z轴，方向遵循右手定则。
[0053]本实施例中的激光摄影测量系统结构简单、方便操作、精度高，是进行相机标定方 法的基础，同时高精度的激光摄影测量系统也是后续进行激光摄影测量而获得高精度测量 结果的保证。
[0化4] 实施例2:
[0055] 在实施例1的基础上，本实施例详细说明利用实施例1所述的激光摄影测量系统进 行相机标定方法的过程。
[0056] 本实施例提供一种激光摄影测量系统的相机标定方法，首先得到在物方坐标系下 各物方标识点16的坐标，然后调整测量装置11，使两台相机正对物方模型拍照，得到24个物 方标识点16的图像数据。通过计算机9上实现相机标定方法的执行程序读出各物方标识点 16的像坐标，利用物方标识点16的物方坐标与像坐标的对应关系，由执行程序计算得到两 台相机和激光测距仪点激光中屯、在物方坐标系中的位姿坐标及内方位元素参数。具体标定 方法如下：
[0057] 步骤1:确定各物方标识点16在物方坐标系下的坐标值(Xci，Yci，Zci) (i = 1，2… 24);
[0058] 步骤2:将测量装置11通过螺纹锁紧件固定在Ξ脚架10上，然后将计算机9通过数 据线8与测量装置相连，启动计算机9上控制相机与激光测距仪点激光中屯、的驱动及控制程 序，通过控制程序调整测量系统位姿，使激光测距仪轴线水平；
[0059] 步骤3:调整物方模型的位置，使其正对着测量装置11，第一物方板13在最前面；
[0060] 步骤4:如图2、图3所示，调整测量装置11的位姿，使点激光12与物方模型的第一物 方板13垂直。具体步骤如下：
[0061] (1)通过计算机9上的控制程序，打开测量装置11的激光测距仪5，激光测距仪5发 射的点激光12，通过肉眼观察使点激光12近似垂直于物方模型的第一物方板13，点激光12 打在第一物方板13的化位置，如图3所示；
[00创 （2似0102线为基准位置，通过旋转支架3，使0102绕点激光中屯、01点逆时针转过目1 到OiA的位置，通过激光测距得到化A的长度La;
[0063] (3) W化化线为基准位置，通过旋转支架3,使化化绕测量装置中屯、化点顺时针转过 θι至化iB的位置，通过激光测距得到化B的长度Lb;
[0064] (4)比较La与Lb的大小，如果La〉Lb，转动水平支架3使〇2右移向B靠近;如果La<Lb，转 动水平支架使化左移向A靠近；
[0065] (5)如果La = Lb，则点激光12在水平方向垂直于第一物方板13;否则重复步骤(2)~ (4)，最终使点激光12与物方模型的第一物方板13在水平方向上垂直；
[0066] (6)保证点激光12在竖直方向上与第一物方板13垂直的方法同上述步骤，只需要 将方向改为在竖直方向上微调点激光12的位置。
[0067] 步骤5:通过计算机9上的控制程序，打开相机4和相机7，对物方模型拍照，通过程 序分别获取各个物方标识点16在像坐标系下的实际像坐标(xoi，yoi) (i = 1，2· · · 24);
[0068] 步骤6:给每台相机的9个参数赋初值，9个参数包括相机的线坐标(Xs，Ys，Zs)、角坐 标(r，w，K)、像主点相对于影像中屯、（即照片中屯、点）的位置x〇、y〇和焦距f，各参数初始化赋 值如下：
[0069] 口二w='A=.v(尸.!',)=0，30</.:? 120;
[0070] Xs = ( Σ Xci) /24，Ys = ( Σ Yci) /24，Zs = ( Σ Zci) /24;
[0071] 步骤7:根据公式(1)，利用相机角坐标元素值界、*、1〇确定旋转矩阵R中各元素的值：
[0072]
[0073] 步骤8:利用基于空间后方交会原理编写的程序计算得到理想情况下物方标识点 在相机的像坐标系下的理论像坐标，根据公式(2)计算各物方标识点(Χι，Υι，Ζι)在像坐标系 下的理论像坐标(xi，yi)(i = l，2···24)，
[0076] 步骤9:然后根据旋转矩阵R、各物方标识点的理论像坐标(Χι，γι)α = 1，2···24^Ρ 实际像坐标(x()l，y()l)(i = l，2…24)按照最小二乘法计算误差方程式系数矩阵Al(i = l，2… 24);
[0077] 步骤10:利用物方标识点像坐标的实际值(xoi，y日i)减去理论像坐标值(xi，y i)构成 误差矩阵以。=1，2-'24):
[007引 [lx, ；Ly]T= [x0i-Xi，y0i-yi]T;
[0079] 步骤11:根据各物方标识点系数矩阵Αια = 1，2···24)和误差矩阵^α = 1，2···24) 分别构造成两个大矩阵4和1^:4=[41 Α2 ...... A24]T，L=[Li L2 ...... L24]T;
[0080] 步骤 12:构造法方程，X=(aTa)-iATl=[A^ ΔΚ Δ& Aw ΔκΔ/細ο Δ础，求解法方 程便可从公式右端得到相机各参数的修正量；
[0081 ] 步骤13:判断Ξ个角坐标作，W，Κ)对应的修正量(Δ口，Δ W，Δ Κ)是否都小于规定的 阔值A，如果满足条件则计算结束;如果不满足条件，通过解算出来的9个相机参数的修正 量修改相机参数值，然后重复步骤7~12;如果迭代解算次数超过20次仍然不满足条件，说 明硬件摄影测量系统有问题。
[0082] 步骤14:两台相机都按照上述步骤进行解算，不同之处在于物方标识点的像坐标 分别取自两台相机各自拍到的图像数据，最终解算得到两台相机的9个参数。由于激光测距 仪点激光中屯、位于两台相机的对称中屯、，因此其线坐标为两台相机线坐标的平均值。
[0083] 实际在拍摄照片进行相机参数解算时，为了避免随机误差，每个相机拍摄3~5张 照片，并解算每张照片的数据，通过分析数据，剔除偏差较大的数据，对剩下的几组数据取 平均值得到相机的9个参数。
[0084] 本发明在实施过程中借鉴了现有技术，限于篇幅，未对现有技术部分进行详细描 述;凡是本发明未提及的技术部分，均可W采用现有技术实现。
[0085] W上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技 术人员来说，本发明可W有各种更改和变化。凡在本发明的精神和原则之内，所做的任何修 改，等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种激光摄影测量系统，包括测量装置，其特征在于，所述测量装置包括底座、支架 及计算机;所述支架设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激 光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相 机，所述计算机分别与相机、激光测距仪、线激光发射器电连接。2. 如权利要求1所述的激光摄影测量系统，其特征在于，所述支架的底部设置有水平分 度盘。3. 如权利要求1所述的激光摄影测量系统，其特征在于，所述激光摄影测量系统还包括 标定装置，所述标定装置包括物方标识点和三个物方板，三个物方板的外形均为正方形，三 个物方板前后依次同轴平行设置且三个物方板的边长依次增大，在每个物方板上设置八个 所述的物方标识点，八个物方标识点按如下方式分布:在每条边的中点设置一个物方标识 点，在相邻两条边的交点设置一个物方标识点。4. 如权利要求3所述的激光摄影测量系统，其特征在于，所述三个物方板的尺寸分别 为：100mm X 100mm X 10mm，200mm X 200mm X 10mm，300mm X 300mm X 1 Omnin5. 如权利要求3所述的激光摄影测量系统，其特征在于，所述物方标识点采用直径5mm 的圆形薄片。6. 如权利要求3所述的激光摄影测量系统，其特征在于，所述标定装置还包括导轨，三 个物方板通过滑动支架设置在导轨上。7. -种如权利要求1-6任一项所述的激光摄影测量系统的相机标定方法，包括以下步 骤， (1) 将测量装置放置在标定装置的正前方，建立物方坐标系，并记录各物方标识点在物 方坐标系下的物方坐标(Xc^Ycd,Zed)(i = l，2···24); (2) 通过计算机控制相机对三个物方板拍照，得到物方标识点的图像数据，通过计算机 内部执行程序识别出图像数据中物方标识点在相机的像坐标系下的实际像坐标(XQ 1J01) (i = l,2···24); (3) 对每台相机的9个相机参数赋初值，9个相机参数包括相机的线坐标(XS，YS，ZS)、角 坐标0，W，4、:像主点相对于影像中心的位置XtKyo和焦距f，各参数初始化赋值如下：利用相机角坐标?，W，A·)的初始值确定旋转矩阵R:(4) 根据公式(2)计算各物方标识点(X1, Y130(1 = 1,2-44)在像坐标系下的理论像坐 标(xi，yi)(i = l，2···24)，然后根据旋转矩阵R、理论像坐标(Xi，yi) (i = I，2…24)和实际像坐标(xoi，y〇i) (i = I， 2···24)按照最小二乘法计算得到误差方程式系数矩阵仏（1 = 1，2 - 24)，并利用物方标识点 的实际像坐标(xoi，y〇i) (i = 1，2…24)减去理论像坐标值(Xi，yi) (i = 1，2··· 24)构成误差矩 阵 Li(i = l，2."24): Li=[lx，ly] =[XOi_Xi，y〇i_yi]; (5) 根据各物方标识点的误差方程式系数矩阵A1Q = I, 2···24)和误差矩阵L1Q = I, 2··· 24)分别构造两大矩阵A和L = A=U1 A2……L2……L24]T; 然后构造法方程，(6) 求解法方程（3)便可从公式右端得到各相机参数的修正量，判断三个角坐标(… w.，功对应的修正量(?, Aw，Aa:).是否都小于规定的阈值Δ，若满足条件则计算结束;若 不满足条件，通过解算出来的9个相机参数的修正量修改相机参数值，然后重复步骤(3)~ (6);若迭代解算次数超过20次仍然不满足条件，则说明硬件摄影测量系统有问题。8.如权利要求7所述的相机标定方法，其特征在于，步骤（1)中，在物方板的中心建立物 方坐标系XYZ，其中水平向右为X轴正方向，竖直向上为Y轴正方向，垂直于物方板为Z轴。
【专利摘要】本发明涉及一种激光摄影测量系统及相机标定方法，包括测量装置和标定装置，测量装置包括底座、支架及计算机；所述支架转动设置在底座上，在支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧对称设置两台所述的相机，所述计算机分别与相机、激光测距仪、线激光发射器电连接；标定装置包括三个物方板和均布设置的物方标识点。通过测量装置和标定装置，利物方板和物方标识点，通过相机拍照以及后续编写的相应程序，可以有效实现测量系统相机的精确标定，且该相机标定方法步骤简单，易于操作，为后续的激光摄影测量提供了可靠的基础。
【IPC分类】G01C25/00, G01C11/02
【公开号】CN105486289
【申请号】CN201610065915
【发明人】李学艺, 吕永刚, 曾庆良, 张鑫, 钟佩思, 王宁宁, 赵丹丹, 张庆雪
【申请人】山东科技大学
【公开日】2016年4月13日
【申请日】2016年1月31日
最新回复(0)