CN113791529B - 一种基于衍射模糊成像原理的无串扰全息3d显示方法 - Google Patents
一种基于衍射模糊成像原理的无串扰全息3d显示方法 Download PDFInfo
- Publication number
- CN113791529B CN113791529B CN202110930000.1A CN202110930000A CN113791529B CN 113791529 B CN113791529 B CN 113791529B CN 202110930000 A CN202110930000 A CN 202110930000A CN 113791529 B CN113791529 B CN 113791529B
- Authority
- CN
- China
- Prior art keywords
- light field
- crosstalk
- diffraction
- wave
- spectrum
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims abstract description 9
- 230000003287 optical effect Effects 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000000644 propagated effect Effects 0.000 claims 1
- 239000004576 sand Substances 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 241000282405 Pongo abelii Species 0.000 description 1
- 230000003796 beauty Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000009916 joint effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
- 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
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N13/00—Stereoscopic video systems; Multi-view video systems; Details thereof
- H04N13/10—Processing, recording or transmission of stereoscopic or multi-view image signals
- H04N13/106—Processing image signals
- H04N13/122—Improving the 3D impression of stereoscopic images by modifying image signal contents, e.g. by filtering or adding monoscopic depth cues
- H04N13/125—Improving the 3D impression of stereoscopic images by modifying image signal contents, e.g. by filtering or adding monoscopic depth cues for crosstalk reduction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N13/00—Stereoscopic video systems; Multi-view video systems; Details thereof
- H04N13/10—Processing, recording or transmission of stereoscopic or multi-view image signals
- H04N13/106—Processing image signals
- H04N13/128—Adjusting depth or disparity
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Signal Processing (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Operations Research (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Holo Graphy (AREA)
Abstract
本发明提出一种基于衍射模糊成像原理的无串扰全息3D显示方法。该方法包括三个步骤:第一步,对于一个3D物体,根据阿贝二次成像理论和菲涅尔衍射原理计算出物体的模糊光场分布,并计算物体的二次衍射模糊成像条件;第二步,基于二次衍射模糊成像特点建立不同深度面光场的串扰关系,计算出串扰光场;第三步,对不同深度面之间的光场串扰,通过叠加光栅相位使串扰光场的空间频谱形成窗口矩阵的特点,使串扰光场以窗口矩阵的形式与目标光场分离,从而生成复振幅全息图,实现无串扰的全息3D显示效果。
Description
一、技术领域
本发明涉及全息3D显示技术,更具体地说,本发明涉及一种基于衍射模糊成像原理的无串扰全息3D显示方法。
二、背景技术
根据全息技术的基本原理,全息显示可分为全息图的记录和再现两个步骤。在全息图的记录过程中,利用了光的干涉原理,以干涉条纹的形式记录物体的振幅和相位信息;在全息图的再现过程中,利用了光的衍射原理,恢复出与被记录物体相同的波前信息,从而提供了人眼视觉里所需要的全部深度信息。因此,全息显示技术被认为是最理想的3D显示技术之一。然而,复杂图像与全深度控制的全息3D显示技术依然难以实现,其根本原因在于当使用2D存储的全息图描绘复杂3D图像所需的所有信息时,不同深度的全息投影图像之间存在相互影响。由于激光具有很高的相干性,在全息重建过程中,单个像点会以艾里斑的形式重建出来,而相邻像点之间则存在一定的重叠区域,这部分区域会发生干涉从而引入串扰光,影响观看效果。虽然国内外学者提出了不少降低串扰的方法,如波前整形法、添加随机相位因子等,但这些方法只能在一定程度上提高全息3D显示的质量,无法完全消除不同深度图像之间的串扰。
三、发明内容
本发明提出一种基于衍射模糊成像原理的无串扰全息3D显示方法。如附图1所示,该方法包括三个步骤:第一步,对于一个3D物体,根据阿贝二次成像理论和菲涅尔衍射原理计算出物体的模糊光场分布,并计算物体的二次衍射模糊成像条件;第二步,基于二次衍射模糊成像特点建立不同深度面光场的串扰关系,计算出串扰光场,得知一个平面对另一个平面的串扰其实是该平面的空间频谱在另一个平面上的二次衍射模糊像;第三步,对不同深度面之间的光场串扰,通过叠加光栅相位使串扰光场的空间频谱形成窗口矩阵的特点,使串扰光场以窗口矩阵的形式与目标光场分离,从而生成复振幅全息图,使目标光场实现无串扰的全息3D显示效果。
在步骤一中,如附图2所示,物波的中心位于坐标原点处,其沿z轴方向进行光场传播,O(ξ,η)代表物波的初始光场分布,接着在物波上叠加焦距为zs的透镜相位信息E(xk,yk,zk)表示物波叠加了透镜之后在距离为zk的衍射光场分布,根据菲涅尔衍射原理,E(xk,yk,zk)与O(ξ,η)的关系为:
其中j表示虚数符号,λ表示波长,表示衍射距离为z=zk的菲涅尔正衍射。当衍射距离为透镜的焦距,即zk=zs时,菲涅尔衍射像是聚焦的;当zk≠zs时,菲涅尔衍射像是离焦的,此时E(xk,yk,zk)为O(ξ,η)的模糊像。
在焦距为zs的透镜相位的作用下,物波衍射到透镜的焦面时得到物波的频谱图像。然后该频谱图像作为新的波源会发生二次衍射。将物波O(ξ,η)的高度记为LH=mdξ,其中,m和dξ分别表示物波在ξ方向的像素数和像素尺寸,根据衍射计算得出二次衍射模糊成像条件为:当z∈[zm,∞)且z≠zs;当时,z∈[zm,z′m]且z≠zs。其中,
当衍射距离z满足上述两个条件中的任意一个条件,则物波的菲涅尔衍射图像就是该物波的模糊像。
其中,Os(-fx,-fy)是光场的空间频谱Os(fx,fy)的倒立像,因此,串扰其实是投影光场的空间频谱叠加上焦距为zs的透镜相位后在衍射距离为z=zk处的菲涅尔衍射光场,即串扰其实是空间频谱光波的二次衍射模糊像。
在步骤三中,为了消除空间频谱光波的二次衍射模糊像对目标光场的影响,在投影光场上添加光栅相位进行卷积,使得目标光场只含高频信号,目标光场用下式表示:
其中,I(xs,ys,zs)表示投影光场的光强分布,dfx和dfy分别表示fx和fy的采样间隔,M和N表示全息图的分辨率,-M≤m≤M,-N≤n≤N。δ表示Dirac函数。用于抵消菲涅尔衍射时在衍射面上产生的二次相位包络,是一个光栅相位,将投影光场中的大量低频信号卷积到高频区域内,从而使其空间频谱形成窗口矩阵的特点。由于投影光场的空间频谱信息被转移到偏离频谱中心的位置,使得中间窗口位置处无空间频谱信息,因此空间频谱光波的二次衍射模糊像主要分布在偏离投影面中心的位置,其不对位于投影面中心的目标光场产生串扰。
对公式(5)进行逆菲涅尔变换,求得重建光场的菲涅尔全息图,最终所得到的全息图复振幅分布H表示为:
四、附图说明
附图1为本发明的一种基于衍射模糊成像原理的无串扰全息3D显示方法流程示意图。
附图2为本发明的物体的模糊成像过程示意图。
附图3为本发明的无串扰全息3D显示的仿真对比结果图。附图3(a)-(b)为本发明的无串扰全息3D显示结果,附图3(c)-(d)为叠加随机相位时的全息3D显示结果。
上述各附图中的图示标号为:
(1)叠加透镜相位后的物波、(2)透镜的焦平面、(3)模糊像。
应该理解上述附图只是示意性的,并没有按比例绘制。
五、具体实施方式
下面详细说明本发明提出的一种基于衍射模糊成像原理的无串扰全息3D显示方法的实施例,对本发明进行进一步的描述。有必要在此指出的是,以下实施例只用于本发明做进一步的说明,不能理解为对本发明保护范围的限制,该领域技术熟练人员根据上述发明内容对本发明做出一些非本质的改进和调整,仍属于本发明的保护范围。
为了实现无串扰的全息3D显示,使用位于两个不同深度平面的两幅图像“美女”、“猩猩”作为被记录物体,其分辨率均为200×200,对应的投影深度分别为7.68cm和18.44cm。接着,将这两幅投影图像的复振幅信息提取出来,并分别叠加透镜相位。令m=M=1000,n=M=1000,设置光波的波长为532nm,叠加光栅相位使两个投影深度面的目标光场与其串扰光场相互分离,根据公式(11)生成分辨率为1000×1000的全息图,全息图的像素大小为dξ=dη=6.4μm。当使用平面波照射该全息图时,其在深度为7.68cm和18.44cm处的仿真重建结果分别如附图3(a)和附图3(b)所示,可以看到,对应深度面的再现像被清晰地重现出来。此时的串扰光场已经完全和目标光场分离开,这说明在光栅相位的作用下,投影图像的空间频谱被改变,从而作为空间频谱模糊像的串扰也会随之改变,最终实现了串扰光场与目标图像的分离。因此,采用本发明所提出的方法能有效地消除串扰的影响,两个平面的平均标准误差值约为0.06。同时,为了进一步说明本发明中消除光场串扰的效果,设置了一组对照组与本发明进行了仿真效果对比。当给每个投影图像施加一个随机相位,即为随机相位时,串扰被离散成了随机散斑,其在深度为7.68cm和18.44cm处的结果如附图3(c)和附图3(d)所示,此时的两个平面的平均标准误差值约为0.44。因此,当为随机相位时,重建光场会受到串扰的影响。因此,本发明的方法能实现高质量的全息3D显示效果。
Claims (1)
1.一种基于衍射模糊成像原理的无串扰全息3D显示方法,其特征在于,该方法包括三个步骤:第一步,对于一个3D物体,根据阿贝二次成像理论和菲涅尔衍射原理计算出物体的模糊光场分布,并计算物体的二次衍射模糊成像条件;第二步,基于二次衍射模糊成像特点建立不同深度面光场的串扰关系,计算出串扰光场,得知一个平面对另一个平面的串扰其实是该平面的空间频谱在另一个平面上的二次衍射模糊像;第三步,对不同深度面之间的光场串扰,通过叠加光栅相位使串扰光场的空间频谱形成窗口矩阵的特点,使串扰光场以窗口矩阵的形式与目标光场分离,从而生成复振幅全息图,使目标光场实现无串扰的全息3D显示效果;
在步骤一中,物波的中心位于坐标原点处,其沿z轴方向进行光场传播,O(ξ,η)代表物波的初始光场分布,接着在物波上叠加焦距为zs的透镜相位信息E(xk,yk,zk)表示物波叠加了透镜之后在距离为zk的衍射光场分布,根据菲涅尔衍射原理,E(xk,yk,zk)与O(ξ,η)的关系为:
其中j表示虚数符号,λ表示波长,表示衍射距离为z=zk的菲涅尔正衍射,当衍射距离为透镜的焦距,即zk=zs时,菲涅尔衍射像是聚焦的;当zk≠zs时,菲涅尔衍射像是离焦的,此时E(xk,yk,zk)为O(ξ,η)的模糊像;
在焦距为zs的透镜相位的作用下,物波衍射到透镜的焦面时得到物波的频谱图像,然后该频谱图像作为新的波源会发生二次衍射,将物波O(ξ,η)的高度记为LH=mdξ,其中,m和dξ分别表示物波在ξ方向的像素数和像素尺寸,根据衍射计算得出二次衍射模糊成像条件为:当时,z∈[zm,∞)且z≠zs;当时,z∈[zm,z′m]且z≠zs;其中,
当衍射距离z满足上述两个条件中的任意一个条件,则物波的菲涅尔衍射图像就是该物波的模糊像;
其中,Os(-fx,-fy)是光场的空间频谱Os(fx,fy)的倒立像,因此,串扰其实是投影光场的空间频谱叠加上焦距为zs的透镜相位后在衍射距离为z=zk处的菲涅尔衍射光场,即串扰其实是空间频谱光波的二次衍射模糊像;
在步骤三中,为了消除空间频谱光波的二次衍射模糊像对目标光场的影响,在投影光场上添加光栅相位进行卷积,使得目标光场只含高频信号,目标光场用下式表示:
其中,I(xs,ys,zs)表示投影光场的光强分布,dfx和dfy分别表示fx和fy的采样间隔,M和N表示全息图的分辨率,-M≤m≤M,-N≤n≤N,δ表示Dirac函数,用于抵消菲涅尔衍射时在衍射面上产生的二次相位包络,是一个光栅相位,将投影光场中的大量低频信号卷积到高频区域内,从而使其空间频谱形成窗口矩阵的特点,由于投影光场的空间频谱信息被转移到偏离频谱中心的位置,使得中间窗口位置处无空间频谱信息,因此空间频谱光波的二次衍射模糊像主要分布在偏离投影面中心的位置,其不对位于投影面中心的目标光场产生串扰;
最终所得到的全息图复振幅分布H表示为:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110930000.1A CN113791529B (zh) | 2021-08-13 | 2021-08-13 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
PCT/CN2021/115786 WO2023015621A1 (zh) | 2021-08-13 | 2021-08-31 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110930000.1A CN113791529B (zh) | 2021-08-13 | 2021-08-13 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113791529A CN113791529A (zh) | 2021-12-14 |
CN113791529B true CN113791529B (zh) | 2022-07-08 |
Family
ID=79181646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110930000.1A Active CN113791529B (zh) | 2021-08-13 | 2021-08-13 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113791529B (zh) |
WO (1) | WO2023015621A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115291490B (zh) * | 2022-01-24 | 2023-11-03 | 四川大学 | 基于光学计算的拓展距离全息显示方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002048800A1 (en) * | 2000-12-13 | 2002-06-20 | Maxsoft Co., Ltd | Phase-encoded multiplexing method using a pseudo-random phase code in a holographic memory system |
CN110363843A (zh) * | 2019-06-18 | 2019-10-22 | 安徽大学 | 基于相移压缩菲涅尔全息的三维图像光学重建方法及系统 |
CN113009710A (zh) * | 2019-12-20 | 2021-06-22 | 杜尔利塔斯有限公司 | 用于在多个平面上形成图像的投影仪 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4421380A (en) * | 1980-10-06 | 1983-12-20 | Mcgrew Stephen P | Full-color hologram |
WO2001023965A1 (en) * | 1999-09-27 | 2001-04-05 | University Of South Florida Division Of Patents And Licencing | Digital interference holographic microscope and methods |
DE102006041875A1 (de) * | 2006-09-01 | 2008-03-13 | Seereal Technologies S.A. | Holographisches Projektionssystem mit Mikrospiegeln zur Lichtmodulation |
US8243353B1 (en) * | 2008-04-07 | 2012-08-14 | Applied Science Innovations, Inc. | Holography-based device, system and method for coded aperture imaging |
TWI655522B (zh) * | 2017-08-08 | 2019-04-01 | 國立臺灣師範大學 | 結構光照明數位全像之方法與裝置 |
CN108519729A (zh) * | 2018-04-24 | 2018-09-11 | 浙江师范大学 | 一种大尺寸高分辨率彩色菲涅尔全息制作方法与显示系统 |
CN111830811A (zh) * | 2020-05-22 | 2020-10-27 | 清华大学 | 基于衍射场叠加的高清三维全息显示方法及其实现装置 |
-
2021
- 2021-08-13 CN CN202110930000.1A patent/CN113791529B/zh active Active
- 2021-08-31 WO PCT/CN2021/115786 patent/WO2023015621A1/zh unknown
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002048800A1 (en) * | 2000-12-13 | 2002-06-20 | Maxsoft Co., Ltd | Phase-encoded multiplexing method using a pseudo-random phase code in a holographic memory system |
CN110363843A (zh) * | 2019-06-18 | 2019-10-22 | 安徽大学 | 基于相移压缩菲涅尔全息的三维图像光学重建方法及系统 |
CN113009710A (zh) * | 2019-12-20 | 2021-06-22 | 杜尔利塔斯有限公司 | 用于在多个平面上形成图像的投影仪 |
Also Published As
Publication number | Publication date |
---|---|
WO2023015621A1 (zh) | 2023-02-16 |
CN113791529A (zh) | 2021-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110363843B (zh) | 基于相移压缩菲涅尔全息的三维图像光学重建方法及系统 | |
JP3471556B2 (ja) | 位相シフトディジタルホログラフィ装置 | |
Soroko | Holography and coherent optics | |
CN113777902B (zh) | 一种基于随机梯度下降算法的曲面全息噪声抑制方法 | |
CN112506019B (zh) | 基于克罗内克积插值的离轴数字全息成像重建方法 | |
CN110032053B (zh) | 一种基于有效利用空间光调制器的全息散斑噪声抑制方法 | |
Zhou et al. | Study on aberration suppressing methods in digital micro-holography | |
Kelly et al. | Digital holographic capture and optoelectronic reconstruction for 3D displays | |
CN103941567A (zh) | 三维数据虚拟分层全息图方法 | |
CN113791529B (zh) | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 | |
CN111240177A (zh) | 一种基于分层像素扫描算法的全息散斑噪声抑制方法 | |
Wang et al. | Zero-order term suppression in off-axis holography based on deep learning method | |
JP2010522350A (ja) | 虚像のないディジタルホログラムの記録・再構成装置及び方法 | |
CN101957171A (zh) | 一种可有效抑制零级和共轭像的同轴数字全息方法 | |
CN117850185B (zh) | 一种基于二元相位调制的离轴数字全息零级项抑制方法 | |
CN117806143A (zh) | 一种基于matlab的离轴数字全息图像再现方法 | |
JP4111614B2 (ja) | 顕微鏡における複素信号検出方法 | |
CN111580372B (zh) | 一种斑点噪声被抑制的大视区计算全息显示方法 | |
CN110018630B (zh) | 一种基于视区特性提高彩色计算全息再现像质量的方法 | |
CN109375489B (zh) | 一种无透镜全息投影中抑制振铃现象的方法 | |
JPWO2012115042A1 (ja) | デジタルホログラフィ装置、及びデジタルホログラフィによる3次元像再生方法 | |
JP6040469B2 (ja) | デジタルホログラフィ装置 | |
CN111880389A (zh) | 一种消除红外数字全息零级衍射的方法 | |
CN110647023A (zh) | 一种基于部分角谱法的快速全息图生成和高质量再现方法 | |
CN114764220B (zh) | 一种基于离轴数字全息改善散斑自相关重建效果的方法 |
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 |