CN109856944B - 一种提高再现像质量的纯相位全息图生成方法 - Google Patents
一种提高再现像质量的纯相位全息图生成方法 Download PDFInfo
- Publication number
- CN109856944B CN109856944B CN201910139552.3A CN201910139552A CN109856944B CN 109856944 B CN109856944 B CN 109856944B CN 201910139552 A CN201910139552 A CN 201910139552A CN 109856944 B CN109856944 B CN 109856944B
- Authority
- CN
- China
- Prior art keywords
- hologram
- reproduced image
- image
- phase
- phase hologram
- 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
Images
Landscapes
- Holo Graphy (AREA)
Abstract
本发明针对纯相位全息图生成过程中,由于散斑噪声和混叠现象造成的再现像质量不高的问题,提出一种提高再现像质量的纯相位全息图生成方法。该方法通过栅格化的物体采样,充分利用混叠现象,使得全息面上的振幅分布均匀,生成的纯相位全息图丢失信息减少,从而能够完整地再现原来的物体,有效地提高了纯相位全息图的再现像质量。本发明的方法简单、效果好,在目前热门的AR显示和投影显示中有极大的应用潜力。
Description
技术领域
本发明涉及一种全息投影显示领域,特别是纯相位全息图的生成方法。
背景技术
全息显示能提供给观察者需要的所有信息,因此被公认为最有前景的3D显示技术。但是,目前的全息显示技术依然有一些问题亟待解决,其中之一是由于振幅或相位的丢失造成的再现像的质量退化问题。纯相位全息图的再现像质量可以通过向物体添加随机相位来得到一定的改善,然而干涉产生的散斑噪声也是一个需要解决的问题。造成纯相位全息图再现质量下降的另一个原因是混叠现象,通常混叠是需要消除的,但是通过研究分析发现,如果混叠使用得当,也能用于纯相位全息图的再现像质量提高。
发明内容
本发明要解决的技术问题是提供一种提高再现像质量的纯相位全息图生成方法。
解决上述技术问题采用如下技术措施:提高再现像质量的纯相位全息图生成方法按如下步骤进行:
(i)放大采样间距:对采样间距为p、分辨率为M1×N1的图像持分辨率不变,但是将采样间距扩大s倍,得到图像I1,扩大后的采样间距为p1,p1=s×p。
(ii)模拟漫反射照明:向图像I1代表的物体添加一定的随机噪声,得到待衍射物体复振幅分布U1,其过程表示为U1=I1×exp(i×R1),其中,i为虚数单位,R1在区间[0, tπ]上具有均匀概率分布,t≤2。
(iii)全息面衍射场分布计算:以物体复振幅分布U1作为输入,采用点源法,计算全息面上的衍射场分布U2,其过程表示为:U2=PBM(U1| λ, d),其中PBM()表示点源法的衍射场计算,λ表示波长,d表示衍射距离,全息图的分辨率为M×N,像素间距为p0。
(iv)纯相位全息图生成:对得到的全息面衍射场分布U2取相位角,即得到纯相位的全息图H,其过程表示为:H=ang(U2),其中ang()表示对复数取相位角操作。
关键参数条件:在关键的步骤(i)中,扩大的采样间距p1和像素数量M1必须同时满足3个条件:第一,采样间距大于艾里斑直径,即p1≥2d×tan(1.22λ/p0/M),其中tan()表示正切函数;第二,采样间距小于人眼的分辨能力,即p1≤d1×1.22λ/De,其中d1是人眼观看再现像的距离,De表示人眼的瞳孔直径;第三,空间光调制器的空间带宽积限制,即p1×M1≤d×tan[arcsin(λ/p0)],其中arcsin表示反正弦函数。
本发明的有益效果在于:由于间距扩大后的图像相当于一个光栅的效果,能够使得全息面上的衍射场的混叠效应得到充分利用,可以提高生成的纯相位全息图的再现像质量。
附图说明
附图1为本发明的纯相位全息图生成流程图。
附图2为全息图的再现像质量对比,2(a)为不进行步骤(i-iii)生成的全息图的再现像,2(b)本发明方法生成的全息图的再现像。
注:上述附图只是示意性的,并没有按比例绘制。
具体实施方式
下面详细说明本发明一种提高再现像质量的纯相位全息图生成方法的一个典型实施例,对本发明进行进一步的具体描述。有必要在此指出的是,以下实施例只用于本发明做进一步的说明,不能理解为对本发明保护范围的限制,该领域技术熟练人员根据上述本发明内容对本发明做出一些非本质的改进和调整,仍属于本发明的保护范围。
本发明所述方法的具体实施方式如下:
纯相位全息图的生成过程(如图1所示)分如下几个步骤:
(i)放大采样间距:对采样间距为p、分辨率为M1×N1的图像持分辨率不变,但是将采样间距扩大s倍,得到图像I1,扩大后的采样间距为p1,p1=s×p。
(ii)模拟漫反射照明:向图像I1代表的物体添加一定的随机噪声,得到待衍射物体复振幅分布U1,其过程表示为U1=I1×exp(i×R1),其中,i为虚数单位,R1在区间[0, tπ]上具有均匀概率分布,t≤2。
(iii)全息面衍射场分布计算:以物体复振幅分布U1作为输入,采用点源法,计算全息面上的衍射场分布U2,其过程表示为:U2=PBM(U1| λ, d),其中PBM()表示点源法的衍射场计算,λ表示波长,d表示衍射距离,全息图的分辨率为M×N,像素间距为p0。
(iv)纯相位全息图生成:对得到的全息面衍射场分布U2取相位角,即得到纯相位的全息图H,其过程表示为:H=ang(U2),其中ang()表示对复数取相位角操作。
关键参数条件:在关键的步骤(i)中,扩大的采样间距p1和像素数量M1必须同时满足3个条件:第一,采样间距大于艾里斑直径,即p1≥2d×tan(1.22λ/p0/M),其中tan()表示正切函数;第二,采样间距小于人眼的分辨能力,即p1≤d1×1.22λ/De,其中d1是人眼观看再现像的距离,De表示人眼的瞳孔直径;第三,空间光调制器的空间带宽积限制,即p1×M1≤d×tan[arcsin(λ/p0)],其中arcsin表示反正弦函数。
本发明的实例中,点源法的衍射场计算具体表示为:U2=PBM(U1)=-iλ-1∑∑[U1´exp(ikr)/r],其中,k=2π/λ,r为物点到全息面上某点的距离;步骤(ii)中的参数t设为1。
下面结合实施例和附图对本发明的内容进行进一步的解释:
物体的分辨率M1×N1为256×256,像素间距p1为12×8um;全息图的分辨率M×N为1920×1080,像素间距p为8.0um;波长λ为532nm,衍射距离d为500mm,再现像观看距离d1为500mm,物体分层层数为256,分层层间距为3×104um,对应的采样间距扩大倍率s为12。
为了说明本发明方法的效果,全息图的再现像质量对比结果如图2所示。2(a)为不进行步骤(i-ii)生成的全息图的再现像,由于振幅信息的缺失造成再现像无法恢复整个物体;2(b)本发明方法生成的全息图的再现像,由于光栅化的物体能够充分利用混叠现象,使得全息面的振幅分布均匀,生成的纯相位全息图丢失信息减少,从而能够较完整地再现原来的物体;从对比可以看出再现像的质量有大幅的提高。因此,本发明的方法简单、效果好,在目前热门的AR显示和投影显示中有极大的应用潜力。
Claims (1)
1.一种提高再现像质量的纯相位全息图生成方法,其特征是按如下步骤进行: (i)放大采样间距:对采样间距为p、分辨率为M1×N1的图像持分辨率不变,但是将采样间距扩大s倍,得到图像I1,扩大后的采样间距为p1,p1=s×p; (ii)模拟漫反射照明:向图像I1代表的物体添加一定的随机噪声,得到待衍射物体复振幅分布U1,其过程表示为U1=I1×exp(i×R1),其中,i为虚数单位,R1在区间[0, tπ]上具有均匀概率分布,t≤2; (iii)全息面衍射场分布计算:以物体复振幅分布U1作为输入,采用点源法,计算全息面上的衍射场分布U2,其过程表示为:U2=PBM(U1| λ, d),其中PBM()表示点源法的衍射场计算,λ表示波长,d表示衍射距离,全息图的分辨率为M×N,像素间距为p0; (iv)纯相位全息图生成:对得到的全息面衍射场分布U2取相位角,即得到纯相位的全息图H,其过程表示为:H=ang(U2),其中ang()表示对复数取相位角操作; 关键参数条件:在关键的步骤(i)中,扩大的采样间距p1和像素数量M1必须同时满足3个条件:第一,采样间距大于艾里斑直径,即p1≥2d×tan(1.22λ/p0/M),其中tan()表示正切函数;第二,采样间距小于人眼的分辨能力,即p1≤d1×1.22λ/De,其中d1是人眼观看再现像的距离,De表示人眼的瞳孔直径;第三,空间光调制器的空间带宽积限制,即p1×M1≤d×tan[arcsin(λ/p0)],其中arcsin表示反正弦函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910139552.3A CN109856944B (zh) | 2019-02-26 | 2019-02-26 | 一种提高再现像质量的纯相位全息图生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910139552.3A CN109856944B (zh) | 2019-02-26 | 2019-02-26 | 一种提高再现像质量的纯相位全息图生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109856944A CN109856944A (zh) | 2019-06-07 |
CN109856944B true CN109856944B (zh) | 2020-03-06 |
Family
ID=66898917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910139552.3A Active CN109856944B (zh) | 2019-02-26 | 2019-02-26 | 一种提高再现像质量的纯相位全息图生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109856944B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095969B (zh) * | 2018-10-23 | 2021-03-16 | 郑州轻工业学院 | 一种基于单个slm抑制计算全息再现像斑点噪声的方法 |
CN112180707B (zh) * | 2020-09-28 | 2021-11-02 | 四川大学 | 基于球面自衍射模型的球面纯相位全息图生成方法 |
CN112269308B (zh) * | 2020-10-16 | 2021-08-10 | 北京航空航天大学 | 一种基于纯相位复合曲面全息图的全息显示*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102089720A (zh) * | 2008-07-09 | 2011-06-08 | 大日本印刷株式会社 | 全息图作成方法、全息图及使用该全息图的安全介质 |
CN104952034A (zh) * | 2014-03-27 | 2015-09-30 | 香港城市大学 | 复杂全息图向相位全息图的转换 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995455B (zh) * | 2014-05-16 | 2017-01-04 | 北京理工大学 | 在有限相位变化范围的抑制3d全息显示散斑噪声的方法 |
CN104376532B (zh) * | 2014-05-19 | 2018-11-09 | 中国人民解放军装甲兵工程学院 | 一种降低n-lut法计算全息再现像相干噪声的方法 |
-
2019
- 2019-02-26 CN CN201910139552.3A patent/CN109856944B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102089720A (zh) * | 2008-07-09 | 2011-06-08 | 大日本印刷株式会社 | 全息图作成方法、全息图及使用该全息图的安全介质 |
CN104952034A (zh) * | 2014-03-27 | 2015-09-30 | 香港城市大学 | 复杂全息图向相位全息图的转换 |
Non-Patent Citations (4)
Title |
---|
Enhancement on the generation of sampled phase-only holograms;P.W.M. Tsang等;《CHINESE OPTICS LETTERS》;20150610;第13卷(第6期);060901-1至060901-4 * |
Generation of phase-only Fresnel hologram based on down-sampling;P.W.M. Tsang等;《OPTICS EXPRESS》;20141020;第22卷(第21期);25208-25214 * |
Non-iterative Generation of Phase-only Holograms based on Dynamic Down-Sampling;P.W.M. Tsang等;《Digital Holography and 3-D Imaging 2017》;20170601;全文 * |
基于点源模型计算全息图快速生成算法的研究进展;金晓宇等;《激光与光电子学进展》;20181010;第55卷(第10期);48-57 * |
Also Published As
Publication number | Publication date |
---|---|
CN109856944A (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8363295B2 (en) | Computer generated hologram | |
JP4932490B2 (ja) | ホログラフィック装置および方法 | |
JP5502113B2 (ja) | コンピュータ生成ビデオホログラムを計算する方法及びデバイス | |
JP6249977B2 (ja) | 映像ホログラムおよび映像ホログラム再生装置 | |
CN109856944B (zh) | 一种提高再现像质量的纯相位全息图生成方法 | |
Blinder et al. | The state-of-the-art in computer generated holography for 3D display | |
US7768684B2 (en) | 3D display | |
CN108885425B (zh) | 用于编码重建三维物体的复值信号的方法和装置 | |
CN111240177B (zh) | 一种基于分层像素扫描算法的全息散斑噪声抑制方法 | |
Shimobaba et al. | Real-time and low speckle holographic projection | |
Zeng et al. | Off-axis phase-only holograms of 3D objects using accelerated point-based Fresnel diffraction algorithm | |
Chang et al. | Fast calculation of computer generated hologram based on single Fourier transform for holographic three-dimensional display | |
CN110647023B (zh) | 一种基于部分角谱法的快速全息图生成和高质量再现方法 | |
Li et al. | Large-size holographic display method based on effective utilization of two spatial light modulators | |
Chang et al. | 36.2: Speckle Suppression in a Scaled Holographic Display from Single Phase‐Only Computer Generated Hologram | |
CN115291490A (zh) | 基于光学计算的拓展距离全息显示方法 | |
KR20200111585A (ko) | 홀로그램 전파 회절무늬의 업 샘플링을 이용한 광시야각 홀로그래픽 디스플레이 장치 |
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 |