有效视角图像切片嵌合法全息体视图的数值重构 下载: 740次
1 引言
全息体视图打印技术是目前广泛应用的全息技术。该技术将双目视差原理与全息打印技术相结合,先对场景采集的一组二维视角图像进行预处理,而后利用全息打印技术将处理之后的图像记录在全息记录介质中生成全息体视图。对全息体视图进行再现就可观察到记录的三维场景。
全息体视图打印技术经历了近半个世纪的发展,目前已经相当成熟。1969年,Debitetto[1]提出水平视差全息体视图全息打印技术,随后King等[2]针对于Debitetto型水平视差全息体视图提出两步法水平视差全息体视图打印技术,该技术首先制作一个透射式水平视差全息体视图,并将其作为母版,而后将母版上的再现像转印到转移干板上,通过这种方法得到的全息体视图可以得到白光再现的正实像。
在视角图像处理方面,Yamaguchi等[3]提出了一种Lippmann全息体视图处理方法,该方法可一步制作出Lippmann型全息体视图。而后,麻省理工学院(MIT)的Halle等[4-10]提出称为Ultragram的一步法全息体视图打印技术,基于此技术可以实现任意深度、全视差、无畸变的全息体视图打印。Geola公司的Ratcliffe等[11-14]在2003年提出了一步激光直写合成全息体视图打印法。2017年本课题组[15-20]提出了一种称为有效视角图像切片嵌合(EPISM)法,并利用该方法对采样的视角图像进行处理,这种方法的基本原理不同于之前的方法,它是基于全息单元图而非全息单元像素来对视角图像进行处理的,并取得了一系列研究成果。
2012年,Park等[21]提出了针对Halle方法预处理后的图像进行数值重构来模拟全息体视图光学再现像的方法,这种方法针对Ultragram型全息体视图的图像处理思路,给出了相应的数值重构算法,数值重构结果可以很好地模拟光学再现像。2014年,Park等[21-22]进一步针对透镜引入径向畸变的情况给出了相应的数值重构算法,并利用峰值信噪比和结构相似性等图像质量指标对畸变造成的影响进行分析。Park等提出的这种数值重构算法可以在光学实验前对光学再现像进行模拟,这使得在实验前就可以利用模拟的再现像对光学实验结果进行预估,为全息体视图打印提供了很好的辅助。
由于EPISM法是基于全息单元图像而非全息单元像素对视角图像进行处理的,因此需要与之相匹配的算法来实现对EPISM法全息体视图的数值重构。
本文基于EPISM法的基本原理,提出了针对EPISM法的全息体视图数值重构算法,并利用光学实验得到EPISM法的光学再现像,将根据数值重构算法得到的数值重构图像与原始采样视角图像以及光学再现图像进行对比,验证了数值重构算法的有效性,同时对算法存在的问题进行了总结。
2 EPISM法的基本原理
在全息体视图打印中,通常将物光与参考光在全息干板上的孔径区域发生干涉形成的全息单元称为hogel,将两步法全息体视图打印的母版称为H1干板,转印干板称为H2干板。
EPISM法的提出是为了实现通过一步全息打印过程实现两步法的全息打印效果。为了达到该目的,利用在液晶显示器(LCD)上模拟两步法中H2干板上的hogel记录图像的方法,将两步法的记录过程通过一步直接实现。
如
3 EPISM法的数值重构
3.1 EPISM法数值重构的基本原理
光学实验是细致复杂的实验过程,包括实验设计、光路搭建、全息打印、效果观察等多个步骤。因此,在实验前对实验效果进行预模拟会对实验有很大帮助。由于EPISM法采用全息单元图(即有效视角切片)替代全息单元像素作为视角图像预处理的基本单位,因此在这里需要基于EPISM法以及光线追踪法的基本原理给出EPISM法的数值重构方法。
图 1. EPISM法的总体原理图。(a)单个全息单元视角图片有效像素部分的提取;(b)多全息单元有效像素部分的拼接组合
Fig. 1. Overall principle of EPISM method. (a) Extraction of effective perspective image segments of single virtual hogel; (b) synthetic effective perspective image mosaicked by effective image segments of multiple virtual hogels
由于经过EPISM法记录的H2干板上每个hogel所记录的都是该hogel中心所看到的有效视角图像切片嵌合之后得到的合成视角图像,因此再现时在LCD位置还原出的是该hogel记录的合成视角图像。如
由于H2干板再现时还原的是H2干板上每个hogel对应的合成视角图像,因此只在干板制作过程中考虑虚拟H1干板的情况,对EPISM法进行数值重构时只需要考虑实际加载合成视角图像的LCD平面、H2平面以及观察平面的位置关系。
3.2 数值重构的具体算法
如
图 3. 数值重构算法的各个参数(水平方向)
Fig. 3. Parameters for numerical reconstruction (horizontal direction)
3.2.1 有效hogel的确定
在视点
如
图 4. 数值重构时有效hogel的确定
Fig. 4. Determination of effective hogel for numerical reconstruction
若要选择有效hogel,hogel的中心坐标应为
将位于hogel中心左侧的视点记为
将位于hogel中心右侧的视点记为
针对上述关系,hogel要与视点有有效观察图像切片,就需要满足左端点落在
3.2.2 有效观察图像切片的确定以及再现像的数值重构
要确定有效观察图像切片在再现图像上的位置,首先需要确定有效观察图像切片的中心坐标,将观察点
H2的hogel在LCD平面的有效观察图像切片尺寸为
将所有视点
4 实验验证
为了验证所提数值重构算法的正确性,针对茶壶模型利用EPISM法进行光学打印,并与数值重构结果进行比较。
这里采用的茶壶模型长6.4 cm,高3.2 cm,深4 cm,倾斜45°放置。选取H2干板的hogel尺寸为1 cm,H2干板的hogel尺寸为0.25 cm,全息干板的尺寸为8 cm×8 cm,采样相机的视场角为30°,相机到茶壶模型的距离为18.6 cm,H2干板到LCD的距离为11.4 cm。根据EPISM法可以得到,相机采样个数应为91×91=8281个。选取可以完整捕捉到整个茶壶的左下、中心、右上三个视角的采样图片作为参考,如
采样后,利用EPISM法对得到的视角图像切片进行嵌合,得到对应H2干板的8×8=64个hogel的合成图像。如
图 6. EPISM法不同位置的合成视角图像
Fig. 6. Synthetic perspective images from different positions by EPISM method
图 7. 茶壶在三个视角的数值重构图像
Fig. 7. Numerically reconstructed images of teapot from three perspectives
利用所提算法对EPISM法处理过的图像进行数值重构,可以首先对左下、中心、右上的3个位置进行数值重构,结果如
另外利用所提算法,针对EPISM法处理过的合成图像,对指定观察距离和观察角度的再现像进行模拟。在这里选取距离为300 cm,以与水平方向呈-12°、-6°、0°、6°、12°的观察角度(如
通过分析模拟结果发现:选择的视角接近0°时可以获得较好的重构图像,当视角远离0°时,重构图像的质量就会降低;当距离较近时,重构图像的效果并不理想,随着距离增加,重构图像的质量逐渐提高。因此,当选择接近0°视角、距离为300 cm进行数值重构时,可以得到较好的模拟再现像。
图 8. 茶壶在水平方向的数值重构图像
Fig. 8. Numerically reconstructed images of teapot in horizontal direction
图 9. 茶壶在竖直方向上的数值重构图像
Fig. 9. Numerical reconstruction images of teapot in vertical direction
图 10. 茶壶在不同距离时的数值重构图像
Fig. 10. Numerically reconstructed images of teapot at different distances
为了进一步验证实验结果,对EPISM法处理过的合成图像进行全息打印,将数值重构结果与光学再现像进行比较。在这里使用VVX09F035M20型LCD屏加载合成图像,该屏的像素为1920 pixel×1200 pixel,取其中的1000 pixel×1000 pixel作为有效像素部分,对应的有效长度近似为10 cm×10 cm,合成全息立体图打印系统的光路如
如
总体而言,EPISM法数值重构的结果可以较好地对光学再现像进行模拟,从而可以在实验前对实验结果进行预估。
图 12. 数值重构图像与光学再现像的比较
Fig. 12. Comparison between numerically reconstructed images and optically reproduced images
5 结论
本课题组给出了一种重构算法,用以对经EPISM法预处理后得到的全息体视图进行数值重构。通过该重构算法可以模拟EPISM法打印的全息体视图在不同位置、不同角度的全息再现像。通过将数值重构图像与原始采样图像、光学实验得到的再现像进行比较后发现,数值重构图像可以对全息体视图的再现结果进行较好的模拟。采用这种方法可以在实验前对实验的结果进行预估,为光学实验提供参考。由于EPISM法使用了近似替代方法,因此对全息体视图再现像进行模拟时会在细节方面存在些许偏差,无法做到精准的还原,下一步将尝试解决这一问题,实现对EPISM法的光学再现结果进行精准的数值重构。
[3] Yamaguchi M, Ohyama N. Holographic 3-D printer[J]. Proceedings of SPIE, 1990, 1212: 84-92.
[4] Halle M W, Benton S A, Klug M A, et al. The generalized holographic stereogram[J]. Proceedings of SPIE, 1991, 1461: 142-155.
[5] Halle MW. The generalized holographic stereogram[D]. [S.l.]: Massachusetts Institute ofTechnology, 1991.
[6] Klug M A, Halle M W, Hubel P M. Full-color ultragrams[J]. Proceedings of SPIE, 1992, 1667: 110-119.
[11] Ratcliffe DB, Vergnes FM, RodinA, et al. Holographic printer: 7800803[P]. 2003-08-21.
[12] Ratcliffe DB, Vergnes FM, RodinA, et al. Holographic printer:6930811[P]. 2005-08-16.
[13] Ratcliffe DB, RodinA, Hrynkiw L.A method of writing a composite 1-step hologram: US7333252[P].2008-02-19.
[14] Ratcliffe DB, NikolskijA, ZacharovasS, et al. Image capture system for a digital holographic printer: 8154584[P]. 2012-04-10.
[17] 黄应清, 苏健, 陈祎贝, 等. 全息打印技术研究进展[J]. 红外与激光工程, 2018, 47(4): 0406008.
[18] Su J, Yan X P, Huang Y Q, et al. Progress in the Synthetic holographic stereogram printing technique[J]. Applied Sciences, 2018, 8(6): 851-860.
[19] 黄应清, 赵锴, 蒋晓瑜, 等. 用波前平面提高全息体视图成像分辨率[J]. 中国激光, 2016, 43(2): 0209002.
[20] 裴闯, 蒋晓瑜, 闫兴鹏, 等. 采用视角切片法提高全息体视图再现视差图像分辨率[J]. 激光与光电子学进展, 2015, 52(9): 090901.
樊帆, 蒋晓瑜, 王培阳, 陈祎贝, 闫兴鹏. 有效视角图像切片嵌合法全息体视图的数值重构[J]. 中国激光, 2018, 45(12): 1209002. Fan Fan, Xiaoyu Jiang, Peiyang Wang, Yibei Chen, Xingpeng Yan. Numerical Reconstruction of Holographic Stereograms with Effective Perspective Image Segmentation and Mosaicking Method[J]. Chinese Journal of Lasers, 2018, 45(12): 1209002.