三维物体主菲涅尔波带全息成像法

杨涛+周皓+季冬+顾济华

摘 要： 为了提高三维物体计算全息图的生成速度，在主菲涅尔波带法的基础上提出了一种改进算法。首先依据衍射距离将二维主菲涅尔波带转换为两组一维波带；然后依据实际应用情况与显示装置进行尺寸匹配；最后将得到的波带叠加获得三维物体全息图。该算法与改进前相比，内存空间的占用较小，运算速度得到显著提高，生成的全息图经由DMD再现获得高质量的三维物体再现图像。

关键词： 计算全息； 主菲涅尔波带； 光电再现； DMD

中图分类号： TN911?34； O438.1 文献标识码： A 文章编号： 1004?373X（2014）20?0112?04

“Host” Fresnel wave zone holographic imaging method based on

three?dimensional objects

YANG Tao， ZHOU Hao， JI Dong， GU Ji?hua

（Department of Physics Science and Technology， Soochow University， Suzhou 215006， China）

Abstract： In order to improve the generating speed to the three?dimensional objects calculation hologram， an improved algorithm based on the “host” Fresnel wave zone method is proposed. The two?dimensional “host” Fresnel wave zone is converted to two one?dimensional wave zones according to the diffraction distance， and then the desired Fresnel wave zones are matched with the size of the display device according to the actual application. Finally， the obtained Fresnel wave zones are superimposed to achieve three?dimensional objects hologram. Compared to the previous algorithm， the improved algorithm occupies less memory space， its computing speed has been significantly improved and a high quality reconstructed image generated by the DMD can be obtained.

Keywords： CGH； “host” Fresnel wave zone； photoelectric reconstruction； DMD

0 引 言

人类对现实世界的感知中，80%的信息来源于双眼的视觉，而目前占据主流地位的显示技术仍是二维显示技术，虽然二维显示技术能满足大部分日常生活需求，但与三维显示相比，缺乏临场感和视觉冲击感。随着科学技术的飞速发展，三维立体显示技术越来越成熟，逐渐成为研究的热点[1?4]。目前三维立体显示技术主要分两种，一种是基于双目视差的立体显示技术，如眼镜/头盔式三维立体显示和光栅式自由立体显示[5?8]、合成成像立体显示[9?13]、体三维显示[14?15]等各种三维显示技术，这些技术都是通过提供具有一定视差的图像，使观察者形成一定的深度印象，从而达到三维立体显示的目的。另外一种则是显示真实立体物体影像的技术，如全息技术，它通过干涉完整地记录了物体的强度、相位等信息，再现像是真正意义的空间立体图像，再现图像可以直接进行裸眼观察，被大多数人认为是未来三维显示最理想的方案[16?19]。

目前采用全息技术显示三维物体的方法主要有光学全息、数字全息和计算全息。其中计算全息具有噪声低、可重复性好、可实现虚拟物体的实体化等特点。但由于三维物体复杂的数学描述和有限的计算机运算能力，计算效率成为了限制计算全息发展的瓶颈。研究者在不提高计算机硬件的前提下，研究了很多简化三维物体的计算全息算法，以达到提高计算速度的目的：Waters提出通过分别计算单个点光源在全息平面上的复振幅以叠加编码成计算全息图[20]；Leseberg将三维物体各个平面的菲涅尔全息图累和组成具有全息图[21]；张晓洁等人提出了主菲涅尔波带法[22?24]，在点物散法的基础上，通过规整同一深度下的所有采样物点来提高运算速度；李勇等研究者对全息图进行变量分离[25?26]，以达到减少计算量，提高运算效率的目的。

本文改进了主菲涅尔波带法以进一步提高全息图生成效率：将分离变量原理[17?19]与主菲涅尔波带法相结合，从主菲涅尔波带的计算公式中分离出关于两个轴上的独立分量，以独立分量代替主菲涅尔波带，并获取其他采样物点轴上的分量，最后通过计算得到各个菲涅尔波带。与传统主菲涅尔波带法相比，改进后的算法不仅补全了文献[23]中丢失的高频信息，而且计算量进一步减少，运算速度得到较大提高。将生成的计算全息图采用DMD进行光电再现得到了高质量的三维物体再现图像。

1 原 理

1.1 主菲涅尔波带法

若采用平面波垂直照射，对于给定衍射距离上的所有采样物点，它们完整的菲涅尔衍射波带都具有相同的衍射图样，只是菲涅尔波带的中心位置不同，根据这一特性，主菲涅尔波带法只需计算其中一个点的菲涅尔衍射波带，其他采样点的菲涅尔波带可以通过平移操作获得，这个需要计算的点称之为主点，一般取物平面中央位置的采样点，主点的全息图则称之为主菲涅尔波带，所有采样点的菲涅尔波带叠加得到全息图。

在满足奈奎斯特采样定理的前提下，主菲涅尔波带法物平面像素间隔[Δx，Δy]和接收面像素间隔[Δξ，Δη]关系另外还满足：

[Δx≥λzMHΔξ，Δy≥λzNHΔη] （1）

式中，[λ]为光源波长；[MH×NH]为接收面像素数目；[z]为垂直衍射距离。依据人眼的分辨率，[Δx]和[Δy]同样也有个上限值，不同的衍射距离有不同的上限值。本文从计算方便出发，物平面和接收面的像素间隔都选取下限值[27]。

从理论角度出发，主菲涅尔波带大小可以任意设置，假设主菲涅尔波带大小与接收面一致，如图1 所示，全息图的接收面限制在：

[ξ∈[-ξh，ξh]] （2）

[η∈[-ηh，ηh]] （3）

点[O]和[A]分别位于物平面的中心位置和右上角边缘位置，分别与衍射面上主菲涅尔波带中心[O′]和[A′]相对应。平移主菲涅尔波带求[A]点的菲涅尔波带时，主菲涅尔波带中心[O′]平移的总长度[d]，其中[x]轴和[y]轴平移长度分别为：

[dξ=ΔxMO2] （4）

[dη=ΔyNO2] （5）

由图1可见，平移之后，以[A′]为中心的主菲涅尔波带只有很少一部分信息被接受面接收，而大部分高频信息都丢失了。为了获得高质量的再现像，主菲涅尔波带必须设置的足够大，才能保证全息图的信息量，即平移后以[A′]为中心的主菲涅尔波带恰好能完整覆盖接收面，设此时主菲涅尔波带在[ξ]轴和[η]轴的长度分别为[xf]和[yf]，由几何关系可得：

[xf=2（dξ+ξh）] （6）

[yf=2（dη+ηh）] （7）

根据两平面间采样间隔关系，将距离换算成像素数目，主菲涅尔波带的像素数目最小应满足：

[M×N=λzMHΔξ2MO+MH×λzNHΔη2NO+NH] （8）

式中，[M]和[N]的值通常不是整数，常取比其值大的最小整数。

图1 主菲涅尔波带法平移原理图

1.2 改进算法

在主菲涅尔波带法的基础上，利用分离变量原理做了相应改进，基本思路是利用主菲涅尔波带在[ξ]轴和[η]轴上的独立分量代替主菲涅尔波带，以达到减少内存空间占用和计算量，提高计算效率的目的。根据菲涅尔的基本衍射公式，主菲涅尔波带在[ξ]轴和[η]轴上的独立分量[x]和[y]的离散形式分别满足以下公式：

[x=expjk[（m-M2-1）Δξ]22z] （9）

[y=expjk[（n-N2-1）Δη]22z] （10）

[C=exp（jkz）jλz] （11）

式中：[m∈[1，M]]；[n∈[1，N]]的正整数。

改进算法的基本思路如图2所示。按上述公式计算得到主菲涅尔波带的独立分量[x]和[y]，将平移操作的对象从主菲涅尔波带转换为[x]和[y]，[ξ]轴，[η]轴上的接收窗口分别限制在[ξ∈[-ξh，ξh]]和[η∈[-ηh，ηh]]以匹配显示器件的分辨率，按坐标相对位置平移得到其他采样点在[ξ]轴独立分量[MH×1]的独立分量[xn]和[η]轴[1×NH]的独立分量[yn]，通过矩阵乘法即可得到每个采样点的菲涅尔波带，叠加组合成全息图。

计算模拟运算中，以双精度类型将[M×N]的主菲涅尔波带保存到内存中，需要内存连续空间[size=8MN/1 0243]GB，改进的算法以独立分量x和y代替主菲涅尔波带，此过程仅需要消耗内存连续空间[size′=8（M+N）1 0243] GB，与未改进相比，减少了内存占用空间[Δsize=8（MN-M-N）1 0243] GB，共少计算了[（MN-M-N）]个元素的值。当实验条件确定时，改进前后内存占用空间与衍射距离的关系如图3所示。由图3可知，衍射距离越远，改进后的算法在内存占用和计算量上越有优势。

图3 改进前后内存占用空间对比

2 实 验

2.1 计算全息图

国际象棋中的“国王”、“皇后”和全息图的位置关系如图4（a）所示，“皇后”和“国王”分别位于全息图[0.3 m]处和[0.5 m]处，共有12 994个采样物点。实验中采用[532 nm]的半导体绿色激光器作为光源，计算像素数目为[1 024×768]的全息图，以匹配DMD的分辨率。设定原始物图大小为[47 mm（H）×47 mm（V）]，像素数目为[512×512]，接受面像素数数目为[1 024×768]，像素间隔为[3.45 μm（H）×3.45 μm（V）]。

根据公式（8）计算[0.3 m]处的主菲涅尔波带，满足要求的主菲涅尔波带的像素数为[28 672×21 504]，以双精度类型存储，[M×N]的主菲涅尔波带需要[4.59 GB]的内存连续空间，受限于实验室电脑（3 GB内存）内存空间的限制，无法将如此庞大的单个矩阵数组存放到内存中，故未改进的主菲涅尔波带法的程序无法顺利模拟，而改进后的算法，[28 672×1]的独立分量[x]和[1×21 504]独立分量[y]仅占用[0.383 MB]的内存连续空间，完全满足运行环境。改进后算法用时150 s得到如图4（b）所示的全息图，并进行了DMD光电再现实验。

图4 原始物图和计算全息图

2.2 DMD光电再现

DMD再现实验采用如图5所示的光路，对图4（b）的全息图进行了再现，图中采用532 nm的半导体绿色激光器，实验采用的DMD分辨率为1 024（H） ×768 （V），单个微镜的尺寸为[10.8 μm（H）×10.8 μm（V）]，镜面之间的间隙为[1 μm]。激光经扩束准直后，照射到DMD上，全息图通过计算机加载到DMD上，分别在离DMD距离[21?23]为3.5 m和5.85 m处获得如图6所示的再现像。

当接收屏与DMD相距3.5 m时，皇后的像清晰可见而国王模糊不清，如图6（a）所示；当相距5.85 m时，远处的国王像清晰可见而近处的皇后模糊不清，如图6（b）所示，不同衍射距离得到相应位置物体的清晰再现像。

Laser：532 nm的半导体绿色激光器；BE：扩束准直装置；

DMD：数字微镜器件； Camera：数码相机

图5 DMD全息显示系统

图6 DMD光电再现

3 结 语

本文对主菲涅尔波带法进行改进，将二维全息图信息变换到一维空间分别进行运算合成以提高计算效率，依据再现器件要求进行尺寸匹配，得到适合的计算全息图。实验结果表明：本文方法具有快速高效生成三维物体全息图的优点，再现图像质量较高，对于三维物体动态显示具有很好的应用前景。

注：本文通讯作者为顾济华。

参考文献

[1] SON J Y， JAVIDI B， KWACK K D. Methods for displaying three dimensional images [J] . IEEE， 2006， 94（3） ： 502?523.

[2] PASTOOR S， W?PKING M. 3?D displays： A review of current technologies [J]. Displays， 1997， 17（2）： 100?110.

[3] 姜太平，沈春林，谭皓.真三维立体显示技术[J].中国图象图形学报，2003（4）：361? 366.

[4] 李克彬，李世其.3D显示技术的最新研究进展[J].计算机工程， 2003，29（12）：3?4.

[5] 田丰，廖薇，刘锦高.基于 DLP 的自由立体显示系统构建[J].电视技术，2010（3）：31?34.

[6] 于瀛洁，蔡明义，张之江.柱透镜光栅自由立体显示中几何参数间关系分析[J].光电子技术，2010，30（1）：10?15.

[7] 黄涛，朱秋东，王涌天，等.立体显示中的垂轴放大率[J].光子学报，2010，39（1）：116?118.

[8] 周磊，陶宇虹，王琼华，等.立体显示用柱面透镜光栅的设计[J].光子学报，2009，38（1）：30?33.

[9] KLEINBERGER P， KLEINBERGER I. Systems for three?dimensional viewing and projection： U.S， 6252707 [P]. 2001?06?26.

[10] EZRA D， WOODGATE G J， OMAR B A， et al. New autostereoscopic display system [C]// IS&T/SPIE's Symposium on Electronic Imaging： Science & Technology. [S.l.]： International Society for Optics and Photonics， 1995： 31?40.

[11] PASTOOR S， LIU J， RENAULT S. An experimental multimedia system allowing 3?D visualization and eye?controlled interaction without user?worn devices [J]. IEEE Transactions on Multimedia， 1999， 1（1）： 41?52.

[12] SANG X， FAN F C， JIANG C C， et al. Demonstration of a large?size real?time full?color three?dimensional display [J]. Optics Letters， 2009， 34（24）： 3803?3805.

[13] 曾吉勇.折反射全景立体成像[D].成都：四川大学，2003.

[14] 徐韡，刘向东，刘旭，等.基于 3DS MAX 信息源的体三维显示技术[J].浙江大学学报：工学版，2006，39（11）：1723?1726.

[15] 邢建芳，龚华军，沈春林，等.基于数字微镜和旋转扫描技术的体三维显示器[J].光电子·激光，2008，19（8）：1011?1015.

[16] 李俊昌.彩色数字全息波前重建算法概论[J].中国激光， 2011，38（5）：1?8.

[17] 杨上供，杜旭日，陈海云，等.计算机与光学方法相结合的彩色彩虹全息术[J].光电子·激光，2009，20（1）：471?474.

[18] SAKAMOTO Y， AOYAMA T， AOKI Y. Volume rendering for computer generated holograms [C]// 19th Congress of the International Commission for Optics： Optics for the Quality of Life. International Society for Optics and Photonics. [S.l.]： [s.n.]， 2003： 556?557.

[19] WANG H， LI Y， JIN H， et al. Three?dimensional visualization of shape measurement data based on a computer generated hologram [J]. Journal of Optics A： Pure and Applied Optics， 2003， 5（5）： S195.

[20] WATERS J P. Holographic image synthesis utilizing theoretical methods [J]. Applied physics Letters， 2004， 9（11）： 405?407.

[21] LESEBERG D. Computer?generated three?dimensional image holograms [J]. Applied optics， 1992， 31（2）： 223?229.

[22] 张晓洁，刘旭，陈晓西.利用菲涅尔波带法计算三维全息[J].光电工程，2004，31（12）：58?60.

[23] ZHANG X， LIU X， CHEN X. Computer?generated holograms for 3D objects using the Fresnel zone plate [C]// Photonics Asia 2004. [S.l.]： International Society for Optics and Photonics， 2005： 109?115.

[24] 陈晓西，刘旭，林远芳，等.三维物体的计算全息的快速算法[J].浙江大学学报：工学版，2009，43（1）：54?56.

[25] 李勇，许富洋，金洪震，等.一种菲涅尔全息图的快速算法[J].光子学报，2010，39（3）：529?531.

（上接第115页）

[26] 陈慧蓉，付胜豪，王元庆，等.一种计算机制全息图快速运算算法[J].计算机工程与应用，2013，49（18）：142?144.

[27] 张晓洁.基于LCD空间光调制器全息显示的若干问题研究[D].杭州：浙江大学，2007.

[28] 李晨.基于DMD的3D数字全息显示系统研究[D].苏州：苏州大学，2012.

[29] 高燕.DMD数字全息及彩色数字全息技术的研究[D].苏州：苏州大学，2013.

[30] KREIS T， ASWENDT P， HOFLING R. Hologram reconstruction using a digital micromirror device [J]. Optical Engineering， 2001， 40（6）： 926?933.

[19] WANG H， LI Y， JIN H， et al. Three?dimensional visualization of shape measurement data based on a computer generated hologram [J]. Journal of Optics A： Pure and Applied Optics， 2003， 5（5）： S195.

[20] WATERS J P. Holographic image synthesis utilizing theoretical methods [J]. Applied physics Letters， 2004， 9（11）： 405?407.

[21] LESEBERG D. Computer?generated three?dimensional image holograms [J]. Applied optics， 1992， 31（2）： 223?229.

[22] 张晓洁，刘旭，陈晓西.利用菲涅尔波带法计算三维全息[J].光电工程，2004，31（12）：58?60.

[23] ZHANG X， LIU X， CHEN X. Computer?generated holograms for 3D objects using the Fresnel zone plate [C]// Photonics Asia 2004. [S.l.]： International Society for Optics and Photonics， 2005： 109?115.

[24] 陈晓西，刘旭，林远芳，等.三维物体的计算全息的快速算法[J].浙江大学学报：工学版，2009，43（1）：54?56.

[25] 李勇，许富洋，金洪震，等.一种菲涅尔全息图的快速算法[J].光子学报，2010，39（3）：529?531.

（上接第115页）

[26] 陈慧蓉，付胜豪，王元庆，等.一种计算机制全息图快速运算算法[J].计算机工程与应用，2013，49（18）：142?144.

[27] 张晓洁.基于LCD空间光调制器全息显示的若干问题研究[D].杭州：浙江大学，2007.

[28] 李晨.基于DMD的3D数字全息显示系统研究[D].苏州：苏州大学，2012.

[29] 高燕.DMD数字全息及彩色数字全息技术的研究[D].苏州：苏州大学，2013.

[30] KREIS T， ASWENDT P， HOFLING R. Hologram reconstruction using a digital micromirror device [J]. Optical Engineering， 2001， 40（6）： 926?933.

[19] WANG H， LI Y， JIN H， et al. Three?dimensional visualization of shape measurement data based on a computer generated hologram [J]. Journal of Optics A： Pure and Applied Optics， 2003， 5（5）： S195.

[20] WATERS J P. Holographic image synthesis utilizing theoretical methods [J]. Applied physics Letters， 2004， 9（11）： 405?407.

[21] LESEBERG D. Computer?generated three?dimensional image holograms [J]. Applied optics， 1992， 31（2）： 223?229.

[22] 张晓洁，刘旭，陈晓西.利用菲涅尔波带法计算三维全息[J].光电工程，2004，31（12）：58?60.

[23] ZHANG X， LIU X， CHEN X. Computer?generated holograms for 3D objects using the Fresnel zone plate [C]// Photonics Asia 2004. [S.l.]： International Society for Optics and Photonics， 2005： 109?115.

[24] 陈晓西，刘旭，林远芳，等.三维物体的计算全息的快速算法[J].浙江大学学报：工学版，2009，43（1）：54?56.

[25] 李勇，许富洋，金洪震，等.一种菲涅尔全息图的快速算法[J].光子学报，2010，39（3）：529?531.

（上接第115页）

[26] 陈慧蓉，付胜豪，王元庆，等.一种计算机制全息图快速运算算法[J].计算机工程与应用，2013，49（18）：142?144.

[27] 张晓洁.基于LCD空间光调制器全息显示的若干问题研究[D].杭州：浙江大学，2007.

[28] 李晨.基于DMD的3D数字全息显示系统研究[D].苏州：苏州大学，2012.

[29] 高燕.DMD数字全息及彩色数字全息技术的研究[D].苏州：苏州大学，2013.

[30] KREIS T， ASWENDT P， HOFLING R. Hologram reconstruction using a digital micromirror device [J]. Optical Engineering， 2001， 40（6）： 926?933.