1、第十章 非相干光处理,非相干光学处理是指采用非相干光照明的信息处理方法,系统传递和处理的基本物理量是光场的强度分布.,早期的光学处理多属于非相干光学处理,由于光场的非相干性质,输入函数和脉冲响应都只能是非负的实函数.对于大量双极性质的输入和脉冲响应,处理起来比较困难.,激光出现后,相干系统具有一个物理上的频谱平面,可以实现傅里叶变换运算,大大增加了处理的灵活性.又由于全息术的推动,使相干光学处理的研究极为活跃,一度曾使非相干处理技术相形失色.,但是多年的实践表明,相干处理系统的突出问题是相干噪声严重,导致对系统元件提出较高要求,而非相干处理系统由于其装置简单,又没有相干噪声,因而再度受到广泛的
2、重视.,本章第4节介绍一类新的处理方法,它采用非相干光源照明(白光）,但采取了一些提高空间相干性和时间相干性的措施,从而在某种程度上既保留了相干处理系统对复振幅进行运算的能力,又增加了处理的灵活性,已受到愈来愈多的重视.,10.1 相干与非相干光学处理,1O.1.1 相干与非相干光学处理的比较,我们把一张透明图像片作为一个线性系统的输入,当用相干光照明它时,图片上每一点的复振幅均在其输出面上产生相应的复振幅输出.整个输出图像是这些复振幅的线性叠加,即,也就是合成复振幅满足复振幅叠加原则.然而人眼、感光胶片或其它接收器可感知的是光强,即合成振幅绝对值的平方,10.1.2,对于完全非相干系统,输入
3、图像上各点的光振动是互不相关的,每个点源发出的光是完全独立的,或者说是完全随机的,其振幅和初相位均随时间作随机变化.而观察的强度是对时间的平均效应.这样一来应于(10.1.2)式中的第二项是与时间的关的,在非相干情况下其平均值为零,即有,由此可知，非相干处理系统是强度的线性系统，满足强度叠加原理。,因此,相干光处理与非相干光处理系统的基本区别在于,前者满足复振幅相干叠加,后者满足强度叠加原则。,显然,复振幅可取正负或其它复数值.这样一来，相干光处理系统有可能完成加、减、乘、除、微分和卷积积分等多种运算,特别是能利用透镜的傅里叶变换性质，在特定的频谱平面上提供输入信息的空间频谱，在这个频谱面上安
4、放滤波器，可以方便而巧妙地进行频域综合，实现空间滤波。,而在非相干光学处理系统中,光强只能取正值.故相干光学处理信息的能力比非相干光学处理系统要丰富得多.这就是为什么一般采用相干光而不是非相干光进行信息处理的主要原因.,然而,相干光学处理也有几个固有缺点.,(1) 相干噪声和散斑噪声问题,(2) 输入和输出上存在的问题,在光学系统中(如透镜、反射镜和分束器等)不可避免地存在一些缺陷,如气泡、擦痕以及尘埃、指印或霉斑等.当用相干光照明时,这些缺陷将产生衍射,而这些衍射波之间又会互相干涉,从而形成一系列杂乱条纹与图像重叠在一起,无法分开.这就是所谓相干噪声。,（1）相干噪声和散斑噪声问题.,另外,
5、当用激光照明一个漫射体时,物体表面上各点反射的光在空间相遇而发生干涉.由于漫射物体表面的微观起伏与光波长相比是粗糙的,也是无规的,因而这种干涉也是无规的.当用相干光照明漫射物体时,这个物体看上去总是麻麻点点的,这就是散斑噪声.,由于以上两种噪声的存在,因此相干光处理的图像总是斑纹重叠,结果总不令人满意,有时甚至把信号淹没.噪声问题成了相干光信息处理发展的严重障碍.,由于信息是以光场复振幅分布的形式在系统中传递和处理,这就要求把输入图像制成透明片,然后用激光照明.这就排除了直接使用发光二极管(LED)阵列等作为输入信号的可能性。而在许多实际应用中的信号是以这种方式提供的。现在已广为使用的光学与电
6、子学混合处理系统,可以直接使用这类非相干信号。,(2)输入和输出上存在的问题,(3) 激光是单色性极好的光源,因此,相干处理系统原则上只能处理单色图像,对彩色图像的处理几乎无能为力.,10.1.2 非相干光学处理系统的噪声抑制,非相干光学处理系统对噪声的抑制作用,是从通信理论中的多余通道的概念发展而来的.例如发送某个信号用了N个信息通道(如同时用几路电话通道来传送一个电话),那么第i个通道的输出信号为,式中,ni为第i个通道上的噪声,不同通道上的噪声是不同的;s为信号,它对所有的通道都是相同的.这样,总的强度输出信号为,E表示求平均,由于噪声是完全随机的，其信号的平均值为零，,另外，不同噪声之
7、间互不相关，因此有,式中,由上面分析可知，单一通道上的信噪比为,当引入N个通道后，信噪比为,因此，多余通道的引用使信噪比提高了N倍。,关于这一点在光学系统中是容易理解的。,扩展光源引入的多余通道,用三个互不相干的点光源代表单色空间非相干扩展光源.光源放在准直透镜LI的前焦面上.显然,不同点光源发出的光经准直透镜后,将通过不同的路径到达像面.由图可见,不同路径的光所成的像是相互重叠的,也就是不同通道上的信号是相同的.这就使得光学元件上的尘埃或其表面缺陷对图像的影响微不足道。,例如,在图中的第三通道中,由于透镜表面的尘埃挡掉了来自物体某一部位的信息,但它还可以从另外两个通道传到输出面.另外,即使系
8、统内各处都有尘埃或缺陷,但由于不同的路径所通过的光学系统的区域是不相同的,也就是说不同通道上的噪声分布是不相同的,而这些通道上各光场之间互不相干,故输出平面上的噪声是不同通道上噪声的强度相加,最终的结果就是对噪声的平均.因此,用空间非相干扩展光源可提高输出图像的信噪比.,10.2.1 卷积和相关,10.2 基于几何光学的非相干处理系统,实现两个函数的卷积和相关是光学信息处理中最基本的运算,在相干光学处理系统中这些运算是通过两次傅里叶变换和频城乘法运算完成的.非相干处理系统由于没有物理上的频谱平面,故不能按照同样的方法处理.但是从空域来看,卷积和相关运算都包括位移、相乘、积分三个基本步骤,采用非
9、相干成像系统也可以完成这些运算.,10.2 基于几何光学的非相干处理系统,若把强度透射率为tl的一张透明片在强度透过率为t2的另一张透明片上成像,那么在第二张透明片后面每点的光强都正比于乘积tl t2.所以用光电探测器来测量透过两块透明片的总强度时,给出的光电流I为,(10.2.1),上图是实现这一运算的系统.透镜L2将tl以相等大小成像在t2上,而透镜L3将透过t2的一个缩小像投射到探测器上.若使其中一张透明片匀速运动,并把测量的光电流响应作为时间的函数,就可以实现tl和t2的一维卷积.例如,让透明片t2按反射的几何位置放入,使得(10.2.1)式变成,若使t2在x和y方向分别移动x0,y0
10、，则t2(-x,-y)变成为,这时探测器的响应为,显然光电探测器测得的,的值是t1，t2在x=x0,y=y0点的卷积值。,相关运算与卷积运算的区别在于,两个函数之一没有折叠的步骤,所以只要使t2透明片按正向几何位置放入就可实现两者的相关运算.,若使t2沿x和y的负方向移动x0,y0，则t2(x,y)变成t2(x+x0,y+yo),于是光电探测器的响应为,这就是t1，t 2在x=x0,y=y0点的相关值。,10.2.2 无运动元件的卷积和相关运算,图10.2.2实现无运动卷积和相关运算的系统,为了避免机械扫描的麻烦,可以采用图10.2.2所示的系统来实现卷积和相关运算.均匀漫射光源S放在透镜L1
11、的前焦面上,透射率为f(x,y)的透明片紧贴着放在L1之后,在离f(x,y)的距离为d处,并且在透镜L2的前面紧贴放置透明片h(x,y),然后在L2的后焦面上用胶片或二维阵列检测器进行记录.,图10.2.2实现无运动卷积和相关运算的系统,为了解释这个系统的工作原理,考虑由光源上特定一点(-xo,-yo)发出的光,经L1后变成平行光,若把第一张透明片投影到第二张透明片上,则通过L2把光束会聚到探测器的(x0,yo)点,如果假定两个透镜的焦距相同,那么在检测器上的强度分布为,这正是所要求的卷积。,这种系统的优点是简单易行,缺点是对f(x ,y)的空间结构越细,得到的相关值误差就越大.因为从f(x,
12、y)到 h(x,y)完全是按几何投影考虑的,完全忽略了结构的衍射,结构越细,衍射越显著,所以用这个系统处理的图像的分辨率是受到限制的.,这正是所要求的相关。若第一张输入透明片按反射的几何位置放入，则检测器上的强度分布为,以几何光学为基础的非相干处理系统有两个明显的限制:,一个是由于照明的非相干性质,系统传递和处理的物理量只能是非负的强度分布,给处理双极性信号和综合双极性脉冲响应造成困难。,另一个限制是我们在所有分析过程中均忽略了衍射效应,这实际上是限制了系统处理的信息容量.因为信息容量的增大,意味着透明片上的空间结构变得越来越精细,通过透明片的光就越来越多地被衍射,只剩下越来越少的光遵从几何光
13、学定律,所以输出将偏离按几何关系给出的结果.,10.3基于衍射的非相干处理非相干频域综合,在相干处理系统中,可以由直接改变变换透镜后焦面上的振幅透过率来综合所需要的滤波运算.当使用非相干光照明时,频域综合仍然是可能的,因为非相干系统的光瞳函数和光学传递函数之间存在着一个简单的自相关函数关系.下图绘出了典型的非相干空间滤波系统。,类似于相干成像系统,输入与输出强度之间的关系可以表示为,hI是强度点扩散函数PSF，上式的归一化傅里叶变换为,为系统的光学传递函数OTF。非相干空间滤波是改变输入光强频谱中各频率余弦分量的对比和相位关系，只要根据输入输出关系,在频城综合出所需的OTF,就可实现各种形式的
14、滤波.,衍射受限系统的OTF等于光瞳函数(即出射光瞳函数,简称光瞳函数)的归一化自相关函数，即,10.3.3,由(10.3.3)式可知,根据系统所需的OTF设计光瞳函数,频域综合可在光瞳面着手。,10.3.3,相干系统中有一个物理上的实实在在的频谱面, 非相干系统中没有这样直接,光瞳函数与传递函数之间通过自相关相联系.,若光瞳面上放置其它形式的滤波器,P应该等于滤波器的透过率函数。对于滤波器的位置精度要求,不像相干系统那么苛刻.但由所需的传递函数确定光瞳函数的解不是唯一的,如何由OTF确定最简单的光瞳函数的步骤现在还不知道.,下面给出非相干频域综合的两个实例,10.3.1 切趾术,在非相干成像
15、系统中,点物在像面上的响应称为点扩散函数。为说明概念,我们考虑一个单透镜成像系统,若孔径光阑紧贴透镜放置,则孔径光阑也是出瞳. 若孔径是圆形,则点扩散函数为圆孔的傅里叶变换.,艾里斑的中央亮斑占有绝大部分能量,根据瑞利判据,系统的分辨率完全决定于中央亮斑半径.次级亮环的峰值仅是中央峰值的0.017,可以忽略它的影响.,圆孔衍射,：艾里斑直径,瑞利判据,对于两个强度相等的不相干的点光源（物点），一个点光源的衍射图样的主极大刚好和另一点光源衍射图样的第一极小相重合，这时两个点光源（或物点）恰为这一光学仪器所分辨.,光学仪器的分辨本领,（两光点刚好能分辨）,但是,这个分辨率判据仅适合于分辨两个等强度
16、光点的情况.当两个光点强度的差别与艾里斑中央和次级大相当时,次极大的存在将干扰我们判断较弱光点的存在.例如观测天狼星附近很弱的伴星, 就会遇到这种情况.切趾术就是为了使中央亮斑周围的亮环去掉而采取的一种非相干频域的综合技术.,光瞳函数,点扩散函数,MTF,由于光瞳边界透过率呈阶跃变化, 导致次级衍射环的产生。要切去点扩展函数的趾部(次级亮环),应把光瞳的透过率分布改为缓变形式.例如采用高斯型透过率孔径函数(光瞳函数),由于高斯型孔径的夫琅禾费衍射图样仍是高斯型的(即高斯函数的傅里叶变换仍是高斯函数),故点扩散函数仍是高斯型分布,能够满意地消除次级环的影响.,光瞳函数,点扩散函数,MTF,从OT
17、F的观点看,这是增大低频的调制传递函数(MTF)值,削弱高频传递能力的结果.上图比较了切趾前后的光瞳函数、点扩散函数和调制传递函数。,下面介绍一个具体的结构.图10.3.4表示一个望远物镜L,孔径光阑P紧贴物镜放置.被观察的远方物体在其后焦面上产生的像乃是孔径函数的夫琅禾费衍射图样,其强度分布如图10.3.3(b)中的实线所示.为了既不增加孔径P,又使中央亮斑之外的次极大被切掉,可在孔径P上放入一块玻璃制成的很薄的平行平板Q,在其上镀以非均匀的吸收膜层,使它的振幅透过率从中心到边缘逐渐减小,呈高斯分布曲线变化.,如图 (a)的虚线所示这样,孔面上光场的分布就由原来的均匀分布变成了高斯分布，所以
18、后焦面上的衍射斑也就是高斯函数的傅里叶变换了,它仍然是高斯分布,如图 (b)中的虚线所示，中央亮斑的宽度虽然略有变宽，但它的边缘次极大被切掉了。,10.3.2 沃耳特(Wolter)最小强度检出滤波器,图10.3.5 沃耳特最小强度检出滤波器,这是一个在光瞳面上建立适当的相位分布,从而改变系统成像性质的例子.矩形光瞳分成两半,其中一半蒸镀了产生相位差的透明薄膜,如图10.3.5(a)所示.这种情况下,光学系统的点扩散函数为,很窄的暗线,如果用这样的光学系统产生接近于点光源或线光源物体的像，则在像的中心将出现很窄的暗线，用它来测定物体的位置特别有利，这种方法用于摄谱仪，可以求出光谱线的正确位置。
19、,10.4 白光光学信息处理技术,采用相干光源能使光学系统实现许多复杂的信息处理运算,这主要是由于相干光学系统的复振幅处理能力很强.可是,正如盖伯所指出的,相干噪声是光学信息处理的头号敌人,此外,相干光源通常是昂贵的,并且对光学处理的环境要求非常严格。,在光学处理中能否降低对光源相干性的要求，但又同时保持对复振幅的线性运算性质？,为了回答这个问题，人们研究了一类新的光学处理方法，称为白光光学处理。,由于采用宽谱带光源,特别适合于处理彩色图像,近年来受到愈来愈多的重视.将白光光学处理归入非相干光学处理一章,仅仅是从它采用了非相干光源这一角度考虑,我们应该注意到,它与通常所说的非相干光学处理是明显
20、不同的.,白光光学处理采用宽谱带白光光源，但采用微小光源尺寸以提高空间相干性，另一方面在输入平面上引入光栅来提高时间相干性。这样既不存在相干噪声，又在某种程度上保留了相干光学处理系统对复振幅进行运算的能力，运算灵活性好。,10.4.1 白光光学处理的基本原理,白光光学处理器,常用的白光光学处理系统如图下图所示.,其中S是白光点光源或者白光光源照明的小孔,L1为准直透镜,L2和L3是消色差傅里叶变换透镜,P1,P2和P3分别是系统的输入平面、频谱平面和输出平面.这一系统类似于相干光学处理的4f系统.但在白光处理中,通常物函数均用光栅抽样(调制)后才放入输入面上,通过对频谱面上色散的物频谱作处理，
21、实现对物函数的处理。令输入透明片的复振幅透过率为t（x1,y1),与输入透明片紧贴的正弦光栅为,式中,。为光栅频率,并假定物透明片对照明光源中各种波长的光波的振幅透过率相同.,对某一确定的波长,在消色差变换透镜L2后焦面P2的空间频谱为,则经光栅抽样后的复振幅分布为,对某一确定的波长,在消色差变换透镜L2后焦面P2的空间频谱为,上式第一项为零级物谱，而且不同波长的零级物谱的中心位置是相同的；第二和第三项是1级信息谱带，每个谱带中心在,色散为彩虹颜色,对于波长间隔为的两种色光，其一级谱中心在x2轴上的偏移量是,假定信号的空间频带宽度这Wt，则不同波长的物谱能够分离的条件是,不同波长的谱在x2平面
22、上的宽度是与波长有关的,平均宽度为,当,不同波长的物谱能够分离,(10.4.5,从式中可以看出，只要光栅频率0远大于输入信息的带宽wt,可以不考虑各波长频谱间的重叠，从而可以在谱面P2上独立对一种波长或几种离散的波长，像相干处理一样进行滤波操作。,对于某一确定的波长n来说，可以采用滤波函数为,的滤波器进行滤波,为了便于讨论，假定各个滤波器都放在同一衍射级，实现对不同波长的+1级谱的处理，同时挡掉其它衍射级。在输出平面上波长为n的像场复振幅为,在输出平面上波长为n的像场强度分布为,当有N个离散的滤波器同时作用于频谱面时，由于不同波长的光是不相干的，因而输出平面上得到的是不同波长输出的非相干叠加,
23、从上述分析可以看出,白光处理技术的确能够处理复振幅信号,并且由于输出强度是互不相干的窄带光强度之和,因而又能抑制令人讨厌的相干噪声.,应该指出,我们采用的分析方法是对确定波长的处理看作相干光处理;而对不同波长处理后像的叠加又看成是完全非相干的,这在理论上是不严格的,更严格的讨论涉及到部分相干理论.,尽管如此,在很多实际应用中,我们只涉及少数几个分离的波长(例如红、绿、蓝三原色),此时若在信号频谱后加滤色片,还可以进一步改善时间相干性.而且在采用矩形光栅时,由于光栅的多级衍射,在各个频谱上都可以进行滤波操作.对于这一类问题的处理,上述的近似分析己经足够了.实际上,(10.4.5)式的条件对很多应
24、用是过分严格了。,10.4.2 假彩色编码,白光信息处理系统对不同波长的单色光,提供了类似于相干光处理系统的运算能力,采用宽谱带光源使系统可以使用不同的色通道,有利于对图像进行彩色化处理-这里介绍几种图像假彩色编码的方法。,等密度假彩色编码,等空间频率假彩色编码,调制,1.等空间频率假彩色编码,将一复振幅透过率t(x1,yl)的黑白透明片与正交光栅一起放入图10.4.1所示的白光处理系统的输入平面Pl处,为分析简便起见,假定正交光栅在两个正交方向上是相加性的,其振幅透过率可以记为,式中, o、o分别是光栅在xl,y1方向上的空间频率.在频谱面P2上,相应于波长的复振幅正比于,由上述方程可见,沿
25、x2和y2轴共有四个彩虹色信号的一级衍射谱.如下图所示,如图所示安放高、低通滤波器，使某一种颜色的低频成分和另一种颜色的高频成分通过。输出平面上物体的低频和高频信息将呈现出不同的颜色。比如，低频和高频分别为蓝色和红色。如果采用各种带通滤波器让不同的颜色通过，就可以得到色彩更丰富的假彩色图像。,2.等密度假彩色编码,如果在上例中,在P2平面上两个呈彩虹颜色的一级谱处安放如下图所示的波器,其中红色滤波器是一个简单的红滤色片,另一个绿色滤波器是由一个绿滤色片和绿色频带中心位置的相位滤波器组成。于是，在输出平面上形成红色原像和绿色反转像叠加的结果,使得原图像不同密度的区域呈现不同的颜色.原图中密度最小
26、处呈红色，密度最大处呈现绿色，中间部分分别对应粉红、黄、淡绿等颜色，密度相同处出现相同的颜色。,二、 调制,调制是指用不同方位角的光栅分别对输入 图像不同的区域预先进行调制。这样制成的输入 透明片放入4f系统中，若用单色光照明，在频谱面上，与图像各区域相对应的频谱成分出现在不同方位上。用狭缝作滤波器，在不方位角上可以抽取不同区域的像。,当系统用白光点光源照明，不同方位的频谱均呈彩虹颜色。如果滤波器上开一些小孔，在不同方位角上，小孔可选取不同颜色的谱，我们将得到一个假彩色化的输出图像。每个区域对应一种色彩。参见阿贝波特仿真实验,3、 相位调制假彩色编码,将黑白图像假彩色化,并使彩色图像特征鲜明容
27、易识别,一直是图像处理中的一个重要的课题.在已研究的实现图像假彩色的许多方法中,相位调制假彩色编码方法由于其光强利用率高、色饱和度好、操作简便，已获得了较广泛的应用。,该方法分三步进行：光栅调制漂白处理滤波解调它实际上是另一种等密度假彩色编码。参见阿贝波特仿真实验。,1）、光栅调制,将作为输入的黑白图像透明片与周期为d、缝宽为a的罗奇光栅密接，在一张感光底片上均匀曝光。,罗奇光栅透过率,设输入图像的光密度为Di(xi,yi), 入射的曝光量为Ei (xi,yi),则原底片 的强度透过率t (xi,yi) 与Di (xi,yi),的关系为：,底片的光密度与曝光量的关系为,BC曝光不足区,AD曝光
28、过度区, 线性工作区的斜率,D0底片的灰雾密度,输入光强经原片衰减和光栅取样后，到达感光底片的曝光量为,通过底片的光密度为,2）漂白处理,将上述负片进行漂白处理，只要适当控制工艺，就可以满足底片上的光程L（x1,y1)近似正比于底片的光密度，即,相应的相位分布,复振幅透过率为,通过以上步骤，图像的密度信息就转化成相位信息。可将以上结果按矩形相位光栅来处理。上述复振幅透过率可表达为,式中,3）滤波解调,将编码相位光栅放在白光处理系统的输入平面P1上，设入射光强度为I（），则频谱面上的复振幅为相位光栅的傅里叶变换：,对于0级谱,m=0,在P2面上放一个小圆孔,分别让0级和第m级谱通过,则在输出平面
29、P3上的复振幅分布为:,对应的强度为,以上结果表明:对于每一衍射级次,输出图像的强度随波长和光程差的变化而变化。,若光源中只有红绿兰3种色光,则强度输出是3种色光输出的非相干迭加:,这样就得到了随光程差变化的彩色输出.当采用白光光源时,各种色光的非相干迭加用积分表示,输出还是彩色输出.由于编码和漂白过程中,已使光程差随图像的密度而改变,因此,实现了按输入图像密度变化的假彩色编码.其图色度丰富,光强利用率高、图像亮度好,已在遥感、生物医学、气象等图像处理中得到了较为广泛的应用.,2、要设计一个“散焦”的(非相干)空间滤波系统,使得它的传递函数的第一零点落在。周/cm的频率上,假定要进行滤波的数据
30、放在一个直径为L的圆透镜前面2f距离处.问所要求的“误聚焦 距离”为多少？（用f,L和。表示），对于。=10周/mm,和L=2cm, 的数值是多少？,解：,散焦系统的点扩散函数在几何光学近似下（没有考虑衍射）是由于离焦而产生的均匀照相明的圆斑。对于非相干系统而言，像面上的光场是所有脉冲响应对应的强度迭加.,所以系统的光学传递函数为,J1是一阶第一类贝塞尔函数。,按题意。是第一个零点，所以有,因为J1的第一个根约为1.22,所以有,1 图题10.1为一投影式非相干光卷积运算装置,由光源S和散射板D产生均匀的非相干光照明,m(x,y)和O(x,y)是两张透明片,在平面P上可以探测到m(x,y)O(x,y)卷积。 (1)写出此装置的系统点扩散函数. (2)写出P平面上光强分布的表达式。 (3)若m(x,y)的空间宽度为L1,0(x,y)的空是宽度为L2,求卷积的空间宽度.,例题,(1)坐标为（x, y)点的单位强度点源，投射到O(x,y)上而在P平面上造成的分布即为系统的点扩散函数，其中心位置在,利用几何关系,以,为中心的函数O的形式为,再考虑到投影面积每边长的放大倍数,于是系统的点扩散函数为,（2）P平面上点（xd,yd）所接收到的强度响应为物函数m(x,y)与以,为中心的点扩散函数乘积的积分，即,（3）,卷积的空间宽度为,
