CN107991854A - 基于层析法的立体仿真全息物体再现方法 - Google Patents
基于层析法的立体仿真全息物体再现方法 Download PDFInfo
- Publication number
- CN107991854A CN107991854A CN201610944698.1A CN201610944698A CN107991854A CN 107991854 A CN107991854 A CN 107991854A CN 201610944698 A CN201610944698 A CN 201610944698A CN 107991854 A CN107991854 A CN 107991854A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- plane
- hologram
- msup
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000004587 chromatography analysis Methods 0.000 title claims abstract description 13
- 238000001093 holography Methods 0.000 title claims abstract description 10
- 239000007787 solid Substances 0.000 title claims abstract description 10
- 230000003287 optical effect Effects 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 4
- 229910052934 alunite Inorganic materials 0.000 description 2
- 239000010424 alunite Substances 0.000 description 2
- 239000003205 fragrance Substances 0.000 description 2
- YNPNZTXNASCQKK-UHFFFAOYSA-N phenanthrene Chemical compound C1=CC=C2C3=CC=CC=C3C=CC2=C1 YNPNZTXNASCQKK-UHFFFAOYSA-N 0.000 description 2
- KPZTWMNLAFDTGF-UHFFFAOYSA-D trialuminum;potassium;hexahydroxide;disulfate Chemical compound [OH-].[OH-].[OH-].[OH-].[OH-].[OH-].[Al+3].[Al+3].[Al+3].[K+].[O-]S([O-])(=O)=O.[O-]S([O-])(=O)=O KPZTWMNLAFDTGF-UHFFFAOYSA-D 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation 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
- 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/16—Processes or apparatus for producing holograms using Fourier transform
-
- 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)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了基于层析法的立体仿真全息物体再现方法,涉及全息光场技术领域,包括以下步骤:将处于(x,y,z)坐标系中的3D物体沿Z轴方向分层,设物函数为G(x,y,z);获得截平面后,再逐个计算每个截面在全息平面上的复振幅分布;计算完所有的截面Gn(x,y)的复振幅分布Pn(x,y),对其进行叠加运算,就可得到全息平面的总光场P(x,y);对P(x,y)进行编码就可获得三维菲涅尔计算全息图,再现时,将三维菲涅尔计算全息图放在光路中,就可观察到具有深度分布的3D仿真物体。本发明的立体仿真全息物体再现方法,可获得的计算全息图易通过光学方法进行再现而获得较好的3D效果。
Description
技术领域
本发明涉及全息光场技术领域,尤其是基于层析法的立体仿真全息物体再现方法。
背景技术
计算全息处理2D图像技术已经很成熟,应用也十分广泛。然而在实际生活中的物体大多数是3D的,3D物体的数据量比2D图形大很多,直接对3D物体进行编码计算,计算量过于庞大,要解决此问题,一是需要提高计算机硬件运算速度,二则需要对编码方法进行改进。
发明内容
本发明提出的一种基于层析法的立体仿真全息物体再现方法,获得的计算全息图易通过光学方法进行再现而获得较好的3D效果。
本发明的技术方案是这样实现的:
基于层析法的立体仿真全息物体再现方法,包括以下步骤:
步骤1:将处于(x,y,z)坐标系中的3D物体沿Z轴方向分层,设物函数为G(x,y,z),在沿Z方向分层后则可表示成:
式中Gn(x,y)为一系列与(x,y)平面平行的平面,每个平面之间间距为△z,n是沿z方向的抽样平面的个数;
步骤2:获得截平面后,再逐个计算每个截面在全息平面上的复振幅分布,以第n个截面Gn(x,y)为例,在全息平面上的复振幅分布P(x,y)是Gn(x,y)的菲涅尔衍射:
步骤3:计算完所有的截面Gn(x,y)的复振幅分布Pn(x,y),对其进行叠加运算,就可得到全息平面的总光场P(x,y);
步骤4:对P(x,y)进行编码就可获得三维菲涅尔计算全息图,再现时,将三维菲涅尔计算全息图放在光路中,就可观察到具有深度分布的3D仿真物体。
本发明通过提供的基于层析法的立体仿真全息物体再现方法,其有益效果在于:对3D物体的复杂程度不敏感;对各层面菲浬耳衍射计算可引入快速傅立叶变换提高算法的效率;该方法获得的计算全息图易通过光学方法进行再现而获得较好的3D效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为层析法制作菲涅尔全息图;
图2为模拟分层示意图;
图3为E在40cm处,S在60cm处的计算全息图;
图4为E在30cm处,S在40cm处再现像图;
图5为E在40cm处,S在60cm处再现像图;
图6为E在40cm处,S在80cm处再现像图;
图7为E在40cm处,S在1m处再现像图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
层析法是将3D物体沿深度方向分层,将各层面的二维截面图分别在全息面成像的复振幅分布相加,以此作为原三维物体的复振幅分布,再加入参考光,通过菲涅尔近似法计算即可制成一张全物场的全息图。
本实施例的基于层析法的立体仿真全息物体再现方法,包括以下步骤:
步骤1:如图1所示,将处于(x,y,z)坐标系中的3D物体沿Z轴方向分层,设物函数为G(x,y,z),在沿Z方向分层后则可表示成:
式中Gn(x,y)为一系列与(x,y)平面平行的平面,每个平面之间间距为△z,n是沿z方向的抽样平面的个数;
步骤2:获得截平面后,再逐个计算每个截面在全息平面上的复振幅分布,以第n个截面Gn(x,y)为例,在全息平面上的复振幅分布P(x,y)是Gn(x,y)的菲涅尔衍射:
式中zn=z0-n△z,是第n个截面Gn(x,y)到全息平面的距离。要求得Pn(x,y),需要对Gn(x,y)做以下几步计算:
(1)给Gn(x,y)乘上一个二次相位因子
(2)对乘积进行傅里叶变换;
(3)对计算结果再乘以一个相位因子
步骤3:计算完所有的截面Gn(x,y)的复振幅分布P n(x,y),对其进行叠加运算,就可得到全息平面的总光场P(x,y)。
步骤4:对P(x,y)进行编码就可获得三维菲涅尔计算全息图,再现时,将三维菲涅尔计算全息图放在光路中,就可观察到具有深度分布的3D仿真物体。
层析法中分层的多少决定了最终再现图的纵向分辨率,分层越多,分辨率越高。同时,层数的增加也会带来计算量的增大,因此,我们需要在满足分辨率要求的同时,尽量减少分层次数。
下面利用计算机模拟,验证上述方法。如图2所示,将一个物体分为S和E两层,首先,设定焦距为40cm,S层处在焦面上(即距全息面40cm),E层和全息平面距离为30cm,将两层传播到全息面的复振幅进行叠加后,加入参考光生成一幅计算全息图,并在距离40cm处再现,其计算全息图如图3所示,再现像如图4所示;接着再给出一组E层位于焦面上,S平面距全息面距离分别为60cm,80cm,1m的再现像,分别如图5,6,7所示。
从图4-7可以看出,当S层在焦面上,E层离焦时,E层再现像模糊不清,S层再现像很清晰;而当E层在焦面上,S层离焦时,则E层再现像清晰,S层再现像模糊,且3幅再现像中,E层再现像清晰度保持不变,S层偏离焦距平面距离越大,再现像清晰度越差。
通过上述两个层面的再现效果,说明模拟实验中制作的计算全息图满足了三维显示的要求,可以同时再现不同深度的图像。因此,通过层析法,可以进一步实现真实3D物体的显示技术。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.基于层析法的立体仿真全息物体再现方法,其特征在于,包括以下步骤:
步骤1:将处于(x,y,z)坐标系中的3D物体沿Z轴方向分层,设物函数为G(x,y,z),在沿Z方向分层后则可表示成:
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msub>
<mi>G</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
式中Gn(x,y)为一系列与(x,y)平面平行的平面,每个平面之间间距为△z,n是沿z方向的抽样平面的个数;
步骤2:获得截平面后,再逐个计算每个截面在全息平面上的复振幅分布,以第n个截面Gn(x,y)为例,在全息平面上的复振幅分布P(x,y)是Gn(x,y)的菲涅尔衍射:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>jkz</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>j&lambda;z</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mi>j</mi>
<mfrac>
<mi>k</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>z</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<munder>
<mrow>
<mo>&Integral;</mo>
<mo>&Integral;</mo>
</mrow>
<mi>&infin;</mi>
</munder>
<mo>&lsqb;</mo>
<msub>
<mi>U</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mi>j</mi>
<mfrac>
<mi>k</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>z</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>j</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mrow>
<msub>
<mi>&lambda;z</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mi>&xi;</mi>
<mo>+</mo>
<mi>y</mi>
<mi>&eta;</mi>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
<msub>
<mi>d</mi>
<mi>&xi;</mi>
</msub>
<msub>
<mi>d</mi>
<mi>&eta;</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤3:计算完所有的截面Gn(x,y)的复振幅分布Pn(x,y),对其进行叠加运算,就可得到全息平面的总光场P(x,y);
步骤4:对P(x,y)进行编码就可获得三维菲涅尔计算全息图,再现时,将三维菲涅尔计算全息图放在光路中,就可观察到具有深度分布的3D仿真物体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610944698.1A CN107991854A (zh) | 2016-10-26 | 2016-10-26 | 基于层析法的立体仿真全息物体再现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610944698.1A CN107991854A (zh) | 2016-10-26 | 2016-10-26 | 基于层析法的立体仿真全息物体再现方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107991854A true CN107991854A (zh) | 2018-05-04 |
Family
ID=62029041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610944698.1A Pending CN107991854A (zh) | 2016-10-26 | 2016-10-26 | 基于层析法的立体仿真全息物体再现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107991854A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109445262A (zh) * | 2018-12-28 | 2019-03-08 | 天津利信科技有限公司 | 全息数据的存储方法及装置 |
CN111443583A (zh) * | 2020-04-14 | 2020-07-24 | 北京航空航天大学 | 一种基于全息图优化分割计算的快速全息图计算方法 |
CN114895542A (zh) * | 2022-04-27 | 2022-08-12 | 安徽大学 | 一种三维计算全息图非迭代快速生成方法 |
-
2016
- 2016-10-26 CN CN201610944698.1A patent/CN107991854A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109445262A (zh) * | 2018-12-28 | 2019-03-08 | 天津利信科技有限公司 | 全息数据的存储方法及装置 |
CN111443583A (zh) * | 2020-04-14 | 2020-07-24 | 北京航空航天大学 | 一种基于全息图优化分割计算的快速全息图计算方法 |
CN114895542A (zh) * | 2022-04-27 | 2022-08-12 | 安徽大学 | 一种三维计算全息图非迭代快速生成方法 |
CN114895542B (zh) * | 2022-04-27 | 2023-11-21 | 安徽大学 | 一种三维计算全息图非迭代快速生成方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107958475A (zh) | 基于深度学习生成网络的变角度光照层析方法及装置 | |
CN104021587B (zh) | 基于计算全息技术的大场景真三维显示快速生成方法 | |
CN110458957A (zh) | 一种基于神经网络的图像三维模型构建方法及装置 | |
CN104133241B (zh) | 波场分离方法和装置 | |
CN103955127A (zh) | 一种相位调制全视差全息体视图实现方法 | |
CN107991854A (zh) | 基于层析法的立体仿真全息物体再现方法 | |
CN103745497B (zh) | 植物生长建模方法和*** | |
CN102054296A (zh) | 一种局部刚性网格变形方法 | |
CN106384386B (zh) | Lod模型生成中网格处理方法、***以及3d重建方法和*** | |
CN103763543A (zh) | 合成全息图的采集方法 | |
CN103941567A (zh) | 三维数据虚拟分层全息图方法 | |
CN103941569A (zh) | Led多角度全息三维数据再现方法 | |
CN101458822B (zh) | 一种3d模型的计算全息图快速生成方法 | |
KR101573362B1 (ko) | 고속의 신축적인 디지털 홀로그램 생성 방법 및 장치 | |
KR101021127B1 (ko) | 룩업 테이블과 영상의 공간적 중복성을 이용한 컴퓨터 형성 홀로그램 산출 방법 및 그 장치 | |
CN109358481A (zh) | 一种复用体全息透镜再现方法 | |
CN102183761B (zh) | 星载干涉合成孔径雷达数字高程模型重建方法 | |
CN103217888A (zh) | 一种计算机制合成彩色彩虹全息图的制作方法 | |
CN108230378A (zh) | 一种基于光线追踪的计算全息遮挡处理算法 | |
CN107464287A (zh) | 基于多目标优化的曲面重构方法 | |
CN103325086B (zh) | 一种基于四面体坐标系的三维图形的变形方法 | |
Min et al. | Wigner inverse transform based computer generated hologram for large object at far field from its perspective light field | |
CN104299265B (zh) | 一种流体环境影响下的群体行为控制方法 | |
CN110426688A (zh) | 一种基于地形背景目标的sar回波模拟方法 | |
US20110310448A1 (en) | Method for calculating computer generated hologram using look-up table and apparatus thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180504 |
|
WD01 | Invention patent application deemed withdrawn after publication |