游騰飛，鐘舜聰，2*，鐘劍鋒，張秋坤，羅曼婷

(1.福州大學(xué)機(jī)械工程及自動(dòng)化學(xué)院，福建福州350108;2.福建省醫(yī)療器械和生物技術(shù)重點(diǎn)實(shí)驗(yàn)室，福建福州350000)

0 引言

自1991年David Huang 等人［1］提出光學(xué)相干層析技術(shù)(optical coherence tomography，OCT)以來(lái)，該技術(shù)就得到了迅速的發(fā)展。OCT 技術(shù)朝著更高分辨率、實(shí)時(shí)性、功能化和信息特異化方向發(fā)展，同時(shí)，OCT 系統(tǒng)朝著小型化方向發(fā)展，更好地滿(mǎn)足使用需求［2］。因?yàn)轭l域OCT 系統(tǒng)不需要參考臂進(jìn)行深度方向的掃描，相比時(shí)域OCT 系統(tǒng)大大縮短了成像時(shí)間，使得OCT系統(tǒng)實(shí)時(shí)成像成為可能［3-5］。為了再提高成像速度，頻域OCT 系統(tǒng)從一維頻域OCT 系統(tǒng)發(fā)展到二維的頻域OCT 系統(tǒng)［6-7］。因?yàn)槎S系統(tǒng)一次采集成像便可產(chǎn)生二維的樣品信息圖像，相比一維OCT 系統(tǒng)需要進(jìn)行多次采集來(lái)構(gòu)成相同的二維的樣品圖像而言，二維OCT系統(tǒng)在成像速度上大為提高，提高了OCT 系統(tǒng)成像的實(shí)時(shí)性。

頻域OCT 系統(tǒng)中包含CCD 相機(jī)、PZT 移相器、線性電動(dòng)移動(dòng)平臺(tái)等，為此，研究者需要開(kāi)發(fā)出能夠協(xié)調(diào)驅(qū)動(dòng)這些硬件的控制軟件，使得所搭建的OCT 系統(tǒng)形成有機(jī)整體，并對(duì)數(shù)據(jù)進(jìn)行處理和圖像重構(gòu)等工作。

VC++編程能力和代碼的執(zhí)行效率高，對(duì)硬件的控制能力強(qiáng)。但是在圖形處理和數(shù)據(jù)分析運(yùn)算等方面代碼復(fù)雜、編程難度大，不利于需要圖形處理和數(shù)據(jù)分析運(yùn)算的軟件的開(kāi)發(fā)和維護(hù)。Matlab 具有強(qiáng)大的數(shù)據(jù)處理分析能力和繪圖能力，使得其在算法研究等科學(xué)領(lǐng)域中廣泛運(yùn)用［8-12］。但是Matlab 代碼不能離開(kāi)其運(yùn)行環(huán)境而獨(dú)立運(yùn)行，并對(duì)底層的硬件控制能力和效率都不如VC++高。為此，在軟件開(kāi)發(fā)過(guò)程中，要是能結(jié)合VC++和Matlab 各自的優(yōu)點(diǎn)［13］，就可以降低軟件開(kāi)發(fā)難道，提高軟件的開(kāi)發(fā)效率。

本研究采用VC++與Matlab 混合編程技術(shù)開(kāi)發(fā)二維頻域光學(xué)相干層析成像系統(tǒng)。

1 光學(xué)系統(tǒng)的搭建

二維相干層析系統(tǒng)示意圖如圖1所示。由圖1 可見(jiàn)，和一維頻域光學(xué)相干層析系統(tǒng)的組成相類(lèi)似，二維相干層析系統(tǒng)可以分為:光源、樣品臂、參考臂、信號(hào)探測(cè)端等部分。不同之處在于平行光入射到分光鏡之前，會(huì)先經(jīng)過(guò)柱透鏡。經(jīng)過(guò)柱透鏡后，平行光會(huì)被匯聚成一條線。利用這條匯聚的光線照射到待測(cè)樣品表面，那么一次采集便能生成這條光線上待測(cè)樣品的結(jié)構(gòu)信息。相比一維的單點(diǎn)成像需要多次采集而拼構(gòu)成一條二維的待測(cè)樣品的結(jié)構(gòu)信息而言，二維的系統(tǒng)大大提高了成像速度。

圖1 二維相干層析系統(tǒng)示意圖

由圖1 可見(jiàn)，光源發(fā)出的寬帶弱相干光被透鏡聚焦后進(jìn)去光纖傳輸，而后從光纖輸出的光通過(guò)透鏡后準(zhǔn)直成一束平行光。該束平行光通過(guò)柱透鏡后聚焦成一條線，透過(guò)分光鏡后，被分成兩束。其中一束入射到PZT 參考鏡中并被參考鏡發(fā)射回來(lái)形成參考光;另一束入射到待測(cè)樣品并被待測(cè)樣品后向散射或者反射回來(lái)，形成樣品光。參考光和樣品光返回分光鏡后相遇，倘若參考臂和樣品臂光程在一個(gè)相干長(zhǎng)度范圍內(nèi)，那么兩束光將發(fā)生干涉，形成干涉光。該干涉光在被CCD 相機(jī)采集之前，還需要使用光柵將該干涉光按波長(zhǎng)分開(kāi)，最后被CCD 相機(jī)采集并上傳到電腦中進(jìn)行數(shù)據(jù)處理分析和圖像重構(gòu)。

系統(tǒng)中柱透鏡使用型號(hào)為T(mén)HORLABS 公司的LJ1821L1-B，焦距f=50 mm。線性移動(dòng)平臺(tái)使用THORLABS 公司型號(hào)為MTS50/M-Z8 的電動(dòng)平臺(tái)。該電動(dòng)平臺(tái)行程為50 mm，分辨率為29 nm。CCD 相機(jī)使用深圳迪美捷公司型號(hào)為DC140M 相機(jī)，該相機(jī)有效像素為1 360×1 024，在該分辨率下幀率為7.5 fps，像元尺寸為4.65 μm×4.65 μm。驅(qū)動(dòng)PZT 移相器的NI 卡采用NI 公司型號(hào)為NI USB-6009 的NI 卡。該NI 卡具有12 位數(shù)字I/O 輸出端，輸出電壓范圍為0～5 V。

2 系統(tǒng)控制軟件的開(kāi)發(fā)

在硬件系統(tǒng)搭建完畢后，還需要開(kāi)發(fā)出相應(yīng)的控制軟件系統(tǒng)，才能使整套系統(tǒng)正常工作，為后續(xù)對(duì)光學(xué)相干層析的研究提供便利。為獲得高效的代碼執(zhí)行效率、硬件控制能力，以及數(shù)據(jù)分析和圖像處理能力，本研究利用VC ++和Matlab 混合編程技術(shù)，充分利用軟件各自的優(yōu)點(diǎn)，開(kāi)發(fā)高效的控制軟件。

在VC++中調(diào)用Matlab 進(jìn)行混合編程比較常用的方法有:①利用Matlab 引擎;②利用Matlab 自帶的MCC 編譯器;③利用MATCOM 編譯器;④利用Matlab的COM Build 工具等方法。以上方法各具優(yōu)點(diǎn)和不足，根據(jù)實(shí)際需要，選擇合適的方法，才能取得事半功倍的效果。二維光學(xué)相干層析技術(shù)需要進(jìn)行二維甚至三維的顯示以及比較復(fù)雜的圖像處理，如果能夠充分利用Matlab 提供的圖像工具箱和圖像函數(shù)，那么對(duì)軟件的開(kāi)發(fā)將大大減小難度，代碼直觀，利于升級(jí)和維護(hù)。為此，選擇在VC ++中調(diào)用Matlab 引擎的方法相比其他方法具有明顯的優(yōu)勢(shì)。

基于VC ++和Matlab 開(kāi)發(fā)的軟件，所需驅(qū)動(dòng)的硬件設(shè)備包括:①NI 卡，用于控制PZT 移相器產(chǎn)生移相;②線性電動(dòng)移動(dòng)平臺(tái)，用于帶動(dòng)待測(cè)樣品進(jìn)行橫向掃描;③面陣CCD 相機(jī)，用于對(duì)干涉光進(jìn)行采集。這些設(shè)備產(chǎn)商都有提供SDK 工具包，便于根據(jù)需要進(jìn)行二次開(kāi)發(fā)。程序界面如圖2所示。

圖2 控制軟件運(yùn)行界面

如圖2所示，軟件主要包括干涉圖像顯示區(qū)、線性電動(dòng)移動(dòng)平臺(tái)控制界面以及一些輸入框和控制按鈕。干涉圖像顯示區(qū)主要是用于顯示CCD 相機(jī)采集到的干涉信號(hào)(如圖2 中左半部分黑白相間的干涉條文)。CCD 相機(jī)自身可以根據(jù)光強(qiáng)條件自動(dòng)調(diào)節(jié)曝光時(shí)間等參數(shù)，也可以根據(jù)需要人工對(duì)包括曝光時(shí)間、對(duì)比度等參數(shù)進(jìn)行設(shè)定。

電動(dòng)平臺(tái)控制界面用于對(duì)線性電動(dòng)移動(dòng)平臺(tái)的控制，達(dá)到對(duì)樣品進(jìn)行掃描的目的?？梢愿鶕?jù)需要設(shè)定電動(dòng)平臺(tái)的加速度、最大速度以及步進(jìn)距離等運(yùn)動(dòng)參數(shù)。

當(dāng)需要對(duì)待測(cè)樣品進(jìn)行實(shí)時(shí)成像的時(shí)候，本研究調(diào)節(jié)好光路，選擇成像所需的方法，單機(jī)“監(jiān)測(cè)數(shù)據(jù)”按鈕，程序便開(kāi)始驅(qū)動(dòng)CCD 相機(jī)進(jìn)行干涉信號(hào)的采集，驅(qū)動(dòng)PZT 移相器移相以及調(diào)用Matlab 引擎進(jìn)行數(shù)據(jù)處理分析和圖像顯示。

程序流程圖如圖3所示。

圖3(a)為一次實(shí)時(shí)成像的流程圖，在關(guān)閉實(shí)時(shí)成像之前，程序一直按照該流程進(jìn)行重復(fù)的圖像顯示。各硬件初始化完成后，單擊“監(jiān)測(cè)數(shù)據(jù)”按鈕，程序控制NI 卡輸出預(yù)設(shè)定的電壓(電壓程序預(yù)設(shè)定好，也可根據(jù)需要在界面輸入而修改默認(rèn)設(shè)定)，移相器接受電壓而伸長(zhǎng)(移相)，而后驅(qū)動(dòng)CCD 采集在該相位的干涉圖，并傳輸并存儲(chǔ)到計(jì)算機(jī)內(nèi)存中，完成一副干涉圖的采集。如果還需另一個(gè)相位的圖像，則NI 卡重新輸出另一電壓值，驅(qū)動(dòng)PZT 移相器產(chǎn)生另一伸長(zhǎng)量，CCD 相機(jī)再進(jìn)行一次數(shù)據(jù)采集，并存儲(chǔ)到計(jì)算機(jī)內(nèi)存中。如果不需其他相位的干涉圖，則程序調(diào)用Matlab引擎，對(duì)剛保存在計(jì)算機(jī)內(nèi)存的數(shù)據(jù)進(jìn)行處理分析，并完成繪圖進(jìn)行顯示。至此便完成一次實(shí)時(shí)成像。

圖3 控制流程圖

如果需要對(duì)圖像的數(shù)據(jù)進(jìn)行保存，方便日后使用和分析處理，也可以單擊“采集數(shù)據(jù)”按鈕進(jìn)行數(shù)據(jù)采集?？刂屏鞒虉D如圖3所示。采集數(shù)據(jù)之前，先設(shè)定好數(shù)據(jù)所保存的位置，電動(dòng)移動(dòng)平臺(tái)的步距以及所需采集的次數(shù)等。在完成各硬件初始化后，NI 卡輸出電壓，驅(qū)動(dòng)PZT 移相器產(chǎn)生移相，接著CCD 相機(jī)采集待測(cè)樣品某個(gè)位置在該相位下的干涉數(shù)據(jù)，保存到剛才設(shè)定的位置。如果還需其他相位的數(shù)據(jù)采集，則NI 卡再次輸出電壓，PZT 發(fā)生移相器發(fā)生相移，然后CCD 再次采集數(shù)據(jù)并保存數(shù)據(jù)。如果待測(cè)樣品該位置數(shù)據(jù)采集結(jié)束，則驅(qū)動(dòng)電動(dòng)移動(dòng)平臺(tái)運(yùn)動(dòng)一個(gè)步距的距離，重新按照上述再次進(jìn)行另外一個(gè)位置的數(shù)據(jù)的采集。

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

筆者在完成二維頻域光學(xué)相干層析系統(tǒng)的硬件搭建和軟件開(kāi)發(fā)調(diào)試后，進(jìn)行對(duì)樣品進(jìn)行成像，測(cè)試系統(tǒng)的可行性。所使用待測(cè)樣品為透明膠帶。該種透明膠帶由塑料薄膜層和膠水層組成。塑料薄膜層、膠水層以及空氣這三者之間它們的折射率都不相同，因此在光通過(guò)他們的分界面時(shí)，由于折射率的變化，光會(huì)產(chǎn)生反射現(xiàn)象，那么對(duì)該透明膠帶進(jìn)行成像，會(huì)分出兩層三個(gè)界面。

CCD 相機(jī)獲取的干涉信號(hào)除了攜帶樣品結(jié)構(gòu)信息的干涉條文，還包含很強(qiáng)的直流噪聲，以及自相干噪聲等。為了防止直流噪聲以及自相干噪聲等對(duì)最后的結(jié)果產(chǎn)生影響，筆者在對(duì)干涉信號(hào)進(jìn)行傅里葉變換前，需要先去除直流噪聲。

經(jīng)過(guò)對(duì)CCD 采集的干涉信號(hào)去除直流分量后，筆者對(duì)信號(hào)做傅里葉變換，最終便可以形成待測(cè)樣品的結(jié)構(gòu)信息。透明膠帶結(jié)構(gòu)信息圖如圖4所示。

圖4(a)左半部分(或右半部分)可以看出具有3條亮線，其中一條比較明亮，另外兩條比較弱。其中比較亮的一條為透明膠帶不黏的那一面的表面，中間那條表示的透明膠帶的塑料薄膜和膠水薄膜的分層界面，最左邊的一條為膠水和基體的分界面面。圖4(b)左半部分(或右半部分)也可以清晰地看出曲線具有3個(gè)峰，相應(yīng)地對(duì)應(yīng)不同的分層界面?？梢钥闯鲈摻Y(jié)果也是符合透明膠帶的實(shí)際情況的。

圖4 透明膠帶結(jié)構(gòu)信息圖

然而，如果僅僅去除直流分量后做傅里葉變換，傅里葉變換后將會(huì)產(chǎn)生共軛鏡像，如圖4所示，在圖像的左右各有一個(gè)以光程差為零的點(diǎn)為對(duì)稱(chēng)軸相對(duì)稱(chēng)的圖像。假如零光程差的位置在樣品內(nèi)部，則勢(shì)必造成樣品的結(jié)構(gòu)信息圖像與其共軛鏡像相重疊，無(wú)法分辨樣品結(jié)構(gòu)信息的問(wèn)題。

為解決該問(wèn)題方法有:①將零光程差的位置調(diào)節(jié)到樣品外部，這樣就不會(huì)造成混疊;②去除共軛鏡像。

本研究采用第二種方法，采用五步移相法，重新構(gòu)造復(fù)數(shù)函數(shù)［14］，再進(jìn)行傅里葉變換，從而去除共軛鏡像，去除共軛鏡像后透明膠帶結(jié)構(gòu)信息如圖5所示。圖5 中(a)圖和(b)圖中均只有圖4 中對(duì)應(yīng)的(a)圖和(b)圖的左半部分的樣品結(jié)構(gòu)信息了。

對(duì)比圖4 和圖5 可知，五步移相法不僅把共軛鏡像去除，而且還增強(qiáng)了圖像的信噪比。

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

基于VC++和Matlab 的二維頻域光學(xué)相干層析成像系統(tǒng)具有易于操作，便于維護(hù)和升級(jí)的特點(diǎn)。在系統(tǒng)測(cè)試過(guò)程中，筆者對(duì)一層透明膠帶進(jìn)行了成像，成像結(jié)果符合實(shí)際，證明了系統(tǒng)是可行的。

圖5 去除共軛鏡像后透明膠帶結(jié)構(gòu)信息

由于一次成像便能產(chǎn)生二維的樣品結(jié)構(gòu)信息，相比一維頻域光學(xué)相干層析成像系統(tǒng)，二維頻域光學(xué)相干層析具有快速成像的優(yōu)點(diǎn)，并且軟件集成度高，功能全面，界面友好，對(duì)光學(xué)相干層析系統(tǒng)的開(kāi)發(fā)具有促進(jìn)作用。

［1］HUANG D.，SWANSON E.，LIN C.et al.，Optical coherence tomography［J］.Science，1991，254(22):1178-1181.

［2］吳 凌.OCT 三維重建與多普勒成像研究［D］.浙江:浙江大學(xué)光電信息工程學(xué)系，2007.

［3］馬振鶴，王瑞康，張 帆，等.快速高分辨率的頻域光學(xué)相干層析成像系統(tǒng)研究［J］.納米技術(shù)與精密工程，2005，3(3):232-235.

［4］DE BOER J F，CENSE B，PARK B H，et al.Improved signal-to-noise ratio in spectral domain compared with time-domain optical coherence tomography［J］.Optics letters，2003，28(21):2067-2069.

［5］LEITGEB R A，HIZENBERGER C K，F(xiàn)ERCHER A F.Performance of fourier domain versus time domain optical coherence tomography［J］.Optics express，2003，11(8):889-894.

［6］WITTE S，BACLAYON M，PETERMAN E J G.et al.Single-shot two-dimensional full-range optical coherence tomography achieved by dispersion control［J］.Optics express，2009，14(17):11335-11349.

［7］ZHONG Jian-feng，ZHONG Shun-cong，ZHANG Qiu-kun.et al.Two-dimensional optical coherence tomography for real-time structural dynamical characterization［J］.Optics and Lasers in Engineering，2015(66):74-79.

［8］石成英，張利民，陳登科.基于VC++與Matlab 聯(lián)合編程的測(cè)量系統(tǒng)軟件開(kāi)發(fā)研究［J］.計(jì)算機(jī)數(shù)字工程，2013，41(9):1535-1537.

［9］李 卓，李益民.Matlab 與VC ++混合編程技術(shù)在圖像處理中的應(yīng)用研究［J］.電腦知識(shí)與技術(shù)，2011，7(22):5450-5452.

［10］沈莉麗.基于Matlab 與VC ++混合編程的數(shù)字信號(hào)處理系統(tǒng)［J］.信息化研究，2010，36(5):43-45.

［11］董建婷，王 杰.VC ++調(diào)用Matlab 實(shí)現(xiàn)圖像處理的3種常用方法比較［J］.現(xiàn)代電子技術(shù)，2012，35(16):46-49.

［12］許雪開(kāi).VC ++與Matlab 混合編程［J］.機(jī)電工程，2007，24(2):26-27.

［13］鐘舜聰，丁 康，孫智娟.Matlab 與C/C++混合編程在信號(hào)處理中的應(yīng)用［J］.微計(jì)算機(jī)應(yīng)用，2000(5):305

［14］何梓昂.FD-OCT 去復(fù)共軛成像的研究［D］.杭州:浙江大學(xué)光電信息工程學(xué)系，2008.