夏春秋，鐘　興，劉春雨，韓培仙1，，金　光

(1.中國科學院大學，北京 100049；2.中國科學院 長春光學精密機械與物理研究所，吉林 長春 130033；3.長光衛(wèi)星技術有限公司，吉林 長春 130031)

?

高分辨率4f傳像系統(tǒng)分辨率的影響因素分析

夏春秋1，2*，鐘興2，3，劉春雨2，韓培仙1，2，金光2

(1.中國科學院大學，北京 100049；2.中國科學院 長春光學精密機械與物理研究所，吉林 長春 130033；3.長光衛(wèi)星技術有限公司，吉林 長春 130031)

研究了在實際應用中4f傳像系統(tǒng)分辨率的影響因素及其作用機制，用于指導系統(tǒng)的設計和裝調(diào)。首先，在理論上分析了由裝調(diào)誤差引起的像差對系統(tǒng)光學傳遞函數(shù)(MTF)產(chǎn)生的影響，并建立了描述二者關系的數(shù)學模型?；谠撃Ｐ脱芯苛? μm分辨率4f光學系統(tǒng)的裝調(diào)誤差，仿真計算了系統(tǒng)含有透鏡偏心、傾斜、位移誤差對MTF的影響，得到了分辨率隨誤差的變化曲線。實驗結果表明，傾斜在15″以內(nèi)，偏心誤差在0.02 mm，位移在0.01 mm范圍內(nèi)，裝調(diào)誤差均會導致分辨率下降，雖然傾斜在極小范圍內(nèi)微調(diào)時分辨率會提高，但總的變化趨勢不變。繪制的裝調(diào)誤差與分辨率的關系圖，可為高分辨率4f系統(tǒng)的誤差預估提供參考。

信息光學；4f傳像系統(tǒng)；裝調(diào)誤差；數(shù)學模型；曲線分析；高分辨率

1　引　言

光學信息處理是以信息光學為基礎，用傅里葉分析的方法研究光學成像和光學變換的技術[1-2]。4f系統(tǒng)的相干光空間濾波是光學信息處理的重要實現(xiàn)系統(tǒng)之一。它以二維圖像作為媒介來實現(xiàn)圖像的識別、圖像的增強與恢復、圖像的傳輸與變換、功率譜分析以及全息術中的傅里葉全息存儲等，因此在圖像光學小波變換、數(shù)字全息及工件誤差檢測、粒子測量等領域有著廣泛的應用[3-5]。

根據(jù)所需要處理的圖形或者工作內(nèi)容不同，4f系統(tǒng)設計主要考慮待處理輸入面的尺寸和所需的最高空間頻率及鏡頭負擔等問題。目前，4f系統(tǒng)的分辨率為1 μm左右。

隨著光學信息處理器件的發(fā)展，4f系統(tǒng)越來越多地與數(shù)字計算機結合使用，其應用前景也更加廣闊，并且隨著圖像精度需求的提高，分辨率的要求也越來越高。分辨率性能的優(yōu)劣直接影響到光學信息處理結果的好壞，它作為成像質(zhì)量判斷的一個重要指標，在設計、裝調(diào)中具有指導意義。本文詳細討論了影響4f系統(tǒng)分辨率的主要因素，建立了影響高分辨4f傳像系統(tǒng)分辨率的主要因素的數(shù)學模型，基于該模型對1 μm分辨率的4f光學系統(tǒng)進行了全面的對比分析和評價，從而提高了實際裝調(diào)和圖像處理的效率。

2　4f傳像系統(tǒng)分辨率的影響因素分析

圖1所示為4f系統(tǒng)的結構示意圖。其中，S為照明點，光源L為準直透鏡， L1，L2分別為第一、第二傅里葉變換透鏡，P1為物面，P為頻譜面，P2為系統(tǒng)像面。由于4f系統(tǒng)是小像差光學系統(tǒng)，所以采用分辨率板和點列圖很難有效地衡量其分辨率。調(diào)制傳遞函數(shù)(Modulation Transfer Function，MTF)曲線表示不同頻率的正弦強度分布函數(shù)，即(物體)經(jīng)光學系統(tǒng)成像后其對比度的衰減程度。它既與光學系統(tǒng)像差有關，又與光學系統(tǒng)的衍射效果有關，能同時運用于小像差光學系統(tǒng)和大像差光學系統(tǒng)。因此這里采用MTF曲線作為4f系統(tǒng)分辨率的判斷依據(jù)。

圖1　4f系統(tǒng)結構示意圖

光學系統(tǒng)的誤差來源主要有光學、機械兩方面。其中，光學誤差主要包括加工工藝引起的半徑厚度、偏心等偏差；光機加工、裝調(diào)中的誤差可以概括為元件的厚度偏差、位移、偏心和傾斜等。在實際系統(tǒng)中，透鏡系統(tǒng)的結構參數(shù)如曲率半徑、厚度等是確定的，各光學面的曲率半徑和面形、玻璃材料的質(zhì)量等級都能控制在所要求的公差之內(nèi)，即使存在面型誤差，也可以通過微調(diào)間距使系統(tǒng)滿足共軸公差要求。因此跟輔助裝調(diào)有關的變量，如光學面之間的空氣間隔、光學面的傾斜與偏心等成為引起系統(tǒng)像質(zhì)變化的主要因素。

實際實驗得到的與位移、偏心、傾斜等因素相關的結論如下：

(1)光路是否共軸。光通過透鏡不同的地方，因為透鏡的厚度不同，從而使得位相調(diào)制函數(shù)不同，而影響成像效果；

(2)光通過傅氏透鏡的中心后，能否將亮斑完全照射到成像面上，使得圖像成像完整；

(3)物面、像面的擺放位置是否嚴格遵守4f透鏡的要求。

在明確裝調(diào)誤差是4f系統(tǒng)分辨率的主要影響因素之后，通過建立裝調(diào)誤差和MTF的數(shù)學模型來探究4f分辨率影響因素及其影響機制，對高分辨率4f系統(tǒng)的誤差預估工作有指導意義[6]。

3　影響因素與MTF的數(shù)學關系模型

如圖2所示，當成像系統(tǒng)成像時，出瞳的理想面為So，像差W(x，y)的存在導致實際波面S偏離了球面。此時出瞳面上的透過率函數(shù)(即光瞳函數(shù))為：

P(x，y)=P(x，y)exp[jkW(x，y)].

(1)

根據(jù)式(1)，系統(tǒng)的通頻帶范圍仍由光瞳尺寸決定，截止頻率和像差無關。但是像差在通帶內(nèi)引入了與頻率有關的相位畸變，會影響像質(zhì)。而在實際光學系統(tǒng)的MTF曲線中，滿足某觀察對比度的截止頻率代表了系統(tǒng)的分辨能力。

圖2　波像差示意圖

在光學系統(tǒng)中，波像差W(x，y)可用Zernike多項式進行描述：

(2)

表1　Zernike多項式的表述

一般采用Zernike多項式系數(shù)Z5，Z6，Z7，Z8，Z9作為系統(tǒng)成像質(zhì)量的評價標準。當存在上述失調(diào)量時，光學系統(tǒng)會產(chǎn)生像差，用Zernike多項式來表示像差，那么它與失調(diào)量之間的函數(shù)關系為：

(3)

其中：Fj(j=1，2，3，…，n)為系統(tǒng)像差，xi(i=1，2，3，…，n)為各元件位置的結構參量，fj(j=1，2，3，…，n)代表像差與鏡面位置之間的函數(shù)關系[7-8]。

經(jīng)過泰勒展開變換，像差與位置結構參量變化之間的近似線性方程組為：

(4)

式(4)可化為失調(diào)量與像差之間的數(shù)學模型：

A·ΔX=ΔF，

(5)

式中：A為靈敏度矩陣；ΔX為系統(tǒng)中各片鏡面需要調(diào)整的變化量，即失調(diào)量。失調(diào)量通常包括沿垂直于光軸方向的移動量、轉動量及各鏡面之間的軸向間隔(TEDY、TETX、TETY)；ΔF為各校正對像的實測值與光學設計值之差，也就是系統(tǒng)成像質(zhì)量隨失調(diào)量的變化值。因此已知失調(diào)量引起的像差，可以采用波像差理論和MTF建立數(shù)學模型，MTF作為光學調(diào)制傳遞函數(shù)，可以用傅里葉級數(shù)或積分的形式表述，從而得到失調(diào)量對分辨率的影響。

根據(jù)4f系統(tǒng)傅里葉變換的規(guī)則，用瞳孔位置坐標相關的澤尼克多項式表達波像差[8-9]，即有：

(6)

a3(x2+y2)+a4(y2-x2)+a5(xy)+

a6[-2x+3x(x2+y2)]+

a7[-2y+3y(x2+y2)]+

a8[1-6(x2+y2)+6(x2+y2)]2+…，

(7)

式中：an為Zernike多項式的展開系數(shù)，Zn(ε，η)為Zernike多項式。各多項式的物理意義分別對應表1中的2～9項。由于在4f傳像系統(tǒng)的實際使用過程中通常有光闌限制，這種情況下的Zernike高階像差很小，對成像的影響也很小，因此主要討論該公式中Zernike的前8項像差，即對其初級像差進行討論和研究。

本文以x，y方向均存在傾斜像差的4f系統(tǒng)為例，用Zernike多項式將式(7)與裝調(diào)和加工引起的像差聯(lián)系起來，得到系統(tǒng)的出瞳函數(shù)為：

(9)

這時出瞳面上的透過率函數(shù)(即光瞳函數(shù))為：

Hc(fx，fy)=P(λdifx，λdify)=

P0(λdifx，λdify)exp(jkW(λdifx，λdify)).

(10)

其波相差為：

(11)

在理想情況下光瞳是圓形光斑，其表達式為：

(12)

其中r為光瞳半徑。則有：

P(ε，η)=P0(ε，η)·exp[jkW(x，y)]=

(13)

脈沖響應為光瞳函數(shù)的傅里葉變換：

H=F{Δφ}=F{P(ε，η)}=

(14)

在實際脈沖響應的計算過程中，相干脈沖響應不再單純是孔徑的夫瑯禾費衍射圖樣，必須考慮波像差的影響。如式(14)所示，雖然傾斜像差不改變截止頻率只產(chǎn)生像點位移，但更高級次的像差會在通帶內(nèi)引入與頻率有關的相位畸變，使像質(zhì)變壞。因此在分析像差對分辨率的影響時，以MTF函數(shù)作為分析依據(jù)具有對小像差系統(tǒng)敏感的優(yōu)點，同時更加客觀和可靠。MTF為脈沖響應模的歸一化形式，即：

(15)

將式(14)的脈沖響應代入式(15)，根據(jù)MTF曲線的值可以得到某一對比度下系統(tǒng)的空間頻率，即分辨率的倒數(shù)。

4f系統(tǒng)屬于小像差系統(tǒng)，因此Zemax中MTF曲線一般用FFT MTF表示。快速傅里葉變換(Fast Fourier Transform，F(xiàn)FT)是對連續(xù)變換的信息進行采樣然后再根據(jù)離散信號進行頻譜計算。在分析偏心、傾斜、位移誤差對實際裝調(diào)測量的光學系統(tǒng)分辨率的影響時，首先對裝調(diào)允許范圍內(nèi)的誤差進行采樣，然后結合ZEMAX和MATLAB通過FFT MTF計算[10-11]可以得到誤差引起的4f系統(tǒng)某一觀察對比度要求下的空間頻率f，分辨率為空間頻率的倒數(shù)。其流程如圖3所示。

圖3　MTF計算的流程圖

在實際分析過程中，同樣的元件誤差在不同透鏡上作用，引起的分辨率改變值是不同的。為了探究特定種類的誤差對整個系統(tǒng)分辨率的影響，這里采用“均方根”的方法。若一片透鏡的某一種誤差量為a，引起整個系統(tǒng)分辨率的變化值為x1，假設含有某誤差量為a的系統(tǒng)中有k片透鏡，則透鏡L1至Lk每一片透鏡分別帶有誤差a，系統(tǒng)所表現(xiàn)的分辨率誤差即為x1，x2，…，xk。通過均方根的形式推導出當系統(tǒng)含有誤差量a時分辨率的變化量，即分辨率為：

(16)

通過分析該方法設置系統(tǒng)中某透鏡含有的單一誤差，可以求得在特定誤差值下系統(tǒng)分辨率的均值，從而可以設置不同的誤差值來觀察分辨率隨誤差的變化。這對探究高分辨4f系統(tǒng)的分辨率影響因素及實際使用中誤差對系統(tǒng)分辨率的影響程度具有指導意義。具體的計算流程如圖4所示。

圖4　系統(tǒng)分辨率的計算流程圖

4　1 μm分辨率的4f系統(tǒng)分辨率影響因素實例分析

本文以非對稱型兩類高分辨率4f光學系統(tǒng)為例，分析了傾斜、偏心、離軸誤差方向對分辨率的影響。

如圖5所示，該系統(tǒng)是由9片透鏡、1片平板及光闌組合而成的非對稱型傅里葉透鏡。該系統(tǒng)物高視場為0～±25 mm，工作單波長為0.532 mm。其照片和測試圖分別如圖6和圖7所示。

圖5　1 μm高分辨率4f傳像系統(tǒng)

圖6　1 μm高分辨率4f傳像系統(tǒng)實物圖

圖7　1 μm高分辨率4f傳像系統(tǒng)測試圖

在沒有任何元件加工誤差與裝調(diào)誤差時，光學系統(tǒng)的MTF曲線如圖8所示。

圖8　1 μm高分辨率4f系統(tǒng)的MTF曲線圖

選取人眼觀察最低對比度0.03和對比度0.05觀察此時各視場的空間頻率，其結果如表2所示。

表21 μm 4f系統(tǒng)在0.03和0.05對比度下的分辨率

Tab.2Resolutions of 1 μm 4fsystem at contrasts of 0.03 and 0.05 respectively

(μm)

本文重點探究了光學面之間的空氣間隔、光學面的傾斜與偏心這3種影響因素對分辨率的影響。

該1 μm分辨率4f系統(tǒng)已經(jīng)在設計公差范圍內(nèi)安裝完畢。在公差范圍內(nèi)進行裝調(diào)誤差分析，如傾斜在15″以內(nèi)，偏心誤差在0.02 mm，位移在0.01 mm范圍內(nèi)。

在分析過程中，每一片透鏡、平板及光闌產(chǎn)生的傾斜誤差從0°開始每增加3″至15″時分辨率的變化規(guī)律。首先將透鏡L1傾斜3″測得此時系統(tǒng)的分辨率，其次將透鏡L2傾斜3″測得系統(tǒng)的分辨率。依次將L3至透鏡L9和平板1傾斜3″，通過計算得到在傾斜3″的情況下系統(tǒng)在0.03和0.05對比度條件下5個視場中系統(tǒng)的分辨率。當系統(tǒng)含有3″傾斜誤差時，每片透鏡傾斜帶來的分辨率改變量和系統(tǒng)分辨率的均值如表3所示。圖9中坐標軸的橫坐標表示每一片透鏡的傾斜角度，縱坐標表示系統(tǒng)在透鏡傾斜誤差產(chǎn)生時分辨率的均方值。

表31 μm 4f系統(tǒng)在對比度為0.03和0.05下含有3″傾斜量的分辨率

Tab.3Resolutions of 1 μm 4fsystem containing 3″ tilt quantity at contrasts of 0.03 and 0.05 respectively

(μm)

續(xù)表3

通過該方法一次遞增傾斜量至15″，得到系統(tǒng)分辨率和傾斜量的關系曲線，如圖9所示。圖9中坐標軸的橫坐標表示每一片透鏡的傾斜角度，縱坐標表示系統(tǒng)在透鏡傾斜誤差產(chǎn)生時分辨率的均方值。

(a)0.03對比度分辨率變化圖

(b)0.05對比度分辨率變化圖

由圖9(a)和圖9(b)可以看出，總體上傾斜使分辨率下降，雖然在極小范圍內(nèi)分辨率有所提高，但總體趨勢仍為下降。

在實際裝調(diào)過程中，裝調(diào)技術可以保證偏心在0.01 mm以內(nèi)，因此可以計算所有透鏡在0.01 mm以內(nèi)的分辨率。首先在最大偏心為0.01 mm 時計算 MTF曲線中對比度為0.03及0.05時的空間頻率。

根據(jù)測量的均值分辨率繪制分辨率與偏心誤差的關系，如圖10所示。由圖10(a)和圖10(b)可以看出,偏心會導致分辨率下降。

(a)0.03對比度分辨率變化圖

(b)0.05對比度分辨率變化圖

(a)0.03對比度分辨率變化圖

(b)0.05對比度分辨率變化圖

同理，在0.01 mm的間隔位移公差下，根據(jù)測量的分辨率均值可以繪制分辨率與位移誤差的關系圖，如圖11所示。由圖11(a)和圖11(b)可以看出，透鏡間隔位移誤差會導致視場分辨率的降低，而且非對稱型的4f透鏡對透鏡間隔十分敏感。

5　結　論

本文從加工、裝調(diào)誤差入手，建立了4f光學系統(tǒng)的分辨率影響因素的數(shù)學模型，并對非對稱型1 μm分辨率的4f光學系統(tǒng)在傾斜、偏心、離軸等誤差方面進行分析，在裝調(diào)公差范圍之內(nèi)，得到了各影響因素與分辨率在0.03對比度和0.05對比度下的關系曲線。結果表明，傾斜、偏心、離軸下降均會導致分辨率的下降。該研究對高分辨率4f系統(tǒng)的裝調(diào)誤差預估其有參考意義[12]。

[1]劉力.傅里葉變換透鏡的設計[D]. 重慶：重慶大學，2003.

LIU L.TheDesignoftheFourierTransformLens[D]. Chongqing：Chongqing University， 2003.(in Chinese)

[2]殷蔚.傅里葉變換透鏡的設計[D]. 重慶：重慶大學，2009.

YIN W.TheDesignoftheFourierTransformLens[D]. Chongqing： Chongqing University， 2009.(in Chinese)

[3]呂乃光. 傅里葉光學[M].北京：機械工業(yè)出版社，2006.

Lü N G.FourierOptics[M]. Beijing：China Machine Press， 2006.(in Chinese)

[4]李立.光學小波變換的4f系統(tǒng)研究[D]. 重慶：重慶大學，2009.

LI L.Studyonthe4fSystemofOpticalWaveletTransform[D]. Chongqing：Chongqing University， 2009. (in Chinese)

[5]崔龍燮. 4f光學信息處理系統(tǒng)[J]. 華北電力學院學報，1995，22(1)：34-40.

CUI L B.4foptical information processing system [J].JournalofNorthChinaInstituteofElectricPower， 1995， 1.22(1)： 34-40.(in Chinese)

[6]龐志海，樊學武，陳欽芳，等. 失調(diào)三反消像散光學系統(tǒng)像差特性[J]. 紅外與激光工程，2013，42(3)：653-657.

PANG ZH H， FAN X W， CHEN Q F，etal.. Misalignment induced aberration characteristic of TMA optical system [J].InfraredandLaserEngineering，2013，42(3)：653-657.(in Chinese)

[7]吳曉波，安問斗，楊鋼.圖像測量系統(tǒng)中的誤差分析即提高測量精度的途徑[J]. 光學 精密工程，1997，5(1)：133-141.

WU X B， AN W D， YANG G. Error analysing and approaches of improving measuring precision in image measuring system [J].Opt.PrecisionEng.， 1997， 5(1)： 133-141.(in Chinese)

[8]盧智勇，孫建鋒，職亞楠，等.直視合成孔徑激光成像雷達內(nèi)發(fā)射場波前像差對成像的影響[J].光學學報，2014，34(7)： 0728001-10.

LU ZH Y，SUN J F， ZHI Y N，etal.. Influence of inner wave-front aberration on down-looking synthetic aperture imaging ladar imaging[J].ActaOpticaSinica， 2014，34(7)： 0728001-10.(in Chinese)

[9]于樹海， 王建立， 董磊， 等. 基于非均勻周期采樣的傅里葉望遠鏡時域信號采集方法[J]. 中國光學， 2013，6(3)： 395-401.

YU SH H，WANG J L，DONG L，etal.. Time region signal collecting method of Fourier telescopy based on non-uniform periodically sampling[J].ChineseOptics， 2013， 6(3)： 395-401.(in Chinese)

[10]史廣維，張新，張建萍，等.幾種失調(diào)反射光學系統(tǒng)的波像差特性[J]. 紅外與激光工程， 2012， 41(1)：160-166.

SHI G W， ZHANG X， ZHANG J P，etal.. Wavefront aberration characteristic of several misaligned reflective optical systems [J].InfraredandLaserEngineering， 2012， 41(1)：160-166.(in Chinese)

[11]朱猛，李翔宇，龍寧波，等.大視場雙縫載頻散斑干涉成像檢測系統(tǒng)[J]. 光學 精密工程，2014，22(1)：13-17.

ZHU M， LI X Y， LONG N B，etal.. Double-slit based carrier frequency speckle interferometric system with large view in field [J].Opt.PrecisionEng.，2014，22(1)：13-17.(in Chinese)

[12]韓亮，田逢春，徐鑫，等.Mallat 算法的光學實現(xiàn)[J]. 光學 精密工程，2008，16(8)：1490-1499.

HAN L，TIAN F CH， XU X，etal.. Optical implementation of Mallat algorithm [J].Opt.PrecisionEng.， 2008，16(8)：1490-1499.(in Chinese)

夏春秋(1991-)，女，山東濰坊人，博士研究生，2008年于中國海洋大學獲得學士學位，主要研究方向為光學設計，光學系統(tǒng)的分析。E-mail：chunq.xia@hotmail.com

導師簡介：

鐘興(1982-)，男，四川自貢人，研究員，2004年于吉林大學獲得學士學位，2009于中科院長春光學精密機械與物理研究所獲得博士學位，研究方向為光學儀器的設計與分析。E-mail：ciomper@163.com

(版權所有未經(jīng)許可不得轉載)

Analysis of influence factors of resolution inhigh-resolution 4f imaging system

XIA Chun-qiu1， 2*， ZHONG Xing2，3， LIU Chun-yu2， HAN Pei-xian1，2， JIN Guang2

(1.UniversityofChineseAcademyofSciences，Beijing100049，China；2.ChangchunInstituteofOpticsandFineMechanics，ChineseAcademyofSciences，Changchun130033，China；3.ChangchunSatelliteTechnologyCo.LTD，Changchun130031，China)

*Correspondingauthor，E-mail：chunq.xia@hotmail.com

The factors and its mechanism which influence 4fimaging systemic resolution in practice were studied to guide the design and adjustment of the system. Firstly, the influence of the aberration caused by adjustment error on the system optical transfer function (MTF) in theory, whereby the mathematic model for describing the relationship of them was established. Based on the model, the adjustment error of 4foptical imaging system of 1 μm resolution was derived, and the influence of lens eccentricity, incline and displacement on the MTF was simulated, thus the variations of the resolution with the errors were achieved. The result indicates that within the incline of 15″, the eccentricity error is 0.02 mm and the displacement is within 0.01 mm. The adjustment error results in the decline of the resolution. Although the resolution increases while the incline within a minimum scope, the total trend will not change. The relationship figure of adjustment error and resolution can provide the reference for the error estimation of high-resolution 4fimaging system.

information optics; 4ftransmission system; alignment; mathematical model; curve analysis; high resolution

2015-08-03；

2015-10-06.

中科院青促會專項

1004-924X(2016)07-1573-09

TP391.41；O435.2

Adoi:10.3788/OPE.20162407.1573