本发明涉及一种基于单幅光栅投影的三维测量方法和系统,属于三维测量领域。
背景技术:
takeda等人1983年提出了基于光栅投影的三维面形测量方法——傅里叶变换轮廓术(ftp)。傅里叶变换轮廓术是通过向物体表面投射光栅条纹,获得被物体高度调制而发生变形的光栅,利用特定算法对变形条纹图像进行傅里叶变换、滤波、傅里叶逆变换和解包裹处理,提取其中的相位,再通过对测量系统的标定,从而获得物体的三维信息。
近年来,s变换已经应用到光学三维面形测量中光栅投影的相位解调过程,同ftp比较,s变换技术在具备ftp的优点之上,处理速度和频谱处理上有更多优势。
解相位是光栅投影法的基本问题之一。首先通过s变换获得条纹图的相位主值,值域位于(-π,+π]区间,再将相位主值恢复成完整的相位场。目前采用的解包裹方法如洪水法、质量图引导法等,过程都较为复杂,耗时较长,且都需两幅条纹图才能获得绝对相位,
如何快速、简单地进行解包裹,真正实现动态物体测量,满足自动化需求在主动光学三维测量中非常重要。
技术实现要素:
为了解决上述问题,本发明通过提供一种基于单幅光栅投影的三维测量方法和系统。
本发明采用的技术方案一方面为一种基于单幅光栅投影的三维测量方法,包括:
物理光栅离焦投影到待测试物体以形成正弦条纹图,采集待测试物体的变形条纹图;
利用s变换法处理变形条纹图以获得主值相位φ(x,y);
基于消除包裹法进行解包裹,获得绝对相位
对基于远心成像的三维测量系统进行标定以建立成像模型,计算被测物体的三维坐标信息。
优选地,利用s变换法处理变形条纹图的步骤包括:
对获得的变形正弦条纹图h(t)进行一维s变换,s变换系数s(τ,f)公式为:
对获得的包含s变换的复杂矩阵采用平顶汉宁窗加权滤波;
对滤波获得的局部基频分量沿时间轴进行叠加,获得完整基频分量,对其进行傅里叶逆变换后,获得基频复信号,表示为
优选地,所述平顶汉宁窗的数学表达式为:
该平顶汉宁窗以s变换幅值最大处为中心,标记该处为s变换脊,fb为s变换脊的频率,fwlow,fwhigh分别表示从s变换脊分别往低频与高频方向上的延伸宽度,fb+fwhigh,fb-fwlow分别表示高低端截止频率,else表示其他频率。
优选地,利用消除包裹法进行解包裹的步骤包括:
利用欧拉公式将包裹相位φ(x,y)变成复数形式,ejφ(x,y)=cosφ(x,y)+jsinφ(x,y);
令φc(x,y)=ejφ(x,y),取φc(x,y)的第x行,x∈[0,m-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcx(x,y)进行一维傅里叶变换,求得相位的水平方向频谱偏移量μ0;
令φc(x,y)=ejφ(x,y),取φc(x,y)的第y列,y∈[0,n-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcy(x,y)进行一维傅里叶变换,得到相位的竖直方向频谱偏移量ν0;
根据傅里叶变换频移特性,在空间域将相位的频谱移动到原始位置,所述傅里叶变换频移特性表达式为:
其中,(t,z)表示空间变量,(μ,ν)表示频率域变量,则有,
对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位。
优选地,对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位,所述四象限反正切操作的表达式:
本发明采用的技术方案另一方面为一种基于单幅光栅投影的三维测量系统,包括:
光栅模块,用于将物理光栅离焦投影到待测试物体以形成正弦条纹图,采集待测试物体的变形条纹图;
处理模块,用于利用s变换法处理变形条纹图以获得主值相位φ(x,y);
还用于基于消除包裹法进行解包裹,获得绝对相位
测量模块,用于对基于远心成像的三维测量系统进行标定以建立成像模型,计算被测物体的三维坐标信息。
优选地,利用s变换法处理变形条纹图的步骤包括:
对获得的变形正弦条纹图h(t)进行一维s变换,s变换系数s(τ,f)公式为:
对获得的包含s变换的复杂矩阵采用平顶汉宁窗加权滤波;
对滤波获得的局部基频分量沿时间轴进行叠加,获得完整基频分量,对其进行傅里叶逆变换后,获得基频复信号,表示为
优选地,所述平顶汉宁窗的数学表达式为:
该平顶汉宁窗以s变换幅值最大处为中心,标记该处为s变换脊,fb为s变换脊的频率,fwlow,fwhigh分别表示从s变换脊分别往低频与高频方向上的延伸宽度,fb+fwhigh,fb-fwlow分别表示高低端截止频率,else表示其他频率。
优选地,利用消除包裹法进行解包裹的步骤包括:
利用欧拉公式将包裹相位φ(x,y)变成复数形式,ejφ(x,y)=cosφ(x,y)+jsinφ(x,y);
令φc(x,y)=ejφ(x,y),取φc(x,y)的第x行,x∈[0,m-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcx(x,y)进行一维傅里叶变换,求得相位的水平方向频谱偏移量μ0;
令φc(x,y)=ejφ(x,y),取φc(x,y)的第y列,y∈[0,n-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcy(x,y)进行一维傅里叶变换,得到相位的竖直方向频谱偏移量ν0;
根据傅里叶变换频移特性,在空间域将相位的频谱移动到原始位置,所述傅里叶变换频移特性表达式为:
其中,(t,z)表示空间变量,(μ,ν)表示频率域变量,则有,
对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位。
优选地,对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位,所述四象限反正切操作的表达式:
本发明的有益效果为基于单幅条纹图,在s变换法得到主值相位后,对该主值相位两端进行补零,使得频域内获得上采样,根据傅里叶变换频移特性将相位的频谱移动到原始位置,消除载波的相位分量,采用远心成像模型,拟合高度-相位映射关系,对三维测量系统进行标定,获取被测物体的三维坐标信息,步骤简单,测量速度更快。
附图说明
图1所示为基于本发明实施例的一种基于单幅光栅投影的三维测量方法示意图;
图2所示为基于本发明实施例的s变换法获得的主值相位图像;
图3所示为基于本发明实施例的对主值相位图某一行补零后的结果图;
图4所示为基于本发明实施例的硬币的三维轮廓图像。
具体实施方式
以下结合实施例对本发明进行说明。
基于发明的实施例1,如图1所示一种基于单幅光栅投影的三维测量方法,包括:物理光栅离焦投影到待测试物体以形成正弦条纹图,采集待测试物体的变形条纹图;利用s变换法处理变形条纹图以获得主值相位φ(x,y);基于消除包裹法进行解包裹,获得绝对相位
基于发明的实施例1所述的方法,利用s变换法处理变形条纹图的步骤包括:
对获得的变形正弦条纹图h(t)进行一维s变换,s变换系数s(τ,f)公式为:
对滤波获得的局部基频分量沿时间轴进行叠加,获得完整基频分量,对其进行傅里叶逆变换后,获得基频复信号,表示为
基于发明的实施例1所述的方法,所述平顶汉宁窗的数学表达式为:
该平顶汉宁窗以s变换幅值最大处为中心,标记该处为s变换脊,fb为s变换脊的频率,fwlow,fwhigh分别表示从s变换脊分别往低频与高频方向上的延伸宽度,fb+fwhigh,fb-fwlow分别表示高低端截止频率,else表示其他频率。
基于发明的实施例1所述的方法,利用消除包裹法进行解包裹的步骤包括:
利用欧拉公式将包裹相位φ(x,y)变成复数形式,ejφ(x,y)=cosφ(x,y)+jsinφ(x,y);
令φc(x,y)=ejφ(x,y),取φc(x,y)的第x行,x∈[0,m-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的1*2kn大小的矩阵φcx(x,y)进行一维傅里叶变换,求得相位的水平方向频谱偏移量μ0;
令φc(x,y)=ejφ(x,y),取φc(x,y)的第y列,y∈[0,n-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的2km*1大小的矩阵φcy(x,y)进行一维傅里叶变换,得到相位的竖直方向频谱偏移量ν0;
根据傅里叶变换频移特性,在空间域将相位的频谱移动到原始位置,所述傅里叶变换频移特性表达式为:
其中,(t,z)表示空间变量,(μ,ν)表示频率域变量,则有,
对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位。
基于发明的实施例1所述的方法,对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位,所述四象限反正切操作的表达式:
基于远心成像的三维形貌测量系统标定方法为现有技术,包括步骤:
步骤s1:搭建远心三维形貌测量系统,所述测量系统包括:远心投影设备、远心摄像设备、平移台;远心摄像设备的光轴垂直于水平放置的平移台,远心投影设备的光轴与平移台成一夹角,并控制远心摄像设备的光轴与远心投影设备的光轴处于同一平面内;
步骤s2:使平移台处于远心投影设备和远心摄像设备的共同景深范围内,控制远心投影设备对平移台投射正弦条纹图,远心摄像设备采集所述正弦条纹图,选取远心摄像设备图像平面上的任一像素点作为标定用像素点,利用多步相移法解得该标定用像素点的绝对相位值,并记录此时的平移台高度值;
控制平移台在远心投影设备和远心摄像设备的共同景深范围内沿远心摄像设备的光轴方向进行数次位移,在移动平移台至不同高度时,获取此高度下标定用像素点的绝对相位值,并记录相应的平移台高度值;
对获得的平移台的高度值与相应的标定用像素点的绝对相位值进行线性拟合,建立远心成像的三维形貌测量系统中绝对相位值与平移台高度值的转换关系;
步骤s3:通过标定远心摄像设备的参量,将远心摄像设备图像平面上的像素坐标转换为世界坐标。
基于本发明的实施例2,一种基于单幅光栅投影的三维测量系统,包括:
光栅模块,用于将物理光栅离焦投影到待测试物体以形成正弦条纹图,采集待测试物体的变形条纹图;
处理模块,用于利用s变换法处理变形条纹图以获得主值相位φ(x,y);
还用于基于消除包裹法进行解包裹,获得绝对相位
测量模块,用于对基于远心成像的三维测量系统进行标定以建立成像模型,计算被测物体的三维坐标信息。
基于本发明的实施例2所述的系统,利用s变换法处理变形条纹图的步骤包括:
对获得的变形正弦条纹图h(t)进行一维s变换,s变换系数s(τ,f)公式为:
对获得的包含s变换的复杂矩阵采用平顶汉宁窗加权滤波;
对滤波获得的局部基频分量沿时间轴进行叠加,获得完整基频分量,对其进行傅里叶逆变换后,获得基频复信号,表示为
基于本发明的实施例2所述的系统,所述平顶汉宁窗的数学表达式为:
该平顶汉宁窗以s变换幅值最大处为中心,标记该处为s变换脊,fb为s变换脊的频率,fwlow,fwhigh分别表示从s变换脊分别往低频与高频方向上的延伸宽度,fb+fwhigh,fb-fwlow分别表示高低端截止频率,else表示其他频率。
基于本发明的实施例2所述的系统,利用消除包裹法进行解包裹的步骤包括:
利用欧拉公式将包裹相位φ(x,y)变成复数形式,ejφ(x,y)=cosφ(x,y)+jsinφ(x,y);
令φc(x,y)=ejφ(x,y),取φc(x,y)的第x行,x∈[0,m-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcx(x,y)进行一维傅里叶变换,求得相位的水平方向频谱偏移量μ0;
令φc(x,y)=ejφ(x,y),取φc(x,y)的第y列,y∈[0,n-1],对其两端进行补零,所述补零的表达式为:
其中,m表示沿水平方向φc(x,y)的像素值大小,n表示沿竖直方向φc(x,y)的像素值大小,k为一个整数,对获得的矩阵φcy(x,y)进行一维傅里叶变换,得到相位的竖直方向频谱偏移量ν0;
根据傅里叶变换频移特性,在空间域将相位的频谱移动到原始位置,所述傅里叶变换频移特性表达式为:
其中,(t,z)表示空间变量,(μ,ν)表示频率域变量,则有,
对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位。
基于本发明的实施例2所述的系统,对φcs(x,y)进行四象限反正切操作,获得消除了包裹的相位,所述四象限反正切操作的表达式:
基于本发明的实施例3,计算出被测物体的三维坐标信息的方法。在windows操作系统下采用matlab编程工具上实现。本实施例采用5角硬币作为被测物体,最终得到物体的绝对相位分布,并生成三维图像。
本发明中首先调整投影仪和相机的位置,使待测物体处于二者共同景深下,调整使得条纹清晰度、正弦性良好时,进行拍摄。
图2为s变换法获得的主值相位图像,如图3所示为对主值相位图某一行补零后的结果图,在s变换法得到主值相位后,对该主值相位两端进行补零,使得频域内获得上采样,根据傅里叶变换频移特性将相位的频谱移动到原始位置,消除载波的相位分量,得到的即是所求受被测物体高度调制的相位;最后采用远心成像模型,拟合高度-相位映射关系,对三维测量系统进行标定,求得被测物体的三维坐标信息,如图4所示为硬币的三维轮廓图像。具体处理过程如下:
首先将摄像机像面像素坐标转换为物理世界的二维面内坐标。对于摄像机图像上任意一点(u,v),a为摄像机的内参数矩阵,r和t分别表示摄像机坐标系相对于世界坐标系的旋转矩阵和平移向量,[r,t]为外参矩阵,根据物相双远心光路的成像特点,则参考面(z=0)上某点的世界坐标与计算机图像坐标的转换关系为:
采用现有技术如径向约束标定方法等可获得内外参数矩阵,即可获得二维面内的物理世界坐标。
通过精密平移台在z轴方向得到一系列的标准高度,分别为h1,h2,······,hn;对于摄像机图像上任意一点(u,v),利用本方案的方法计算对应的连续相位分布,分别为
······,
对获得的一系列标准高度h1,h2,······,hn,与获得的对应的δφ1(u,v),δφ2(u,v),······,δφn(u,v),进行线性差值拟合,得到高度-相位的映射关系;最后计算出被测物体的三维坐标信息。
以上所述,只是本发明的较佳实施例而已,本发明并不局限于上述实施方式,只要其以相同的手段达到本发明的技术效果,都应属于本发明的保护范围。在本发明的保护范围内其技术方案和/或实施方式可以有各种不同的修改和变化。