近代光信息处理第3章课件

第三章非相干光学信息处理3/14/20231第三章非相干光学信息处理3．1杨氏干涉仪和空间相干性3．2非相干像的形成3．3MTF的测量3．4非相干空间滤波3．5迈克耳孙干涉仪和时间相干性3．6傅里叶变换光谱仪3．7投影显示的消像素技术3．8计算层析技术3．9结论3/14/20232光学信息处理第三章非相干光学信息处理由于廉价的激光器的广泛应用，非相干光学信息处理已变得不那么重要了，与相干光学信息处理相比，非相干光学信息处理的优势很少．现在很少有人试图去建立一个非相干光学信息处理器，例如非相干光学相关器(参见第四章)．尽管如此，大量光学仪器仍是采用非相干光或自然光作为光源的，其中大多数是成像仪器，例如照相机、显微镜、望远镜、投影仪、制版设备等等．应当说，常规意义下的成像，也是光学信息处理的重要应用，在这个意义上，非相干光学信息处理的基本概念仍然有必要加以研究，这些概念已广泛应用于非相干光成像之中．3/14/20233光学信息处理光源相干性如果点光源位于轴外，则干涉条纹也将发生位移，亮纹将在W，V等处出现，如图3.2所示．此时仍然有(a+c)＝(b+d)，和(b+f)＝(a+e+)．

如果S1和S2同时存在，将看到两个独立的干涉图样的非相干叠加，因为S1和S2是非相干的。图3．2离轴点光源杨氏干涉仪3/14/20235光学信息处理空间相干性杨氏干涉仪可以用来研究光波的相干性。通过P1和P2两个小孔是否在屏上产生干涉条纹来确定照明这两点的光波是否相干。若屏上出现高反差的条纹，光波就是相干的；若屏上出现均匀的照明，光波就是非相干的；若屏上出现低反差的条纹，光波就是部分相干的。以P1、P2的位置为函数的相干性表征光波在P1、P2

的相干的程度，称为空间相干性。3/14/20236光学信息处理空间相干性的测量我们可以改变P1和P2的间距来测量空间相干性。间距增大时，发生两个效应，一个是条纹间距的变小，另一个是条纹反差度的下降。条纹反差度决定了空间相干性。如果小孔的间距大于某一极限后屏上的条纹不再出现，则称此极限间距为空间相干宽度(spatialcoherencewidth)．在相干光处理系统中，我们总是假定空间相干宽度大于光学系统的横向特征尺度；在非相干光处理系统中，我们总是假定空间相干宽度为零；而在部分相干光处理系统中，假定空间相干宽度大于零，并小于系统的特征尺度。3/14/20237光学信息处理复振幅的脉冲响应仍是h(x,y)，相应的强度分布为|h(x,y)

|2．由于照明光为非相干光，从各个点光源辐射的光波彼此是不相干的，各点光源的像也是彼此不相干的，输出像是输入平面物体上各点的像的强度叠加，其强度分布为

|g(,)|2=∞-∞|f(x,y)|2|h(-x,-y)

|2

dxdy在频域中:GI(u,v)=FI(u,v)HI(u,v)式中GI，FI和HI分别表示|g|2,|f|2和|h|2的傅里叶变换．|h(x,y)|2又称点扩散函数，记为PSF(pointspreadfunction)，而HI(u,v)则称为非相干成像系统的传递函数，简称光学传递函数，简写为OTF(opticaltransferfunction)．2、非相干光的成像过程(非相干光的照明)3/14/20239光学信息处理2、非相干光的成像过程(非相干光的照明)由于H(u,v)是h(x,y)的傅里叶变换，根据傅里叶变换的法则，|h(x,y)|2的傅里叶变换为H(u,v)的自相关，亦即

HI(u,v)=∞-∞H*(p,q)H(p+u,q+v)

dpdq上式表明OTF是CTF的自相关．OTF通常是复函数，可表为OTF=|OTF|exp(i)=MTFexp(i)记MTF=|OTF|.MTF称为调制传递函数(modulationtransferfunction);而相位则记为PTF=,PTF称为相位传递函数(phasetransferfunction).3/14/202310光学信息处理3.3MTF的测量

非相干成像系统的MTF可以借助于输入平面上的余弦光栅来测量.余弦光栅的光强分布为i(x)=1+cos(2pox)(1)设系统的输出为o(x)=1+mcos(2pox+)

(2)式中反差度即调制度m可如下测出

在频域中，输入函数可表为I(p)=(p)+(p-po)/2+(p+po)/2(4)输出信号可写作o(p)=I(p)OTF(p)=OTF(0)(p)+OTF(po)(p-po)/2+OTF(-po)(p+po)/2通常的归一化手续规定OTF(0)=13/14/202311光学信息处理3.4非相干空间滤波在相干光学信息处理系统(4f系统)中，当我们把相干光源(激光)换成非相干光源(钨丝灯)，傅里叶平面上的傅里叶变换图像就消失了，这一情形与杨氏干涉仪类似．这是否意味着我们不能实现空间滤波?答案是否定的。设想在傅氏平面上设置一小窗口滤波器H(u)，系统的CTF=H(u)，而OTF则是CTF的自相关.图3.3滤波平面上的实窗口函数生成的CTF及OTF3/14/202313光学信息处理CTF是高通滤波器，从u=a到u=a+b，但MTF仍是低通滤波器，从u=-b到u=b与a无关

由一组无规则分布的小孔构成的孔径的作用相当于低通滤波器．这样一个滤波器的截止频率可以由针孔的直径导出，相当于b．如果用照相机去拍摄一个场景，该滤波器可以直接加在镜头上，拍得的照片中即不包含高频分量．日常生活的经验告诉我们：当我们缩小照相机的光圈时，拍得的照片的分辨率(也就是“解析度”)下降，但景深加大．3/14/202314光学信息处理非相干VanderLugt相关器Lohmann指出，VanderLugt相关器也能用在非相干光的情形．相干VanderLugt相关器的输出中，相关项为（参见节4.3(14)式）c(,)=∞-∞f(x,y)

g[x-(-b),y-]

dxdy强度分布为

|c(,)|2=|

∞-∞f(x,y)g[x-(-b),y-]

dxdy|2当输入物体用非相干光时，相关项的强度分布为|c(,)|2=∞-∞|f(x,y)|2|g[x-(-b),y-]

|2

dxdy即|f|2和|g|2

的相关．因而当f与g全同时相关峰出现在(b,0)处，也就是相干光处理器的相关峰位置．3/14/202315光学信息处理3.5迈克耳孙干涉仪和时间相干性迈克耳孙干涉仪见图．当两臂长度相等时(a=b)，相干条纹出现．注意两个反射镜应稍微倾斜一点，否则在屏上看不到条纹．若使得反射镜M2沿光轴方向移远，使b>a，干涉条纹的反差就会下降．当(2b-2a)大于一定长度l后，屏上的条纹消失，变成均匀的亮斑，l称光波的相干长度．相干时间定义为=l

/c(式中c为光速)图迈克耳孙干涉仪3/14/202317光学信息处理3.5迈克耳孙干涉仪和时间相干性在迈克耳孙干涉仪中，两个光束能够形成干涉条纹的前提条件是它们到达屏的时间差不大于，或它们的光程差不大于l；否则就不会产生干涉条纹．相干长度l相当于波列的平均长度．因此对于一个给定的时刻，沿光波传播方向相干性度量体现为时间相干性，在垂直于传播方向的截面中相干性度量体现为空间相干性.图迈克耳孙干涉仪3/14/202318光学信息处理3.6傅里叶变换光谱仪考虑屏上的一个点，称为观察点．该点的相位差取决于两光路的光程差p．由图3.5有p=2(b–a)相干叠加的光强度为I(p,)=S()[1+cos(2p/c)]图迈克耳孙干涉仪式中S()是产生干涉前的光强，称初始光强．S()表征了光波中的频率成分含量，正是我们感兴趣的光谱函数．当光程差为p时，在观察点探测到的总光强为

I(p)=∞oI(p,)d

=∞oS()d+∞oS()cos(2p/c)d3/14/202319光学信息处理s(x)=2∞oS()cos(2x)dS()=∞-∞s(x)cxp(i2x)dx设x=p/c,代入上两式得到s(p/c)=2∞oS()cos(2p/c)d=2[I(p)-A]S()=∞-∞s(p/c)cxp(i2p/c)d(p/c)S()=(2/c)∞-∞[I(p)–A]cxp(i2p/c)dp这里I(p)和A都是可测量.最后，我们看到光波的频率分布S()可以由[I(p)–A]的傅里叶逆变换得到，而[I(p)–A]可以用移动反射镜M2的过程中多次抽样测量的数据来充分逼近．与前面的假设相对应，我们设S(-)=S()，S(-)没有物理意义，我们将它略去．3/14/202321光学信息处理3.7投影显示的消像素技术当使用液晶显示器LCD进行投影成像时，LCD上的像素结构就会出现在投影屏上．1、利用小孔滤波消除像素结构由于像素的周期结构，在频率平面上出现一系列傅里叶频项．这正是原始图形与周期结构函数的乘积经傅里叶变换后，在频率平面上形成的图形的谱与函数列阵的卷积．每一个谱项都只是中心谱项在不同位置的“复现”，因此只要在频率平面上放置小孔滤波器，仅让一个谱项(例如零级谱项)通过，就可以消除像素结构，见因3.6．然而在这一过程中大部分能量都被滤波器拦去，输出像十分暗淡．3/14/202322光学信息处理图3.6利用小孔滤波消除像素结构3/14/202323光学信息处理图3.6利用小孔滤波消除像素结构图3.7利用不同相位延迟的相位滤波消除像素结构3/14/202325光学信息处理3.7投影显示的消像素技术2、利用不同相位延迟的相位滤波消除像素结构只要相位的延迟大于相干长度，各谱项间的相干性就被破坏，结果像素结构就消除了，而强度并不受到影响．对于白光光源，由于频谱很宽，典型的相干长度仅10m左右．上述技术显著地改善了液晶投影显示的成像质量．特别是一些分辨率较低的投影仪，利用这一技术改善了像质，使它的投影像看起来好像是高分辨率的投影仪．3/14/202326光学信息处理F(u,v)的获得如果p(x)用图3.10的方法得到，F(u,0)就可以计算出来(沿u轴方向的傅里叶变换).F(u,0)=∞-∞p(x)

cxp(-i2ux)

dx

为了获得更多的信息，可以把坐标轴(u,v)转动一个小角度得到(u’,v’)，相应地把(x,y)旋转角得到(x’,y’)，而物体不转动，见图3.11．于是得到p(x’)，并采用同样的手续获得F(u’,v’=0)；图3.113/14/202329光学信息处理F(u,v)的获得再把坐标系转过,又获取F(u”,v”=0)，总共旋转180o，就得到傅里叶频域中的一系列数据，参见图3.12．图中旋转坐标系测出p(x)，p(x’)，p(x”)，…，由此计算出F(u,v)，物体保持不动．图3.11图3.123/14/202330光学信息处理断层图像f(x,y)的获得根据F(u,v)的数据集合后，就可以通过傅里叶逆变换得到f(x,y)．有两种不同的处理方法：第一种方法:通过内插，得到在直角坐标系中的F(u,v)，然后进行傅里叶逆变换.f(x,y)=∞-∞F(u,v)

cxp[i2(ux+vy)]

dudv第二种方法:直接在极坐标下进行傅里叶变换，因而免去了在傅里叶域中的内插手续．极坐标系下的傅里叶逆变换为:f(x,y)=2o

∞o

F(cos,sin)

.exp[i2(xcos+ysin)]

dd其中的积分用求和代替．3/14/202331光学信息处理图3.13人体的CT图像脊椎和肋骨呈现出明亮的白色，图像左边的大块的组织是肝脏，右上角可看到胃里的液体和气体．3/14/202332光学信息处理3.9结论对比了非相干光与相干光的成像过程与空间滤波；非相干光学传递函数(OTF)为相干光学传递函数(CTF)的自相关；非相干系统的本性就是低频滤波；简要分析了光的相干性；傅里叶变换光谱仪；计算层析术(CT)；消除液晶显示技术中像素结构的方法．3/14/202333光学信息处理3.10图像的相乘和积分设两张透明片的透过率分别为1(x,y)和2(x,y)，系统1

S是均匀非相干光源，经透镜L1均匀照明xy平面。将两张透明片紧贴置于xy平面上，在平面后便可得到两者的乘积:I(x,y)=k1(x,y)2(x,y)式中k是比例常数。透镜L2的作用是将xy平面上的图像成一缩小像投射在小的光电探测器D上，光电流的数值为：I=k-1(x,y)2(x,y)dxdy光电探测器上得到的便是两个图像的积分运算。3/14/202334光学信息处理但是，如果要适时更换透明片，则采用下图所示的系统更为方便。L2的作用是将x1y1平面以放大率M=1成像于x2y2平面上。应该说明的是，置于x1y1上的透明片应该倒置，形成1(-x1,-y1)，原因是L2成像后将使之坐标反转。系统23/14/202335光学信息处理3.11图像的相关和卷积实现图像相关运算方法:运动法、无运动法1、运动法

采用3.10系统2,1仍然反置。令1在x1方向上位移xo，在y1方向上yo，则D的光电流输出将为I=k-1(x-xo,y-yo)2(x,y)dxdy

因为对于一个实函数而言，其共扼函数与其本身是相同的，用1*代替1，上式可看成是两者之间的相关运算，即1

2在(xo,yo)点的值。若使1沿x方向以速度v匀速移动，则光电探测器将得到两者在y=yo处的一维相关运算，它是一个时间的函数I(vt)。3/14/202336光学信息处理若在x方向每扫描一次，图形就向上移动y1的距离，则得到光电流的一维阵列:Im(vt)=k-1*(x-vt,y-ym)2(x,y)dxdy=1

2

上式是完整二维相关运算，它在Y方向是抽样的.

卷积运算的实现只需把x1y1平面上的1置于正方向，则很容易得到两者的卷积：1*2，这里不再详述。（自己证明）3/14/202337光学信息处理2、无运动法光学系统原理：考虑S面上一点(-xS,-yS)发出的光，经L1后成为平行光透过1照明2，照明光强度分布正比于1[-x+(d/f)xS,-y+(d/f)yS]。经2后由L2聚焦到焦面xDyD上。这里假定L1和L2焦距相等。位于xDyD的探测器测得的强度为

IS=k-1[(d/f)xS-x,(d/f)yS–y]2(x,y)dxdy=

1*2

3/14/202338光学信息处理应该看到，以几何光学为基础的非相干处理系统只能处理光的强度分布，即只能处理非负的实函数，在有些应用中会受到很大的限制。另一方面，由于系统完全是根据几何光学原理设计的，对于细节过于丰富的图像，由于衍射效应其内含的高频信息往往会丢失，使得输出结果引入较大的偏差。因此，以几何光学为基础的非相干光学处理系统只能在保证几何光学定理成立的条件下才能使用。3/14/202339光学信息处理3.12

白光光学信息处理技术采用相干光源能使光学系统实现许多复杂的信息处理运算，但相干光学信息处理的相干噪声较大。此外，相干光源通常是昂贵的，并且对光学处理的环境要求非常严格。

非相干光学处理采用横向扩展的光源，没有空间相干性，若同时采用白光，则时间相干性也减少到很小的程度，因此这种处理方法具有噪声低、结构简单的优点。可是，非相干处理系统没有物理上的频谱平面，因而频域综合就比较困难。由于系统的输入和脉冲响应都只能是非负的实函数，这又大大限制了系统所能完成的运算．3/14/202340光学信息处理于是，人们会提出这样一个问题：在光学处理中能否降低对光源相干性的要求，但又同时保持对复振幅的线性运算性质?为了回答这个问题，人们研究了一类新的光学处理方法，称为白光光学处理。

白光光学处理采用宽谱带白光光源，但采用微小的光源尺寸以提高空间相干性，另一方面在输入平面上引入光栅来提高时间相干性，这样既不存在相干噪声，又在某种程度上保留了相干光学处理系统对复振幅进行运算的能力，运算灵活性好。由于采用宽谱带光源，特别适合于处理彩色图像，近年来受到愈来愈多的重视。将白光光学处理归人非相干光学处理一章，仅仅是从它采用了非相干光源这一角度考虑，我们应该注意到，它与通常所说的非相干光学处理是明显不同的。3/14/202341光学信息处理1.白光光学处理的基本原理白光光学处理系统如图所示．其中S是白光点光源或者白光光源照明的小孔，这一系统类似于相干光学处理的4f系统。但在白光处理中，通常物函数均用光栅抽样(调制)后才放入输入面上，通过对频谱面上色散的物频谱作处理，实现对物函数的处理。3/14/202342光学信息处理令输入透明片的复振幅透过率为t(x1,y1)，与输入透明片紧贴的正弦光栅为

tg(x1)=1+cos(2ox1)式中o为光栅频率，并假定物透明片对照明光源中各种波长的光波的振幅透过率相同。则经光栅抽样后的复振幅分布为

f(x1,y1)=t(x1,y1)[1+cos(2ox1)]对某一确定的波长，在消色差变换透镜L2后焦面P2的空间频谱为F(,)=T(,)*[(,)+(-o,)/2+(+o,)/2]=T(,)+T(-o,)/2+T(+o,)/2(1)3/14/202343光学信息处理

利用P2平面上频率坐标与空间坐标的关系：=x2/f,=y2/f,方程(1)可写为F(x2,y2;)=T(x2/f,y2/f)+T(x2/f-o,y2/f)/2+T(x2/f+o,y2/f)/2(2)从(2)式看到：第一项为零级物谱，而且不同波长的零级物谱的中心位置是相同的；第二项和第三项是±1级信号谱带，每个谱带中心在x2=±fo处，色散为彩虹颜色．对于波长间隔为的两种色光，其一级谱中心在x2轴上的偏移量是x2=fo.假定信号的空间频带宽度为Wt，则不同波长的物谱能够分离的条件是

/>>Wt/

o(3)式中，为两种色光的平均波长．3/14/202344光学信息处理显然，只要光栅频率o远大于输入信号带宽，就可以忽略各波长频谱间的重叠，从而在+1级或—-1级谱面，象相干处理那样，对一系列的波长进行滤波操作．对于某一确定波长n来说，若设滤波函数为Hn(x2/nf-o,y2/nf),则经过滤波和L3的逆傅里叶变换后,如同相干处理那样，在输出平面上波长为n的像场复振幅为gn(x3,y3;n)=F-1{T(x2/nf-o,y2/nf)Hn(x2/nf-o,y2/nf)}忽略与强度分布无关的量，输出面上波长为n的像强度分布为

I(x3,y3;n)=t(x3,y3)*hn(x3,y3;n)2式中，hn是Hn的逆傅里叶变换．3/14/202345光学信息处理实际上滤波器Hn总不可能做到只让n的光波通过，至少包含n的某一波长间隔n的光波都能通过．当然，当n比n小得多时，可以作为准单色处理．考虑到这一点，可以把通过滤波后在像平面上的像强度分布写成

In=n

t(x3,y3)*hn(x3,y3;n)2式中，hn是第n个滤波器的脉冲响应．当有N个离散的滤波器同时作用于频谱面时，由于不同波长的色光是不相干的，因而输出面上得到的是不同波长输出的非相干叠加，即I(x3,y3)=∑n

t(x3,y3)*hn(x3,y3;n)2

3/14/202346光学信息处理从上述分析可以看出，白光处理技术的确能够处理复振幅信号，并且由于输出强度是互不相干的窄带光强度之和，因而又能抑制令人讨厌的相干噪声．应该指出，我们采用的分析方法是对确定波长的处理看作相干光处理，而对不同波长处理后像的叠加又看成是完全非相干的，这在理论上是不严格的，更严格的讨论涉及到部分相干理论．尽管如此，在很多实际应用中，我们只涉及少数几个分离的波长(例如红、绿、蓝三原色)，此时若在信号频谱后加滤色片，还可以进一步改善时间相干性．而且在采用矩形光栅时，由于光栅的多级衍射，在各个频谱上都可以进行滤波操作．对于这一类问题的处理，上述的近似分析已经足够了．实际上，(3)式的条件对很多应用是过份严格了．3/14/202347光学信息处理2.实时假彩色编码

白光信息处理系统对不同波长的单色光，提供了类似于相干光处理系统的运算能力，采用宽带光源使系统可以使用不同的色通道，有利于对图像进行彩色化处理．这里介绍两种图像假彩色编码的方法：等密度假彩色编码和等空间频率假彩色编码．这两种方法都不需要对输入的图像透明片进行预处理，而只需要在白光信息处理系统的频谱面上放置适当的滤波器，就可以在输出平面上直接得到彩色化的图像．由于具有实时处理的特点，因而又称为实时假彩色编码．3/14/202348光学信息处理(1)等空间频率假彩色编码将一复振幅透过率t(x1,y1)的黑白透明片与正交光栅一起放入白光处理系统的输入平面P1处，为分析简便起见，假定正交光栅在两个正交方向上是相加性的，其振幅透过率可以记为tg(x1,y1)=[1+cos(2ox1)/2+cos(2oy1)/2]式中o，o分别是光栅在x1,y1方向上的空间频率．在频谱面P2上，相应于波长的复振幅分布正比于

F(x2,y2;)=T(x2/f,y2/f)+T(x2/f-o,y2/f)/4+T(x2/f+o,y2/f)/4+T(x2/f,y2/f-o)/4+T(x2/f,y2/f+o)/4

由上述方程可见，沿x2和y2轴共有四个彩虹色信号的一级衍射谱．3/14/202349光学信息处理由于空间滤波只有在沿着垂直于颜色弥散的方向上才有效，所以我们用右图所示的一维空间滤波器来进行假彩色化．图中位于x2轴上蓝色谱带处的是一维低通空间滤波器H1(y2/f)，只让y2方向的低频通过；位于y2轴上蓝色谱带处的是一维低通空间滤波器H1(x2/f)，只让x2方向上的低频通过；位于x2轴上红色谱带处的是一维高通空间滤波器H2(y2/f)，只让y2方向的高频通过；位于y2轴上红色谱带处的是一维高通空间滤波器H2(x2/f)，只让x2方向的高频通过．3/14/202350光学信息处理式中Tb和Tr分别是所选择的蓝色及红色彩色信号谱．在输出面P3上，相应的复振幅分布为于是，平面P2上经过滤波后的谱函数可写为3/14/202351光学信息处理如果光栅的空间频率o及o足够高，则上式可近似地表示为

式中b

和r

是信号的蓝色及红色的光谱宽度；h1及h2分别是H1和H2的点扩散函数。上式表明，两个非相干像在输出平面P3合成彩色编码像，像的低频结构呈蓝色，高频结构呈红色。相等的空间频率结构呈现同一颜色，故称为等空间频率编码。3/14/202352光学信息处理等空间频率假彩色编码3/14/202353光学信息处理(2)等密度假彩色编码如果在P2平面上两个呈彩虹颜色的一级谱处安放如图所示的滤波器，其中红色滤波器是一个简单的红滤色片，另一个绿色滤波器是由一个绿滤色片和绿色频带中心位置的相位滤波器组成。于是，在输出平面上形成红色原像和绿色反转像叠加的结果，使得原图像不同密度的区域呈现不同的颜色。3/14/202354光学信息处理分析：谱平面上放置一个红色全通滤波器，一个绿色滤色片并在其中心加一个相位滤波器组成，其数学表达式为:在白光处理的输出平面P3上的复振幅分布为:谱平面上滤波后的频谱分布为:3/14/202355光学信息处理如果光栅频率足够高，则(10)式可近似地写成：

tgn(x3,y3)=tg(x3,y3)–2〈tg(x3,y3)〉式中tgn(x3,y3)是绿色的对比度反转像，这里的〈tg(x3,y3)〉表示tg(x3,y3)的集平均或系综平均。由于像tr和tgn分别来自光源中不同颜色的光谱带，它们之间是非相干的，所以输出面强度分布是

I(x3,y3)=|g(x3,y3;)|2d

=rIr(x3,y3)+gIrn(x3,y3)Ir(x3,y3)是红色正像，Irn(x3,y3)是绿色负像，r和g分别是红色和绿色的光谱宽度。当这两个像重合在一起时就得到了密度假彩色编码的像。原物中密度最小处呈红色，密度最大处呈绿色，中间部分分别对应粉红、黄、浅绿等颜色，密度相同处出现相同的颜色。3/14/202356光学信息处理等密度假彩色编码3/14/202357光学信息处理3.13相位调制假彩色编码在已研究的实现图像彩色化的许多方法中，相位调制彩色编码方法由于其光强利用率高，色饱和度好，噪声低，操作简单，已在遥感、生物医学、气象等图像处理中得到应用。相位调制彩色编码方法经过对信息的调制和解调，实现了空间强度调制信息与空间波长调制信息的转换。这些技术对于理解光学信息处理的实质具有典型的意义，也是白光光学处理的原理性概念的一种具体应用．密度假彩色编码可分为三个步骤：光栅抽样，漂白处理，白光处理系统中滤波解调．3/14/202358光学信息处理1、光栅抽样将周期为a的朗奇光栅与输入图像重叠在一张复制底片上均匀曝光。设输入图像的密度为Di(x,y)，罗奇光栅的透过率为Ts(x,y)=rect(2x/a)*comb(x/a)/a

(1)则经光栅抽样后所得负片的密度分布D(x,y)为

D(x,y)={[D10-Di(x,y)]-D0}rect(2x/a)*comb(x/a)/a+D0式中Do是底片的灰雾密度，D10是可通过改变曝光条件来控制的常数，是底片的反差系数。这样得到的一张矩形级数光栅其底片光密度可以简记为

D(x,y)=D0Ts为0处D(x,y)=D10-Di(x,y)

Ts为1处(2)3/14/202359光学信息处理(2)漂白处理将经抽样所得到的负片进行漂白处理，并适当控制漂白工艺，可以得到近似满足光程差正比于底片密度的效果，即L(x,y)=Lo=CDoTs=0L(x,y)=Lo=C[D10-Di(x,y)]

Ts=1其相位分布

(x,y)=o=2Lo/Ts=0(x,y)=1=2L1/Ts=1复振幅透过率

exp[i(x,y)]=exp(i2CDo/)Ts=0exp[i(x,y)]=exp{i2C[D10-Di(x,y)]/}Ts=13/14/202360光学信息处理或

最后得到编码的相位光栅的振幅透过率为(3)3/14/202361光学信息处理(3)白光信息处理系统中的滤波解调将编码相位光栅放在白光信息处理系统的输入平面上，设入射单色光强度为A()，则频谱平面(x2,y2)上的复振幅为(4)式中x2=f,x2=f,

将(3)式代人(4)式，计算并化简后得(5)3/14/202362光学信息处理当m=0时，有其对应的强度只与相位差和有关，即如滤波器分别让零级频谱和m级频谱通过，则在输出平面上的复振幅分布为当m≠0时，有(6)3/14/202363光学信息处理若用d表示与相位差相对应的光程差，则(6)式可改写为该式表明对于每一个衍射级次，输出图像的强度随波长和光程差而变化。图(a)和(b)分别为零级和1级的输出强度随光程差d而变化的曲线，其中A()=1，作为参变量。3/14/202364光学信息处理

假如用图中画出的红、绿、蓝三种色光r，g，b

照相，则强度输出是三种色光输出的非相干叠加:I(d)=I(d,r)+I(d,g)+I(d,b)

即得到随d而变化的彩色输出。当采用白光光源时，各色光的非相干叠加变成下列积分：

I(d)=I(d,)d仍然是随d而变化的彩色输出。由于在编码和漂白处理中，已使光程差随输入密度而改变，因此得到按输入图像密度变化的假彩色编码。这种编码方法输出图像的色度丰富，饱和度也很好，在低衍射级次(包括零极)输出的情况下，也能得到彩色化效果很好的输出图像，因此光强度利用率高，图像亮度好。3/14/202365光学信息处理3.14彩色胶片资料存贮

彩色胶片的资料存贮是胶片工业中长期没有得到解决的问题，其主要原因是彩色胶片所用的染料不稳定而造成逐渐褪色。虽然有一些技术可以用来保存彩色像，但都存在着不少问题。最常用的方法是用三个基色滤色片将彩色像保持在三个不同的黑白底片上，然后用三个基色幻灯机来合成彩色像。在这过程中三个像要非常精确地重合，此外这个方法还有两个主要缺点：一是每个彩色片需用的存贮体积是黑白片的三倍；二是彩色像合成系统，非常昂贵和复杂。3/14/202366光学信息处理回顾1906年依伏斯首先报导了从黑白透明片恢复出彩色像的技术。他使用了一个用衍射现象来产生彩色像的幻灯器，其中使用了不同频率或不同方位的光栅。1969年缪勒用了相似的方法即用三色光栅来进行彩色像的编码。在解码时，用三个准单色光源来恢复彩色像。马可夫斯基(1972年)、克罗松和基纳雷(1978年)、杨振寰和母国光等人(1980年)的工作。3/14/202367光学信息处理彩色胶片记录在黑白透明片上的技术记录过程是将彩色胶片通过伦奇光栅依次用三次接触法记录在黑白底片上，每次记录时通过不同的基色滤色片(即红、蓝、绿三色)，而对应的伦奇光栅的方位角分别为0o位置(红色滤色片)、60o位置(蓝色滤色片)、120o位置(绿色滤色片)。3/14/202368光学信息处理如果这三次曝光很准确地记录在照相底片上，这就是一张空间编码的黑白透明片，相应的强度分布可用下式表示

po是光栅的空间频率，Tn表示负像的强度透射率，而Tr、Tb、及Tg分别是在底片上红、蓝、绿三色的曝光量，(x，y)、(x’，y’)及(x”，y”)分别对应于光栅在0o、60o及120o位置时的坐标系统，n1是底片的伽玛值，K1为常数。3/14/202369光学信息处理

可利用接触曝光来得到适合于处理的正像透明片。多重光栅空间调制的正像透明片的强度透过率为

其中K2为比例常数，Tp为正片密度透过率，n2为第二张透明片的伽玛值。为了得到像振幅的线性透过率，我们取n1

n2=2，这样振幅透过率为3/14/202370光学信息处理为了从所得到的黑白编码透明片再现逼真的彩色像，我们采用白光处理技术。如图所示，将该透明片放置在白光处理器的输入平面P1处，于是在空间频谱面P2上对波长为的光强分布为3/14/202371光学信息处理

为简单起见，略去上式中的常数因子，则

其中Tr(p,q),Tb(p’,q’)及Tg(p”,q”)分别是Tr(x,y),Tb(x’,y’)及Tg(x”,t”)的傅里叶变换。而p=(2/f),q=(2/f),p’=(2/f)’,q=(2/f)’,p”=(2/f)”,q