傅里叶变换光学.docx
《傅里叶变换光学.docx》由会员分享,可在线阅读,更多相关《傅里叶变换光学.docx(24页珍藏版)》请在冰豆网上搜索。
![傅里叶变换光学.docx](https://file1.bdocx.com/fileroot1/2022-10/9/27379361-8921-4a9b-b62c-dc2b6cfdfc80/27379361-8921-4a9b-b62c-dc2b6cfdfc801.gif)
傅里叶变换光学
中山大学光信息专业实验报告:
傅里叶光学变换系统
实验人:
何杰勇(11343022)合作人:
徐艺灵组号B13
一、实验目的和内容
1、了解透镜对入射波前的相位调制原理。
2、加深对透镜复振幅、传递函数、透过率等参量的物理意义的认识。
3、观察透镜的傅氏变换(FT)图像,观察4f系统的反傅氏变换(IFT)图像,并进行比较。
4、在4f系统的变换平面(T)插入各种空间滤波器,观察各种试件相应的频谱处理图像。
二、实验原理
1、透镜的FT性质及常用函数与图形的关学频谱分析
透镜由于本身厚度的不同,使得入射光在通过透镜时,各处走过的光程差不同,即所受时间延迟不同,因而具有相位调制能力。
图1为简化分析,假设任意点入射光线在透镜中的传播距离等于改点沿光轴方向透镜的厚度,并忽略光强损失,即通过透镜的光波振幅分布不变,仅产生位相的变化,且其大小正比于透镜在该点的厚度。
设原复振幅分布为
的光通过透镜后,其复振幅分布受到透镜的位相调制,附加了一个位相因子
后变为
:
图1
(1)
若对于任意一点(x,y)透镜的厚度为
,透镜的中心厚度为
。
光线由该点通过透镜时在透镜中的距离为
,空气空的距离为
-
,透镜折射率为n,则该点的总的位相差为:
(2)
(2)中的k=2π/λ,为入射光波波数。
用位相延迟因子
来表示即为:
(3)
由此可见只要知道透镜的厚度函数
就可得出其相位调制。
在球面镜傍轴区域,用抛物面近似球面,可以得到球面透镜的厚度函数为:
(4)
其中
、
是构成透镜的两个球面的曲率半径。
公式(4)对双凹、双凸、或凹凸透镜都成立。
引入焦距f,其定义为:
(5)
代入(3)得:
(6)
式(6)即是透镜位相调制的表达式,它表明复振幅
通过透镜时,透镜各点都发生位相延迟。
从式(6)容易看出第一项位相因子
仅表示入射光波的常量位相延迟,不影响位相的空间分布,即波面形状,所以在运算过程中可以略去。
第二项
是具有调制作用的因子,它表明光波通过透镜的位相延迟与该点到透镜中心的距离的平方成正比。
而且与透镜的焦距有关。
当考虑透镜孔径后,有:
(7)
其中的
为透镜的光瞳函数,表达式为:
(8)
2、透镜的傅里叶变换性质
在单色平面波垂直照射下,夫琅和斐衍射光场的复振幅分布正比于衍射屏透射系数的傅里叶变换。
衍射图像的强度分布正比于衍射屏的功率谱分布。
一般情况下,我们是将夫朗和斐衍射图像成像到透镜的像方焦平面出,这就是说,作为成像元件的透镜,就相当于傅里叶变换器。
如图2所示,设单位振幅的单色平面光垂直照射一透射系数为
的衍射屏,与衍射屏相距Z处放置一焦距为f的薄透镜L,先观察其像方平面L的光场分布。
为了讨论方便,这里我们忽略透镜材料的吸收、散射、透镜表面的反射以及透镜孔径大小等因素的影响。
图2透镜的傅里叶变换性质
设
、
、
、
分别表示衍射屏后、透镜输入平面、输出平面以及像方平面出光波场的复振幅分布。
由于透镜的相位调制特性,输出平面与输入平面出光波场之间的关系由下式决定:
(9)
而从透镜输出平面到像方焦平面,光波相当于经历一次菲涅耳衍射。
夫朗和斐近似下观察到平面上的衍射光场复振幅:
=
(10)
式中u和v分别表示
和
方向的空间频率。
于是由(9)和(10)式,透镜像方焦平面上的光波场复振幅
分布应具有如下形式:
=
(
)(11)
在单位振幅的平面波垂直照射下,透镜衍射屏的光波场复振幅分布
即等于衍射屏的透射系数
,故其频谱分布为:
(12)
该频谱分量从衍射屏传播到透镜的输入平面处,产生一个相位延迟
,即有:
(13)
在傍轴条件下
具有如下的形式:
(14)
由此可以得到透镜输入平面处光波场的频谱分布为:
(15)
代入(11)得透镜像方焦平面处的广场分布为:
=
(
)(16)
从上式可以看到,在单色平面波垂直照射下,透镜像方焦平面处的光场除了一个常数因子外和一个二次因子外,其余的反应了衍射屏透射系数得傅里叶变换。
经过进一步的分析我们可以得到在用透镜对二维关学图像进行傅里叶变换时,若将图像放置在透镜的物方焦平面上,则在透镜的像方焦平面上得到输入图像准确的傅里叶变换。
若将输入图像放置在透镜与其像方焦平面之间,则像方焦平面上频谱图样的大小可随衍射屏到像方焦平面的距离的变化而改变;并且当输入图像紧贴透镜后放置时可获得最大的频谱图样。
而对于球面波照射时,傅里叶变换平面将不是在透镜的像方平面。
而是光源的共轭像平面上。
3.透镜孔径的衍射与滤波特性
由于孔径的衍射效应,任何具有有限大小通过光孔径的光学成像系统,均不存在如几何光学中所说的理想像点。
所谓共轭像点,实际上是由系统孔径引起的,以物点的几何像点为中心的夫琅和斐衍射图样的中央亮斑——艾里斑。
其次,透镜有限大小的通光孔径,也限制了衍射屏函数的较高频率成分(具有较大入射倾角的平面波分量)的传播。
这可以从图3可以看出:
图3:
透镜孔径引起渐晕效应
透过衍射屏的基频平面波分量1可以全部通过透镜,具有较高(空间)频率的平面波分量2只能部分通过,而高频平面波分量3则完全不能通过。
这样,在透镜像方焦平面上的光波场中就缺少了衍射屏透射光场中部分高频成分,因此,所得衍射屏函数的频谱将不完整。
这种现象称为衍射的渐晕效应。
由此可将,从光信息处理角度来讲,透镜孔径的有限大小,使得系统存在着有限大小的通频宽带和截止频率;从光学成像的角度来讲,则使得系统存在着一个分辨极限。
4.相干光学图像处理系统(4f系统)
用夫琅和斐衍射来实现图像的频谱分解,最重要的意义是为空间滤波创造了条件,由于衍射场就是屏函数的傅里叶频谱面,空间频率(u,v)与衍射场点位置(
)一一对应,使得人们可见从改变频谱入手来改造图像,进行信息处理。
为此,设计了图4所示的图像处理系统。
图44f图像处理系统
在此系统中,两个透镜
、
成共焦组合,
的前焦面(x,y)为物平面O,图像由此输入,
的后焦面
为像平面I,图像在此输出。
共焦平面(
)称为变换平面T,在此可以安插各种结构和性能的屏(即空间滤波器)。
当平行光照射在物平面上时,整个OTI系统成为相干成像系统。
由于变换平面上空间滤波器的作用,使输出图像得以改造,所以OTI系统又是一个相干光学信息处理系统。
这里先研究它的成像问题。
我们将相干光学系统的成像过程看作两步:
第一步,从O面到T面,使第一次夫琅和斐衍射,它起分频作用。
第二步,从T面到I面,再次夫琅和斐衍射,起合成作用,即综合频谱输出图像。
在这样的两步中,变换平面T处于关键地位,若在此处设置光学滤波器,就能起到选频作用。
要想作到图像的严格复原,T面必须完全畅通无阻。
此处的4f系统每次衍射都是从焦面到焦面,这就保证了复振幅的变换是纯粹的傅里叶变换。
如果光波能够自由通过变换平面,即连续两次的傅里叶变换,函数的形式基本复原,只是自变量变号,
即图像倒置。
在有源滤波器的情况下,
这里为滤波器的透过率函数,这也是我们进行滤波实验的依据。
5.空间滤波实验
要从输入图像中提取或排除某种信息,就要事先研究这类信息的频谱特征,然后针对它制备相应的空间滤波器置于变换平面,经过第二次衍射合成后,就可以达到预期的效果,光信息处理的原理也就是基于如此。
三、实验仪器与装置图
实验仪器:
激光器、准直系统、傅里叶透镜、傅里叶变换试件、频谱处理器、CCD光电接收器;
实验装置图:
如图5
图5实验装置图
四、实验内容
1.根据傅里叶变换光路装置简图摆好光路,打开激光电源,调整光路。
2.开启电脑,运行csylaser软件。
调节光路中各器件的位置,以得到样品较为清晰的傅里叶变换图像(根据所用样品,最终应得到“米”字图像)。
并将图像保存,作为原始数据。
3.根据反傅里叶变换光路装置简图(4f系统)摆好光路,调节器件位置,以得到样品最为尖锐的反傅里叶变换图像,并保存。
在调节时,主要是调节CCD的位置,傅里叶透镜的位置摆放好不要轻易乱动。
4.在频谱处理器的位置加上带有狭缝的滤波片,将激光依次透过狭缝,观察不同的狭缝对于光波的透过作用的不同,保存图像,并分析。
5.关闭激光器和电脑电源,整理好仪器。
实验结束。
四、实验数据记录与分析
1.观察样品的傅里叶频谱图图样。
图6所示为样品的原图样,图7为其频谱图:
图6样品示意图
图7样品傅里叶变换频谱图
由图可知,样品经过傅里叶变换得到的频谱图
理论验证:
用MatLab程序编辑一个二维矩阵做出一个图6所示的图像,使其发光部分值为1,不发光部分为0。
如图8(a)所示。
图像在焦平面上的频谱图为图像经过了一次二维傅里叶变换,再将频谱搬移到中心,得到的频谱图如图8(b),可以看到,理论和实际得到的图像很相似,都为“米”字型。
图8MatLab的模拟图(a)原始图像;(b)图像的频谱图
对比分析:
结合图7与图8(b)可以看出,所测样品的傅里叶变换图像,类似一个“米”字。
分析可知,中心的十字经变换后仍为十字形,而顶角处的三角形则经过傅里叶变换后变为“×”形,因为最终样品的傅里叶变换图像为“米”字形。
图7中的“米”字不是很清晰,分析其原因,主要是由于激光器所发激光的光强太大,导致CCD过曝光,使中心十字过亮而“×”形不明显。
图质稍有模糊,分析原因,除了CCD的采样分辨率太小之外,可能是由于光路没有调节至完全共轴,或者CCD没有恰好在透镜焦点上,导致了实际上没有在焦平面上获取图像。
由数学理论分析可知,频谱可以认为表示的是图像的衬比度的变化程度。
图中处于中心位置的是零频位置,也就是图像的直流成分,可以理解为光强没有发生变化的表象;越远离中心的频谱自然指的是衬比度变化的部分,越远离中心,衬比度的变化程度越大。
又由阿贝成像原理,我们可以知道物体在焦面上成的像,其实是图像在透镜经过夫琅禾费衍射所形成的像,于是水平部分指的应是竖直方向上的衬比度变化,竖直部分指的是水平方向上的衬比度变化。
而图中的“X”部分就是箭头的斜边部分的衬比度了。
故我们实验中得到的频谱图与理论得到的很符合。
改进方法:
降低激光光强,是CCD尽量不产生过曝;
调整光路,使CCD尽量正对傅里叶透镜的焦面,得到尽量更好的图样;
更换焦距较大的透镜,由于焦距太小,导致对光的折射角度过大,使对焦点的确定变得很难。
故可以试着更换焦距较大的透镜。
2.观察样品的反傅里叶频谱图图样。
摆好并调节好光路后,可得样品的反傅里叶变换图像如下图所示:
图9样品的反傅立叶变换图像
由上图可以看出,待测样品的反傅里叶变换图像就是它本身,但是图像发生了反转。
调节光路,当获得最为清晰,且边缘最为尖锐的图像时,为较理想的图像。
理论验证:
用MatLab软件编辑程序将由步骤1得到的图像的频谱图进行一次傅里叶变换,得打如图10的图样。
可知,源图像的频谱图经过傅里叶变换后得到的图样形状和原图一致,但是发生了翻转,与实验得到的结果相同。
图10基于MatLab的频谱傅里叶变换图样
对比分析:
实验所得图像和理论图相比,只有部分图像,原因是CCD所采图像不全,其根本的原因是没有调好4F系统导致在CCD上得到的是放大后的像,于是只能得到部分图像。
实验得到的图像边缘模糊,不光滑,其可能原因是CCD没有刚好落在透镜焦点位置,也可能是因为两傅里叶透镜没有完全共轴,或者各器件之间的距离与4F系统所要求的距离有偏差。
实验得到的图像的亮度不完全均匀,原因可能是激光的光点没有全部照在CCD上,使亮度不均匀。
改进方法:
保证光路准直后,检查射