[发明专利]一种基于小波变换的改进窗口傅里叶三维测量法

摘要

本发明公开了一种基于小波变换的改进窗口傅里叶三维测量法，其主要目的在于求解出精确的条纹图像的相位分布，并由相位分布得到物体的三维形貌信息。其实现步骤为：将一幅黑白正弦调制条纹图像投影到被测物体上并用CCD进行采集；对采集到的变形条纹图像逐行进行小波变换，提取小波变换脊，然后计算并去除求取脊过程中相位二阶导所带来的误差，最后得到精确的窗口大小矩阵。将精确窗口矩阵代入窗口傅里叶变换，经过滤波等步骤求得变形条纹图的相对相位信息。建立条纹图质量图，然后采用洪水算法相位展开，得到条纹图像的绝对相位分布。根据相位高度转换公式由绝对相位分布获取被测物体的三维信息。

主权项

1.一种基于小波变换的改进窗口傅里叶三维测量法，其特征在于，包括如下步骤：步骤1)，将黑白条纹图像投影到被测物体表面，使用CCD摄像头对被测物体表面进行拍摄，得到一幅高度为c、宽度为r的变形条纹图像g(x，y)：g(x，y)＝A(x，y)+B(x，y)cos[2πf₀x+φ(x，y)]    (1)式其中，A(x，y)是背景光强分布，B(x，y)是物体表面反射率，f₀是正弦条纹频率，φ(x，y)是待求的相对相位分布图，(x，y)表示变形条纹图像的二维坐标；步骤2)，对所述变形条纹图g(x，y)逐行做小波变换，得到变形条纹图像每一点的近似尺度因子，具体过程如下：步骤2.1)，将y视为常数y₁，采用一维小波变换对所述变形条纹图像g(x，y)的第y₁行进行处理，处理过程为：Wy1(a1,b)=∫-∞+∞g(x,y1)M*a,b(x)dx]]>  (2)式其中，y₁表示某一行的行号；a是尺度因子，取值范围为10到90，每隔0.2取一个值；b是平移因子，取值范围为1到条纹的宽度r，每隔1取一个值，单位为像素；经一维小波变换获得的是一个400行r列的二维复数矩阵，称第y₁行小波变换矩阵，a₁是矩阵中元素的行标号；其中，M^*_a，b(x)为M_a，b(x)的共轭函数；根据式(3)的小波函数M(x)得到M^*_a，b(x)表达式如式(4)：M(x)=1(fb2π)1/4exp(2πifcx)exp(-x22fb2)]]>  (3)式M*a,b(x)=1a(fb2π)1/4exp(-2πifcx-ba)exp[-(x-b)22fb2a2];]]>  (4)式其中，f_b是小波函数的带宽，f_c是小波函数的中心频率，i为单位复数，x为自变量；步骤2.2)，求取条纹图像g(x，y)的近似尺度因子分布图f_a1(x，y₁)：计算二维复数矩阵的对应的模矩阵搜索模矩阵第x列中值最大的元素，并将模矩阵的第x列中值最大元素的行标号赋值给a_max，则条纹图像的近似尺度因子分布图f_a1(x，y₁)在坐标(x，y₁)处数值a_rx为a_rx＝10+0.2×a_max；遍历条纹图像所有坐标点，求得条纹图像的近似尺度因子分布图f_a1(x，y)；步骤3)，对经所述步骤2.1)小波变换后的变形条纹图去除小波变换脊误差，得到变形条纹图像每一点的准确尺度因子，具体过程如下：步骤3.1)，去除变形条纹图像每一个坐标点的脊误差，处理过程为：W(x，y₁)＝W₀+W₁+W₂+ε₀    (5)式其中W(x，y₁)为变形条纹图中任意一坐标点的小波变换简化形式，其中ε₀表示由变形条纹图中每一点的相位二阶导所带来的误差，其表达式为：(6)式其中，t为中间计算变量，表示相位二阶导，对近似尺度因子逐行求导求出求出最后求出对应尺度因子a的ε₀，得到ε₀为400行的一维复数数组在坐标(x，y₁)的坐标点上的误差数组ε₀(x，y₁)；W₀，W₁，W₂分别为小波变换计算过程中的简化表达式，形式如下：W0(a,b)=A(x,y)aπ(fb2π)1/4exp(-2fb2fc2π2)]]>  (7)式(8)式(9)式其中，为相位、表示相位导数，W(x，y₁)中逐点减去误差ε₀(x，y₁)，求出精确小波变换脊的值W_ε(x，y₁)＝W₀+W₁+W₂，求出每一行精确小波变换脊矩阵W_400×r，矩阵W_400×r为400行r列的复数矩阵；步骤3.2)，求取条纹图像的精确尺度因子分布图f_a2(x，y)：求出W_400×r的对应的模矩阵C_400×r(a₁，b)，搜索模矩阵C_400×r(a₁，b)第x列中值最大的元素，并将模矩阵C_400×r(a₁，b)的第x列中值最大元素的行标号赋值给a_amax，则条纹图像的近似尺度因子分布图f_a2(x，y₁)在坐标(x，y₁)处数值a_acr为：a_acr＝10+0.2×a_amax；遍历条纹图像所有坐标点，求得条纹图像的精确尺度因子分布图f_a2(x，y)；步骤4)，对变形条纹图逐行做窗口傅里叶变换，求出变形条纹图相对相位分布图，具体过程如下：步骤4.1)，将y视为常数，采用一维窗口傅立叶变换对变形条纹图像g(x，y)的每行进行处理，一维窗口傅立叶变换过程为：WF(b,ξ)=∫-∞+∞g(x)Wδ(x-b)exp(-jξx)dx]]>  (10)式其频域表达式为：WF(fs,y)=Σn=0n=+∞Pn(fs-nf0,y)]]>  (11)式其中，g(x)为一行条纹图像；WF(b，ξ)表示一维窗口傅立叶变换，ξ表示频域计算因子，δ表示一维窗口傅立叶的窗口尺度因子，δ取值为一维窗口的位置(x-b，y)对应的精确尺度因子分布图f_a2(x，y)相应位置上的点的值，W_δ(x-b)表示窗口函数，表达式为：Wδ(x-b)=δ-12πexp[-(b-x)22δ2]]]>  (12)式其中，n表示阶次，依次取值为0，1，2，……，无穷；P_n(f_s-nf₀，y)表示任一点傅里叶变换后对应的n阶频谱；f_s表示频域的变量；步骤4.2)，对傅里叶变换后的频谱滤波并提取相位信息，具体过程如下：对P_n(f_s-nf₀，y)进行滤波并提取出一阶频谱P₁(f_s-f₀，y)，再对P₁(f_s-f₀，y)进行逆傅里叶变换，得到含有相位信息的B(x，y)exp[iφ(x，y)]，计算B(x，y)exp[iφ(x，y)]的角度值即可得到含有物体高度信息的变形条纹图相对相位值φ(x，y)，得到的相位值是介于0-2π之间；遍历条纹图像所有坐标点，得到变形条纹图的相对相位分布图φ_A(x，y)；步骤5)，建立条纹图像质量图，采用洪水算法展开相对相位分布图φ_A(x，y)，得到实际相位图具体过程如下：步骤5.1)，利用相对相位分布图中的相位梯度建立质量图，质量图可以按照以下公式计算：Δ(x,y)=Wrap2{φA(x+1,y)-φA(x,y)}+Wrap2{φA(x,y+1)-φA(x,y)}]]>  (13)式其中W_rap{}是包裹函数，当值大于2π或小于-2π时将其减去或加上2π；步骤5.2)，在相对相位分布图中央找到质量值最高的坐标点，选择该点作为相位展开的起始点，将该点放入一个空的栈中；步骤5.3)，判断栈是否为空；如果是，则相位展开过程结束，得到展开相位结果并进入步骤6)；如果否，弹出栈顶的点，展开该点的四邻点中没有处理的坐标点，并将这些未处理的点入栈；步骤5.4)，按照质量图中的质量值将栈中所有点排序，质量最高的点放在栈顶，转到步骤5.3)继续处理；步骤6)，读取最终的展开相位结果根据经典光栅投影的从展开相位结果到物体高度h(x，y)的转换公式，最终求得测量物体的三维信息；其中转换公式h(x，y)为：(14)式其中，l、d是测量系统的几何参数，l是投影仪到测量平面的距离，d是CCD摄像头到投影仪的距离；表示相位变化量，为展开相位结果，为初始相位结果，ω₀为投影光栅的角频率。