技术领域
[0001] 本
发明属于三维信息重构的领域。基于黑白正弦光栅投影，采用小波变换法分析
变形条纹图像求解相对
相位，建立
质量图指导相位展开，得到精确绝对相位的过程。
背景技术
[0002] 三维测量技术能够描述物体的三维特征和获取物体表面的三维信息，在产品检测和加工控制、医疗领域、文物保护领域、航空航天领域、文化领域等有着广泛的应用。
[0003] 光栅投影法是一种重要的三维测量技术，通过向物体表面投射正弦光栅，将物体的高度信息以相位的形式隐含在光栅中，利用CCD获得物体表面的光栅条纹图像，并使用特定
算法对条纹图像进行处理，提取其中的相位，从而建立物体的三维信息。常用的求解条纹图像相位的方法有
相移法、傅里叶变换法、窗口傅里叶变换法、小波变换法等。
[0004] 傅里叶变换法由于只需采集一幅条纹图像即可完成相位的测量，可以实现动态测量，应用非常广泛。但是傅里叶变换是一种全局的
信号分析工具，不能提取局部信号的特征，在变换中存在的
频谱混叠问题影响着相位测量的
精度。近年来，随着小波分析理论的迅速发展，小波变换被引入光学三维测量领域，使用小波变换分析条纹图像进而实现物体的三维测量，即小波变换轮廓术。小波分析是一种分析非平稳信号的有效工具，相比于已经在
信号处理领域得到广泛应用的传统傅里叶变换法，小波变换在分析信号的局部特征时更有优势。小波变换具有多
分辨率分析的特性，不仅能够像傅里叶变换那样得到信号的整体特征，而且能够对局部信号的细节进行分析，具有较好的时间/空间局部性。所以采用小波变换求解条纹图像的相位避免了傅里叶变换中出现的频谱混叠问题，测量精度更高。
[0005] 小波变换法求解得到的相位值都是介于0-2π的，所以需要相位展开过程。为了实现动态测量，一般不采用添加辅助图像增加信息协助相位展开。简单的扫描线展开相位方法运算速度非常快，但鲁棒性不高，容易出现错误传递现象。质量图指导相位展开算法通过对条纹图像建立一幅反映每一
像素点可靠性的质量图选取一条最佳的相位展开路径，这种方法鲁棒性较高，能够相对精确地展开相位，但运算时间较长，不适用于实时测量。
[0006] 质量指导相位展开算法的关键步骤是质量图的建立。质量图主要有以下几种建立方法。相位梯度法，每个点与四邻的相位梯度的最大值作为每个点的质量值，质量值越大，说明该点的质量越差。表面向量质量图法，使用每一点的法向量和CCD投影方向的负单位向量的内积作为该点的质量值，法向量可由该点和四邻点的相位值计算得出。以上两种方法没有涉及到小波变换的参数，并不适用于小波变换轮廓术。小波变换脊幅值法，根据小波变换脊处的模值表示局部信号和小波函数相似程度的原理，利用小波变换脊的模值建立质量图，指导初始相位的展开，这种方法充分利用了小波变换得到的矩阵，但这种质量图建立方法没有考虑小波变换脊处的尺度因子，当脊处的尺度因子过大或者过小时，模值不能准确说明局部信号的质量好坏；并且当局部信号的幅值变化时，小波变换矩阵中小波脊处的幅值也会变化，故不能准确描述局部信号的正弦性和可靠性。
发明内容
[0007] 技术问题：针对小波变换轮廓术中相位展开的精确性和实时性问题，本发明的目的在于提高小波变换相位展开过程的速度和求解绝对相位的精度。本方法充分利用了小波变换矩阵中的有用信息，结合了直接相位展开速度快和洪
水相位展开法精度高的优点，只需投影一幅条纹图像，实现对较复杂物体表面三维信息的获取，具有较高的精确性，同时实时性较好。
[0008] 技术方案：一种基于小波变换的三维测量方法，具体步骤如下：
[0009] 步骤1：将黑白条纹图像投影到被测物体表面，使用CCD对被测物体表面进行拍摄，得到一幅高度为c、宽度为r的变形条纹图像g(x，y)：
[0010] g(x，y)＝A(x，y)+B(x，y)cos[2πf0x+φ(x，y)]
[0011] 其中，A(x，y)是背景光强分布，B(x，y)是物体表面反射率，f0是正弦条纹
频率，φ(x，y)是待求的相对相位分布，(x，y)表示变形条纹图像的二维坐标，[0012] 步骤2：对变形条纹图像逐行做小波变换，得到变形条纹图像的相对相位分布图，具体过程如下：
[0013] 步骤2.1：将y视为常数，采用一维连续小波变换对变形条纹图像g(x，y)的第y行进行处理，处理过程为：
[0014]
[0015] 其中，a是尺度因子，取值范围为10到50，每隔0.2取一个值，b是平移因子，取值范围为1到条纹图像宽度r，每隔1取一个值，单位为像素，获得的W(a1，b)是一个200行r列的二维复数矩阵，a1是矩阵W(a1，b)中元素的行标号，称第y行小波变换矩阵，M*a，b(x)是Ma，b(x)的共轭函数，M(x)是小波函数，表达式如下：
[0016]
[0017] 其中，fb是小波函数的带宽，fc是小波函数的中心频率，i为单位复数，[0018] 步骤2.2：求取条纹图像的最佳尺度因子分布图ar(x，y)和相对相位分布图φ(x，y)，求取ar(x，y)和φ(x，y)在坐标(x，y)处数值的方法如下：
[0019] 求出W(a1，b)的对应的模矩阵A(a1，b)和
角度矩阵 搜索矩阵A(a1，b)第x列中值最大的元素，求出该元素在矩阵A(a1，b)中的行标号为amax，则arx＝10+0.2×amax，arx为条纹图像的最佳尺度因子分布图ar(x，y)在坐标(x，y)处数值，矩阵 中第x列中行标号为amax的元素数值为条纹图像的相对相位分布图φ(x，y)在坐标(x，y)处数值，[0020] 遍历条纹图像所有坐标点，求得条纹图像的最佳尺度因子分布图ar(x，y)和相对相位分布图φ(x，y)，
[0021] 步骤3：建立质量图Q(x，y)，
[0022] 步骤3.1：对频率为f0的一维正弦信号做小波变换：
[0023]
[0024] 求出二维复数矩阵W1(a1，b)的模矩阵A1(a1，b)，搜索每列元素中数值最大的元素，记录其行标号，求得这些行标号的平均值 求得最优尺度因子ar1：
[0025]
[0026] 步骤3.2：求质量图Q(x，y)：
[0027] Q(x，y)＝|ar(x，y)-ar1|
[0028] 步骤4：根据质量图Q(x，y)将相对相位分布图φ(x，y)分成两个部分，[0029] 步骤4.1：建立高度c，宽度r的0值矩阵D，
[0030] 步骤4.2：遍历质量图Q(x，y)的每一个点，求出质量图Q(x，y)出现频率最高的元素数值Q1，得到
阈值T＝1.05×Q1，
[0031] 步骤4.3：遍历质量图Q(x，y)的每一个点，当Q(x，y)中的数值大于阈值T时，矩阵D中相应
位置元素置1，
[0032] 步骤5：展开相对相位分布图φ(x，y)中D矩阵相应位置元素值为0的点的相位，[0033] 步骤5.1：建立高度c，宽度r的0值矩阵S，
[0034] 步骤5.2：在相对相位分布图φ(x，y)的中心像素的20×20领域内选择一个D矩阵相应位置元素值为0的点作为展开的起始点，起始点的绝对相位值记为该点在相对相位分布图φ(x，y)中的数值，沿着行方向向图像两侧展开D矩阵相应位置元素值为0的点的相位，每展开一个像素点，该点在S矩阵相应位置元素值置1，在直接展开过程中如遇到某些点的同一行的前一点并未处理，可以将该点的上一行或者下一行的相邻点视为前一点然后展开，每行都如此运算，
[0035] 每一点具体展开过程如下：
[0036]
[0037] 式中， 是前一点的绝对相位值， 是当前点的相对相位值， 是求得的当前点的绝对相位值，round是取整函数，
[0038] 步骤6：展开相对相位分布图φ(x，y)中S矩阵相应位置元素值为0的点的相位，[0039] 步骤6.1：在S矩阵中标记0值的连通域，对于每个连通域做如下处理：
[0040] 步骤6.2：选择连通域中的任一边界点做为起始点，入当前为空的堆栈，[0041] 步骤6.3：查找起始点四邻点中S矩阵相应位置元素值为0的点，如果四邻点中找不到这样的点，跳到步骤6.5；如果可以找到这样的点，展开这些点的相位，将这些点在S矩阵中相应位置元素值置1，并将这些点存入堆栈，
[0042] 步骤6.4：对堆栈中的点按照其在条纹图像质量图Q(x，y)中相应位置的数值大小进行排序，质量值最大点放置于堆栈的顶部，
[0043] 步骤6.5：将堆栈的顶部的点取出做为起始点，判断堆栈是否为空，如果为空，则该连通域展开过程结束；如果不为空，转到步骤6.3，
[0044] 步骤7：经过步骤5和步骤6处理，得到条纹图像的绝对相位图φ1(x，y)，根据经典光栅投影的相位到高度的转换公式，最终求得测量物体的三维信息。
[0045] 有益效果：与
现有技术相比，本发明具有如下优点：首先，本发明采用小波变换法求解条纹图像的相对相位，相比于傅里叶变换三维测量法，避免了频谱混叠现象，可以得到更加精确的相对相位结果，相比于相移法，本发明只需投影一幅黑白条纹图像，可以实现动态测量；其次，本发明建立质量图的方法充分利用了小波变换矩阵中的信息，相比其他方法，在相位展开中使用本发明方法所建立的质量图，能够得到一条可靠性更高的相位展开路径，降低了相位展开中出错的概率，可以得到更加精确的绝对相位结果；最后，本发明利用质量图将相对相位图分成两个部分，对这两个部分分别采用不同的算法进行相位展开，相比于传统的洪水相位展开算法，极大减少了相位展开所需的运算时间，同时保持了较高的鲁棒性。综上，本发明能够快速准确地获取被测物体的三维高度信息，具有较好的实时性和鲁棒性。
附图说明
[0047] 图2是步骤6中采用洪水算法展开“质量”较低的点的相位的具体过程流程图。
[0048] 图3是CCD采集到的以
泡沫塑料为被测物体的变形条纹图像。
[0049] 图4是采用小波变换法求解得到的相对相位分布图。
[0050] 图5是得到的最佳尺度因子分布图。
[0051] 图6是根据质量图计算的阈值对最佳尺度因子分布图做阈值分割得到的二值图像。
[0052] 图7是步骤5过后的展开相位图像。
[0053] 图8是最终得到的绝对相位分布图。
[0054] 图9是根据经典的光栅投影的相位到高度的转换公式，得到的表示物体三维信息的点
云图像。
具体实施方式
[0055] 下面结合附图对本发明的具体实施方式做进一步说明。在windows
操作系统下选用VC++6.0作为编程工具，对CCD采集到的变形条纹图像进行处理。该实例采用泡沫塑料作为被测物体，最终得到比较精确的含有护板三维信息的全场绝对相位分布。
[0056] 图1是本发明的整个过程的流程图。
[0057] 图2是步骤6采用洪水算法展开“质量”较低的点的相位的具体过程流程图。
[0058] 针对小波变换轮廓术中的精确性和快速性的问题，本发明首先采用小波变换求解采集到的变形条纹图像的相对相位分布，利用小波变换脊处的尺度因子建立质量图指导相对相位图的质量展开，在此
基础上，为了提高相位展开算法的速度，根据质量图将相对相位分布分成质量好坏的两个部分，质量普遍较好的部分采用扫描线法直接相位展开，质量普遍较差的部分采用质量指导选取最佳路径的方法。这种方法省去了全局质量指导算法的一些不必要的步骤，提高了相位展开算法的效率，而且保持了较高的测量精度。在得到精确的绝对相位分布后，根据经典光栅投影的相位到高度转换公式，最终得到测量物体的三维信息。
[0059] 本发明基于小波变换的三维测量方法，具体实施步骤如下：
[0060] 步骤1：将黑白条纹图像投影到被测物体表面，使用CCD对被测物体表面进行拍摄，得到一幅高度为c、宽度为r的变形条纹图像g(x，y)：
[0061] g(x，y)＝A(x，y)+B(x，y)cos[2πf0x+φ(x，y)]
[0062] 其中，A(x，y)是背景光强分布，B(x，y)是物体表面反射率，由于环境光是不变的，所以这A(x，y)、B(x，y)两个参数可以当做常数处理，f0是正弦条纹频率，可取f0＝32，φ(x，y)是待求的相对相位分布，(x，y)表示变形条纹图像的二维坐标，图3为变形条纹图像，
[0063] 步骤2：对变形条纹图像逐行做小波变换，得到变形条纹图像的相对相位分布图，具体过程如下：
[0064] 步骤2.1：将y视为常数，采用一维连续小波变换对变形条纹图像g(x，y)的第y行进行处理，处理过程为：
[0065]
[0066] 其中，a是尺度因子，取值范围为10到50，每隔0.2取一个值，b是平移因子，取值范围为1到条纹图像宽度r，每隔1取一个值，单位为像素，获得的W(a1，b)是一个200行r列的二维复数矩阵，a1是矩阵W(a1，b)中元素的行标号，称第y行小波变换矩阵，M*a，b(x)是Ma，b(x)的共轭函数，M(x)是小波函数，表达式如下：
[0067]
[0068] 其中，fb是小波函数的带宽，fc是小波函数的中心频率，i为单位复数，本发明采用fb＝2，fc＝0.845，进行小波变换，
[0069] 步骤2.2：求取条纹图像的最佳尺度因子分布图ar(x，y)和相对相位分布图φ(x，y)，求取ar(x，y)和φ(x，y)在坐标(x，y)处数值的方法如下：
[0070] 求出W(a1，b)的对应的模矩阵A(a1，b)和角度矩阵
[0071]
[0072]
[0073] 搜索矩阵A(a1，b)第x列中值最大的元素，称为条纹图像(x，y)坐标处的点的小波变换脊，求出该元素在矩阵A(a1，b)中的行标号为amax，则arx＝10+0.2×amax，这个式子实际上是根据行标号为amax寻找小波变换脊处的尺度因子的数值，arx为条纹图像的最佳尺度因子分布图ar(x，y)在坐标(x，y)处数值，矩阵 中第x列中行标号为amax的元素数值为条纹图像的相对相位分布图φ(x，y)在坐标(x，y)处数值，
[0074] 遍历条纹图像所有坐标点，求得条纹图像的最佳尺度因子分布图ar(x，y)和相对相位分布图φ(x，y)，实验效果分别见图5和图4，
[0075] 步骤3：建立质量图Q(x，y)，建立质量图的理论基础是小波变换矩阵的模值代表了条纹图像的局部信号和小波函数的相似程度，由于本发明中选用的小波函数是morlet函数，具有较好的正弦性，而根据步骤2.2，小波变换矩阵的模值在小波变换脊处达到最大，所以可以采用小波变换脊处的尺度因子判断条纹图像局部信号的正弦性好坏。当条纹图像某一点的最佳尺度因子过大时，说明该像素点附近的局部信号和低频率的正弦信号比较吻合，可能是无条纹区域，故该点的“质量”即可靠性可以判定为低；当该点的最佳尺度因子过小时，说明该像素点附近的局部信号与高频率的正弦信号比较吻合，可能是受到噪声干扰，该点的“质量”也可以判定为低。故首先应该找到这样一个数值，当像素点的最佳尺度因子等于该数值时，像素点的“质量”最好，步骤3.1为找到这样一个数值的过程，步骤3.2为建立质量图的过程，质量图中每一点质量的含义实际上是该点的尺度因子与最优尺度因子的“距离”。本步骤中求取的质量图用于在步骤5中将相对相位分布图分成两个部分，在步骤6中指导相位展开过程，得到相位展开的最佳路径，
[0076] 步骤3.1：对频率为f0的一维正弦信号做小波变换：
[0077]
[0078] 求出二维复数矩阵W1(a1，b)的模矩阵A1(a1，b)，搜索每列元素中数值最大的元素，记录其行标号，求得这些行标号的平均值 求得最优尺度因子ar1：
[0079]
[0080] 步骤3.2：求质量图Q(x，y)：
[0081] Q(x，y)＝|ar(x，y)-ar1|
[0082] 本发明所建立的质量图中，质量值越大，说明该像素点的可靠性越低，[0083] 步骤4：根据质量图Q(x，y)将相对相位分布图φ(x，y)分成两个部分，[0084] 步骤4.1：建立高度c，宽度r的0值矩阵D，D矩阵和条纹图像大小相同，用于区分根据质量图将相对相位图分成的两个部分，
[0085] 步骤4.2：遍历质量图Q(x，y)，统计得到取值最多的质量值Q1，得到阈值T＝1.05×Q1，
[0086] 步骤4.3：遍历质量图Q(x，y)，当Q(x，y)中的数值大于阈值T时，矩阵D中相应位置元素置1，图6为D阵的图像模式，其中白色区域代表质量较好优先展开的点，黑色区域为需要采用洪水算法展开的点的区域，
[0087] 步骤5：展开相对相位分布图φ(x，y)中D矩阵相应位置元素值为0的点的相位，该步骤展开的点都是“质量”较好的点，且质量值相差不多，按照传统的洪水相位展开算法，这些点仍然要参与排序，但实际上排序后和排序前的相位展开路径对于精度的影响不大，故这一部分采用最简单的扫描线相位展开算法即可得到比较精确的结果，而运算且速度非常快，
[0088] 步骤5.1：建立高度c，宽度r的0值矩阵S，S矩阵与条纹图像大小相同，用于标记相对相位分布图φ(x，y)的点是否已经展开，
[0089] 步骤5.2：在相对相位分布图φ(x，y)的中心像素的20×20领域内选择一个D矩阵相应位置元素值为0的点作为展开的起始点，起始点的绝对相位值记为该点在相对相位分布图φ(x，y)中的数值，沿着行方向向图像两侧展开D矩阵相应位置元素值为0的点的相位，每展开一个像素点，该点在S矩阵相应位置元素值置1，在直接展开过程中如遇到某些点的同一行的前一点并未处理，可以将该点的上一行或者下一行的相邻点视为前一点然后展开，如果这两点也没有处理，留待步骤6处理，每行都如此运算，
[0090] 每一点具体展开过程如下：
[0091]
[0092] 式中， 是前一点的绝对相位值， 是当前点的相对相位值， 是求得的当前点的绝对相位值，round是取整函数，图7为经过步骤5过程后的相位分布图，[0093] 步骤6：展开相对相位分布图φ(x，y)中S矩阵相应位置元素值为0的点的相位，这一部分为“质量”较差的点，为达到较高展开精度故采用洪水相位展开算法处理，洪水算法的流程图如图2所示，由于这部分点数量较少，所以运算时间也较快，[0094] 步骤6.1：在S矩阵中标记0值的连通域，二值矩阵的一个连通域是指二值矩阵中相互连通且取值相同的元素的集合，如图6中每个封闭的黑色区域都是一个连通域，对于每个连通域做如下处理：
[0095] 步骤6.2：选择连通域中的任一边界点做为起始点，入当前为空的堆栈，[0096] 步骤6.3：查找起始点四邻点中S矩阵相应位置元素值为0的点，如果四邻点中找不到这样的点，跳到步骤6.5；如果可以找到这样的点，展开这些点的相位，将这些点在S矩阵中相应位置元素值置1，并将这些点存入堆栈，
[0097] 步骤6.4：对堆栈中的点按照其在条纹图像质量图Q(x，y)中相应位置的数值大小进行排序，质量值最大点放置于堆栈的顶部，
[0098] 步骤6.5：将堆栈的顶部的点取出做为起始点，判断堆栈是否为空，如果为空，则该连通域展开过程结束；如果不为空，转到步骤6.3，
[0099] 步骤7：经过步骤5和步骤6处理，已经遍历了条纹图像的相对相位分布图φ(x，y)的所有点，得到条纹图像的绝对相位图φ1(x，y)，如图8所示为被测物体的绝对相位图，在个人计算机上，基于VC++平台，本发明展开相位运算时间为27s，而全场相位展开运算时间近90min，可见本发明方法实时性较好，根据经典光栅投影的相位到高度转换公式，最终求得测量物体的三维信息，图9为将物体的三维信息以点云的形式表现出来。
[0100] 由条纹图像的绝对相位图φ1(x，y)到物体高度的公式如下：
[0101]
[0102] 其中，l、d是测量系统的几何参数，l是投影仪到测量平面的距离，d是CCD摄像头到投影仪的距离， 表示相位变化量， 为展开相位结果，为初始相位结果，由测量参考面决定，ω0为投影光栅的角频率，由系统标定可得。
