近代光信息處理第3章課件

第三章非相干光學(xué)信息處理3/14/20231第三章非相干光學(xué)信息處理3．1楊氏干涉儀和空間相干性3．2非相干像的形成3．3MTF的測量3．4非相干空間濾波3．5邁克耳孫干涉儀和時間相干性3．6傅里葉變換光譜儀3．7投影顯示的消像素技術(shù)3．8計算層析技術(shù)3．9結(jié)論3/14/20232光學(xué)信息處理第三章非相干光學(xué)信息處理由于廉價的激光器的廣泛應(yīng)用，非相干光學(xué)信息處理已變得不那么重要了，與相干光學(xué)信息處理相比，非相干光學(xué)信息處理的優(yōu)勢很少．現(xiàn)在很少有人試圖去建立一個非相干光學(xué)信息處理器，例如非相干光學(xué)相關(guān)器(參見第四章)．盡管如此，大量光學(xué)儀器仍是采用非相干光或自然光作為光源的，其中大多數(shù)是成像儀器，例如照相機、顯微鏡、望遠鏡、投影儀、制版設(shè)備等等．應(yīng)當(dāng)說，常規(guī)意義下的成像，也是光學(xué)信息處理的重要應(yīng)用，在這個意義上，非相干光學(xué)信息處理的基本概念仍然有必要加以研究，這些概念已廣泛應(yīng)用于非相干光成像之中．3/14/20233光學(xué)信息處理光源相干性如果點光源位于軸外，則干涉條紋也將發(fā)生位移，亮紋將在W，V等處出現(xiàn)，如圖3.2所示．此時仍然有(a+c)＝(b+d)，和(b+f)＝(a+e+)．

如果S1和S2同時存在，將看到兩個獨立的干涉圖樣的非相干疊加，因為S1和S2是非相干的。圖3．2離軸點光源楊氏干涉儀3/14/20235光學(xué)信息處理空間相干性楊氏干涉儀可以用來研究光波的相干性。通過P1和P2兩個小孔是否在屏上產(chǎn)生干涉條紋來確定照明這兩點的光波是否相干。若屏上出現(xiàn)高反差的條紋，光波就是相干的；若屏上出現(xiàn)均勻的照明，光波就是非相干的；若屏上出現(xiàn)低反差的條紋，光波就是部分相干的。以P1、P2的位置為函數(shù)的相干性表征光波在P1、P2

的相干的程度，稱為空間相干性。3/14/20236光學(xué)信息處理空間相干性的測量我們可以改變P1和P2的間距來測量空間相干性。間距增大時，發(fā)生兩個效應(yīng)，一個是條紋間距的變小，另一個是條紋反差度的下降。條紋反差度決定了空間相干性。如果小孔的間距大于某一極限后屏上的條紋不再出現(xiàn)，則稱此極限間距為空間相干寬度(spatialcoherencewidth)．在相干光處理系統(tǒng)中，我們總是假定空間相干寬度大于光學(xué)系統(tǒng)的橫向特征尺度；在非相干光處理系統(tǒng)中，我們總是假定空間相干寬度為零；而在部分相干光處理系統(tǒng)中，假定空間相干寬度大于零，并小于系統(tǒng)的特征尺度。3/14/20237光學(xué)信息處理復(fù)振幅的脈沖響應(yīng)仍是h(x,y)，相應(yīng)的強度分布為|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，F(xiàn)I和HI分別表示|g|2,|f|2和|h|2的傅里葉變換．|h(x,y)|2又稱點擴散函數(shù)，記為PSF(pointspreadfunction)，而HI(u,v)則稱為非相干成像系統(tǒng)的傳遞函數(shù)，簡稱光學(xué)傳遞函數(shù)，簡寫為OTF(opticaltransferfunction)．2、非相干光的成像過程(非相干光的照明)3/14/20239光學(xué)信息處理2、非相干光的成像過程(非相干光的照明)由于H(u,v)是h(x,y)的傅里葉變換，根據(jù)傅里葉變換的法則，|h(x,y)|2的傅里葉變換為H(u,v)的自相關(guān)，亦即

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

dpdq上式表明OTF是CTF的自相關(guān)．OTF通常是復(fù)函數(shù)，可表為OTF=|OTF|exp(i)=MTFexp(i)記MTF=|OTF|.MTF稱為調(diào)制傳遞函數(shù)(modulationtransferfunction);而相位則記為PTF=,PTF稱為相位傳遞函數(shù)(phasetransferfunction).3/14/202310光學(xué)信息處理3.3MTF的測量

非相干成像系統(tǒng)的MTF可以借助于輸入平面上的余弦光柵來測量.余弦光柵的光強分布為i(x)=1+cos(2pox)(1)設(shè)系統(tǒng)的輸出為o(x)=1+mcos(2pox+)

(2)式中反差度即調(diào)制度m可如下測出

在頻域中，輸入函數(shù)可表為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通常的歸一化手續(xù)規(guī)定OTF(0)=13/14/202311光學(xué)信息處理3.4非相干空間濾波在相干光學(xué)信息處理系統(tǒng)(4f系統(tǒng))中，當(dāng)我們把相干光源(激光)換成非相干光源(鎢絲燈)，傅里葉平面上的傅里葉變換圖像就消失了，這一情形與楊氏干涉儀類似．這是否意味著我們不能實現(xiàn)空間濾波?答案是否定的。設(shè)想在傅氏平面上設(shè)置一小窗口濾波器H(u)，系統(tǒng)的CTF=H(u)，而OTF則是CTF的自相關(guān).圖3.3濾波平面上的實窗口函數(shù)生成的CTF及OTF3/14/202313光學(xué)信息處理CTF是高通濾波器，從u=a到u=a+b，但MTF仍是低通濾波器，從u=-b到u=b與a無關(guān)

由一組無規(guī)則分布的小孔構(gòu)成的孔徑的作用相當(dāng)于低通濾波器．這樣一個濾波器的截止頻率可以由針孔的直徑導(dǎo)出，相當(dāng)于b．如果用照相機去拍攝一個場景，該濾波器可以直接加在鏡頭上，拍得的照片中即不包含高頻分量．日常生活的經(jīng)驗告訴我們：當(dāng)我們縮小照相機的光圈時，拍得的照片的分辨率(也就是“解析度”)下降，但景深加大．3/14/202314光學(xué)信息處理非相干VanderLugt相關(guān)器Lohmann指出，VanderLugt相關(guān)器也能用在非相干光的情形．相干VanderLugt相關(guān)器的輸出中，相關(guān)項為（參見節(jié)4.3(14)式）c(,)=∞-∞f(x,y)

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

dxdy強度分布為

|c(,)|2=|

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

dxdy|2當(dāng)輸入物體用非相干光時，相關(guān)項的強度分布為|c(,)|2=∞-∞|f(x,y)|2|g[x-(-b),y-]

|2

dxdy即|f|2和|g|2

的相關(guān)．因而當(dāng)f與g全同時相關(guān)峰出現(xiàn)在(b,0)處，也就是相干光處理器的相關(guān)峰位置．3/14/202315光學(xué)信息處理3.5邁克耳孫干涉儀和時間相干性邁克耳孫干涉儀見圖．當(dāng)兩臂長度相等時(a=b)，相干條紋出現(xiàn)．注意兩個反射鏡應(yīng)稍微傾斜一點，否則在屏上看不到條紋．若使得反射鏡M2沿光軸方向移遠，使b>a，干涉條紋的反差就會下降．當(dāng)(2b-2a)大于一定長度l后，屏上的條紋消失，變成均勻的亮斑，l稱光波的相干長度．相干時間定義為=l

/c(式中c為光速)圖邁克耳孫干涉儀3/14/202317光學(xué)信息處理3.5邁克耳孫干涉儀和時間相干性在邁克耳孫干涉儀中，兩個光束能夠形成干涉條紋的前提條件是它們到達屏的時間差不大于，或它們的光程差不大于l；否則就不會產(chǎn)生干涉條紋．相干長度l相當(dāng)于波列的平均長度．因此對于一個給定的時刻，沿光波傳播方向相干性度量體現(xiàn)為時間相干性，在垂直于傳播方向的截面中相干性度量體現(xiàn)為空間相干性.圖邁克耳孫干涉儀3/14/202318光學(xué)信息處理3.6傅里葉變換光譜儀考慮屏上的一個點，稱為觀察點．該點的相位差取決于兩光路的光程差p．由圖3.5有p=2(b–a)相干疊加的光強度為I(p,)=S()[1+cos(2p/c)]圖邁克耳孫干涉儀式中S()是產(chǎn)生干涉前的光強，稱初始光強．S()表征了光波中的頻率成分含量，正是我們感興趣的光譜函數(shù)．當(dāng)光程差為p時，在觀察點探測到的總光強為

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

=∞oS()d+∞oS()cos(2p/c)d3/14/202319光學(xué)信息處理s(x)=2∞oS()cos(2x)dS()=∞-∞s(x)cxp(i2x)dx設(shè)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的過程中多次抽樣測量的數(shù)據(jù)來充分逼近．與前面的假設(shè)相對應(yīng)，我們設(shè)S(-)=S()，S(-)沒有物理意義，我們將它略去．3/14/202321光學(xué)信息處理3.7投影顯示的消像素技術(shù)當(dāng)使用液晶顯示器LCD進行投影成像時，LCD上的像素結(jié)構(gòu)就會出現(xiàn)在投影屏上．1、利用小孔濾波消除像素結(jié)構(gòu)由于像素的周期結(jié)構(gòu)，在頻率平面上出現(xiàn)一系列傅里葉頻項．這正是原始圖形與周期結(jié)構(gòu)函數(shù)的乘積經(jīng)傅里葉變換后，在頻率平面上形成的圖形的譜與函數(shù)列陣的卷積．每一個譜項都只是中心譜項在不同位置的“復(fù)現(xiàn)”，因此只要在頻率平面上放置小孔濾波器，僅讓一個譜項(例如零級譜項)通過，就可以消除像素結(jié)構(gòu)，見因3.6．然而在這一過程中大部分能量都被濾波器攔去，輸出像十分暗淡．3/14/202322光學(xué)信息處理圖3.6利用小孔濾波消除像素結(jié)構(gòu)3/14/202323光學(xué)信息處理圖3.6利用小孔濾波消除像素結(jié)構(gòu)圖3.7利用不同相位延遲的相位濾波消除像素結(jié)構(gòu)3/14/202325光學(xué)信息處理3.7投影顯示的消像素技術(shù)2、利用不同相位延遲的相位濾波消除像素結(jié)構(gòu)只要相位的延遲大于相干長度，各譜項間的相干性就被破壞，結(jié)果像素結(jié)構(gòu)就消除了，而強度并不受到影響．對于白光光源，由于頻譜很寬，典型的相干長度僅10m左右．上述技術(shù)顯著地改善了液晶投影顯示的成像質(zhì)量．特別是一些分辨率較低的投影儀，利用這一技術(shù)改善了像質(zhì)，使它的投影像看起來好像是高分辨率的投影儀．3/14/202326光學(xué)信息處理F(u,v)的獲得如果p(x)用圖3.10的方法得到，F(xiàn)(u,0)就可以計算出來(沿u軸方向的傅里葉變換).F(u,0)=∞-∞p(x)

cxp(-i2ux)

dx

為了獲得更多的信息，可以把坐標(biāo)軸(u,v)轉(zhuǎn)動一個小角度得到(u’,v’)，相應(yīng)地把(x,y)旋轉(zhuǎn)角得到(x’,y’)，而物體不轉(zhuǎn)動，見圖3.11．于是得到p(x’)，并采用同樣的手續(xù)獲得F(u’,v’=0)；圖3.113/14/202329光學(xué)信息處理F(u,v)的獲得再把坐標(biāo)系轉(zhuǎn)過,又獲取F(u”,v”=0)，總共旋轉(zhuǎn)180o，就得到傅里葉頻域中的一系列數(shù)據(jù)，參見圖3.12．圖中旋轉(zhuǎn)坐標(biāo)系測出p(x)，p(x’)，p(x”)，…，由此計算出F(u,v)，物體保持不動．圖3.11圖3.123/14/202330光學(xué)信息處理斷層圖像f(x,y)的獲得根據(jù)F(u,v)的數(shù)據(jù)集合后，就可以通過傅里葉逆變換得到f(x,y)．有兩種不同的處理方法：第一種方法:通過內(nèi)插，得到在直角坐標(biāo)系中的F(u,v)，然后進行傅里葉逆變換.f(x,y)=∞-∞F(u,v)

cxp[i2(ux+vy)]

dudv第二種方法:直接在極坐標(biāo)下進行傅里葉變換，因而免去了在傅里葉域中的內(nèi)插手續(xù)．極坐標(biāo)系下的傅里葉逆變換為:f(x,y)=2o

∞o

F(cos,sin)

.exp[i2(xcos+ysin)]

dd其中的積分用求和代替．3/14/202331光學(xué)信息處理圖3.13人體的CT圖像脊椎和肋骨呈現(xiàn)出明亮的白色，圖像左邊的大塊的組織是肝臟，右上角可看到胃里的液體和氣體．3/14/202332光學(xué)信息處理3.9結(jié)論對比了非相干光與相干光的成像過程與空間濾波；非相干光學(xué)傳遞函數(shù)(OTF)為相干光學(xué)傳遞函數(shù)(CTF)的自相關(guān)；非相干系統(tǒng)的本性就是低頻濾波；簡要分析了光的相干性；傅里葉變換光譜儀；計算層析術(shù)(CT)；消除液晶顯示技術(shù)中像素結(jié)構(gòu)的方法．3/14/202333光學(xué)信息處理3.10圖像的相乘和積分設(shè)兩張透明片的透過率分別為1(x,y)和2(x,y)，系統(tǒng)1

S是均勻非相干光源，經(jīng)透鏡L1均勻照明xy平面。將兩張透明片緊貼置于xy平面上，在平面后便可得到兩者的乘積:I(x,y)=k1(x,y)2(x,y)式中k是比例常數(shù)。透鏡L2的作用是將xy平面上的圖像成一縮小像投射在小的光電探測器D上，光電流的數(shù)值為：I=k-1(x,y)2(x,y)dxdy光電探測器上得到的便是兩個圖像的積分運算。3/14/202334光學(xué)信息處理但是，如果要適時更換透明片，則采用下圖所示的系統(tǒng)更為方便。L2的作用是將x1y1平面以放大率M=1成像于x2y2平面上。應(yīng)該說明的是，置于x1y1上的透明片應(yīng)該倒置，形成1(-x1,-y1)，原因是L2成像后將使之坐標(biāo)反轉(zhuǎn)。系統(tǒng)23/14/202335光學(xué)信息處理3.11圖像的相關(guān)和卷積實現(xiàn)圖像相關(guān)運算方法:運動法、無運動法1、運動法

采用3.10系統(tǒng)2,1仍然反置。令1在x1方向上位移xo，在y1方向上yo，則D的光電流輸出將為I=k-1(x-xo,y-yo)2(x,y)dxdy

因為對于一個實函數(shù)而言，其共扼函數(shù)與其本身是相同的，用1*代替1，上式可看成是兩者之間的相關(guān)運算，即1

2在(xo,yo)點的值。若使1沿x方向以速度v勻速移動，則光電探測器將得到兩者在y=yo處的一維相關(guān)運算，它是一個時間的函數(shù)I(vt)。3/14/202336光學(xué)信息處理若在x方向每掃描一次，圖形就向上移動y1的距離，則得到光電流的一維陣列:Im(vt)=k-1*(x-vt,y-ym)2(x,y)dxdy=1

2

上式是完整二維相關(guān)運算，它在Y方向是抽樣的.

卷積運算的實現(xiàn)只需把x1y1平面上的1置于正方向，則很容易得到兩者的卷積：1*2，這里不再詳述。（自己證明）3/14/202337光學(xué)信息處理2、無運動法光學(xué)系統(tǒng)原理：考慮S面上一點(-xS,-yS)發(fā)出的光，經(jīng)L1后成為平行光透過1照明2，照明光強度分布正比于1[-x+(d/f)xS,-y+(d/f)yS]。經(jīng)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光學(xué)信息處理應(yīng)該看到，以幾何光學(xué)為基礎(chǔ)的非相干處理系統(tǒng)只能處理光的強度分布，即只能處理非負的實函數(shù)，在有些應(yīng)用中會受到很大的限制。另一方面，由于系統(tǒng)完全是根據(jù)幾何光學(xué)原理設(shè)計的，對于細節(jié)過于豐富的圖像，由于衍射效應(yīng)其內(nèi)含的高頻信息往往會丟失，使得輸出結(jié)果引入較大的偏差。因此，以幾何光學(xué)為基礎(chǔ)的非相干光學(xué)處理系統(tǒng)只能在保證幾何光學(xué)定理成立的條件下才能使用。3/14/202339光學(xué)信息處理3.12

白光光學(xué)信息處理技術(shù)采用相干光源能使光學(xué)系統(tǒng)實現(xiàn)許多復(fù)雜的信息處理運算，但相干光學(xué)信息處理的相干噪聲較大。此外，相干光源通常是昂貴的，并且對光學(xué)處理的環(huán)境要求非常嚴(yán)格。

非相干光學(xué)處理采用橫向擴展的光源，沒有空間相干性，若同時采用白光，則時間相干性也減少到很小的程度，因此這種處理方法具有噪聲低、結(jié)構(gòu)簡單的優(yōu)點?？墒牵窍喔商幚硐到y(tǒng)沒有物理上的頻譜平面，因而頻域綜合就比較困難。由于系統(tǒng)的輸入和脈沖響應(yīng)都只能是非負的實函數(shù)，這又大大限制了系統(tǒng)所能完成的運算．3/14/202340光學(xué)信息處理于是，人們會提出這樣一個問題：在光學(xué)處理中能否降低對光源相干性的要求，但又同時保持對復(fù)振幅的線性運算性質(zhì)?為了回答這個問題，人們研究了一類新的光學(xué)處理方法，稱為白光光學(xué)處理。

白光光學(xué)處理采用寬譜帶白光光源，但采用微小的光源尺寸以提高空間相干性，另一方面在輸入平面上引入光柵來提高時間相干性，這樣既不存在相干噪聲，又在某種程度上保留了相干光學(xué)處理系統(tǒng)對復(fù)振幅進行運算的能力，運算靈活性好。由于采用寬譜帶光源，特別適合于處理彩色圖像，近年來受到愈來愈多的重視。將白光光學(xué)處理歸人非相干光學(xué)處理一章，僅僅是從它采用了非相干光源這一角度考慮，我們應(yīng)該注意到，它與通常所說的非相干光學(xué)處理是明顯不同的。3/14/202341光學(xué)信息處理1.白光光學(xué)處理的基本原理白光光學(xué)處理系統(tǒng)如圖所示．其中S是白光點光源或者白光光源照明的小孔，這一系統(tǒng)類似于相干光學(xué)處理的4f系統(tǒng)。但在白光處理中，通常物函數(shù)均用光柵抽樣(調(diào)制)后才放入輸入面上，通過對頻譜面上色散的物頻譜作處理，實現(xiàn)對物函數(shù)的處理。3/14/202342光學(xué)信息處理令輸入透明片的復(fù)振幅透過率為t(x1,y1)，與輸入透明片緊貼的正弦光柵為

tg(x1)=1+cos(2ox1)式中o為光柵頻率，并假定物透明片對照明光源中各種波長的光波的振幅透過率相同。則經(jīng)光柵抽樣后的復(fù)振幅分布為

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光學(xué)信息處理

利用P2平面上頻率坐標(biāo)與空間坐標(biāo)的關(guān)系：=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光學(xué)信息處理顯然，只要光柵頻率o遠大于輸入信號帶寬，就可以忽略各波長頻譜間的重疊，從而在+1級或—-1級譜面，象相干處理那樣，對一系列的波長進行濾波操作．對于某一確定波長n來說，若設(shè)濾波函數(shù)為Hn(x2/nf-o,y2/nf),則經(jīng)過濾波和L3的逆傅里葉變換后,如同相干處理那樣，在輸出平面上波長為n的像場復(fù)振幅為gn(x3,y3;n)=F-1{T(x2/nf-o,y2/nf)Hn(x2/nf-o,y2/nf)}忽略與強度分布無關(guān)的量，輸出面上波長為n的像強度分布為

I(x3,y3;n)=t(x3,y3)*hn(x3,y3;n)2式中，hn是Hn的逆傅里葉變換．3/14/202345光學(xué)信息處理實際上濾波器Hn總不可能做到只讓n的光波通過，至少包含n的某一波長間隔n的光波都能通過．當(dāng)然，當(dāng)n比n小得多時，可以作為準(zhǔn)單色處理．考慮到這一點，可以把通過濾波后在像平面上的像強度分布寫成

In=n

t(x3,y3)*hn(x3,y3;n)2式中，hn是第n個濾波器的脈沖響應(yīng)．當(dāng)有N個離散的濾波器同時作用于頻譜面時，由于不同波長的色光是不相干的，因而輸出面上得到的是不同波長輸出的非相干疊加，即I(x3,y3)=∑n

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

3/14/202346光學(xué)信息處理從上述分析可以看出，白光處理技術(shù)的確能夠處理復(fù)振幅信號，并且由于輸出強度是互不相干的窄帶光強度之和，因而又能抑制令人討厭的相干噪聲．應(yīng)該指出，我們采用的分析方法是對確定波長的處理看作相干光處理，而對不同波長處理后像的疊加又看成是完全非相干的，這在理論上是不嚴(yán)格的，更嚴(yán)格的討論涉及到部分相干理論．盡管如此，在很多實際應(yīng)用中，我們只涉及少數(shù)幾個分離的波長(例如紅、綠、藍三原色)，此時若在信號頻譜后加濾色片，還可以進一步改善時間相干性．而且在采用矩形光柵時，由于光柵的多級衍射，在各個頻譜上都可以進行濾波操作．對于這一類問題的處理，上述的近似分析已經(jīng)足夠了．實際上，(3)式的條件對很多應(yīng)用是過份嚴(yán)格了．3/14/202347光學(xué)信息處理2.實時假彩色編碼

白光信息處理系統(tǒng)對不同波長的單色光，提供了類似于相干光處理系統(tǒng)的運算能力，采用寬帶光源使系統(tǒng)可以使用不同的色通道，有利于對圖像進行彩色化處理．這里介紹兩種圖像假彩色編碼的方法：等密度假彩色編碼和等空間頻率假彩色編碼．這兩種方法都不需要對輸入的圖像透明片進行預(yù)處理，而只需要在白光信息處理系統(tǒng)的頻譜面上放置適當(dāng)?shù)臑V波器，就可以在輸出平面上直接得到彩色化的圖像．由于具有實時處理的特點，因而又稱為實時假彩色編碼．3/14/202348光學(xué)信息處理(1)等空間頻率假彩色編碼將一復(fù)振幅透過率t(x1,y1)的黑白透明片與正交光柵一起放入白光處理系統(tǒng)的輸入平面P1處，為分析簡便起見，假定正交光柵在兩個正交方向上是相加性的，其振幅透過率可以記為tg(x1,y1)=[1+cos(2ox1)/2+cos(2oy1)/2]式中o，o分別是光柵在x1,y1方向上的空間頻率．在頻譜面P2上，相應(yīng)于波長的復(fù)振幅分布正比于

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光學(xué)信息處理由于空間濾波只有在沿著垂直于顏色彌散的方向上才有效，所以我們用右圖所示的一維空間濾波器來進行假彩色化．圖中位于x2軸上藍色譜帶處的是一維低通空間濾波器H1(y2/f)，只讓y2方向的低頻通過；位于y2軸上藍色譜帶處的是一維低通空間濾波器H1(x2/f)，只讓x2方向上的低頻通過；位于x2軸上紅色譜帶處的是一維高通空間濾波器H2(y2/f)，只讓y2方向的高頻通過；位于y2軸上紅色譜帶處的是一維高通空間濾波器H2(x2/f)，只讓x2方向的高頻通過．3/14/202350光學(xué)信息處理式中Tb和Tr分別是所選擇的藍色及紅色彩色信號譜．在輸出面P3上，相應(yīng)的復(fù)振幅分布為于是，平面P2上經(jīng)過濾波后的譜函數(shù)可寫為3/14/202351光學(xué)信息處理如果光柵的空間頻率o及o足夠高，則上式可近似地表示為

式中b

和r

是信號的藍色及紅色的光譜寬度；h1及h2分別是H1和H2的點擴散函數(shù)。上式表明，兩個非相干像在輸出平面P3合成彩色編碼像，像的低頻結(jié)構(gòu)呈藍色，高頻結(jié)構(gòu)呈紅色。相等的空間頻率結(jié)構(gòu)呈現(xiàn)同一顏色，故稱為等空間頻率編碼。3/14/202352光學(xué)信息處理等空間頻率假彩色編碼3/14/202353光學(xué)信息處理(2)等密度假彩色編碼如果在P2平面上兩個呈彩虹顏色的一級譜處安放如圖所示的濾波器，其中紅色濾波器是一個簡單的紅濾色片，另一個綠色濾波器是由一個綠濾色片和綠色頻帶中心位置的相位濾波器組成。于是，在輸出平面上形成紅色原像和綠色反轉(zhuǎn)像疊加的結(jié)果，使得原圖像不同密度的區(qū)域呈現(xiàn)不同的顏色。3/14/202354光學(xué)信息處理分析：譜平面上放置一個紅色全通濾波器，一個綠色濾色片并在其中心加一個相位濾波器組成，其數(shù)學(xué)表達式為:在白光處理的輸出平面P3上的復(fù)振幅分布為:譜平面上濾波后的頻譜分布為:3/14/202355光學(xué)信息處理如果光柵頻率足夠高，則(10)式可近似地寫成：

tgn(x3,y3)=tg(x3,y3)–2〈tg(x3,y3)〉式中tgn(x3,y3)是綠色的對比度反轉(zhuǎn)像，這里的〈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分別是紅色和綠色的光譜寬度。當(dāng)這兩個像重合在一起時就得到了密度假彩色編碼的像。原物中密度最小處呈紅色，密度最大處呈綠色，中間部分分別對應(yīng)粉紅、黃、淺綠等顏色，密度相同處出現(xiàn)相同的顏色。3/14/202356光學(xué)信息處理等密度假彩色編碼3/14/202357光學(xué)信息處理3.13相位調(diào)制假彩色編碼在已研究的實現(xiàn)圖像彩色化的許多方法中，相位調(diào)制彩色編碼方法由于其光強利用率高，色飽和度好，噪聲低，操作簡單，已在遙感、生物醫(yī)學(xué)、氣象等圖像處理中得到應(yīng)用。相位調(diào)制彩色編碼方法經(jīng)過對信息的調(diào)制和解調(diào)，實現(xiàn)了空間強度調(diào)制信息與空間波長調(diào)制信息的轉(zhuǎn)換。這些技術(shù)對于理解光學(xué)信息處理的實質(zhì)具有典型的意義，也是白光光學(xué)處理的原理性概念的一種具體應(yīng)用．密度假彩色編碼可分為三個步驟：光柵抽樣，漂白處理，白光處理系統(tǒng)中濾波解調(diào)．3/14/202358光學(xué)信息處理1、光柵抽樣將周期為a的朗奇光柵與輸入圖像重疊在一張復(fù)制底片上均勻曝光。設(shè)輸入圖像的密度為Di(x,y)，羅奇光柵的透過率為Ts(x,y)=rect(2x/a)*comb(x/a)/a

(1)則經(jīng)光柵抽樣后所得負片的密度分布D(x,y)為

D(x,y)={[D10-Di(x,y)]-D0}rect(2x/a)*comb(x/a)/a+D0式中Do是底片的灰霧密度，D10是可通過改變曝光條件來控制的常數(shù)，是底片的反差系數(shù)。這樣得到的一張矩形級數(shù)光柵其底片光密度可以簡記為

D(x,y)=D0Ts為0處D(x,y)=D10-Di(x,y)

Ts為1處(2)3/14/202359光學(xué)信息處理(2)漂白處理將經(jīng)抽樣所得到的負片進行漂白處理，并適當(dāng)控制漂白工藝，可以得到近似滿足光程差正比于底片密度的效果，即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復(fù)振幅透過率

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

最后得到編碼的相位光柵的振幅透過率為(3)3/14/202361光學(xué)信息處理(3)白光信息處理系統(tǒng)中的濾波解調(diào)將編碼相位光柵放在白光信息處理系統(tǒng)的輸入平面上，設(shè)入射單色光強度為A()，則頻譜平面(x2,y2)上的復(fù)振幅為(4)式中x2=f,x2=f,

將(3)式代人(4)式，計算并化簡后得(5)3/14/202362光學(xué)信息處理當(dāng)m=0時，有其對應(yīng)的強度只與相位差和有關(guān)，即如濾波器分別讓零級頻譜和m級頻譜通過，則在輸出平面上的復(fù)振幅分布為當(dāng)m≠0時，有(6)3/14/202363光學(xué)信息處理若用d表示與相位差相對應(yīng)的光程差，則(6)式可改寫為該式表明對于每一個衍射級次，輸出圖像的強度隨波長和光程差而變化。圖(a)和(b)分別為零級和1級的輸出強度隨光程差d而變化的曲線，其中A()=1，作為參變量。3/14/202364光學(xué)信息處理

假如用圖中畫出的紅、綠、藍三種色光r，g，b

照相，則強度輸出是三種色光輸出的非相干疊加:I(d)=I(d,r)+I(d,g)+I(d,b)

即得到隨d而變化的彩色輸出。當(dāng)采用白光光源時，各色光的非相干疊加變成下列積分：

I(d)=I(d,)d仍然是隨d而變化的彩色輸出。由于在編碼和漂白處理中，已使光程差隨輸入密度而改變，因此得到按輸入圖像密度變化的假彩色編碼。這種編碼方法輸出圖像的色度豐富，飽和度也很好，在低衍射級次(包括零極)輸出的情況下，也能得到彩色化效果很好的輸出圖像，因此光強度利用率高，圖像亮度好。3/14/202365光學(xué)信息處理3.14彩色膠片資料存貯

彩色膠片的資料存貯是膠片工業(yè)中長期沒有得到解決的問題，其主要原因是彩色膠片所用的染料不穩(wěn)定而造成逐漸褪色。雖然有一些技術(shù)可以用來保存彩色像，但都存在著不少問題。最常用的方法是用三個基色濾色片將彩色像保持在三個不同的黑白底片上，然后用三個基色幻燈機來合成彩色像。在這過程中三個像要非常精確地重合，此外這個方法還有兩個主要缺點：一是每個彩色片需用的存貯體積是黑白片的三倍；二是彩色像合成系統(tǒng)，非常昂貴和復(fù)雜。3/14/202366光學(xué)信息處理回顧1906年依伏斯首先報導(dǎo)了從黑白透明片恢復(fù)出彩色像的技術(shù)。他使用了一個用衍射現(xiàn)象來產(chǎn)生彩色像的幻燈器，其中使用了不同頻率或不同方位的光柵。1969年繆勒用了相似的方法即用三色光柵來進行彩色像的編碼。在解碼時，用三個準(zhǔn)單色光源來恢復(fù)彩色像。馬可夫斯基(1972年)、克羅松和基納雷(1978年)、楊振寰和母國光等人(1980年)的工作。3/14/202367光學(xué)信息處理彩色膠片記錄在黑白透明片上的技術(shù)記錄過程是將彩色膠片通過倫奇光柵依次用三次接觸法記錄在黑白底片上，每次記錄時通過不同的基色濾色片(即紅、藍、綠三色)，而對應(yīng)的倫奇光柵的方位角分別為0o位置(紅色濾色片)、60o位置(藍色濾色片)、120o位置(綠色濾色片)。3/14/202368光學(xué)信息處理如果這三次曝光很準(zhǔn)確地記錄在照相底片上，這就是一張空間編碼的黑白透明片，相應(yīng)的強度分布可用下式表示

po是光柵的空間頻率，Tn表示負像的強度透射率，而Tr、Tb、及Tg分別是在底片上紅、藍、綠三色的曝光量，(x，y)、(x’，y’)及(x”，y”)分別對應(yīng)于光柵在0o、60o及120o位置時的坐標(biāo)系統(tǒng)，n1是底片的伽瑪值，K1為常數(shù)。3/14/202369光學(xué)信息處理

可利用接觸曝光來得到適合于處理的正像透明片。多重光柵空間調(diào)制的正像透明片的強度透過率為

其中K2為比例常數(shù)，Tp為正片密度透過率，n2為第二張透明片的伽瑪值。為了得到像振幅的線性透過率，我們?nèi)1

n2=2，這樣振幅透過率為3/14/202370光學(xué)信息處理為了從所得到的黑白編碼透明片再現(xiàn)逼真的彩色像，我們采用白光處理技術(shù)。如圖所示，將該透明片放置在白光處理器的輸入平面P1處，于是在空間頻譜面P2上對波長為的光強分布為3/14/202371光學(xué)信息處理

為簡單起見，略去上式中的常數(shù)因子，則

其中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