CN109458955B - 基于平面度约束的离轴圆条纹投影测量零相位点求解方法 - Google Patents
基于平面度约束的离轴圆条纹投影测量零相位点求解方法 Download PDFInfo
- Publication number
- CN109458955B CN109458955B CN201811574048.8A CN201811574048A CN109458955B CN 109458955 B CN109458955 B CN 109458955B CN 201811574048 A CN201811574048 A CN 201811574048A CN 109458955 B CN109458955 B CN 109458955B
- Authority
- CN
- China
- Prior art keywords
- phase point
- zero phase
- circular fringe
- zero
- point
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
基于平面度约束的离轴圆条纹投影测量零相位点求解方法,先按照测量需求,搭建离轴圆条纹投影三维测量***;然后采用离轴圆条纹投影三维测量***对一个标准平面进行测量,采样得到所需圆条纹图,并处理得到圆条纹图相位信息;再构建离轴圆条纹投影三维测量***的虚拟高度求解模型;然后根据虚拟高度求解模型,建立对标准平面进行测量时平面度与零相位点取值的关联关系;最后综合圆条纹图相位信息、虚拟高度求解模型以及平面度与零相位点位置关联关系,建立零相位点搜索方案,实现对零相位点的求解;本发明能够实现对离轴圆条纹投影轮廓术零相位点的高精度求解,有助于促进圆条纹投影轮廓术测量性能的提高。
Description
技术领域
本发明属于三维测量技术领域,具体涉及基于平面度约束的离轴圆条纹投影测量零相位点求解方法。
背景技术
三维测量在质量检测、快速制造、逆向工程、模式识别等领域有着越来越广泛的需求,对于三维测量技术的研究也日益迫切。圆条纹投影轮廓术是发明人所在团队近期提出的一种新型光学三维测量技术,为机械零部件等的高精度三维测量提供了一种新选择。但该技术实际测量精度受限于测量噪声、相位解调精度等因素,特别是在共轴布置时存在对噪声的放大效应,会显著降低测量精度。鉴于此,为提高测量精度,需要采用离轴布置方式实现圆条纹投影三维测量。但此时采样条纹图零相位点缺失,阻碍了三维重建的实现。
由于圆条纹投影轮廓术是一种新型光学三维测量技术,相关研究薄弱,针对上述问题的研究仍旧欠缺,亟需推进。
发明内容
为了克服上述现有技术的额确定,本发明的目的在于提供基于平面度约束的离轴圆条纹投影测量零相位点求解方法,显著提高圆条纹投影轮廓术的三维测量性能。
为了达到上述目的,本发明采取的技术方案为:
基于平面度约束的离轴圆条纹投影测量零相位点求解方法,包括以下步骤:
步骤1:按照测量需求,搭建离轴圆条纹投影三维测量***;
步骤2:采用所搭建离轴圆条纹投影三维测量***对一个标准平面进行测量,采样得到所需圆条纹图,并处理得到圆条纹图相位信息;
步骤3:构建离轴圆条纹投影三维测量***的虚拟高度求解模型;
步骤4:根据虚拟高度求解模型,建立对标准平面进行测量时平面度与零相位点取值的关联关系;
步骤5:综合圆条纹图相位信息、虚拟高度求解模型以及平面度与零相位点位置关联关系,建立零相位点搜索方案,实现对零相位点的求解。
所述的步骤1中的离轴圆条纹投影三维测量***,由投影***和配置远心镜头的相机构成,投影***光轴与相机光轴需要保持平行且在空间上偏离一定距离,在相机视场位于投影***视场内的前提下使投影***光轴与相机光轴尽可能远离。
所述的步骤2中的基于离轴圆条纹投影三维测量***的测量过程与共轴圆条纹投影轮廓术一致,设求得的圆条纹图相位分布为Φ(xc,yc),(xc,yc)为采样圆条纹图的图像坐标。
所述的步骤3中虚拟高度求解模型为:
其中:μc为相机像元尺寸,(xc,yc)为采样圆条纹图的图像坐标,(xc0,yc0)为采样圆条纹图“零相位点”的图像坐标,β为远心镜头放大倍率,Φ(xc,yc)为采样圆条纹图在点(xc,yc)的相位。
所述的采样圆条纹图“零相位点”实际上为一个虚拟点,因为不存在于采样圆条纹图上,因而需要求解。
所述的步骤4中对标准平面进行测量时平面度与零相位点取值的关联关系为,所选取零相位点距其真实零相位点越近,重建得到轮廓的平面度越高,反之亦然,也就是说当且仅当选用的零相位点位于其真实位置时,重建得到轮廓的平面度最高。
所述的步骤5中的零相位点搜索方案为,根据仪器布置及成像***参数的先验知识,确定零相位点初步分布区域;按照测量精度要求,遍历搜索初步分布区域的每一点作为潜在零相位点;对每一个遍历的零相位点,借助公式(1)求解被测标准平面对应的高度分布;对每一个零相位点对应的标准平面求解高度分布进行平面拟合,得到每次拟合的均方根误差RMSE;搜索得到平面度最高的点即为零相位点真实值。
对应的拟合平面为:
其中:aq,bq,和cq是平面拟合求得的系数;平面拟合的RMSE为:
其中:row和column是表征遍历总点数的参量。
当遍历得到所有点的RMSE后,其最小值对应的点即为实际零相位点:
其中:Q为零相位点总遍历个数。
本发明具有以下有益效果:通过建立虚拟高度求解模型及平面度约束,提出了一种适用于离轴圆条纹投影轮廓术的零相位点高精度求解方法;本发明方法能够以极小的成本显著提高圆条纹投影轮廓术三维测量精度,促进其在工业测量中的应用推广。
附图说明
图1为本发明实施例离轴圆条纹投影三维测量***的示意图。
图2为实施例采用图1离轴圆条纹投影三维测量***进行测量时所投影的圆条纹图。
图3为实施例图2圆条纹图由图1离轴圆条纹投影三维测量***模拟投影到标准平面后,进而模拟采样得到的圆条纹图。
图4为实施例图3所示采样圆条纹图处理所得相位。
图5为实施例进行零相位点搜索的初步分布区域及搜索路径示意图。
图6为实施例对图3标准平面选用不同零相位点进行虚拟高度重建时得到的三维点云分布。
图7为实施例测量中的采样圆条纹图。
图8为实施例重建得到的标准球三维轮廓。
具体实施方式
下面结合附图和实施例对本发明作详细描述。
基于平面度约束的离轴圆条纹投影测量零相位点求解方法,包括以下步骤:
步骤1:按照测量需求,搭建离轴圆条纹投影三维测量***,如图1所示,可见,离轴圆条纹投影三维测量***与传统共轴圆条纹投影测量***类似,由投影***和配置远心镜头的相机构成,且投影***光轴与相机光轴要保持平行;不同的是,离轴圆条纹投影三维测量***的投影***光轴与相机光轴需要在空间上偏离一定距离,在相机视场位于投影***视场内的前提下使投影***光轴与相机光轴尽可能远离,这通常会使得投影光轴位于相机视场外,参照图2,图2为采用图1所示的离轴圆条纹投影三维测量***进行测量时所投影的圆条纹图;
步骤2:采用所搭建离轴圆条纹投影三维测量***对一个标准平面进行测量,采样得到所需的一系列圆条纹图,并处理得到圆条纹图相位信息,参照图3和图4,图3对一个标准平面进行模拟测量时得到的采样圆条纹图,图4为处理图3后得到的圆条纹图相位Φ(xc,yc),(xc,yc)为采样圆条纹图的图像坐标;
模拟测量的优势在于,能够预先确定待求零相位点的实际值,同时避免其它因素的干扰;在本模拟过程中,相机分辨率为1280×960,虚拟零相位点实际值为(xc0,yc0)=(4960,1);
步骤3:构建离轴圆条纹投影三维测量***的虚拟高度求解模型,如下所示:
其中:μc为相机像元尺寸,(xc,yc)为采样圆条纹图的图像坐标,(xc0,yc0)为采样圆条纹图“零相位点”的图像坐标,β为远心镜头放大倍率,Φ(xc,yc)为采样圆条纹图在点(xc,yc)的相位;
由公式(1)可推得,零相位点(xc0,yc0)是进行圆条纹投影轮廓术三维重建的关键参量;但是,图2所示投影圆条纹图的零相位点(也就是图2中的圆环中心)并没有出现在图3中,这意味着无法像共轴式圆条纹投影轮廓术一样从采样圆条纹图中借助相位信息直接求得零相位点坐标,因而需要采用专门方法对其进行求解;
步骤4:根据虚拟高度求解模型,建立对标准平面进行测量时平面度与零相位点取值的关联关系;通过对公式(1)进行分析可得,对标准平面进行测量时,测得三维点云的平面度与零相位点取值的关联关系为:所选取零相位点距其真实零相位点越近,重建得到轮廓的平面度越高,反之亦然,也就是说当且仅当选用的零相位点位于其真实位置时,重建得到轮廓的平面度最高;
步骤5:综合圆条纹图相位信息、虚拟高度求解模型以及平面度与零相位点位置关联关系,建立零相位点搜索方案,实现对零相位点的便捷、高精度求解。
零相位点的搜索方案为:
5.1:根据仪器布置及成像***参数等先验知识,确定零相位点初步分布区域;
5.2:按照测量精度要求,遍历搜索初步分布区域的每一点作为潜在零相位点;对每一个遍历的零相位点,借助公式(1)求解被测标准平面对应的高度分布;
5.3:对每一个零相位点对应的测量高度分布进行平面拟合,得到每次拟合的均方根误差(RMSE);
5.4:搜索得到平面度最高的点即为零相位点真实值。
参照图5,图5为进行零相位点搜索时初步分布区域及搜索路径示意图,图中的十字线中心表示零相位点真实位置。实际求解中,搜索范围大小、搜索点密度等参数需要按照测量精度要求选定;设第q个零相位点遍历点的像素坐标为对应的求解高度分布可如下求解:
对应的拟合平面为:
其中:aq,bq,和cq是平面拟合求得的系数;平面拟合的RMSE为:
其中:row和column是表征遍历总点数的参量。
参照图6,图6为综合图4圆条纹图相位及公式(2)求解得到的分别对应于零相位点(xc0,yc0)=(4960,11)、(xc0,yc0)=(4960,1)的标准平面形貌。由图可见,当零相位点偏离其真实位置时,重建得到的轮廓形貌呈现明显弯曲(亦即平面度显著降低),而当零相位点位于其真实位置时重建得到的轮廓呈一个平面(亦即平面度高)。因而,平面度约束是一种有效的零相位点求解指标。
当遍历得到所有点的RMSE后,其最小值对应的点即为实际零相位点:
其中:Q为零相位点总遍历个数。
利用本发明方法,搭建了实际离轴圆条纹投影三维测量***,对一个标准球进行了测量,参照图7和图8,图7为测量中的采样圆条纹图,图8为重建得到的标准球三维轮廓。
以上所述仅为本发明的一个实施例,并不用以限制本发明,凡在本发明基础上所做的任何修改、等同替换及拓展等,均应包含在本发明的保护范围内。
Claims (2)
1.基于平面度约束的离轴圆条纹投影测量零相位点求解方法,其特征在于,包括以下步骤:
步骤1:按照测量需求,搭建离轴圆条纹投影三维测量***;
步骤2:采用所搭建离轴圆条纹投影三维测量***对一个标准平面进行测量,采样得到所需圆条纹图,并处理得到圆条纹图相位信息;
步骤3:构建离轴圆条纹投影三维测量***的虚拟高度求解模型;
步骤4:根据虚拟高度求解模型,建立对标准平面进行测量时平面度与零相位点取值的关联关系;
步骤5:综合圆条纹图相位信息、虚拟高度求解模型以及平面度与零相位点位置关联关系,建立零相位点搜索方案,实现对零相位点的求解;
所述的步骤2中的基于离轴圆条纹投影三维测量***的测量过程与共轴圆条纹投影轮廓术一致,设求得的圆条纹图相位分布为Φ(xc,yc),(xc,yc)为采样圆条纹图的图像坐标;
所述的步骤3中虚拟高度求解模型为:
其中:μc为相机像元尺寸,(xc,yc)为采样圆条纹图的图像坐标,(xc0,yc0)为采样圆条纹图“零相位点”的图像坐标,β为远心镜头放大倍率,Φ(xc,yc)为采样圆条纹图在点(xc,yc)的相位;
所述的采样圆条纹图“零相位点”实际上为一个虚拟点,因为它不存在于采样圆条纹图上,因而需要求解;
所述的步骤4中对标准平面进行测量时平面度与零相位点取值的关联关系为,所选取零相位点距其真实零相位点越近,重建得到轮廓的平面度越高,反之亦然,也就是说当且仅当选用的零相位点位于其真实位置时,重建得到轮廓的平面度最高;
所述的步骤5中的零相位点搜索方案为,根据仪器布置及成像***参数的先验知识,确定零相位点初步分布区域;按照测量精度要求,遍历搜索初步分布区域的每一点作为潜在零相位点;对每一个遍历的零相位点,借助公式(1)求解被测标准平面对应的高度分布;对每一个零相位点对应的标准平面求解高度分布进行平面拟合,得到每次拟合的均方根误差RMSE;搜索得到平面度最高的点即为零相位点真实值;
设第q个零相位点遍历点的像素坐标为对应的求解高度分布如下求解:
对应的拟合平面为:
其中:aq,bq,和cq是平面拟合求得的系数;平面拟合的RMSE为:
其中:row和column是表征遍历总点数的参量;
当遍历得到所有点的RMSE后,其最小值对应的点即为实际零相位点:
其中:Q为零相位点总遍历个数。
2.根据权利要求1所述的基于平面度约束的离轴圆条纹投影测量零相位点求解方法,其特征在于:所述的步骤1中的离轴圆条纹投影三维测量***,由投影***和配置远心镜头的相机构成,投影***光轴与相机光轴需要保持平行且在空间上偏离一定距离,在相机视场位于投影***视场内的前提下使投影***光轴与相机光轴尽可能远离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811574048.8A CN109458955B (zh) | 2018-12-21 | 2018-12-21 | 基于平面度约束的离轴圆条纹投影测量零相位点求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811574048.8A CN109458955B (zh) | 2018-12-21 | 2018-12-21 | 基于平面度约束的离轴圆条纹投影测量零相位点求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109458955A CN109458955A (zh) | 2019-03-12 |
CN109458955B true CN109458955B (zh) | 2020-01-14 |
Family
ID=65614277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811574048.8A Active CN109458955B (zh) | 2018-12-21 | 2018-12-21 | 基于平面度约束的离轴圆条纹投影测量零相位点求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109458955B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110006364B (zh) * | 2019-03-18 | 2020-12-29 | 南京师范大学 | 基于圆条纹径向空间载波相移的三维实时显微测量方法 |
CN113074657B (zh) * | 2021-03-16 | 2022-02-01 | 华中科技大学 | 一种基于虚拟平面采样的三坐标相位映射表标定方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004109229A2 (en) * | 2003-06-11 | 2004-12-16 | Solvision | 3d and 2d measurement system and method with increased sensitivity and dynamic range |
CN102721376B (zh) * | 2012-06-20 | 2014-12-31 | 北京航空航天大学 | 一种大视场三维视觉传感器的标定方法 |
CN104729429B (zh) * | 2015-03-05 | 2017-06-30 | 深圳大学 | 一种远心成像的三维形貌测量***标定方法 |
CN105953749B (zh) * | 2016-06-21 | 2018-12-18 | 西安交通大学 | 一种光学三维形貌测量方法 |
CN106950035B (zh) * | 2017-02-23 | 2019-04-12 | 西安交通大学 | 液晶光栅干涉测量空间光调制器相位调制特性的装置及方法 |
-
2018
- 2018-12-21 CN CN201811574048.8A patent/CN109458955B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109458955A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103530880B (zh) | 基于投影高斯网格图案的摄像机标定方法 | |
CN102183214B (zh) | 一种大口径非球面镜结构光检测方法 | |
CN105783786B (zh) | 一种基于结构光视觉的零件倒角测量方法与装置 | |
CN111369607B (zh) | 一种基于图片解析的预制构件拼装匹配方法 | |
CN109556540A (zh) | 一种非接触式基于3d图像的物体平面度检测方法、计算机 | |
CN102376089A (zh) | 一种标靶校正方法及*** | |
CN110223355B (zh) | 一种基于双重极线约束的特征标志点匹配方法 | |
CN107576281A (zh) | 一种测量管路弯曲半径的方法和装置 | |
Wu et al. | A novel high precise laser 3D profile scanning method with flexible calibration | |
CN109544628A (zh) | 一种指针式仪表的准确读数识别***及方法 | |
CN109458955B (zh) | 基于平面度约束的离轴圆条纹投影测量零相位点求解方法 | |
CN109448043A (zh) | 平面约束下的立木高度提取方法 | |
CN112013788A (zh) | 基于叶片局部前缘曲线特征标定转动中心的方法 | |
CN109141266A (zh) | 一种钢结构测量方法及*** | |
CN103424087A (zh) | 一种大尺度钢板三维测量拼接***及方法 | |
CN106767405B (zh) | 相位映射辅助三维成像***快速对应点匹配的方法及装置 | |
CN109242909A (zh) | 一种面向高精度二维尺寸测量的线阵相机标定算法 | |
Wang et al. | A direct calibration method for line structured light measurement system based on parallel lines | |
CN111968182A (zh) | 一种双目相机非线性模型参数的标定方法 | |
Qin et al. | A novel hierarchical iterative hypothesis strategy for intrinsic parameters calibration of laser structured-light weld vision sensor | |
CN115661446A (zh) | 基于深度学习的指针式仪表示数自动读取***及方法 | |
CN113390394B (zh) | 一种带有摄影比例尺的光束法平差算法 | |
CN112414316B (zh) | 一种应变片敏感栅尺寸参数测量方法 | |
TWI444586B (zh) | 零件形位公差檢測系統及方法 | |
CN110260817B (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 |