CN115291490B - 基于光学计算的拓展距离全息显示方法 - Google Patents
基于光学计算的拓展距离全息显示方法 Download PDFInfo
- Publication number
- CN115291490B CN115291490B CN202210076577.5A CN202210076577A CN115291490B CN 115291490 B CN115291490 B CN 115291490B CN 202210076577 A CN202210076577 A CN 202210076577A CN 115291490 B CN115291490 B CN 115291490B
- Authority
- CN
- China
- Prior art keywords
- hologram
- diffraction
- distance
- calculation
- optical
- 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
- 238000000034 method Methods 0.000 title claims abstract description 81
- 230000003287 optical effect Effects 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 title claims abstract description 23
- 230000008569 process Effects 0.000 claims abstract description 37
- 230000008901 benefit Effects 0.000 claims abstract description 5
- 230000009466 transformation Effects 0.000 claims abstract description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 9
- 238000007796 conventional method Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
-
- 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/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- 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/16—Processes or apparatus for producing holograms using Fourier transform
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Holo Graphy (AREA)
Abstract
本发明提出一种基于光学计算的扩展距离全息显示方法,该发明针对传统菲尼尔卷积法的计算全息图生成过程中,由于衍射距离增加造成的重建质量急剧下降和由于卷积运算的三个傅里叶变换造成的计算速度慢的问题。由于该方法在全息图生成的过程中减少了傅里叶变换的数值计算,从而能够提升全息图计算生成速度,相比传统方法计算生成速度要快8倍以上。此外,该方法得到的全息图再现时衍射距离较远,与传统方法相比,具有长距离显示的优势。本发明的方法计算生成全息图的速度快、再现质量高,在目前热门的真三维显示中有极大的应用潜力。
Description
技术领域
本发明涉及一种全息显示领域,特别是全息图的生成和再现方法。
背景技术
全息显示能提供给观看者需要的所有信息,因此被公认为最有前景的三维显示技术。但是,目前的全息显示技术依然有一些问题亟待解决。其中之一:传统菲尼尔卷积法虽然广泛用于计算生成全息图,但是它适用于近场衍射和全息重建质量随着再现距离的增加而急剧退化。同时,传统菲尼尔卷积法的全息图生成速度受到卷积操作中三次快速傅里叶变换的限制,这也是另一个亟待解决的问题。然而,也有许多应用场景需要更长的再现距离。
发明内容
本发明要解决的技术问题是提供一种全息图生成快速,且能在扩展的再现距离上实现高质量重建的方法。解决上述技术问题采用如下技术措施:基于光学计算的拓展距离全息显示方法按如下步骤进行: (i)基于光学计算的菲尼尔卷积法计算全息面衍射场分布:首先,对物体U0进行随机相位调制得到复振幅分布U1,其过程表示为U1=U0´exp(i´φ),其中,i为虚数单位,φ为在[0,2π]间分布的随机相位,特别地,如果随机相位的相位皆为0,则是不添加随机相位,不添加随机相位是添加随机相位的一种特殊情况;然后,对于复振幅分布为U1的物平面衍射场经过距离为z的衍射过程,采用基于光学计算的菲尼尔卷积法计算全息面衍射场分布U2,其过程表示为U2=FT(U1)´FT(hFrz),其中,FT为傅里叶变换,hFrz为衍射过程的点扩散函数。 (ii)全息面衍射场分布编码为计算全息图:将得到的全息面衍射场分布U2进行编码,得到用于空间光调制器件上加载并进行显示的全息图Holo,其过程表示为Holo=Encode(U2),其中Encode()表示对复振幅的编码过程函数。 (iii)全息图的再现:首先,将得到的全息图Holo加载到位于光学傅里叶透镜一边焦平面的空间光调制器上;然后,在光学傅里叶透镜的另一边焦平面设置零级光滤波器;最后,再经过距离为z的衍射过程,就得到高质量的光学再现;全息图再现的数字模拟过程表示为U0’=FrT{DCF[FT(Holo)], z},其中,FT表示光学傅里叶变换过程,DCF表示零级光滤波函数,FrT{(), z}表示距离为z的衍射过程。本发明的有益效果在于:由于该方法在全息图生成的过程中减少了傅里叶变换的数值计算,从而能够提升全息图计算生成速度,相比传统方法计算生成速度要快8倍以上;此外,该方法得到的全息图再现时衍射距离较远,与传统方法相比,具有长距离显示的优势。
附图说明
附图1为本发明的全息图生成与再现流程图。
附图2为全息图计算生成时间与传统方法的对比,图2(a)为两种方法在不同全息图分辨率下的计算生成时间对比,2(b)为本发明方法与传统方法的时间比值,即加速倍率。
附图3为全息图再现质量与传统方法的对比,3(a-d)本发明方法和3(e-h)传统方法生成的全息图在不同重建距离下的再现结果。
注:上述附图只是示意性的,并没有按比例绘制。
具体实施方式
下面详细说明本发明基于光学计算的拓展距离全息显示方法的一个典型实施例,对本发明进行进一步的具体描述。有必要在此指出的是,以下实施例只用于本发明做进一步的说明,不能理解为对本发明保护范围的限制,该领域技术熟练人员根据上述本发明内容对本发明做出一些非本质的改进和调整,仍属于本发明的保护范围。
本发明所述方法的具体实施方式如下:
基于光学计算的拓展距离全息显示方法(如图1所示)分如下几个步骤:
(i)基于光学计算的菲尼尔卷积法计算全息面衍射场分布:首先,对物体U0进行随机相位调制得到复振幅分布U1,其过程表示为U1=U0´exp(i´φ),其中,i为虚数单位,φ为在[0,2π]间分布的随机相位,特别地,如果随机相位的相位皆为0,则是不添加随机相位,不添加随机相位是添加随机相位的一种特殊情况;然后,对于复振幅分布为U1的物平面衍射场经过距离为z的衍射过程,采用基于光学计算的菲尼尔卷积法计算全息面衍射场分布U2,其过程表示为U2=FT(U1)´FT(hFrz),其中,FT为傅里叶变换,hFrz为衍射过程的点扩散函数。 (ii)全息面衍射场分布编码为计算全息图:将得到的全息面衍射场分布U2进行编码,得到用于空间光调制器件上加载并进行显示的全息图Holo,其过程表示为Holo=Encode(U2),其中Encode()表示对复振幅的编码过程函数。 (iii)全息图的再现:首先,将得到的全息图Holo加载到位于光学傅里叶透镜一边焦平面的空间光调制器上;然后,在光学傅里叶透镜的另一边焦平面设置零级光滤波器;最后,再经过距离为z的衍射过程,就得到高质量的光学再现;全息图再现的数字模拟过程表示为U0’=FrT{DCF[FT(Holo)], z},其中,FT表示光学傅里叶变换过程,DCF表示零级光滤波函数,FrT{(), z}表示距离为z的衍射过程。
本发明的实例中,基于光学计算的菲尼尔卷积法计算全息面衍射场分布时,衍射过程的点扩散函数hFrz,具体表示为hFrz(x, y)=exp(ikz)×exp{ik[x2+y2]/(2z)}/(iλz),其中,k=2π/λ,z为衍射距离,λ为波长,x和y分别为空域横轴和纵轴坐标。
本发明的实例中,全息衍射场的编码过程函数Encode(),具体为取相位函数。
下面结合实施例和附图对本发明的内容进行进一步的解释:
全息图计算生成过程中,物体的分辨率为1024×1024,像素间距为8.0um,波长λ为671nm,衍射距离z为100、400、900、1650 mm;全息图的分辨率为1024×1024,像素间距为8.0um;计算平台为Intel Core i5 7400M 中央处理器, 8 GB 内存和微软 Windows 10 操作系统,程序语言为python 3.8。本发明方法,OCFC法,的计算生成时间与传统方法,FR-CV法,的计算生成时间对比如图2所示,计算生成时间结果显示,本发明的方法比传统方法快8倍以上。在不同重建距离下的再现结果如图3所示,3(a-d)和3(e-h)分别为本发明方法和传统方法生成的全息图在不同重建距离下的再现结果,从结果对比来看,本发明的方法再现质量较高。
因此,本发明的方法计算生成全息图的速度快、再现质量高,在目前热门的真三维显示中有极大的应用潜力。
Claims (1)
1.基于光学计算的拓展距离全息显示方法,其特征是按如下步骤进行:
(i)基于光学计算的菲尼尔卷积法计算全息面衍射场分布:首先,对物体U0进行随机相位调制得到复振幅分布U1,其过程表示为U1=U0´exp(i´φ),其中,i为虚数单位,φ为在[0,2π]间分布的随机相位,如果随机相位的相位皆为0,则是不添加随机相位,不添加随机相位是添加随机相位的一种特殊情况;然后,对于复振幅分布为U1的物平面衍射场经过距离为z的衍射过程,采用基于光学计算的菲尼尔卷积法计算全息面衍射场分布U2,其过程表示为U2=FT(U1)´FT(hFrz),其中,FT为傅里叶变换,hFrz为衍射过程的点扩散函数,具体表示为hFrz(x, y)=exp(ikz)×exp{ik[x2+y2]/(2z)}/(iλz),其中,k=2π/λ,z为衍射距离,λ为波长,x和y分别为空域横轴和纵轴坐标;
(ii)全息面衍射场分布编码为计算全息图:将得到的全息面衍射场分布U2进行编码,得到用于空间光调制器件上加载并进行显示的全息图Holo,其过程表示为Holo=Encode(U2),其中Encode()表示对复振幅的编码过程函数;
(iii)全息图的再现:首先,将得到的全息图Holo加载到位于光学傅里叶透镜一边焦平面的空间光调制器上;然后,在光学傅里叶透镜的另一边焦平面设置零级光滤波器;最后,再经过距离为z的衍射过程,就得到高质量的光学再现;全息图再现的数字模拟过程表示为U0’=FrT{DCF[FT(Holo)], z},其中,FT表示光学傅里叶变换过程,DCF表示零级光滤波函数,FrT{(), z}表示距离为z的衍射过程;
该方法在全息图生成的过程中减少了傅里叶变换的数值计算,能够提升全息图计算生成速度;该方法得到的全息图再现时衍射距离较远,具有长距离显示的优势。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210076577.5A CN115291490B (zh) | 2022-01-24 | 2022-01-24 | 基于光学计算的拓展距离全息显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210076577.5A CN115291490B (zh) | 2022-01-24 | 2022-01-24 | 基于光学计算的拓展距离全息显示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115291490A CN115291490A (zh) | 2022-11-04 |
CN115291490B true CN115291490B (zh) | 2023-11-03 |
Family
ID=83820855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210076577.5A Active CN115291490B (zh) | 2022-01-24 | 2022-01-24 | 基于光学计算的拓展距离全息显示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115291490B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10301466A (ja) * | 1997-02-27 | 1998-11-13 | Nippon Telegr & Teleph Corp <Ntt> | 計算機ホログラム合成方法および装置およびこの方法を記録した記録媒体 |
US7277209B1 (en) * | 1997-11-20 | 2007-10-02 | European Community (Ec) | Computer-assisted holographic method and device |
JP2011141776A (ja) * | 2010-01-08 | 2011-07-21 | Pegasus Net Kk | 光学的情報読取装置及びホログラム型2次元コード |
CN102419478A (zh) * | 2011-11-30 | 2012-04-18 | 华侨大学 | 长距离近似无衍射光束的产生装置 |
WO2012069810A1 (en) * | 2010-11-22 | 2012-05-31 | Light Blue Optics Ltd | Spatial light modulators |
JP2015036799A (ja) * | 2013-08-15 | 2015-02-23 | 国立大学法人北海道大学 | 複素振幅画像再生装置および複素振幅画像再生方法、ならびに、散乱位相画像生成装置および散乱位相画像生成方法 |
KR20150033501A (ko) * | 2013-09-23 | 2015-04-01 | 한국전자통신연구원 | 광시야각 홀로그래픽 디스플레이 장치 |
CN108519729A (zh) * | 2018-04-24 | 2018-09-11 | 浙江师范大学 | 一种大尺寸高分辨率彩色菲涅尔全息制作方法与显示系统 |
CN109375489A (zh) * | 2018-10-19 | 2019-02-22 | 四川大学 | 一种无透镜全息投影中振铃现象的抑制方法 |
CN109901370A (zh) * | 2019-01-22 | 2019-06-18 | 四川大学 | 纯相位全息图和单随机相位编码的光学图像加解密方法 |
CN110647023A (zh) * | 2019-09-16 | 2020-01-03 | 四川大学 | 一种基于部分角谱法的快速全息图生成和高质量再现方法 |
CN111897197A (zh) * | 2020-08-18 | 2020-11-06 | 四川大学 | 基于双相位编码的傅里叶相位全息图生成方法 |
CN113468554A (zh) * | 2021-06-02 | 2021-10-01 | 北京理工大学 | 基于深度学习的高保真全息轨道角动量通信方法 |
CN113777902A (zh) * | 2021-09-01 | 2021-12-10 | 北京航空航天大学 | 一种基于随机梯度下降算法的曲面全息噪声抑制方法 |
CN113791529A (zh) * | 2021-08-13 | 2021-12-14 | 北京航空航天大学 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITRM20030398A1 (it) * | 2003-08-13 | 2005-02-14 | Consiglio Nazionale Ricerche | Metodo per la modifica della risoluzione spaziale |
US20180188686A1 (en) * | 2016-12-30 | 2018-07-05 | Electronics And Telecommunications Research Institute | Method and apparatus for hologram resolution transformation |
TWI664586B (zh) * | 2018-03-22 | 2019-07-01 | National Taiwan Normal University | 透過瑕疵複數光場之透明基板瑕疵檢測方法 |
US11137719B2 (en) * | 2018-12-11 | 2021-10-05 | University Of North Carolina At Chapel Hill | Methods, systems, and computer readable media for improved digital holography and display incorporating same |
-
2022
- 2022-01-24 CN CN202210076577.5A patent/CN115291490B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10301466A (ja) * | 1997-02-27 | 1998-11-13 | Nippon Telegr & Teleph Corp <Ntt> | 計算機ホログラム合成方法および装置およびこの方法を記録した記録媒体 |
US7277209B1 (en) * | 1997-11-20 | 2007-10-02 | European Community (Ec) | Computer-assisted holographic method and device |
JP2011141776A (ja) * | 2010-01-08 | 2011-07-21 | Pegasus Net Kk | 光学的情報読取装置及びホログラム型2次元コード |
WO2012069810A1 (en) * | 2010-11-22 | 2012-05-31 | Light Blue Optics Ltd | Spatial light modulators |
CN102419478A (zh) * | 2011-11-30 | 2012-04-18 | 华侨大学 | 长距离近似无衍射光束的产生装置 |
JP2015036799A (ja) * | 2013-08-15 | 2015-02-23 | 国立大学法人北海道大学 | 複素振幅画像再生装置および複素振幅画像再生方法、ならびに、散乱位相画像生成装置および散乱位相画像生成方法 |
KR20150033501A (ko) * | 2013-09-23 | 2015-04-01 | 한국전자통신연구원 | 광시야각 홀로그래픽 디스플레이 장치 |
CN108519729A (zh) * | 2018-04-24 | 2018-09-11 | 浙江师范大学 | 一种大尺寸高分辨率彩色菲涅尔全息制作方法与显示系统 |
CN109375489A (zh) * | 2018-10-19 | 2019-02-22 | 四川大学 | 一种无透镜全息投影中振铃现象的抑制方法 |
CN109901370A (zh) * | 2019-01-22 | 2019-06-18 | 四川大学 | 纯相位全息图和单随机相位编码的光学图像加解密方法 |
CN110647023A (zh) * | 2019-09-16 | 2020-01-03 | 四川大学 | 一种基于部分角谱法的快速全息图生成和高质量再现方法 |
CN111897197A (zh) * | 2020-08-18 | 2020-11-06 | 四川大学 | 基于双相位编码的傅里叶相位全息图生成方法 |
CN113468554A (zh) * | 2021-06-02 | 2021-10-01 | 北京理工大学 | 基于深度学习的高保真全息轨道角动量通信方法 |
CN113791529A (zh) * | 2021-08-13 | 2021-12-14 | 北京航空航天大学 | 一种基于衍射模糊成像原理的无串扰全息3d显示方法 |
CN113777902A (zh) * | 2021-09-01 | 2021-12-10 | 北京航空航天大学 | 一种基于随机梯度下降算法的曲面全息噪声抑制方法 |
Non-Patent Citations (3)
Title |
---|
"Autofocusing algorithm for a digital holographic imaging system using convolutional neural networks";Son,K;《Japanese journal of applied physics》;第57卷(第9期);09SB02-1-09SB02-6 * |
"Holographic display with optical computational fresnel convolution to broaden distance";Jun Wang;《Optics Express》;第30卷(第3期);4288-4301 * |
"holographic display with optical computational fresnel convolution to broaden distance";jun wang;《optics express》;4288-4301 * |
Also Published As
Publication number | Publication date |
---|---|
CN115291490A (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | 4K-DMDNet: diffraction model-driven network for 4K computer-generated holography | |
US9541899B2 (en) | Fast generation of pure phase digital holograms | |
US8890931B2 (en) | Fast generation of holograms | |
CN111240177B (zh) | 一种基于分层像素扫描算法的全息散斑噪声抑制方法 | |
Arai et al. | An accelerated hologram calculation using the wavefront recording plane method and wavelet transform | |
US9823623B2 (en) | Conversion of complex holograms to phase holograms | |
CN109856944B (zh) | 一种提高再现像质量的纯相位全息图生成方法 | |
CN115291490B (zh) | 基于光学计算的拓展距离全息显示方法 | |
CN110647023B (zh) | 一种基于部分角谱法的快速全息图生成和高质量再现方法 | |
CN115690252A (zh) | 一种基于卷积神经网络的全息图重建方法及系统 | |
CN115373247B (zh) | 基于双向相位补偿的锥面全息图快速生成方法 | |
Chen et al. | Generation of phase-only holograms based on aliasing reuse and application in holographic see-through display system | |
CN114237002B (zh) | 像面全息图的生成系统 | |
CN109946943B (zh) | 一种基于光束整形抑制散斑噪声的全息显示系统 | |
Han et al. | Spherical crown holographic display with planar SLM | |
JP2012008207A (ja) | ルックアップテーブルと画像の時間的重複性を用いた3次元動画の計算機合成ホログラムの算出方法及びその装置 | |
Shibukawa et al. | Digital image diffusion technique for suppressing interpage crosstalk in holographic data storage | |
KR102546710B1 (ko) | 디지털 홀로그램 표시 장치 및 디지털 홀로그램 영상 재생 방법 | |
KR20230140187A (ko) | 디지털 홀로그램 변환 및 렌더링 방법 그리고 장치 | |
Sang et al. | Optimized real-time complex-amplitude holographic display based on computer-generated holograms | |
JP4028908B2 (ja) | ホログラム作成装置、ホログラム作成方法およびホログラフィ表示装置 | |
Sui et al. | Noise suppression for double-phase hologram by resampled macro-pixel encoding | |
JP3272752B2 (ja) | 白色光再生ホログラム記録装置 | |
CN114995083B (zh) | 基于平面空间光调制器的柱面全息显示方法 | |
US11768463B2 (en) | Digital hologram display apparatus and displaying method of digital holographic image |
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 |