李 丹，孫云飛

(1.吉林建筑科技學(xué)院，吉林 長春 130000；2.吉林師范大學(xué)信息技術(shù)學(xué)院，吉林 四平 136000)

1 引言

隨著第一臺激光器的誕生，相干光源逐漸進(jìn)入到大眾視野，由于其獨特的特性在工業(yè)生產(chǎn)、農(nóng)業(yè)發(fā)展、醫(yī)學(xué)發(fā)展、航天航空、軍事等領(lǐng)域，得到了迅速地發(fā)展以及廣泛應(yīng)用。尤其是在國防方面，激光武器超強(qiáng)的抗干擾能力、超高精度以及無可比擬的速度，是其它武器所不具備的。對激光技術(shù)的研究具有非常深遠(yuǎn)的意義。自相位共軛技術(shù)[1]是探究激光應(yīng)用過程中不可跨越的課題，自相位共軛現(xiàn)象在1972年被B.Ya.Zeldovich首次發(fā)現(xiàn)后，逐漸走進(jìn)相關(guān)專家學(xué)者的研究范疇內(nèi)，但大多數(shù)研究內(nèi)容為如何控制時間波形[2]、能量反射率的計算以及相位共軛保真度計算等方面，因硬件設(shè)備限制，無法精準(zhǔn)地獲得散射光束空間內(nèi)的光強(qiáng)分布[3]特性，而相位共軛技術(shù)中，光強(qiáng)分布能直接影響光束的質(zhì)量與能量，在多個應(yīng)用方面均起到了決定性的作用。

近幾年，在計算機(jī)網(wǎng)絡(luò)技術(shù)的幫助下，相關(guān)專家學(xué)者對光強(qiáng)分布的檢測展開研究，本文參考了多項課題內(nèi)容，提出一種散射光束空間內(nèi)光強(qiáng)分布檢測算法。首先，憑借CCD高速攝像機(jī)可靠性高、輸出效果好、耐振動[4]以及光譜響應(yīng)范圍寬等一系列優(yōu)點，精準(zhǔn)采集散射光束空間內(nèi)光斑圖像，由于采集到的原始圖像中含有較多的雜質(zhì)和干擾噪聲，為了得到清晰的檢測結(jié)果，對原始圖像進(jìn)行自適應(yīng)濾波處理操作[5]，得到光束光斑最真實的光強(qiáng)分布情況。然后，將散射光束空間內(nèi)的高斯光束分為若干條平行細(xì)光束[6]，通過分析這些細(xì)光束的傳輸方向和傳輸方式，建立唯像模型，得到各自的光強(qiáng)分布值；最后，將這些細(xì)光束的光強(qiáng)值結(jié)合在一起，得到總的光強(qiáng)分布值。仿真中，在入射光能量不同的環(huán)境下對光強(qiáng)分布進(jìn)行檢測，結(jié)果驗證了本文方法具有較高檢測精度和理想檢測結(jié)果。

2 CCD圖像采集與處理

本文利用CCD高速攝像機(jī)采集散射光束空間內(nèi)光斑圖樣。首先，使CCD攝像機(jī)的感光平面與光軸始終保持在垂直狀態(tài)下，保證在感光范圍內(nèi)可以采集到有效光斑圖像；然后，調(diào)節(jié)中性密度片，使進(jìn)入感光區(qū)域內(nèi)的光達(dá)到理想狀態(tài)，然后進(jìn)行多次、重復(fù)的圖像采集；最后，在所采集的圖樣中，挑選最佳圖樣數(shù)據(jù)，三維圖像如圖1(a)所示，二維圖像如圖1(b)所示。從圖中可以看出，大量的噪聲和雜質(zhì)影響了圖樣的質(zhì)量，無法實現(xiàn)精準(zhǔn)地檢測。本文借助MATLAB軟件中的圖像工具箱來消除噪聲。將中值濾波的窗口設(shè)置為3×4大小，利用其獨有的強(qiáng)濾除特性，實現(xiàn)對原始圖像的二維中值濾波處理。與此同時，結(jié)合自適應(yīng)濾波技術(shù)，共同作用于圖像中存在的噪聲和其它雜質(zhì)，完成預(yù)處理后的圖像如圖1(c)所示。

圖1 CCD攝像機(jī)圖像采集與處理

對比圖1中的三幅圖像，可以看出，原始圖像中的噪聲干擾基本被處理干凈，噪聲的平均值也逐漸下降到一個很低的水平。觀察圖1中的三幅圖像，可以發(fā)現(xiàn)其曲線變化規(guī)律與高斯光束曲線變化規(guī)律大致相同，因此，本文將散射空間內(nèi)的激光束看作是一組高斯光束。

3 散射光束空間光強(qiáng)分布檢測方法

為了實現(xiàn)光強(qiáng)分布的精準(zhǔn)檢測，本文首先對散射空間內(nèi)的高斯光束進(jìn)行分解處理[7]，得到多條平行細(xì)光束，通過計算平行細(xì)光束在傳輸過程中直射光和散射光的光強(qiáng)值，得到每條細(xì)光束距離z處的光強(qiáng)分布情況。再將每一條細(xì)光束的光強(qiáng)分布值相加，所得結(jié)果即為散射空間內(nèi)光強(qiáng)的分布檢測結(jié)果。

3.1 構(gòu)建平行細(xì)光束傳輸過程唯像模型

圖2為高斯光束劃分后的平行細(xì)光束傳輸示意圖。其中，Z軸代表細(xì)光束的傳輸方向，原點O表示散射空間內(nèi)光束的發(fā)出地。以Z軸和原點O為基礎(chǔ)，構(gòu)建一個柱坐標(biāo)。當(dāng)光束傳輸?shù)骄嚯xz處時，在橫截面XO′Y上隨機(jī)選取一點M(其坐標(biāo)為(r，φ，z)，其中，φ表示線段O′M與X軸所形成的角的大小)在橫截面上，光斑的光強(qiáng)通常是以圓對稱的形式分布的，因此，可以忽略φ的影響，直接對分布在X軸上的光強(qiáng)進(jìn)行檢測即可。

圖2 平行細(xì)光束傳輸示意圖

在平行細(xì)光束的傳輸過程中，由于傳輸介質(zhì)和傳輸環(huán)境的不同，使得細(xì)光束的光強(qiáng)分為兩種不同的類型：一種是直射光強(qiáng)，該種類型的光強(qiáng)表示光束自發(fā)出后一直到距離z處始終保持一種直射狀態(tài)，表達(dá)示如式(1)所示

Id(z)=I0(0)e-cz

(1)

式中，I0(0)表示光束剛發(fā)出時的光強(qiáng)，滿足條件I0(0)=1；c表示的是空氣衰減系數(shù)值。

另一部分為散射光強(qiáng)，表示光束發(fā)出后，經(jīng)過傳輸介質(zhì)以及傳輸環(huán)境的散射作用，使在固定傳輸空間內(nèi)的光束傳輸方向發(fā)生了改變[8]，導(dǎo)致最終的光強(qiáng)作用在傳輸空間以外。綜上所述，利用式(2)計算得到散射角θp以內(nèi)的散射光通量值

Φs(θp，z)=Fs(θp，z)(Φ(π，z)-e-cz)

(2)

式中，Φs(θp，z)表示當(dāng)傳輸距離為z時，角θp以內(nèi)的散射光通量值；Fs(θp，z)表示Φs(θp，z)在總的散射光通量(Φ(π，z)-e-cz)中所占的比例大小；Φ(π，z)表示在距離z處，直射光強(qiáng)與散射光強(qiáng)的光通量總值。

通過式(2)，可以計算得到θp和角(θp+dθ)以內(nèi)的散射光通量值[9]，再與圖2相結(jié)合，可以得到當(dāng)傳輸距離為z時、平行細(xì)光束與Z軸所形成的夾角在方向P上的散射光強(qiáng)，計算公式如式(3)所示

(3)

3.2 散射光束空間光強(qiáng)分布的數(shù)學(xué)計算模型

當(dāng)只考慮平行細(xì)光束在傳輸過程中的光強(qiáng)分布特性[10]時，可將θ0看作是遠(yuǎn)場發(fā)散角(也就是半角)，那么在傳輸距離為z、傳輸方向為θi時，點(xi，z)處的光強(qiáng)計算公式如式(4)所示

(4)

根據(jù)上文的計算與分析，在橫截面[11]上任意取一點Q，坐標(biāo)記為(xQ，z)。通過式(1)和式(4)可計算得到點Q處的直射光強(qiáng)為

Id(θQ，z)=I0(θQ，z)exp(-cz/cosθQ)

(5)

式中，Id(θQ，z)代表的是當(dāng)平行細(xì)光束的初始光強(qiáng)為θQ時、在θQ方向上傳輸了z/cosθQ的距離后，匯聚在點Q處的直射光強(qiáng)分布值。當(dāng)c的值為0時，Id(θQ，z)的值可以通過式(4)所示的光強(qiáng)計算公式計算得到[12]。

通過計算式(3)，可以得到在傳輸方向θij≠θQ上、初始光強(qiáng)為I0(θij)的細(xì)光束在點Q處的散射光強(qiáng)大小值為Is(|θij-θQ|，z)，其中，|θij-θQ|代表的是方向θij與方向θQ之間所形成的散射角大小。

綜上所述，將平行細(xì)光束在傳輸方向θij≠θQ上、點Q處的散射光強(qiáng)綜合進(jìn)行相加計算，即可得到細(xì)光束的總散射光強(qiáng)，計算公式如式(6)所示

(6)

那么散射光束空間內(nèi)，高斯光束在點Q處的總光強(qiáng)分布情況計算公式為

It(θQ，z)=Id(θQ，z)+Is(θQ，z)

θij≠θQθQ≤θ

(7)

通過上述公式計算得到散射光束空間光強(qiáng)分布結(jié)果，實現(xiàn)散射光束空間光強(qiáng)分布檢測。

4 仿真研究

在提出了光強(qiáng)分布檢測算法后，需要進(jìn)行相關(guān)的仿真來驗證方法在實際應(yīng)用中是否合理有效。因此，在實驗室中架構(gòu)檢測裝置，展開仿真驗證。

4.1 仿真裝置

考慮到激光獨有的特性：可以將光束匯聚在某處，形成非常小的匯聚點，得到一個類似理想點光源的入射光；再加上激光超強(qiáng)的相干性，但是極易受到光學(xué)元件、空氣中的灰塵以及其它細(xì)微雜質(zhì)的影響。因此，本文設(shè)計的實驗裝置圖如圖3所示。

圖3 實驗裝置圖

觀察圖3，調(diào)節(jié)Q點的裝置由Nd：YAG激光器來實現(xiàn)，它的主要作用是負(fù)責(zé)輸出基膜的偏振光。輸出的偏振光光束波長在1064nm上下波動，最高能量可達(dá)到40mJ，重復(fù)振動的頻率幅值為1Hz，入射脈沖的寬度高達(dá)7ns。為了防止Stokes光產(chǎn)生的散射光反射到激光器，在CCD1與CCD2之間、ED1與ED2之間還設(shè)置了隔離器。反射鏡的主要作用是反射光束，使光束可以按照規(guī)定進(jìn)行到能量計ED1和ED2中，以此完成對入射光能量的檢測。裝置還配備了聚焦透鏡，使反射光可以很好地聚焦在布里淵介質(zhì)池內(nèi)，介質(zhì)池選擇的是FC-72，介質(zhì)池的長度控制在60cm以內(nèi)。

裝置中CCD1與CCD2的主要作用是用來實時測量入射光的能量變化情況以及獲取入射光的光斑圖像。透射光光斑則由CCD3實現(xiàn)測量。在三臺CCD高速攝像機(jī)前方，都會適當(dāng)放置衰減片，以防過強(qiáng)的激光損壞攝像設(shè)備。為此，實驗選取的高速攝像機(jī)為MTV-1881EX型號，獲取到的圖像通過連接器傳送至計算機(jī)。CCD攝像機(jī)在整個實驗過程中為固定機(jī)位，進(jìn)入光能量的不同由衰減片來控制，CCD攝像機(jī)則在不斷變化的環(huán)境下實現(xiàn)光束光斑圖像的采集。為了實驗結(jié)果更具有說服力，在不同的入射光能量下，均進(jìn)行10次圖像的采集，綜合對比之下再進(jìn)行后續(xù)的計算和檢測。

4.2 實驗測試

圖4為利用CCD攝像機(jī)采集到的散射光束空間內(nèi)光斑圖像，當(dāng)入射能量處在不斷變化的情況下，通過調(diào)節(jié)衰減片的衰減率，保證CCD光強(qiáng)不會陷入飽和狀態(tài)，因此，圖像中光強(qiáng)之間不可進(jìn)行比較。

圖4 入射能量不同下的光斑圖像

利用本文方法對上述6幅圖像中的任意一幅進(jìn)行入射光斑采集以及圖像處理測試，測試結(jié)果如圖5所示。

圖5 圖像采集與處理測試結(jié)果

圖5為對10次圖像采集結(jié)果處理后的光強(qiáng)分布曲線圖。為了驗證本文方法對光強(qiáng)分布檢測的相對誤差，以10次結(jié)果作為一個統(tǒng)計量。通過觀察圖5可以看出，10條曲線的波動幅度相差較小，沒有出現(xiàn)較為明顯的曲線波動。這說明本文方法的圖像采集與處理效果是非?？煽康模瑸楹罄m(xù)進(jìn)行光強(qiáng)分布的檢測做好了充分的準(zhǔn)備工作。

接下來利用本文方法對經(jīng)過處理后的光斑圖像進(jìn)行光強(qiáng)分布的檢測，結(jié)果如圖6所示。

圖6 不同入射光下的光強(qiáng)分布檢測結(jié)果

從圖6中可以看出，利用本文方法檢測到的光強(qiáng)分布曲線對比原始圖像光強(qiáng)分布曲線，二者之間有著相同的變化幅度，在多個地方都出現(xiàn)了重合的現(xiàn)象，擬合程度高，這說明利用本文方法得到的光強(qiáng)分布檢測結(jié)果，具有較高的準(zhǔn)確性，檢測結(jié)果的可靠性較高。

5 結(jié)論

由于傳統(tǒng)方法在檢測光強(qiáng)分布時，無法根據(jù)激光的特性進(jìn)行精度檢測，存在較大的誤差，因此，本文提出了濾波去噪和高斯光束劃分下散射光束空間光強(qiáng)分布的檢測方法。首先，利用CCD高速攝像機(jī)獲取光束的光斑圖像，對其進(jìn)行預(yù)處理，得到其光斑變化規(guī)律；然后，將散射光束空間內(nèi)的高斯光束分為若干條平行細(xì)光束，再對每一條細(xì)光束進(jìn)行光強(qiáng)分布的計算，得到單條細(xì)光束的光強(qiáng)分布值；最后，將這些光強(qiáng)分布值結(jié)合在一起，得到最終的光強(qiáng)分布檢測結(jié)果。通過仿真，結(jié)果驗證了本文方法具有理想的光強(qiáng)分布檢測結(jié)果。為激光技術(shù)的繼續(xù)深入研究，提供了一種科學(xué)而又準(zhǔn)確的檢測方法。