CN112097689B - 一种3d结构光系统的标定方法 - Google Patents
一种3d结构光系统的标定方法 Download PDFInfo
- Publication number
- CN112097689B CN112097689B CN202010951187.9A CN202010951187A CN112097689B CN 112097689 B CN112097689 B CN 112097689B CN 202010951187 A CN202010951187 A CN 202010951187A CN 112097689 B CN112097689 B CN 112097689B
- Authority
- CN
- China
- Prior art keywords
- image
- projector
- coordinates
- point
- camera
- 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
Links
Images
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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring 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/254—Projection of a pattern, viewing through a pattern, e.g. moiré
-
- 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring 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/2536—Measuring 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 using several gratings with variable grating pitch, projected on the object with the same angle of incidence
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种3D结构光系统的标定方法,属于计算机视觉技术领域。本发明提供了一种3D结构光系统的标定方法,该方法使用精度更高,标定流程更复杂的包括圆点图案的标定板,使用精度更高的多频相移码来寻找相机和投影传感器的对应关系,局部提升相位精度,能够更精确确定传感器之间位置关系,同时采用灵活的排序算法对圆心位置进行排序,标定板基本可以摆放为各种姿态,提取更多的标定数据。
Description
技术领域
本发明涉及计算机视觉技术领域,尤其涉及一种3D结构光系统的标定方法。
背景技术
目前,3D结构光成像技术在工业上主要应用在精密测量,缺陷检测等场景,涵盖包含表面组装技术、汽车、航空、半导体、医疗、制药、食品加工等生产制造领域。和二维检测相比,3D检测技术提供了多一维的高度信息,在实际测量过程中,不仅可以精确定位被测物,精确测量尺寸,颜色等二维信息,还能够获取物体表面轮廓信息,这样就能测量分析物体表面的平面度、坡度、弯曲度和瑕疵等复杂信息。3D结构光测量具有非接触、速度快、精度高、抗干扰能力强等一系列优点,是3D检测行业精密测量的一个非常重要的方向。
3D结构光标定算法,决定了3D测量系统的检测精度,该算法主要是确立相机、结构光投影的内部参数以及他们构成系统的外部参数,正是通过确定这些系统模型参数才能恢复出被测物的3D深度信息。
目前结构光3D系统标定原理都采用张正友标定法,但是在具体实施的时候,方法却不尽相同,当然,精度差异也非常大,提供精度较高的标定方法逐渐成为一个研究的领域。
在中国专利申请文献CN110443856A,公开了一种3D结构光模组标定方法,该方法包括如下步骤:S0、制作标定板,在标定板上打印棋盘格,所述棋盘格铺满整个标定板;S1、拍摄标定照片,放置3D结构光模组于所述标定板平面的正前方,固定所述3D结构光模组和所述标定板的相对位置,依次拍摄所述标定板的红外照片、点阵投射到所述标定板的点阵投射照片,将所述标定板的红外照片和所述点阵投射照片作为一组标定照片;S2、拍摄标定照片组,更换所述3D结构光模组和所述标定板的相对位置,重复步骤S1拍摄多组标定照片;S3、标定3D结构光模组,采用标定算法通过多组所述标定照片对红外相机进行标定,然后对点阵投射器进行标定,得到所述点阵投射器相对于所述红外相机的几何关系。所述标定板为白色背景的标定板,所述棋盘格的形状为正方形,所述棋盘格中黑色格子的灰度可调;所述标定算法为张正友平面标定法;标定过程包括:S31、提取像素坐标,提取所述标定板的红外照片的棋盘格点像素坐标,并设定对应的世界坐标数值;S32、标定红外相机,计算棋盘格点世界坐标到像素坐标的单应性映射,根据旋转矩阵的约束关系计算得到对应的红外相机内参和外参;S33、计算转换关系,计算点阵投射的光斑坐标与所述点阵投射器的点阵几何关系之间的转换关系,将红外照片的棋盘格点像素坐标转换成所述点阵投射器的像素坐标,采用张正友平面标定算法计算所述点阵投射器的内参和外参。该方案采用棋盘格标定板,需要检测标定板的角点,作为Mark点,在实际标定应用中,标定板位姿变化摆放时,受光照以及成像质量不均匀等影响,棋盘格标定板的角点提取精度很容易受影响,采用圆的标定板,可以很好的进行圆的提取,不易受干扰。
在中国专利申请文献CN110827357A,公开了一种组合图案标定板以及结构光相机参数标定方法,标定板上的图案包括识别标志和码块区域;所述识别标志为黑白边界区分的闭合的四边形,其边界内为黑色像素,边界外部为白色像素;所述码块区域位于四边形边界内部,码块由NxN的黑色和白色方形色块组成的编码信息。标定方法包括如下步骤:S1:利用权利要求1中的多个组合图案标定板,各组合图案标定板具有图案不同的码块区域;S2:采集分布排列的多个组合图案标定板的图像,并进行预处理;S3:计算图像每个像素的局部梯度;将梯度辐值和方向一致的像素区域设置为连通域;S4:将连通域中梯度大于设定阈值的连通域作为标定图案的边缘,边缘像素点使用直线拟合成不同的线段;S5:遍历所有的线段,检测四个依次邻接的线段能否组成完整的四边形,并检测所有具有完整四边形的区域是否包含有效的编码区域;S6:将检测到的编码与编码库的图案进行匹配,计算编码对应的图案在标定板中的坐标,以及四边形四条线段的交点在图像中对应的坐标;S7:根据图像中的像素坐标和对应的标定板坐标,对结构光相机参数进行标定。该方案中标定板位姿不易倾斜过大,过大会造成图案匹配不上,只适宜小幅度倾斜,这样标定出来的参数就不能充分表述系统的模型。
现有技术至少存在以下不足:
1.确立相机传感器和投影仪传感器对应关系精度不够。
2.所获取标定图像状态有限,无法充分描述系统模型。
发明内容
为解决现有技术中存在的精度不够问题,本发明提供了一种3D结构光系统的标定方法,该方法使用精度更高,标定流程更复杂的包括圆点图案的标定板,使用精度更高的多频相移码来寻找相机和投影传感器的对应关系,局部提升相位精度,能够更精确确定传感器之间位置关系,同时采用灵活的自动梯形校正算法对圆点进行排序,标定板基本可以摆放为各种姿态,提取更多的标定数据。
本发明提供了一种3D结构光系统的标定方法,包括如下步骤:
制作标定板步骤,制作系统标定板并烧录,所述标定板上的图案为多个圆形靶标,相邻的所述圆形靶标的圆心距一致,多个所述圆形靶标排列为行和列,所述圆形靶标的大小不同;
标定板图像采集步骤:
设定一个标定板姿态,分别在标定板上投影水平和竖直两个方向的多个频率的正弦光栅条纹图案,另外投影一张没有光栅条纹的光,并在每次投影的同时用相机采集该姿态下投影各种图案的标定板图像;
调整标定板姿态,重复标定板图像采集步骤,直至标定板姿态数量达到预设的目标值,得到多个标定板图像;
标定步骤,根据采集到的多个所述标定板图像,对结构光系统进行标定,具体包括如下步骤:
提取标定板图像中所有圆形靶标的圆心位置;
根据提取的所有圆形靶标的圆心位置对所有圆心进行排序;
采用双线性插值算法得到每个圆心坐标的绝对相位值;
利用得到的每个圆心坐标的绝对相位值,得到各圆心在相机图像中坐标对应的在投影仪图像中坐标;
根据得到的相机及投影仪图像中的圆心坐标,利用张正友标定方法得到相机及投影仪的内参、畸变系数和外参;
利用同一标定板姿态下的相机和投影仪的外参,计算每个标定板姿态下相机和投影仪坐标系之间的转换矩阵;
根据得到的每个标定板姿态下的转换矩阵,得到标定的相机和投影仪坐标系转换矩阵。
优选地,还包括,在采用双线性插值算法得到每个圆心坐标的绝对相位值之前,利用多频多步相移算法,得到标定板平面相位场在水平和垂直方向的绝对相位值。
优选地,所述提取标定板图像中所有圆形靶标的圆心位置具体包括以下步骤:
对相机采集到的没有打光栅条纹的标定板图像进行二值化;
对二值化后的标定板图像进行blob分析,得到标定板圆的质心位置,实现圆形靶标的粗定位;
通过Canny边缘提取,得到圆形靶标的像素边缘位置;
通过对圆形靶标的边缘像素点进行灰度插值,然后再进行参数拟合,得到圆的亚像素边缘轮廓;
对得到的各个圆的亚像素边缘轮廓点进行最小二乘椭圆拟合,得到各椭圆的中心位置,将各椭圆的中心位置作为各圆形靶标的圆心位置。
优选地,对所有圆心进行排序具体包括以下步骤:
根据得到的标定板各椭圆的中心位置,确定采集的标定板图案中不规则四边形的四个顶点;
将四个顶点按顺序连线成四边形;
当所述四边形旋转角度不超过预设值时,确定各圆点的位置坐标,并根据各圆点的位置坐标顺序,确定各圆心的顺序;
当所述四边形旋转角度超过预设值时,按照如下方法进行排序:
确定所述四边形的每条边上得所有圆点,将两条相对的边上相对的两个圆点进行连线;
确定所有连线形成的各个四边形内部所有的圆点,并确定各圆点的顺序。
优选地,当所述四边形旋转角度不超过预设值时,确定各圆心顺序具体包括以下步骤:
确定标定板图像中五大圆形靶标的位置坐标;
由这些已知阵列坐标的大圆点,根据摄影几何的直线投影不变性,对位于阵列其他位置的小圆点在图像上的坐标进行预估;
通过blob分析得到的各圆形靶标的坐标,计算各点到预估的所有点的距离,距离最小的点为对应的位置的点,根据坐标顺序进行排序。
优选地,确定五大圆形靶标的位置坐标具体包括如下步骤:
选择五个半径相等的圆形靶标,所述五个半径相等的圆形靶标的半径大于标定板图案中其它圆形靶标的半径;
计算各点之间的距离,得到距离最大的两个圆点,分别为F1和F2,以及距离最小的两个圆点,分别为N1和N2,则剩下的一点确定为五大圆形靶标中第五大圆点C5;
分别计算F1、F2到N1、N2的距离之和;
距离和较小的点为五大圆形靶标中第一大圆点C1,所述第一大圆点位于图案水平中线与垂直中线左侧的垂直中线交叉点;
距离和较大的点为五大圆形靶标中第二大圆点C2,所述第二大圆点位于图案水平中线与垂直中线右侧的垂直中线交叉点;
计算N1、N2与C1的距离,距离较小的点为五大圆形靶标中第三大圆点C3,较大的点为五大圆形靶标中第四大圆点C4。
优选地,预估小圆点坐标具体包括如下步骤:
用两个已知的大圆点C1和C2确定一条水平直线,计算与该直线平行的第二行两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板水平方向的两个分量;
用两个已知的大圆点C3和C4确定一条垂直直线,与该直线平行的第二列两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板垂直方向的两个分量;
根据已知大圆点C1的坐标,以及标定板水平和垂直方向各自的两个分量,预估得到所有圆形靶标的坐标,并存储所有圆形靶标位置坐标。
优选地,采用双线性插值算法得到每个圆心坐标在水平和垂直方向的绝对相位值,具体包括如下步骤:
根据提取得到的各圆心亚像素位置坐标,利用双线性插值得到该圆心对应的相位主值,插值公式如下:
φ(i+u,j+v)=(1-u)(1-v)φ(i,j)+(1-u)vφ(i,j+1)+u(1-v)φ(i+1,j)+uvφ(i+1,j+1);
其中,
φ为相位主值;
i,j分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的整数部分;
u,v分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的小数部分;
相位主值通过如下方法获得:
在采集的投影了正弦条纹码的图像中,多步相移法光栅光强分布函数为:
Ik(x,y)=I’(x,y)+I”(x,y)cos(φ+(k-1)π/2)
其中,
k取值为1-n,其中n为多步相移法中的步数;
x,y为像素坐标;
I(x,y)为第k幅相移图像中(x,y)位置的像素点的灰度值;
I'(x,y)为图像的平均灰度;
I”(x,y)为图像的灰度调制;
φ为相位主值;
通过光强分布函数,得到每个像素点对应的相位主值:
其中,
k取值为1-n,其中n为多步相移法中的步数;
x,y为像素坐标;
Ik为该点像素灰度值;
N为光栅条纹周期个数;
φ(x,y)为相位主值。
优选地,利用得到的每个圆心坐标在水平和垂直方向的绝对相位值,通过如下公式得到各圆心在相机图像中坐标对应的在投影仪图像中坐标:
其中,
up为圆心c点在投影仪图像中u方向的坐标;
vp为圆心c点在投影仪图像中v方向的坐标;
N为光栅条纹周期个数;
W为投影仪水平方向的分辨率;
H为投影仪垂直方向的分辨率;
Φu(uc,vc)为圆心c点垂直方向绝对相位值;
Φv(uc,vc)为圆心c点水平方向绝对相位值;
优选地,还包括,根据标定得到的相机和投影仪内参、畸变系数和外参,以及相机与投影仪之间的坐标系转换矩阵,建立相机和投影仪与世界坐标系的映射关系,具体包括如下步骤:
通过多频外差法,得到绝对相位灰度图;
对于空间任一点p[xw,yw,zw],其投影在相机中的图像坐标为p(uc,vc),根据得到的拍摄图像上的每个像素点的绝对相位值,得到对应的投影仪图像中的一条直线,相机图像与投影仪图像的对应关系为:
其中:
N为光栅条纹周期个数;
W为投影仪水平方向的分辨率;
Φ(uc,vc)为该像素点的绝对相位值;
up为该像素点对应的投影仪图像中垂直方向的坐标。
根据上述相机图像与投影仪图像的对应关系,由如下成像原理公式,得到唯一的点p的三维坐标(XW,YW,ZW):
SC[uc,vc,1]=ACMC[XW,YW,ZW,1]
SP[up,vp,1]=APMP[XW,YW,ZW,1]
其中,
AC和AP分别为相机和投影仪的内参;
MC和MP分别为相机和投影仪的外参;
SC和SP分别为相机和投影仪的比例因子;
(uc,vc)和(up,vp)分别为相机和投影仪的图像坐标,两者均使用预先标定出的系统畸变参数进行畸变矫正;
(XW,YW,ZW)为点p的唯一三维坐标。
与现有技术相对比,本发明的有益效果如下:
(1)本发明采用圆靶标标定板,通过提取圆的亚像素轮廓进行椭圆拟合,得到圆心的亚像素位置,精度更高,同时标定过程中,标定板随意摆放,造成照明亮度不一致,圆的目标特征比较明显,也能够适应光照变化,准确提取圆心位置;
(2)本发明利用双线性插值得到每个圆心坐标所对应的绝对相位值,然后利用绝对相位值计算出对应的投影仪图像坐标。
附图说明
图1是本发明优选实施方式的标定方法流程图;
图2是本发明优选实施方式的标定步骤流程图;
图3是本发明优选实施方式的带角度的标定图像中圆点图案的排序示意图。
具体实施方式
下面结合附图1-3,对本发明的具体实施方式作详细的说明。
本发明提供了一种3D结构光系统的标定方法,包括如下步骤:
制作标定板步骤,制作系统标定板并烧录,所述标定板上的图案为多个圆形靶标,相邻的所述圆形靶标的圆心距一致,多个所述圆形靶标排列为行和列,所述圆形靶标的大小不同;
标定板图像采集步骤,
设定一个标定板姿态,分别在标定板上投影水平和竖直两个方向的多个频率的正弦光栅条纹图案,另外投影一张没有光栅条纹的光,并在每次投影的同时用相机采集该姿态下投影各种图案的标定板图像;
调整标定板姿态,重复标定板图像采集步骤,直至标定板姿态数量达到预设的目标值,得到多个标定板图像;
标定步骤,根据采集到的多个所述标定板图像,对结构光系统进行标定,具体包括如下步骤:
提取标定板图像中所有圆形靶标的圆心位置;
根据提取的所有圆形靶标的圆心位置对所有圆心进行排序;
采用双线性插值算法得到每个圆心坐标的绝对相位值;
利用得到的每个圆心坐标的绝对相位值,得到各圆心在相机图像中坐标对应的在投影仪图像中坐标;
根据得到的相机及投影仪图像中的圆心坐标,利用张正友标定方法得到相机及投影仪的内参、畸变系数和外参;
利用同一标定板姿态下的相机和投影仪的外参,计算每个标定板姿态下相机和投影仪坐标系之间的转换矩阵;
根据得到的每个标定板姿态下的转换矩阵,得到标定的相机和投影仪坐标系转换矩阵。
将得到的所有转换矩阵中对应位置的参数叠加求和并取平均,作为最终的投影仪和相机坐标系的转换矩阵。
标定得到的参数包括:相机内参和畸变系数,投影仪内参和畸变系数以及相机和投影仪坐标系之间的刚性变换关系矩阵。
作为优选实施方式,还包括,在采用双线性插值算法得到每个圆心坐标的绝对相位值之前,利用多频多步相移算法,得到标定板平面相位场在水平和垂直方向的绝对相位值。具体可采用十步相移码和三频外差法,得到水平和垂直方向的绝对相位值,包括如下步骤:
在标定图像采集步骤中,对标定板分别投影的水平和竖直方向不同频率的正弦光栅条纹图案,所述正弦光栅条纹图案采用三种不同波长,波长分别为:λ1,λ2,λ3,其对应相位主值分别为:φ1,φ2,φ3;每个投影图案采集十张图像,用于十步相移码计算;
采用三频外差法计算绝对相位值:
将第一波长λ1的第一相位φ1与第二波长λ2的第二相位φ2进行叠加,得到相位φ12,其中相位φ12所对应的波长为λ12,叠加计算得到的λ12为:
将第二波长λ2的第二相位φ2与第三波长λ3的第三相位φ3进行叠加,得到相位φ23,其中相位φ23所对应的波长为λ23,叠加计算得到的λ23为:
将λ12和λ23进行叠加,得到最终叠加的波长λ123。
通过以下公式求取叠加的绝对相位值:
式中,
ni为被测物体表面某点在对应光栅图中的条纹级数,ni包含整数部分Ni和小数部分Δni;
φi为对应光栅的包裹相位;
Φi为对应光栅的绝对相位;
由光栅12和光栅23叠加产生光栅123,选择λ1,λ2,λ3,使得光栅123的波长λ123覆盖全场,使N123=0,并由此得到光栅123的绝对相位。
作为优选实施方式,所述提取标定板图像中所有圆形靶标的圆心位置具体包括以下步骤:
对相机采集到的没有打光栅条纹的标定板图像进行二值化;
对二值化后的标定板图像进行blob分析,得到标定板圆的质心位置,实现圆形靶标的粗定位;
通过Canny边缘提取,得到圆形靶标的像素边缘位置;
通过对圆形靶标的边缘像素点进行灰度插值,然后再进行参数拟合,得到圆的亚像素边缘轮廓;
对得到的各个圆的亚像素边缘轮廓点进行最小二乘椭圆拟合,得到各椭圆的中心位置,将各椭圆的中心位置作为各圆形靶标的圆心位置。
作为优选实施方式,对所有圆心进行排序具体包括以下步骤:
根据得到的标定板各椭圆的中心位置,确定采集的标定板图案中不规则四边形的四个顶点;
将四个顶点按顺序连线成四边形;
当所述四边形旋转角度不超过预设值时,确定各圆点的位置坐标,并根据各圆点的位置坐标顺序,确定各圆心的顺序;
当所述四边形旋转角度超过预设值时,按照如下方法进行排序:
确定所述四边形的每条边上得所有圆点,将两条相对的边上相对的两个圆点进行连线;
确定所有连线形成的各个四边形内部所有的圆点,并确定各圆点的顺序。
作为优选实施方式,当所述四边形旋转角度不超过预设值时,确定各圆心顺序具体包括以下步骤:
确定标定板图像中五大圆形靶标的位置坐标;
由这些已知阵列坐标的大圆点,根据摄影几何的直线投影不变性,对位于阵列其他位置的小圆点在图像上的坐标进行预估;
通过blob分析得到的各圆形靶标的坐标,计算各点到预估的所有点的距离,距离最小的点为对应的位置的点,根据坐标顺序进行排序。
作为优选实施方式,确定五大圆形靶标的位置坐标具体包括如下步骤:
选择五个半径相等的圆形靶标,所述五个半径相等的圆形靶标的半径大于标定板图案中其它圆形靶标的半径;
计算各点之间的距离,得到距离最大的两个圆点,分别为F1和F2,以及距离最小的两个圆点,分别为N1和N2,则剩下的一点确定为五大圆形靶标中第五大圆点C5;
分别计算F1、F2到N1、N2的距离之和;
距离和较小的点为五大圆形靶标中第一大圆点C1,所述第一大圆点位于图案水平中线与垂直中线左侧的垂直中线交叉点;
距离和较大的点为五大圆形靶标中第二大圆点C2,所述第二大圆点位于图案水平中线与垂直中线右侧的垂直中线交叉点;
计算N1、N2与C1的距离,距离较小的点为五大圆形靶标中第三大圆点C3,较大的点为五大圆形靶标中第四大圆点C4。
在大小上,五个大圆半径相等,且稍大于标定板其余圆的半径;在排布上,距离最大的两个圆与标定板图案中心圆属同一行,分别位于图案中间圆的左侧和右侧,与中心圆均相隔两个小圆,距离最小的两个圆在同一行且位置相邻,其中一圆与图案中心圆属同一列,且相隔一小圆,当该两圆位于图案水平中线下方时,其中一大圆在垂直中心线的右侧,最后一大圆同样与图案中心圆属于同一列,在图案水平中心线的另一侧,且同样相隔一小圆。
在标定过程中,标定板的位姿是随意摆放的,确定五个大圆的位置,可以得到沿图案水平方向和垂直方向的方向向量,从而得到图案中所有点的大概位置,通过blob分析和椭圆拟合提取的真实圆心,并和预估的圆心做圆心距比较,距离最小的圆心便是对应该位置的圆。
作为优选实施方式,预估小圆点坐标具体包括如下步骤:
用两个已知的大圆点C1和C2确定一条水平直线,计算与该直线平行的第二行两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板水平方向的两个分量;
用两个已知的大圆点C3和C4确定一条垂直直线,与该直线平行的第二列两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板垂直方向的两个分量;
根据已知大圆点C1的坐标,以及标定板水平和垂直方向各自的两个分量,预估得到所有圆形靶标的坐标,并存储所有圆形靶标位置坐标。
作为优选实施方式,采用双线性插值算法得到每个圆心坐标在水平和垂直方向的绝对相位值,具体包括如下步骤:
根据提取得到的各圆心亚像素位置坐标,利用双线性插值得到该圆心对应的相位主值,插值公式如下:
φ(i+u,j+v)=(1-u)(1-v)φ(i,j)+(1-u)vφ(i,j+1)+u(1-v)φ(i+1,j)+uvφ(i+1,j+1);
其中,
φ为相位主值;
i,j分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的整数部分;
u,v分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的小数部分;
相位主值通过如下方法获得:
在采集的投影了正弦条纹码的图像中,多步相移法光栅光强分布函数为:
Ik(x,y)=I’(x,y)+I”(x,y)cos(φ+(k-1)π/2)
其中,
k取值为1-n,其中n为多步相移法中的步数;
x,y为像素坐标;
I(x,y)为第k幅相移图像中(x,y)位置的像素点的灰度值;
I'(x,y)为图像的平均灰度;
I”(x,y)为图像的灰度调制;
φ为相位主值;
通过光强分布函数,得到每个像素点对应的相位主值:
其中,
k取值为1-n,其中n为多步相移法中的步数;
x,y为像素坐标;
Ik为该点像素灰度值;
N为光栅条纹周期个数;
φ(x,y)为相位主值。
当采用十步相移法时,在采集的投影了正弦条纹码的图像中,十步相移法光栅光强分布函数为:
Ik(x,y)=I’(x,y)+I”(x,y)cos(φ+(k-1)π/2)
其中,
k取值为1-10;
x,y为像素坐标;
I(x,y)为第k幅相移图像中(x,y)位置的像素点的灰度值;
I'(x,y)为图像的平均灰度;
I”(x,y)为图像的灰度调制;
φ为相位主值;
通过光强分布函数,得到每个像素点对应的相位主值:
其中,
k取值为1-10;
x,y为像素坐标;
Ik为该点像素灰度值;
N为光栅条纹周期个数;
φ(x,y)为相位主值。
通过上面方法获得的相位主值对应的像素点是整数坐标,根据提取得到的各圆心亚像素位置坐标,利用双线性插值得到该圆心对应的相位主值,插值公式如下:
φ(i+u,j+v)=(1-u)(1-v)φ(i,j)+(1-u)vφ(i,j+1)+u(1-v)φ(i+1,j)+uvφ(i+1,j+1)
其中,
φ为相位主值;
i,j分别为用于插值的相机图像中的点在垂直和水平方向的坐标;
u,v分别为用于插值的投影仪图像中的点在垂直和水平方向的坐标;
作为优选实施方式,利用得到的每个圆心坐标在水平和垂直方向的绝对相位值,通过如下公式得到各圆心在相机图像中坐标对应的在投影仪图像中坐标:
其中,
up为圆心c点在投影仪图像中u方向的坐标;
vp为圆心c点在投影仪图像中v方向的坐标;
N为光栅条纹周期个数;
W为投影仪水平方向的分辨率;
H为投影仪垂直方向的分辨率;
Φu(uc,vc)为圆心c点垂直方向绝对相位值;
Φv(uc,vc)为圆心c点水平方向绝对相位值;
作为优选实施方式,还包括,根据标定得到的相机和投影仪内参、畸变系数和外参,以及相机与投影仪之间的坐标系转换矩阵,建立相机和投影仪与世界坐标系的映射关系,具体包括如下步骤:
通过多频外差法,得到绝对相位灰度图;此处可采用三频外差法;
对于空间任一点p[xw,yw,zw],其投影在相机中的图像坐标为p(uc,vc),根据得到的拍摄图像上的每个像素点的绝对相位值,得到对应的投影仪图像中的一条直线,相机图像与投影仪图像的对应关系为:
其中:
N为光栅条纹周期个数;
W为投影仪水平方向的分辨率;
Φ(uc,vc)为该像素点的绝对相位值;
up为该像素点对应的投影仪图像中垂直方向的坐标。
根据上述相机图像与投影仪图像的对应关系,由如下成像原理公式,得到唯一的点p的三维坐标(XW,YW,ZW):
SC[uc,vc,1]=ACMC[XW,YW,ZW,1]
SP[up,vp,1]=APMP[XW,YW,ZW,1]
其中,
AC和AP分别为相机和投影仪的内参;
MC和MP分别为相机和投影仪的外参;
SC和SP分别为相机和投影仪的比例因子;
(uc,vc)和(up,vp)分别为相机和投影仪的图像坐标,两者均使用预先标定出的系统畸变参数进行畸变矫正;
(XW,YW,ZW)为点p的唯一三维坐标。
由相机和投影仪图像的对应关系表达式知,在投射水平或垂直单向条纹的情况下,相机中一像素坐标点对应的投影仪图像中一水平或垂直方向的坐标。由垂直的绝对相位值可以在投影仪图像中确定一条垂直的对应线,即确定投影仪图像对应点在垂直方向的坐标值,作为该点列坐标;同理水平的绝对相位值可以确定一条水平的对应线,即确定投影仪图像对应点在水平方向的坐标值,作为该点行坐标。
实施例1
根据本发明的一个具体实施方案,下面对本发明中的圆心排序方法进行详细说明。标定板上的圆形靶标图案如图3所示。圆点排序包括如下步骤:
首先确定中间五大圆形靶标的位置坐标:
1.计算各点之间的距离,求出距离最大的两点(设为F1和F2)和最小的两点(设为N1和N2),则剩下的一点即可确定为点第72点(如图3所示);
2.分别计算F1、F2到N1、N2的距离之和,则和较小的点为第47点,较大的点为第53点。
3.计算N1、N2和第47点的距离,距离较小的点为第27点,较大的点为第28点。
然后,由这些已知阵列坐标的大圆点,根据摄影几何的直线投影不变性,对位于阵列其他位置的小圆点在图像上的点进行预估。主要步骤如下:
1.用第47点和第53点两个已知的点确定一条直线,计算第二行两个已知点之间的距离,以及两个点连线的X和Y方向分量,从而得到标定板水平方向的两个分量;
2.同理,由第28点和第72点两点确定垂直方向的两个分量;
3.由已知的第47点的坐标,以及标定板水平和垂直方向的各自的两个分量,预判得到第1点的坐标,从而由第1点的坐标,以及标定板水平和垂直方向的各自的两个分量,预判得到所有99点的坐标(之前已优先提取得到大圆的坐标,可以不用预判,直接存储其坐标位置);
4.通过blob分析得到的所有点的质心坐标,计算各点到预判的99点的距离,距离最小的点为对应的位置的点,从而实现排序,得到如图的排序效果。
实施例2
根据本发明的一个具体实施方案,下面详细说明本发明中标定步骤中的具体实施。
1.通过对不打光栅条纹的图像处理,得到每个圆圆心的亚像素坐标;
2.对于水平和垂直的光栅图像,通过十步相移和多频外差法得到每个像素点的水平和垂直方向的相位主值;
3.双线性插值分别得到该点的水平和垂直方向的绝对相位值;
4.通过相机图像和投影仪图像的对应关系公式得到该点的投影仪图像坐标;
5.已知摄像机的针孔模型:
在标定中,空间三维世界坐标为(XW,YW,ZW),也就是标定板的物理坐标,(u,v)为提取的圆心的相机图像坐标或投影仪图像坐标,mij为对应性矩阵的元素。
标定板图案中,圆心距一致,可以建立平面坐标系,通过得到的每一组相机图像的圆心坐标位置和投影仪图像的圆心坐标位置,以及标定板的物理坐标,可以通过张正友标定算法分别得到相机和投影仪的内参,畸变系数,同时,对于每一个位姿下的标定板,张正友标定法可以分别得到对应的相机和投影仪一组外参,利用内参,畸变系数和一组外参可以确定一组单应性矩阵,即上述M矩阵。
利用得到的每一个姿态下的相机的外参矩阵,乘以同一位姿下投影仪外参矩阵的逆,可以得到一新的矩阵,该矩阵就是相机和投影仪的关系矩阵,通过求取每一姿态下的关系矩阵,然后将得到的关系矩阵中的(i,j)位置的元素加和取平均,作为相机和投影仪的转换矩阵。
6.以图3中的点47为例,第5步中的针孔模型包含三个方程,整理消去ZC,可以得到两个关于mij的线性方程,即:
m11XW+m12YW+m13ZW+m14-uXWm31-uYWm32-uZWm33=um34
m21XW+m22YW+m23ZW+m24-vXWm31-vYWm32-vZWm33=vm34
通过张正友标定法分别得到的相机和投影仪单应性矩阵,从而可以得到四个线性方程组,由于只有XW,YW,ZW三个变量,从而可以联立方程得到该点空间三维坐标。
因为外参和转换关系矩阵只包含旋转和平移的刚性关系,所以不会改变物体形状和尺寸。通过外参矩阵,可将空间三维坐标转换至相机坐标系下,利用相机和投影仪的转换关系矩阵就可以得到相机坐标系一点在投影仪下的位置坐标。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均包含在本发明的保护范围之内。
Claims (7)
1.一种3D结构光系统的标定方法,其特征在于,包括如下步骤:
制作标定板步骤,制作系统标定板并烧录,所述标定板上的图案为多个圆形靶标,相邻的所述圆形靶标的圆心距一致,多个所述圆形靶标排列为行和列,所述圆形靶标的大小不同;
标定板图像采集步骤:
设定一个标定板姿态,分别在标定板上投影水平和竖直两个方向的多个频率的正弦光栅条纹图案,另外投影一张没有光栅条纹的图案,并在每次投影的同时用相机采集该姿态下投影各种图案的标定板图像;
调整标定板姿态,重复标定板图像采集步骤,直至标定板姿态数量达到预设的目标值,得到多个标定板图像;
标定步骤,根据采集到的多个所述标定板图像,对结构光系统进行标定,具体包括如下步骤:
提取标定板图像中所有圆形靶标的圆心位置,具体包括以下步骤:
对相机采集到的没有打光栅条纹的标定板图像进行二值化;
对二值化后的标定板图像进行blob分析,得到标定板圆的质心位置,实现圆形靶标的粗定位;
通过Canny边缘提取,得到圆形靶标的像素边缘位置;
通过对圆形靶标的边缘像素点进行灰度插值,然后再进行参数拟合,得到圆的亚像素边缘轮廓;
对得到的各个圆的亚像素边缘轮廓点进行最小二乘椭圆拟合,得到各椭圆的中心位置,将各椭圆的中心位置作为各圆形靶标的圆心位置;
根据提取的所有圆形靶标的圆心位置对所有圆心进行排序,具体包括以下步骤:
根据得到的标定板各椭圆的中心位置,确定采集的标定板图案中不规则四边形的四个顶点;
将四个顶点按顺序连线成四边形;
当所述四边形旋转角度不超过预设值时,确定各圆点的位置坐标,并根据各圆点的位置坐标顺序,确定各圆心的顺序;
当所述四边形旋转角度超过预设值时,按照如下方法进行排序:
确定所述四边形的每条边上的 所有圆点,将两条相对的边上相对的两个圆点进行连线;
确定所有连线形成的各个四边形内部所有的圆点,并确定各圆点的顺序;
利用十步相移和多频外差法,得到标定板平面相位场在水平和垂直方向的绝对相位值;
采用双线性插值算法得到每个圆心坐标的绝对相位值;
利用得到的每个圆心坐标的绝对相位值,得到各圆心在相机图像中坐标对应的在投影仪图像中坐标;
根据得到的相机及投影仪图像中的圆心坐标,利用张正友标定方法得到相机及投影仪的内参、畸变系数和外参;
利用同一标定板姿态下的相机和投影仪的外参,计算每个标定板姿态下相机和投影仪坐标系之间的转换矩阵;
根据得到的每个标定板姿态下的转换矩阵,得到标定的相机和投影仪坐标系转换矩阵。
2.根据权利要求1所述的3D结构光系统的标定方法,其特征在于,当所述四边形旋转角度不超过预设值时,确定各圆心顺序具体包括以下步骤:
确定标定板图像中五大圆形靶标的位置坐标;
由这些已知阵列坐标的大圆点,根据摄影几何的直线投影不变性,对位于阵列其他位置的小圆点在图像上的坐标进行预估;
通过blob分析得到的各圆形靶标的坐标,计算各点到预估的所有点的距离,距离最小的点为对应的位置的点,根据坐标顺序进行排序。
3.根据权利要求2所述的3D结构光系统的标定方法,其特征在于,确定五大圆形靶标的位置坐标具体包括如下步骤:
选择五个半径相等的圆形靶标,所述五个半径相等的圆形靶标的半径大于标定板图案中其它圆形靶标的半径;
计算各点之间的距离,得到距离最大的两个圆点,分别为F1和F2,以及距离最小的两个圆点,分别为N1和N2,则剩下的一点确定为五大圆形靶标中第五大圆点C5;
分别计算F1、F2到N1、N2的距离之和;
距离和较小的点为五大圆形靶标中第一大圆点C1,所述第一大圆点位于图案水平中线与垂直中线左侧的垂直中线交叉点;
距离和较大的点为五大圆形靶标中第二大圆点C2,所述第二大圆点位于图案水平中线与垂直中线右侧的垂直中线交叉点;
计算N1、N2与C1的距离,距离较小的点为五大圆形靶标中第三大圆点C3,较大的点为五大圆形靶标中第四大圆点C4。
4.根据权利要求2所述的3D结构光系统的标定方法,其特征在于,预估小圆点坐标具体包括如下步骤:
用两个已知的大圆点C1和C2确定一条水平直线,计算与该直线平行的第二行两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板水平方向的两个分量;
用两个已知的大圆点C3和C4确定一条垂直直线,与该直线平行的第二列两个已知点之间的距离,以及所述两个已知点连线的X和Y方向分量,得到标定板垂直方向的两个分量;
根据已知大圆点C1的坐标,以及标定板水平和垂直方向各自的两个分量,预估得到所有圆形靶标的坐标,并存储所有圆形靶标位置坐标。
5.根据权利要求1所述的3D结构光系统的标定方法,其特征在于,采用双线性插值算法得到每个圆心坐标在水平和垂直方向的绝对相位值,具体包括如下步骤:
根据提取得到的各圆心亚像素位置坐标,利用双线性插值得到该圆心对应的相位主值,插值公式如下:
φ(i+u,j+v)=(1-u)(1-v)φ(i,j)+(1-u)vφ(i,j+1)+u(1-v)φ(i+1,j)+uvφ(i+1,j+1);
其中,
φ为相位主值;
i,j分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的整数部分;
u,v分别为计算得到的圆心坐标在图像坐标系中垂直和水平方向的小数部分;
相位主值通过如下方法获得:
在采集的投影了正弦条纹码的图像中,多步相移法光栅光强分布函数为:
Ik(x,y)=I′(x,y)+I″(x,y)cos(φ+(k-1)π/2)
其中,
k取值为1-n,其中n为多步相移法中的步数,此处取10;
x,y为像素坐标;
I(x,y)为第k幅相移图像中(x,y)位置的像素点的灰度值;
I'(x,y)为图像的平均灰度;
I″(x,y)为图像的灰度调制;
φ为相位主值;
通过光强分布函数,得到每个像素点对应的相位主值:
其中,
k取值为1-n,其中n为多步相移法中的步数,此处取10;
x,y为像素坐标;
Ik为该点像素灰度值;
N为光栅条纹周期个数;
φ(x,y)为相位主值。
7.根据权利要求1所述的3D结构光系统的标定方法,其特征在于,还包括,根据标定得到的相机和投影仪内参、畸变系数和外参,以及相机与投影仪之间的坐标系转换矩阵,建立相机和投影仪与世界坐标系的映射关系,具体包括如下步骤:
通过多频外差法,得到绝对相位灰度图;
对于空间任一点p[xw,yw,zw],其投影在相机中的图像坐标为p(uc,vc),根据得到的拍摄图像上的每个像素点的绝对相位值,得到对应的投影仪图像中的一条直线,相机图像与投影仪图像的对应关系为:
其中:
N为光栅条纹周期个数;
W为投影仪水平方向的分辨率;
Φ(uc,vc)为该像素点的绝对相位值;
up为该像素点对应的投影仪图像中垂直方向的坐标;
根据上述相机图像与投影仪图像的对应关系,由如下成像原理公式,得到唯一的点p的三维坐标(XW,YW,ZW):
SC[uc,vc,1]=ACMC[XW,YW,ZW,1]
SP[up,vp,1]=APMP[XW,YW,ZW,1]
其中,
AC和AP分别为相机和投影仪的内参;
MC和MP分别为相机和投影仪的外参;
SC和SP分别为相机和投影仪的比例因子;
(uc,vc)和(up,vp)分别为相机和投影仪的图像坐标,两者均使用预先标定出的系统畸变参数进行畸变矫正;
(XW,YW,ZW)为点p的唯一三维坐标。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010951187.9A CN112097689B (zh) | 2020-09-11 | 2020-09-11 | 一种3d结构光系统的标定方法 |
PCT/CN2020/130543 WO2022052313A1 (zh) | 2020-09-11 | 2020-11-20 | 3d结构光系统的标定方法、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010951187.9A CN112097689B (zh) | 2020-09-11 | 2020-09-11 | 一种3d结构光系统的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112097689A CN112097689A (zh) | 2020-12-18 |
CN112097689B true CN112097689B (zh) | 2022-02-22 |
Family
ID=73752106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010951187.9A Active CN112097689B (zh) | 2020-09-11 | 2020-09-11 | 一种3d结构光系统的标定方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112097689B (zh) |
WO (1) | WO2022052313A1 (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111811433B (zh) * | 2020-07-15 | 2022-03-08 | 河北工业大学 | 基于红蓝正交条纹的结构光系统标定方法及装置和应用 |
CN112991462A (zh) * | 2021-03-15 | 2021-06-18 | 扬州大学 | 一种基于圆点图的摄像机标定方法 |
CN113989386B (zh) * | 2021-10-27 | 2023-05-30 | 武汉高德智感科技有限公司 | 一种红外相机标定方法及系统 |
CN114882110B (zh) * | 2022-05-10 | 2024-04-12 | 中国人民解放军63921部队 | 适用于微纳卫星自组装的相对位姿测量及靶标设计方法 |
CN115082557B (zh) * | 2022-06-29 | 2024-03-15 | 中交第二航务工程局有限公司 | 一种基于双目视觉的塔柱吊装相对姿态测量方法 |
CN115201796B (zh) * | 2022-07-26 | 2024-04-30 | 白犀牛智达(北京)科技有限公司 | 一种车辆传感器的外参校正方法 |
CN115307576B (zh) * | 2022-08-02 | 2024-04-09 | 清华大学 | 一种结构光测量中的阶梯边界补偿方法及装置 |
CN115930784B (zh) * | 2023-01-09 | 2023-08-25 | 广州市易鸿智能装备有限公司 | 一种视觉检测系统的点检方法 |
CN116182702B (zh) * | 2023-01-31 | 2023-10-03 | 桂林电子科技大学 | 一种基于主成分分析的线结构光传感器标定方法及系统 |
CN116563388B (zh) * | 2023-04-28 | 2024-05-07 | 神力视界(深圳)文化科技有限公司 | 标定数据采集方法、装置、电子设备及存储介质 |
CN116277979B (zh) * | 2023-05-24 | 2023-09-08 | 南京铖联激光科技有限公司 | 一种dlp打印机光机畸变矫正方法 |
CN116817794B (zh) * | 2023-06-27 | 2024-02-13 | 浙江大学 | 一种基于结构光的水下高精度三维成像装置及方法 |
CN117029705B (zh) * | 2023-06-27 | 2024-03-22 | 苏州瑞威盛科技有限公司 | 一种基于非接触式3d视觉的齿轮跨棒距测量系统及方法 |
CN117017496B (zh) * | 2023-09-28 | 2023-12-26 | 真健康(北京)医疗科技有限公司 | 柔性体表定位装置及穿刺手术导航定位系统 |
CN117451657A (zh) * | 2023-10-23 | 2024-01-26 | 南京审计大学 | 一种烟气干扰环境下远距离长波红外线激光测量仪及测量方法 |
CN117369197B (zh) * | 2023-12-06 | 2024-05-07 | 深圳市安思疆科技有限公司 | 3d结构光模组、成像系统及获得目标物体深度图的方法 |
CN117911542A (zh) * | 2024-03-19 | 2024-04-19 | 杭州灵西机器人智能科技有限公司 | 一种标定板、标定板识别方法、系统、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105393A (zh) * | 2006-07-13 | 2008-01-16 | 周波 | 投射多频光栅的物体表面三维轮廓的视觉测量方法 |
CN102853783A (zh) * | 2012-09-18 | 2013-01-02 | 天津工业大学 | 一种高精度多波长三维测量方法 |
CN110132431A (zh) * | 2019-03-29 | 2019-08-16 | 黑龙江科技大学 | 一种图像灰度区间扩展的多频外差光栅绝对相位计算方法 |
CN110595387A (zh) * | 2019-08-01 | 2019-12-20 | 佛山市南海区广工大数控装备协同创新研究院 | 一种基于多频率结构光的三维重建系统标定方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7061628B2 (en) * | 2001-06-27 | 2006-06-13 | Southwest Research Institute | Non-contact apparatus and method for measuring surface profile |
JP2006003212A (ja) * | 2004-06-17 | 2006-01-05 | Konica Minolta Sensing Inc | 位相計測システム |
CN100342207C (zh) * | 2004-12-17 | 2007-10-10 | 北京航空航天大学 | 一种微小构件通孔内表面三维形貌检测装置及其标定和使用方法 |
CN103234482B (zh) * | 2013-04-07 | 2015-12-02 | 哈尔滨工程大学 | 一种基于正弦光栅的结构光测量系统标定方法 |
CN108717715B (zh) * | 2018-06-11 | 2022-05-31 | 华南理工大学 | 一种用于弧焊机器人的线结构光视觉系统自动标定方法 |
-
2020
- 2020-09-11 CN CN202010951187.9A patent/CN112097689B/zh active Active
- 2020-11-20 WO PCT/CN2020/130543 patent/WO2022052313A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105393A (zh) * | 2006-07-13 | 2008-01-16 | 周波 | 投射多频光栅的物体表面三维轮廓的视觉测量方法 |
CN102853783A (zh) * | 2012-09-18 | 2013-01-02 | 天津工业大学 | 一种高精度多波长三维测量方法 |
CN110132431A (zh) * | 2019-03-29 | 2019-08-16 | 黑龙江科技大学 | 一种图像灰度区间扩展的多频外差光栅绝对相位计算方法 |
CN110595387A (zh) * | 2019-08-01 | 2019-12-20 | 佛山市南海区广工大数控装备协同创新研究院 | 一种基于多频率结构光的三维重建系统标定方法 |
Non-Patent Citations (3)
Title |
---|
"基于数学光栅投影的结构光三维测量技术与系统研究";李中伟;《中国博士学位论文全文数据库信息科技辑》;20091115;第I138-1317页第I138-52页:正文第24-57页:正文第27-44页 * |
"结构光三维形貌与变形测量系统标定研究";刘海军;《中国优秀硕士学位论文全文数据库信息科技辑》;20150215(第02期);第I138-1317页:正文第27-44页 * |
刘海军."结构光三维形貌与变形测量系统标定研究".《中国优秀硕士学位论文全文数据库信息科技辑》.2015,(第02期), * |
Also Published As
Publication number | Publication date |
---|---|
WO2022052313A1 (zh) | 2022-03-17 |
CN112097689A (zh) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112097689B (zh) | 一种3d结构光系统的标定方法 | |
US8090194B2 (en) | 3D geometric modeling and motion capture using both single and dual imaging | |
CN110276808B (zh) | 一种单相机结合二维码测量玻璃板不平度的方法 | |
US9367952B2 (en) | 3D geometric modeling and 3D video content creation | |
CN110686599B (zh) | 基于彩色格雷码结构光的三维测量方法、系统、装置 | |
Xu et al. | Real-time 3D shape inspection system of automotive parts based on structured light pattern | |
WO2007015059A1 (en) | Method and system for three-dimensional data capture | |
CN113205593B (zh) | 一种基于点云自适应修复的高反光表面结构光场三维重建方法 | |
US20130127998A1 (en) | Measurement apparatus, information processing apparatus, information processing method, and storage medium | |
CN109716348B (zh) | 独立地处理多个所关注的区域 | |
CN109974625B (zh) | 一种基于色相优化灰度的彩色物体结构光三维测量方法 | |
CN110163918A (zh) | 一种基于射影几何的线结构光标定方法 | |
US20160025591A1 (en) | Automated deflectometry system for assessing reflector quality | |
CN112815843B (zh) | 一种3d打印过程中工件表面打印偏差的在线监测方法 | |
CN106952262B (zh) | 一种基于立体视觉的船板加工精度分析方法 | |
CN110763204A (zh) | 一种平面编码靶标及其位姿测量方法 | |
CN113505626A (zh) | 一种快速三维指纹采集方法与系统 | |
CN114972633A (zh) | 交叉激光线约束下的快速扫描点云插值方法 | |
CN112614215B (zh) | 一种三维测量表面温度场重建方法 | |
JP2004077290A (ja) | 3次元形状計測装置および方法 | |
CN114998571B (zh) | 基于固定尺寸标志物的图像处理及颜色检测方法 | |
Usamentiaga | Easy rectification for infrared images | |
CN115841517A (zh) | 一种基于dic双圆交比的结构光标定方法和装置 | |
TW201913574A (zh) | 應用於物件檢測之攝影機與雷射測距儀的數據融合方法 | |
Abzal et al. | Development of an automatic map drawing system for ancient bas-reliefs |
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 |