外差式三频不等步相移解相位方法
【专利摘要】本发明涉及一种外差式三频不等步相移解相位方法，本方法选取三个频率，首先采用四步相移法求得中间频率的包裹相位及平均光强，然后在光强已知的情况下采用二步相移法求得另外两个频率对应的包裹相位，最后根据外差方法解得真实相位，另外提出将等值线正余弦滤波方法应用于光栅投影三维测量中，利用条纹方向信息对包裹相位在等相位线上进行滤波，本算法利用较少投影图片进行准确测量，同时充分利用条纹方向信息对包裹相位图进行滤波，有效去除包裹相位图中的不一致点，外差式三频不等步相移解相位方法用在相位法进行三维形状测量中可有效提高测量效率和测量准确度。
【专利说明】外差式三频不等步相移解相位方法
【技术领域】
[0001] 本发明属于图像处理【技术领域】，特别涉及一种外差式三频不等步相移解相位方 法，能够高效精确获取物体的三维信息。
【背景技术】
[0002] 光栅投影三维测量技术是一种利用相位信息进行三维轮廓精确测量的方法，目前 获取相位的方法主要有相移法与傅里叶变换轮廓法。相移法因精度高、分辨率高且对物体 表面反射率变化不太敏感而被广泛应用。该方法通过采集多帧相移条纹图来计算包含被测 物表面形变信息的包裹相位值。由于包裹相位值在[-^n]区间内，因此为了获得连续 相位映射，需要将包裹相位展开，得到真实的相位值。相位解包裹方法可分为基于单频率和 基于多频率的解包裹方法，基于单频率的解包方法测量时只获取一组包裹相位值，故其在 解包范围、解包精度及误差容忍度方面均存在不足，而基于多频率的解包方法测量时获取 多组包裹相位值，增大了解包范围、提高了解包精度、降低了误差影响。基于多频率的解包 裹方法主要有格雷（Gray)码法和多频相移方法。格雷码法对光照和被测物表面明暗比较 敏感，经常出现像素边界解码不正确，容易受到噪声干扰且测量范围受格雷码条纹数限制。 多频相移方法通过增加不同频率条纹来提高测量精度及测量范围，增强适应性。目前多频 相移方法已经成为工业上应用最广的目标物高精度三维形貌重建方法，但多频相移方法需 要拍摄不同频率的多幅图像。
[0003] 基于光栅投影的三维测量技术，在读取采用余弦编码的条纹图样时，其图像灰度 分布应为余弦分布。但由于投影系统的非线性输出、摄像设备中的电子噪声干扰和环境噪 声等因素的影响，导致在编码图样的灰度分布中存在大量噪声，这将导致在相位恢复时产 生误差，影响不连续物体表面的三维重建。在使用优化算法减小相位误差的过程中，首先需 要对投影系统进行非线性校正，但校正过程需要大量时间，且要求投影仪的伽马值相对固 定。其次，需要对相位进行误差补偿，基于样条拟合的方法在相位计算时需要经过多次迭 代，运算量较大；基于查找表的方法要将预先计算出的相位误差值存储在一个查找表里用 于补偿误差，计算相位误差的过程需要大量时间且对该方法对环境要求比较高；基于图像 滤波的方法对相位图像进行普通频域或空域滤波，可以极大地消除噪声，降低相位误差，但 对含有高度跳变的曲面，会损失物体边界跳变信息。
【发明内容】
[0004] 本发明的目的就是为克服现有技术的不足，提供一种外差式三频不等步相移解相 位方法，利用方法中的关键技术实现在减少拍摄图像的数量获得准确的测量精度同时能够 降低噪声和物体高度跳变对精度的影响。
[0005] 本发明是通过这样的技术方案实现的：一种外差式三频不等步相移解相位方法， 利用数字投影仪、摄像机、PC机组成的测量平台实现条纹的投影照明和条纹图的数据采集 和处理，其特征在于，包括下列步骤：
[0006] (1)投影三个不同频率共8幅条纹图，并采集条纹图像，其中采用中间频率f2投射 四步相移条纹，频率和f3投射二步相移条纹；
[0007](2)采用不等步相移解相原理求解频率f\、f2、f3对应的包裹相位 4 4 Z内(^7)。
[0008] (3)分别对求得的包裹相位值进行等值线正余弦滤波；
[0009] (4)采用外差原理求出频率f2变形条纹图像的真实相位;
[0010] (5)利用双目相位立体匹配方法计算物体三维表面点云数据。
[0011] 本发明的实质性特点是，选取三个频率余弦波，对每个频率的余弦波进行不同步 数的空间相移，只需投射并拍摄八幅条纹图即可获得较高质量的真实相位图。同时在对投 影系统进行非线性校正的基础上，针对外差法差频叠加包裹相位噪声的问题，提出在对包 裹相位图进行解包裹前采用等值线正余弦滤波方法对其进行降噪处理，本发明与现有技术 相比较具有如下优点：
[0012] 1.利用较少投影图片进行准确测量。
[0013] 2.充分利用条纹方向信息对包裹相位图进行滤波，有效去除包裹相位图中的不一 致点。
【专利附图】
【附图说明】
[0014] 图1 :本发明的测量方法流程图；
[0015] 图2 :外差原理示意图；
[0016] 图3 :条纹的等值线和法向曲线示意图。
【具体实施方式】
[0017]本发明的流程图如图1所示，首先投影并采集三个不同频率共8幅条纹图（中间 频率f2投射四步相移条纹，频率和f3投射二步相移条纹），然后采用不等步相移解相 原理计算背景光强及频率f2的包裹相位值，并在背景光强已知的情况下求得频率 和f3对应的包裹相位值jWfKv」和4抑dvj，之后对三个包裹相位分别进行等值线正余 弦滤波。其次采用外差原理将4差频叠加得到等效频率为f12的相位函 数9 12 (X，y)，将4 和4 差频叠加得到等效频率为f23的相位函数0 23 (x，y)， 并分别对其进行正余弦滤波。最后再通过二次外差解得含物体表面形变信息的相位函数 9 123(x，y)，利用双目相位匹配方法得到物体三维表面点云数据。
[0018] 下面结合附图，对本发明技术方案的具体实施过程加以说明。
[0019] 1.三个不同频率共8幅条纹图投影和采集
[0020] 计算机模拟生成8幅波长XpX2和X3分别为15、17、20(单位：像素/周期，投 影条纹宽度W= 1020像素）的满足三频不等步相移编码方案的灰度条纹图。使用投影仪 将8幅条纹图投射到被测物表面，用摄像机进行采集。
[0021] 2.不等步相移法解相位
[0022] 在光栅投影三维测量中，投射余弦编码条纹到被测物体表面，编码条纹被物体表 面信息调制形成变形条纹，假设投影光强满足标准余弦分布，则变形条纹图像的光强分布 函数为：
[0023] I(x,y)=a(x,v) + b(x,y)cos[270c + h(p(x,y)\ (1)
[0024]I(x，y)为摄像机拍摄的变形条纹图像光强，a(x，y)为背景光的强度，b(x，y)为调 制振幅，f为条纹频率，为待求解的相位函数。式中含有三个未知量，因此要解得同 一个待测点的相位值至少要获取三个光强函数值。
[0025] 三频率外差不等步相移解相位的方法是对三种不同频率条纹（三个频率分别记 为&、&和&)采用不同步数相移算法求解得到各自频率的包裹相位值。频率&采用标准 四步相移条纹，将余弦光栅移动4次，每次移动步长a=Ji/2,这样可以采集到的4幅图像 来求得物体的相位值，在这里设In(x，y)代表第n幅图像光强，则：
[0026] /"(x,v) =a(x,y)+b(x,v)cos[2^/；x+A(p2(x,y)+(n-1)^/2] (n=l, 2, 3, 4) (2)
[0027] 利用这四幅变形条纹图像求得相位为：
[0028]
【权利要求】
1. 一种外差式三频不等步相移解相位方法，利用数字投影仪、摄像机、PC机组成的测 量平台实现条纹的投影照明和条纹图的数据采集和处理，其特征在于，包括下列步骤： (1) 投影三个不同频率共8幅条纹图，并采集条纹图像，其中采用中间频率4投射四步 相移条纹，频率和f3投射二步相移条纹； (2) 采用不等步相移解相原理求解频率f\、f2、f3对应的包裹相位4内仏))、a ^(p<(x,y)\ (3) 分别对求得的包裹相位值进行等值线正余弦滤波； (4) 采用外差原理求出频率f2变形条纹图像的真实相位O2 ; (5) 利用双目相位立体匹配方法计算物体三维表面点云数据。
2. 根据权利要求1所述的外差式三频不等步相移解相位方法，其特征在于，步骤（2) 中， 在光栅投影三维测量中，投射余弦编码条纹到被测物体表面，编码条纹被物体表面信 息调制形成变形条纹，假设投影光强满足标准余弦分布，则变形条纹图像的光强分布函数 为： I(x,y) = a(x,y) + h{x, y)cos[27ifa + h(p{x,y)\ q) I(x，y)为摄像机拍摄的变形条纹图像光强，a(x，y)为背景光的强度，b(x，y)为调制振 幅，f为条纹频率，」批V,.V〉为待求解的相位函数；式中含有三个未知量，因此要解得同一个 待测点的相位值至少要获取三个光强函数值； 三频率外差不等步相移解相位的方法是对三种不同频率条纹采用不同步数相移算法 求解得到各自频率的包裹相位值；三个频率分别记为f\、f2和f3 ; 频率4采用标准四步相移条纹，将余弦光栅移动4次，每次移动步长a = /2,这样 可以采集到的4幅图像来求得物体的相位值，在这里设In(x，y)代表第n幅图像光强，则： I11 (x, v) = ?(x,.v) + b{x,y)cos[2rf2x + Acp2(x,y) + (?-\)nfl\ (n=l, 2, 3, 4) (2) 利用这四幅变形条纹图像求得相位为：
假设相机和投影仪有相当大的景深，同时被测物表面反射是线性的，则在4幅图像中 同一位置像素点（i，j)处的平均光强度系数a (x，y)是个常数； a(x, y) = (I1 (x, y)+I2(x, y)+I3(x, y)+I4(x, y))/4 (4) 频率和f3都采用二步相移方法，相移步长a = 3 Ji /2,这样可以各采集到正弦和余 弦两幅条纹图，光强函数分别为 /, (x, j') = a(x,r) + ^(^,_v)sin(2^/；x + ⑶ /6(x,_v) = ^(x, v) + Z)(x,]〇cos(2^/；x + A^(x,j')) (g) /7(x,_v) = i3(x, v) + v)sin(2^3x + A^3(x,_v)) ⑵ /8(x,_v) = a(x, v) + /?(x, v)cos(2^x + A^3(x,j')) (8) 利用式（4)求得a (x，y)，代入式（9)和式（10)可求频率&和&变形条纹图像的包裹 相位值4和d WfxJ)
其中In(x, y)代表第n幅图像光强，a(x, y)是平均光强度系数，arctan为求解反正切 函数。
3. 根据权利要求1所述的外差式三频不等步相移解相位方法，其特征在于，步骤（3) 中， 直接使用传统滤波方法对实际包裹相位图进行滤波，去除噪声的同时会损失相位周期 跳变信息；由于正、余弦图像的灰度值是连续变化的，对正、余弦图像滤波不会损害原图像 的周期跳变信息部分，本算法提出采用正余弦滤波方法对光栅投影三维测量过程中的包裹 相位图进行滤波，先计算相位图的正、余弦结果图，然后再将这两幅图像根据条纹方向信息 在等值线窗口内进行中值滤波，去除噪声，最后将经过滤波后的正、余弦图像做除法运算， 得到相位图； 等值线是与条纹当前走向局部平行的线，法向曲线是与条纹当前走向局部垂直的线， 如附图3所示；为确定条纹等值线和法向曲线先求取条纹方向； 设图像上某一点的坐标为PiO^yi),此点的条纹方向为Qi，由条纹方向图的定义 可知，方向图的方向值9是条纹等值线的切线方向；在当前点的邻域内可以近似用切 线代替条纹等值线，因此，由该点根据条纹方向图可以得到等值线上该点的两个相邻点 f*i+l (Xi+1，5^+1)和 Pi-I (X i-1， Yi-i)：
同理，继续向两边扩展，可分别得到Pi+2(xi+2，yi+2)和Pp 2O^yiJ等；由此可以得到 一条过当前点PiO^yi)且沿条纹方向的曲线，即为条纹等值线，曲线上各点的相位值与点 PiUi, Yi)的相位值相等，即在条纹的等值线上条纹的相位保持不变； 对于条纹图来说，沿着条纹切线方向的图像灰度差异最小，而垂直于条纹切线方向的 图像灰度差异最大；所以可以计算每一点沿着各个方向的灰度差异，取得最小值时对应的 方向即为条纹方向；梯度法是一种常用的求取条纹方向的方法； 由于数字图像只在整像素上取值，而由上述方法得到的点（Xi，yi)并不一定是整数点； 为减少条纹方向图只取整数像素点上的值所引起的误差，可以取与所得非整数点（Xi，yi) 相邻的四个整像素点的条纹方向值的线性插值作为该点的条纹方向值，以减小离散误差。
4. 根据权利要求1所述的外差式三频不等步相移解相位方法，其特征在于，步骤（4) 中， 外差原理是将两种不同频率f\、f2条纹对应的相位函数和经过差频叠加 得到一种新的等效频率为f12的相位函数912 (X)，不同频率f\、f2和f12条纹对应的波长分 别为X i、X 2和入12 ; 等效波长A12经过计算可表示为：
利用外差原理将包裹在区间[-^ n]内的相位图像展开，得到全场范围内连续的相 位图像；假设不同频率对应的波长分别为\2和A3值，则由A 1和\ 2计算得到等效波 长入12, A2和八计算得到等效波长入23，入12和入23计算得到等效波长入 123;为了实现全 场范围内连续的真实相位图像，等效波长A123应该大于等于投影仪投射的行像素数W(投 影仪投射坚条纹），即入123彡W; 不同频率条纹包裹相位展开的计算公式为： = InNi + (/ = L2J) (13) 式中，Oi为真实相位值，Ni为条纹整数级数，2奶为A i对应的相位函数； 对于待测物体形面上同一点不同条纹波长A a和入b (a分别等于1，2, 12 ;对应的b分 别等于2, 3, 23)，差频叠加得到等效波长A ab的相位函数0 ab(x)，测得的条纹相位值满足：
^2=(^ 2-12+ ^ 2-23) /2 (18) 式中Round为取整函数，9 123(x, y)为等效波长：V 123对应的相位，^V12为利用等效波 长入12及它的相位Q12 (x，y)得到的频率f2变形条纹图像的真实相位，O2^23为利用等效波 长入23及它的相位Q23 (x，y)得到的频率f2变形条纹图像的真实相位；为减小误差，真实相 位CD2为①2_12和①2_23的平均值。
5.根据权利要求1所述的外差式三频不等步相移解相位方法，其特征在于，步骤（5) 中，基于相位的立体匹配是根据空间同一点在左右摄像机像平面上，像点的局部相位等条 件和极线约束条件进行立体视觉匹配，找出物体上任意一点对应于左右摄像机图像的像点 坐标对；然后根据空间点对应于两台摄像机图像的像点点对和摄像机系统的标定结果，利 用空间点、像点和摄像机光心在一条直线上的原理，通过空间直线求交点的方法，求解空间 点的三维坐标，最终可得到三维点云。
【文档编号】G06T7/00GK104330052SQ201410675447
【公开日】2015年2月4日 申请日期:2014年11月21日 优先权日:2014年11月21日
【发明者】耿磊, 刘文超, 肖志涛, 张芳, 吴骏, 李月龙, 杨振杰, 叶琨, 甘鹏 申请人:天津工业大学