1、 第六章第六章 空间滤波空间滤波空间滤波的目的空间滤波的目的:是通过有意识地改变像的频谱:是通过有意识地改变像的频谱,使像产生所希望的变换使像产生所希望的变换 空间滤波和光学信息处理紧密联系在一起,空间滤波是光学信息处理的空间滤波和光学信息处理紧密联系在一起,空间滤波是光学信息处理的常用手段之一。常用手段之一。光学信息处理光学信息处理是一个更为宽广的领域,它主要是用光学方法实现对输入是一个更为宽广的领域,它主要是用光学方法实现对输入信息的各种变换或处理。信息的各种变换或处理。18731873年年阿贝阿贝(AbbeAbbe)提出二次成像理论。阿贝于提出二次成像理论。阿贝于18931893年、波特
2、年、波特(Porter)(Porter)于于 19061906年为验证这一理论所做的实验,科学地说明了成像质量与系统传递年为验证这一理论所做的实验,科学地说明了成像质量与系统传递的空间频谱之间的关系。的空间频谱之间的关系。19351935年策尼克年策尼克(ZenikeZenike)提出的相衬显微镜是空提出的相衬显微镜是空间滤波技术早期最成功的应用。间滤波技术早期最成功的应用。50 50年代,艾里亚斯年代,艾里亚斯(Elias)(Elias)及其同事发表的经典论文及其同事发表的经典论文光学和通信理论光学和通信理论和和光学处理的傅里叶方法光学处理的傅里叶方法为光学信息处理提供了有力的数学工具。为光
3、学信息处理提供了有力的数学工具。6060年代由于激光的出现和全息术的重大发展,光学信息处理进入了蓬勃年代由于激光的出现和全息术的重大发展,光学信息处理进入了蓬勃发展的新时期。发展的新时期。1946 1946年杜费年杜费(DuffieuxDuffieux)把光学成像系统看作线性滤波器,成功地用傅里叶把光学成像系统看作线性滤波器,成功地用傅里叶方法分析成像过程,发表了方法分析成像过程,发表了傅里叶变换及其在光学中的应用傅里叶变换及其在光学中的应用的著名论的著名论著。著。6.1 6.1 空间滤波的理论基础空间滤波的理论基础 P1P1P2P2P3P3 阿贝研究显微镜成像问题时,提出了一种不同于几何光学
4、的新观点:他阿贝研究显微镜成像问题时,提出了一种不同于几何光学的新观点:他将物看成是不同空间频率信息的集合,相干成像过程分两步完成:第一步是将物看成是不同空间频率信息的集合,相干成像过程分两步完成:第一步是入射光场经物平面入射光场经物平面P1P1发生夫琅发生夫琅和和费衍射,在透镜后焦面费衍射,在透镜后焦面P2P2上形成一系列衍射上形成一系列衍射斑;第二步是各衍射斑作为新的次波源发出球面次波,在像面上互相叠加,斑;第二步是各衍射斑作为新的次波源发出球面次波,在像面上互相叠加,形成物体的像。将显微镜成像过程看成是上述两步成像的过程是波动光学的形成物体的像。将显微镜成像过程看成是上述两步成像的过程是
5、波动光学的观点。后来人们称其为阿贝成像理论。阿贝成像理论不仅用傅里叶变换阐述观点。后来人们称其为阿贝成像理论。阿贝成像理论不仅用傅里叶变换阐述了显微镜成像的机理,更重要的是首次引入频谱的概念启发人们用改造频谱了显微镜成像的机理,更重要的是首次引入频谱的概念启发人们用改造频谱的手段来改造信息。的手段来改造信息。阿贝阿贝-波特实验是对阿贝成像原理最好的验证和演示。这项实验的一般波特实验是对阿贝成像原理最好的验证和演示。这项实验的一般做法做法如图所示如图所示:用平行:用平行相干光束相干光束照明一张细丝网格,在成像透镜的后照明一张细丝网格,在成像透镜的后焦面上出现周期性网格的傅里叶频谱焦面上出现周期性
6、网格的傅里叶频谱,由这些傅里叶频谱分量的再组合由这些傅里叶频谱分量的再组合从而在像平面上再现网格的像。若把各种遮挡物从而在像平面上再现网格的像。若把各种遮挡物(如光圈、狭缝、小光如光圈、狭缝、小光屏屏)放在频谱平面上,就能以不同方式改变像的频谱,从而在像平面上放在频谱平面上,就能以不同方式改变像的频谱,从而在像平面上得到由改变后的频谱分量重新组合得到的对应的像。得到由改变后的频谱分量重新组合得到的对应的像。图图(a)(a)是实验装置图。图是实验装置图。图 (b)(b)是使用一条水平狭缝时透过的频谱,对应的像是使用一条水平狭缝时透过的频谱,对应的像如图如图(c)(c)所示,它只包括网格的垂直结构
7、所示,它只包括网格的垂直结构.如果将狭缝旋转如果将狭缝旋转9090,则透过的则透过的频谱和对应的像如图频谱和对应的像如图(d)(d)、图、图(e)(e)所示。若在焦面上放一个可变光圈所示。若在焦面上放一个可变光圈,开始时开始时光圈缩小,只能通过光圈缩小,只能通过轴上的轴上的傅里叶分量;然后逐渐加大光圈,就可以看傅里叶分量;然后逐渐加大光圈,就可以看到网格的像怎样由傅里叶分量一步步综合出来。如果去掉光圈换上一个到网格的像怎样由傅里叶分量一步步综合出来。如果去掉光圈换上一个小光屏挡住零级频谱,则可以看到网格像的对比度反转。这些实验以其小光屏挡住零级频谱,则可以看到网格像的对比度反转。这些实验以其简
8、单的装置十分明确地演示了阿贝成像原理,对空间滤波的作用给出了简单的装置十分明确地演示了阿贝成像原理,对空间滤波的作用给出了直观的说明,为光学信息处理的概念奠定了基础。直观的说明,为光学信息处理的概念奠定了基础。带有灰尘的网格带有灰尘的网格(灰尘的频谱特性)灰尘的频谱特性)频谱频谱小孔阵列滤波器小孔阵列滤波器输出像输出像圆孔只让云雾状光斑通过圆孔只让云雾状光斑通过输出灰尘的像输出灰尘的像用不同直径圆孔滤波得到的像用不同直径圆孔滤波得到的像 可见:孔越小,传播到像面的高频成分就越少,像的可见:孔越小,传播到像面的高频成分就越少,像的细节损失得越多。细节损失得越多。拼合的月球表面照片拼合的月球表面照
9、片拼合的月球表面照片的频谱拼合的月球表面照片的频谱在频谱中挡住水平亮线在频谱中挡住水平亮线经滤波后月球表面的像经滤波后月球表面的像(不含拼合的垂直线结构不含拼合的垂直线结构)现在我们以一维光栅为例,用傅里叶分析的手段讨论空间滤波过程,现在我们以一维光栅为例,用傅里叶分析的手段讨论空间滤波过程,以便以便更透彻更透彻地了解改变系统频谱对像结构的影响。为简明起见,采用最典地了解改变系统频谱对像结构的影响。为简明起见,采用最典型的相干滤波系统,通常称为型的相干滤波系统,通常称为4f4f系统系统,如图所示。图中,如图所示。图中L1L1是准直透镜,是准直透镜,L2L2和和L3L3为傅里叶变换透镜,焦距均为
10、为傅里叶变换透镜,焦距均为f f。P1P1、P2P2和和P3P3分别是分别是物面物面、频谱面频谱面和和像面像面,并且,并且P3P3平面采用反演坐标系。平面采用反演坐标系。6.1.2 6.1.2 空间滤波的傅里叶分析空间滤波的傅里叶分析物面物面频谱面频谱面像面像面 设光栅常数为设光栅常数为d d,缝宽为缝宽为a a,光栅沿,光栅沿x x方向的宽度为方向的宽度为L L,则它的透过率为则它的透过率为:Lxrectdxcombdaxrectxt11111)(在在P P2 2平面上的光场分布应正比于平面上的光场分布应正比于dmLdamdaLTmsincsinc)(物面物面频谱面频谱面像面像面dLdadL
11、daLdaL1csin)(csin1csin)(csin)(csin假定假定L dL d,可以避免各级频谱的重叠,可以避免各级频谱的重叠。下面讨论在频谱面上放置不同下面讨论在频谱面上放置不同的滤波器时,在输出面上像场的变化情况:的滤波器时,在输出面上像场的变化情况:(1 1)滤波器是一个适当宽度的狭缝,只允许零级谱通过,也就是说只让)滤波器是一个适当宽度的狭缝,只允许零级谱通过,也就是说只让上式上式 中第一项通过,则狭缝后的透射光场为中第一项通过,则狭缝后的透射光场为)(csin)()(LdaLHT式中式中H H()是狭缝的透过函数。在输出平面上的场分布为)是狭缝的透过函数。在输出平面上的场分
12、布为 )()()(13 HTFxg Lxrectda3空间滤波的全部过程如图下图所示。空间滤波的全部过程如图下图所示。dmLdamdaLTmsincsinc)(物体物体物体的谱物体的谱滤波函数滤波函数滤波后的谱滤波后的谱输出像输出像 一维光栅经滤波的像(透过零级)一维光栅经滤波的像(透过零级)(2 2)狭缝加宽能允许零级和正、负一级频谱通过,这时透射的频谱包括上)狭缝加宽能允许零级和正、负一级频谱通过,这时透射的频谱包括上式中的前三项,即:式中的前三项,即:dLdadLdaLdaL1csin)(csin1csin)(csin)(csin)()(HT于是输出平面上的场分布为于是输出平面上的场分布
13、为 )()()(13 HTFxg dxjLxrectdacLxrectda3332expsin dxjLxrectdac332expsin dxdacLxrectda332cossin21 此时像与物的周期相同,但由于高频信息的丢失,像的结构变成余弦振幅此时像与物的周期相同,但由于高频信息的丢失,像的结构变成余弦振幅光栅,空间滤波过程如图。光栅,空间滤波过程如图。sincsinc(a/da/d)起什么作用?)起什么作用?滤波函数滤波函数滤波后的谱滤波后的谱输出像输出像(3 3)滤波面放置双缝,只允许正、负二级谱通过,这时系统透射的频谱为)滤波面放置双缝,只允许正、负二级谱通过,这时系统透射的频
14、谱为 dcLdcLdacdaL2sin2sin)2(sin )()(HT输出平面上的场分布输出平面上的场分布 )()()(13 HTFxg dxLxrectdacda334cos)2(sin2 在这种情况下,在这种情况下,像的结构是余弦振幅光栅。,像的结构是余弦振幅光栅。如下图所示。如下图所示。强度分布如何?强度分布如何?(4 4)频谱面上放置不透光的小圆屏,挡住零级谱,而让其余频率成分通过,)频谱面上放置不透光的小圆屏,挡住零级谱,而让其余频率成分通过,这时透射频谱可表示为这时透射频谱可表示为)()(HT)(sin)(LcdaLT 像面上的光场分布正比于像面上的光场分布正比于 )()(13
15、TFxg )(sin1 LcdaLF)()(33Lxrectdaxt Lxrectdxcombdaxrect3331)(3Lxrectda 当当a=d/2a=d/2,即缝宽等于缝的间隙时,直流分量为,即缝宽等于缝的间隙时,直流分量为1/21/2,像场的复振幅分布仍,像场的复振幅分布仍为光栅结构,并且周期与物相同,但强度分布是均匀的,即实际上看不见为光栅结构,并且周期与物相同,但强度分布是均匀的,即实际上看不见条纹,如下图所示。条纹,如下图所示。当当ad/2,ad/2,直流分量大直流分量大于于1/21/2。去掉零级谱以。去掉零级谱以后像场分布如右图所后像场分布如右图所示。对应物体上亮的示。对应物
16、体上亮的部分变暗、暗的部分部分变暗、暗的部分变亮,实现了变亮,实现了对比反对比反转转。上述讨论说明了。上述讨论说明了利用空间滤波技术,利用空间滤波技术,可以改变成像系统中可以改变成像系统中像场的光分布。像场的光分布。例例 1 1:在在 4 f4 f 系 统 中系 统 中,在在 x x1 1y y1 1平 面 上 放 置 一 正 弦 光 栅平 面 上 放 置 一 正 弦 光 栅,其 振 幅 透 过 率 为其 振 幅 透 过 率 为 t(xt(x)=t)=t0 0+t+t1 1cos(2cos(20 0 x x1 1)(1)(1)在频谱面的中央设置一小圆屏挡住光栅的零级谱,求像的强度分布及可见在频
17、谱面的中央设置一小圆屏挡住光栅的零级谱,求像的强度分布及可见度度 (2)(2)移动小圆屏,挡住光栅的移动小圆屏,挡住光栅的+1+1级谱,像面的强度分布和可见度又如何?级谱,像面的强度分布和可见度又如何?(1)(1)设用振幅为设用振幅为1 1的单色平面波垂直照明物平面,频谱面上的零级斑对应于的单色平面波垂直照明物平面,频谱面上的零级斑对应于物平面上与物平面上与t t0 0项相联系的直流信息,所以挡住零级斑相当于完全通过系统项相联系的直流信息,所以挡住零级斑相当于完全通过系统的物信息为:的物信息为:u(xu(x1 1,y,y1 1)=t)=t1 1cos(2cos(20 0 x x1 1)解解:按
18、一般程序应先求出按一般程序应先求出t(xt(xl l)的频谱,然后求出滤波后的频谱,再作逆傅里的频谱,然后求出滤波后的频谱,再作逆傅里叶变换而求得像。但也可按如下方式考虑:叶变换而求得像。但也可按如下方式考虑:所以输出的信息成为所以输出的信息成为u ui i(x(x3 3,y,y3 3)=t)=t1 1cos(2cos(20 0 x x3 3)物面物面频谱面频谱面像面像面输出图像强度成为输出图像强度成为)2(cos),(30221233xtyxuIii )4cos(1213021xt 可见度为可见度为12/2/2121minmaxminmax ttIIII(2 2)如果挡住)如果挡住+1+1级
19、谱,输出强度又如何变化呢?级谱,输出强度又如何变化呢?我们可以把输入图像的物信息展开我们可以把输入图像的物信息展开)2exp(21)2exp(21)(10110101xjtxjttxt 谱平面上的谱平面上的+1+1级谱与上式中第二项对应,挡住级谱与上式中第二项对应,挡住+1+1级谱后完全通过的物信息为级谱后完全通过的物信息为)2exp(21)(101010 xjttxu 此时输出信息为此时输出信息为)2exp(21),(301033xjttyxui 输出图像的强度分布为输出图像的强度分布为)2cos(41),(),(3010212023333xttttyxuyxIii 除直流分量以外,其交流成
20、分的空间频率仍旧是除直流分量以外,其交流成分的空间频率仍旧是 0 0,但条纹的可见度降为,但条纹的可见度降为202110minmaxminmax4/ttttIIII 例例2 2:在:在4f4f系统中,在系统中,在x x1 1y y1 1平面上有两个图像,它们的中心在平面上有两个图像,它们的中心在x x1 1轴上,距离坐标轴上,距离坐标原点分别为原点分别为a a和和-a-a。今在频谱面上放置一正弦光栅,其振幅透过率为。今在频谱面上放置一正弦光栅,其振幅透过率为)2cos(1),(aT 试证明在像面中心可得到两个图像的相加。试证明在像面中心可得到两个图像的相加。解:设用平面波垂直照射物平面,则解:
21、设用平面波垂直照射物平面,则x x1 1y y1 1平面上两个像的复振幅分布为平面上两个像的复振幅分布为),(),(),(11211111yaxuyaxuyxu 其频谱为其频谱为),(U滤波函数滤波函数)2cos(1),(aT 可看成是系统的传递函数。因此像的复振幅为可看成是系统的传递函数。因此像的复振幅为 ),(),(),(133 HUFyxui ),(),(3333yxhyxu ),(33yxh式中式中是是),(H的点扩散函数,即的点扩散函数,即 )2cos(1),(133 aFyxh ),(21),(21),(333333yaxyaxyx ),(21),(21),(),(),(33333
22、33333yaxyaxyxyxuyxui 像的复振幅为像的复振幅为 ),2(),2(21),(),(),(),(21332331332331332331yaxuyaxuyaxuyaxuyxuyxu 由上式可知在像面中心得到图像由上式可知在像面中心得到图像u u1 1和和u u2 2的相加。的相加。例例3 3:在:在4f4f成像系统中,为了在像面上得到输入图像的成像系统中,为了在像面上得到输入图像的图像,试问在频图像,试问在频谱面上应该使用怎样的滤波器?谱面上应该使用怎样的滤波器?解:设输入图像的复振幅为解:设输入图像的复振幅为u u0 0(x(x1 1),其频谱为,其频谱为U U0 0(),因
23、此有),因此有 dxjUxu)2exp()()(1010 又设输出像的复振幅为又设输出像的复振幅为u ui i(x(x3 3),),由题知由题知:dxjUdxdxudxdxui)2exp()()()(3033033 dxjUj)2exp()(230 透过变换平面的频谱应为透过变换平面的频谱应为)(2)()()(000 UjUTU 2)(jT 滤波函数为滤波函数为可取正、负两值。为实现负值,可将两块模片叠合,一块是振幅模片,可取正、负两值。为实现负值,可将两块模片叠合,一块是振幅模片,其透过率为其透过率为 2)(1 T另一块是另一块是相位模片相位模片,做成在,做成在的正范围和负范围中,其相位差为
24、的正范围和负范围中,其相位差为 的相位的相位掩模,其透过率函数为掩模,其透过率函数为)(2 Tjj 0 0 2j=2 jj 微分运算的微分运算的j2j2滤波器滤波器6.2 6.2 系统与滤波器系统与滤波器 空间滤波系统需要完成从空域到频域,又从频域还原到空域的两次空间滤波系统需要完成从空域到频域,又从频域还原到空域的两次傅里叶变换,以及在频域的乘法运算。典型的滤波系统是三透镜系统傅里叶变换,以及在频域的乘法运算。典型的滤波系统是三透镜系统,这种系统,这种系统简称为简称为4f4f系统。系统。输入和输出的关系为输入和输出的关系为 ),(),(),(),(),(333300133yxhyxgHGFy
25、xg 式中式中g g0 0(x(x3 3,y,y3 3)物体的几何像,物体的几何像,h h是是H H的逆傅里叶变换,的逆傅里叶变换,称为滤波器的脉冲称为滤波器的脉冲响应。响应。6.2.16.2.1空间滤波系统空间滤波系统 在光学处理系统中,空间滤波器是位于空间频率平面上的一种在光学处理系统中,空间滤波器是位于空间频率平面上的一种模片模片,用于实现对输入信息的某种变换。空间波波器的透过率一般是复函数用于实现对输入信息的某种变换。空间波波器的透过率一般是复函数:),(exp),(),(jAH 6.2.2 6.2.2 空间滤波器空间滤波器根据透过率函数的性质根据透过率函数的性质,空间滤波器可以分为以
26、下几种空间滤波器可以分为以下几种:1.1.二元振幅滤波器二元振幅滤波器 这种滤波器的复振幅透过率是这种滤波器的复振幅透过率是0 0或或1 1。由二元滤波所作用的频率区间又。由二元滤波所作用的频率区间又可细分为可细分为:低通滤波器:低通滤波器:它只允许位于频谱面中心及其附近的低频分量通过,它只允许位于频谱面中心及其附近的低频分量通过,可以用来滤掉高频噪声。可以用来滤掉高频噪声。高通滤波器高通滤波器:它阻挡低频分量而允许高频通过,可以实现图像的衬度反它阻挡低频分量而允许高频通过,可以实现图像的衬度反转或边缘增强。转或边缘增强。带通滤波器带通滤波器:它只允许特定区间的空间频谱通过,可以去除随机噪声。
27、它只允许特定区间的空间频谱通过,可以去除随机噪声。方向滤波器方向滤波器:它阻挡它阻挡(或允许或允许)特定方向上的频谱分量通过,可以突出某特定方向上的频谱分量通过,可以突出某些方向性特征。些方向性特征。上述四种二元振幅滤波器的形状如下图所示。上述四种二元振幅滤波器的形状如下图所示。四种二元振幅滤波器四种二元振幅滤波器 方向滤波器 印刷电路中掩模疵点的检查 组合照片上接缝的去除 地震记录中强信号的提取 这种滤波器仅改变各频率成分的相对振幅分布,而不改变其相位分布。这种滤波器仅改变各频率成分的相对振幅分布,而不改变其相位分布。2.2.振幅滤波器振幅滤波器 它只改变空间频谱的相位,不改变它的振幅分布。
28、由于不衰减入射它只改变空间频谱的相位,不改变它的振幅分布。由于不衰减入射光场的能量,具有很高的衍射效率。这种滤波器通常用真空镀膜的方法光场的能量,具有很高的衍射效率。这种滤波器通常用真空镀膜的方法得到,但由于工艺方法的限制,要得到复杂的相位变化是很困难的。得到,但由于工艺方法的限制,要得到复杂的相位变化是很困难的。3.3.相位滤波器相位滤波器 这种滤波器对各种频率成分的振幅和相位都同时起调制作用,滤波这种滤波器对各种频率成分的振幅和相位都同时起调制作用,滤波函数是函数是复函数。复函数。它的应用很广泛,但难于制造。它的应用很广泛,但难于制造。19631963年范德拉格特用年范德拉格特用全全息方法
29、息方法综合出复数空间滤波器,综合出复数空间滤波器,19651965年罗曼和布劳恩用年罗曼和布劳恩用计算全息技术计算全息技术制制作成复数滤波器,从而克服了制作空间滤波器的重大障碍。作成复数滤波器,从而克服了制作空间滤波器的重大障碍。4.4.复数滤波器复数滤波器 在一般情况下在一般情况下,用显微镜只能观察物体的亮暗变化用显微镜只能观察物体的亮暗变化,不能辨别物体不能辨别物体相位相位的变化的变化。最初相位物体最初相位物体(如细菌标本如细菌标本)的观察必须采用的观察必须采用染色法染色法,但染但染色的同时会杀死细菌色的同时会杀死细菌,改变标本的原始结构改变标本的原始结构,从而不能在显微镜下如实从而不能在
30、显微镜下如实研究标本的生命过程研究标本的生命过程。19351935年策尼克提出的相衬显微镜年策尼克提出的相衬显微镜,利用相位滤波利用相位滤波器将物体的相位变化转换成可以观察到的光的强弱变化器将物体的相位变化转换成可以观察到的光的强弱变化。这种转换通常这种转换通常又称为幅相变换又称为幅相变换。一、策尼克相衬显微镜一、策尼克相衬显微镜P1P1P2P2P3P36.3 6.3 空间滤波应用空间滤波应用P1P1P2P2P3P3如图透明相位物体放在如图透明相位物体放在P1P1面上,其复振幅透过率为面上,其复振幅透过率为 ),(exp),(1111yxjyxt 假定相移假定相移 11,则可忽略,则可忽略 2
31、 2及更高级的项,于是复振幅透过率可以写为及更高级的项,于是复振幅透过率可以写为:),(1),(1111yxjyxt P1P1P2P2P3P3物光波实际上可看作两部分,强的直接透射光和由于相位起伏造成的弱物光波实际上可看作两部分,强的直接透射光和由于相位起伏造成的弱衍射光,一个普通的显微镜对上述物体所生成的像,其强度可以写成衍射光,一个普通的显微镜对上述物体所生成的像,其强度可以写成1),(1),(21111 yxjyxI 策尼克认识到策尼克认识到,衍射光之所以观察不到衍射光之所以观察不到,是由于它与很强的本底之间是由于它与很强的本底之间相差相差90900 0,只有改变这两部分之间的相位正交关
32、系,只有改变这两部分之间的相位正交关系,才能使两部分光才能使两部分光叠加时产生干涉叠加时产生干涉,从而产生可观察的像强度变化从而产生可观察的像强度变化。直接透射光在谱面直接透射光在谱面上将会聚成轴上的一个焦点上将会聚成轴上的一个焦点,而衍射光由于包含较高的空间频率而在而衍射光由于包含较高的空间频率而在谱面上较为分散谱面上较为分散。由于由于,因此因此可以简单地在谱面放置相位滤波器可以简单地在谱面放置相位滤波器。使零频的相位使零频的相位相对于其它频率的相对于其它频率的相位改变相位改变/2/2。滤波函数为滤波函数为:),(H0,j其它,1),(H0,j其它,1滤波后的频谱滤波后的频谱),(),(),
33、(),(jjFH像面复振幅分布为像面复振幅分布为),(),(3333yxjjyxg 像强度分布为像强度分布为23333),(),(yxjjyxI ),(2133yx 于是像的强度和相移成线性关系于是像的强度和相移成线性关系。即可通过像的强度观察物体的位相分即可通过像的强度观察物体的位相分布。取正号时布。取正号时,相位最大的部位光强最强相位最大的部位光强最强,叫做叫做正相衬正相衬:取负号时取负号时,相位相位值大的部位光强弱值大的部位光强弱,叫做负相衬。叫做负相衬。由于直接透射光相对于衍射光太强由于直接透射光相对于衍射光太强,像的对比度很低像的对比度很低。如果使零级衍如果使零级衍射光产生相移的同时
34、射光产生相移的同时,受到部分衰减受到部分衰减,可以提高像衬度可以提高像衬度,更有利于观察更有利于观察。相衬显微镜是空间滤波技术早期最成功的应用,根据该原理设计的显微镜相衬显微镜是空间滤波技术早期最成功的应用,根据该原理设计的显微镜为生物实验室广泛采用,对于研究有机体的生命机能提供了有力的工具,为生物实验室广泛采用,对于研究有机体的生命机能提供了有力的工具,策尼克因此在策尼克因此在19531953年获诺贝尔奖年获诺贝尔奖。硅藻、纤维质和细菌的相位显微镜照片硅藻、纤维质和细菌的相位显微镜照片二、补偿滤波器二、补偿滤波器 提高光学系统的成像质量始终是光学工作者所追求的目标。提高光学系统的成像质量始终
35、是光学工作者所追求的目标。5050年代初期,麦尔查年代初期,麦尔查(MarechaMarecha)认为照片中的缺陷是由于成像系认为照片中的缺陷是由于成像系统的统的传递函数传递函数中存在相应缺陷引起的中存在相应缺陷引起的,因而如果能在频谱平面上因而如果能在频谱平面上放置适当的滤波器放置适当的滤波器,使得滤波器的滤波函数补偿原来系统传递函使得滤波器的滤波函数补偿原来系统传递函数的缺陷数的缺陷,则两者的乘积产生一个较为满意的频率响应则两者的乘积产生一个较为满意的频率响应,于是照于是照片的质量将得到部分改善片的质量将得到部分改善。假定成像缺陷是由于成像系统严重离焦引起的假定成像缺陷是由于成像系统严重离
36、焦引起的,则在几何光则在几何光学近似下学近似下,离焦系统的脉冲响应是一个均匀的圆形光斑离焦系统的脉冲响应是一个均匀的圆形光斑。即点扩即点扩散函数为散函数为式中式中a a为圆形光斑半径为圆形光斑半径,1/a1/a2 2是归一化因子是归一化因子。为求为求得得相应的传递函数相应的传递函数,可可将上式作傅里叶将上式作傅里叶-贝塞尔变换贝塞尔变换,即即 arcircarh21)(aJaH21)(1上式所表达的传递函数的高频损失严重上式所表达的传递函数的高频损失严重,且在且在某中间频率区域传递函数的符号发生反转某中间频率区域传递函数的符号发生反转。)(),(),(HFG 加加 补偿滤波器后补偿滤波器后)(
37、)(),(),(HHFG 5050年代初期年代初期,巴黎大学研究所麦尔查等人采用下图(巴黎大学研究所麦尔查等人采用下图(a)a)所示的组合滤所示的组合滤波器波器,放在放在4f4f系统的频谱面上补偿这个带缺陷的传递函数系统的频谱面上补偿这个带缺陷的传递函数。其中吸收板其中吸收板用来衰减很强的低频峰值用来衰减很强的低频峰值,以便提高像的对比以便提高像的对比,突出细节突出细节。相移板使相移板使H H的第一个负瓣相移的第一个负瓣相移,以纠正以纠正对比反转对比反转。图图 (b)(b)表示原来的以及补偿后表示原来的以及补偿后的传递函数的传递函数,输出图像的像质量而获得改善输出图像的像质量而获得改善。相移板
38、相移板吸收吸收(a)(a)(b)(b)例:用例:用CRT(CRT(阴级射线管阴级射线管)记录一帧图像透明片记录一帧图像透明片。设扫描点之间的间隔为设扫描点之间的间隔为0.2mm0.2mm,图像最高空间频率为图像最高空间频率为10mm10mm-1-1。如欲完全去掉离散扫描点如欲完全去掉离散扫描点,得得到一帧到一帧连续灰阶图像连续灰阶图像,空间滤波器的形状和尺寸应当如何设计空间滤波器的形状和尺寸应当如何设计?输出输出图像的分辨率如何图像的分辨率如何(设傅里叶变换物镜的焦距设傅里叶变换物镜的焦距f=l000mm,f=l000mm,=632.8nm)=632.8nm)。解:扫描点的表达式为解:扫描点的
39、表达式为 mnnyymxxyxf),(),(010111 mnnyymxxyx),(1010100),(100100yyxxcombyx mnnyymxxyxf),(),(010111 其频谱为其频谱为),(),(00 yxcombF mnnyymxxyx),(1010100),(100100yyxxcombyx mnnfyymfxx),(2020 mnynfyxmfxyx),(1020200 由此可见,点状结构的频谱仍然是点状结构,其频谱出现的位置为由此可见,点状结构的频谱仍然是点状结构,其频谱出现的位置为0202ynfyxmfx 点状结构的频谱是高频,所以采用低通滤波将其滤点状结构的频谱是高频,所以采用低通滤波将其滤 掉,低通滤波器的圆孔掉,低通滤波器的圆孔半径为半径为mmxfxr164.3)0,1(02 能传递的最高空间频率为能传递的最高空间频率为mmx/510max线输出图像的分辨率为输出图像的分辨率为5 5线线 /mm/mm。