基于NProd相關(guān)系數(shù)和人工免疫的激光圖像匹配算法

強鉦捷，張合新，孟 飛，張騰飛

（第二炮兵工程大學(xué)控制科學(xué)與工程系，西安 710025）

0 引言

現(xiàn)代制導(dǎo)技術(shù)對制導(dǎo)的性能要求越來越高，計算速度要求越來越快。對于同源激光圖像，其匹配基準(zhǔn)圖像和實時匹配圖像的成像機理和灰度特征相近，因此可采用基于灰度匹配方法。這種匹配方法制導(dǎo)精度高，但抗干擾性和實時性較差[1]。為此，本文針對激光圖像的同源灰度匹配方法存在的不足，提出一種NProd相關(guān)系數(shù)和人工免疫遺傳算法相結(jié)合的灰度配準(zhǔn)匹配方法。應(yīng)用NProd相關(guān)系數(shù)作為相似性度量，補償干擾的影響；利用人工免疫算法的魯棒性和強大并行搜索能力提高匹配速度和匹配穩(wěn)定性。

1 NProd相關(guān)系數(shù)及相關(guān)曲面

1.1 NProd系數(shù)

采用NProd相關(guān)系數(shù)[2]作為算法相似性度量，即計算圖像中每一點的NProd相關(guān)系數(shù)值，并將所有的NProd系數(shù)構(gòu)成相關(guān)曲面。

NProd系數(shù)的數(shù)學(xué)表達(dá)如下：

設(shè)F={Xi,j}為m×m的模板圖，G={Yi,j}為n×n的基準(zhǔn)圖，且m＜n，則點（u,v）的NProd相關(guān)系數(shù)為[3]

式中RNProd就是模板在點（u,v）處的NProd系數(shù)，且當(dāng)RNProd=1時，相似性度量值取得極大值，表明該點的匹配程度最高。

1.2 相關(guān)曲面特征分析

相關(guān)曲面以匹配坐標(biāo)點（x,y）作為平面坐標(biāo)，以NProd相關(guān)系數(shù)值作為縱坐標(biāo)繪制而成，如圖1所示。

圖1 相關(guān)曲面示意圖Fig.1 Schematic diagram of correlation surface

在圖1中，最高峰的位置坐標(biāo)代表圖像匹配度的最佳位置，次高峰是與最高峰最為接近的匹配點。對于激光主動成像圖像，往往存在噪聲干擾造成噪聲最高峰下降、次高峰上升，使得圖像存在一個或多個次高峰的情況，這說明在圖像中有一個或多個相似的匹配點，會對圖像匹配的精度產(chǎn)生影響。所以對于這種情況，應(yīng)首先進(jìn)行圖像去噪處理，增強圖像的峰值特征，其次再求取多峰值函數(shù)的最優(yōu)解。而人工免疫算法的高度并行性在處理多峰值求解問題時具有收斂速度快，魯棒性強等優(yōu)勢。因此，選擇人工免疫算法作為匹配算法的搜索策略。

2 人工免疫算法

2.1 免疫系統(tǒng)基本原理

人體的免疫可分為非特異性免疫和特異性免疫。非特異性免疫指的是人體本身所固有的免疫能力，它在首次遇到特定抗原時就能夠?qū)⑵湎麥?，是人體的第一道防線，通常人們也將非特異性免疫稱為先天性免疫或者固有免疫；特異性免疫指的是，當(dāng)非特異性免疫無法完全消滅入侵的抗原時，T淋巴細(xì)胞與B淋巴細(xì)胞將被激活，并通過一系列的分化過程形成效應(yīng)細(xì)胞，實現(xiàn)機體的免疫功能。通常人們也將特異性免疫稱作獲得性免疫或者適應(yīng)性免疫。特異性免疫通?？梢苑譃閮刹剑菏紫葹槌醮螒?yīng)答，主要是完成學(xué)習(xí)、適應(yīng)、記憶抗原特征的工作。其次是二次應(yīng)答。當(dāng)機體再次受到同種抗原入侵時，二次應(yīng)答可快速針對抗原產(chǎn)生對應(yīng)抗體消滅抗原[7]。所以特異性免疫具有適應(yīng)性、魯棒性和自學(xué)習(xí)的特點，是本文研究的重點。

2.2 匹配算法設(shè)計

STEP1：問題識別

將激光雷達(dá)圖像的匹配問題轉(zhuǎn)化成人工免疫算法求解問題，就是確定抗體和抗原。一般可將問題的目標(biāo)函數(shù)作為抗原，問題的解看做抗體。因此，匹配問題用人工遺傳算法的思想來描述就是：求和抗原最佳匹配的抗體。結(jié)合激光圖像目標(biāo)識別問題，將待匹配的實時圖作為目標(biāo)函數(shù)，即抗原；將基準(zhǔn)圖上最佳匹配的點作為問題的解，即抗體。

STEP2：產(chǎn)生初始抗體

以二進(jìn)制碼的形式來表示圖像的灰度值，并在解空間的按照如下規(guī)則產(chǎn)生初始抗體群：

1）設(shè)定種群大小

2）將模板窗口分成若干個區(qū)域

3）在每個區(qū)域中任意選取若干點，并計算每點的NProd值

4）選擇親和度（NProd值）大的點作為初始種群。

STEP3：計算親和度

采用式（1）來計算匹配點的NProd相關(guān)系數(shù)值作為抗體的親和度值。親和度越大則匹配匹配程度越高，當(dāng)NProd的值為最大值1時，證明結(jié)果是最為理想的狀態(tài)。

STEP4：提取疫苗

為了使抗體種群能夠朝著最優(yōu)方向進(jìn)化，我們將上一代種群pi中的最優(yōu)個體保留作為抗體疫苗amax，并將其應(yīng)用于pi+1代的種群中作為參考標(biāo)準(zhǔn)來修正pi+1代種群中的抗體。

STEP5：建立記憶庫

為了避免重復(fù)計算，可以記錄抗體的親和度值建立記憶庫。這樣當(dāng)再次計算該抗體親和度時就可以直接調(diào)用了。另外，記憶庫還具有記憶各次抽取疫苗的功能，也就是記錄每次計算種群中的最佳個體。

STEP6：抗體選擇

設(shè)種群中親和度最高的抗體為amax，則單個抗體ai與amax之間距離di表示二者相似度，di值越大表明抗體之間差異越明顯。

di表達(dá)式如下

若di值較為集中，則應(yīng)增大變異概率Pc保證抗體群多樣性；若di值比較分散，則減小變異概率Pc，保證種群進(jìn)化方向。

STEP7：抗體的抑制與促進(jìn)

對親和度低的個體進(jìn)行刪除，以保證抗體種群進(jìn)化程度；同時以一定比例克隆最佳抗體amax用以代替被刪除的抗體。

STEP8：抗體交叉

采用兩種交叉方式：

1）在種群內(nèi)部選取擁有較高親和度的抗體，并對其實施均勻交叉，以獲得最佳的抗體；

2）利用抗體和疫苗實施免疫試驗，利用疫苗來檢驗抗體，并對其進(jìn)行適當(dāng)?shù)男薷?。交叉方式見圖2。

圖2 抗體交叉Fig.2 Cross-reaction of antibodies

STEP9：抗體變異

按照設(shè)定的變異概率Pc進(jìn)行變異操作，過程如下圖3所示。

圖3 抗體變異Fig.3 Variation of antibodies

STEP10：更新記憶庫

STEP11：算法終止

設(shè)定閾值T，若上一代種群的平均NProd值與當(dāng)前種群的NProd值差異小于T時，且NProd值大于給定的系數(shù)CNProd，則可判定終止。否則返回STEP3重復(fù)操作。

匹配算法流程如圖4所示。

圖4 人工免疫算法流程圖Fig.4 Schematic diagram of artificial immune algorithm

3 仿真驗證

3.1 性能評價標(biāo)準(zhǔn)

為了全面檢驗算法在激光圖像目標(biāo)識別中的表現(xiàn)，本文通過匹配的準(zhǔn)確性、實時性和魯棒性三方面來檢驗。具體參考標(biāo)準(zhǔn)有：

1）準(zhǔn)確性指標(biāo)——匹配概率Pc

式中：nR表示正確匹配次數(shù)；nT表示匹配試驗的總次數(shù)。

2）實時性指標(biāo)——匹配時間TM

匹配時間采用多次匹配的平均時間值，定義為：

匹配的可靠性定義為式（3）。

式中：Ttotal為多次匹配的時間總和，nT為匹配的次數(shù)。

算法的魯棒性由匹配裕度RMM和匹配適應(yīng)度RMA綜合衡量。

由式（7）可得，RMA與Pc和RMM成正比。當(dāng)SSNRmax一定時，RMA的值越大，Pc就越高，也就是說算法擁有更強的畸變適應(yīng)性。

3.2 仿真結(jié)果

仿真實驗一：首先采用分辨率和信噪比較高的仿真激光雷達(dá)圖像進(jìn)行仿真實驗。通過將人工免疫算法（AIA）和均方誤差算法（MSE）、平均絕對差算法（MAD）、歸一化積相關(guān)算法（NProd）、標(biāo)準(zhǔn)遺傳算法（SGA）進(jìn)行對比，來檢驗本文算法在匹配概率、匹配時間和魯棒性上的效果。其中：圖5（a）為256×256的匹配基準(zhǔn)圖，圖5（b）直方圖為均衡化圖像，圖5（c）為匹配相關(guān)曲面圖。

式中：Std（?）表示求圖像的標(biāo)準(zhǔn)差；X、Y分別表示去均值后的基準(zhǔn)圖與實時圖。

匹配適應(yīng)度

圖5 仿真激光雷達(dá)圖像匹配Fig.5 Simulation LiDAR images matching

表1 匹配算法數(shù)據(jù)仿真結(jié)果Tab.1Arguments of simulation LiDAR images matching

由仿真結(jié)果可以看出，SGA的匹配時間最短，但由于SGA算法在受到噪聲干擾時容易陷入次優(yōu)解或漫游狀態(tài)，因此匹配概率也最低。MAD和MSE都是基于圖像灰度參數(shù)計算進(jìn)行比較匹配的，因此匹配概率高于SGA，但抗干擾性不理想，且二者匹配時間較長。NProd算法在匹配概率和魯棒性上都有比較理想的表現(xiàn)，但匹配時間過長。而本文算法采用NProd相關(guān)系數(shù)作為相似性度量，具有較好的匹配概率和魯棒性，并且充分利用人工免疫算法在并行運算搜索上的優(yōu)勢，彌補了算法在實時性上的不足，因此匹配效果最好。

此外，在仿真中，還對本算法和SGA算法在收斂性和收斂速度方面的表現(xiàn)進(jìn)行了對比，如圖6所示。

由圖6看出，SGA在進(jìn)行匹配時很早就陷入次優(yōu)解，進(jìn)入漫游狀態(tài)；而本文算法的種群是不斷朝向最優(yōu)方向進(jìn)化的，且在20代后得到最優(yōu)解，具有較快的收斂速度。

仿真實驗二：為進(jìn)一步驗證本文算法魯棒性和抗干擾能力方面的優(yōu)勢，實驗二采用信噪比較低的真實激光雷達(dá)圖像進(jìn)行匹配仿真。圖7（a）為256×256的匹配基準(zhǔn)圖，圖7（b）為直方圖均衡化圖像，圖7（c）為實時匹配窗口，圖7（d）為匹配相關(guān)曲面圖。

表2 激光雷達(dá)圖像匹配數(shù)據(jù)Tab.2 Arguments of LiDAR images matching

從表2可以看出，由于仿真二采用的真實激光雷達(dá)圖像信噪比相對較低，匹配算法受到不同程度的影響，在匹配概率與魯棒性兩項數(shù)據(jù)上都有不同幅度的下降，具體下降幅度參見表3。

表3 下降幅度比較Tab.3 Comparison of arguments descend range

通過對比表3中各算法在匹配精度和魯棒性上的下降幅度，可以得到結(jié)論：當(dāng)匹配圖像質(zhì)量下降和受到干擾的情況下，幾種典型的灰度匹配算法受到的干擾影響較大，在匹配概率上出現(xiàn)大幅度的下滑；SGA算法在強干擾下，過早地進(jìn)入了漫游狀態(tài)，匹配概率下降最大；而本算法在這兩方面的下降幅度最小，仍能保持較高的匹配性能，這進(jìn)一步證明了本算法具備較強的抗干擾能力與適應(yīng)性。

4 結(jié)論

本文采用具有較高準(zhǔn)確性和較強魯棒性的NProd相關(guān)系數(shù)作為匹配算法的相似性度量，并對人工免疫算法進(jìn)行了改進(jìn)，利用其強大的并行計算和良好收斂性作為搜索策略，保證了匹配算法的實時性，提出了一種基于NProd相關(guān)系數(shù)和人工免疫算法的激光圖像匹配算法。最后通過對激光仿真圖像和真實激光雷達(dá)圖像兩組仿真實驗，證明本文算法能夠結(jié)合NProd相關(guān)系數(shù)法和人工免疫算法優(yōu)點，具有較好的匹配精度，且算法魯棒性和實時性較強。

[1]杜文超,孟小芬,劉欽,等.基于SAR圖像的快速景象匹配方法[J].雷達(dá)科學(xué)與技術(shù),2014,（1）:39-43.

[2]Shi-qi Huang,Dai-zhi Liu.A novel method speckl enoise reduciton and ship target detection in SAR images[J].Pattern Recongniton,2009:1533-1542.

[3]劉毅飛,張旭明,丁明躍.歸一化互相關(guān)灰度圖像匹配的多核信號處理器實現(xiàn)[J].計算機應(yīng)用,2011,12:3334-3336.

[4]閔波.基于獨特型網(wǎng)絡(luò)的人工免疫模型研究與應(yīng)用[D].西安:西安電子科技大學(xué),2010.

[5]肖人彬,王磊.人工免疫系統(tǒng)：原理、模型、分析及展望[J].計算機學(xué)報,2002,25（12）:1281-1293.

[6]何珍梅,徐雪松.人工免疫系統(tǒng)研究綜述[J].華東交通大學(xué)學(xué)報,2007,04:79-83.

[7]金章贊,廖明宏,肖剛.否定選擇算法綜述[J].通信學(xué)報,2013,01:159-170.

[8]楊小岡.導(dǎo)彈多源圖像匹配定位方法與仿真技術(shù)研究[D].西安：第二炮兵工程學(xué)院博士論文,2006.