基于細(xì)菌覓食算法的波前編碼相位掩模的優(yōu)化研究

陸 宇，張星煦，趙廷玉

(浙江理工大學(xué)理學(xué)院物理系,杭州 310018)

?

基于細(xì)菌覓食算法的波前編碼相位掩模的優(yōu)化研究

陸 宇，張星煦，趙廷玉

(浙江理工大學(xué)理學(xué)院物理系,杭州 310018)

為延拓光學(xué)成像系統(tǒng)的景深，獲得對(duì)物距不敏感的波前編碼光學(xué)成像系統(tǒng)，需求得最佳的相位掩模參數(shù)。以經(jīng)典的三次相位掩模為研究對(duì)象，采用細(xì)菌覓食算法，優(yōu)化了焦距為200.0 mm、入瞳直徑為25.4 mm的雙膠合波前編碼系統(tǒng)。仿真和實(shí)驗(yàn)結(jié)果表明：未加相位掩模的普通光學(xué)系統(tǒng)景深僅為30.0 mm，而插入優(yōu)化后的相位掩模的波前編碼系統(tǒng)景深可達(dá)80.0 mm，大大拓展了景深。細(xì)菌覓食算法優(yōu)化波前編碼系統(tǒng)的相位掩模有著運(yùn)算速度快、易跳出局部最優(yōu)解的優(yōu)點(diǎn)。

細(xì)菌覓食算法；波前編碼；相位掩模；優(yōu)化；景深

0 引 言

光學(xué)成像系統(tǒng)的景深延拓是應(yīng)用光學(xué)領(lǐng)域的熱點(diǎn)問(wèn)題。Dowski等提出的波前編碼成像技術(shù)常用于延拓非相干光學(xué)系統(tǒng)的景深[1]。波前編碼光學(xué)成像技術(shù)是在光學(xué)系統(tǒng)中引入自由曲面形式的相位掩模對(duì)物體進(jìn)行光學(xué)編碼，使得系統(tǒng)對(duì)物距的變化不敏感，即不同物距下光學(xué)傳遞函數(shù)(optical transfer function,OTF)和點(diǎn)擴(kuò)散函數(shù)(point spread function,PSF)保持不變[2]。因此物距不同的點(diǎn)，在成像面上會(huì)得到模糊程度一致的圖像。利用數(shù)字濾波技術(shù)對(duì)模糊圖像處理可獲得最終的清晰圖像。波前編碼技術(shù)在不影響成像分辨率的前提下，可大大拓展光學(xué)系統(tǒng)的景深，增加圖像的信息量，并且有效抑制了離焦相關(guān)的像差，減少光學(xué)元件個(gè)數(shù)，降低成本。這項(xiàng)技術(shù)特別適用于實(shí)時(shí)顯示立體動(dòng)態(tài)物體，在顯微成像、手機(jī)攝像、紅外空空導(dǎo)彈等方面具有廣泛的應(yīng)用前景。

自由曲面形式的相位掩模的優(yōu)化設(shè)計(jì)是實(shí)現(xiàn)波前編碼系統(tǒng)的關(guān)鍵技術(shù)，而合適的優(yōu)化算法是相位掩模優(yōu)化設(shè)計(jì)的重要組成部分。目前，常用于優(yōu)化相位掩模的算法有遺傳算法[3-4]、模擬退火算法[5]等。傳統(tǒng)優(yōu)化算法從單個(gè)初始值開始迭代尋求最優(yōu)解，算法簡(jiǎn)單但容易陷入局部最優(yōu)解。模擬退火算法在傳統(tǒng)算法的基礎(chǔ)上引入一定的概率突跳特性，在解空間中隨機(jī)尋找目標(biāo)函數(shù)的全局最優(yōu)解。該算法從較高的初始溫度開始，以一定的步徑降低系統(tǒng)溫度至預(yù)設(shè)的最低溫度結(jié)束。但使用該算法時(shí)產(chǎn)生的新解與初始解與步徑大小有關(guān)：當(dāng)步徑較小時(shí)，該算法難以跳出局部最優(yōu)解；當(dāng)步徑較大時(shí)，求取最優(yōu)解的過(guò)程中數(shù)值波動(dòng)較大，得到的最優(yōu)解不穩(wěn)定；所以需要多次調(diào)試參數(shù)才能獲得最優(yōu)解。而遺傳算法不依賴于梯度信息，但可能因?yàn)樗阉鞣秶^小或變異太少導(dǎo)致早熟收斂，容易陷入局部最優(yōu)解，而無(wú)法得到全局最優(yōu)解。本文使用細(xì)菌覓食算法[6-7]對(duì)波前編碼系統(tǒng)的相位掩模參數(shù)進(jìn)行優(yōu)化，該算法具有較好的魯棒性，且在兼顧跳出局部最優(yōu)解的同時(shí)保證了算法的收斂速度，節(jié)約了光學(xué)系統(tǒng)設(shè)計(jì)中優(yōu)化所用的時(shí)間。

1 基于細(xì)菌覓食算法的相位掩模優(yōu)化

1.1 波前編碼系統(tǒng)理論基礎(chǔ)

波前編碼由光學(xué)編碼和數(shù)字圖像處理兩部分組成，系統(tǒng)框圖如圖1所示。圖1中相位掩膜、透鏡和CCD組成了光學(xué)編碼部分，通過(guò)添加相位掩膜，使得CCD探測(cè)到的圖像對(duì)物距不敏感，再對(duì)圖像進(jìn)行復(fù)原，獲得最終的清晰圖像。

圖1 波前編碼系統(tǒng)框圖

波前編碼理論以光學(xué)傳遞函數(shù)為理論指導(dǎo)，采用兩次變換法[8]可以推導(dǎo)光學(xué)傳遞函數(shù)，分析其對(duì)物距的不敏感性。波前編碼系統(tǒng)的光瞳函數(shù)P(x,y,ψ)可以表示為：

(1)

(2)

其中：λ為光的波長(zhǎng)；D為光學(xué)系統(tǒng)的入瞳直徑；f為焦距；do和di分別為物距和像距。

在傅里葉光學(xué)理論中，光學(xué)系統(tǒng)的點(diǎn)擴(kuò)散函數(shù)表征了光學(xué)系統(tǒng)在空間域上的成像特性，它與光瞳函數(shù)的關(guān)系可以表示為：

(3)

其中：κ為常數(shù)，在實(shí)際計(jì)算時(shí)通常忽略不計(jì)[9]；F表示傅里葉變換。

光學(xué)傳遞函數(shù)是點(diǎn)擴(kuò)散函數(shù)的傅里葉變換，其表達(dá)式為：

(4)

調(diào)制傳遞函數(shù)(modulation transfer function,MTF)是輸出像與輸入像的對(duì)比度之比，可表示為光學(xué)傳遞函數(shù)的模：

(5)

從上面的分析可知，光學(xué)傳遞函數(shù)和點(diǎn)擴(kuò)散函數(shù)兩者都表示光學(xué)系統(tǒng)的成像能力，前者在頻域中表示，后者在空間域中表示。從式(4)可知，兩者之間存在傅立葉變換關(guān)系。因此，如果一個(gè)系統(tǒng)的光學(xué)傳遞函數(shù)對(duì)物距不敏感，則其點(diǎn)擴(kuò)散函數(shù)也不敏感；反之亦然。考慮到光學(xué)傳遞函數(shù)是復(fù)函數(shù)，從計(jì)算的簡(jiǎn)便性出發(fā)，光學(xué)上通常用調(diào)制傳遞函數(shù)替代光學(xué)傳遞函數(shù)，即調(diào)制傳遞函數(shù)和點(diǎn)擴(kuò)散函數(shù)都可以用來(lái)表示對(duì)物距的敏感性。

1.2 評(píng)價(jià)函數(shù)

本文用調(diào)制傳遞函數(shù)來(lái)表示光學(xué)系統(tǒng)的成像能力。波前編碼光學(xué)系統(tǒng)所成的像對(duì)物距不敏感，即處于不同物距下的物體在像面的模糊程度一致，不同物距下的調(diào)制傳遞函數(shù)或點(diǎn)擴(kuò)散函數(shù)相同。與此同時(shí)，光學(xué)系統(tǒng)生成的中間圖像的模糊程度不能超過(guò)圖像處理的極限，系統(tǒng)在低頻段的MTF值不能太低，這樣模糊圖像可實(shí)現(xiàn)復(fù)原。

本文使用費(fèi)希爾信息量[10-11]評(píng)價(jià)調(diào)制傳遞函數(shù)的離焦不變性，綜合考慮數(shù)字圖像處理中圖像的復(fù)原質(zhì)量，對(duì)調(diào)制傳遞函數(shù)加以限制，得出評(píng)價(jià)函數(shù)如下：

(6)

其中：u為空間頻率，M(u,ψ)為光學(xué)系統(tǒng)的MTF，th為系統(tǒng)可接受的最小MTF。

1.3 基于細(xì)菌覓食算法的相位掩膜板優(yōu)化

在構(gòu)建評(píng)價(jià)函數(shù)的基礎(chǔ)上，需要優(yōu)化搜索算法求解最佳的相位掩膜的參數(shù)。為了保證對(duì)參數(shù)全局優(yōu)化和收斂速度，本文選擇了細(xì)菌覓食算法進(jìn)行優(yōu)化。

細(xì)菌覓食算法是模擬大腸桿菌覓食行為的一種新型算法。近年來(lái)已被用于解決實(shí)際中的數(shù)學(xué)問(wèn)題[12-13]。細(xì)菌覓食包含趨化、復(fù)制和驅(qū)散三種行為方式。趨化是指細(xì)菌通過(guò)翻轉(zhuǎn)和前進(jìn)，向營(yíng)養(yǎng)豐富的地方聚集。翻轉(zhuǎn)是指細(xì)菌向隨機(jī)方向移動(dòng)，前進(jìn)是指根據(jù)翻轉(zhuǎn)后細(xì)菌的評(píng)價(jià)函數(shù)值，決定細(xì)菌接下來(lái)的是否按照原方向繼續(xù)移動(dòng)。復(fù)制是根據(jù)評(píng)價(jià)函數(shù)值，選擇較差的細(xì)菌繼承較好細(xì)菌的位置，加快尋優(yōu)速度。細(xì)菌通過(guò)趨化和復(fù)制可獲得局部最優(yōu)解。驅(qū)散是讓細(xì)菌以一定概率出現(xiàn)在搜索空間的任意位置，可跳出局部最優(yōu)解，從而尋找全局最優(yōu)解。

細(xì)菌覓食算法求解相位板的最佳參數(shù)的流程如圖2所示，步驟如下：

a) 初始化參數(shù)。設(shè)置相位板參數(shù)的搜索范圍(αmin,αmax)，細(xì)菌總數(shù)S，趨化步數(shù)Nc，單項(xiàng)運(yùn)動(dòng)最大次數(shù)Ns，復(fù)制次數(shù)Nre，驅(qū)散次數(shù)Ned，驅(qū)散概率Ped。

b) 隨機(jī)產(chǎn)生S組初始解P0，不滿足限制條件MTF>th的解重新生成。

c) 驅(qū)散次數(shù)加1，判斷是否達(dá)到驅(qū)散結(jié)束條件。是，輸出解P0；否，進(jìn)行下一步。

d) 復(fù)制次數(shù)加1，判斷是否達(dá)到復(fù)制結(jié)束條件。是，評(píng)價(jià)函數(shù)值較大的S/2組細(xì)菌死亡，剩下的S/2組細(xì)菌一分為二，保留細(xì)菌總數(shù)S，跳至步驟c；否，進(jìn)行下一步。

e) 趨化次數(shù)加1，判斷是否達(dá)到趨化結(jié)束條件。是，轉(zhuǎn)步驟f；否，轉(zhuǎn)步驟g。

f)S組解以一定概率隨機(jī)出現(xiàn)在解空間內(nèi)，轉(zhuǎn)步驟d。

g) 計(jì)算解P0的評(píng)價(jià)函數(shù)J0(1∶S)。

h) 產(chǎn)生新解P1(1∶S),并計(jì)算其評(píng)價(jià)函數(shù)J1(1∶S)。

i) 比較J0和J1中各組值，將評(píng)價(jià)函數(shù)值較好的解保存至P(0∶S)，再將最優(yōu)值賦予P0(0∶S)。

j) 運(yùn)動(dòng)步數(shù)加1，判斷是否達(dá)到單項(xiàng)最大運(yùn)動(dòng)次數(shù)。是，轉(zhuǎn)步驟e；否，轉(zhuǎn)步驟g。

圖2 基于細(xì)菌覓食算法的波前編碼系統(tǒng)相位板參數(shù)優(yōu)化流程

2 優(yōu)化結(jié)果與分析

針對(duì)焦距為200.0mm，通光孔徑為25.4mm的雙膠合透鏡(GCL-010606,大恒光電)光學(xué)系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)。以經(jīng)典的三次相位掩模為例，在光學(xué)系統(tǒng)的光瞳面插入相位掩模，使用細(xì)菌覓食算法進(jìn)行優(yōu)化設(shè)計(jì)?？紤]到圖像的大部分能量都集中在低頻段，取umax=70線對(duì)/mm，為確保模糊圖像可被復(fù)原，設(shè)閾值th=0.1。初始參數(shù)取解空間為[0,100]，細(xì)菌總數(shù)S=8，趨化步數(shù)Nc=3，單項(xiàng)運(yùn)動(dòng)最大次數(shù)Ns=3，復(fù)制次數(shù)Nre=3，驅(qū)散次數(shù)Ned=3，驅(qū)散概率Ped=0.3。使用細(xì)菌覓食算法進(jìn)行多次優(yōu)化，得到的相位板參數(shù)值的波動(dòng)在2%以內(nèi)，較為穩(wěn)定。多次優(yōu)化取平均值，得出最佳相位板參數(shù)：α=66.21。圖3為插入優(yōu)化后的三次相位掩膜的波前編碼系統(tǒng)在不同物距下的MTF曲線?？梢?，插入優(yōu)化后的相位掩模后，雖然光學(xué)系統(tǒng)的MTF值下降，但各種物距下的MTF曲線基本重合，并且MTF沒(méi)有零點(diǎn)，意味著低頻與高頻的圖像信息沒(méi)有丟失。

圖3 三次相位板波前編碼MTF曲線

本文仿真了傳統(tǒng)的光學(xué)系統(tǒng)，即只有雙膠合透鏡沒(méi)有相位掩膜板的情況作為對(duì)比。光學(xué)系統(tǒng)在不同物距下的調(diào)制傳遞函數(shù)曲線如圖4所示。由圖4可見，該系統(tǒng)對(duì)1000.0 mm的物距清晰成像，隨著物距do偏離量的增加，系統(tǒng)的MTF曲線變化明顯，MTF值迅速下降并出現(xiàn)零點(diǎn)，造成了圖像信息丟失。

圖4 相位掩模加入前的傳統(tǒng)光學(xué)系統(tǒng)MTF曲線

3 實(shí)驗(yàn)結(jié)果

為了驗(yàn)證上述的仿真結(jié)果，本文搭建了實(shí)驗(yàn)裝置，如圖5所示。相位掩模采用Holoeye空間光調(diào)制器(spatial liquid modulation，SLM)(型號(hào)：SLM-0017)實(shí)現(xiàn)，SLM前后的偏振片是為了實(shí)現(xiàn)相位調(diào)制而采用[14]。從左往右依次為CCD相機(jī)、偏振片、雙膠合透鏡、SLM-0017、偏振片、被測(cè)物。被測(cè)書本傾斜放置，使得每行字與光學(xué)系統(tǒng)之間的距離不同，便于比較不同物距下的成像質(zhì)量。

1.CCD相機(jī)；2.偏振片；3.雙膠合透鏡；4.空間光調(diào)制器；5.偏振片；6.被測(cè)物圖5 實(shí)驗(yàn)實(shí)物照片

考慮PSF容易測(cè)得，下面用PSF評(píng)價(jià)實(shí)際光學(xué)系統(tǒng)的成像特征，像面捕捉到的傳統(tǒng)光學(xué)系統(tǒng)(未加相位掩模)與波前編碼系統(tǒng)(加入相位掩模)的PSF如圖6和圖7所示。可以發(fā)現(xiàn)，傳統(tǒng)光學(xué)系統(tǒng)在對(duì)焦時(shí)，其PSF為一個(gè)銳利的點(diǎn)，但當(dāng)離焦時(shí)，PSF退化為彌散斑，且隨著物距偏離量的增大，彌散斑越來(lái)越大，導(dǎo)致成像性能惡化，直至無(wú)法分辨圖像。但對(duì)于插入相位掩膜的波前編碼光學(xué)成像系統(tǒng)，在對(duì)焦時(shí)，PSF為彌散斑(形狀和大小由相位掩膜板形式?jīng)Q定)，雖然會(huì)導(dǎo)致模糊，但隨著物距偏離量的增大，PSF幾乎不變，這樣便可以使用同一種濾波器進(jìn)行復(fù)原，得到清晰圖像。

圖6 傳統(tǒng)光學(xué)系統(tǒng)不同物距下的PSF

圖7 波前編碼系統(tǒng)不同物距下的PSF

圖8為傳統(tǒng)光學(xué)系統(tǒng)和波前編碼系統(tǒng)的實(shí)際成像比較。傳統(tǒng)光學(xué)系統(tǒng)實(shí)際成像如圖8(a)所示，從圖中可以看出，被測(cè)書本中間部分的文字清晰可見但隨著物距偏離量的增大，圖像的模糊程度越來(lái)越大，上下兩端的文字已經(jīng)完全模糊，無(wú)法進(jìn)行分辨。而使用插入優(yōu)化后的相位掩膜板的波前編碼系統(tǒng)，圖像的模糊程度一致，如圖8(b)所示。導(dǎo)致圖像模糊的PSF如圖7所示，以此作為反卷積核，使用簡(jiǎn)單的維納濾波算法[15]可復(fù)原出較為清晰的最終圖像如圖8(c)所示，可以看到整幅圖像都清晰可見，實(shí)現(xiàn)了大景深成像。未加相位掩模的光學(xué)系統(tǒng)景深為30.0 mm，而加入相位掩模的雙膠合系統(tǒng)的清晰成像的景深延拓至80.0 mm，實(shí)現(xiàn)了景深的延拓。

圖8 傳統(tǒng)光學(xué)系統(tǒng)與波前編碼系統(tǒng)成像的比較

4 細(xì)菌覓食算法和模擬退火算法的對(duì)比

本文分別使用細(xì)菌覓食算法和經(jīng)典的模擬退火算法優(yōu)化相位掩模并對(duì)比優(yōu)化結(jié)果。后者在局部尋優(yōu)的基礎(chǔ)上，在當(dāng)前解的鄰域內(nèi)隨機(jī)產(chǎn)生一個(gè)新解，在溫度參數(shù)的控制下，除了接受優(yōu)化解之外，還根據(jù)Metropolis準(zhǔn)則，以一定的概率接受惡化解。這使得該算法具有跳出局部最優(yōu)解的能力。但產(chǎn)生的新解大小與預(yù)設(shè)的步徑或領(lǐng)域有關(guān)，預(yù)設(shè)值較小時(shí)算法跳不出局部最優(yōu)解；當(dāng)該預(yù)設(shè)值較大時(shí)可能導(dǎo)致跳過(guò)全局最優(yōu)解或多次所得最優(yōu)解的數(shù)值波動(dòng)較大。該算法需要多次調(diào)試找到合適的參數(shù)才能快速準(zhǔn)確的得到全局最優(yōu)解。表1中給出初始溫度T0，截止溫度T1，步徑m與最佳參數(shù)α和所用時(shí)間t的關(guān)系。

表1 模擬退火算法不同參數(shù)的優(yōu)化比較

由表1可知，T1和T0不變時(shí)，當(dāng)m較大時(shí)，最優(yōu)解α波動(dòng)較大；當(dāng)m較小時(shí)，α難以跳出局部最優(yōu)解。當(dāng)m不變，溫度差較大時(shí)，程序運(yùn)行時(shí)間增加；當(dāng)T0較小時(shí)，陷入局部最優(yōu)解。需多次調(diào)整才可快捷準(zhǔn)確的使用?？傮w而言，模擬退火算法容易陷入局部最優(yōu)解，而較難得到全局最優(yōu)解，如表1的第三列和最后兩列，優(yōu)化得到的參數(shù)均超過(guò)40，和最優(yōu)解66.21相差很大。

而細(xì)菌覓食算法求解相位掩膜的參數(shù)魯棒性強(qiáng)，使用簡(jiǎn)單，產(chǎn)生的新解為解空間的隨機(jī)數(shù)，與步徑無(wú)關(guān)。算法運(yùn)行時(shí)間與細(xì)菌總數(shù)，趨化步數(shù)，單項(xiàng)運(yùn)動(dòng)最大次數(shù)，復(fù)制次數(shù)，驅(qū)散次數(shù)等有關(guān)，參數(shù)的設(shè)置較易。細(xì)菌算法用時(shí)和優(yōu)化參數(shù)如表2所示。

表2 細(xì)菌覓食算法不同參數(shù)的優(yōu)化比較

由表2可知，盡管細(xì)菌覓食算法使用的初始參數(shù)不同，但優(yōu)化結(jié)果基本在全局最優(yōu)解附近波動(dòng)。也就是說(shuō)，使用細(xì)菌覓食算法進(jìn)行優(yōu)化時(shí)，各運(yùn)算步驟較小時(shí)便可求得較穩(wěn)定的最優(yōu)解，且多次優(yōu)化的最優(yōu)解α的波動(dòng)較小。由于該算法是在解空間內(nèi)全局搜索，所用極易跳出局部最優(yōu)解。

5 結(jié) 論

本文提出的基于細(xì)菌覓食算法的波前編碼系統(tǒng)相位板的優(yōu)化，以費(fèi)希爾信息量評(píng)價(jià)系統(tǒng)MTF的離焦不敏感性，優(yōu)化得出了三次相位板的最佳參數(shù)。細(xì)菌覓食算法具有良好的魯棒性，在對(duì)搜索范圍進(jìn)行合適的調(diào)整后，便可快速有效地進(jìn)行運(yùn)算，具有群體智能算法并行搜索、易跳出局部極小值等優(yōu)點(diǎn)，以較短的時(shí)間求得相位板參數(shù)全局最優(yōu)解。對(duì)焦距f=200.0 mm，入瞳直徑D=25.4 mm的雙膠合成像系統(tǒng)進(jìn)行成像實(shí)驗(yàn)，插入優(yōu)化后的相位掩膜板，不同物距的點(diǎn)擴(kuò)散函數(shù)具有極高的相似性。實(shí)驗(yàn)結(jié)果表明像面探測(cè)到的圖像在不同物距下模糊程度具有一致性，使用同一個(gè)的維納濾波器進(jìn)行圖像處理可獲得最終的清晰圖像，最終將雙膠合系統(tǒng)的清晰成像的景深延拓至80.0 mm，是未加相位掩模的傳統(tǒng)系統(tǒng)的2.7倍，實(shí)現(xiàn)了景深延拓的目的。

[1] DOWSKI E R, CATHEY W T. Extended depth of field through wave-front coding[J]. Applied Optics,1995,34(11):1859-66.

[2] LIU J, MIAO E, SUI Y, et al. Optimized non-integer order phase mask to extend the depth of field of an imaging system[J]. Optics Communications,2016,374:92-96.

[3] 黃薇薇,葉子,張文宇,等.基于遺傳算法的波前編碼相位板參數(shù)優(yōu)化[J].光學(xué)儀器,2007,29(4):17-22.

[4] 閆鋒.基于調(diào)制傳遞函數(shù)相似度的波前編碼系統(tǒng)優(yōu)化技術(shù)[J].光學(xué)學(xué)報(bào),2013(4):68-73.

[5] 張文字,陳燕萍,趙廷玉,等.一種基于OTF穩(wěn)定性的波前編碼相位板優(yōu)化方法[J].光子學(xué)報(bào),2007,36(11):2017-2021.

[6] PASSINO K M. Biomimicry of bacterial foraging for distributed optimization and control[J]. IEEE Control Systems,2002,22(3):52-67.

[7] ELATTAR E E. Adaptive bacterial foraging and genetic algorithm for unit commitment problem with ramp rate constraint[J]. International Transactions on Electrical Energy Systems,2015,26(7):1555-1569.

[8] 莊松林.光學(xué)傳遞函數(shù)[M].北京:機(jī)械工業(yè)出版社,1981:276-282.

[9] 陳燕平.景深延拓波前編碼系統(tǒng)理論及其應(yīng)用研究[D].杭州:浙江大學(xué),2007:6-7.

[10] PRASAD S. Fisher-information-based analysis of a phase-diversity-speckle imaging system[J]. Journal of the Optical Society of America A Optics Image Science & Vision,2004,21(11):pp.2073-88.

[11] MOTKA L, STOKLASA B, D’ANGELO M, et al. Optical resolution from Fisher information[J]. The European Physical Journal Plus,2016,131(5):1-13.

[12] 李珺,黨建武.細(xì)菌覓食算法求解高維優(yōu)化問(wèn)題[J].計(jì)算機(jī)應(yīng)用研究,2016,33(4):1024-1027.

[13] 王飛,邊會(huì)媛,韓雪,等.基于細(xì)菌覓食算法的砂礫巖巖性識(shí)別方法[J].地球科學(xué)與環(huán)境學(xué)報(bào),2016,38(2):277-284.

[14] 陳棟.基于液晶空間光調(diào)制的波前可控技術(shù)研究[D].南京:南京理工大學(xué),2011:17-19.

[15] 胡小平,陳國(guó)良,毛征宇,等.離焦模糊圖像的維納濾波復(fù)原研究[J].儀器儀表學(xué)報(bào),2007,28(3):479-482.

(責(zé)任編輯： 康 鋒)

The Optimization of the Phase Mask in a Wavefront Coding System Based on Bacterial Foraging Algorithm

LUYu,ZHANGXingxu,ZHAOTingyu

(Department of Physics, School of Sciences, Zhejiang Sci-Tech University, Hangzhou 310018, China)

In order to extend the depth of field of an optical system and acquire a wavefront coding system insensitive to the object distance, the optimal parameters of the phase mask need to be obtained. A doublet system with a classical cubic phase mask is optimized based on bacterial foraging algorithm. The focal length is 200.0 mm and entrance pupil diameter is 25.4 mm in the doublet system. The simulation and experimental results show that the depth of field in traditional optical system is only 30.0 mm, while it is 80.0 mm in the wavefront coding system by adding the optimized phase mask. It greatly extends the depth of field. Optimization using bacterial foraging algorithm have a fast speed and easy to jump out of local optimal solution.

bacterial foraging algorithm; wavefront coding; phase mask; optimization; depth of field

10.3969/j.issn.1673-3851.2017.07.023

2017-02-26 網(wǎng)絡(luò)出版日期： 2017-05-24

國(guó)家自然科學(xué)基金項(xiàng)目(61205158)；浙江省自然科學(xué)基金項(xiàng)目(LY15F050013)

陸宇(1991-)，男，安徽蚌埠人，碩士研究生，主要從事光學(xué)設(shè)計(jì)方面的研究。

趙廷玉，E-mail：zhaotingyu@zstu.edu.cn

TS195.644

A

1673- 3851 (2017) 04- 0604- 07