張寶武，霍劍鋒，張明月，劉媛媛，余桂英，姚蘆鷸

(1.中國(guó)計(jì)量大學(xué) 計(jì)量測(cè)試工程學(xué)院，浙江 杭州 310018；2.比薩大學(xué)物理系，意大利 比薩 56127)

基于VirtualLab Fusion的TDOCT系統(tǒng)合成光源仿真

張寶武1,2，霍劍鋒1，張明月1，劉媛媛1，余桂英1，姚蘆鷸1

(1.中國(guó)計(jì)量大學(xué) 計(jì)量測(cè)試工程學(xué)院，浙江 杭州 310018；2.比薩大學(xué)物理系，意大利 比薩 56127)

基于VirtualLab Fusion構(gòu)建了時(shí)域光學(xué)相干層析技術(shù)(time domain optical coherence tomography，TDOCT)光路系統(tǒng)，導(dǎo)入不同顆數(shù)發(fā)光二極管(light emitting diode，LED)的合成光譜，從參考鏡的位置變化中獲得干涉信號(hào)，以此分析出合成光源的半高寬和相干度，以及對(duì)應(yīng)TDOCT的縱向分辨率.結(jié)果表明，合成光源的最大強(qiáng)度值隨著LED顆數(shù)從1到10的增加先增加后趨向不變；而半高寬隨之不斷增加；TDOCT系統(tǒng)的縱向分辨率隨之不斷提高.當(dāng)LED增加至3顆的時(shí)候，合成光源已經(jīng)可以使TDOCT系統(tǒng)縱向分辨率提高至9.8 μm，這與一顆超輻射發(fā)光二極管(super luminescent light-emitting diode，SLD)給出的分辨率級(jí)別相當(dāng).

光學(xué)相干層析技術(shù)；VirtualLab Fusion軟件；發(fā)光二極管；合成光源

光學(xué)相干層析技術(shù)(OCT)因具有高分辨率、非接觸、無(wú)損傷、速度快，以及能夠?qū)铙w生物進(jìn)行實(shí)時(shí)成像等特點(diǎn)，于是被迅速應(yīng)用于生物醫(yī)學(xué)的影像檢測(cè)，并正在被拓展到工業(yè)測(cè)量等領(lǐng)域.1991年Huang等首次提OCT概念[1]，應(yīng)用它進(jìn)行了人眼視網(wǎng)膜的微結(jié)構(gòu)和冠狀動(dòng)脈臂的成像，從而推動(dòng)了OCT技術(shù)在生物醫(yī)學(xué)領(lǐng)域的應(yīng)用.1997年Colston等人用OCT技術(shù)取得了離體豬前磨牙的牙本質(zhì)以及牙周組織圖像[2]，證明了OCT技術(shù)在牙科領(lǐng)域的應(yīng)用潛力[3].另外，OCT系統(tǒng)可以對(duì)人體的表皮和真皮超細(xì)微結(jié)構(gòu)[4]，胃腸道表層結(jié)構(gòu)[5]，對(duì)冠狀動(dòng)脈組織[6]等進(jìn)行高分辨成像，以觀察病變組織的情況，從而實(shí)現(xiàn)對(duì)疾病的診斷.近年來(lái)，隨著OCT技術(shù)的發(fā)展，又產(chǎn)生了許多新的研究熱點(diǎn)：一方面多種方式控制OCT光源系統(tǒng)的偏振、相位和光譜等信號(hào)，發(fā)展出眾多功能型OCT系統(tǒng)[7]，另一方面擴(kuò)展OCT的應(yīng)用領(lǐng)域[8].

對(duì)于OCT系統(tǒng)來(lái)說(shuō)，光源相干長(zhǎng)度是其縱向分辨率的重要決定因素.一般要求其光源具有較寬的頻譜寬度、輸出功率高、穩(wěn)定性好、易于耦合等特點(diǎn).為此，如何改進(jìn)光源，擴(kuò)大其光譜寬度以便獲得足夠短的相干時(shí)間，從而實(shí)現(xiàn)減小相干長(zhǎng)度的目的成為眾多研究者的關(guān)注目標(biāo).目前已研制出并應(yīng)用于OCT中的寬光譜光源有超輻射發(fā)光二極管(SLD)[9]、超短脈沖飛秒激光[10].其中，SLD光譜寬度在40～70 nm，OCT系統(tǒng)縱向分辨率可達(dá)5～20 μm；藍(lán)寶石鎖模飛秒激光的光譜寬度在350 nm，OCT系統(tǒng)縱向分辨率高達(dá)1 μm.不過(guò)，這些光源價(jià)格昂貴，而且功率越高，光源噪聲越大.經(jīng)研究，單個(gè)普通發(fā)光二極管(LED)也能作為OCT系統(tǒng)的光源，并使OCT在一定的需求上取得較好的性能[11].另外，LED具有價(jià)格便宜、噪聲小和安全等特點(diǎn).但是，單個(gè)LED的不足之處在于其光譜寬度很窄(光譜半高寬度大約為30 nm).為了促進(jìn)LED光源在OCT中的使用，研究者提出了多個(gè)LED光源組合合成寬帶光源的應(yīng)用思路[11]，以此來(lái)提高OCT的分辨率和性價(jià)比.

為了驗(yàn)證多顆LED光源合成理論分析的正確性，在系統(tǒng)加工之前需要對(duì)其進(jìn)行虛擬仿真.在光學(xué)仿真方面，Matlab, ZEMAX, Tracepro等一直是主流工具，但是其缺點(diǎn)在于：第一，使用者需要知道光線追跡理論，并且還需要進(jìn)行數(shù)值計(jì)算的方法設(shè)計(jì)和分析，使用起來(lái)比較復(fù)雜，不適應(yīng)光學(xué)系統(tǒng)的快速可視化的設(shè)計(jì)仿真；第二，編制這樣的仿真程序工作量較大，難度較高，因此只限于少數(shù)人員使用.

VirtualLab Fusion[12]是由德國(guó)LightTrans公司研發(fā)的一款模塊化光學(xué)建模軟件，它模塊化的拖拽建模方式避免了繁冗復(fù)雜的理論推導(dǎo)和仿真模型的可視化編程.另外，它著眼于光源開始的整個(gè)光學(xué)系統(tǒng)，讓系統(tǒng)整體的虛擬運(yùn)行，從而避免了仿真對(duì)象的孤立性.這種全系統(tǒng)化的建模仿真思想為相關(guān)的光學(xué)設(shè)計(jì)和仿真提供了一種嶄新的技術(shù)方案.

本文基于VirtualLab Fusion，以時(shí)域光學(xué)相干層析技術(shù)(TDOCT)為對(duì)象，對(duì)不同多顆LED的合成光源進(jìn)行了仿真，復(fù)現(xiàn)了文獻(xiàn)[11]的結(jié)果，驗(yàn)證了理論的正確性和方法的可行性.

1 多顆LED合成光源TDOCT縱向分辨率

利用多顆LED組合光源的TDOCT系統(tǒng)如圖1.其中光源部分由N顆LED光源組成，每顆LED發(fā)出的光分別經(jīng)過(guò)獨(dú)立的光纖輸送到同一個(gè)N：1型的波分復(fù)用器，合成后形成一個(gè)干涉長(zhǎng)度短并具有高度部分相干的寬帶光源.這個(gè)合成光經(jīng)過(guò)一個(gè)光纖邁克爾遜干涉儀以后，最后被光電探測(cè)器接收.如果經(jīng)參考臂的光程與經(jīng)樣品臂的光程匹配到光源的相干長(zhǎng)度之內(nèi)，樣品的反射光就與參考臂上返回的參考光相干涉.然后，光電探測(cè)器的信號(hào)再由處理電路進(jìn)行放大、濾波、選頻、解調(diào)，最后再經(jīng)過(guò)數(shù)據(jù)采集與處理由計(jì)算機(jī)讀取、處理和顯示，完成樣品內(nèi)部結(jié)構(gòu)信息檢測(cè).

1 多顆LED組合光源的TDOCT系統(tǒng)Figure 1 TDOCT system with synthetic optical source of LED

由文獻(xiàn)[13]推導(dǎo)得，當(dāng)同時(shí)使用N個(gè)互不相干的高斯型LED光源合成時(shí)，干涉儀輸出的相干光強(qiáng)交流部分為

(1)

式(1)中，Δλ為合成光源的光譜半高寬；λ0為合成光源的中心頻率；LC定義為單顆LED的相干長(zhǎng)度；ΔL為參考光和信號(hào)光之間的光程差.另外，為了方便，(1)式中已經(jīng)設(shè)定每個(gè)光源的相干長(zhǎng)度均相等，且光強(qiáng)相同，每個(gè)光源的波長(zhǎng)為λi.從(1)式，可以看出，由多顆LED光源合成后的相干交流部分相當(dāng)于一個(gè)中心頻率為λ0的光源.光強(qiáng)交變部分除具有高斯函數(shù)外，還多出一個(gè)函數(shù)對(duì)高斯函數(shù)調(diào)制，影響高斯函數(shù)的包絡(luò)形狀.

由于TDOCT為背向光的探測(cè)技術(shù)，因此光要在組織中往返一次，所以，使用高斯光譜光源的TDOCT系統(tǒng)的縱向分辨率為合成光源相干長(zhǎng)度的一半.由此可見，OCT的縱向分辨率由合成光源的相干長(zhǎng)度決定，相干長(zhǎng)度越長(zhǎng)，分辨率越低.

2 基于VirtualLab Fusion的多顆LED合成光源仿真

2.1 TDOCT的VirtualLab Fusion的建模

基于VirtualLab Fusion構(gòu)建的TDOCT模型光路流程圖如圖2，(a)為元件平面布局；(b)為元件三維布局.其中“Gaussian Wave”代表系統(tǒng)的高斯型光源；“Ideal Beam Splitter”為理想光束分束器，這里設(shè)置為50/50非偏振型；M1為理想平面鏡，設(shè)定為OCT系統(tǒng)的參考鏡，M2為理想平面鏡，設(shè)置為TDOCT系統(tǒng)樣品；“Virtual Screen”為虛擬屏，用以顯示某一個(gè)截面上的光強(qiáng)分布；“Value Monitoring”為光強(qiáng)探測(cè)器，用來(lái)探測(cè)虛擬屏上某一點(diǎn)的光強(qiáng)變化.器件之間的相對(duì)距離可以通過(guò)單擊元件下面的坐標(biāo)來(lái)設(shè)置.圖中的光路連接方式即表示光在系統(tǒng)中的走向，光源發(fā)出的光經(jīng)過(guò)理想分束器后分別照射到參考鏡和樣品鏡上，之后，經(jīng)過(guò)兩個(gè)鏡反射的兩路光再匯合.這個(gè)匯合的光可以分別由虛擬屏和光強(qiáng)探測(cè)器接收.

圖2 VirtualLab Fusion對(duì)TDOCT系統(tǒng)的建模Figure 2 System of TDOCT in Virtuallab Fusion

完成TDOCT的建模以后，就可以通過(guò)改變初始條件來(lái)實(shí)現(xiàn)對(duì)系統(tǒng)的仿真.因?yàn)閂irtualLab Fusion自身沒有多顆LED合成光源，因此，我們首先需要通過(guò)光源編輯窗口導(dǎo)入已保存的特定光譜.對(duì)于單顆LED來(lái)說(shuō)，其光譜近似為(2)式表示的高斯型光譜函數(shù)

(2)

式(2)中，Δv0為光源譜寬，I0為光源強(qiáng)度系數(shù).利用頻率和波長(zhǎng)之間的關(guān)系即可獲得波長(zhǎng)表示的高斯型光譜函數(shù).基于此，運(yùn)用MATLAB得到不同顆數(shù)LED合成后的光源光譜，如圖3所示.其中LED顆數(shù)用N表示，第一顆LED光源的中心波長(zhǎng)為850 nm，其后相鄰兩顆LED光源之間的中心波長(zhǎng)相距10 nm，且每顆LED光源的半高寬為30 nm，強(qiáng)度以中心波長(zhǎng)850 nm對(duì)應(yīng)的強(qiáng)度為歸一化條件.從圖中可以看出，隨著LED顆數(shù)的增加，合成光源的強(qiáng)度隨之增加，譜寬隨之展寬.為了更細(xì)致地考察這種變化趨勢(shì)，我們獲取了圖3中每條曲線的最大值和半高寬，它們隨LED顆數(shù)的變化情況如圖4所示.從圖4中我們可以看出，當(dāng)LED顆數(shù)增加至6顆以后，合成光源的最大強(qiáng)度值不再增加，而半高寬持續(xù)增大.

圖3 不同顆數(shù)LED組合后的光源光譜Figure 3 Spectrum of synthetic optical source with different numbers of LED

圖4 合成光源的強(qiáng)度最大值和譜寬隨LED顆數(shù) 的變化情況Figure 4 Maximum intensity and FWHM of the synthetic optical source with different numbers of LED

然后，將圖3不同顆數(shù)LED的合成光源依次導(dǎo)入圖2(a)的光源中，如圖5.圖5(a)為中心波長(zhǎng)為850 nm，半高寬為30 nm的單顆LED能量波長(zhǎng)分布(部分)；圖5(b)為VirtualLab Fusion生成的合成光源，強(qiáng)度以中心波長(zhǎng)850 nm對(duì)應(yīng)的強(qiáng)度為歸一化條件.

2.2 合成光源干涉信號(hào)仿真

建模結(jié)束后，我們首先對(duì)比了850 nm單頻光源和和單顆LED光源情況下虛擬屏給出的干涉，如圖6.其中(a)為850 nm單頻光源的情況，(b)為中心波長(zhǎng)為850 nm單顆LED的情況.從中可以看出，當(dāng)光源為理想單頻光時(shí)，TDOCT將給出非常好的干涉圖形；當(dāng)光源是寬頻光源時(shí)，TDOCT將給出對(duì)比度不高的干涉條紋.

圖5 光源基本參數(shù)編輯窗口Figure 5 Parameter window for light source

圖6 不同光源下TDOCT給出的干涉圖形Figure 6 Interference pattern of TDOCT with different optical source

隨后，我們利用“Parameter Run”功能，實(shí)現(xiàn)參考鏡M1位置變化時(shí)合成光源干涉信號(hào)的仿真.其中參考鏡M1在其初始位置前后平移距離設(shè)置為-15～15 μm，步數(shù)為100.這也就是說(shuō)，TDOCT的兩條光的光程差有-15～15 μm的變化量.不同顆數(shù)LED合成光源后的干涉信號(hào)隨光程差的變化如圖7所示.從中可以看出，隨著LED顆數(shù)的增加，干涉信號(hào)所覆蓋的光程差縮短，即干涉長(zhǎng)度減小，相應(yīng)地縱向分辨率就會(huì)提高.圖8給出了不同顆數(shù)LED合成后，相應(yīng)TDOCT縱向分辨率的變化情況.從中可以看出，隨著LED顆數(shù)的增加，TDOCT的縱向分辨率隨之提高.當(dāng)LED為3顆的時(shí)候，縱向分辨率為9.8 μm，這已經(jīng)可ΔL以達(dá)到單個(gè)SLD光源的級(jí)別了.也就是說(shuō)，多顆LED合成光源可以用來(lái)代替單顆SLD光源.當(dāng)然，對(duì)于TDOCT來(lái)說(shuō)，光譜寬度越大，相干長(zhǎng)度越短，縱向分辨率越高；但是，不能無(wú)限制的增加LED顆數(shù)，因?yàn)殡S著LED顆數(shù)的增加，總光強(qiáng)會(huì)保持不變，這樣會(huì)增加干涉信號(hào)的背景信息，降低信噪比.所以，具體LED顆數(shù)還要根據(jù)實(shí)驗(yàn)來(lái)確定.

圖7 不同顆數(shù)LED組合時(shí)TDOCT的干涉波形Figure 7 Interference of the TDOCT with synthetic optical source of different numbers of LED

圖8 TDOCT縱向分辨率隨LED顆數(shù)變化的情況Figure 8 Vertical resolution of TDOCT varying with the number of LED

3 結(jié) 語(yǔ)

本文基于最新光學(xué)仿真軟件VirtualLab Fusion，以TDOCT為對(duì)象，仿真研究了多顆LED合成光源的參數(shù).研究結(jié)果表明，LED在一定顆數(shù)組合下所合成的光源能夠達(dá)到SLD光源的光譜半高寬和相干長(zhǎng)度，并且可以給出與SLD級(jí)別相當(dāng)?shù)腡DOCT的縱向分辨率.這種仿真結(jié)果是基于VirtualLab Fusion的系統(tǒng)虛擬運(yùn)行環(huán)境下獲得的.在驗(yàn)證多顆LED合成光源可行性的同時(shí)也為相關(guān)光學(xué)建模問題提供了一種整體虛擬仿真新思路.

(致謝：非常感謝訊技光電科技(上海)有限公司提供的VirtualLab Fusion密鑰和多年來(lái)的技術(shù)支持)

[1] CHAN K P, YAMADA M, DEVARAJ B, et al. Optical imaging through highly scattering media by use of heterodyne detection in the 1.3 μm wavelength region[J]. Optics Letters,1995,20(5):492-494.

[2] COLSTON B W, EVERETT M J, DA SILVA L B, et al. Optical coherence tomography for diagnosing periodontal disease[J]. SPIE,1997,2973(2):216-220.

[3] 石博雅，孟卓，劉鐵根,等.OCT系統(tǒng)對(duì)人體牙齒組織的非失真成像深度的研究[J].光學(xué)學(xué)報(bào)，2014,34(2):197-203. SHI Boya, MENG Zhuo, LIU Tiegen, et al. Non-distorted imaging depth of optical coherence tomography system in human dental tissues[J]. Acta Optica Sinica.,2014,34(2):197-203.

[4] ALEX A, POVAZAY B, HOFER B, et al. Multispectral in vivo three-dimensional optical coherence tomography of human skin[J]. Journal of Biomedical Optics,2010,15(2):431-440.

[5] OSIAC E, SAFTOIU A, GHEONEA D I, et al. Optical coherence tomography and Doppler optical coherence tomography in the gastrointestinal tract[J]. World Journal of Gastroenterology,2011,130(1):15-20.

[6] LIU Zhongguo，WANG Ruikang. Stripe motion artifact suppression in phase-resolved OCT blood flow images of the human eye based on the frequency rejection filter[J]. Chinese Optics Letters，2013，11(3)，51-55.

[7] WANG Yi，CHEN Xiaodong, LENG Ji，et al. Common-path OCT using optical delay stair[J]. Chinese Optics Letters，2015，13(1):55-56.

[8] 蘇亞，孟卓，于海民,等.OCT無(wú)創(chuàng)檢測(cè)技術(shù)的人體血糖平衡延遲時(shí)間研究[J].激光技術(shù)，2015,39(1):19-22. SU Ya, MENG Zhuo, YU Haimin, et al. Study on blood glucose lag time in noninvasive measurement using optical coherence tomography[J]. Laser Technology,2015,39(1):19-22.

[9] HUANG D, SWANSON E A, LIU C P, et al. Optical coherence tomography[J]. Science,1991,254(5035):1178-1181.

[10] DREXLER W, MORGNER U, KARTNER F X, et al. In vivo ultrahigh-resolution optical coherence tomography[J]. Optics Letters.,1999,24(17):1221-1223.

[11] 劉新文，王惠南.提高OCT分辨率的多LED合成新光源[J].南京航空航天大學(xué)學(xué)報(bào)，2005，37(4)：530-534. LIU Xinwen, WANG Huinan, Synthesis investigation of several LED Sources for OCT resolution improment[J]. Journal of Nanjing University of Aeronautics &Astronautics,2005,37(4):530-534.

[12] 王馳，畢書博，王利,等.超小自聚焦光纖探頭研究用場(chǎng)追跡數(shù)值模擬技術(shù)[J].物理學(xué)報(bào)，2013,62(2):340-349. WANG Chi, BI Shubo, WANG Li, et al. Field-tracing based numerical simulation technique for the investigation of ultra-small self-focusing optical fiber probe[J]. Acta Physica Sinica,2013,62(2):340-349.

[13] 姚蘆鷸. 提高OCT分辨率若干關(guān)鍵技術(shù)的研究[D]. 杭州：中國(guó)計(jì)量大學(xué)，2016. YAO Luyu. Study on several key techniques for improving resolution of optical coherence tomography[D]. Hangzhou: China Jiliang University,2016.

Simulations of synthetic light source of TDOCT based on VirtualLab Fusion

ZHANG Baowu1,2, HUO Jianfeng1, ZHANG Mingyue1, LIU Yuanyuan,YU Guiying1, YAO Luyu1

(1.College of Metrology and Measurement Engineering, China Jiliang University, Hangzhou 310018, China;2.Department of Physics, University of Pisa, Pisa 56127, Italy)

The light path system of time domain optical coherence tomography(TDOCT)was simulated with the software of VirtualLab fusion. The light source was synthesized with multi LED lights. The interference signal of the TDOCT was obtained by moving the reference mirror. The full width at half maximum (FWHM) and coherence length of the synthetic light source, the corresponding vertical resolutions of the TDOCT were analyzed. The results showed that the maximum intensity of the synthetic light source increased at the beginning and then kept stable with the LED lights increased from 1 to 10, while the vertical resolution of the TDOCT kept increasing. When the synthetic light source was mixed with 3 LED lights, the vertical resolution of the TDOCT was 9.8μm, which was similar to that by using one SLD.

OCT; VirtualLab Fusion; LED; synthetic light source

2096-2835(2016)04-0452-06

10.3969/j.issn.2096-2835.2016.04.017

2016-07-05 《中國(guó)計(jì)量大學(xué)學(xué)報(bào)》網(wǎng)址：zgjl.cbpt.cnki.net

國(guó)家自然科學(xué)基金資助項(xiàng)目(No.11404312)，國(guó)家留學(xué)基金資助項(xiàng)目(No.201408330449)，浙江省“儀器科學(xué)與技術(shù)”重中之重學(xué)科開放基金資助項(xiàng)目(No.JL150508)，廣西高校光電信息處理重點(diǎn)實(shí)驗(yàn)開放基金資助項(xiàng)目(No.KFJJ2014-03)，廣西自動(dòng)檢測(cè)技術(shù)與儀器重點(diǎn)實(shí)驗(yàn)開放基金資助項(xiàng)目(No.YQ15204).

TH741

A