全息图的数字化频域滤波及数值再现研究

更新时间:2023-07-24 14:21:01 阅读量: 实用文档 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

【摘要】采用计算机对普通离轴计算全息图及博奇型修正离轴参考光计算全息图进行数字化滤波操作,可在频域将零级及孪生像消除,从而得到单一的清晰实像或虚像。利用快速傅里叶变换算法计算菲涅耳衍射积分,实现了计算全息图的设计制作、频域滤波、再现过程的全数字化,计算机数

第35卷 第6期 电 子 科 技 大 学 学 报 Vol.35 No.6 2006年12月 Journal of University of Electronic Science and Technology of China Dec. 2006 ·光电子学工程与应用·

全息图的数字化频域滤波及数值再现研究

仇 宇

(绵阳师范学院学生工作部 四川 绵阳 621000 )

【摘要】采用计算机对普通离轴计算全息图及博奇型修正离轴参考光计算全息图进行数字化滤波操作,可在频域将零级及孪生像消除,从而得到单一的清晰实像或虚像。利用快速傅里叶变换算法计算菲涅耳衍射积分,实现了计算全息图的设计制作、频域滤波、再现过程的全数字化,计算机数值模拟结果表明该方法具有再现图像信噪比高、操作过程简便、计算速度快、灵活性强等特点。

关 键 词 计算机制全息图; 傅里叶变换; 数字滤波; 数值再现

中图分类号 O438.1 文献标识码 A

Digital Frequency Filtering and Numerical Reconstruction of Computer Generated Hologram QIU Yu (Students’ Administration Department, Mianyang Normal University Mianyang Sichuan 621000)

Abstract Digital filtering method for frequency domain is employed to process conventional and modified off-axis reference beam computer generated holograms, then the spatial spectra of zero order image and one twin image can be eliminated, and single clear virtual image or real image can be reconstructed numerically. In this paper, fast Fourier transform algorithm is used to calculate the Fresnel diffraction integration, and the total digitization of design, fabrication, frequency filtering, and reconstruction procedures are achieved successfully. The numeric simulation by computer shows that this method has some advantages such as high signal-to-noise ratio, convenient operation procedures, fast computing speed, and high flexibility.

Key words computer generated hologram; Fourier transform; digital filtering; numerical reconstruction

传统的光学全息图制作需有真实的记录物体存在,还需要激光光源、防震平台、光学元件、记录材料等硬件条件,对记录环境条件(振动、噪音等)也有较高的要求,均给全息图的制作质量与可重复性带来较多问题。不用光学干涉记录而用人工的方法,即使用计算机与绘图仪,根据物波的数学描述,将其振幅和相位以一定的编码方式记录下来形成全息图,这就是计算全息图;采用光敏电子元件(Charge Couple Device,CCD)代替普通照相干版,并将记录的全息图输入计算机,再通过一定的算法重建物体的再现像,在计算机屏幕上显示出来,这就是数字全息术。计算全息和数字全息在技术上突破了光学全息的某些限制和不足,有利于实现光学技术与计算机技术、电子技术的有机结合,有利于改善光学技术的处理精度、可靠性、技术复杂性等,是光学全息和未来全息技术的重要发展方向[1]。

由于CCD的分辨率(通常在100 lines/mm左右)低于传统的感光记录介质大约一个数量级,限制了数字全

使得同轴及准同轴记录的数字全息图的再现像受到直透光和共轭息图记录时的参物光夹角只能在1°左右[2],

像的严重干扰,特别是零级像,占据了大部分能量而在图像的当中形成一个大又亮的亮斑,使得重建图像的质量不高。目前,消除数字全息图的零级像及共轭像的方法主要有两种:(1) 采用相移技术来改变参考光束的相位并记录多幅相应的全息图,再通过计算机数字相减技术,能很好地去除直透光和共轭像的影响,获得高质量的数字图像[3];(2) 利用离轴数字全息图技术对全息图进行傅里叶变换,通过在频谱面上将彼此分离的零级像及共轭像的频谱滤掉,可再现得到所需的单一图像[4]。本文将离轴数字全息图的数字滤波技术收稿日期:2006 05 29

作者简介:仇 宇(1965 – ),男,本科,讲师,主要从事网络技术方面的研究.

【摘要】采用计算机对普通离轴计算全息图及博奇型修正离轴参考光计算全息图进行数字化滤波操作,可在频域将零级及孪生像消除,从而得到单一的清晰实像或虚像。利用快速傅里叶变换算法计算菲涅耳衍射积分,实现了计算全息图的设计制作、频域滤波、再现过程的全数字化,计算机数

第6期 仇 宇 : 全息图的数字化频域滤波及数值再现研究 971 应用于计算全息图中,可实现全息图制作、存储、处理和再现全过程的数字化。

1 修正离轴参考光计算全息图的编码及制作

普通菲涅耳离轴光学全息图的记录光路如图1所示,全

息图记录平面xy位于物平面xiyi的菲涅耳衍射区内,假设

物体的振幅透过率函数为U(xi, yi),则传播距离z到达全息图

平面的物光波为U(xi, yi)的菲涅耳衍射光场[5]:

O(x,y)=1

iλzeikz∫∫ +∞

∞U(xi,yi)eikx xi)2+(y yi)2]2zdxidyi (1)

将式(1)的指数因子展开并进行整理可得:

O(x,y)=1

iλzeeikzik22(x+y)2z 图1 菲涅耳离轴光学全息图的记录光路 ∫∫ +∞

∞U(xi,yi)eik2xi+yi2) i2π(fx+fy)xiyi2zedxidyi=1

iλzeeikzik22(x+y)2zik2(xi+yi2) 2z U(xi,yi)e (2)

式中 fx=x/λz,fy=y/λz,显然菲涅耳衍射积分可看成是乘积U(xi,yi)eik2xi+yi2)2z的傅里叶变换,采用二维快速傅里叶变换算法(FFT)可简化式(1)的计算,提高了计算速度。

设O(x,y)=A(x,y)ei (x,y),平面参考光波R(x,y)=Rei2παx,α为平面参考波的载频,则在全息图平面上的干涉光强分布为:

I(x,y)=O(x,y)+R(x,y)=O(x,y)+R2+R*(x,y)O(x,y)+O*(x,y)R=

A(x,y)+R+2RA(x,y)cos[2παx (x,y)]

222222 (3) 式中 O(x,y)+R2=A(x,y)+R2表示零级像,再现时表现为沿再现光波传播方向的零级衍射光晕(干扰噪

声),作用是产生偏置分量,使I(x,y)成为非负实函数;R*(x,y)O(x,y)包含原始物光波的复振幅信息,再现时生成物体的虚像;O*(x,y)R包含物光波的共轭波函数,再现时生成物体的共轭像(通常为实像)。式(3)为采用光学照相方法记录的全息图的透过率函数,其中,O(x,y)+R2=A(x,y)+R2在再现时为产生多 余的衍射像,增加了全息图记录时的带宽要求。根据抽样定理,计算全息图的抽样点数正比于全息图函数

要求的抽样点数较多,制作时很不方便。为此博奇(Burch)的总带宽,因此,直接采用式(3)制作计算全息图时,

利用计算机设计的灵活性,提出了博奇型修正离轴参考光计算全息图的编码方法[6]。该方法通过采用常数 偏置项来代替式(3)中的A(x,y)+R2,重新构造新的全息图函数为: 222

I(x,y)=0.5{1+A(x,y)cos[2παx (x,y)]}

(4)

由此可将全息图的总抽样点数降低为原先式(3)所要求的25%,全息图再现时的噪声分量只剩下中间的平面波分量,达到了降低带宽、减少抽样数、提高再现图像质量等诸多优点。

基于以上的分析,本文以字母A作为待记录的物体,如图2所示,图像分辨率为256×256像素。按照式(3)、

(4)的光强分布,采用计算机模拟技术,设计了离轴计算全息图及博奇型修正离轴参考光计算全息图,物平面与全息图平面的距离z选为190 mm,其结果如图3、4所示。

图2 记录物体的图像 图3 普通离轴计算全息图 图4 博奇型修正离轴参考光计算全息图

【摘要】采用计算机对普通离轴计算全息图及博奇型修正离轴参考光计算全息图进行数字化滤波操作,可在频域将零级及孪生像消除,从而得到单一的清晰实像或虚像。利用快速傅里叶变换算法计算菲涅耳衍射积分,实现了计算全息图的设计制作、频域滤波、再现过程的全数字化,计算机数

972 电 子 科 技 大 学 学 报 第35卷 2 离轴计算全息图的数字滤波

由以上分析可知,离轴全息图的再现像由于受到离轴参考光角度的调制,其零级像、原始像及共轭像在空间是分离的,原始像及共轭像对称地分布在零级像两边。若对全息图进行傅里叶变换,为:

F(μ,ν)= {I(x,y)}=∫∫+∞

∞I(x,y)e i2π(μx+νy)dxdy (5)

则在频谱面上即可得到分离的零级衍射像、原始像、共轭像的频谱。图5、6所示为采用快速傅里叶变换算法得到的离轴计算全息图及博奇型修正离轴参考光计算全息图的频谱图,由于采用常数项偏置,博奇型修正离轴参考光计算全息图的零级频谱只剩下一个中央亮点,成功消除了普通离轴全息图零级频谱中由于物函数的自相关项引起的弥散光斑。

针对在频谱平面上彼此分离的零级像、原始像、共轭像的频谱,采用窗口滤波函数rrect(μ,ν)对全息图的频谱进行滤波操作h(x,y)=F(μ,ν)rrect(μ,ν),方便地滤掉不需要的频谱分量,保留所需共轭像或原始像的频谱。

图5 普通离轴计算全息图的频谱图 图6 博奇型修正离轴参考光计算全息图的频谱图

3 离轴计算全息图的数值再现

对上面滤波获得的单一像的频谱进行逆傅里叶变换,得到的复振幅分布h(x,y)中只包含共轭像或原始像的信息,利用计算机对h(x,y)进行数值再现,对其再次进行菲涅耳衍射积分计算,传播z0距离后到达像平面上的复振幅分布为:

1ikz0U′(X,Y)=eikz0∫∫ +∞ ∞h(x,y)R′(x,y)eik[(X x)2+(Y y)2]2z0dxdy (6)

再现距离z0、再现光波R′(x,y)的选取将决定再现出的单一像是物体的共轭像或原始像,如采用记录时参考光波的共轭波作为再现光波,在全息图平面后z0=z处再现出清晰的共轭像(通常为实像),而采用原始参考

将计算得到的再现结果直接光波作为再现光波,在全息图平面前z0= z处再现出清晰的原始像(通常为虚像)。

输出到计算机屏幕等数字显示设备上,可实现物体的数值再现。

本文以博奇型修正离轴参考光计算全息图为例,说明滤波前后再现像的效果。图7所示为未进行滤波操作的再现像,正如前面所分析的那样,零级像占据了大部分能量,在图像的当中形成一个大的亮斑,对再现的共轭像形成明显的干扰(图像上叠加了周期性的直条纹),使得重建图像的质量不高。图8、9所示为滤波操作后分别保留共轭像及原始像的频谱并进行数值再现的图像,可见零级像及孪生像均得到了彻底的消除,且数值再现后的单一图像清晰、无干扰,信噪比明显提高。

图7 博奇型修正离轴参考 图8 数字滤波后仅保留共 图9 数字滤波后仅保留原

光计算全息图的再现像 轭像频谱后的再现像 始像频谱后的再现像

(下转第976页)

本文来源:https://www.bwwdw.com/article/8inm.html

微信扫码分享

《全息图的数字化频域滤波及数值再现研究.doc》
将本文的Word文档下载到电脑,方便收藏和打印
推荐度:
点击下载文档
下载全文
范文搜索
下载文档
Top