CN106813596B - 一种自标定阴影莫尔三维轮廓测量方法 - Google Patents
一种自标定阴影莫尔三维轮廓测量方法 Download PDFInfo
- Publication number
- CN106813596B CN106813596B CN201710037589.6A CN201710037589A CN106813596B CN 106813596 B CN106813596 B CN 106813596B CN 201710037589 A CN201710037589 A CN 201710037589A CN 106813596 B CN106813596 B CN 106813596B
- Authority
- CN
- China
- Prior art keywords
- phase
- measurement
- bar graph
- striped
- measuring
- 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
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/2504—Calibration devices
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
本发明属于光学精密测量技术领域,具体涉及一种自标定阴影莫尔三维轮廓测量方法,依次包括下述步骤:第一步:将被测物体置于测量空间,开启光源,并使用摄像机记录产生的条纹图;第二步:运用精密位移台以预定义的距离等步长移动光栅两次,并记录每次移动后产生的条纹图;第三步:对采集的条纹图进行预处理:运用条纹降噪的方法去除图像中的高频噪音信号;第四步:对条纹图信号的相位进行解调处理;第五步:运用去包裹方法对上面提取的测量相位进行相位展开;第六步:根据解算的相位和测量灵敏度利用下式进行相位‑高度映射,完成本次测量任务。本发明的测量方法仅需要测量三帧条纹图,节约了测量时间,本发明测量方法的计算过程准确、方便。
Description
技术领域
本发明涉及光学精密测量技术领域,具体涉及一种应用三帧条纹图的相移阴影莫尔三维轮廓测量方法。
背景技术
阴影莫尔三维轮廓测量技术是一种全场、非接触的光学三维轮廓测量方法,自1970年提出以来,已广泛应用于人体测量、晶圆表面检测、BGA共面度测试、自动生产、航空工业等领域。
早期的阴影莫尔技术主要基于莫尔等高线法。莫尔等高线法简单直观,但不能有效利用条纹间距的灰度级信息,其测量灵敏度仅限一个条纹间距。
基于傅立叶的阴影莫尔方法通过倾斜被测物体为测量视场引入载频,然后通过傅里叶变换技术获得测量相位。然而该方法已受衍射影响:若倾斜角较大,受衍射的影响(假定主光栅的周期≥0.1mm),此时产生的莫尔条纹图对比度会随物体和光栅间距的增加而减小。这使得条纹图光强分布不均匀;若倾斜角太小,引入的载频过小,这样在一些实验中,又不能有效进行不同频率信号的分离。
相比之下,相移法阴影莫尔技术测量方法简单,且是一种自动的解调方法。基本原理是:通过在时域改变实验参数,引起测量光场的相位变化,同时使用摄像机记录并存储该变化,进而运用计算机程序逐点地分析,实现测量相位自动解调。相移技术可以充分利用条纹图中的灰度信息,提高了测量分辨力,其过程主要包括相位解调、相位去包裹、测量灵敏度标定和相位-高度映射等。
传统上,测量灵敏度通常依赖已知高度分布标定块进行标定。该过程需要制作标定块,因此费时费力。Dirckx等首次提出了一种仅依赖莫尔测量装置自身的标定方法。在Dirckx方法中,通过移动物体,使得测量视场产生时域光强变化,进而通过分析某一观测点光强与物体移动距离曲线的周期获得即测量灵敏度。Dirckx方法简单、易行,但周期确定麻烦,而且随着应用领域对阴影莫尔技术测量精度要求的提高,传统的阴影技术需要发展成为高灵敏度阴影莫尔技术。在其系统中,细光栅(光栅周期范围为:0.2mm:0.025mm)得到了应用。此时,测量装置的动态范围受泰伯效应影响明显减小。因此应用Dirckx方法时,物体极易移出测量装置的动态范围。最近,提出了一种基于主量分析法的阴影莫尔测量灵敏度标定方法能方便的确定光强变化周期。然而该方法需要相位解包裹,而且在应用中要求引入的相移量分布在[0,2π]范围内。因此缺乏灵活性。
在相位解调方面,受阴影莫尔技术相位-高度非线性影响,相移阴影莫尔视场无法引入一个全局的相移。如果应用传统的三步相移技术,但需要忽略相位-高度非线性关系,近似实现全场相移。从而致使相移技术无法得到精确解。
近年来,基于较少帧条纹图的相移算法得到了较大的发展。其中Vargas等提出的克莱姆正交化二帧算法,解调速度快,测量精度得到了广泛的应用。然而该方法的精度受条纹图背景去除精度影响较大。最近Wang等提出了一种正则化条纹图差的方法对其进行了改进,并通过实验证明了该方法的有效性。须提及的是,尽管Wang的算法可以减少条纹图背景对解调结果的影响,但是该方法的解调结果相对测量相位有一个位移。因此在高度测量场合产生了问题。
发明内容
一种自标定阴影莫尔三维轮廓测量方法,以改善现有的相移阴影莫尔物体三维轮廓光学测量装置工件测量的精度,节约测量时间,计算过程准确、方便。
为了达到上述目的,本发明的解决方案是:
一种自标定阴影莫尔三维轮廓测量方法,其特征在于:依次包括下述步骤:
第一步:将被测物体置于测量空间,开启光源,并使用摄像机记录产生的条纹图;
第二步:运用精密位移台以预定义的距离等步长移动光栅两次,并记录每次移动后产生的条纹图;
第三步:对采集的条纹图进行预处理:运用条纹降噪的方法去除图像中的高频噪音信号;
第四步:对条纹图信号的相位进行解调处理:
设A(x,y)为背景,B(x,y)为调制项,φ(x,y)为相位,δ为相移,n为时域条纹序号,则摄取的条纹图光强分别可表述为:
In(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+nδ],(n=0,1,2) (1)
为了表述清晰,下面的推导过程省略坐标项(x,y)。
首先应用运用条纹互减去除背景项,得到的新条纹表述为:
进一步对新条纹进行加、减运算得:
显然(3)式振幅不等,因此对其进行正则化,并约掉常数项得:
式中代表了矩阵范数运算,其中,m、n为条纹图大小。采集的条纹图中条纹数目大于1,根据三角函数积分计算特点知,上述矩阵范数近似相等,因此可得:
进而,得:
θ=φ+δ/2=arctan(IN1/IN0) (6)
将(6)式带入(3)式,得:
δ=arctan[Is·tan(IN1/IN0)/Ia] (7)
因此得测量相位为:
φ=θ-δ (8)
上式中参数分别表示:l为光源-CCD平面与参考光栅的距离、d为光源与CCD的距离、Δl为相移时光栅移动距离和p光栅常数;
第五步:运用去包裹方法对上面提取的测量相位进行相位展开;
第六步:根据解算的相位和测量灵敏度利用下式进行相位-高度映射,完成本次测量任务。
设测量系统灵敏度为s,根据阴影莫尔相位-高度映射关系得:
上述测量方法的步骤2中所说的预定距离为0.001mm。
本发明提出一种基于自标定、正交相位解调的阴影莫尔三维传感技术。该方法通过电动位移台移动光栅两次,获得三帧相移条纹图,运用条纹间差运算去除条纹图背景,产生两帧新条纹图,然后加、减运算的方法产生的新条纹,并利用克莱姆正交化方法对新条纹进行规范化处理,进而解算出相移。与现有技术相比,本发明的优点是:
1、传统方法需要采集较多帧(如,16帧)相移条纹图,然后根据给定点光强变化的周期确定测量灵敏度,其过程费时费力,本专利的方法仅需要三帧条纹图,因此节约了测量时间。
2、本专利方法中灵敏度的确定采用解析的方法,灵敏度确定过程采用具体的公式计算,而不是传统标定方法中设计程序寻找特征点的数值分析方法,因此计算过程准确、方便。
3、本专利的方法采用分析相移条纹图的方法标定测量灵敏度,因此可在相位解调过程中同时解算测量灵敏度,无须额外进行灵敏度标定工作及设计相移装置,具有自标定的特点。
4、该方法将为精密零件表面轮廓三维快速表征开拓新途径,有望实现IC电路硅基底表面平整度检测,可以大大简化相移阴影莫尔三维轮廓测量的技术水平。
5、应用要求低:该方法不需要光源的一致性,对环境要求低,适用于工业场合的测量需要。可以自动、高精度完成被测物体表面形貌测量,为实现工况下关键零部件变形检测等打下基础。
附图说明:
图1是本发明测量装置结构图;
图2是本发明测量原理图;
图3是本发明进行测量时的虚拟仪器界面示意图,其中:图3(a),摄像机实时图;图3(b),记录的相移条纹图;图3(c),包裹相图;图3(d),去包裹相图;图3(e),测量结果;图3(f),截取部分测量结果;
图中,1-线光源,2-CCD摄像机,3-电动位移台,4-光栅,5-被测物品。
具体实施方式
实施例1:
参见图1和图2,本发明所用装置为常规的装置,包括CCD摄像机2、线光源1、光栅4和电动位移台3,CCD摄像机2的光轴垂直于光栅面,CCD摄像机2和线光源1的光源中心置于距离光栅4所在平面L的距离,线光源1与CCD摄像机1中心距离为d。
测量时,被测物品5设置于测量位置上。一种自标定阴影莫尔三维轮廓测量方法,依次包括下述步骤:
第一步:将被测物体置于测量空间,开启光源,并使用摄像机记录产生的条纹图,具体地说,就是将工件置于测量空间,通过虚拟仪器界面进行初始化,如调节摄像机的亮度、饱和度等,同时开启光源照射测量光栅,在物体表面产生条纹图,然后控制摄像机摄取并存储光场变化,便得到了用于解调物体表面形貌的相移条纹图;
第二步:运用精密位移台以预定义的距离等步长移动光栅两次,并记录每次移动后产生的相移条纹图,本实施例中所采用的预定义距离为0.001mm;
第三步:对采集的条纹图进行预处理:运用条纹降噪的方法去除图像中的高频噪音信号;
第四步:对条纹图信号的相位进行解调处理:
假设A(x,y)为背景,B(x,y)为调制项,φ(x,y)为相位,δ为相移,n为时域条纹序号,则摄取的条纹图光强分别可表述为:
In(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+nδ],(n=0,1,2) (1)
为了表述清晰,下面的推导过程省略坐标项(x,y)。
首先应用运用条纹互减去除背景项,得到的新条纹表述为:
进一步对新条纹进行加、减运算得:
显然(3)式振幅不等,因此对其进行正则化,并约掉常数项得:
式中代表了矩阵范数运算,其中,m、n为条纹图大小。一般地,采集的条纹图中条纹数目大于1,根据三角函数积分计算特点知,上述矩阵范数近似相等。因此可得:
进而,得:
θ=φ+δ/2=arctan(IN1/IN0) (6)
将(6)式带入(3)式,得:
δ=arctan[Is·tan(IN1/IN0)/Ia] (7)
因此得测量相位为:
设测量系统灵敏度为s,根据阴影莫尔相位-高度映射关系得:
上式中参数分别表示:l为光源-CCD平面与参考光栅的距离、d为光源与CCD的距离、Δl为相移时光栅移动距离和p光栅常数。;
第五步:运用去包裹方法对上面提取的测量相位进行相位展开;
第六步:根据解算的相位和测量灵敏度进行相位-高度映射,完成本次测量任务。
本发明的方法可采用Labview平台将光栅位置调整、光源控制、条纹图产生、采集和数据处理等环节集成在虚拟仪器环境中,构建阴影莫尔三维轮廓虚拟测量仪器。整个测量过程集成起来,可有效实现测量过程控制与操作自动化、便捷化。
测量时,首先将物体置于测量空间,先通过虚拟仪器界面进行初始化,如调节摄像机的亮度、饱和度等,然后开启单色光源照射测量光栅,物体表面产生变形光栅,并与测量光栅交叠产生莫尔条纹图。通过移动光栅引起测量视场莫尔条纹图的相位变化,然后运用摄像机记录两次视场条纹图变化,解调时应用相移技术提取测量视场的相位及引入的相移,进而根据相移解算测量灵敏度,将绝对相位映射至高度空间。
测量系统设计的自动工作流程为:
开始→开启光源→摄像机初始化→采集初始条纹图→记录条纹图→移动光栅产生相移→记录条纹图→再次移动光栅相移→记录条纹图→条纹预处理→信息处理→输出结果。
虚拟仪器界面参见图3,图中:(a),摄像机实时图;(b),记录的相移条纹图;(c),包裹相图;(d),去包裹相图;(e),测量结果;(f),截取部分测量结果。
Claims (2)
1.一种自标定阴影莫尔三维轮廓测量方法,其特征在于:依次包括下述步骤:
第一步:将被测物体置于测量空间,开启光源,并使用摄像机记录产生的条纹图;
第二步:运用精密位移台以预定义的距离等步长移动光栅两次,并记录每次移动后产生的条纹图;
第三步:对采集的条纹图进行预处理:运用条纹降噪的方法去除图像中的高频噪音信号;
第四步:对条纹图信号的相位进行解调处理:
设A(x,y)为背景,B(x,y)为调制项,φ(x,y)为相位,δ为相移,n为时域条纹序号,则摄取的条纹图光强分别可表述为:
In(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+nδ],(n=0,1,2) (1)
为了表述清晰,下面的推导过程省略坐标项(x,y);
首先应用运用条纹互减去除背景项,得到的新条纹表述为:
进一步对新条纹进行加、减运算得:
显然(3)式振幅不等,因此对其进行正则化,并约掉常数项得:
式中代表了矩阵范数运算,其中,m、n为条纹图大小;采集的条纹图中条纹数目大于1,根据三角函数积分计算特点知,上述矩阵范数近似相等,因此可得:
进而,得:
θ=φ+δ/2=arctan(IN1/IN0) (6)
将(6)式带入(3)式,得:
δ=arctan[Is·tan(IN1/IN0)/Ia] (7)
因此得测量相位为:
φ=θ-δ (8)
第五步:运用去包裹方法对第四步得到的测量相位进行相位展开;
第六步:根据解算的相位和测量灵敏度利用下式进行相位-高度映射,完成本次测量任务;
设测量系统灵敏度为s,根据阴影莫尔相位-高度映射关系得:
上式中参数分别表示:l为光源-CCD平面与参考光栅的距离、d为光源与CCD的距离、△l为相移时光栅移动距离和p光栅常数。
2.根据权利要求1所述的自标定阴影莫尔三维轮廓测量方法,其特征在于:所述测量方法的第二步中所说的预定距离为0.001mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710037589.6A CN106813596B (zh) | 2017-01-18 | 2017-01-18 | 一种自标定阴影莫尔三维轮廓测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710037589.6A CN106813596B (zh) | 2017-01-18 | 2017-01-18 | 一种自标定阴影莫尔三维轮廓测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106813596A CN106813596A (zh) | 2017-06-09 |
CN106813596B true CN106813596B (zh) | 2019-01-18 |
Family
ID=59112297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710037589.6A Active CN106813596B (zh) | 2017-01-18 | 2017-01-18 | 一种自标定阴影莫尔三维轮廓测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106813596B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107727002B (zh) * | 2017-09-29 | 2020-03-13 | 西安工业大学 | 阴影莫尔参数标定用结构、标定方法及失调误差消除方法 |
CN108254086B (zh) * | 2017-12-29 | 2022-01-28 | 西安工业大学 | 一种随机光纤点衍射测量方法 |
CN108775874A (zh) * | 2018-06-04 | 2018-11-09 | 西安工业大学 | 一种应用随机相移技术的阴影三维轮廓测量装置和方法 |
CN109341584A (zh) * | 2018-11-02 | 2019-02-15 | 西安工业大学 | 一种齿轮齿面三维形貌表征方法 |
CN113409367B (zh) * | 2021-07-08 | 2023-08-18 | 西安交通大学 | 条纹投影测量点云逐点加权配准方法、设备和介质 |
CN114018176A (zh) * | 2021-10-27 | 2022-02-08 | 华中科技大学 | 一种投影图像处理模块、三维重构方法及其系统 |
CN114623782B (zh) * | 2022-03-16 | 2024-07-09 | 合肥因赛途科技有限公司 | 一种无低通滤波的数字莫尔三维测量方法 |
CN114858094B (zh) * | 2022-04-02 | 2023-11-10 | 西安工业大学 | 一种针对动态物体的学习型单帧莫尔轮廓测量方法 |
CN115790452B (zh) * | 2023-02-06 | 2023-05-16 | 西安知象光电科技有限公司 | 一种齿轮齿面三维形貌莫尔表征与测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1206845A (zh) * | 1998-06-22 | 1999-02-03 | 西安交通大学 | 变光源阴影莫尔测量方法及其装置 |
CN1900651A (zh) * | 2006-07-27 | 2007-01-24 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
US20070086020A1 (en) * | 2005-10-19 | 2007-04-19 | The University Of Maryland | Shadow moire using non-zero talbot distance |
US20080239288A1 (en) * | 2007-04-02 | 2008-10-02 | Korea Advanced Institute Of Science And Technology | 3D shape measurement apparatus and method using stereo moire technique |
CN105157616A (zh) * | 2015-07-31 | 2015-12-16 | 西安工业大学 | 一种阴影莫尔轮廓测量装置、其标定方法和测量方法 |
-
2017
- 2017-01-18 CN CN201710037589.6A patent/CN106813596B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1206845A (zh) * | 1998-06-22 | 1999-02-03 | 西安交通大学 | 变光源阴影莫尔测量方法及其装置 |
US20070086020A1 (en) * | 2005-10-19 | 2007-04-19 | The University Of Maryland | Shadow moire using non-zero talbot distance |
CN1900651A (zh) * | 2006-07-27 | 2007-01-24 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
US20080239288A1 (en) * | 2007-04-02 | 2008-10-02 | Korea Advanced Institute Of Science And Technology | 3D shape measurement apparatus and method using stereo moire technique |
CN105157616A (zh) * | 2015-07-31 | 2015-12-16 | 西安工业大学 | 一种阴影莫尔轮廓测量装置、其标定方法和测量方法 |
Non-Patent Citations (4)
Title |
---|
Shadow moiré technology based fast method for the measurement of surface topography;Hong Zhao等;《APPLIED OPTICS》;20131231;第52卷(第33期);第7874-7881页 |
Three-frame self-calibration phase shift algorithm using the Gram–Schmidt orthonormalization approach;HUBING DU等;《Applied Optics》;20160831;第55卷(第24期);第6713-6717页 |
一种多距离融合的大深度测量范围相移阴影叠栅轮;艾永旭等;《光学学报》;20161231;第36卷(第12期);第1212001-1-1212001-7页 |
基于正则化主量分析的相移阴影叠栅解调技术;杜虎兵;《激光与光电子学进展》;20161231;第100501-1-100501-7 |
Also Published As
Publication number | Publication date |
---|---|
CN106813596A (zh) | 2017-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106813596B (zh) | 一种自标定阴影莫尔三维轮廓测量方法 | |
Su et al. | Reliability-guided phase unwrapping algorithm: a review | |
Zappa et al. | Static and dynamic features of Fourier transform profilometry: A review | |
CN106092158B (zh) | 物理参数估计方法、装置和电子设备 | |
CN102183214A (zh) | 一种大口径非球面镜结构光检测方法 | |
CN108426539A (zh) | 一种基于双频正交光栅投影检测反射元件面形的方法 | |
CN111397634B (zh) | 星敏感器固定端面热变形的高分辨干涉检测装置及方法 | |
CN113175894A (zh) | 一种物体表面三维形貌白光干涉测量装置及方法 | |
CN103344209B (zh) | 一种反射镜零重力面形测试方法 | |
CN105157616A (zh) | 一种阴影莫尔轮廓测量装置、其标定方法和测量方法 | |
CN107192337A (zh) | 基于单缝衍射原理利用ccd测量微小位移的方法 | |
Martínez et al. | Industrial calibration procedure for confocal microscopes | |
WO2023124867A1 (zh) | 基于伪Wigner-Ville分布的白光扫描干涉三维重建方法 | |
Dvoynishnikov et al. | Optoelectronic method of contactless reconstruction of the surface profile of complexly shaped three-dimensional objects | |
CN204854655U (zh) | 一种阴影莫尔轮廓测量装置 | |
CN107560638B (zh) | 一种星敏感器变温标定参数的获取方法 | |
CN108489422B (zh) | 一种变频相移最小二乘迭代叠加面形分离的方法 | |
Patiño et al. | A Gaussian filtering method to reduce directionality on high-density point clouds digitized by a conoscopic holography sensor | |
Zhu et al. | A method for measuring the guideway straightness error based on polarized interference principle | |
Macher et al. | Detection of visually perceptible sink marks on high gloss injection molded parts by phase measuring deflectometry | |
Yu | Image processing methods based on physical models | |
Somthong et al. | Average surface roughness evaluation using 3-source photometric stereo technique | |
Du et al. | Calibration of the high sensitivity shadow moiré system using random phase-shifting technique | |
CN113091648A (zh) | 一种三维空间信息采集分析处理的方法及装置 | |
Lei et al. | Vibration testing parameters measured by sampling moire method |
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 |