CN110057295A - 一种免像控的单目视觉平面距离测量方法 - Google Patents
一种免像控的单目视觉平面距离测量方法 Download PDFInfo
- Publication number
- CN110057295A CN110057295A CN201910275119.2A CN201910275119A CN110057295A CN 110057295 A CN110057295 A CN 110057295A CN 201910275119 A CN201910275119 A CN 201910275119A CN 110057295 A CN110057295 A CN 110057295A
- Authority
- CN
- China
- Prior art keywords
- point
- camera
- plane
- angle
- distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/022—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by means of tv-camera scanning
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种免像控的单目视觉平面距离测量方法，包括：测量系统集成阶段：将相机和激光测距仪固连构成直接定向视觉测量系统；测量系统检校阶段：对测量系统进行集成检校，确定相机内参、畸变像差以及相机与激光测距仪之间的偏心角，修正相机的非线性畸变、实测高程和姿态；平面距离测量阶段：基于相机在倾斜视角下像点、光心和物点共线的中心投影成像原理建立像平面和物平面之间的“像点‑距离”的变换关系，最终实现免像控的单目视觉平面距离测量。本发明具有较低的硬件成本和操作复杂度。本发明无需在待测平面布设像控点，可在数分钟内完成测点布设，大大降低了野外工作量；也不依赖于场景中的先验知识，具有普适性。
Description
技术领域
本发明涉及平面距离测量领域，特别是涉及一种免像控的单目视觉平面距离测量方法。
背景技术
视觉测量是随着计算机视觉技术的发展而形成的一门非接触式测量技术，相比传统测量手段具有量测范围大、获取测量信息多、受工作环境限制小、不干扰测量对象等优点。视觉测量系统根据所用摄像机数目的多少，可以分为单目视觉测量、双目视觉测量、多目视觉测量和结构光视觉测量等。其中，利用单台相机拍摄的单幅图像进行物体几何尺寸或空间位置测量的技术在机器视觉领域被称为单目视觉测量。由于具有结构简单、易于标定、操作便捷等优点，被广泛应用于工业检测及测绘工程中的平面测量。然而，单幅图像由于成像过程中深度信息的缺失，造成三维世界坐标到二维图像坐标的映射关系是不可逆的，通常需要附加共面约束条件来建立单应矩阵描述的像方平面和物方平面间可逆的坐标变换关系。主要采用的方法是基于点对应求解单应矩阵的直接线性变换法，即：使用待测平面中4个以上坐标己知的像控点解算单应矩阵，再利用单应性矩阵解算待测像点在物平面中的坐标，进而计算两点之间的距离；其测量精度对像控点的精度、数量及分布条件等因素较为敏感，并且在水面等复杂施测环境中像控点的布设较为困难，往往限制了该方法的使用。其他方法利用视场中已知的距离或规则的几何结构作为先验知识解算单应矩阵，例如：Wang等提出了一种利用基于直线对应的单应矩阵解算方法，由于直线可通过多点拟合得到，因此测量精度高高于点对应的方法；丁伟利等提出一种利用矩形结构获取测量信息的方法，仅需利用场景中易于测量的矩形一边长度作为输入，就能够对场景中不同位置进行测量；王美珍等提出了基于交比射影变换不变量的单幅图像平面几何信息提取方法，利用场景中几何结构计算平面中的灭点灭线，并利用待测平面中事先测得距离的线段构建交比来计算待求线段的物理距离，避免了布设像控点的繁琐；然而，这些方法需要利用图像中的几何信息，只适用于特定场景，缺乏普适性。
发明内容
发明目的：本发明的目的是提供一种免像控的单目视觉平面距离测量方法，具有普适性。
技术方案：本发明所述的免像控的单目视觉平面距离测量方法，包括：
测量系统集成阶段：将相机和激光测距仪固连构成直接定向视觉测量系统；
测量系统检校阶段：对测量系统进行集成检校，确定相机内参、畸变像差以及相机与激光测距仪之间的偏心角，修正相机的非线性畸变、实测高程和姿态；
平面距离测量阶段：基于相机在倾斜视角下像点、光心和物点共线的中心投影成像原理建立像平面和物平面之间的“像点-距离”的变换关系，最终实现免像控的单目视觉平面距离测量。
进一步，所述将相机和激光测距仪固连构成直接定向视觉测量系统的具体过程为：采用一根连接杆将相机和激光测距仪固连构成直接定向视觉测量系统，连接杆上设有三个共线分布的螺孔，位于中间的螺孔用于连接三脚架，位于两侧的螺孔分别连接相机和激光测距仪，将相机与激光测距仪的三轴调至平行。
进一步，位于两侧的螺孔中的一个螺孔通过双头连接件连接相机，位于两侧的螺孔中的另一个螺孔通过双头连接件连接激光测距仪。
进一步，所述对测量系统进行集成检校的过程采用一块网格尺寸为(DX,DY)、网格数为m×n的平面棋盘格标定板作为参照物，其中，DX为网格的长度，DY为网格的宽度，m为横向网格数，n为纵向网格数；对测量系统进行集成检校的过程具体包括以下步骤：
S11：拍摄9幅不同姿态下的棋盘格标定板图像，采用张正友法检校相机，得到相机内参(fx,fy,Cx,Cy)和畸变像差(k1,k2,p1,p2)，其中Cx表示像主点的横坐标，Cy表示像主点的纵坐标，fx表示相机在像平面x轴方向上的等效焦距，fy表示相机在像平面y轴方向上的等效焦距，k1表示第一径向畸变系数，k2表示第二径向畸变系数，p1表示第一切向畸变系数，p2表示第二切向畸变系数；
S12：将平面棋盘格标定板水平放置，采用三脚架固定测量系统并拍摄一幅完整平面棋盘格标定板的畸变图像I(ud,vd)；(ud,vd)为畸变图像坐标；
S13：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、横滚角的测量值以及激光测距仪到平面棋盘格标定板的斜距DL；
S14：为了校正图像的非线性畸变，首先根据相机内参描述的尺度缩放关系将无畸变图像坐标(u,v)转换为无畸变相机坐标(x,y)：
然后根据像差模型建立无畸变相机坐标(x,y)到畸变相机坐标(x′,y′)的映射关系：
接下来根据式(1)的反变换关系将(x′,y′)变换到畸变图像坐标(ud,vd)：
最后根据式(1)～(3)建立的(u,v)到(ud,vd)的映射关系，采用最邻近插值法完成畸变图像I(ud,vd)到畸变校正图像I(u,v)的灰度映射；
S15：在畸变校正图像I(u,v)中检测平面棋盘格标定板上所有(m-1)×(n-1)个内角点的亚像素坐标，其中，第i个内角点的亚像素坐标为(ui,vi)，i＝1,2,…,(m-1)×(n-1)；
S16：令相机俯仰角中间值ωk＝ωL-2+0.1×k、相机横滚角中间值
S17：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O′为光心在平面棋盘格标定板上的投影点，计算光心O到平面棋盘格标定板的高程H：
H＝DL·sinωL+d·cosωk (4)
式(4)中，d为测量基准点到相机平面的投影距离；
S18：对于平面棋盘格标定板上的任意点P，其在畸变校正后的原始像平面上对应像点p的坐标用(u1,v1)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角中间值进行校正，得到的像平面上像点p′的坐标
其中，w1为图像的宽度，h1为图像的高度；
S19：将相机透镜平面的光心O在平面棋盘格标定板上的垂足点C作为参考点，计算物点P到点C的水平距离M和垂直距离N，过程如下：定义p″为p′在的像平面主纵线上的投影点，P″为物点P在平面棋盘格标定板主纵线O′C上的投影点，根据三角关系得到O到P″的距离L为：
L＝H/sinβ (6)
其中，β表示射线OP″和主纵线O′C间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(8)得到：
f＝s·(fx+fy)/2 (8)
根据三角形的相似性，有：
P″P/P″O＝p′p″/p″O (9)
代入对应参数，得：
同理，有：
S20：对于像平面中任意两点p1和p2，首先根据步骤S19分别求得其物点相对于垂足点C的距离(M1,N1)和(M2,N2)，然后代入下式：
其中，M1为点p1相对于垂足点C的水平距离，N1为点p1相对于垂足点C的垂直距离，M2为点p2相对于垂足点C的水平距离，N2为点p2相对于垂足点C的垂直距离，D(p1,p2)为点p1和点p2之间的距离；
S21：根据式(12)分别求得平面棋盘格标定板中相邻水平内角点之间的距离D(pa,b,pa+1,b)和相邻垂直内角点之间的距离D(pa,b,pa,b+1)作为平面棋盘格标定板网格尺寸的测量值，pa,b为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b的内角点，pa+1,b为平面棋盘格标定板的内角点中横向编号为a+1、纵向编号为b的内角点，pa,b+1为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b+1的内角点，a＝1,2,…,(m-1)，b＝1,2,…,(n-1)；根据下式统计平面棋盘格标定板上所有网格距离测量值的总体均方根误差：
S22：将RMSE取得最小值所对应的ωk作为相机俯仰角的估计值ωE，将RMSE取得最小值所对应的作为相机横滚角的估计值并根据式(14)进行计算：
式(14)中，为相机横滚角估计值，ωE为相机俯仰角估计值，为相机横滚偏心角，Δω为相机俯仰偏心角。
进一步，所述平面距离测量阶段具体包括以下步骤：
S31：采用三脚架固定测量系统，通过调节合适的方位角使相机视场覆盖完整的测量区域且光轴指向测量区域中心，调节水平使横滚角近似为0°；
S32：拍摄一幅图像并对图像进行畸变校正，消除径向和切向畸变并使校正后的像主点位于图像中心；
S33：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、激光测距仪横滚角的测量值和以及激光测距仪到物平面的斜距DL1；
S34：根据测量系统检校阶段得到的相机横滚偏心角和相机俯仰偏心角Δω，通过式(15)得到相机的横滚角和俯仰角ωC；
S35：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O″为光心在物平面上的投影点，计算光心O到物平面的高程H1：
H1＝DL1·sinωL+d·cosωC (16)
S36：对于物平面上的任意点Q，其在畸变校正后的原始像平面上对应像点q的坐标用(u2,v2)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角进行校正，得到的像平面上像点q′的坐标
其中，w2为图像的宽度，h2为图像的高度；
S37：将相机透镜平面的光心O在物平面上的垂足点C′作为参考点，计算物点Q到点C′的水平距离M′和垂直距离N′，过程如下：定义q″为q′在的像平面主纵线上的投影点，Q″为物点Q在物平面主纵线O″C′上的投影点，根据三角关系得到O到Q″的距离L1为：
L1＝H1/sinβ1 (18)
其中，β1表示射线OQ″和主纵线O″C′间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(20)得到：
f＝s·(fx+fy)/2 (20)
根据三角形的相似性，有：
Q″Q/Q″O＝q′q″/q″O (21)
代入对应参数，得：
同理，有：
S38：对于像平面中任意两点q1和q2，首先根据步骤S37分别求得其物点相对于垂足点C′的距离(M1′,N1′)和(M2′,N2′)，然后代入下式：
其中，M1′为点q1相对于垂足点C′的水平距离，N1′为点q1相对于垂足点C′的垂直距离，M2′为点q2相对于垂足点C′的水平距离，N2′为点q2相对于垂足点C′的垂直距离，D(q1,q2)为点q1和点q2之间的距离。
有益效果：本发明公开了一种免像控的单目视觉平面距离测量方法，与现有技术相比，具有如下的有益效果：
1)本发明仅需一台相机和一台具有倾角测量功能的激光测距仪，而无需全站仪等专业测绘仪器，具有较低的硬件成本和操作复杂度；
2)本发明无需在待测平面布设像控点，可在数分钟内完成测点布设，大大降低了野外工作量；也不依赖于场景中的先验知识，具有普适性。
附图说明
图1为本发明具体实施方式中直接定向视觉测量系统的结构示意图；
图2为本发明具体实施方式中测量系统检校阶段的流程图；
图3为本发明具体实施方式中相机在倾斜视角下的中心投影成像模型示意图。
具体实施方式
本具体实施方式公开了一种免像控的单目视觉平面距离测量方法，包括：
测量系统集成阶段：将相机和激光测距仪固连构成直接定向视觉测量系统；
测量系统检校阶段：对测量系统进行集成检校，确定相机内参、畸变像差以及相机与激光测距仪之间的偏心角，修正相机的非线性畸变、实测高程和姿态；
平面距离测量阶段：基于相机在倾斜视角下像点、光心和物点共线的中心投影成像原理建立像平面和物平面之间的“像点-距离”的变换关系，最终实现免像控的单目视觉平面距离测量。
将相机和激光测距仪固连构成直接定向视觉测量系统的具体过程为：采用一根连接杆将相机和激光测距仪固连构成直接定向视觉测量系统，如图1所示，连接杆1上设有三个共线分布的1/4″螺孔，位于中间的螺孔用于连接三脚架，位于两侧的螺孔分别通过双头连接件4连接相机2和激光测距仪3，将相机2与激光测距仪3的三轴调至平行。激光测距仪3的测量距离大于100m，测距精度优于3mm；内置两轴倾角测量模块3-3，俯仰角的测量范围为-90°～+90°，横滚角的测量范围为-45°～+45°，测角精度优于0.1°。此外，图1中还示出了相机底座2-1、相机光轴2-2、激光测距仪底座3-1、激光测距仪光轴3-2。
对测量系统进行集成检校的过程如图2所示，采用一块网格尺寸为(DX,DY)、网格数为m×n的平面棋盘格标定板作为参照物，其中，DX为网格的长度，DY为网格的宽度，m为横向网格数，n为纵向网格数；对测量系统进行集成检校的过程具体包括以下步骤：
S11：拍摄9幅不同姿态下的棋盘格标定板图像，采用张正友法检校相机，得到相机内参(fx,fy,Cx,Cy)和畸变像差(k1,k2,p1,p2)，其中Cx表示像主点的横坐标，Cy表示像主点的纵坐标，fx表示相机在像平面x轴方向上的等效焦距，fy表示相机在像平面y轴方向上的等效焦距，k1表示第一径向畸变系数，k2表示第二径向畸变系数，p1表示第一切向畸变系数，p2表示第二切向畸变系数；
S12：将平面棋盘格标定板水平放置，采用三脚架固定测量系统并拍摄一幅完整平面棋盘格标定板的畸变图像I(ud,vd)；(ud,vd)为畸变图像坐标；
S13：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、横滚角的测量值以及激光测距仪到平面棋盘格标定板的斜距DL；
S14：为了校正图像的非线性畸变，首先根据相机内参描述的尺度缩放关系将无畸变图像坐标(u,v)转换为无畸变相机坐标(x,y)：
然后根据像差模型建立无畸变相机坐标(x,y)到畸变相机坐标(x′,y′)的映射关系：
接下来根据式(1)的反变换关系将(x′,y′)变换到畸变图像坐标(ud,vd)：
最后根据式(1)～(3)建立的(u,v)到(ud,vd)的映射关系，采用最邻近插值法完成畸变图像I(ud,vd)到畸变校正图像I(u,v)的灰度映射；
S15：在畸变校正图像I(u,v)中检测平面棋盘格标定板上所有(m-1)×(n-1)个内角点的亚像素坐标，其中，第i个内角点的亚像素坐标为(ui,vi)，i＝1,2,…,(m-1)×(n-1)；
S16：令相机俯仰角中间值ωk＝ωL-2+0.1×k、相机横滚角中间值
S17：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O′为光心在平面棋盘格标定板上的投影点，计算光心O到平面棋盘格标定板的高程H：
H＝DL·sinωL+d·cosωk (4)
式(4)中，d为测量基准点到相机平面的投影距离；
S18：对于平面棋盘格标定板上的任意点P，其在畸变校正后的原始像平面上对应像点p的坐标用(u1,v1)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角中间值进行校正，得到的像平面上像点p′的坐标
其中，w1为图像的宽度，h1为图像的高度；
S19：将相机透镜平面的光心O在平面棋盘格标定板上的垂足点C作为参考点，计算物点P到点C的水平距离M和垂直距离N，过程如下：定义p″为p′在的像平面主纵线上的投影点，P″为物点P在平面棋盘格标定板主纵线O′C上的投影点，根据三角关系得到O到P″的距离L为：
L＝H/sinβ (6)
其中，β表示射线OP″和主纵线O′C间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(8)得到：
f＝s·(fx+fy)/2 (8)
根据三角形的相似性，有：
P″P/P″O＝p′p″/p″O (9)
代入对应参数，得：
同理，有：
S20：对于像平面中任意两点p1和p2，首先根据步骤S19分别求得其物点相对于垂足点C的距离(M1,N1)和(M2,N2)，然后代入下式：
其中，M1为点p1相对于垂足点C的水平距离，N1为点p1相对于垂足点C的垂直距离，M2为点p2相对于垂足点C的水平距离，N2为点p2相对于垂足点C的垂直距离，D(p1,p2)为点p1和点p2之间的距离；
S21：根据式(12)分别求得平面棋盘格标定板中相邻水平内角点之间的距离D(pa,b,pa+1,b)和相邻垂直内角点之间的距离D(pa,b,pa,b+1)作为平面棋盘格标定板网格尺寸的测量值，pa,b为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b的内角点，pa+1,b为平面棋盘格标定板的内角点中横向编号为a+1、纵向编号为b的内角点，pa,b+1为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b+1的内角点，a＝1,2,…,(m-1)，b＝1,2,…,(n-1)；根据下式统计平面棋盘格标定板上所有网格距离测量值的总体均方根误差：
S22：将RMSE取得最小值所对应的ωk作为相机俯仰角的估计值ωE，将RMSE取得最小值所对应的作为相机横滚角的估计值并根据式(14)进行计算：
式(14)中，为相机横滚角估计值，ωE为相机俯仰角估计值，为相机横滚偏心角，Δω为相机俯仰偏心角。
相机在倾斜视角下的中心投影成像模型示意图如图3所示。平面距离测量阶段具体包括以下步骤：
S31：采用三脚架固定测量系统，通过调节合适的方位角使相机视场覆盖完整的测量区域且光轴指向测量区域中心，调节水平使横滚角近似为0°；
S32：拍摄一幅图像并对图像进行畸变校正，消除径向和切向畸变并使校正后的像主点位于图像中心；
S33：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、激光测距仪横滚角的测量值和以及激光测距仪到物平面的斜距DL1；
S34：根据测量系统检校阶段得到的相机横滚偏心角和相机俯仰偏心角Δω，通过式(15)得到相机的横滚角和俯仰角ωC；
S35：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O″为光心在物平面上的投影点，计算光心O到物平面的高程H1：
H1＝DL1·sinωL+d·cosωC (16)
S36：对于物平面上的任意点Q，其在畸变校正后的原始像平面上对应像点q的坐标用(u2,v2)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角进行校正，得到的像平面上像点q′的坐标
其中，w2为图像的宽度，h2为图像的高度；
S37：将相机透镜平面的光心O在物平面上的垂足点C′作为参考点，计算物点Q到点C′的水平距离M′和垂直距离N′，过程如下：定义q″为q′在的像平面主纵线上的投影点，Q″为物点Q在物平面主纵线O″C′上的投影点，根据三角关系得到O到Q″的距离L1为：
L1＝H1/sinβ1 (18)
其中，β1表示射线OQ″和主纵线O″C′间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(20)得到：
f＝s·(fx+fy)/2 (20)
根据三角形的相似性，有：
Q″Q/Q″O＝q′q″/q″O (21)
代入对应参数，得：
同理，有：
S38：对于像平面中任意两点q1和q2，首先根据步骤S37分别求得其物点相对于垂足点C′的距离(M1′,N1′)和(M2′,N2′)，然后代入下式：
其中，M1′为点q1相对于垂足点C′的水平距离，N1′为点q1相对于垂足点C′的垂直距离，M2′为点q2相对于垂足点C′的水平距离，N2′为点q2相对于垂足点C′的垂直距离，D(q1,q2)为点q1和点q2之间的距离。
Claims (5)
1.一种免像控的单目视觉平面距离测量方法，其特征在于：包括：
测量系统集成阶段：将相机和激光测距仪固连构成直接定向视觉测量系统；
测量系统检校阶段：对测量系统进行集成检校，确定相机内参、畸变像差以及相机与激光测距仪之间的偏心角，修正相机的非线性畸变、实测高程和姿态；
平面距离测量阶段：基于相机在倾斜视角下像点、光心和物点共线的中心投影成像原理建立像平面和物平面之间的“像点-距离”的变换关系，最终实现免像控的单目视觉平面距离测量。
2.根据权利要求1所述的免像控的单目视觉平面距离测量方法，其特征在于：所述将相机和激光测距仪固连构成直接定向视觉测量系统的具体过程为：采用一根连接杆将相机和激光测距仪固连构成直接定向视觉测量系统，连接杆上设有三个共线分布的螺孔，位于中间的螺孔用于连接三脚架，位于两侧的螺孔分别连接相机和激光测距仪，将相机与激光测距仪的三轴调至平行。
3.根据权利要求2所述的免像控的单目视觉平面距离测量方法，其特征在于：位于两侧的螺孔中的一个螺孔通过双头连接件连接相机，位于两侧的螺孔中的另一个螺孔通过双头连接件连接激光测距仪。
4.根据权利要求1所述的免像控的单目视觉平面距离测量方法，其特征在于：所述对测量系统进行集成检校的过程采用一块网格尺寸为(DX,DY)、网格数为m×n的平面棋盘格标定板作为参照物，其中，DX为网格的长度，DY为网格的宽度，m为横向网格数，n为纵向网格数；对测量系统进行集成检校的过程具体包括以下步骤：
S11：拍摄9幅不同姿态下的棋盘格标定板图像，采用张正友法检校相机，得到相机内参(fx,fy,Cx,Cy)和畸变像差(k1,k2,p1,p2)，其中Cx表示像主点的横坐标，Cy表示像主点的纵坐标，fx表示相机在像平面x轴方向上的等效焦距，fy表示相机在像平面y轴方向上的等效焦距，k1表示第一径向畸变系数，k2表示第二径向畸变系数，p1表示第一切向畸变系数，p2表示第二切向畸变系数；
S12：将平面棋盘格标定板水平放置，采用三脚架固定测量系统并拍摄一幅完整平面棋盘格标定板的畸变图像I(ud,vd)；(ud,vd)为畸变图像坐标；
S13：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、横滚角的测量值以及激光测距仪到平面棋盘格标定板的斜距DL；
S14：为了校正图像的非线性畸变，首先根据相机内参描述的尺度缩放关系将无畸变图像坐标(u,v)转换为无畸变相机坐标(x,y)：
然后根据像差模型建立无畸变相机坐标(x,y)到畸变相机坐标(x′,y′)的映射关系：
接下来根据式(1)的反变换关系将(x′,y′)变换到畸变图像坐标(ud,vd)：
最后根据式(1)～(3)建立的(u,v)到(ud,vd)的映射关系，采用最邻近插值法完成畸变图像I(ud,vd)到畸变校正图像I(u,v)的灰度映射；
S15：在畸变校正图像I(u,v)中检测平面棋盘格标定板上所有(m-1)×(n-1)个内角点的亚像素坐标，其中，第i个内角点的亚像素坐标为(ui,vi)，i＝1,2,…,(m-1)×(n-1)；
S16：令相机俯仰角中间值ωk＝ωL-2+0.1×k、相机横滚角中间值
S17：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O′为光心在平面棋盘格标定板上的投影点，计算光心O到平面棋盘格标定板的高程H：
H＝DL·sinωL+d·cosωk (4)
式(4)中，d为测量基准点到相机平面的投影距离；
S18：对于平面棋盘格标定板上的任意点P，其在畸变校正后的原始像平面上对应像点p的坐标用(u1,v1)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角中间值进行校正，得到的像平面上像点p′的坐标
其中，w1为图像的宽度，h1为图像的高度；
S19：将相机透镜平面的光心O在平面棋盘格标定板上的垂足点C作为参考点，计算物点P到点C的水平距离M和垂直距离N，过程如下：定义p″为p′在的像平面主纵线上的投影点，P″为物点P在平面棋盘格标定板主纵线O′C上的投影点，根据三角关系得到O到P″的距离L为：
L＝H/sinβ (6)
其中，β表示射线OP″和主纵线O′C间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(8)得到：
f＝s·(fx+fy)/2 (8)
根据三角形的相似性，有：
P″P/P″O＝p′p″/p″O (9)
代入对应参数，得：
同理，有：
S20：对于像平面中任意两点p1和p2，首先根据步骤S19分别求得其物点相对于垂足点C的距离(M1,N1)和(M2,N2)，然后代入下式：
其中，M1为点p1相对于垂足点C的水平距离，N1为点p1相对于垂足点C的垂直距离，M2为点p2相对于垂足点C的水平距离，N2为点p2相对于垂足点C的垂直距离，D(p1,p2)为点p1和点p2之间的距离；
S21：根据式(12)分别求得平面棋盘格标定板中相邻水平内角点之间的距离D(pa,b,pa+1,b)和相邻垂直内角点之间的距离D(pa,b,pa,b+1)作为平面棋盘格标定板网格尺寸的测量值，pa,b为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b的内角点，pa+1,b为平面棋盘格标定板的内角点中横向编号为a+1、纵向编号为b的内角点，pa,b+1为平面棋盘格标定板的内角点中横向编号为a、纵向编号为b+1的内角点，a＝1,2,…,(m-1)，b＝1,2,…,(n-1)；根据下式统计平面棋盘格标定板上所有网格距离测量值的总体均方根误差：
S22：将RMSE取得最小值所对应的ωk作为相机俯仰角的估计值ωE，将RMSE取得最小值所对应的作为相机横滚角的估计值并根据式(14)进行计算：
式(14)中，为相机横滚角估计值，ωE为相机俯仰角估计值，为相机横滚偏心角，Δω为相机俯仰偏心角。
5.根据权利要求4所述的免像控的单目视觉平面距离测量方法，其特征在于：所述平面距离测量阶段具体包括以下步骤：
S31：采用三脚架固定测量系统，通过调节合适的方位角使相机视场覆盖完整的测量区域且光轴指向测量区域中心，调节水平使横滚角近似为0°；
S32：拍摄一幅图像并对图像进行畸变校正，消除径向和切向畸变并使校正后的像主点位于图像中心；
S33：采用激光测距仪测量得到激光测距仪俯仰角的测量值ωL、激光测距仪横滚角的测量值和以及激光测距仪到物平面的斜距DL1；
S34：根据测量系统检校阶段得到的相机横滚偏心角和相机俯仰偏心角Δω，通过式(15)得到相机的横滚角和俯仰角ωC；
S35：定义O为相机透镜平面的光心，o为光心在像平面上的投影点，O″为光心在物平面上的投影点，计算光心O到物平面的高程H1：
H1＝DL1·sinωL+d·cosωC (16)
S36：对于物平面上的任意点Q，其在畸变校正后的原始像平面上对应像点q的坐标用(u2,v2)表示；首先将图像坐标系的原点由左上角平移至o点，然后采用坐标旋转对相机横滚角进行校正，得到的像平面上像点q′的坐标
其中，w2为图像的宽度，h2为图像的高度；
S37：将相机透镜平面的光心O在物平面上的垂足点C′作为参考点，计算物点Q到点C′的水平距离M′和垂直距离N′，过程如下：定义q″为q′在的像平面主纵线上的投影点，Q″为物点Q在物平面主纵线O″C′上的投影点，根据三角关系得到O到Q″的距离L1为：
L1＝H1/sinβ1 (18)
其中，β1表示射线OQ″和主纵线O″C′间的夹角，有：
其中，s表示已知的相机的像元尺寸，f表示相机的实际焦距，f通过式(20)得到：
f＝s·(fx+fy)/2 (20)
根据三角形的相似性，有：
Q″Q/Q″O＝q′q″/q″O (21)
代入对应参数，得：
同理，有：
S38：对于像平面中任意两点q1和q2，首先根据步骤S37分别求得其物点相对于垂足点C′的距离(M1′,N1′)和(M2′,N2′)，然后代入下式：
其中，M1′为点q1相对于垂足点C′的水平距离，N1′为点q1相对于垂足点C′的垂直距离，M2′为点q2相对于垂足点C′的水平距离，N2′为点q2相对于垂足点C′的垂直距离，D(q1,q2)为点q1和点q2之间的距离。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910275119.2A CN110057295B (zh) | 2019-04-08 | 2019-04-08 | 一种免像控的单目视觉平面距离测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910275119.2A CN110057295B (zh) | 2019-04-08 | 2019-04-08 | 一种免像控的单目视觉平面距离测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110057295A true CN110057295A (zh) | 2019-07-26 |
| CN110057295B CN110057295B (zh) | 2020-12-25 |
Family
ID=67318523
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910275119.2A Active CN110057295B (zh) | 2019-04-08 | 2019-04-08 | 一种免像控的单目视觉平面距离测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110057295B (zh) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110458893A (zh) * | 2019-07-29 | 2019-11-15 | 武汉光庭信息技术股份有限公司 | 高级辅助驾驶视觉感知传感器的横滚角标定方法及系统 |
| CN110530441A (zh) * | 2019-08-21 | 2019-12-03 | 河海大学 | 一种全天候在线视频测流系统 |
| CN110763306A (zh) * | 2019-09-30 | 2020-02-07 | 中国科学院西安光学精密机械研究所 | 一种基于单目视觉的液位测量系统及方法 |
| CN111192235A (zh) * | 2019-12-05 | 2020-05-22 | 中国地质大学（武汉） | 一种基于单目视觉模型和透视变换的图像测量方法 |
| CN111595254A (zh) * | 2020-06-04 | 2020-08-28 | 中国人民解放军陆军装甲兵学院 | 一种透镜阵列与lcd显示屏轴向间距测量方法及系统 |
| CN111666876A (zh) * | 2020-06-05 | 2020-09-15 | 北京百度网讯科技有限公司 | 用于检测障碍物的方法、装置、电子设备和路侧设备 |
| CN111982072A (zh) * | 2020-07-29 | 2020-11-24 | 西北工业大学 | 一种基于单目视觉的目标测距方法 |
| CN112556589A (zh) * | 2020-12-03 | 2021-03-26 | 厦门大学 | 一种高度测量方法和高度缺陷判定方法 |
| CN112802129A (zh) * | 2021-04-13 | 2021-05-14 | 之江实验室 | 一种基于单目视觉的焊接安全距离测量方法 |
| CN113050108A (zh) * | 2021-03-23 | 2021-06-29 | 湖南盛鼎科技发展有限责任公司 | 电子界址视觉测量系统及测量方法 |
| CN113074635A (zh) * | 2021-03-29 | 2021-07-06 | 哈尔滨市科佳通用机电股份有限公司 | 一种标定螺栓及利用该标定螺栓检测螺栓松脱的方法 |
| CN113804916A (zh) * | 2021-09-17 | 2021-12-17 | 河海大学 | 一种基于最大流速先验信息的频域时空图像测速法 |
| CN114034227A (zh) * | 2021-11-03 | 2022-02-11 | 泰山玻璃纤维邹城有限公司 | 玻璃纤维拉丝e距离的测量方法 |
| WO2023284202A1 (zh) * | 2021-07-13 | 2023-01-19 | 北京金博星指纹识别科技有限公司 | 一种图像分辨率归一化处理方法及装置 |
| CN116141674A (zh) * | 2022-12-12 | 2023-05-23 | 合肥工业大学 | 一种基于视觉测距的3d打印机基板自动调平装置及方法 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113124819B (zh) * | 2021-06-17 | 2021-09-10 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种基于平面镜的单目测距方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204388803U (zh) * | 2015-01-19 | 2015-06-10 | 长春师范大学 | 结构光视觉系统中激光器和相机的组合固定装置 |
| CN106197270A (zh) * | 2016-07-07 | 2016-12-07 | 大连理工大学 | 一种便携式龙门视觉测量装置 |
-
2019
- 2019-04-08 CN CN201910275119.2A patent/CN110057295B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204388803U (zh) * | 2015-01-19 | 2015-06-10 | 长春师范大学 | 结构光视觉系统中激光器和相机的组合固定装置 |
| CN106197270A (zh) * | 2016-07-07 | 2016-12-07 | 大连理工大学 | 一种便携式龙门视觉测量装置 |
Non-Patent Citations (1)
| Title |
|---|
| 张振等: "基于物像尺度变换的河流水面流场定标方法", 《仪器仪表学报》 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110458893A (zh) * | 2019-07-29 | 2019-11-15 | 武汉光庭信息技术股份有限公司 | 高级辅助驾驶视觉感知传感器的横滚角标定方法及系统 |
| CN110458893B (zh) * | 2019-07-29 | 2021-09-24 | 武汉光庭信息技术股份有限公司 | 高级辅助驾驶视觉感知传感器的横滚角标定方法及系统 |
| CN110530441A (zh) * | 2019-08-21 | 2019-12-03 | 河海大学 | 一种全天候在线视频测流系统 |
| CN110763306A (zh) * | 2019-09-30 | 2020-02-07 | 中国科学院西安光学精密机械研究所 | 一种基于单目视觉的液位测量系统及方法 |
| CN111192235A (zh) * | 2019-12-05 | 2020-05-22 | 中国地质大学（武汉） | 一种基于单目视觉模型和透视变换的图像测量方法 |
| CN111595254A (zh) * | 2020-06-04 | 2020-08-28 | 中国人民解放军陆军装甲兵学院 | 一种透镜阵列与lcd显示屏轴向间距测量方法及系统 |
| CN111595254B (zh) * | 2020-06-04 | 2021-09-21 | 中国人民解放军陆军装甲兵学院 | 一种透镜阵列与lcd显示屏轴向间距测量方法及系统 |
| CN111666876A (zh) * | 2020-06-05 | 2020-09-15 | 北京百度网讯科技有限公司 | 用于检测障碍物的方法、装置、电子设备和路侧设备 |
| CN111666876B (zh) * | 2020-06-05 | 2023-06-09 | 阿波罗智联(北京)科技有限公司 | 用于检测障碍物的方法、装置、电子设备和路侧设备 |
| CN111982072A (zh) * | 2020-07-29 | 2020-11-24 | 西北工业大学 | 一种基于单目视觉的目标测距方法 |
| CN112556589A (zh) * | 2020-12-03 | 2021-03-26 | 厦门大学 | 一种高度测量方法和高度缺陷判定方法 |
| CN112556589B (zh) * | 2020-12-03 | 2022-01-07 | 厦门大学 | 一种高度测量方法和高度缺陷判定方法 |
| CN113050108A (zh) * | 2021-03-23 | 2021-06-29 | 湖南盛鼎科技发展有限责任公司 | 电子界址视觉测量系统及测量方法 |
| CN113050108B (zh) * | 2021-03-23 | 2024-01-09 | 湖南盛鼎科技发展有限责任公司 | 电子界址视觉测量系统及测量方法 |
| CN113074635A (zh) * | 2021-03-29 | 2021-07-06 | 哈尔滨市科佳通用机电股份有限公司 | 一种标定螺栓及利用该标定螺栓检测螺栓松脱的方法 |
| CN112802129A (zh) * | 2021-04-13 | 2021-05-14 | 之江实验室 | 一种基于单目视觉的焊接安全距离测量方法 |
| WO2023284202A1 (zh) * | 2021-07-13 | 2023-01-19 | 北京金博星指纹识别科技有限公司 | 一种图像分辨率归一化处理方法及装置 |
| CN113804916A (zh) * | 2021-09-17 | 2021-12-17 | 河海大学 | 一种基于最大流速先验信息的频域时空图像测速法 |
| CN114034227A (zh) * | 2021-11-03 | 2022-02-11 | 泰山玻璃纤维邹城有限公司 | 玻璃纤维拉丝e距离的测量方法 |
| CN114034227B (zh) * | 2021-11-03 | 2024-02-06 | 泰山玻璃纤维邹城有限公司 | 玻璃纤维拉丝e距离的测量方法 |
| CN116141674A (zh) * | 2022-12-12 | 2023-05-23 | 合肥工业大学 | 一种基于视觉测距的3d打印机基板自动调平装置及方法 |
| CN116141674B (zh) * | 2022-12-12 | 2024-04-12 | 合肥工业大学 | 一种基于视觉测距的3d打印机基板自动调平装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110057295B (zh) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110057295A (zh) | 一种免像控的单目视觉平面距离测量方法 | |
| CN102376089B (zh) | 一种标靶校正方法及系统 | |
| US7479982B2 (en) | Device and method of measuring data for calibration, program for measuring data for calibration, program recording medium readable with computer, and image data processing device | |
| CN108510551B (zh) | 一种远距离大视场条件下相机参数的标定方法及系统 | |
| CN104240262B (zh) | 一种用于摄影测量的相机外参数标定装置及标定方法 | |
| CN103278138B (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
| JP3728900B2 (ja) | キャリブレーション方法及び装置、並びにキャリブレーション用データ生成方法 | |
| Lucchese | Geometric calibration of digital cameras through multi-view rectification | |
| CN109029299A (zh) | 舱段销孔对接转角的双相机测量装置及测量方法 | |
| CN109855603B (zh) | 一种聚焦测量方法及终端 | |
| CN106595700A (zh) | 基于三点坐标测量的靶道空间基准标定方法 | |
| CN110375648A (zh) | 棋盘格靶标辅助的单台相机实现的空间点三维坐标测量方法 | |
| CN102169573A (zh) | 高精度的宽视场镜头实时畸变矫正方法及系统 | |
| CN103782232A (zh) | 投影仪及其控制方法 | |
| US20130113897A1 (en) | Process and arrangement for determining the position of a measuring point in geometrical space | |
| CN109163657A (zh) | 一种基于双目视觉三维重建的圆形目标位姿检测方法 | |
| CN109900205B (zh) | 一种高精度的单线激光器和光学相机的快速标定方法 | |
| CN110378969A (zh) | 一种基于3d几何约束的汇聚式双目相机标定方法 | |
| CN109934878A (zh) | 一种基于相机坐标系的线性标定系统及其方法 | |
| CN108469254A (zh) | 一种适用于仰视和俯视位姿的大视场多视觉视频测量系统全局标定方法 | |
| CN110736447B (zh) | 一种集成式图像采集设备竖直方向水平位置检校方法 | |
| CN109916304A (zh) | 镜面/类镜面物体三维测量系统标定方法 | |
| CN111879354A (zh) | 一种无人机精细化测量系统 | |
| CN109212497A (zh) | 一种空间六自由度车载雷达天线位姿偏差测量及对接方法 | |
| CN109919856A (zh) | 基于双目视觉的沥青路面构造深度检测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |