趙忠超 楊旭鋒 許天旭 何九如 弓巧俠 杜艷麗董林 袁斌 馬鳳英

(鄭州大學(xué)物理工程學(xué)院,鄭州 450001)

1 引 言

全息術(shù)[1]是利用光的干涉和衍射原理記錄并再現(xiàn)物體真實三維空間信息的技術(shù).傳統(tǒng)全息術(shù)采用相干光源照明,不僅對系統(tǒng)的穩(wěn)定性要求極高,而且還會引入嚴(yán)重的相干噪聲,限制了其在眾多領(lǐng)域中的應(yīng)用.自1961年Mertz等[2]將全息術(shù)推廣到非相干領(lǐng)域后,基于各種不同分光技術(shù)的非相干全息術(shù)得到了飛速發(fā)展[3?8].此外,隨著計算機處理速度的提高和電荷耦合器(CCD)的問世,Goodman等[9]首次提出了利用計算機再現(xiàn)全息圖的思想.1971年,Huang[10]提出了數(shù)字全息的概念.1994 年,Schnar和Jüptner[11,12]利用CCD記錄并實現(xiàn)了全息圖再現(xiàn)的完全數(shù)字化.自此,數(shù)字全息術(shù)[13?17]受到了越來越多的關(guān)注.數(shù)字全息技術(shù)省卻了干板復(fù)雜的化學(xué)濕處理過程,在純相位物體定量成像方面具有顯著的優(yōu)勢.

近年來,國內(nèi)外眾多研究小組對非相干光源照明下的數(shù)字全息進(jìn)行了研究.Kim[4,18]利用邁克耳孫干涉裝置在自然光下實現(xiàn)了彩色物體的遠(yuǎn)距離全息記錄和再現(xiàn),并應(yīng)用于自適應(yīng)光學(xué)領(lǐng)域.Wan等[3]對非相干數(shù)字全息成像系統(tǒng)的研究進(jìn)展和成像特性分別進(jìn)行了介紹,并基于三角干涉光路提出了離軸非相干傅里葉彩色數(shù)字全息術(shù)(incoherent Fourier triangular color digital holography,IFTCH).Rosen和Brooker[19]于2007年提出了一種基于SLM分光的FINCH技術(shù),該技術(shù)具有光路簡單、既無時間掃描也無空間掃描、獲得全息圖速度快、分辨率高,且易于與現(xiàn)有成熟光學(xué)系統(tǒng)匹配的特點,因此一經(jīng)提出便被應(yīng)用于熒光顯微[20]、白光反射全息成像[21]和彩色3D成像[22]等領(lǐng)域.2012年,Bouchal和Bouchal[23]采用FINCH系統(tǒng)實現(xiàn)了圖像邊緣信息的提取,具體方法為:在SLM上加載螺旋相位調(diào)制的雙透鏡相位掩模,記錄物體的全息圖并數(shù)值再現(xiàn)后得到邊緣提取的圖像信息.但文獻(xiàn)[23]中并未給出系統(tǒng)螺旋點擴散函數(shù)的數(shù)學(xué)表達(dá)式,也未驗證系統(tǒng)對相位物體的邊緣提取與識別能力.

本文從波動光學(xué)角度出發(fā),給出了基于純相位SLM的FINCH系統(tǒng)螺旋點擴散函數(shù)表達(dá)式;實驗研究了系統(tǒng)的脈沖響應(yīng),得到了與點擴散函數(shù)模擬結(jié)果一致的結(jié)論,證明了點擴散函數(shù)表達(dá)式的正確性;采用該系統(tǒng)對分辨率板成像,在幾乎不損失分辨率的情況下實現(xiàn)了圖像邊緣提取;對非染色的洋蔥細(xì)胞進(jìn)行全息拍攝,驗證了系統(tǒng)對相位物體的邊緣提取與識別能力.

2 理論分析

圖1為FINCH記錄光路原理圖,zs是物平面與透鏡L之間的距離,d是空間光調(diào)制器到透鏡L的距離,zh是SLM到CCD的距離.

圖1 FINCH系統(tǒng)記錄光路原理圖Fig.1.Schematic diagram of FINCH system.

將SLM所有像素隨機分成兩等份,分別加載兩個透鏡(焦距分別為fd1和fd2)的相位.為實現(xiàn)螺旋相位調(diào)制,在其中一個透鏡上疊加螺旋相位mφ.SLM的振幅反射系數(shù)為

式中λ為波長,θ為相移常數(shù),m為拓?fù)浜蓴?shù),φ=arctan(y/x)為SLM平面的方位角,B和B′均為常數(shù).

假設(shè)點源位于透鏡L(焦距f0)焦平面附近的光軸上,坐標(biāo)為(0,0,?zs).傍軸條件下,其發(fā)出的球面波在透鏡L前表面的復(fù)振幅可近似表示為

式中A為球面半徑為1處的振幅,利用數(shù)學(xué)關(guān)系可得

令σ2=2z2s,經(jīng)過透鏡L后復(fù)振幅為

式中

令

再經(jīng)SLM反射后變成兩束自相干的光波,傳播距離zh后在CCD表面的復(fù)振幅分別為

和公式

兩束相干光束在CCD表面發(fā)生干涉,強度分布即點擴散函數(shù)(PSF)為

式中δ為螺旋光束傳輸中產(chǎn)生的相位延遲,γ為CCD平面的方位角,C和C1為關(guān)于(x2+y2)的函數(shù),

c.c.是等號右邊第二項的復(fù)共軛.由點擴散函數(shù)表達(dá)式(9)可知CCD平面上記錄的點源全息圖具有螺旋結(jié)構(gòu),對應(yīng)拓?fù)浜蔀閙.具有不同相位常數(shù)的點源全息經(jīng)線性處理消除零級像和孿生像后得到的復(fù)值全息具有螺旋相位,螺旋相位中心存在一個相位奇點,此處相位不確定,振幅為零,因此相位奇點處光強為零(螺旋相位波前在該處干涉相消),即再現(xiàn)后得到中心光場被抑制的空心點像.

設(shè)非相干光照射下物體表面強度為g(xs,ys,?zs),其在CCD上記錄的干涉圖案為三維物體表面所有點源全息圖的非相干疊加,

其中MT為橫向放大率.

相位常數(shù)分別為θ1,θ2,θ3時,CCD記錄的全息圖分別為H1(x,y),H2(x,y),H3(x,y).再現(xiàn)時,將3次記錄的全息圖進(jìn)行線性疊加得到的復(fù)全息HF(x,y)為

利用角譜衍射算法[24]模擬衍射過程實現(xiàn)數(shù)值再現(xiàn),得到再現(xiàn)像

式中k=2π/λ,F和F?1分別表示正逆傅里葉變換,fx和fy表示與x,y對應(yīng)的頻域坐標(biāo).

3 模擬與實驗分析

3.1 點擴散函數(shù)

采用(9)式模擬FINCH系統(tǒng)的螺旋點擴散函數(shù),相位常數(shù)為0°、拓?fù)浜蔀?時的仿真結(jié)果如圖2(a)所示,對應(yīng)的模擬再現(xiàn)像如圖2(c)所示.為驗證點擴散函數(shù)表達(dá)式的正確性,以直徑20μm的小孔作為點源進(jìn)行實驗,圖2(b)為相移常數(shù)為0°時CCD記錄的小孔全息圖(拓?fù)浜蒻=1),圖2(d)為對應(yīng)的再現(xiàn)像.比較圖2(a)和圖2(b)及圖2(c)和圖2(d),并考慮到CCD響應(yīng)靈敏度等因素的限制,模擬結(jié)果與實驗結(jié)果基本一致.

圖2 (a)θ為0?時模擬點源全息圖;(b)θ為0?時實驗點源全息圖;(c)模擬點源再現(xiàn)像;(d)實驗點源再現(xiàn)像Fig.2.(a)The simulated holograms with the phase factor equals zero;(b)the experimental holograms with the phase factor equals zero;(c)simulated reconstructed image;(d)experimental reconstructed image.

3.2 振幅物體的全息記錄與再現(xiàn)

搭建如圖3所示的反射式非相干數(shù)字全息記錄系統(tǒng),測試物體為USAF1951分辨率板.反射式空間光調(diào)制器像素大小為8μm,像素數(shù)為1920×1080,CCD像素大小為4.54μm,像素數(shù)為2750×2200,實驗中只使用了1024×1024個像素點.其中光源為白光連續(xù)光譜光源,BS為分束器,準(zhǔn)直透鏡L2的焦距為250 mm,偏振片P的偏振方向與SLM的液晶長軸方向一致,窄帶濾光片BF的帶寬為20 nm,中心波長為632.8 nm.分辨率板到透鏡L2的距離為330 mm,透鏡L2到SLM的距離為150 mm,SLM和CCD之間的距離為zh,為獲得最佳干涉效果,其具體數(shù)值需與SLM上加載的透鏡焦距大小相匹配.實驗中取fd1=245 mm,fd2=255 mm,zh=250 mm,拓?fù)浜蒻=1.

圖3 非相干數(shù)字全息實驗光路Fig.3.Experimental set-up of incoherent digital holography.

為驗證系統(tǒng)對振幅物體的邊緣提取與識別特性,需要在SLM上分別加載雙透鏡相位掩模和螺旋相位調(diào)制的雙透鏡掩模,并在兩種相位調(diào)制模式下,改變相位常數(shù)θ分別為0°,120°和240°時各拍攝三張全息圖后進(jìn)行數(shù)字重建.圖4(a)和圖4(c)分別為θ=0°時雙透鏡相位掩模及對應(yīng)模式下記錄的全息圖,圖4(b)和圖4(d)分別為θ=0°時螺旋相位調(diào)制雙透鏡相位掩模及對應(yīng)模式下記錄的全息圖.在計算機中對所記錄不同相位常數(shù)的全息圖進(jìn)行處理及再現(xiàn),獲得的再現(xiàn)像如圖5所示.其中圖5(a)為雙透鏡相位掩模調(diào)制下分辨率板的再現(xiàn)像,圖5(b)為螺旋相位掩模調(diào)制下分辨率板的再現(xiàn)像,圖5(c)中的實線和虛線分別對應(yīng)圖5(a)和圖 5(b)中豎線部分的歸一化強度分布曲線.對比圖5(a)和圖5(b)可知,當(dāng)SLM加載螺旋相位調(diào)制的雙透鏡掩模(m=1)時,可實現(xiàn)圖像各向同性的邊緣提取,再現(xiàn)像背景噪聲小,可在不犧牲分辨率的情況下實現(xiàn)圖像的邊緣提取與識別.

圖4 θ為0?時(a)雙透鏡相位掩模,(b)螺旋相位調(diào)制下的雙透鏡相位掩模,(c)雙透鏡模式下記錄的全息圖,(d)螺旋相位調(diào)制模式下記錄的全息圖Fig.4.(a)Double lens phase mask with the phase factor equals zero,(b)holograms captured under double lens masks with the phase factor equals zero,(c)spiral phase mask with the phase factor equals zero,(d)holograms captured under spiral phase mask with the phase factor equals zero.

圖5 分辨率板再現(xiàn)像 (a)雙透鏡相位掩模調(diào)制;(b)螺旋相位掩模調(diào)制;(c)(a)和(b)中豎線部分的歸一化強度分布Fig.5. The reconstructed images of resolution target:(a)Double lens phase mask modulation;(b)spiral phase mask modulation;(c)normalized intensity curves of the vertical line in(a)and(b).

3.3 純相位物體的全息記錄與再現(xiàn)

為驗證螺旋相位調(diào)制的FINCH系統(tǒng)對純相位物體的成像特性,采用非染色洋蔥表皮細(xì)胞作為記錄物體.為獲得實驗結(jié)果,需在圖3所示實驗裝置中記錄物體前加一個20倍的顯微鏡頭,并采用透射式照明光路,其他元件及參數(shù)均不變.圖6(a)和圖6(b)分別為相位常數(shù)為0°時雙透鏡相位掩模和螺旋調(diào)制的雙透鏡相位掩模兩種模式下記錄的全息圖,圖6(c)和圖6(d)分別為圖6(a)和圖6(b)兩種模式下對應(yīng)的再現(xiàn)圖.可以看出,螺旋相位調(diào)制模式得到的再現(xiàn)像有效地抑制了細(xì)胞內(nèi)部的光場,具有更清晰的邊緣信息.該特性在非標(biāo)記、動態(tài)實時細(xì)胞形態(tài)分析方面有潛在的應(yīng)用價值.

圖6 (a)θ為0?時雙透鏡相位掩模下記錄的全息圖;(b)θ為0?時螺旋相位掩模下記錄的全息圖;(c)雙透鏡相位掩模下的再現(xiàn)像;(d)螺旋相位掩模下的再現(xiàn)像Fig.6. Holograms recorded under(a)double lens phase mask and(b)spiral phase mask with the phase factor equals zero;reconstructed images under(c)double lens phase mask and(d)spiral phase mask.

4 結(jié) 論

從波動光學(xué)角度出發(fā),分析了螺旋相位調(diào)制下FINCH系統(tǒng)的記錄與再現(xiàn)原理,給出了系統(tǒng)螺旋點擴散函數(shù)的表達(dá)式,模擬結(jié)果與實驗結(jié)果相吻合.采用螺旋相位調(diào)制的FINCH系統(tǒng)對分辨率板成像,在幾乎不損失分辨率的情況下實現(xiàn)了圖像邊緣提取的效果.對非染色洋蔥細(xì)胞成像,獲得了相比于非螺旋調(diào)制模式更清晰的細(xì)胞邊緣信息.該系統(tǒng)對光路的穩(wěn)定性和光源的相干性等方面要求不高,可以在非相干光照明下實現(xiàn)對振幅或相位物體的邊緣提取,在活細(xì)胞的分裂、形變監(jiān)測及醫(yī)學(xué)檢驗等領(lǐng)域有重要應(yīng)用前景.

[1]Gabor D 1948Nature161 777

[2]Mertz L,Young N O 1962Optical Instruments and TechniquesLondon,United Kingdom,July 11–14,1961 p305

[3]Wan Y,Wang D,Man T 2014Opt.Express22 8565

[4]Kim M K 2013Opt.Express21 9636

[5]Naik D N,Pedrini G,Osten W 2013Opt.Express21 3990

[6]Rosen J,Brooker G 2007Opt.Express15 2244

[7]Cochran G 1966J.Opt.Soc.Am.56 1513

[8]Stroke G W,Restrick R C 1965Appl.Phys.Lett.7 229

[9]Goodman J W,Lawrence R W 1967Appl.Phys.Lett.11 77

[10]Huang T S 1971Proc.IEEE59 1335

[11]Schnars U,Jüptner W P O 2002Meas.Sci.Technol.13 R85

[12]Schnars U,Jüptner W 1994Appl.Opt.33 179

[13]Yang J,Wu X C,Wu Y C,Yao L C,Chen L H,Qiu K Z,Cen K F 2015Acta Phys.Sin.64 114209(in Chinese)[陽靜,吳學(xué)成,吳迎春,姚龍超,陳玲紅,邱坤贊,岑可法2015物理學(xué)報64 114209]

[14]Li J C,Peng Z J,Fu Y C 2011Chin.Phys.Lett.28 064201

[15]Zhang Q X,Lü X X,Yu Q T,Liu G Y 2009Chin.Phys.B18 2764

[16]Lu X W,Li J Z,Chen H Y 2010Chin.Phys.Lett.27 104209

[17]Li J C 2012Acta Phys.Sin.61 134203(in Chinese)[李俊昌2012物理學(xué)報61 134203]

[18]Kim M K 2012Opt.Lett.37 2694

[19]Rosen J,Brooker G 2007Opt.Lett.32 912

[20]Rosen J,Brooker G 2008Nat.Photon.2 190

[21]Rosen J,Brooker G 2012Adv.Opt.Techn.1 151

[22]Shi X,Zhu W F,Yuan B,Du Y L,Gong Q X,Guo M T,Liang E J,Ma F Y 2015Chin.J.Las.42 265(in Chinese)[石俠,朱五鳳,袁斌,杜艷麗,弓巧俠,郭茂田,梁二軍,馬鳳英2015中國激光42 265]

[23]Bouchal P,Bouchal Z 2012Opt.Lett.37 2949

[24]Li J C,Song Q X,Pascal P,Gui J B,Lou Y L 2014Chin.J.Las.41 81(in Chinese)[李俊昌,宋慶和,Picart Pascal,桂進(jìn)斌,樓宇麗2014中國激光41 81]