搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

面向纯相位型全息显示的自适应混合约束迭代算法

高乾程 何泽浩 刘珂瑄 韩超 曹良才

引用本文:
Citation:

面向纯相位型全息显示的自适应混合约束迭代算法

高乾程, 何泽浩, 刘珂瑄, 韩超, 曹良才

Adaptive mixed-constraint Gerchberg-Saxton algorithm for phase-only holographic display

Gao Qian-Cheng, He Ze-Hao, Liu Ke-Xuan, Han Chao, Cao Liang-Cai
PDF
HTML
导出引用
  • 纯相位型计算全息图舍弃全息面上的振幅信息, 降低了全息显示的图像质量. 迭代法通过多次正、逆向的波前传播计算, 优化了纯相位型计算全息图的相位轮廓, 提升了全息再现的图像质量. 然而, 传统迭代法可能存在迭代发散、收敛速度偏慢以及再现质量受限等问题. 本文提出了一种基于角谱传播模型的自适应混合约束迭代算法, 通过在迭代中加入自适应的反馈机制和频带约束策略, 增加了迭代优化的自由度, 提升了迭代的收敛性和全息再现的图像质量. 仿真和实验结果表明, 提出的算法能够在较短的计算时间内获得更高重建质量的纯相位型计算全息图, 对于高质量的全息显示具有重要应用价值.
    At present, spatial light modulators are incapable of modulating both the amplitude and phase of the wavefront simultaneously. Therefore, when a spatial light modulator is used for holographic display, it is necessary to encode the complex amplitude of the object wave into amplitude-only or phase-only computer-generated-hologram. The phase-only holographic display has attracted much attention of researchers due to its characteristics of high diffraction efficiency and no conjugate image. However, current optimization algorithms for generating phase-only hologram have the problems of iterative divergence, slow convergence speed, and poor reconstruction quality, which is difficult to satisfy the requirements for high-quality holographic display. In this work, we propose an accurate adaptive mixed constraint Gerchberg-Saxton algorithm by constraining the frequency bandwidth in the hologram plane and adaptively constraining the amplitude of the reconstructed image in the image plane based on the angular spectrum propagation theory. Firstly, we use the angular spectrum propagation model without paraxial approximation to simulate the forward and backward propagation of the light wave for ensuring the accuracy of the wavefront propagation. Secondly, dividing the image plane into signal area and noise area according to the spatial distribution of target image, and different adaptive feedback strategies are set up for the two regions based on the optimized effect of the phase-only hologram. The adaptive feedback strategy is established, which can accelerate the convergence speed of the proposed algorithm and enhance the hologram of freedom of the optimization. Finally, the frequency bandwidth constraint strategy is introduced in the hologram plane to optimize the edge pixels and compensate for the frequency information loss of the phase-only computer-generated hologram, which improves the reconstruction quality of the hologram. After 100 iterations, the average correlation coefficient of the holographic reconstructed image of the proposed algorithm is about 0.9857, and the average peak signal-to-noise ratio is about 31.02 dB. The correlation coefficient and the peak signal-to-noise ratio of the reconstructed images of the proposed algorithm are better than those of the Gerchberg-Saxton algorithm with only frequency bandwidth constraint strategy, and the proposed algorithm has clearer texture and better display effect. The results of numerical simulations and optical experiments show the feasibility and effectiveness of the proposed method. The proposed adaptive mixed constraint Gerchberg-Saxton algorithm is a promising technology for high-quality holographic display.
      通信作者: 韩超, hanchaozh@126.com ; 曹良才, clc@tsinghua.edu.cn
    • 基金项目: 国家自然科学基金区域创新发展联合基金 (批准号: U22A2079 ) 、安徽工程大学检测技术与节能装置安徽省重点实验室开放基金(批准号: DTESD2020A06)、2021年度安徽高校研究生科学研究项目(批准号: YJS20210447)和芜湖市科技计划项目(批准号: 2021cg21)资助的课题.
      Corresponding author: Han Chao, hanchaozh@126.com ; Cao Liang-Cai, clc@tsinghua.edu.cn
    • Funds: Project supported by the Regional Innovation and Development Joint Funds of the National Natural Science Foundation of China (Grant No. U22A2079), the Open Research Fund of Anhui Key Laboratory of Detection Technology and Energy Saving Devices, Anhui Polytechnic University, China (Grant No. DTESD2020A06), the 2021 Anhui University Graduate Scientific Research Project, China (Grant No. YJS20210447), and the Science and Technology Planning Project of Wuhu City, China (Grant No. 2021cg21).
    [1]

    Goodman J W, Lawrence R W 1967 Appl. Phys. Lett. 11 77Google Scholar

    [2]

    Gabor D 1948 Nature 161 777Google Scholar

    [3]

    Zhang H, Zhao Y, Cao L C, Jin G F 2014 Chin. Opt. Lett. 12 060002Google Scholar

    [4]

    He Z H, Sui X M, Jin G F, Cao L C 2019 Appl. Opt. 58 A74Google Scholar

    [5]

    Ma Q G, Cao L C, He Z H, Zhang S D 2019 Chin. Opt. Lett. 17 111001Google Scholar

    [6]

    贾甲, 王涌天, 刘娟, 李昕, 谢敬辉 2012 激光与光电子学进展 49 050002Google Scholar

    Jia J, Wang Y T, Liu J, Li X, Xie J H 2012 Lasers Optoelectron. Prog. 49 050002Google Scholar

    [7]

    Chang C L, Bang Kiseung, Wetzstein Gordon, Lee Byoungho, Gao L 2020 Optica 7 1563Google Scholar

    [8]

    Zhang H, Xie J H, Liu J, Wang Y T 2009 Appl. Opt. 48 5834Google Scholar

    [9]

    王迪, 侯页好, 黄倩, 郑义微, 王琼华 2022 中国激光 49 1909001Google Scholar

    Wang D, Hou Y H, Huang Q, Zheng Y W, Wang Q H 2022 Chin. J. Lasers 49 1909001Google Scholar

    [10]

    王一同, 周宏强, 闫景逍, 合聪, 黄玲玲 2021 中国激光 48 1918004

    Wang Y T, Zhou H Q, Yan J X, He C, Huang L L 2021 Chin. J. Lasers 48 1918004

    [11]

    范爽, 张亚萍, 王帆, 高云龙, 钱晓凡, 张永安, 许蔚, 曹良才 2018 物理学报 67 094203Google Scholar

    Fan S, Zhang Y P, Wang F, Gao Y L, Qian X F, Zhang Y A, Xu W, Cao L C 2018 Acta Phys. Sin. 67 094203Google Scholar

    [12]

    席思星, 于娜娜, 王晓雷, 朱巧芬, 董昭, 王微, 刘秀红, 王华英 2019 物理学报 68 110502Google Scholar

    Xi S X, Yu N N, Wang X L, Zhu Q F, Dong Z, Wang W, Liu X H, Wang H Y 2019 Acta Phys. Sin. 68 110502Google Scholar

    [13]

    曹良才, 何泽浩, 刘珂瑄, 隋晓萌 2022 红外与激光工程 51 20210935Google Scholar

    Cao L C, He Z H, Liu K X, Sui X M 2022 Inf. Laser. Eng. 51 20210935Google Scholar

    [14]

    何泽浩, 隋晓萌, 曹良才, 金国藩 2021 中国激光 48 1209002Google Scholar

    He Z H, Sui X M, Cao L C, Jin G F 2021 Chin. J. Lasers 48 1209002Google Scholar

    [15]

    He Z H, Sui X M, Jin G F, Chu D P, Cao L C 2021 Opt. Express 29 119Google Scholar

    [16]

    Liu K X, He Z H, Cao L C 2021 Chin. Opt. Lett. 19 50501Google Scholar

    [17]

    Sui X M, He Z H, Jin G F, Chu D P, Cao L C 2021 Opt. Express 29 2597Google Scholar

    [18]

    Liu K, He Z H, Cao L C 2022 Appl. Phys. Lett. 120 061103Google Scholar

    [19]

    Wu J C, Liu K X, Sui X M, Cao L C 2021 Opt. Lett. 46 2908Google Scholar

    [20]

    Gerchberg R W, Saxton W O 1972 Optik 35 237

    [21]

    Zhou P C, Bi Y, Sun M Y, Wang H, Li F, Qi Y 2014 Appl. Opt. 53 G209Google Scholar

    [22]

    Chang C L, Xia J, Yang L, Lei W, Yang Z M, Chen J H 2015 Appl. Opt. 54 6994Google Scholar

    [23]

    Chen L Z, Zhang H, He Z H, Wang X Y, Cao L C, Jin G F 2020 Appl. Sci. 10 3652Google Scholar

    [24]

    Wu Y, Wang J, Chen C, Liu C J, Jin F M, Chen N 2021 Opt. Express 29 1412Google Scholar

  • 图 1  AMCGS算法中的信号区与噪声区

    Fig. 1.  Signal area and noise area on the reconstruction plane in the AMCGS algorithm.

    图 2  AMCGS算法原理图

    Fig. 2.  Schematic of AMCGS algorithm.

    图 3  目标图像(a)—(d)及其AMCGS反馈策略仿真重建结果(e)—(h) (a), (e) 狒狒; (b), (f) 船; (c), (g) 芭芭拉; (d), (h) 女人

    Fig. 3.  Target images (a)–(d) and corresponding simulation results by the proposed AMCGS feedback (e)–(h): (a), (e) Baboon;(b), (f) ship; (c), (g) barbara; (d), (h) women.

    图 4  AMCGS算法与RSGS算法仿真重建结果CC值随迭代次数的变化曲线 (a) 狒狒; (b) 船; (c) 芭芭拉; (d) 女人

    Fig. 4.  CCs under different iteration numbers of the reconstructions by the AMCGS algorithm and RSGS algorithm: (a) Baboon; (b) ship; (c) barbara; (d) women.

    图 5  AMCGS算法与RSGS算法仿真重建结果PSNR值随迭代次数的变化曲线 (a) 狒狒; (b) 船; (c) 芭芭拉; (d) 女人

    Fig. 5.  PSNRs under different iteration numbers of the reconstructions by the AMCGS algorithm and RSGS algorithm: (a) Baboon; (b) ship; (c) barbara; (d) women.

    图 6  全息再现仿真结果对比 (a)—(d) 传统随机相位法; (e)—(h) AMCGS算法

    Fig. 6.  Comparison of holographic reconstructions simulation results: (a)–(d) Traditional random-phase method; (e)–(h) AMCGS algorithm.

    图 7  光学实验系统结构

    Fig. 7.  Schematic of optical experimental setup.

    图 8  光学实验结果对比 (a)—(d) 传统随机相位法; (e)—(h) AMCGS算法

    Fig. 8.  Comparison of optics experimental results: (a)–(d) Conventional random phase method; (e)–(h) AMCGS algorithm.

    表 1  传统GS算法, RSGS算法以及AMCGS算法迭代100次所需时间

    Table 1.  Time required for 100 iterations of GS algorithm, RSGS algorithm and AMCGS algorithm.

    算法类型时间/s
    BaboonBarbaraShipWomen
    传统GS算法3.022.982.963.12
    RSGS算法10.0510.0310.0610.10
    AMCGS算法10.4710.1110.0810.26
    下载: 导出CSV
  • [1]

    Goodman J W, Lawrence R W 1967 Appl. Phys. Lett. 11 77Google Scholar

    [2]

    Gabor D 1948 Nature 161 777Google Scholar

    [3]

    Zhang H, Zhao Y, Cao L C, Jin G F 2014 Chin. Opt. Lett. 12 060002Google Scholar

    [4]

    He Z H, Sui X M, Jin G F, Cao L C 2019 Appl. Opt. 58 A74Google Scholar

    [5]

    Ma Q G, Cao L C, He Z H, Zhang S D 2019 Chin. Opt. Lett. 17 111001Google Scholar

    [6]

    贾甲, 王涌天, 刘娟, 李昕, 谢敬辉 2012 激光与光电子学进展 49 050002Google Scholar

    Jia J, Wang Y T, Liu J, Li X, Xie J H 2012 Lasers Optoelectron. Prog. 49 050002Google Scholar

    [7]

    Chang C L, Bang Kiseung, Wetzstein Gordon, Lee Byoungho, Gao L 2020 Optica 7 1563Google Scholar

    [8]

    Zhang H, Xie J H, Liu J, Wang Y T 2009 Appl. Opt. 48 5834Google Scholar

    [9]

    王迪, 侯页好, 黄倩, 郑义微, 王琼华 2022 中国激光 49 1909001Google Scholar

    Wang D, Hou Y H, Huang Q, Zheng Y W, Wang Q H 2022 Chin. J. Lasers 49 1909001Google Scholar

    [10]

    王一同, 周宏强, 闫景逍, 合聪, 黄玲玲 2021 中国激光 48 1918004

    Wang Y T, Zhou H Q, Yan J X, He C, Huang L L 2021 Chin. J. Lasers 48 1918004

    [11]

    范爽, 张亚萍, 王帆, 高云龙, 钱晓凡, 张永安, 许蔚, 曹良才 2018 物理学报 67 094203Google Scholar

    Fan S, Zhang Y P, Wang F, Gao Y L, Qian X F, Zhang Y A, Xu W, Cao L C 2018 Acta Phys. Sin. 67 094203Google Scholar

    [12]

    席思星, 于娜娜, 王晓雷, 朱巧芬, 董昭, 王微, 刘秀红, 王华英 2019 物理学报 68 110502Google Scholar

    Xi S X, Yu N N, Wang X L, Zhu Q F, Dong Z, Wang W, Liu X H, Wang H Y 2019 Acta Phys. Sin. 68 110502Google Scholar

    [13]

    曹良才, 何泽浩, 刘珂瑄, 隋晓萌 2022 红外与激光工程 51 20210935Google Scholar

    Cao L C, He Z H, Liu K X, Sui X M 2022 Inf. Laser. Eng. 51 20210935Google Scholar

    [14]

    何泽浩, 隋晓萌, 曹良才, 金国藩 2021 中国激光 48 1209002Google Scholar

    He Z H, Sui X M, Cao L C, Jin G F 2021 Chin. J. Lasers 48 1209002Google Scholar

    [15]

    He Z H, Sui X M, Jin G F, Chu D P, Cao L C 2021 Opt. Express 29 119Google Scholar

    [16]

    Liu K X, He Z H, Cao L C 2021 Chin. Opt. Lett. 19 50501Google Scholar

    [17]

    Sui X M, He Z H, Jin G F, Chu D P, Cao L C 2021 Opt. Express 29 2597Google Scholar

    [18]

    Liu K, He Z H, Cao L C 2022 Appl. Phys. Lett. 120 061103Google Scholar

    [19]

    Wu J C, Liu K X, Sui X M, Cao L C 2021 Opt. Lett. 46 2908Google Scholar

    [20]

    Gerchberg R W, Saxton W O 1972 Optik 35 237

    [21]

    Zhou P C, Bi Y, Sun M Y, Wang H, Li F, Qi Y 2014 Appl. Opt. 53 G209Google Scholar

    [22]

    Chang C L, Xia J, Yang L, Lei W, Yang Z M, Chen J H 2015 Appl. Opt. 54 6994Google Scholar

    [23]

    Chen L Z, Zhang H, He Z H, Wang X Y, Cao L C, Jin G F 2020 Appl. Sci. 10 3652Google Scholar

    [24]

    Wu Y, Wang J, Chen C, Liu C J, Jin F M, Chen N 2021 Opt. Express 29 1412Google Scholar

计量
  • 文章访问数:  229
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-08-25
  • 修回日期:  2022-10-14
  • 上网日期:  2022-11-11
  • 刊出日期:  2023-01-20

面向纯相位型全息显示的自适应混合约束迭代算法

  • 1. 安徽工程大学, 高端装备先进感知与智能控制教育部重点实验室, 芜湖 241000
  • 2. 清华大学精密仪器系, 精密测试技术及仪器国家重点实验室, 北京 100084
  • 通信作者: 韩超, hanchaozh@126.com ; 曹良才, clc@tsinghua.edu.cn
    基金项目: 国家自然科学基金区域创新发展联合基金 (批准号: U22A2079 ) 、安徽工程大学检测技术与节能装置安徽省重点实验室开放基金(批准号: DTESD2020A06)、2021年度安徽高校研究生科学研究项目(批准号: YJS20210447)和芜湖市科技计划项目(批准号: 2021cg21)资助的课题.

摘要: 纯相位型计算全息图舍弃全息面上的振幅信息, 降低了全息显示的图像质量. 迭代法通过多次正、逆向的波前传播计算, 优化了纯相位型计算全息图的相位轮廓, 提升了全息再现的图像质量. 然而, 传统迭代法可能存在迭代发散、收敛速度偏慢以及再现质量受限等问题. 本文提出了一种基于角谱传播模型的自适应混合约束迭代算法, 通过在迭代中加入自适应的反馈机制和频带约束策略, 增加了迭代优化的自由度, 提升了迭代的收敛性和全息再现的图像质量. 仿真和实验结果表明, 提出的算法能够在较短的计算时间内获得更高重建质量的纯相位型计算全息图, 对于高质量的全息显示具有重要应用价值.

English Abstract

    • 真实世界中的景物通常是具有深度信息的三维物体[1]. 然而, 传统显示技术仅能呈现二维平面图像, 深度信息的缺失限制了显示画面的真实性[2]. 探索三维呈现方式、提升画面展示的真实性, 是显示技术的重要发展方向[3]. 全息三维显示提供了人眼可以感知的所有种类的深度线索[4], 是三维内容的理想呈现方式, 近年来受到广泛关注[5]. 最早出现的光学全息显示技术使用相干条纹实现三维场景的记录和再现[6], 对记录介质、光学器件和实验系统的要求较为严格[7], 且多数情况下只能实现静态场景的三维显示[8]. 伴随着计算机技术发展而产生的计算全息显示技术使用数值计算实现了相干条纹的获取[9], 使用空间光调制器(spatial light modulator, SLM)实现了条纹的加载和重建[10], 极大地降低了全息显示的实现难度[11]. 同时, 可刷新SLM的快速发展[12], 为动态三维全息显示的实现奠定器件基础[13].

      当前, 动态SLM通常难以实现振幅和相位的同时调制. 为了匹配动态SLM的调制特性, 物波函数的复振幅需要被编码为振幅型或相位型的计算全息图(computer-generated hologram, CGH)[14] . 纯相位型CGH能量利用率高、共轭像易于消除, 是动态三维信息的理想承载方式之一[15]. 不过, 纯相位型CGH舍弃了全息面上的振幅信息, 降低了全息三维显示的图像质量. 提升纯相位型CGH显示质量的常用方法包括误差扩散法[16]、双相位编码法[17]、迭代编码法[18]及深度学习法[19]等. 迭代编码法通过多次正、逆向的波前传播计算, 可以优化全息面上的相位轮廓, 进而提升纯相位型CGH的重建质量. Gerchberg-Saxton(GS)算法是一种典型的迭代优化算法[20], 但常面临收敛速度受限的问题. 加权GS方法通过在像平面引入比例反馈因子的方式提升了迭代计算的收敛速度[21], 但部分结果仍存在迭代发散、收敛速度偏慢以及重建质量受限等问题, 需要进一步优化提升.

      本文提出了一种基于角谱传播的自适应混合约束迭代(adaptive mixed constraint Gerchberg-Saxton, AMCGS)算法. 该算法基于自适应的反馈机制和频带约束策略, 在迭代过程中为像平面和全息平面的不同区域设置自适应的反馈因子, 提升了迭代优化过程的收敛速度和稳定性, 在较短的计算时间内获得了更高重建质量的纯相位型CGH, 可以应用在高质量的全息显示场景中.

    • 纯相位型CGH的迭代优化方法需要基于特定的波前传播模型进行设计, 如傅里叶变换模型、菲涅耳衍射模型和角谱传播模型等. 其中, 傅里叶变换模型模拟夫琅禾费衍射, 可用于远场衍射分布的计算, 但对于三维场景的适用性较弱, 需要进行针对性优化; 菲涅耳衍射模型的物面采样间距与传输距离相关, 增加了计算难度. 角谱传播模型可以精确模拟平面波集的光学传播过程, 物面采样间距与传输距离无关. 因此, 本文提出的AMCGS算法基于角谱传播模型设计. 平面波集由全息面传播到像平面的过程定义为正向角谱传播过程, 其数学形式可描述为

      $ U\left( {x,y} \right) = {\Im ^{ - 1}}\left\{ {\Im \left[ {u\left( {x,y} \right)} \right] \cdot H\left( {{f_x},{f_y}} \right)} \right\}, $

      式中, $ \Im $表示傅里叶变换, $ {\Im ^{ - 1}} $表示逆傅里叶变换; $ u\left( {x, y} \right) $表示全息面的复振幅分布, $ U\left( {x, y} \right) $表示像平面的复振幅分布; $ H\left( {{f_x}, {f_y}} \right) $表示与传播距离$ z $相关的传递函数, 其具体形式为

      $\begin{split} & H\left({f}_{x},{f}_{y}\right) \\ =\;& \left\{ \begin{array}{*{20}{l}} &{\rm{exp}}\left(\text{i}\dfrac{2\text{π}}{\lambda }z\sqrt{1-{\lambda }^{2}{f}_{x}^{2}-{\lambda }^{2}{f}_{y}^{2}}\right),\\ & \qquad\;\; {\lambda }^{2}{f}_{x}^{2}+{\lambda }^{2}{f}_{y}^{2} < 1\text{, } \\& 0,\qquad 其他\text{, }\end{array} \right. \end{split}$

      其中, $ \lambda $表示照明波长, $ {f_x} $$ {f_y} $分别表示沿$ x $轴与$ y $轴方向的空间频率. 同理, 平面波集由像平面传播到全息面的过程定义为逆向角谱传播过程.

      为了获得全息面上的复振幅分布, 提取复振幅$ U\left( {x, y} \right) $的相位, 并使用逆向角谱传播模型:

      $ \begin{split}u'\left( {x,y} \right) =\;& {\Im ^{ - 1}}\{ \Im \left\{ {A\left( {x,y} \right) \cdot \arg \left[ {U\left( {x,y} \right)} \right]} \right\} \\ &\times {H_{ - z}}\left( {{f_x},{f_y}} \right) \}, \end{split}$

      式中, $ \arg $表示提取相位, $ A\left( {x, y} \right) $表示像平面上受到约束的振幅, ${H_{{{ - }}z}}\left( {{f_x}, {f_y}} \right)$表示与传播距离$ z $相关的逆向传递函数, $ u'\left( {x, y} \right) $表示经逆向传播获取的全息面复振幅. 通过多次受约束的正、逆向角谱传播, 全息面上的相位轮廓可以得到优化, 进而得到重建质量更高的纯相位型CGH.

    • 传统GS算法的迭代约束对所有采样点赋予相同的权值. 这一反馈策略限制迭代优化的自由度, 降低了迭代过程的收敛速度和稳定性. 在像平面上设置信号区和噪声区, 并为不同区域设置不同的反馈权值, 可以增加迭代优化的自由度[22], 进而加快迭代优化的收敛速度[23]. AMCGS算法对于信号区与噪声区的设置如图1所示. 在迭代过程中, 目标图像需要扩展到原先尺寸的2倍. 由原图像构成的中心区域为信号区, 由采样点扩展获得的边缘区域为噪声区.

      图  1  AMCGS算法中的信号区与噪声区

      Figure 1.  Signal area and noise area on the reconstruction plane in the AMCGS algorithm.

      对于原始采样点数为N×N的目标图像, AMCGS算法的流程如图2所示. 为了按照图1所示方法在像平面构建出信号区和噪声区, 像平面的图像采样点数首先设置为2N×2N [24]. 考虑到角谱传播模型的全息面和像平面采样点数相同, 因此在第一次迭代中, 采样点数为2N×2N的均匀振幅和随机相位构成了全息面的初始复振幅分布. 使用正向角谱传播模型即(1)式, 可获得初始复振幅分布在像平面上对应的重建复振幅分布. 在执行逆向角谱传播前, 按照如下策略对像平面的重建振幅进行约束:

      图  2  AMCGS算法原理图

      Figure 2.  Schematic of AMCGS algorithm.

      $ {A}_{\text{cons}}=\left\{ \begin{array}{*{20}{l}} &2{A}_{\text{T}}-\omega \cdot {A}_{\text{R}},&{A}_{\text{T}}\in S,~{A}_{\text{R}}\in S\text{, }\\ &{A}_{\text{R}},\text{}&{A}_{\text{T}}\in N,~{A}_{\text{R}}\in N\text{, }\end{array}\right.$

      $ \omega = {{\text{e}}^{{A_{\text{T}}} - {A_{\text{R}}} - \ln 2}}, $

      式中, $ S $为信号区, $ N $为噪声区, $ {A_{\text{T}}} $为目标图像的振幅分布, $ {A_{\text{R}}} $为像平面的重建振幅分布, $ {A_{{\text{cons}}}} $为经过自适应混合约束后的振幅, $ \omega $为自适应振幅约束系数. $ \omega $会随着迭代过程变化, 自适应地优化迭代步长以调整迭代进程.

      提取复振幅(1)式的相位, 与约束振幅(4)式共同构成像平面的约束复振幅. 使用逆向角谱传播模型(3)式, 可获得约束复振幅分布在全息面上的复振幅分布. 提取全息面复振幅(3)式的相位, 即可获得采样点数为2N×2N的纯相位型CGH. 由图2可以看出, 噪声区的引入增加了纯相位型CGH的采样点数. 对于原始采样点数为N×N的目标图像, 未经处理的CGH采样点数为2N×2N. 然而, 在全息再现过程中, 上述纯相位型CGH边缘像素对应的重建区域主要为噪声区. 这些边缘像素不仅对目标物体全息再现的贡献较小, 还限制了系统的空间带宽积利用率和光能利用率. 为了避免空间带宽积的浪费, 第二次迭代前需要对纯相位型CGH边缘像素进行优化. AMCGS算法首先基于频带约束策略对采样点数为2N×2N的CGH进行了裁切, 仅保留了其中心区域的N×N个像素; 其次将采样点数为N×N的CGH再次拓展为2N×2N, 以保证其在像平面的重建结果仍旧可以构建出信号区和噪声区. 为了避免光能的浪费, 在上述采样点拓展过程中, 纯相位型CGH的边缘像素取值均为0. 第二次及以后的迭代过程均使用采样点数为2N×2N的均匀振幅和补零拓展后的纯相位型CGH构成全息面的复振幅. 重复自适应迭代过程(1)式—(5)式, 直至达到预设的循环次数, 即可获得目标图像经过优化的纯相位型CGH.

    • 除了可以提升迭代优化的自由度, AMCGS算法还可以拓展角谱传播模型的有效距离. 简便起见, 以一维角谱传播模型为例, 讨论角谱传播模型的有效距离. 此时, 角谱传递函数即(2)式可简化为

      $ H({f}_{x}) = \left\{ \begin{array}{*{20}{l}} &\exp\Big(\text{i}\dfrac{2\text{π}}{\lambda }z\sqrt{1-{\lambda }^{2}{f}_{x}^{2}}\Big),& {\lambda }^{2}{f}_{x}^{2} < 1,\\&0,& 其他.\end{array}\right.$

      提取一维角谱传递函数的相位, 可以得到

      $ \phi \left( {{f_x}} \right) = \frac{{2{\text{π }}}}{\lambda }z\sqrt {1 - {\lambda ^2}f_x^2} . $

      相位函数$ \phi \left( {{f_x}} \right) $的局域空间频率为

      $ {L_f}({f_x}) = \frac{1}{{2{\text{π }}}}\frac{{\partial \phi ({f_x})}}{{\partial {\text{π }}}} = - \frac{{z{f_x}}}{{\sqrt {{\lambda ^{ - 2}} - {f_x}^2} }}. $

      为了避免角谱传播模型中的混叠误差, 一维采样过程中, 频率采样间隔需要满足奈奎斯特采样定理:

      $ \Delta {f_x^{- 1}} = 2{L_{{\text{1D}}}} \geqslant 2\left| {{L_f}\left( {{f_x}} \right)} \right|, $

      式中, Δfx为频率采样间隔, L1D为一维全息平面的总宽度. 考虑到(8)式的单调性, 为了使(9)式始终成立, 仅需保证(9)式在fx取最大值时成立即可.此时, 重建距离z需要满足:

      $ {\left. {{f_x}} \right|_{\max }} = {{{{N_{{\text{1D}}}}} \mathord{\left/ {\vphantom {{{N_{{\text{1D}}}}} {2L}}} \right. } {2L}}_{{\text{1D}}}} \text{, } $

      $ \begin{split}z \leqslant {z_{\max }} =\;& \frac{{{L_{{\text{1D}}}}\sqrt {4{\lambda ^{ - 2}}L_{{\text{1D}}}^2 - N_{{\text{1D}}}^2} }}{{{N_{{\text{1D}}}}}}\\ =\;& {N_{{\text{1D}}}}\Delta {x_{{\text{1D}}}}\sqrt {4{\lambda ^{ - 2}}\Delta x_{{\text{1D}}}^2 - 1} ,\end{split} $

      式中, N1D为一维采样点数, Δx1D为一维采样间隔, zmax为一维角谱传播模型的有效传播距离. AMCGS算法在构建信号区和噪声区的过程中将目标图像采样点数扩展到原始尺寸的2倍, 因此其有效距离也将扩展为传统角谱传播模型的2倍.

    • 本文用峰值信噪比(peak signal-to-noise ratio, PSNR)和相关系数(correlation coefficient, CC)评价全息再现的图像质量. 其中, PSNR (单位: dB)的数学形式为

      $\begin{split}& {\text{PSNR}} = \\ &10 \cdot \lg \left( {\frac{{I_{\max }^2}}{{\displaystyle\frac{1}{{m \times n}}\sum\limits_{i = 0}^{m - 1} {\sum\limits_{j = 0}^{n - 1} {{{\big\| {R( {i,j} ) - K( {i,j} )} \big\|}^2}} } }}} \right), \end{split}$

      式中, $ {I_{\max }} $表示图像的最大灰阶值, RK分别表示目标图像与重建图像的振幅分布. PSNR取值越高, 重建图像的质量相对越好. CC的数学形式为

      $ {\text{CC}} = \frac{\displaystyle{\sum {(X - \overline X )(Y - \overline Y )} }}{\displaystyle{\sqrt {\sum {{{(X - \overline X )}^2}\sum {{{(Y - \overline Y )}^2}} } } }}, $

      式中, $ X $$ \overline X $分别表示目标图像的振幅分布和振幅均值, $ Y $$ \overline Y $分别表示重建图像的振幅分布和振幅均值. CC取值越大, 两幅图像的相关性越强.

      本文采用图3(a)(d)所示的狒狒(baboon), 船(ship), 芭芭拉(barbara), 女人(women)四幅图像作为目标图像. 上述目标图像的原始采样点数均为512×512, 补零填充后采样点数扩展至1024×1024. 图3(e)(h)展示了没有加入频带约束策略时, 不同迭代次数下全息再现结果的PSNR值. 当迭代次数小于20时, 目标图像全息再现结果的PSNR值随着迭代次数的增加快速提升; 当迭代次数大于40时, 目标图像全息再现结果的PSNR值渐渐趋于稳定. 仿真结果表明, 本文设计的反馈策略具有稳定的收敛性, 且全息再现质量相比非迭代算法提升显著.

      图  3  目标图像(a)—(d)及其AMCGS反馈策略仿真重建结果(e)—(h) (a), (e) 狒狒; (b), (f) 船; (c), (g) 芭芭拉; (d), (h) 女人

      Figure 3.  Target images (a)–(d) and corresponding simulation results by the proposed AMCGS feedback (e)–(h): (a), (e) Baboon;(b), (f) ship; (c), (g) barbara; (d), (h) women.

      为了提升空间带宽积和光能的利用率, 将频带约束引入上述反馈策略, 得到完整的AMCGS算法. 同时, 将加入频带约束的传统GS算法命名为RSGS算法, 作为对照组进行仿真实验. 针对图3(a)(d)的4个目标图像, AMCGS算法和RSGS算法重建结果的CC值随迭代次数变化的曲线, 在图4(a)(d)中给出. 可以看出, 随着迭代次数的增加, AMCGS算法再现结果的CC值比RSGS算法更大, 其再现图像与原始图像的相关性更高. 在迭代100次时, AMCGS算法全息再现结果的CC值平均约为0.9857; 在迭代300次时, AMCGS算法全息再现结果的CC值平均约为0.9934, 此时AMCGS算法已趋近于收敛. AMCGS算法100次迭代得到的全息再现图像已经可以满足人眼视觉质量要求. PSNR作为评价全息再现质量的又一定量指标, 本文基于计算成本与全息再现质量间的平衡, 绘制了AMCGS算法和RSGS算法迭代100次PSNR值的变化曲线, 如图5(a)(d)所示. 在迭代100次时, AMCGS算法全息再现结果的PSNR平均值约为31 dB, 相比RSGS算法平均提升了约5 dB, 得到的全息再现结果已可以应用于高质量全息显示. 上述结果均证实本文提出的AMCGS算法能获得更好的全息再现质量. 大量增加迭代次数以测试AMCGS算法再现结果的PSNR值随迭代次数变化的关系, 迭代2×104次时再现图像的PSNR值为41.47 dB, 迭代105次时再现图像的PSNR值为42.7 dB, 仅提升1.2 dB左右. 2×104次迭代后PSNR值随迭代次数变化曲线趋于收敛, 大幅度增加迭代次数对再现图像PSNR值的提升较为有限.

      图  4  AMCGS算法与RSGS算法仿真重建结果CC值随迭代次数的变化曲线 (a) 狒狒; (b) 船; (c) 芭芭拉; (d) 女人

      Figure 4.  CCs under different iteration numbers of the reconstructions by the AMCGS algorithm and RSGS algorithm: (a) Baboon; (b) ship; (c) barbara; (d) women.

      图  5  AMCGS算法与RSGS算法仿真重建结果PSNR值随迭代次数的变化曲线 (a) 狒狒; (b) 船; (c) 芭芭拉; (d) 女人

      Figure 5.  PSNRs under different iteration numbers of the reconstructions by the AMCGS algorithm and RSGS algorithm: (a) Baboon; (b) ship; (c) barbara; (d) women.

      将传统随机相位法和AMCGS算法分别应用于4个目标图像, 传统随机相位法对应的全息再现结果, 如图6(a)(d)所示, 迭代100次, AMCGS算法对应的全息再现结果如图6(e)(h)所示. 对于不同的目标图像, AMCGS算法均可以取得良好的全息再现质量, 相应的全息再现图像对比度高, 图像的细节清晰可见, 更符合人眼视觉质量要求. 使用CC作为评价指标, AMCGS算法的全息再现结果较传统随机相位法平均提升了约0.45; 使用PSNR作为评价指标, AMCGS算法的全息再现结果较传统随机相位法平均提升了约20 dB.

      图  6  全息再现仿真结果对比 (a)—(d) 传统随机相位法; (e)—(h) AMCGS算法

      Figure 6.  Comparison of holographic reconstructions simulation results: (a)–(d) Traditional random-phase method; (e)–(h) AMCGS algorithm.

      本文仿真实验所采用的计算平台为Matlab2021b, Windows 10操作系统和Core I5-11400H (2.70 GHz). 100次迭代后, AMCGS算法可以达到符合人眼视觉质量的全息再现结果. 基于计算成本与全息再现质量间的平衡, 针对4个目标图像, 采样点数设置为512×512. 传统GS算法, RSGS算法以及AMCGS算法完成迭代所需时间见表1. 从表1可以看出, AMCGS算法相较于无需增加采样数的传统GS算法, 其迭代时间较长, 然而传统GS算法难以用于高质量全息显示场景; AMCGS算法相较于同样增加采样数的RSGS算法, 二者迭代时间相差无几, 结合图5可知, 100次迭代后AMCGS算法全息再现结果的PSNR值相比于RSGS算法平均提升了约5 dB, 且AMCGS算法具有更大的优化空间.

      算法类型时间/s
      BaboonBarbaraShipWomen
      传统GS算法3.022.982.963.12
      RSGS算法10.0510.0310.0610.10
      AMCGS算法10.4710.1110.0810.26

      表 1  传统GS算法, RSGS算法以及AMCGS算法迭代100次所需时间

      Table 1.  Time required for 100 iterations of GS algorithm, RSGS algorithm and AMCGS algorithm.

    • 为了进一步探索AMCGS算法在光学实验中的重建质量, 本文构建了如图7所示的光学系统. 波长为532 nm的相干光束通过衰减片、空间滤波器、准直透镜和起偏器后成为对SLM调制敏感的平面照明波. 纯相位型CGH加载在SLM上, 用以调制平面照明波. 实验中使用的SLM型号为HoloEye Gaea-2Vis, 其像素间距为3.74 μm, 采样点数为3840×2160. 经过纯相位型CGH调制的重建光波经过4f系统滤波后, 由Cannon 500D相机采集.

      图  7  光学实验系统结构

      Figure 7.  Schematic of optical experimental setup.

      光学重建结果如图8所示. 对于不同的目标物体, AMCGS算法的光学重建结果的对比度较好、散斑噪声较弱且亮度分布均匀, 与数值仿真中取得的结果相似. 实验重建结果证明了AMCGS算法的实际有效性.

      图  8  光学实验结果对比 (a)—(d) 传统随机相位法; (e)—(h) AMCGS算法

      Figure 8.  Comparison of optics experimental results: (a)–(d) Conventional random phase method; (e)–(h) AMCGS algorithm.

    • 本文提出了一种基于角谱传播模型的AMCGS算法. 该算法将像平面分割为信号区和噪声区, 通过在不同区域加入不同的自适应反馈机制, 提升了迭代优化的稳定性和收敛速度; 在全息面使用频带约束策略, 提升了空间带宽积和光能的利用率, 提高了全息再现的图像质量. 实验表明, 针对不同的目标物体, AMCGS算法全息再现结果的PSNR值相比RSGS算法平均提升了约5 dB, CC值相比RSGS算法平均提升了约0.025, 且重建图像具有更高的对比度、更低的散斑噪声和更均匀的亮度分布. AMCGS算法对于高质量的全息显示具有重要的应用价值.

参考文献 (24)

目录

    /

    返回文章
    返回