CN112611320A - 一种基于梯度光强的多频外差面结构光三维重建方法 - Google Patents
一种基于梯度光强的多频外差面结构光三维重建方法 Download PDFInfo
- Publication number
- CN112611320A CN112611320A CN202011313378.9A CN202011313378A CN112611320A CN 112611320 A CN112611320 A CN 112611320A CN 202011313378 A CN202011313378 A CN 202011313378A CN 112611320 A CN112611320 A CN 112611320A
- Authority
- CN
- China
- Prior art keywords
- images
- frequency
- image
- gradient
- structured light
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005259 measurement Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000004069 differentiation Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 17
- 230000010363 phase shift Effects 0.000 description 7
- 238000000691 measurement method Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/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/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é
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
Description
技术领域
本发明属于面结构光三维测量领域,具体涉及一种基于梯度光强的多频外差面结构光三维重建方法。
背景技术
随着现代科技的不断发展,三维测量方法中的重建技术的在各行各业的应用需求不断的增加。该技术的主要思想是通过一定的技术方法获得待测物体的三维形貌信息,重建物体三维模型。近年来三维重建技术在工业检测、医疗、文物恢复等领域得到越来越多的应用。
依据不同的实现形式,三维测量技术根据是否与被测物体接触主要可分为:接触式与非接触式。接触式三维测量技术具有较高精度,但因其需要和被测物体表面接触而具有破坏性。相比于接触式测量,非接触式测量有着无破坏性,高效率,工作距离大等优点。而在众多的非接触式三维测量技术中,基于面结构光的三维视觉重建拥有高精度、操作简单、实时性强等特点,其中多频外差相移技术(Multi-frequency Heterodyne Phase-shiftingTechnology)的面结构光三维重建技术更是受到许多研究者的重视,是视觉重建领域的研究热点。
对于三维测量技术中越来越高的速度要求,多频外差相移技术依然存在一些问题,为了从采集图像中获取被测物体的准确的相位信息,传统的多频外差相移技术通常需要投射一系列不同频率的相移图像计算包裹相位,利用多频外差原理进行相位展开计算,这导致投射的条纹图像增多,采集图像的时间增加,测量过程受到干扰的可能性加大,结构光三维重建的速度也随之降低。大量待投射图像限制了多频外差相移技术在高速、高精度实时结构光三维重建技术中的应用。
发明内容
为实现本发明目的而采用的技术方案是这样的,一种基于梯度光强的多频外差面结构光三维重建方法,包括如下步骤:
1)搭建结构光测量系统。所述结构光测量系统包括工业相机、投影仪、操作终端和被测物体。所述投影仪和工业相机可通信的连接于所述操作终端。所述工业相机和投影仪分别布置在被测物体的前侧。
2)调整工业相机和投影仪的位置,待位置调整完毕后,对整个结构光测量系统进行标定,保存标定参数。
3)在所述操作终端确定m组待生成的结构光条纹图不同的频率为fi,(i=1,2,...,m)。
10)所述操作终端根据最高频率采集图像的绝对相位φ1(x,y)与步骤2)获得的标定参数计算被测物体表面点的三维坐标数据,完成被测物体三维重建。
进一步,在步骤4)中,单个条纹频率fi对应一组结构光条纹图像。每组结构光条纹图像包含有两幅结构光条纹图像和共需投射2m幅结构光条纹图像,将2m幅结构光条纹图像记为投射图像每一组结构光条纹图像具有如下表达形式:
式(2)至(3)中,表示单个条纹频率对应的该组中第1幅投射图像中(x,y)点处的灰度值大小。表示单个条纹频率对应的该组中第2幅投射图像中(x,y)点处的灰度值大小。Ap代表图像的直流分量,Bp代表图像的幅值。φi(x,y)为该幅图像中的绝对相位。
进一步,在步骤6)中,单个条纹频率fi所对应的一组采集图像的梯度光强具有如下表达形式:
Ii1'(x,y)=-Bcsinφi(x,y)dφi (4)
Ii2'(x,y)=Bccosφi(x,y)dφi (5)
式(4)至(5)中,Iij'(x,y)表示单个条纹频率对应的该组中第j幅采集图像中(x,y)点处的梯度光强值大小。Bc代表采集图像的灰度调制。dφi为采集图像中相位的微分。
进一步,在步骤7)中,单个条纹频率对应的每组采集图像包裹相位计算公式如下:
进一步,所述工业相机和投影仪的光轴相交,投影区域相互重叠,且被测物体正好位于工业相机和投影仪的公共视场内。
本发明的技术效果是毋庸置疑的,其相比传统多频外差相移技术的有益效果在于:
1、本发明方法不需要投射一系列不同频率的相移图像,仅有两幅采集图像就可计算得到包裹相位,大幅减少了投射图像数量,采集时间和计算速度都有大幅度提升,测量效率显著提高。
2、本发明方法提高了投射条纹图像的周期数,能够有效提高三维重建的精度,相比传统的多频外差相移方法,可靠性和精度均有所保证。
3、本发明方法在不改变硬件条件下,仅靠算法实现测量效率和测量精度,相比传统多频外差技术,该方法具有简单易操作,实时性强的优点,更适合高速三维重建技术。
附图说明
图1为本发明的基于梯度光强的多频外差面结构光三维重建方法的测量系统示意图;
图2为本发明的基于梯度光强的多频外差面结构光三维重建方法的流程示意图;
图3为本发明的该基于梯度光强的多频外差面结构光三维重建方法所生成的条纹图像示意图。
图中:工业相机1、投影仪2、操作终端3和被测物体4。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
本实施例公开一种基于梯度光强的多频外差面结构光三维重建方法,该方法适用于多种频率外差测量方法,参见图1和图2,具体实现步骤如下:
1)搭建结构光测量系统。所述结构光测量系统包括工业相机1、投影仪2、操作终端3和被测物体4。所述投影仪2和工业相机1可通信的连接于所述操作终端3。所述工业相机1和投影仪2分别布置在被测物体4的前侧。所述工业相机1和投影仪2的光轴相交,投影区域相互重叠,且被测物体4正好位于工业相机1和投影仪2的公共视场内。
2)调整工业相机1和投影仪2的位置,待位置调整完毕后,对整个结构光测量系统进行标定,获得投影仪2-工业相机1系统的三维坐标-相位关系,保存标定参数。
3)在所述操作终端3确定m组正弦光栅条纹图像不同的频率为fi,(i=1,2,...,m)。在本实施例中选取传统三种条纹频率和基于摩尔序列的外差相位展开方法,确定条纹频率数m=3。根据相关研究,优选三种不同的外差频率解包裹相位,频率为:f1=70、f2=64、f3=59。则本实施例需采集三组不同频率下的图像,每组图像共有两幅竖直条纹光栅图像。设定每幅待投射图像的宽度为W=1024,高度为H=768。
4)依据本实施例设定生成的结构光条纹光栅表达式中的待定参数分别为:
A、B(1)
式(1)至(2)中,A为图像的直流分量,B为图像的幅值,φ(x,y)为图像中的绝对相位,p为节距,p=W/f,代表一个条纹周期内像素点总数。则对应于pi(i=1,2,...,m),节距的一组竖直条纹光栅(x,y)处光强的表达式为:
根据上述步骤3)中预设f1=70、f2=64、f3=59,可对应生成三组条纹光栅图像,每组均包含两幅竖直的条纹图像。对应于本实施例生成的六幅结构光条纹图像,如附图3所示。
5)投影仪2将步骤2中生成的六幅结构光条纹图像投射至被测物体4,同时通过该工业相机2同步采集经物体高度调制后的条纹图像,并将六幅采集图像存储于该控制终端3。
6)计算步骤5)中被采集图像的梯度光强。对应于fi,(i=1,2,3)频率的一组竖直采集图像(x,y)处梯度光强的表达式为:
I1'(x,y)=-Bcsinφi(x,y)dφi (5)
I2'(x,y)=Bccosφi(x,y)dφi (6)
式(5)至(6)中,Bc代表竖直采集图像的灰度调制,dφi为竖直采集图像中相位的微分。根据上述该步骤3中预设的三种条纹频率,共得到三组采集图像的梯度光强值。
7)计算步骤6)中每组竖直采集图像的梯度光强对应的包裹相位。本实施例所得各组采集图像的梯度光强值为:
对应f1的一组竖直采集图像的梯度光强值:I11'、I12';对应f2的一组竖直采集图像的梯度光强值:I21'、I22';对应f2的一组竖直采集图像的梯度光强值:I31'、I32';
根据上述假设可解出各组采集图像的包裹相位:
对应f1的一组竖直采集图像的包裹相位为:
对应f2的一组竖直采集图像的包裹相位为:
对应f2的一组竖直采集图像的包裹相位为:
8)依据步骤7)解出的各组竖直采集图像对应的包裹相位,通过基于摩尔序列的外差相位展开算法得到最高频率f1=70对应的该组竖直采集图像的相位级数整数部分,即计算条纹级数k1;
10)根据步骤9)中获得的最高频率采集图像的绝对相位信息与预先标定得到的投影仪-相机系统的三维坐标-相位关系计算待测物体表面点的三维坐标数据,完成待测物体三维重建。
本实施例提供的一种基于梯度光强的多频外差面结构光三维重建方法,通过生成一组基于梯度光强的多频面结构光条纹图像,用投影仪2将不同频率的投射图像投射到被测物体4上。首先计算被采集图像的梯度光强,接着通过梯度光强得到包裹相位信息。再根据多频外差的相位展开原理进行绝对相位的求解。减少了投射图像张数,增加了条纹图像的周期数,提高了三维重建的速度和精度。
实施例2:
本实施例提供一种较为基础的实现方式,一种基于梯度光强的多频外差面结构光三维重建方法,该方法适用于多种频率外差测量方法,参见图1和图2,包括如下步骤:
1)搭建结构光测量系统。所述结构光测量系统包括工业相机1、投影仪2、操作终端3和被测物体4。所述投影仪2和工业相机1可通信的连接于所述操作终端3。所述工业相机1和投影仪2分别布置在被测物体4的前侧。
2)调整工业相机1和投影仪2的位置,待位置调整完毕后,对整个结构光测量系统进行标定,获得投影仪2-工业相机1系统的三维坐标-相位关系,保存标定参数。
3)在所述操作终端3确定m组待生成的结构光条纹图不同的频率为fi,(i=1,2,...,m)。
5)所述操作终端3将投射图像通过投影仪2分别投射到被测物体4表面,并通过工业相机1采集图像,获得经过被测物体4调制后的m组采集图像,并将m组采集图像记为总采集图像并输出至操作终端3,每组待投射结构光条纹图像对应的均与工业相机1采集的一组采集图像对应。
根据外差原理,两个高频率的条纹图像可以合成低频率的外差条纹图像。m个不同频率的采集图像,经过多级外差运算合成的最后一级外差条纹图像的节距(一个条纹周期T包含的像素点的个数)覆盖图像全场,即最后一级合成外差条纹图像的频率为1。
10)所述操作终端3根据最高频率采集图像的绝对相位φ1(x,y)与步骤2)获得的标定参数,标定得到的投影仪-相机系统的三维坐标-相位关系计算被测物体4表面点的三维坐标数据,完成被测物体4三维重建。
实施例3:
本实施例主要步骤同实施例2,进一步,在步骤4)中,单个条纹频率fi对应一组结构光条纹图像。每组结构光条纹图像包含有两幅结构光条纹图像和共需投射2m幅结构光条纹图像,将2m幅结构光条纹图像记为投射图像每一组结构光条纹图像具有如下表达形式:
式(2)至(3)中,表示单个条纹频率对应的该组中第1幅投射图像中(x,y)点处的灰度值大小。表示单个条纹频率对应的该组中第2幅投射图像中(x,y)点处的灰度值大小。Ap代表图像的直流分量,Bp代表图像的幅值。φi(x,y)为该幅图像中的绝对相位。
实施例4:
本实施例主要步骤同实施例3,进一步,在步骤6)中,单个条纹频率fi所对应的一组采集图像的梯度光强具有如下表达形式:
Ii1'(x,y)=-Bcsinφi(x,y)dφi (4)
Ii2'(x,y)=Bccosφi(x,y)dφi (5)
式(4)至(5)中,Iij'(x,y)表示单个条纹频率对应的该组中第j幅采集图像中(x,y)点处的梯度光强值大小。Bc代表采集图像的灰度调制。dφi为采集图像中相位的微分。
实施例5:
本实施例主要步骤同实施例4,进一步,在步骤7)中,单个条纹频率对应的每组采集图像包裹相位计算公式如下:
实施例6:
本实施例主要步骤同实施例2,进一步,所述工业相机1和投影仪2的光轴相交,投影区域相互重叠,且被测物体4正好位于工业相机1和投影仪2的公共视场内。
Claims (5)
1.一种基于梯度光强的多频外差面结构光三维重建方法,其特征在于,包括如下步骤:
1)搭建结构光测量系统;所述结构光测量系统包括所述工业相机(1)、投影仪(2)、操作终端(3)和被测物体(4);所述投影仪(2)和工业相机(1)可通信的连接于所述操作终端(3);所述工业相机(1)和投影仪(2)分别布置在被测物体(4)的前侧;
2)调整工业相机(1)和投影仪(2)的位置,待位置调整完毕后,对整个结构光测量系统进行标定,保存标定参数;
3)在所述操作终端(3)确定m组待生成的结构光条纹图不同的频率为fi,(i=1,2,...,m);
5)所述操作终端(3)将投射图像通过投影仪(2)分别投射到被测物体(4)表面,并通过工业相机(1)采集图像,获得经过被测物体(4)调制后的m组采集图像,并将m组采集图像记为总采集图像并输出至操作终端(3);
10)所述操作终端(3)根据最高频率采集图像的绝对相位φ1(x,y)与步骤2)获得的标定参数计算被测物体(4)表面点的三维坐标数据,完成被测物体(4)三维重建。
3.根据权利要求2所述的一种基于梯度光强的多频外差面结构光三维重建方法,其特征在于,在步骤6)中,单个条纹频率fi所对应的一组采集图像的梯度光强具有如下表达形式:
Ii1'(x,y)=-Bcsinφi(x,y)dφi (4)
Ii2'(x,y)=Bccosφi(x,y)dφi (5)
式(4)至(5)中,Iij'(x,y)表示单个条纹频率对应的该组中第j幅采集图像中(x,y)点处的梯度光强值大小;Bc代表采集图像的灰度调制;dφi为采集图像中相位的微分;。
5.根据权利要求1所述的一种基于梯度光强的多频外差面结构光三维重建方法,其特征在于:所述工业相机(1)和投影仪(2)的光轴相交,投影区域相互重叠,且被测物体(4)正好位于工业相机(1)和投影仪(2)的公共视场内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011313378.9A CN112611320A (zh) | 2020-11-21 | 2020-11-21 | 一种基于梯度光强的多频外差面结构光三维重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011313378.9A CN112611320A (zh) | 2020-11-21 | 2020-11-21 | 一种基于梯度光强的多频外差面结构光三维重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112611320A true CN112611320A (zh) | 2021-04-06 |
Family
ID=75225025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011313378.9A Pending CN112611320A (zh) | 2020-11-21 | 2020-11-21 | 一种基于梯度光强的多频外差面结构光三维重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112611320A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113658321A (zh) * | 2021-07-26 | 2021-11-16 | 浙江大华技术股份有限公司 | 三维重建方法、系统及相关设备 |
-
2020
- 2020-11-21 CN CN202011313378.9A patent/CN112611320A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113658321A (zh) * | 2021-07-26 | 2021-11-16 | 浙江大华技术股份有限公司 | 三维重建方法、系统及相关设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10911672B2 (en) | Highly efficient three-dimensional image acquisition method based on multi-mode composite encoding and epipolar constraint | |
CN104061879B (zh) | 一种连续扫描的结构光三维面形垂直测量方法 | |
Su et al. | Dynamic 3-D shape measurement method based on FTP | |
CN109186496B (zh) | 一种基于移动最小二乘法的三维面形测量方法 | |
CA2528791A1 (en) | Full-field three-dimensional measurement method | |
CN104279981B (zh) | 一种基于条纹反射的镜面/类镜面物体绝对面形的测量方法及装置 | |
CN113251952B (zh) | 一种光栅平动结构光三维测量系统及三维测量方法 | |
CN102261896A (zh) | 一种基于相位测量的物体三维形貌测量方法及系统 | |
CN105066904B (zh) | 基于相位梯度阈值的流水线产品三维面型检测方法 | |
CN106032976A (zh) | 基于波长选择的三条纹投影相位展开方法 | |
CN110174079A (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
CN109631798A (zh) | 一种基于π相移方法的三维面形垂直测量方法 | |
CN110672038A (zh) | 一种基于双频相移条纹投影的快速三维测量方法 | |
CN111947599B (zh) | 基于学习的条纹相位恢复和散斑相关的三维测量方法 | |
CN108332684A (zh) | 一种基于结构光照明显微技术的三维轮廓测量方法 | |
CN109523627A (zh) | 一种基于泰勒指数表达式型面结构光的三维重建方法 | |
CN110207620A (zh) | 一种通过不同频率确定数字光栅投影结构光级数的三维重建方法 | |
CN105300317A (zh) | 基于正弦和三角波条纹投影的三维测量方法 | |
CN112504165A (zh) | 一种基于双边滤波优化的复合立体相位展开方法 | |
CN109751969A (zh) | 一种利用正反格雷码线移光栅的三维非接触扫描方法 | |
CN112611320A (zh) | 一种基于梯度光强的多频外差面结构光三维重建方法 | |
CN112212806B (zh) | 一种基于相位信息导向的立体相位展开方法 | |
CN111947600A (zh) | 基于相位级次代价滤波的鲁棒立体相位展开方法 | |
CN106840039B (zh) | 一种基于复合光栅投影的快速三维测量方法 | |
CN114739322B (zh) | 一种三维测量方法、设备和存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210406 |