基于剪切原理消除载频的方法.pdf

《基于剪切原理消除载频的方法.pdf》由会员分享，可在线阅读，更多相关《基于剪切原理消除载频的方法.pdf（14页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 102288203 A (43)申请公布日 2011.12.21 CN 102288203 A *CN102288203A* (21)申请号 201110193880.5 (22)申请日 2011.07.12 G01D 3/028(2006.01) G01B 11/25(2006.01) (71)申请人 昆明理工大学 地址 650093 云南省昆明市五华区学府路 253 号 (72)发明人 饶帆 钱晓凡 (54) 发明名称 基于剪切原理消除载频的方法 (57) 摘要 本发明涉及一种数字全息和光学投影栅的三 维轮廓术再现光场中消除载频的方法， 具体地说 是一种基于剪切原。

2、理消除载频方法， 属于光学技 术领域。将数字全息图输入到计算机中进行衍射 计算， 得到 KM 个数据点上的二维全息再现复光 场， 然后在计算机中实现光场剪切干涉， 通过剪切 算法计算二维全息再现复光场相位沿 x 轴方向和 y 轴方向的梯度， 减少欠采样的发生， 最后利用再 现复光场相位的梯度减去其平均值得到消除载频 后的二维全息再现复光场相位梯度， 通过最小二 乘建立 KM 矩形网格上的离散泊松方程， 再用离 散余弦变换求解计算出相位值， 得到相位值的最 小二乘解。 该方法有利于降低相位解包裹的难度， 能够有效地消除载频， 使重构的相位值更接近待 测的真实相位值。 (51)Int.Cl. (1。

3、9)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 1 页 说明书 4 页 附图 8 页 CN 102288216 A1/1 页 2 1. 一种基于剪切原理消除载频的方法， 其特征在于按以下步骤进行 ： （1） 将数字全息图输入到计算机中进行衍射计算， 得到KM个数据点上的二维全息再 现复光场 A ； （2） 在计算机中进行光场剪切， 对典型的KM个数据点上的二维全息再现复光场 A， 分 别沿x轴方向和y轴方向作s个像素的平移 （剪切） ， 创建新的剪切光场 Bx、 By； （3） 将二维全息再现复光场 A 分别与得到的剪切光场 Bx、 By相除， 得到新光场 Cx、 Cy， 。

4、计 算新光场 Cx、 Cy的相位， 得到二维全息再现复光场 A 相位沿x轴方向和y轴方向的梯度 Dx、 Dy； （4） 再分别将 Dx、 Dy减去其本身的平均值， 得到复光场 A 相位沿x轴方向和y轴方向的 新梯度 Ex、 Ey， 新梯度 Ex、 Ey已经不含载频， 从而实现对载频的消除 ； （5） 利用消除载频后的相位梯度 Ex、 Ey， 通过最小二乘法建立KM矩形网格上的二维离 散泊松方程， 再用离散余弦变换求解计算出相位值的最小二乘解， 即为消除载频后的再现 复光场相位。 权 利 要 求 书 CN 102288203 ACN 102288203 A CN 102288216 A1/4 页。

5、 3 基于剪切原理消除载频的方法 技术领域 0001 本发明涉及一种数字全息和光学投影栅的三维轮廓术再现光场中消除载频的方 法， 具体地说是一种基于剪切原理消除载频方法， 属于光学技术领域。 背景技术 0002 利用数字全息的再现光场可以得到物体的相位分布， 作为一种光学测量的手段， 因其具有非接触、 高灵敏度、 实时和全场等优点， 已经在断裂力学、 形貌测量、 生物粒子监 视、 高精度的温度和压强的测量、 构件缺陷检测、 力学和微电子学中的无损检测等众多领域 得到广泛的应用。为了克服数字全息再现过程中孪生像和零级像的影响， 通常需要使用离 轴全息， 因此再现光场的相位中包含有载频的相位， 直。

6、接影响到真实相位的测量， 必须消除 载频相位才能重构物体真实相位。 0003 光学投影栅的三维轮廓术是将条纹投影到三维物体的表面， 用被三维物体轮廓调 制的形变条纹获得包裹的相位信息， 再经过相位展开 （解包裹） 得到物体三维形貌的方法。 在该轮廓术中， 为了提高精度， 通常要使投影栅具有较高的空间频率 （载频） ， 除非在频域中 移频， 重构的相位中一般也含有载频的相位， 同样必须消除载频才能重构物体真实三维形 貌。 发明内容 0004 本发明的目的是克服现有技术中存在的问题， 提供一种能较好地消除载频给解包 裹带来的困难， 同时避免人为查找载频带来的误差， 继而获得比传统方法更好的重构相位。

7、 的方法。 0005 本发明所采用的技术方案的步骤是 ： （1） 将数字全息图输入到计算机中进行衍射计算， 得到KM个数据点上的二维全息再 现复光场 A ； （2） 在计算机中进行光场剪切， 对典型的KM个数据点上的二维全息再现复光场 A， 分 别沿x轴方向和y轴方向作s个像素的平移 （剪切） ， 创建新的剪切光场 Bx、 By； 数据点KM 和像素s的数量根据实际需要确定 ； （3） 分别将二维全息再现复光场 A 与得到的剪切光场 Bx、 By相除， 得到新光场 Cx、 Cy， 计 算新光场 Cx、 Cy的相位， 得到二维全息再现复光场 A 相位沿x轴方向和y轴方向的梯度 Dx、 Dy 新光。

8、场 C 的相位恰好是二维全息再现复光场 A 相位x轴方向和y轴方向的梯度 ； （4） 再分别将 Dx、 Dy减去其本身的平均值， 得到复光场 A 相位沿x轴方向和y轴方向的 新梯度 Ex、 Ey， 新梯度 Ex、 Ey已经不含载频， 从而实现对载频的消除 （即从 Dx、 Dy中减去其本 身的平均值， 在相位解包裹前完成消除载频） ； （5） 利用消除载频后的相位梯度 Ex、 Ey， 通过最小二乘法建立KM矩形网格上的二维 离散泊松 （Poisson） 方程 ： 其中 ,为 A 的待求相位， 而 说 明 书 CN 102288203 ACN 102288203 A CN 102288216 A2。

9、/4 页 4 可以通过 Ex、 Ey计算得到 ， 再用离散余弦变换 （DCT） 求解计算出相位值的最小二乘 解 （计算二维全息再现复光 场相位） ， 即为消除载频后的再现复光场相位。 0006 二维离散泊松 （Poisson） 方程是数学中形如的微分方 程， 离散余弦变换 （DCT） 是与傅里叶变换相关的一种变换， 它类似于离散傅里叶变换， 但是 只使用实数。离散余弦变换相当于一个长度大概是它两倍的离散傅里叶变换， 这个离散傅 里叶变换是对一个实偶函数进行的。离散余弦变换经常被信号处理和图像处理使用， 用于 对信号和图像进行有损数据压缩。用离散余弦变换求解离散泊松方程的原理是 ： 将方程两 边。

10、作离散余弦变换， 通过函数导数 （梯度）的余弦变换与函数u(x,y) 的余 弦变换之间的关系可以建立新方程， 求解该方程可以得到待求函数u(x,y)。 0007 本发明将剪切原理引入到数字全息和光学投影栅三维轮廓术中， 消除了载频对真 实相位的影响。通过将数字全息再现光场， 在计算机中作人为的平移， 实现沿x轴方向和y 轴方向的剪切， 再从相位梯度中减去其平均值， 实现了对载频的消除， 即在相位解包裹前完 成消除载频。 0008 本发明针对 CCD、 CMOS 的分辨率只有 100 200 线 / 毫米， 且有效光敏面尺寸小， 限制了光电转换器件的空间带宽积， 导致数字全息再现光场相位的空间频。

11、率一般较高， 常 常出现包裹相位欠采样， 使包裹相位中存在大量不可靠数据点， 给消除载频带来了不小困 难， 甚至载频无法彻底的消除等问题， 可以很好地减小欠采样对消除载频的影响， 继而获得 比传统方法更好的相位解包裹结果。 附图说明 0009 图 1 是本发明再现光场剪切原理示意图 ； 图 2 是本发明实施例 1 采集数字全息图光路图， 图中 ： 1- 数字全息重构原始再现光场， 2- 剪切 （平移） 后的光场， 3- 剪切量 s， 4- 坐标轴， 5- 剪切方向， 6- 激光器， 7- 分束镜， 8- 全 反镜， 9- 显微物镜， 10- 针孔滤波器， 11- 准直透镜， 12- 待测物， 。

12、13- 分束镜， 14- 汇聚透镜， 15-CMOS， 16- 计算机， 17- 全反镜， 18- 显微物镜， 19- 针孔滤波器， 20- 准直透镜 ； 图 3 是本发明实施例 1 蜡烛燃烧前消载频后的干涉条纹 ； 图 4 是本发明实施例 1 蜡烛燃烧后消载频后的干涉条纹 ； 图 5 是本发明实施例 1 实验所得干涉图 ； 图 6 是本发明实施例 1 计算所得的解包裹相位图 ； 图 7 是本发明实施例 1 计算所得的干涉图。 0010 图 8 是本发明本发明实施例 2 光学投影栅三维轮廓术的记录光路示意图， 图中 ： 21- 计算机， 22- 投影机， 23- 数字相机， 24- 待测三维物。

13、体 （比如人的头部） ； 图 9 是本发明本发明实施例 2 记录的投影栅变形检测图 ； 图 10 是本发明本发明实施例 2 用传统的四步相移法得到的包裹相位 （将其乘以虚数 j=， 并取 e 指数， 可得到复光场 A） ； 图 11 是本发明本发明实施例 2 计算得到的不含载频新梯度 Ex； 说 明 书 CN 102288203 ACN 102288203 A CN 102288216 A3/4 页 5 图 12 是本发明本发明实施例 2 计算得到的不含载频新梯度 Ey； 图 13 是本发明本发明实施例 2 计算得到的三维物体相位 （三维形貌。已经去除载频） ； 图 14 是同一个三维物体没有。

14、去除载频时得到的相位 ； 图 15 是本发明本发明实施例 2 最终输出的三维物体形貌。 具体实施方式 0011 下面结合附图和实施例， 对本发明作进一步说明， 但本发明的内容不限于所述范 围。 0012 实施例 1 ： 消除数字全息再现光场中的载频。 0013 首先， 用所示的数字全息图的记录光路记录一幅数字全息图 （参见图 2） 。激光器 6 发出的激光束通过分束镜 7 分为两束， 其中一束经显微物镜 9 汇聚并通过针孔滤波器 10 后， 经准直透镜 11 变为平行光， 照射到待测物体 12 上， 透射光通过分束镜 13 作为物光经汇 聚透镜 14 汇聚到 CMOS 15 上 ； 另一束经全。

15、反镜 8、 17 和显微物镜 18 汇聚并通过针孔滤波器 19后， 照射到分束镜13上， 作为参考光经汇聚透镜14到达CMOS 15上。 参、 物光在CMOS 15 上干涉， 光强转换为电信号传输到计算机 16 中保存为数字全息图。 0014 然后， 按以下步骤消除数字全息光场中的载频 ： （1） 利用上述数字全息实验装置记录一幅激光通过蜡烛后的数字全息图， 并将数字全 息图输入到计算机中进行衍射计算， 得到 450350 个数据点上的二维全息再现复光场 A。 0015 （2） 在计算机中进行光场剪切， 分别将该二维全息再现复光场 A 沿x轴方向和y轴 方向作 s=1 个像素的平移剪切， 创建。

16、新的剪切光场 Bx、 By； （3） 分别将二维全息再现复光场 A 与得到的剪切光场 Bx、 By相除， 得到新光场 Cx、 Cy， 计 算新光场 Cx、 Cy的相位， 得到二维全息再现复光场 A 相位沿 x轴方向和y轴方向的梯度 Dx、 Dy； （4） 再分别将相位 Dx、 Dy减去其本身的平均值， 得到复光场 A 相位沿x轴方向和y轴方 向的新梯度 Ex、 Ey， 实现对载频的消除 ； （5） 利用消除载频后的相位梯度 Ex、 Ey， 通过最小二乘法建立KM（450350） 矩形网 格上的离散泊松方程， 再用离散余弦变换求解计算出相位值的最小二乘解， 即为再现复光 场相位。 0016 在本。

17、例中， 比较计算所得的干涉图 （图 7） 与实验所得干涉图 （图 5） ， 两者有很好的 吻合度。 0017 实施例 2 ： 消除光学投影栅三维轮廓术再现光场中的载频。 0018 参见图8， 首先， 将计算机21中预先储存的四幅相位相差/2的投影栅图， 用投影 仪 22 投射到待测三维物体 24 上， 投影栅变形检测图通过数字相机 23 采集后储存到计算机 21 中， 利用传统的四步相移算法得到物体的包裹相位图。 0019 接着， 按以下步骤进行 ： （1） 将包裹相位乘以虚数 j=， 并取 e 指数， 得到的 12801024 个数据点上的二维 复光场 A。 0020 （2） 在计算机中进行。

18、光场剪切， 分别将该二维复光场 A 沿x轴方向和y轴方向作 s=1 个像素的平移剪切， 创建新的剪切光场 Bx、 By； 说 明 书 CN 102288203 ACN 102288203 A CN 102288216 A4/4 页 6 （3） 分别将二维复光场 A 与得到的剪切光场 Bx、 By相除， 得到新光场 Cx、 Cy， 计算新光场 Cx、 Cy的相位， 得到二维复光场 A 相位沿x轴方向和y轴方向的梯度 Dx、 Dy； （4） 再分别将相位 Dx、 Dy减去其本身的平均值， 得到复光场 A 相位沿x轴方向和y轴方 向的新梯度 Ex、 Ey， 实现对载频的消除 ； （5） 利用消除载频。

19、后的相位梯度 Ex、 Ey， 通过最小二乘法建立KM（12801024） 矩形 网格上的离散泊松方程， 再用离散余弦变换求解计算出相位值的最小二乘解， 即为复光场 相位。 0021 在本例中， 比较图 13 和图 14， 可以看到消除载频后才能得到物体正确的三维形貌 （相位） 。 说 明 书 CN 102288203 ACN 102288203 A CN 102288216 A1/8 页 7 图 1 图 2 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A2/8 页 8 图 3 图 4 说 明 书 附 图 CN 102288203 AC。

20、N 102288203 A CN 102288216 A3/8 页 9 图 5 图 6 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A4/8 页 10 图 7 图 8 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A5/8 页 11 图 9 图 10 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A6/8 页 12 图 11 图 12 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A7/8 页 13 图 13 图 14 说 明 书 附 图 CN 102288203 ACN 102288203 A CN 102288216 A8/8 页 14 图 15 说 明 书 附 图 CN 102288203 ACN 102288203 A 。