CN112100562A - 基于相位补偿的可见光波段球面计算全息图快速生成方法 - Google Patents
基于相位补偿的可见光波段球面计算全息图快速生成方法 Download PDFInfo
- Publication number
- CN112100562A CN112100562A CN202010835995.9A CN202010835995A CN112100562A CN 112100562 A CN112100562 A CN 112100562A CN 202010835995 A CN202010835995 A CN 202010835995A CN 112100562 A CN112100562 A CN 112100562A
- Authority
- CN
- China
- Prior art keywords
- calculation
- spherical
- diffraction
- hologram
- phase compensation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 86
- 238000004364 calculation method Methods 0.000 title claims abstract description 58
- 230000003287 optical effect Effects 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 abstract description 12
- 238000001093 holography Methods 0.000 abstract description 4
- 238000013459 approach Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
- G03H2001/0825—Numerical processing in hologram space, e.g. combination of the CGH [computer generated hologram] with a numerical optical element
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Holo Graphy (AREA)
Abstract
本发明提出一种基于相位补偿的可见光波段球面计算全息图快速生成方法。该方法将球面计算全息图的生成分成两步,首先计算物平面到波前记录面的衍射过程,再采用相位补偿法,计算得到与波前记录面相截的球面计算全息面的衍射场,即可得到球面计算全息图。采用相位补偿法大大降低了采样定理要求的采样点数,从而实现了可见光波段球面计算全息图的快速生成。该方法相比传统点源法,计算生成速度提高了近万倍;与卷积法或球面波谱法相比,需要的采样点数降低了近万倍;该方法重建图像质量高;提出的方法首次解决了平面到球面的衍射计算问题,同时为球面计算全息的光学再现效果验证提供了可行的实践途径。
Description
技术领域
本发明涉及一种全息显示技术,特别是球面全息图的计算生成方法。
背景技术
全息显示作为一种最理想的真三维显示技术,一直受到极大的关注。而球面计算全息因其可提高大视场角而成为近期的一个研究热点。但是球面计算全息图的生成时间耗时巨大,最新的球面计算全息图生成技术采用卷积法或者球谐波谱法大大提高了其计算生成速度,但是这些方法在可见光波长下的满足采样定理的采样点数十分庞大,达到4×1010之大,计算时间也接近1分钟。因而为方便计算,通常将波长限制在远红外波谱,对应的采样点数大大减少,可采用正常的512×256的采样点数,其计算速度大约几十毫秒。波长的局限性制约了球面全息的应用与发展。因此,实现可见光波段下正常采样点数的球面计算全息图的快速生成方法是一个亟待解决的研究问题。
发明内容
本发明针对上述球面计算全息图的生成速度问题,特别是在可见光波段,由于采样点数庞大造成计算速度慢的技术问题,提出一种基于相位补偿的可见光波段球面计算全息图快速生成方法。该方法包括球面计算全息图的计算生成过程和重建过程两个部分。
所述的计算生成过程如图1所示,具体描述为:
步骤一,先计算从为平面的物面到与球面相截的波前记录平面的衍射过程,记物面衍射场为Uo(xo,yo,z1),波前记录平面的衍射场为Uw(xw,yw,z2),其衍射过程表示为Uw=DT(Uo),其中DT(·)表示面到面的衍射计算过程。
步骤二,对波前记录面的衍射场的振幅保持不变,对其相位进行补偿,得到目标球面的衍射场,即球面计算全息图,记为Uh(xh,yh,zh)或Uh(θh,φh,R),R为球面半径,其过程表示为:Uh=Uw×exp(-j×k×Dpc),其中j为虚数单位,波数k=2π/λ,Dpc表示相位补偿角,计算过程为:Dpc=R×cosθh×cosφh-z2,其中sinθh=xh/(R×cosθh),sinφh=yh/R。
所述的重建过程为计算生成过程的逆过程,数值计算过程具体描述为:
步骤一,采用相位补偿法计算全息面到波前记录面的衍射场,表示为:Uw=Uh×exp(j×k×Dpc);步骤二,计算从波前记录面到物面的衍射场,其过程表示为Uo=IDT(Uw),其中IDT(·)表示面到面逆衍射计算过程。
光学过程则直接将球面计算全息图加载到特定的空间光调制器上,在物平面上可得到重建的物体像。
该方法的有益效果在于:采用相位补偿法大大降低了采样定理要求的采样点数,从而实现了可见光波段球面计算全息图的快速生成。该方法相比传统点源法,计算生成速度提高了近万倍;与卷积法或球面波谱法相比,需要的采样点数降低了近万倍;该方法重建图像质量高;提出的方法首次解决了平面到球面的衍射计算问题,同时为球面计算全息的光学再现效果验证提供了可行的实践途径。
附图说明
附图1为本发明的球面计算全息图的生成方法示意图。
附图2为本发明的模拟与光学实验采用的测试图像。
附图3为本发明的模拟实验结果。
附图4为本发明的光学实验结果。
具体实施方式
下面详细说明本发明一种基于相位补偿的可见光波段球面计算全息图快速生成方法的一个典型实施例,对该方法进行进一步的具体描述。有必要在此指出的是,以下实施例只用于该方法做进一步的说明,不能理解为对该方法保护范围的限制,该领域技术熟练人员根据上述该方法内容对该方法做出一些非本质的改进和调整,仍属于本发明的保护范围。
本发明提出一种基于相位补偿的可见光波段球面计算全息图快速生成方法,该方法包括球面计算全息图的计算生成过程和重建过程两个部分。
所述的计算生成过程如图1所示,具体描述为:
步骤一,先计算从为平面的物面到与球面相截的波前记录平面的衍射过程,记物面衍射场为Uo(xo,yo,z1),波前记录平面WRP的衍射场为Uw(xw,yw,z2),其衍射过程表示为Uw=DT(Uo),其中DT(·)表示面到面的衍射计算过程。
步骤二,对波前记录面的衍射场的振幅保持不变,对其相位进行补偿,得到目标球面的衍射场,即球面计算全息图,记为Uh(xh,yh,zh)或Uh(θh,φh,R),R为球面半径,其过程表示为:Uh=Uw×exp(-j×k×Dpc),其中j为虚数单位,波数k=2π/λ,Dpc表示相位补偿角,计算过程为:Dpc=R×cosθh×cosφh-z2,其中sinθh=xh/(R×cosθh),sinφh=yh/R。
所述的重建过程为计算生成过程的逆过程,数值计算过程具体描述为:
步骤一,采用相位补偿法计算全息面到波前记录面的衍射场,表示为:Uw=Uh×exp(j×k×Dpc);步骤二,计算从波前记录面到物面的衍射场,其过程表示为Uo=IDT(Uw),其中IDT(·)表示面到面逆衍射计算过程。
光学过程则直接将球面计算全息图加载到特定的空间光调制器上,在物平面上可得到重建的物体像。
所述的面到面的衍射计算过程可采用但不限于角谱衍射计算公式或者菲尼尔衍射计算公式。
在本发明的实例中,采用256*256的灰度图如图2所示。波长λ和波前记录面边长Lo分别为532nm和5.3mm。图3,图4分别为模拟结果和光学结果。图3(a),(c),(e)和图3(b),(d),(f)分别为提出的方法生成的全息图和采用点源法模拟重建的结果。图4为不同球面张角下的光学重建结果。结果表明,该方法可以有效可行,且重建质量高。
Claims (2)
1.基于相位补偿的可见光波段球面计算全息图快速生成方法,其特征在于,该方法包括球面计算全息图的计算生成过程和重建过程两个部分;所述的计算生成过程具体描述为:步骤一,先计算从为平面的物面到与球面相截的波前记录平面的衍射过程,记物面衍射场为Uo(xo,yo,z1),波前记录平面的衍射场为Uw(xw,yw,z2),其衍射过程表示为Uw=DT(Uo),其中DT(·)表示面到面的衍射计算过程;步骤二,对波前记录面的衍射场的振幅保持不变,对其相位进行补偿,得到目标球面的衍射场,即球面计算全息图,记为Uh(xh,yh,zh)或Uh(θh,φh,R),R为球面半径,其过程表示为:Uh=Uw×exp(-j×k×Dpc),其中j为虚数单位,波数k=2π/λ,Dpc表示相位补偿角,计算过程为:Dpc=R×cosθh×cosφh-z2,其中sinθh=xh/(R×cosθh),sinφh=yh/R;所述的重建过程为计算生成过程的逆过程,数值计算过程具体描述为:步骤一,采用相位补偿法计算全息面到波前记录面的衍射场,表示为:Uw=Uh×exp(j×k×Dpc);步骤二,计算从波前记录面到物面的衍射场,其过程表示为Uo=IDT(Uw),其中IDT(·)表示面到面逆衍射计算过程;光学过程则直接将球面计算全息图加载到特定的空间光调制器上,在物平面上可得到重建的物体像。
2.根据权利要求1所述的面到面的衍射计算过程可采用但不限于角谱衍射计算公式或者菲尼尔衍射计算公式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010835995.9A CN112100562B (zh) | 2020-08-19 | 2020-08-19 | 基于相位补偿的可见光波段球面计算全息图快速生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010835995.9A CN112100562B (zh) | 2020-08-19 | 2020-08-19 | 基于相位补偿的可见光波段球面计算全息图快速生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100562A true CN112100562A (zh) | 2020-12-18 |
CN112100562B CN112100562B (zh) | 2023-09-15 |
Family
ID=73753919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010835995.9A Active CN112100562B (zh) | 2020-08-19 | 2020-08-19 | 基于相位补偿的可见光波段球面计算全息图快速生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100562B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115167087A (zh) * | 2022-05-27 | 2022-10-11 | 四川大学 | 基于随机梯度下降优化算法的球面全息再现质量提升方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010150296A1 (en) * | 2009-06-26 | 2010-12-29 | Consiglio Nazionale Delle Ricerche | Method for controlling depth-of-focus in 3d image reconstructions, in particular for synthesizing three-dimensional dynamic scene for three- dimensional holography display, and holographic apparatus utilizing such a method |
CN104182996A (zh) * | 2014-04-17 | 2014-12-03 | 中国人民解放军装甲兵工程学院 | 一种数字基元全息图的压缩存储与快速恢复方法 |
CN105116705A (zh) * | 2015-09-29 | 2015-12-02 | 北京工业大学 | 多波长自适应数字全息成像系统及方法 |
CN106097333A (zh) * | 2016-06-08 | 2016-11-09 | 上海交通大学 | 微波全息图的生成与重构方法及其系统 |
US20200192287A1 (en) * | 2018-12-11 | 2020-06-18 | The University Of North Carolina At Chapel Hill | Methods, systems, and computer readable media for improved digital holography and display incorporating same |
CN115373247A (zh) * | 2022-07-29 | 2022-11-22 | 四川大学 | 基于双向相位补偿的锥面全息图快速生成方法 |
-
2020
- 2020-08-19 CN CN202010835995.9A patent/CN112100562B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010150296A1 (en) * | 2009-06-26 | 2010-12-29 | Consiglio Nazionale Delle Ricerche | Method for controlling depth-of-focus in 3d image reconstructions, in particular for synthesizing three-dimensional dynamic scene for three- dimensional holography display, and holographic apparatus utilizing such a method |
CN104182996A (zh) * | 2014-04-17 | 2014-12-03 | 中国人民解放军装甲兵工程学院 | 一种数字基元全息图的压缩存储与快速恢复方法 |
CN105116705A (zh) * | 2015-09-29 | 2015-12-02 | 北京工业大学 | 多波长自适应数字全息成像系统及方法 |
CN106097333A (zh) * | 2016-06-08 | 2016-11-09 | 上海交通大学 | 微波全息图的生成与重构方法及其系统 |
US20200192287A1 (en) * | 2018-12-11 | 2020-06-18 | The University Of North Carolina At Chapel Hill | Methods, systems, and computer readable media for improved digital holography and display incorporating same |
CN115373247A (zh) * | 2022-07-29 | 2022-11-22 | 四川大学 | 基于双向相位补偿的锥面全息图快速生成方法 |
Non-Patent Citations (5)
Title |
---|
RUOXUE YANG 等: "Fast Diffraction Calculation for Spherical Computer-Generated Hologram Using Phase Compensation Method in Visible Range" * |
SERGIO DE NICOLA 等: "Wave front reconstruction of Fresnel off-axis holograms with compensation of aberrations by means of phase-shifting digital holography" * |
赵锴 等: "深度和波长补偿算法生成彩色计算全息图" * |
邵晓鹏;刘飞;李伟;杨力铭;杨思原;刘佳维;: "计算成像技术及应用最新进展" * |
马韬: "多层衍射光学元件设计理论及其在混合光学系统中的应用" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115167087A (zh) * | 2022-05-27 | 2022-10-11 | 四川大学 | 基于随机梯度下降优化算法的球面全息再现质量提升方法 |
CN115167087B (zh) * | 2022-05-27 | 2023-12-26 | 四川大学 | 基于随机梯度下降优化算法的球面全息再现质量提升方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112100562B (zh) | 2023-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Acceleration method for computer generated spherical hologram calculation of real objects using graphics processing unit | |
CN112100562B (zh) | 基于相位补偿的可见光波段球面计算全息图快速生成方法 | |
CN101452581B (zh) | 一种计算全息动画的快速生成方法 | |
CN104182996A (zh) | 一种数字基元全息图的压缩存储与快速恢复方法 | |
Ma et al. | A novel feed-forward neural network-based method for fast hologram generation | |
Xu et al. | Calculation of computer-generated hologram (CGH) from 3D object of arbitrary size and viewing angle | |
CN112180707B (zh) | 基于球面自衍射模型的球面纯相位全息图生成方法 | |
Tong et al. | A Deep‐Learning Approach for Low‐Spatial‐Coherence Imaging in Computer‐Generated Holography | |
CN115690252A (zh) | 一种基于卷积神经网络的全息图重建方法及系统 | |
CN110647023B (zh) | 一种基于部分角谱法的快速全息图生成和高质量再现方法 | |
CN115113508A (zh) | 基于光衍射神经网络的全息显示散斑抑制方法 | |
Zhang et al. | Research on three-dimensional computer-generated holographic algorithm based on conformal geometry theory | |
Dong et al. | Vision transformer-based, high-fidelity, computer-generated holography | |
Kang et al. | Quality comparison and acceleration for digital hologram generation method based on segmentation | |
Kovachev et al. | Reconstruction of computer generated holograms by spatial light modulators | |
Qi et al. | Accurate complex modulation by the iterative spatial cross-modulation method | |
Zhao et al. | Fast calculation method for full-color computer-generated hologram of real objects captured by depth camera | |
Bo | Deep learning approach for computer-generated holography | |
Nomura et al. | Electoronic hologram generation using high quality color and depth information of natural scene | |
CN108628138B (zh) | 基于人眼视觉特性的全息头盔显示的处理方法 | |
Xia et al. | Phase-only stereoscopic hologram calculation based on Gerchberg–Saxton iterative algorithm | |
Orzó et al. | Gpu implementation of volume reconstruction and object detection in digital holographic microscopy | |
Hwang et al. | Off-axis complex hologram encoding method for holographic display with amplitude-only modulation | |
Wang et al. | Advances in computer-generated hologram for rendering and acceleration based on the point-polygon hybrid method | |
Zeng et al. | Physically Guided Generative Adversarial Network for Holographic 3D Content Generation from Multi-view Light Field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |