基于频域合成的计算机制半周视彩色彩虹全息 下载: 881次
1 引言
全息显示技术[1-3]是一种最符合人眼观看的显示技术,其中周视全息可以实现物体360°的显示,是一种极具震撼力的全息三维显示方法。早期人们采用光学方法制作柱面全息[4],后来的计算机技术发展促进了计算柱面全息技术[5]的进步。王典民等[6]在1990年借助彩虹全息的原理实现了一种平板型周视彩虹全息技术,采用复杂的光路实现了单色周视彩虹全息的记录和再现。随着全息直写系统的发展,利用计算全息的方法有可能实现平板型周视彩虹全息,钱惠国等[7]在2005年提出了计算机制平板型周视全息,比较仔细地介绍了光场的传输特性,给出了一个简单三维模型的再现结果。毋东元等[8]采取了比较巧妙的参考光设置方法,将计算全息与光学全息相结合来进行周视全息的制作,给出了初步的结果。然而两种方法所实现的全息图均为单色全息,实际光学再现效果并不理想。2014年,李志光等[9]在前人的研究基础上,实现了计算机与光学相结合彩色平板周视全息显示。他们采用计算机与光学技术相结合的方法首先制作三基色环带菲涅耳全息图,之后再进行光学拍摄得到平板周视全息。光学全息拍摄对实验条件要求严格,需要较多的实验设备,操作复杂。平板型周视彩虹全息图进行360°显示时,由于其具有对称结构,眼睛在某个位置观察,不仅可以看到与该视角对应的正确三维像,同时也会看到与该位置对称的观察窗口处的共轭像,影响观看效果[10]。可以采用体全息的方式通过光学平台拍摄[11]周视彩虹全息来去除共轭像干扰。对于平面型的计算全息,很难去除共轭像的干扰。
实现大侧视角的平板型周视全息图,全息图的空间频率极高[12],在目前全息输出系统输出能力有限的条件下,很难实现,为了避免共轭光对观看的影响,本文提出了一种基于频域合成的计算机制半周视彩色彩虹全息算法。首先分析了半周视彩色彩虹全息图的基本原理,总结了半环形观察窗口与频域的对应关系,使用三维物体特定方向的视图序列,通过将三个颜色分量线性插值,并利用二维傅里叶变换将其衍射到频域,组合形成半周视彩色彩虹全息图的物光频谱分布,通过二维傅里叶逆变换将物光频谱衍射到全息面,得到全息面上的物光分布,取其实部,加上一个偏置分量,得到半周视彩色彩虹全息图。所计算的全息图再现的三维彩色显示的效果逼真,颜色绚丽,可以多人同时观赏。通过所提出的算法,避免了光学拍摄的复杂性,可以在笔记本电脑上实现全息图的计算。
2 基本原理
2.1 半周视彩色彩虹全息计算原理
半周视彩色彩虹全息再现原理如
图 1. 半周视彩色彩虹全息图再现原理示意图
Fig. 1. Schematic of the reconstruction of half-circle view-able color rainbow hologram
图 2. 半周视彩色彩虹全息图物光传播原理图
Fig. 2. Principle of object light propagation of half-circle view-able color rainbow hologram
观察窗口为半环状分布,具有轴对称结构,以一个截面来分析物光的传播特性。如
式中
物体上顶点
式中Δ
全息图上的最大空间频率由
根据空域与频域关系可知,当物体宽度为
最终所计算的全息图尺寸
则频域内的坐标
由上述分析可知,观察窗口位置确定后,由观察窗口所决定的频谱范围就确定了。
对于二维的情况,物体尺寸为
设红色分量每一块的分辨率为(
式中(
频域内的分块数分别为:
式中红色分量的分块数是
图 3. 半周视彩色彩虹全息图频谱合成示意图
Fig. 3. Schematic of the frequency synthesization of half-circle view-able color rainbow hologram
三维物体视图投影角度(
根据投影角度算出投影图像,通过线性插值得到分辨率为
式中,
对于狭缝所决定的频域范围外的区域块,局部频谱数据为零。对所有局部频域进行叠加,可以得到整个频域内的分布,红色分量的物光频谱可以表示为:
频域分块决定了可以观察到视图的数量和分辨率。分块尺寸大,则每个视图的分辨率高,总的可观看的视图数量减少,视图数量与视图分辨率的选择需要取舍。假设两个视图之间的视角为Δ
同理可以在频域内合成绿色和蓝色分量的物光频谱
图 4. 频谱合成示意图。(a) Fred(fx,fy);(b) Fgreen(fx,fy);(c) Fblue(fx,fy);(d) F(fx,fy)
Fig. 4. Diagram of frequency synthesization. (a) Fred(fx,fy); (b) Fgreen(fx,fy); (c) Fblue(fx,fy); (d) F(fx,fy)
则全息面上的物光分布可以表示为:
式中ifft2(*)表示对*进行二维傅里叶逆变换,
对于大尺寸、高采样频率的全息图,由于巨大的数据量,直接对
(13)式可以改写为:
式中ifft(*,1)表示对*进行行方向一维傅里叶逆变换,ifft(*,2)表示对*进行列方向一维傅里叶逆变换。
全息图可以表示为:
式中real(*)表示对*取实部,
通过(15)式可以得到计算机制合成半周视彩色彩虹全息图,之后通过全息输出系统[13-15]进行输出,显影定影漂白后,得到可用于白光再现的全息图。
2.2 三维模型投影算法
三维点云模型数据可以表示为
对于一个确定的视角(
式中
对
设置图像采样间隔为Δ
根据投影关系,某一个点在矩阵中的坐标位置为:
式中round(*)指的是对*取整。
判断该位置的三个通道数据是否都为零。如果不都为零,说明该处已经被之前的物点填充,当前的物点是被遮挡的一个物点。如果都为零,则将当前物点的颜色值赋值给矩阵当前位置。对三维物体所有物点进行计算,可以得到该视角下的投影图像。
3 实验及结果分析
实验中采用的三维模型为通过建模制作的三维人像,投影图像的尺寸设置为47 mm×47 mm,物体距离半环形观察窗口的距离为
图 6. 全息图的(a)局部放大图及(b)显示装置
Fig. 6. (a) Local magnification and (b) display equipment of hologram
4 结论
根据光路可逆原理,分析了半周视彩色彩虹全息的物光频谱分布。通过频域合成的方法得到半周视彩色彩虹全息图的物光频谱,之后逆衍射到全息面得到物光复振幅,然后编码为全息图。所计算的全息图侧视角较小,为15.25°,当全息输出系统进一步发展,能够输出更高频率的全息图时,则可以实现更大侧视角的半周视彩色彩虹全息显示。
[1] 王迪, 唐文华, 王君, 等. 大视角的全息真彩色显示方法[J]. 激光与光电子学进展, 2015, 52(6): 060901.
[2] 吴琼, 王辉, 施逸乐, 等. 彩色反射全息图颜色再现定量分析[J]. 中国激光, 2016, 43(11): 1109001.
[3] 赵锴, 黄应清, 蒋晓瑜, 等. 深度和波长补偿算法生成彩色计算全息图[J]. 中国激光, 2016, 43(6): 0609007.
[6] 王典民, 哈流柱, 王民草. 周视彩虹全息术[J]. 光学学报, 1990, 10(11): 996-998.
[7] 钱惠国, 王辉, 金洪震. 计算机制作平板型周视彩虹全息[J]. 光电子·激光, 2005, 16(3): 372-375.
[8] 毋东元, 王辉, 谢裕涛, 等. 计算机与光学全息术结合制作大侧视角周视全息图[J]. 光子学报, 2013, 42(12): 1454-1459.
[9] 李志光, 吴琼, 王辉, 等. 计算集成平板周视彩色全息图[J]. 中国激光, 2014, 41(11): 1109003.
[10] 李志光. 计算集成平板彩色周视全息图的理论与技术研究[D]. 金华: 浙江师范大学, 2015.
Li ZG. Theoretical and technique research of computer integrated circular-viewing plane color hologram[D]. Jinhua: Zhejiang Normal University, 2015.
[11] 李勇, 王辉, 应朝福, 等. 大视角计算机制彩虹全息图的研究[J]. 光子学报, 2003, 32(5): 595-597.
[12] 拜晓慧, 万霖, 任亚萍, 等. 利用反射全息实现计算全息三维显示[J]. 光子学报, 2012, 41(5): 591-595.
[13] 杨鑫, 李勇, 王辉, 等. 基于频域滤波的计算机制真彩色彩虹全息图模拟再现[J]. 中国激光, 2012, 39(8): 0809001.
[14] 杨上供, 杜旭日, 陈海云, 等. 计算机与光学方法相结合的彩色彩虹全息术[J]. 光电子·激光, 2009, 20(1): 74-77.
[15] 施逸乐, 王辉, 吴琼, 等. 彩色计算全息颜色匹配的研究[J]. 光子学报, 2013, 42(1): 104-109.
Article Outline
杨鑫, 王辉, 李勇, 李志光. 基于频域合成的计算机制半周视彩色彩虹全息[J]. 中国激光, 2018, 45(3): 0309001. Yang Xin, Wang Hui, Li Yong, Li Zhiguang. Computer Generated Half-Circle View-Able Color Rainbow Hologram Based on Frequency Domain Synthesis[J]. Chinese Journal of Lasers, 2018, 45(3): 0309001.