CN111536902B - 一种基于双棋盘格的振镜扫描系统标定方法 - Google Patents

一种基于双棋盘格的振镜扫描系统标定方法 Download PDF

Info

Publication number
CN111536902B
CN111536902B CN202010322971.3A CN202010322971A CN111536902B CN 111536902 B CN111536902 B CN 111536902B CN 202010322971 A CN202010322971 A CN 202010322971A CN 111536902 B CN111536902 B CN 111536902B
Authority
CN
China
Prior art keywords
coordinate system
camera
galvanometer
calibration
coordinate
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.)
Active
Application number
CN202010322971.3A
Other languages
English (en)
Other versions
CN111536902A (zh
Inventor
杨树明
王腾
李述胜
袁野
胡鹏宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202010322971.3A priority Critical patent/CN111536902B/zh
Publication of CN111536902A publication Critical patent/CN111536902A/zh
Application granted granted Critical
Publication of CN111536902B publication Critical patent/CN111536902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/2504Calibration devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/002Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/2518Projection by scanning of the object
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30204Marker
    • G06T2207/30208Marker matrix

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Laser Beam Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于双棋盘格的振镜扫描系统标定方法,使用的标定装置为两块规格相同、以一定角度相互固定的棋盘格标定板。该方法首先使用摄像机拍摄若干不同姿态下的标定装置图像,使用图像分割算法对同一图像中的两块标定板进行分割,再根据张正友标定法标定得到摄像机的内参数与外参数。然后保持标定装置的最后一个姿态不变,通过控制振镜镜片转动与摄像机曝光采集图像,获得不同光平面下投射到两块标定板上的激光条纹图像。通过图像处理与计算得到的高精度标定点拟合得到光平面方程,最终获得振镜坐标系与相机坐标系之间的刚体变换,实现振镜扫描系统的整体标定。本发明所需的标定装置容易获得,且标定步骤简单,速度较快,精度较高。

Description

一种基于双棋盘格的振镜扫描系统标定方法
技术领域
本发明属于光学三维扫描系统标定技术领域,具体涉及一种基于双棋盘格的振镜扫描系统标定方法。
背景技术
光学三维测量具有高精度、高效率和非接触等特点,目前已经广泛应用于质量检测,尺寸测量、逆向工程等诸多领域。其中线结构光扫描技术因为其精度高,结构简单,稳定性高等优势,目前已成为应用最广泛的三维测量方法之一。
传统的线结构光扫描系统由光源、摄像机和机械扫描装置组成,根据激光三角发测量原理计算得到单线的三维坐标,通过机械扫描装置移动被测量对象或移动光源完成对物体三维形貌的扫描测量。传统的机械扫描装置体积重量大,且由于惯性较大导致扫扫描速度慢,以上缺点限制了其在工业领域的广泛应用。
振镜扫描作为一种矢量扫描方式,具有体积小,扫描速度快、定位精度高、转动惯量小等优点。使用高精度的振镜扫描器件取代传统的机械扫描装置使得线结构光测量系统具有自扫描的特性,提高了测量系统的集成度。振镜扫描系统通过精确控制反光镜片转动改变激光方向实现高速、高分辨率扫描,同时可以提供很高的定位精度与重复精度,从而实现对物体的高速高精度测量。
振镜扫描系统的标定是确保测量精度的关键步骤,但是关于标定方法的研究较少,现有的标定方法实际操作过程复杂、标定速度较慢,且往往需要依赖高精度的标定装置,难以满足实际生产中现场快速标定的要求。
发明内容
针对现有标定方法存在的问题,本发明的目的在于提供一种操作简单灵活,标定速度快、精度较高、成本较低的基于双棋盘格的振镜扫描系统标定方法。
为达到上述目的,本发明采用如下技术方案来实现的:
一种基于双棋盘格的振镜扫描系统标定方法,包括以下步骤:
1)将两块规格完全相同的棋盘格标定板铰接连接成V字形,形成标定装置,确保两块棋盘格标定板夹角能够在设定范围内变化,且能够在任意角度保持固定;
2)将步骤1)中所述的标定装置固定在某一夹角下,以不通过姿态放置在摄像机视场内拍摄若干张图像,结合张正友标定法对相机镜头进行标定,获得内参数与外参数;
3)驱动振镜镜片转动,拍摄得到不同角度下的激光条纹图像,利用激光条纹的像素坐标与步骤2)中得到的内参数与外参数完成激光平面的拟合标定,实现对整个振镜扫描系统的标定。
本发明进一步的改进在于,所述步骤2)中,拍摄得到的图像先进行图像分割再进行相机镜头的标定,分割后的每幅图像只包含一幅完整的棋盘格阵列图案,且分割后图像的像素坐标系保持不变。
本发明进一步的改进在于,所述步骤2)中,利用张正友标定法标定出摄像机镜头的等效焦距fx与fy,相机的主点坐标u0与v0,镜头畸变系数k1与k2,以及相机坐标系分别与每个姿态下的棋盘格上建立的世界坐标系之间刚体变换矩阵[R1|T1],[R2|T2],....,[Rn|Tn]。
本发明进一步的改进在于,所述步骤3)中,分别建立以下坐标系:以图像左上角为原点,像素的行、列方向为坐标轴方向建立像素坐标系OP-uv,以镜头光轴与图像平面的交点为原点,像素的行、列方向为坐标轴方向建立图像坐标系O1-xy,以镜头光心为原点,像素的行、列方向分别为x轴y轴方向,根据右手系规则确定z轴方向建立相机坐标系OC-XCYCZC,分别以上下两块棋盘格的左上角点为原点,棋盘格的行、列方向为x轴y轴方向,根据右手系规则确定z轴方向建立各自的世界坐标系OU-XUYUZU和OD-XDYDZD,以振镜转轴为y轴方向,y轴与相机坐标系OC-XC-ZC平面的交点为原点,振镜控制电压U=0时的激光平面法线方向为x轴方向,根据右手系规则确定z轴方向建立振镜坐标系OG-XGYGZG,系统标定过程的本质为确定振镜坐标系与相机坐标系之间的相对位置关系。
本发明进一步的改进在于,所述步骤3)中,计算出分别位于两个棋盘格各自世界坐标系下的激光条纹坐标,然后统一转换为相机坐标系下的坐标:
对于固定在上方的棋盘格,根据相机成像模型及世界坐标系的定义有:
Figure BDA0002462141210000031
式中,ρ为空间比例因子,(uU,vU)为激光条纹的像素坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,(XU,YU,ZU)为其在世界坐标系下的坐标,(XUC,YUC,ZUC)为其在相机坐标系下的坐标,RU和TU分别为相机坐标系与上方棋盘格世界坐标系之间的旋转矩阵和平移向量;
Figure BDA0002462141210000032
则MU为4×4矩阵:
Figure BDA0002462141210000033
将MU带入(1)式中求解可得:
Figure BDA0002462141210000041
式中,r11、r12…r33、t1、t2、t3为MU矩阵的元素;
对于固定在下方的棋盘格,令
Figure BDA0002462141210000042
则MD为4×4矩阵:
Figure BDA0002462141210000043
同理计算可得:
Figure BDA0002462141210000044
式中,(uD,vD)为激光条纹的像素坐标,(XDC,YDC,ZDC)为其在相机坐标系下的坐标,r11、r12…r33、t1、t2、t3为MD矩阵的元素。
本发明进一步的改进在于,所述步骤3)中,通过计算双棋盘格上的激光条纹标并拟合光平面计算出振镜坐标系原点以及三个坐标轴在相机坐标系中的位置,进而标定出相机坐标系与振镜坐标系之间的刚体变换矩阵:
Figure BDA0002462141210000045
式中,(XC,YC,ZC)为任意一点在相机坐标系下的坐标,(XG,YG,ZG)为任意一点在振镜坐标系下的坐标,RG为3×3旋转矩阵,TG为3×1的平移向量,0T为1×3的零矩阵。
本发明进一步的改进在于,所述步骤3)中,建立振镜坐标系与像素坐标系之间的转换关系,并结合光平面约束条件得到空间三维坐标与像素二维坐标直接的映射关系,完成系统的整体标定:
根据相机模型与权利要求6中的转换关系(6)可得:
Figure BDA0002462141210000051
式中,ρ为空间比例因子,(u,v)为空间任意一点对应的像素坐标,(XG,YG,ZG)为空间任意一点在振镜坐标系下的坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,RG和TG分别为振镜坐标系与相机坐标系之间的旋转矩阵和平移向量;
根据振镜激光扫描的工作原理,在振镜坐标系下给出光平面方程的约束条件:
Figure BDA0002462141210000052
式中,XG和ZG为光面面上某点在振镜坐标系下的坐标,θ为振镜相对XG轴的机械转动角度,U为振镜扫描系统的控制电压,k为控制电压与振镜机械转动角度之间的比例系数,单位为°/V;进一步整理可得约束条件为:
XG cos2kU-ZG sin2kU=0 (9)
联立式(7)和式(9)得到振镜坐标系下的三维坐标与像素坐标之间的映射关系:
Figure BDA0002462141210000053
本发明至少具有如下有益的技术效果:
本发明提供的一种基于双棋盘格的振镜扫描系统标定方法,将摄像机标定与整个系统的标定有机的结合在一起,利用两块相互固定的的棋盘格标定板回去高精度标志点对光平面进行标定,进一步完成对整个振镜扫描测量系统的标定。相比与其他标定方法,标定器件简单易获得,标定步骤更少,速度更快,且精度较高。
进一步,本发明的标定过程中,每张标定图片包含两张不同角度的棋盘格图案,通过算法对图像进行分割再作为标定算法的输入可以将标定图片数量减少一半,提高了标定的操作速度,在工业化大批量标定时具有明显的优势。
进一步,本发明将两块棋盘格标定板上的激光条纹的三维坐标统一到相机坐标系中,从而可以直接计算出振镜坐标系与相机坐标系之间的转换关系,简化了标定的计算步骤。
进一步,本发明利用两块以一定夹角固定的棋盘格标定出了相机坐标系与振镜坐标系之间的刚体变换矩阵,在完成对相机的标定之后,只需保持标定装置的位姿不变,改变振镜控制电压使得激光平面与双标定板相交,拍摄获得若干不同电压下的多幅激光折线条纹图像,计算并拟合出对应的光平面方程。对光平面的交线进行最小二乘优化,从而精确的标定出振镜坐标系与相机坐标系直接的刚体变换矩阵。
进一步,本发明根据振镜的工作原理以及振镜坐标系的定义,建立了控制电压与激光平面约束条件之间的联系,结合标定得到的振镜坐标系与像素坐标系之间的转换关系,便可得到像素坐标和控制电压与空间三维坐标之间的单一映射关系。在进行系统标定之后,通过上述映射关系便可快速地计算得到被测对象的三维形貌信息。
附图说明
图1是本发明的双棋盘格标定装置的示意图,其中图1中(a)为左视图,图1中(b)为轴侧图。
图2是本发明的坐标系建立的示意图。
图3是本发明振镜坐标系标定的原理图。
具体实施方式
为进一步阐述本发明的目的、操作步骤以及优势,下面结合附图详细说明本发明的实施方式:
本发明提供的一种基于双棋盘格的振镜扫描系统标定方法,包括以下步骤:
步骤一:制作双棋盘格标定装置,使用可在任意位置固定的机械铰链将两块规格相同的棋盘格标定板按图1所示固定连接,两标定板之间的角度可通过旋转铰链改变。以下约定两棋盘格标定板名称分别为上标定板和下标定板。
步骤二:为更加明确的说明摄相机标定与系统标定的过程,根据相机成像模型与振镜扫描系统工作原理建立以下六个坐标系,分别为:像素坐标系、图像坐标系、相机坐标系、上标定板世界坐标系、下标定板世界坐标系和振镜坐标系。以上坐标系的建立如图2所示:
(1)像素坐标系OP-uv:以图像左上角为原点OP,像素的行、列方向分别为u,v坐标轴方向建立二维直角坐标系,坐标轴单位为像素。
(2)图像坐标系O1-xy:以镜头光轴与图像平面的交点为原点O1,平行于像素的行、列方向为x,y坐标轴方向建立二维直角坐标系,坐标轴单位为毫米,用来描述图像中一点的物理位置。
(3)相机坐标系OC-YCZC:以镜头光心为原点OC,平行于像素的行、列方向分别为XC轴和YC轴方向,根据右手坐标系规则确定ZC轴方向建立三维直角坐标系,坐标轴单位为毫米。
(4)上标定板世界坐标系OU-XUYUZU:以上方棋盘格标定板的左上角点为原点OU,棋盘格的行、列方向分别为XU轴和YU轴方向,根据右手系规则确定ZU轴方向建立三维直角坐标系,坐标轴单位为毫米。
(5)下标定板世界坐标系OD-XDYDZD:以下方棋盘格标定板的左上角点为原点OD,棋盘格的行、列方向分别为XD轴和YD轴方向,根据右手系规则确定ZD轴方向建立三维直角坐标系,坐标轴单位为毫米。
(6)振镜坐标系OG-XGYGZG:以振镜转轴为YG轴方向,YG与相机坐标系OC-XC-ZC平面的交点为原点OG,振镜控制电压U=0时的激光平面法线方向为XG轴方向,根据右手系规则确定ZG轴方向建立三维直角坐标系,坐标轴单位为毫米。
步骤三:建立坐标系之间的转换关系。根据相机成像模型,像素坐标系与图像坐标系之间的转换关系由原点坐标偏移和坐标单位之间的转换确定:
Figure BDA0002462141210000081
式中,(u,v)为像素坐标系中的点,(x,y)为图像坐标系中的点,dx和dy为单个像素在x和y方向上的物理尺寸,单位为毫米,(u0,v0)为图像坐标系原点在像素坐标系中的坐标,即相机的主点坐标。
根据相机的小孔成像模型与透视投影原理,图像坐标与相机坐标之间的转换关系为:
Figure BDA0002462141210000082
式中,(x,y)为图像坐标系中的点,(XC,YC,YC)为相机坐标系中的点,f为相机镜头的物理焦距,单位为毫米,ρ为空间比例因子。
以相机坐标系为纽带,建立上述坐标系之间的转换关系。联立式(1)和式(2)得到相机坐标系与像素坐标系之间的转换关系:
Figure BDA0002462141210000083
式中,ρ为空间比例因子,(u,v)为像素坐标系中的点,(u0,v0)为图像坐标系原点在像素坐标系中的坐标,fx=f/dx和fy=f/dy为相机镜头的等效焦距,单位为1,(XC,YC,YC)为相机坐标系中的点。
根据坐标系之间的刚体变换原理,上标定板坐标系与相机坐标系之间的转化关系为:
Figure BDA0002462141210000091
式中,(XU,YU,ZU)为点在上标定板世界坐标系下的坐标,(XC,YC,ZC)为点在相机坐标系下的坐标,RU为3×3旋转矩阵,TU为3×1的平移向量,0T为1×3的零矩阵,RU和TU可通过相机标定计算得出。
同理可得下标定板坐标系与相机坐标系之间的转化关系为:
Figure BDA0002462141210000092
式中,(XD,YD,ZD)为点在下标定板世界坐标系下的坐标,(XC,YC,ZC)为点在相机坐标系下的坐标,RD为3×3旋转矩阵,TD为3×1的平移向量,0T为1×3的零矩阵,RD和TD可通过相机标定计算得出。
根据坐标系之间的刚体变换原理,振镜坐标系与相机坐标系之间的转化关系为:
Figure BDA0002462141210000093
式中,(XC,YC,ZC)为点在相机坐标系下的坐标,(XG,YG,ZG)为点在振镜坐标系下的坐标,RG为3×3旋转矩阵,TG为3×1的平移向量,0T为1×3的零矩阵。RG和TG即为待标定的参数,通过系统标定计算得出。
步骤四:摄相机标定。将步骤一中的双棋盘格标定板放置在相机视场中,在不同位置和不同角度拍摄10张图像,每张图像均包含两幅棋盘格图案且曝光正常。使用图像处理算法对采集到的图像进行分割,分割后的两张图像分别只包含一副完整的棋盘格图案并保持原始的图像尺寸不变。对图片进行从1到20的编号,其中奇数序号图片均为上标定板图案,偶数序号图片均为下标定板图案。
将得到的20张图片作为张正友标定法的输入进行计算,计算得到相机的内参矩阵K、畸变参数向量D和20组外参矩阵[R1|T1],[R2|T2],....,[R20|T20]。其中K即为相机坐标系与像素坐标系之间的转换矩阵,D为包含三阶径向畸变参数与二阶切向畸变参数的5×1向量,每组外参矩阵[R|T]对应相机坐标与不同姿态下的棋盘格世界坐标系之间的刚体变换矩阵。
步骤五:获取相机坐标系下的标志点坐标。计算不同角度的光平面下激光条纹的三维坐标并将其转换到相机坐标系下:
(1)获取系统标定所需的图像。保持相机标定时标定装置的最后一个位姿不变,打开激光器,改变控制电压使得振镜转动一定角度,反射激光同时照射在两个标定板上,调整相机曝光采集此时的激光条纹图像,将振镜控制电压依次设置为U1,U2,…,Un(n≥4)并分别拍摄n幅图像。
(2)计算激光条纹三维坐标。如图3所示,操作(1)中采集的图像中激光条纹为一条折线,折线的两段分别位于上下两块标定板上。首先利用相机标定时计算出的畸变参数对图像进行畸变矫正,然后利用霍夫变换确定折线条纹的顶点并分割折线条纹,最后使用中心线提取算法分别提取上下标定板上的条纹像素坐标(uU,vU)和(uD,vD)。对于上标定板的上激光条纹,其在上标定板世界坐标系中有ZU=0,结合步骤三中的式(3)和式(4)有:
Figure BDA0002462141210000101
式中,ρ为空间比例因子,(u0,v0)为图像坐标系原点在像素坐标系中的坐标,fx和fy为相机镜头的等效焦距,R19和T19分别为标定板最后一个位姿下,上标定板坐标系相对相机坐标系的旋转矩阵和平移向量,(XU,YU,ZU)为点在上标定板世界坐标系下的坐标。内参矩阵与外参矩阵均已在相机标定时计算得到,因此可以计算出不同控制电压下激光条纹的三维坐标(XU,YU,ZU=0)。
同理可以计算得到不同控制电压下,下标定板上激光条纹的三维坐标(XD,YD,ZD=0)。
(3)将激光条纹三维坐标统一转换到相机坐标系。上述计算得到的激光条纹三维坐标分别位于空间中两相交直线上,因此可以唯一的确定一个激光平面方程。但操作(2)中计算得到的三维坐标分别位于上下标定板各自的世界坐标系中,因此需要将其统一转换到相机坐标系中:
对于上标定板的激光条纹,坐标转换可通过步骤三中式(4)完成,令:
Figure BDA0002462141210000111
式中,RU=R19和TU=T19分别为标定板最后一个位姿下,上标定板坐标系相对相机坐标系的旋转矩阵和平移向量。
将MU带入式(4)中,并带入三维坐标(XU,YU,ZU=0)的具体计算结果,得:
Figure BDA0002462141210000112
式中,(uU,vU)为激光条纹的像素坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,(XU,YU,ZU)为激光条纹在世界坐标系下的坐标,(XUC,YUC,ZUC)为激光条纹在相机坐标系下的坐标,r11、r12…r33、t1、t2、t3为MU矩阵的元素。
对于固定在下方的棋盘格,令:
Figure BDA0002462141210000121
式中,RD=R20和TD=T20分别为标定板最后一个位姿下,下标定板坐标系相对相机坐标系的旋转矩阵和平移向量。
同理计算可得:
Figure BDA0002462141210000122
式中,(uD,vD)为激光条纹的像素坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,(XU,YU,ZU)为激光条纹在世界坐标系下的坐标,(XDC,YDC,ZDC)为其在相机坐标系下的坐标,r11、r12…r33、t1、t2、t3为矩阵MD的元素。
步骤六:标定振镜坐标系与相机坐标系的转换关系。
(1)拟合得到不同控制电压下的激光平面。当控制电压为U1时,设对应激光平面在相机坐标系下的方程为:
A1x+B1y+C1z+D1=0 (11)
将步骤五中计算所得的激光条纹坐标(XUC,YUC,ZUC)、(XDC,YDC,ZDC)带入式(11)中并进行最小二乘拟合,得到该光平面的单位法向量n1=(i1,j1,k1)。同理可计算得到控制电压为U2,U3,…,Un(n≥4)时对应光平面的法向量n1=(i2,j2,k2),n1=(i3,j3,k3),…,nn=(in,jn,kn)。
(2)设振镜坐标系XG轴在相机坐标系下的单位方向向量为nx=(ax,bx,cx),根据振镜坐标系的定义以及控制电压与转动角度之间的关系可知:
Figure BDA0002462141210000131
式中,U1、U2、U3为振镜控制电压,k为控制电压与振镜转角的比例系数,求解式(12)可得出nx
设振镜坐标系YG轴在相机坐标系下的单位方向向量为ny=(ay,by,cy),由定义知,ny垂直于所有光平面的法向量,因此有:
ni·ny=0(i=1,2,...,n) (13)
由于n≥4,因此式(13)为一超定方程组,求解该方程组便可得出ny
设振镜坐标系ZG轴在相机坐标系下的单位方向向量为nz=(az,bz,cz),根据三坐标轴互相垂直可得nz=nx×ny
(3)理想情况下,所有激光平面应相交于YG轴所在的直线,但由于激光器安装误差,实际情况为YG轴上存在一点P(x0,y0,z0)到所有光平面的距离最近,联立光平面方程建立方程组:
Aix+Biy+Ciz+Di=0,i=1,2,...,n (14)
由于n≥4,因此式(14)为一超定方程组,该方程组的最小二乘解即为点坐标。
根据YG轴的方向向量ny和轴上一点P可得YG轴所在直线在相机坐标系下的方程:
Figure BDA0002462141210000132
设振镜坐标系的原点在相机坐标系中的坐标为OG(xG,yG,zG),由定义知,原点OG为YG轴与相机坐标系中OC-XC-ZC平面的交点,故有yG=0,联立YG轴所在直线方程得OG坐标为:
Figure BDA0002462141210000133
式中,(x0,y0,z0)为YG轴上P点的坐标,(ay,by,cy)YG轴的单位方向向量。
至此完成了振镜坐标系原点坐标以及坐标轴方向向量的计算,再根据步骤三中式(6)可标定出振镜坐标系与相机坐标系之间的转换矩阵为:
Figure BDA0002462141210000141
步骤七:结合光平面约束条件完成整体标定。联立步骤三中式(3)和式(6)可的振镜坐标系与像素坐标系之间的映射关系:
Figure BDA0002462141210000142
式中,ρ为空间比例因子,(u,v)为像素坐标系中的点,(u0,v0)为图像坐标系原点在像素坐标系中的坐标,fx和fy为相机镜头的等效焦距,(XG,YG,ZG)为点在振镜坐标系下的坐标,RG为3×3旋转矩阵,TG为3×1的平移向量,0T为1×3的零矩阵。
由振镜工作原理知,振镜控制电压U与机械转角θ之间为线性关系,设其比例系数为k,则有θ=kU。又根据反射定律,激光平面转动的角度
Figure BDA0002462141210000143
与机械转角θ之间的关系为
Figure BDA0002462141210000144
因此可得出振镜坐标系下的激光平面约束条件为:
XG cos2kU-ZG sin2kU=0 (19)
联立式(18)和式(19)便可完成系统的整体标定:
Figure BDA0002462141210000145
若已知一组振镜控制电压和对应的激光条纹像素坐标,便可通过求解式(20)得到振镜坐标系下对应的三维坐标。
以上结合附图对本发明的具体实施方式作了说明,但这些说明不能被理解为限制了本发明的范围,本发明的保护范围由随附的权利要求书限定,任何在本发明权利要求基础上的改动都是本发明的保护范围。

Claims (4)

1.一种基于双棋盘格的振镜扫描系统标定方法,其特征在于,包括以下步骤:
1)将两块规格完全相同的棋盘格标定板铰接连接成V字形,形成标定装置,确保两块棋盘格标定板夹角能够在设定范围内变化,且能够在任意角度保持固定;
2)将步骤1)中所述的标定装置固定在某一夹角下,以不同的姿态放置在摄像机视场内拍摄若干张图像,结合张正友标定法对相机镜头进行标定,获得内参数与外参数;拍摄得到的图像先进行图像分割再进行相机镜头的标定,分割后的每幅图像只包含一幅完整的棋盘格阵列图案,且分割后图像的像素坐标系保持不变;利用张正友标定法标定出摄像机镜头的等效焦距fx与fy,相机的主点坐标u0与v0,镜头畸变系数k1与k2,以及相机坐标系分别与每个姿态下的棋盘格上建立的世界坐标系之间刚体变换矩阵[R1|T1],[R2|T2],....,[Rn|Tn];
3)驱动振镜镜片转动,拍摄得到不同角度下的激光条纹图像,利用激光条纹的像素坐标与步骤2)中得到的内参数与外参数完成激光平面的拟合标定,实现对整个振镜扫描系统的标定;分别建立以下坐标系:以图像左上角为原点,像素的行、列方向为坐标轴方向建立像素坐标系OP-uv,以镜头光轴与图像平面的交点为原点,像素的行、列方向为坐标轴方向建立图像坐标系O1-xy,以镜头光心为原点,像素的行、列方向分别为x轴y轴方向,根据右手系规则确定z轴方向建立相机坐标系OC-XCYCZC,分别以上下两块棋盘格的左上角点为原点,棋盘格的行、列方向为x轴y轴方向,根据右手系规则确定z轴方向建立各自的世界坐标系OU-XUYUZU和OD-XDYDZD,以振镜转轴为y轴方向,y轴与相机坐标系OC-XC-ZC平面的交点为原点,振镜控制电压U=0时的激光平面法线方向为x轴方向,根据右手系规则确定z轴方向建立振镜坐标系OG-XGYGZG,系统标定过程的本质为确定振镜坐标系与相机坐标系之间的相对位置关系。
2.根据权利要求1中所述的一种基于双棋盘格的振镜扫描系统标定方法,其特征在于,所述步骤3)中,计算出分别位于两个棋盘格各自世界坐标系下的激光条纹坐标,然后统一转换为相机坐标系下的坐标:
对于固定在上方的棋盘格,根据相机成像模型及世界坐标系的定义有:
Figure FDA0002930264130000021
式中,ρ为空间比例因子,(uU,vU)为激光条纹的像素坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,(XU,YU,ZU)为其在世界坐标系下的坐标,(XUC,YUC,ZUC)为其在相机坐标系下的坐标,RU和TU分别为相机坐标系与上方棋盘格世界坐标系之间的旋转矩阵和平移向量;
Figure FDA0002930264130000022
则MU为4×4矩阵:
Figure FDA0002930264130000023
将MU带入(1)式中求解可得:
Figure FDA0002930264130000024
式中,r11、r12…r33、t1、t2、t3为MU矩阵的元素;
对于固定在下方的棋盘格,令
Figure FDA0002930264130000025
则MD为4×4矩阵:
Figure FDA0002930264130000031
同理计算可得:
Figure FDA0002930264130000032
式中,(uD,vD)为激光条纹的像素坐标,(XDC,YDC,ZDC)为其在相机坐标系下的坐标,r11、r12…r33、t1、t2、t3为MD矩阵的元素。
3.根据权利要求2中所述的一种基于双棋盘格的振镜扫描系统标定方法,其特征在于,所述步骤3)中,通过计算双棋盘格上的激光条纹标并拟合光平面计算出振镜坐标系原点以及三个坐标轴在相机坐标系中的位置,进而标定出相机坐标系与振镜坐标系之间的刚体变换矩阵:
Figure FDA0002930264130000033
式中,(XC,YC,ZC)为任意一点在相机坐标系下的坐标,(XG,YG,ZG)为任意一点在振镜坐标系下的坐标,RG为3×3旋转矩阵,TG为3×1的平移向量,0T为1×3的零矩阵。
4.根据权利要求3中所述的一种基于双棋盘格的振镜扫描系统标定方法,其特征在于,所述步骤3)中,建立振镜坐标系与像素坐标系之间的转换关系,并结合光平面约束条件得到空间三维坐标与像素二维坐标直接的映射关系,完成系统的整体标定:
根据相机模型与权利要求3中的转换关系(6)可得:
Figure FDA0002930264130000041
式中,ρ为空间比例因子,(u,v)为空间任意一点对应的像素坐标,(XG,YG,ZG)为空间任意一点在振镜坐标系下的坐标,fx与fy为摄像机镜头的等效焦距,(u0,v0)为相机的主点坐标,RG和TG分别为振镜坐标系与相机坐标系之间的旋转矩阵和平移向量;
根据振镜激光扫描的工作原理,在振镜坐标系下给出光平面方程的约束条件:
Figure FDA0002930264130000042
式中,XG和ZG为光面面上某点在振镜坐标系下的坐标,θ为振镜相对XG轴的机械转动角度,U为振镜扫描系统的控制电压,k为控制电压与振镜机械转动角度之间的比例系数,单位为°/V;进一步整理可得约束条件为:
XGcos2kU-ZGsin2kU=0 (9)
联立式(7)和式(9)得到振镜坐标系下的三维坐标与像素坐标之间的映射关系:
Figure FDA0002930264130000043
CN202010322971.3A 2020-04-22 2020-04-22 一种基于双棋盘格的振镜扫描系统标定方法 Active CN111536902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010322971.3A CN111536902B (zh) 2020-04-22 2020-04-22 一种基于双棋盘格的振镜扫描系统标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010322971.3A CN111536902B (zh) 2020-04-22 2020-04-22 一种基于双棋盘格的振镜扫描系统标定方法

Publications (2)

Publication Number Publication Date
CN111536902A CN111536902A (zh) 2020-08-14
CN111536902B true CN111536902B (zh) 2021-03-09

Family

ID=71980112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010322971.3A Active CN111536902B (zh) 2020-04-22 2020-04-22 一种基于双棋盘格的振镜扫描系统标定方法

Country Status (1)

Country Link
CN (1) CN111536902B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202100024218A1 (it) * 2021-09-21 2023-03-21 Flexsight S R L Apparato automatizzato di calibrazione di una pluralità di sensori di immagine appartenenti ad uno scanner 3d e metodo di calibrazione di detti sensori di immagine mediante detto apparato automatizzato di calibrazione

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112132906B (zh) * 2020-09-22 2023-07-25 西安电子科技大学 一种深度相机与可见光相机之间的外参标定方法及系统
CN114638899B (zh) * 2020-12-16 2024-09-27 宇通客车股份有限公司 利用相机提取棋盘格标定板特征点及对应的联合标定方法
CN112700480B (zh) * 2020-12-29 2022-07-12 河北工业大学 一种面向小尺寸物体旋转扫描的点云快速配准方法及应用
CN112884845B (zh) * 2021-01-14 2023-03-31 山东浪潮科学研究院有限公司 一种基于单摄像机室内机器人障碍物定位的方法
CN112927133B (zh) * 2021-02-07 2022-04-26 湖南桥康智能科技有限公司 一种基于一体化标定参数的图像空间投影拼接方法
CN113177989B (zh) * 2021-05-07 2024-07-19 深圳云甲科技有限公司 一种口内扫描仪标定方法及装置
CN113393538B (zh) * 2021-07-14 2024-05-28 南京苏润科技发展有限公司 一种基于双棋盘格标定的闸门门缝检测方法
CN113655466B (zh) * 2021-07-16 2023-09-29 国家石油天然气管网集团有限公司 一种基于结构光三角测量的标定系统和标定方法
CN113744344B (zh) * 2021-08-18 2023-09-08 富联裕展科技(深圳)有限公司 激光设备的标定方法、装置、设备及存储介质
CN113989369A (zh) * 2021-10-21 2022-01-28 霍夫纳格智能科技(嘉兴)有限公司 一种用于激光加工系统高精度的标定方法及装置
CN116109492A (zh) * 2021-11-10 2023-05-12 梅卡曼德(北京)机器人科技有限公司 条纹光强不均匀的改善方法、装置、设备、介质及相机
CN114608806A (zh) * 2022-02-15 2022-06-10 德特克(上海)数字技术有限责任公司 一种基于双目相机的激光振镜扫描系统的标定方法
CN114923410B (zh) * 2022-05-09 2024-05-14 一汽解放汽车有限公司 一种纵梁孔位在线检测方法及装置
CN116485918B (zh) * 2023-06-25 2023-09-08 天府兴隆湖实验室 一种标定方法、系统及计算机可读存储介质
CN117553682B (zh) * 2023-11-03 2024-07-26 苏州福斯特万电子科技有限公司 五轴点胶设备的转台中心标定方法及系统
CN117606747B (zh) * 2023-11-22 2024-06-18 北京天翔睿翼科技有限公司 激光振镜系统的高精度化标定方法
CN117422776A (zh) * 2023-12-14 2024-01-19 西北工业大学青岛研究院 一种非参数畸变模型的水下视觉相机标定装置及方法
CN118134992A (zh) * 2024-01-27 2024-06-04 西安工业大学 一种非规则物体尺寸测量方法及装置
CN118067015B (zh) * 2024-04-18 2024-07-09 钛玛科(北京)工业科技有限公司 一种基于嵌入式线阵相机的测宽系统
CN118397108B (zh) * 2024-06-26 2024-08-27 中国海洋大学 一种用于水下声光信息结合的标定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07120224A (ja) * 1993-10-22 1995-05-12 Serutetsuku Syst:Kk 形鋼の寸法測定装置
CN102941410A (zh) * 2012-10-18 2013-02-27 北京航空航天大学 一种点扫描三维测量系统振镜标定方法
CN106247979A (zh) * 2016-08-17 2016-12-21 四川川大智胜软件股份有限公司 一种基于振镜的三维扫描仪
CN107014312A (zh) * 2017-04-25 2017-08-04 西安交通大学 一种振镜式线激光扫描三维测量系统的整体标定方法
CN110375648A (zh) * 2019-08-05 2019-10-25 华南农业大学 棋盘格靶标辅助的单台相机实现的空间点三维坐标测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07120224A (ja) * 1993-10-22 1995-05-12 Serutetsuku Syst:Kk 形鋼の寸法測定装置
CN102941410A (zh) * 2012-10-18 2013-02-27 北京航空航天大学 一种点扫描三维测量系统振镜标定方法
CN106247979A (zh) * 2016-08-17 2016-12-21 四川川大智胜软件股份有限公司 一种基于振镜的三维扫描仪
CN107014312A (zh) * 2017-04-25 2017-08-04 西安交通大学 一种振镜式线激光扫描三维测量系统的整体标定方法
CN107014312B (zh) * 2017-04-25 2019-05-24 西安交通大学 一种振镜式线激光扫描三维测量系统的整体标定方法
CN110375648A (zh) * 2019-08-05 2019-10-25 华南农业大学 棋盘格靶标辅助的单台相机实现的空间点三维坐标测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
结冰风洞中冰形生长的光学三维测量研究;龚晓亮;《中国优秀博士论文全文集工程科技Ⅱ辑》;20170430(第4期);文章第2.2.3节-2.2.5节、2.3整节、附图2-13、2-15 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202100024218A1 (it) * 2021-09-21 2023-03-21 Flexsight S R L Apparato automatizzato di calibrazione di una pluralità di sensori di immagine appartenenti ad uno scanner 3d e metodo di calibrazione di detti sensori di immagine mediante detto apparato automatizzato di calibrazione

Also Published As

Publication number Publication date
CN111536902A (zh) 2020-08-14

Similar Documents

Publication Publication Date Title
CN111536902B (zh) 一种基于双棋盘格的振镜扫描系统标定方法
CN107014312B (zh) 一种振镜式线激光扫描三维测量系统的整体标定方法
CN106408556B (zh) 一种基于一般成像模型的微小物体测量系统标定方法
CN108921901B (zh) 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法
CN106056620B (zh) 线激光相机测量系统标定方法
CN110966935B (zh) 基于标志点的偏折测量系统一体化几何标定方法
CN106871787B (zh) 大空间线扫描成像三维测量方法
CN103530880B (zh) 基于投影高斯网格图案的摄像机标定方法
CN109099883A (zh) 高精度大视场机器视觉测量与标定装置及方法
CN109615661A (zh) 光场相机内参数标定装置及方法
CN110455225B (zh) 基于结构光视觉的矩形花键轴同轴度及键位置度测量方法
CN109272555B (zh) 一种rgb-d相机的外部参数获得及标定方法
CN110455198B (zh) 基于线结构光视觉的矩形花键轴键宽及直径测量方法
CN113205603A (zh) 一种基于旋转台的三维点云拼接重建方法
CN109141226A (zh) 单相机多角度的空间点坐标测量方法
CN115861445B (zh) 一种基于标定板三维点云的手眼标定方法
CN112229323B (zh) 基于手机单目视觉的棋盘格合作目标的六自由度测量方法及其应用
CN110501026B (zh) 基于阵列星点的相机内方位元素标定装置及方法
CN112082511A (zh) 一种基于转台的大物体快速标定测量方法
CN110136068A (zh) 基于双侧远心镜头相机之间位置标定的音膜球顶装配系统
CN115187612A (zh) 一种基于机器视觉的平面面积测量方法、装置及系统
CN110827359A (zh) 基于棋盘格三面体的相机与激光器外参检校方法及装置
CN112700480B (zh) 一种面向小尺寸物体旋转扫描的点云快速配准方法及应用
CN114170321A (zh) 一种基于测距的相机自标定方法及系统
CN110108203A (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