CN116182783A - 一种散斑纹理球形靶标及其位姿估计方法 - Google Patents

一种散斑纹理球形靶标及其位姿估计方法 Download PDF

Info

Publication number
CN116182783A
CN116182783A CN202310218505.4A CN202310218505A CN116182783A CN 116182783 A CN116182783 A CN 116182783A CN 202310218505 A CN202310218505 A CN 202310218505A CN 116182783 A CN116182783 A CN 116182783A
Authority
CN
China
Prior art keywords
spherical target
digital
camera
speckle pattern
target
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.)
Pending
Application number
CN202310218505.4A
Other languages
English (en)
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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN202310218505.4A priority Critical patent/CN116182783A/zh
Publication of CN116182783A publication Critical patent/CN116182783A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提出一种散斑纹理球形靶标及其位姿估计方法,属于航洋定位技术领域,散斑纹理球形靶标为三维球体,球体表面均匀覆盖黑白相间的随机散斑图案,方法是将散斑纹理球形靶标刚体连接到待测物体上;建立相机小孔成像模型;建立与散斑纹理球形靶标相同纹理的数字球形靶标;通过相机对散斑纹理球形靶标进行成像,确认数字球形靶标的位置,固定数字球形靶标的球心,按一定规则旋转数字球形靶标,将数字球形靶标图像与实际球形靶标图像作对比,求取像素级汉明距离,在经过多次旋转之后,从中求取出最小的汉明距离,最小汉明距离对应的数字球形靶标图像与实际球形靶标图像是一致的,获得了相机和散斑纹理球形靶标的位置关系,从而获得待测物体的位姿。

Description

一种散斑纹理球形靶标及其位姿估计方法
技术领域
本发明属于航洋定位技术领域,具体地,涉及一种散斑纹理球形靶标及其位姿估计方法。
背景技术
目前,在位姿监测应用上,简单几何图形特征的靶标使用广泛,这类靶标多为十字形、圆形、方形等对称图形,识别算法简单,但识别精度低,且识别结果不具有唯一性。还有诸如棋盘格标定板式的标定靶。采用棋盘格标定板进行标定时,为了将标定板图像上检测到的角点坐标与其世界坐标对应上,必须保证全部角点都能同时被检测到。同时,当图像畸变太大或者有一定的离焦模糊时,角点检测的精度会明显下降。然后是基于棋盘格标定板的变种,如将棋盘格的方块改为圆阵列。虽然可以在一定程度上解决棋盘格标定板存在的问题。但是其需要拍摄大量的图片来标定相机,且对光线的要求高。相对复杂的靶标设计通常由几种简单的几何图形组合而成,能提升一定的识别精度,且后两个靶标的识别还具有唯一性,但对识别出来的图形的完整性有一定要求。
在测量目标物体的位姿时,通常首先需要在物体上部署由多个特征点组成的靶标,再通过相机对靶标进行成像从而测算目标的位姿。这种靶标在应用时经常出现下列问题:
(1)靶标相对相机角度过大时,成像畸变严重。此时从照片中提取靶标特征点信息误差较大,严重影响了位姿估计的精度。
(2)在监测目标位姿时,靶标随着目标随机运动,常出现靶标特征点被遮挡、移除相机视野等问题,导致靶标特征点缺失,使目标位姿监测任务失败。
发明内容
针对以上问题,本发明的目的是提供了一种用于目标位姿估计的散斑纹理球形靶标,并提供了基于该靶标的目标位姿估计方法,散斑纹理球形靶标的球体表面印有随机散斑图案,在进行目标姿态估计时,通过相机对散斑纹理球形靶标进行成像,对比散斑纹理球形靶标成像图片和模拟数字球形靶标成像图片来确定靶标的位姿。
本发明为实现上述目的所采用的技术方案是:
本发明提出了一种散斑纹理球形靶标,该靶标是一个半径为R的三维球体,球体表面均匀覆盖黑白相间的随机散斑图案,且任意两个视角下球形靶标散斑图案均不重复。
本发明还提出了一种散斑纹理球形靶标的位姿估计方法,该方法包括如下步骤:
步骤1、将散斑纹理球形靶标通过刚体连接方式固定到无法通过相机直接测量位姿的待测物体上;
在测量位姿之前,使用三维扫描仪对散斑纹理球形靶标与待测物体整体进行三维扫描,建立散斑纹理球形靶标与待测物体整体的计算机三维模型;
步骤2、建立相机小孔成像模型,根据相机小孔成像模型,得到散斑纹理球形靶标上任意点P从相机系下坐标转换为图像平面I下坐标的过程表示为式(1):
Figure BDA0004115837750000021
式中,K为相机的内参数矩阵,其中f表示相机焦距,dx和dy分别表示图像平面I上每个像素在相机的像平面C的x轴和y轴方向上的物理大小;u0和v0分别表示二维像平面坐标系O1-xy的原点O1在图像平面I上的坐标;(XC,YC,ZC)表示任意点P在相机坐标系O-XCYCZC下的坐标,u和v分别表示点P在图像平面I上的坐标;建立相机的计算机数字模型;
步骤3、使用三维扫描仪对散斑纹理球形靶标进行三维扫描,以建立与散斑纹理球形靶标相同纹理的数字球形靶标;
步骤4、通过相机对散斑纹理球形靶标进行成像,根据拍摄的散斑纹理球形靶标图像,确认数字球形靶标的位置,固定数字球形靶标的球心,按一定规则旋转数字球形靶标,每旋转一次靶标,通过模拟相机拍摄数字球形靶标,将数字球形靶标图像与实际球形靶标图像作对比,求取像素级汉明距离,在经过多次旋转之后,获得一系列汉明距离,从中求取出最小的汉明距离,最小汉明距离下对应的数字球形靶标图像与实际球形靶标图像是一致的,因此便获得了相机和散斑纹理球形靶标的位置关系,也即获得了待测物体的位姿。
进一步,在步骤2中,相机小孔成像模型的建立过程如下:
以相机光心O为原点,建立三维相机坐标系O-XCYCZC,在相机的像平面C中,以点O1为原点,建立二维像平面坐标系O1-xy;在图像平面I中以图像左上角Oi为原点,建立二维图像平面坐标系Oi-uv;点P(XC,YC,ZC)为散斑纹理球形靶标上任意一点,线段Pp经过相机光心O与相机的像平面C交于点p(xC,yC)。
进一步,步骤4包括如下过程:
通过相机对散斑纹理球形靶标进行成像,根据拍摄的散斑纹理球形靶标图像,确认数字球形靶标的位置;
固定数字球形靶标的球心,按如下规则旋转数字球形靶标:
a、相机正对数字球形靶标的x轴方向;
b、将数字球形靶标沿着z轴顺时针根据实际需求转动一个角度α,α<2°,每转动一次拍摄一张数字球形靶标照片,并记录转动角度,直至旋转一周;
c、将数字球形靶标绕y轴顺时针转动一个角度α,α<2°,然后重复步骤b;
d、将数字球形靶标绕y轴逆时针转动一个角度α,α<2°,然后重复步骤b;
至此,便获得了数字球形靶标的整个正面纹理图像,将所有数字球形靶标图像与实际球形靶标图像一一作对比,计算每一次旋转的像素级汉明距离d,图像对比之前,将RGB图像转为灰度图像,则只需要比较两幅图像的一个灰度通道,如式(2)所示;
Figure BDA0004115837750000041
式中i表示第i个像素,l表示像素的总数量,
Figure BDA0004115837750000042
表示实际球形靶标图像中像素i的灰度值,/>
Figure BDA0004115837750000043
表示数字球形靶标图像中像素i的灰度值;
然后从所有的d中求取出最小的汉明距离dmin,最小汉明距离dmin下对应的数字球形靶标图像与实际球形靶标图像是一致的,因此便获得了相机和散斑纹理球形靶标的位置关系,也即获得了待测物体的位姿。
通过上述设计方案,本发明可以带来如下有益效果:本发明提出的一种散斑纹理球形靶标及其位姿估计方法,散斑纹理球形靶标的球体表面印有随机散斑图案,在进行目标姿态估计时,通过相机对散斑纹理球形靶标进行成像,对比散斑纹理球形靶标成像图片和模拟数字球形靶标成像图片来确定靶标的位姿,不需要比对整个靶标来确认两张图片是否一致,对识别精度有很高提升,鲁棒性也有很高提升。相对于一些对称靶标,散斑纹理球形靶标上散斑的随机性不会让拍摄的图片出现歧义,能唯一确认靶标的位姿。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明示意性实施例及其说明用于理解本发明,并不构成本发明的不当限定,在附图中:
图1为散斑纹理球形靶标主视图;
图2为小孔成像模型图;
图3为数字球形靶标旋转示意图;
图4为采用本发明的方法进行水下仪器位姿监测的实施图;
图5为数字球形靶标旋转到某个角度时与实际球形靶标图像对比之实际球形靶标图像示意图;
图6为数字球形靶标旋转到某个角度时与实际球形靶标图像对比之某个角度下的数字球形靶标图;
图7为数字球形靶标与实际球形靶标匹配成功之实际球形靶标图像示意图;
图8为数字球形靶标与实际球形靶标匹配成功之最小的汉明距离对应的数字球形靶标图;
图9为数字球形靶标与实际球形靶标匹配成功的示意图。
具体实施方式
一种散斑纹理球形靶标,该靶标是一个半径为R的三维球体,球体表面均匀覆盖黑白相间的随机图像。当任意两个视角下散斑纹理球形靶标散斑图案的归一化互相关系数均小于阈值0.95时,证明任意视角下随机散斑图案是唯一的。
散斑纹理球形靶标的球体主视图如图1所示。
Figure BDA0004115837750000061
Figure BDA0004115837750000062
Figure BDA0004115837750000063
Figure BDA0004115837750000064
Figure BDA0004115837750000065
式中:I1和I2分别表示两幅灰度图像;m和n分别表示图像像素点的行数和列数;i和j表示第i行第j列的像素点;μ1和μ2分别表示I1和I2灰度均值;δ1和δ2分别表示I1和I2灰度标准差;p(I1,I2)表示I1和I2的归一化互相关系数。
散斑的设计:假设散斑纹理球形靶标的半径为R,要保证球面上任意两个视角下图案内容均不重复。这一步可以通过球体靶标数字模型两个视角下的图像归一化互相关系数来验证。
散斑纹理球形靶标的使用方法:
1、将散斑纹理球形靶标通过刚体连接固定到无法通过相机直接测量位姿(位置和姿态)的物体(例如水下物体等)上面。由于是刚性连接,所以散斑纹理球形靶标的位姿与被测物体的位姿是唯一对应的。然后在被测物体工作之前,使用三维扫描仪对散斑纹理球形靶标与待测物体整体进行三维扫描,建立散斑纹理球形靶标与待测物体整体的计算机三维模型。当确定了散斑纹理球形靶标的位姿后,通过计算机三维模型,即可确定待测物体的位姿;
2、建立相机小孔成像模型,如图2所示,以相机光心O为原点,建立三维相机坐标系O-XCYCZC;在相机的像平面C中,以点O1为原点,建立二维像平面坐标系O1-xy;在图像平面I中以图像左上角Oi为原点,建立二维图像平面坐标系Oi-uv;点P(XC,YC,ZC)为散斑纹理球形靶标上任意一点,线段Pp经过相机光心O与相机的像平面C交于点p(xC,yC)。根据小孔成像模型,散斑纹理球形靶标上任意点P从相机系下坐标转换为图像平面I下坐标的过程可以表示为式(1):
Figure BDA0004115837750000071
式中,K为相机的内参数矩阵,其中f表示相机焦距,dx和dy分别表示图像平面I上每个像素在相机的像平面C的x轴和y轴方向上的物理大小;u0和v0分别表示二维像平面坐标系O1-xy的原点O1在图像平面I上的坐标;(XC,YC,ZC)表示任意点P在相机坐标系O-XCYCZC下的坐标,u和v分别表示任意点P在图像平面I上的坐标;
通过上述标定方法对相机的焦距、像素等内参数进行标定,获取相机内参数矩阵,根据内参数矩阵,建立相机的计算机数字模型;
3、使用三维扫描仪对散斑纹理球形靶标进行三维扫描,以建立和散斑纹理球形靶标相同纹理的数字球形靶标;
4、通过相机对散斑纹理球形靶标进行成像,根据拍摄的散斑纹理球形靶标图像,确认数字球形靶标的位置;固定数字球形靶标的球心,按图3所示旋转数字球形靶标:
a、相机正对数字球形靶标的x轴方向;
b、将数字球形靶标沿着z轴顺时针根据实际需求转动一个角度α(α<2°),每转动一次拍摄一张数字球形靶标照片,并记录转动角度,直至旋转一周;
c、将数字球形靶标绕y轴顺时针转动一个角度α(α<2°),然后重复步骤b;
d、将数字球形靶标绕y轴逆时针转动一个角度α(α<2°),然后重复步骤b;
至此,便获得了数字球形靶标的整个正面纹理图像,将所有数字球形靶标图像与实际球形靶标图像一一作对比,计算每一次旋转的像素级汉明距离d,如式(2)所示。
Figure BDA0004115837750000081
式中i表示第i个像素,l表示像素的总数量,
Figure BDA0004115837750000082
表示实际球形靶标图像中像素i的灰度值,/>
Figure BDA0004115837750000083
表示数字球形靶标图像中像素i的灰度值;
然后从所有的d中求取出最小的汉明距离dmin,最小汉明距离dmin下对应的数字球形靶标图像与实际球形靶标图像是一致的,因此便获得了相机和散斑纹理球形靶标的位置关系,也即获得了待测物体的位姿。
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。除非另作定义,此处使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。为了避免混淆本发明的实质,公知的方法、过程和流程并没有详细叙述。
该散斑纹理球形靶标可放在水面浮体上,下面刚性连接水下仪器,可以实时观测水下仪器的位姿,如图4所示。
1、散斑纹理球形靶标的球体半径R为2m,球体表面全部覆盖如图1所示的随机散斑图案。
2、根据小孔成像模型,标定相机内参数,获得相机内参数矩阵K为:
Figure BDA0004115837750000091
3、求解汉明距离:通过matlab软件读取每组图像,然后求解像素级汉明距离,如图5和图6为比对过程中的实际球形靶标图像示意图和某一次旋转角度下的数字球形靶标图像,二者之间的汉明距离为191937;
4、确认水下仪器位姿
经过旋转步骤之后,得到一系列汉明距离,求解出最小的汉明距离。确认最小汉明距离对应的散斑纹理球形靶标位姿,如图7和图8所示,图7与图8之间的像素汉明距离为0。为了能直观地表示数字球形靶标和实际球形靶标图像一致性,将两张图片中圆形区域做减法,得到图9所示的圆形区域为纯白色的结果图(图像一致结果为0,0为白色)。根据最小汉明距离对应的散斑纹理球形靶标位姿确认水下仪器的位姿。

Claims (4)

1.一种散斑纹理球形靶标,其特征在于:该靶标是一个半径为R的三维球体,球体表面均匀覆盖黑白相间的随机散斑图案,且任意两个视角下球形靶标散斑图案均不重复。
2.一种散斑纹理球形靶标的位姿估计方法,其特征在于,该方法包括如下步骤:
步骤1、将散斑纹理球形靶标通过刚体连接方式固定到无法通过相机直接测量位姿的待测物体上;
在测量位姿之前,使用三维扫描仪对散斑纹理球形靶标与待测物体整体进行三维扫描,建立散斑纹理球形靶标与待测物体整体的计算机三维模型;
步骤2、建立相机小孔成像模型,根据相机小孔成像模型,散斑纹理球形靶标上任意点P从相机系下坐标转换为图像平面I下坐标的过程表示为式(1):
Figure FDA0004115837740000011
式中,K为相机的内参数矩阵,其中f表示相机焦距,dx和dy分别表示图像平面I上每个像素在相机的像平面C的x轴和y轴方向上的物理大小;u0和v0分别表示二维像平面坐标系O1-xy的原点O1在图像平面I上的坐标;(XC,YC,ZC)表示任意点P在相机坐标系O-XCYCZC下的坐标,u和v分别表示任意点P在图像平面I上的坐标;建立相机的计算机数字模型;
步骤3、使用三维扫描仪对散斑纹理球形靶标进行三维扫描,以建立与散斑纹理球形靶标相同纹理的数字球形靶标;
步骤4、通过相机对散斑纹理球形靶标进行成像,根据拍摄的靶标图像,确认数字球形靶标的位置,固定数字球形靶标的球心,按一定规则旋转数字球形靶标,每旋转一次靶标,通过模拟相机拍摄数字球形靶标,将数字球形靶标图像与实际球形靶标图像作对比,求取像素级汉明距离,在经过多次旋转之后,获得一系列汉明距离,从中求取出最小的汉明距离,最小汉明距离下对应的数字球形靶标图像与实际球形靶标图像是一致的,因此便获得了相机和球形靶标的位置关系,也即获得了待测物体的位姿。
3.根据权利要求2所述的散斑纹理球形靶标的位姿估计方法,其特征在于:在步骤2中,相机小孔成像模型的建立过程如下:
以相机光心O为原点,建立三维相机坐标系O-XCYCZC,在相机的像平面C中,以点O1为原点,建立二维像平面坐标系O1-xy;在图像平面I中以图像左上角Oi为原点,建立二维图像平面坐标系Oi-uv;点P(XC,YC,ZC)为散斑纹理球形靶标上任意一点,线段Pp经过相机光心O与相机的像平面C交于点p(xC,yC)。
4.根据权利要求2所述的散斑纹理球形靶标的位姿估计方法,其特征在于:步骤4包括如下过程:
通过相机对散斑纹理球形靶标进行成像,根据拍摄的靶标图像,确认数字球形靶标的位置;
固定数字球形靶标的球心,按如下规则旋转数字球形靶标:
a、相机正对数字球形靶标的x轴方向;
b、将数字球形靶标沿着z轴顺时针根据实际需求转动一个角度α,α<2°,每转动一次拍摄一张数字球形靶标照片,并记录转动角度,直至旋转一周;
c、将数字球形靶标绕y轴顺时针转动一个角度α,α<2°,然后重复步骤b;
d、将数字球形靶标绕y轴逆时针转动一个角度α,α<2°,然后重复步骤b;
至此,便获得了数字球形靶标的整个正面纹理图像,将所有数字球形靶标图像与实际球形靶标图像一一作对比,计算每一次旋转的像素级汉明距离d,图像对比之前,将RGB图像转为灰度图像,则只需要比较两幅图像的一个灰度通道,如式(2)所示;
Figure FDA0004115837740000031
式中i表示第i个像素,l表示像素的总数量,
Figure FDA0004115837740000033
表示实际球形靶标图像中像素i的灰度值,/>
Figure FDA0004115837740000032
表示数字球形靶标图像中像素i的灰度值;
然后从所有的d中求取出最小的汉明距离dmin,最小汉明距离dmin下对应的数字球形靶标图像与实际球形靶标图像是一致的,因此便获得了相机和散斑纹理球形靶标的位置关系,也即获得了待测物体的位姿。
CN202310218505.4A 2023-03-09 2023-03-09 一种散斑纹理球形靶标及其位姿估计方法 Pending CN116182783A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310218505.4A CN116182783A (zh) 2023-03-09 2023-03-09 一种散斑纹理球形靶标及其位姿估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310218505.4A CN116182783A (zh) 2023-03-09 2023-03-09 一种散斑纹理球形靶标及其位姿估计方法

Publications (1)

Publication Number Publication Date
CN116182783A true CN116182783A (zh) 2023-05-30

Family

ID=86442322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310218505.4A Pending CN116182783A (zh) 2023-03-09 2023-03-09 一种散斑纹理球形靶标及其位姿估计方法

Country Status (1)

Country Link
CN (1) CN116182783A (zh)

Similar Documents

Publication Publication Date Title
CN110555889B (zh) 一种基于CALTag和点云信息的深度相机手眼标定方法
CN109146980B (zh) 基于单目视觉的优化的深度提取和被动测距方法
CN109035320B (zh) 基于单目视觉的深度提取方法
CN110276808B (zh) 一种单相机结合二维码测量玻璃板不平度的方法
CN106548489B (zh) 一种深度图像与彩色图像的配准方法、三维图像采集装置
CN111272102A (zh) 一种线激光扫描三维测量标定方法
CN113012234B (zh) 基于平面变换的高精度相机标定方法
CN111383194B (zh) 一种基于极坐标的相机畸变图像校正方法
CN112132907B (zh) 一种相机标定方法、装置、电子设备及存储介质
JP2000227309A (ja) 3次元位置姿勢センシング装置
CN111243033A (zh) 一种优化双目相机外参数的方法
CN112614188B (zh) 一种基于交比不变性的点阵式标定板及其识别方法
CN114998448B (zh) 一种多约束双目鱼眼相机标定与空间点定位的方法
CN110763204A (zh) 一种平面编码靶标及其位姿测量方法
CN116912333B (zh) 一种基于作业围栏标定杆的相机姿态自标定方法
CN111699513B (zh) 一种标定板、内参数标定方法、机器视觉系统及存储装置
CN109859137A (zh) 一种广角相机非规则畸变全域校正方法
CN107680035B (zh) 一种参数标定方法和装置、服务器及可读存储介质
CN114359405A (zh) 一种离轴沙姆3d线激光相机的标定方法
CN113329179A (zh) 拍摄对位方法、装置、设备及存储介质
CN115147499A (zh) 标定参数确定方法、混合标定板、装置、设备和介质
CN112446926B (zh) 一种激光雷达与多目鱼眼相机的相对位置标定方法及装置
CN116182783A (zh) 一种散斑纹理球形靶标及其位姿估计方法
CN113483669B (zh) 一种基于立体靶标的多传感器位姿标定方法及装置
CN115937326A (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