光學(xué)鏡頭成像質(zhì)量研究

席瑞駿，肖雙江，楊 慧，劉 華

(東北師范大學(xué) 物理學(xué)院，吉林 長(zhǎng)春 130024)

攝影愛(ài)好者在測(cè)試鏡頭性能時(shí)，主要參考對(duì)象有：OECF光電轉(zhuǎn)換函數(shù)、信噪比、色彩準(zhǔn)確性(色差)、分辨力(空間頻率響應(yīng))、均勻性(亮度、色彩)等[1]，而MTF曲線或者彌散斑圖像是2種較好的綜合考量標(biāo)準(zhǔn). 目前，調(diào)制傳遞函數(shù)(MTF)可以通過(guò)點(diǎn)源法、矩形脈沖法、刃邊法、矩形靶標(biāo)法、高分辨率圖像法[2]、頻譜比較法[3]、對(duì)彌散斑進(jìn)行傅里葉變換等方法得到. 本文采用頻譜對(duì)比法，即用不同頻率的矩形光柵測(cè)量成像系統(tǒng)的MTF和觀察成像系統(tǒng)的彌散斑來(lái)評(píng)價(jià)鏡頭的成像質(zhì)量，對(duì)比2種方法所得結(jié)論的差異，與ZEMAX中理論圖像相比較，直觀展示彗差和球差相關(guān)性質(zhì)，并分析實(shí)驗(yàn)誤差. 本文以2個(gè)透鏡為例，透鏡1是主要像差為彗差的透鏡，透鏡2是主要像差為球差的透鏡. 這2個(gè)透鏡的其他像差對(duì)于成像的影響相對(duì)較小，在本實(shí)驗(yàn)中可忽略不計(jì).

1 MTF和彌散斑

1.1 MTF

MTF成像曲線圖是由鏡頭的生產(chǎn)廠家在極為客觀嚴(yán)謹(jǐn)?shù)臏y(cè)試環(huán)境下測(cè)量的，是鏡頭成像品質(zhì)權(quán)威且客觀的技術(shù)參考依據(jù). 通過(guò)MTF可以分析鏡頭的分辨率和對(duì)比度，MTF成像曲線下的面積越大說(shuō)明鏡頭的綜合性能越好，并且還可以反映出成像系統(tǒng)在各個(gè)頻率段的成像素質(zhì).

光學(xué)系統(tǒng)的調(diào)制傳遞函數(shù)表示為給定空間頻率下像和物調(diào)制度之比[4]：

(1)

其中物調(diào)制度M定義為

(2)

1.2 彌散斑

點(diǎn)光源(即星點(diǎn))經(jīng)過(guò)光學(xué)系統(tǒng)后在像面前后不同截面上所成的衍射像的光強(qiáng)分布即為彌散斑[5]. 由于點(diǎn)擴(kuò)散函數(shù)(PSF)是準(zhǔn)確描述彌散斑圖像信息的函數(shù)，而MTF又可以通過(guò)對(duì)PSF進(jìn)行傅里葉變換得到，所以通過(guò)觀察彌散斑，可以基本判斷出成像系統(tǒng)的像差情況，并且判斷出MTF曲線的基本趨勢(shì).

2 實(shí)驗(yàn)方法

2.1 模擬實(shí)驗(yàn)

利用ZEMAX模擬像差分別只有彗差、球差的單透鏡組成的成像系統(tǒng)1和成像系統(tǒng)2，直接得出成像系統(tǒng)的彌散斑圖像和MTF曲線. 成像系統(tǒng)1彌散斑圖像如圖1所示，MTF曲線如圖2所示，成像系統(tǒng)2彌散斑圖像如圖3所示， MTF曲線如圖4所示.

圖1 成像系統(tǒng)1彌散斑圖像

圖2 成像系統(tǒng)1 MTF曲線

圖3 成像系統(tǒng)2彌散斑圖像

圖4 成像系統(tǒng)2 MTF曲線

曲線中T代表子午方向(切向)的MTF曲線，S代表弧矢方向(徑向)的MTF曲線. 通過(guò)在ZEMAX中改變像差的大小和透鏡的孔徑，可以分析出它們對(duì)彌散斑圖像和MTF曲線的影響.

2.2 驗(yàn)證實(shí)驗(yàn)

2.2.1 實(shí)驗(yàn)過(guò)程

1.4.3 安全性 用藥前后檢測(cè)患者的血肌酐（Scr）、尿素氮（BUN）、胱抑素C水平。上述指標(biāo)均采用定量免疫分析法以TDx-FLx型熒光偏振免疫分析儀（美國(guó)Abbott公司）檢測(cè)，試劑盒均由美國(guó)Abbott公司提供。當(dāng)患者Scr水平較用藥前基線水平連續(xù)兩次增長(zhǎng)5 mg/L或增長(zhǎng)率≥50%時(shí)，則判定其腎功能發(fā)生異常[9]。

通過(guò)具體操作實(shí)驗(yàn)，驗(yàn)證從ZEMAX模擬軟件中得出的結(jié)論是否與實(shí)驗(yàn)吻合. 實(shí)驗(yàn)主要裝置有：白色LED光源(P>1 W)，準(zhǔn)直管組件(f=400 mm，通光口徑Φ50 mm，光潔度Ⅳ級(jí)透鏡)，透鏡1(Φ23 mm，f=90 mm，K9材料，光潔度Ⅳ級(jí))，透鏡2(Φ23 mm，f=100 mm，K9材料，光潔度Ⅳ級(jí)，光闌可調(diào))，國(guó)標(biāo)分辨力板A1，CMOS相機(jī)(1.3×106Pixel，分辨率1 280×1 024，黑白逐行掃描，像素大小5.2 μm×5.2 μm)和精密光學(xué)導(dǎo)軌. 以上實(shí)驗(yàn)儀器均由北京杏林睿光科技有限公司生產(chǎn).

首先拍攝彌散斑圖像. 在準(zhǔn)直管中，放有只有中心點(diǎn)透光的玻璃板，光線經(jīng)過(guò)透鏡1和透鏡2在CMOS上成像. 通過(guò)觀察CMOS拍攝的圖像可以明顯觀察到彗差和球差的彌散斑圖像. 改變透鏡1的光軸與成像系統(tǒng)光軸之間的水平夾角和透鏡2的光闌大小，觀察彗差和球差的彌散斑發(fā)生的變化.

透鏡1彌散斑圖像如圖5所示，透鏡2彌散斑圖像如圖6所示.

圖5 透鏡1彌散斑圖像

圖6 透鏡2彌散斑圖像

測(cè)量MTF曲線的實(shí)驗(yàn)裝置如圖7所示，從左到右依次為白色LED光源、國(guó)標(biāo)分辨力板A1、透鏡和CMOS相機(jī). 通過(guò)調(diào)節(jié)透鏡對(duì)國(guó)標(biāo)分辨力板A1上不同空間頻率的矩形光柵進(jìn)行成像，使用相機(jī)采集程序保存圖像. 在使用透鏡1時(shí)，改變透鏡1的光軸與成像系統(tǒng)光軸之間的水平夾角，選取3個(gè)不同角度得到透鏡1的MTF曲線. 在使用透鏡2時(shí)，改變透鏡的光闌得到3個(gè)不同通光口徑的MTF曲線.

圖7 實(shí)驗(yàn)裝置

2.2.2 圖像處理

對(duì)不同空間頻率的矩形光柵采集的圖像，編寫(xiě)Matlab程序讀取圖像灰度值的極大值和極小值，利用圖像處理軟件ImageJ處理圖像. 對(duì)于同一頻率不同方向的4個(gè)矩形光柵，由于彗差在子午方向和弧矢方向產(chǎn)生的畸變不同，球差在各個(gè)方向產(chǎn)生的畸變相同，因此對(duì)透鏡1的圖像取左上、右下光柵圖像進(jìn)行處理；透鏡2圖像只取左上光柵圖像進(jìn)行處理. 分別選取每張圖像中間的3個(gè)區(qū)域，得到3組灰度最大值及最小值，取其平均值，代入式(1)和式(2)得到該空間頻率矩形光柵的像對(duì)比度和MTF值. 以矩形光柵的頻率為橫坐標(biāo)，MTF值為縱坐標(biāo)，畫(huà)出MTF曲線，如圖8所示和圖9所示.

圖8 透鏡1的MTF曲線

圖9 透鏡2的MTF曲線

3 結(jié)果分析

3.1 彌散斑分析

3.1.1 對(duì)透鏡1的分析

利用ZEMAX得出不同彗差大小的成像系統(tǒng)的彌散斑圖像，發(fā)現(xiàn)彗差越大，彌散斑圖像的面積越大. 在實(shí)驗(yàn)中，拍攝不同偏轉(zhuǎn)角度對(duì)應(yīng)的彌散斑圖像，發(fā)現(xiàn)偏轉(zhuǎn)角度越大，面積越大. 由于彗差是軸外像差，偏轉(zhuǎn)角度越大則彗差越大，與模擬結(jié)果一致.

3.1.2 對(duì)透鏡2的分析

利用ZEMAX模擬透鏡的彌散斑，得出不同球差對(duì)應(yīng)的成像系統(tǒng)的彌散斑圖像，發(fā)現(xiàn)球差的大小與彌散斑圖像的面積不是簡(jiǎn)單的正比關(guān)系，與實(shí)驗(yàn)得到的結(jié)論吻合.

3.2 MTF分析

實(shí)驗(yàn)得到的MTF曲線都低于模擬得到的MTF曲線，這是多方面因素導(dǎo)致的. 首先，國(guó)標(biāo)分辨力板A1中的矩形光柵的空間頻率范圍是從3.13 mm-1開(kāi)始，而在ZEMAX模擬中是從0開(kāi)始；其次，ZEMAX模擬的是僅有彗差或者球差這種像差的理想成像系統(tǒng)，并且沒(méi)有其他任何誤差，而在實(shí)驗(yàn)中，透鏡1或2不只存在一種像差，只是以某種像差為主. 數(shù)據(jù)處理的過(guò)程中，默認(rèn)物的調(diào)制度M為1，而實(shí)際上國(guó)標(biāo)分辨力板A1的調(diào)制度并不為1. 其他因素將在誤差分析中討論. 因此需要比較的是理論與實(shí)際MTF曲線在趨勢(shì)上的一致性.

3.2.1 對(duì)透鏡1的分析

首先分析當(dāng)透鏡無(wú)偏轉(zhuǎn)角度時(shí)， MTF曲線與模擬MTF圖在趨勢(shì)上一致，模擬曲線中T代表子午方向(切向)對(duì)應(yīng)右下矩形光柵，S代表弧矢方向(徑向)對(duì)應(yīng)左上矩形光柵. 在實(shí)驗(yàn)中，用分辨力板左上和右下方向的光柵測(cè)出的MTF曲線彼此的趨勢(shì)雖一致，量值并不相同，當(dāng)頻率在3.13～7.87 mm-1時(shí)，弧矢方向的MTF值大于子午方向的MTF值，這與ZEMAX模擬結(jié)果一致，驗(yàn)證了彗差透鏡在不同方向的彗差大小不相同.

接下來(lái)分析不同偏轉(zhuǎn)角度對(duì)MTF的影響，并與ZEMAX的模擬結(jié)果相對(duì)比. 在ZEMAX中，對(duì)成像系統(tǒng)1中的透鏡分別進(jìn)行0°，10°和20°的偏轉(zhuǎn)，利用ZEMAX的模擬成像功能，發(fā)現(xiàn)隨著偏轉(zhuǎn)角度的增大，同一光柵頻率下的對(duì)比度降低. 在實(shí)驗(yàn)中，觀察圖9發(fā)現(xiàn)，當(dāng)光柵頻率為3.13～7.87 mm-1時(shí)，彗差透鏡偏轉(zhuǎn)角度測(cè)出的MTF曲線低于無(wú)偏轉(zhuǎn)角度的MTF曲線，與模擬結(jié)果一致，當(dāng)光柵頻率大于7.87 mm-1時(shí)，由于CMOS相機(jī)的分辨率較低以及MTF值小導(dǎo)致相對(duì)誤差變大等因素，有無(wú)偏轉(zhuǎn)角度對(duì)MTF值影響不大.

3.2.2 對(duì)透鏡2的分析

觀察圖10發(fā)現(xiàn)，透鏡的MTF曲線，與模擬MTF曲線在趨勢(shì)上一致；利用ZEMAX改變同一成像系統(tǒng)中透鏡前的通光口徑，發(fā)現(xiàn)隨著口徑的增大，MTF曲線下降. 在實(shí)驗(yàn)中，當(dāng)增大透鏡前通光口徑時(shí)，對(duì)應(yīng)的MTF曲線下降，成像質(zhì)量降低，與模擬結(jié)果一致，驗(yàn)證了球差的性質(zhì).

3.3 誤差分析

3.3.1 操作誤差

在實(shí)驗(yàn)測(cè)量前，僅借助輔助棒難以將各器件嚴(yán)格調(diào)節(jié)至共軸等高，這會(huì)使各元件光軸與成像系統(tǒng)光軸不在同一條直線上，從而降低成像質(zhì)量.

3.3.2 儀器誤差

在實(shí)驗(yàn)中，光源是白色LED光源，光強(qiáng)可能有不均勻性，使經(jīng)過(guò)矩形光柵前的光強(qiáng)分布不同，從而使得圖像處理過(guò)程中獲取像素的灰度值發(fā)生偏差[4].

在使用CMOS拍攝時(shí)，一方面，不可避免地存在許多噪點(diǎn)，因此對(duì)圖像進(jìn)行處理時(shí)，獲取了不同區(qū)域灰度值的最大值及最小值并取平均值，盡量減少噪點(diǎn)帶來(lái)的影響[6-7]；另一方面，在CMOS中像素構(gòu)造方式是矩陣排列，因此在對(duì)有周期性圖案成像時(shí)就有可能產(chǎn)生疊紋效應(yīng). 疊紋效應(yīng)是由于2組周期性的圖案在分布方向上相差一個(gè)角度造成的.

在本實(shí)驗(yàn)中，成像的光線在CMOS感光區(qū)域與透鏡之間有來(lái)回反射的光線，形成二次成像，并與第一次成像相互疊加，影響成像質(zhì)量[8]. 當(dāng)圖像對(duì)比度較低時(shí)，由于數(shù)據(jù)處理方式的限制，圖像的噪點(diǎn)會(huì)大幅度影響計(jì)算出的對(duì)比度，因此圖9和圖10中，當(dāng)對(duì)比度較低時(shí)，圖像出現(xiàn)一定范圍的波動(dòng).

4 結(jié)束語(yǔ)

通過(guò)實(shí)驗(yàn)繪制MTF曲線，并觀察彌散斑，理論上用ZEMAX軟件模擬MTF曲線這3種方法具體展示了球差和彗差在不同通光口徑、不同偏轉(zhuǎn)角度對(duì)成像質(zhì)量的影響. 其中，MTF曲線可以直觀地反映出成像系統(tǒng)的質(zhì)量高低，并且成像系統(tǒng)的MTF曲線也體現(xiàn)了系統(tǒng)的對(duì)比度和分辨率的傳遞能力，而利用彌散斑圖像的方法雖然直觀,但帶有主觀性，不能作定量評(píng)價(jià). MTF曲線的測(cè)試對(duì)測(cè)試設(shè)備和測(cè)試方法的要求都比較高，本文結(jié)合實(shí)驗(yàn)室現(xiàn)有的條件，選擇物像頻譜對(duì)比法測(cè)透鏡的MTF曲線，該方法操作簡(jiǎn)單，成本低，可作為基礎(chǔ)實(shí)驗(yàn)來(lái)開(kāi)設(shè)，加深學(xué)生對(duì)于光學(xué)像差的理解. 但該方法也有值得改進(jìn)之處，如得到的MTF曲線不夠光滑，可通過(guò)合適的數(shù)據(jù)插值擬合算法(牛頓插值法)進(jìn)行處理[4]，對(duì)CMOS獲取的圖像進(jìn)行空間域?yàn)V波[6]且可深入探討球差的大小、離焦程度、波長(zhǎng)等因素對(duì)彌散斑的影響[5].