郭 軍，馬金鳳

（廣州海洋地質(zhì)調(diào)查局資料處理研究所，廣東 廣州 510760）

隨著多波束聲納、側掃聲納、淺層剖面儀等水下遙測工具的不斷出現(xiàn)，人們對海底底質(zhì)的認識也不斷清晰起來，借助水下遙測工具獲得的海底圖像，技術人員可辨識不同類別的海底質(zhì)地，同時也可發(fā)現(xiàn)海底的人工目標。

由于不同類型的底質(zhì)在聲圖上構成的紋理結構是不同的［1－2］。紋理是海底表面結構粗糙程度的直接反應，是屬于特定類型的底質(zhì)特有的屬性，能將其與不同類型的底質(zhì)區(qū)分開來。本文采用灰度共生矩陣來統(tǒng)計紋理，利用不同的權矩陣對其進行加權處理可得到紋理特征統(tǒng)計量。

近年來神經(jīng)網(wǎng)絡的廣泛應用為海底底質(zhì)分類提供了新的方法，如采用BP神經(jīng)網(wǎng)絡、Hopflied神經(jīng)網(wǎng)絡及小波神經(jīng)網(wǎng)絡等對底質(zhì)進行監(jiān)督分類［3］。然而在實際應用中，上述網(wǎng)絡還存在著分類精度依賴于網(wǎng)絡訓練樣本，且訓練速度慢、不易收斂等缺點。自組織競爭神經(jīng)網(wǎng)絡是一種“無師自通”的網(wǎng)絡，具有自組織、自適應的學習能力。本文采用其對海底底質(zhì)進行自動分類。

1 特征矢量提取

底質(zhì)分類中采用的特征矢量有3種：直接的反向散射強度、對反向散射強度譜分析后提取的參數(shù)、影像紋理分析［4－5］。實際應用中，由于反向散射強度模型與測量值之間存在著一定的偏差，可能導致不同聲納系統(tǒng)獲得的結果不一致［6］。若只采用前兩種矢量，可能會影響分類精度。而影像紋理不直接采用像素絕對大小，其作為特征矢量更具有普遍性，本文將其作為底質(zhì)分類的輸入特征矢量。

1.1 共生矩陣

灰度共生矩陣反映圖像灰度關于方向、相鄰間隔及變換幅度的綜合信息［7］。本文采用反差、方差、熵、逆差距、角二階距和灰度相關6個紋理統(tǒng)計特征，公式為

式中：ux，uy，σx，σy分別為px，py的均值和均方差；u為均值；

1.2 K－L變換

設有向量集

其數(shù)學期望為E（X），U是X的協(xié)方差矩陣C的特征向量按其特征根由大到小的順序排列而構成的變換矩陣。那么

為主成分分析算法。

設λi為變換矩陣U的特征值，將其按λ1≥λ2≥…λn排列，令為前d個主成分的累積貢獻率，反映了前d個主成分綜合原始7個特征矢量信息的能力。究竟采用幾個主成分來代替原始的7個特征矢量才合適呢？通常的做法是取較小的d，使得前d個主成份的累積貢獻率不低于某一個水平（通常為85%）。本文采用V作為對聲圖特征信息進行主成分分析后保留多少主分量的衡量準則。

2 自組織競爭神經(jīng)網(wǎng)絡

自組織神經(jīng)網(wǎng)絡是一種沒有導師監(jiān)督學習方式的神經(jīng)網(wǎng)絡［8］，該網(wǎng)絡具有自組織、自適應改變網(wǎng)絡參數(shù)與結構的特點，能夠自行分析、比較樣本的內(nèi)在規(guī)律，并對具有共同特征的樣本進行正確的分類，其網(wǎng)絡結構如圖1所示。

圖1 自組織競爭神經(jīng)網(wǎng)絡結構

本文采用的自組織競爭神經(jīng)網(wǎng)絡具體設計步驟如下：

1）初始化網(wǎng)絡。設置網(wǎng)絡的連接權值wij（i＝1，2，…，N，j＝1，2，…，M），其中N為輸入層神經(jīng)元個數(shù)，M為競爭層神經(jīng)元個數(shù)，且滿足約束條件

2）給輸入層賦初值，作數(shù)據(jù)的歸一化。

3）計算競爭層神經(jīng)元狀態(tài)，xi輸入樣本向量的第i個元素，則有

4）根據(jù)競爭機制，具有最大加權值的神經(jīng)元k贏得競爭勝利輸出為

5）修正競爭后的權值，對所有的輸入層神經(jīng)元i有

6）若各連接權值的調(diào)整量趨近于0，則結束，否則返回步驟2）。

3 實驗分析

實驗用圖來源于加利福尼亞大學海圖實驗室提供的位于加利福尼北部Mackerricher州立保護區(qū)內(nèi)的側掃聲納圖像。根據(jù)已有的數(shù)據(jù)顯示，該地區(qū)的海底底質(zhì)大致分為5類：泥土、沙地、巖石、沙礫、泥沙。實驗中，在上述5個不同的特征區(qū)域提取5個樣本圖像數(shù)據(jù)，如圖2所示，定義樣本的尺寸大小為128×128像素，將每個樣本圖像分成16小塊，尺寸大小為32×32像素，其中8塊用于網(wǎng)絡的訓練和學習，剩余8塊用于網(wǎng)絡的測試分類。原始樣本圖像的灰度級為256級，由此得出的共生矩陣過大，增加了運算量。在不影響紋理的情況下，對樣本圖像的灰度級做適當?shù)淖儞Q以減少灰度級，再計算共生矩陣。本文將樣本圖像的灰度級降至16級，計算樣本的共生矩陣后，對提取出的特征矢量和均值進行K－L變換已達到降維和減少冗余信息的目的，然后對樣本圖像構建自組織競爭神經(jīng)網(wǎng)路，在該網(wǎng)絡中對樣本圖像實現(xiàn)自動分類。

圖2 海底底質(zhì)樣本

實驗中計算出子樣本的均值、反差、方差、熵、逆差距、角二階距和灰度相關共有7個特征統(tǒng)計量，計算結果見表1、表2。

表1 沙地樣本特征統(tǒng)計量

表2 泥沙樣本特征統(tǒng)計量

實驗以沙地為例，顯示了4個樣本的數(shù)據(jù)，計算結果如圖3所示。從圖3中可以看出，前4個主成分的累計貢獻率為94.80%，遠遠超過了85%，這說明前4個主成分完全可以將原有特征矢量的主要信息反映出來。為此本文將特征矢量經(jīng)過K－L變換后的維數(shù)設置為4，即每個樣本數(shù)據(jù)可以用4個主要的特征矢量表示。沙地和泥沙經(jīng)過K－L變換后的特征矢量見表3、表4。

圖3 沙地主成分累積貢獻率

表3 沙地樣本主成分分析特征向量

表4 泥沙樣本主成分分析特征向量

實驗共有16×5＝80組樣本數(shù)據(jù)，其中前40組作為網(wǎng)絡訓練集數(shù)據(jù)，后40組作為測試集數(shù)據(jù)。采用一個80×4的矩陣，即每個樣本中包含4個元素，由于要區(qū)分的類別數(shù)目為5，因此設置自組織競爭神經(jīng)網(wǎng)絡的競爭層神經(jīng)元數(shù)目為5。為了加快學習速度，將學習速率設置為0.1，訓練次數(shù)為500次。實驗結果見表5、表6、表7。從實驗結果可以看出，自組織競爭神經(jīng)網(wǎng)絡成功地對40個訓練樣本進行聚類，其中分類標識為1的樣本為沙石，標識為2的樣本為泥土，標識為3的樣本為巖石，標識為4的樣本為沙地，標識為5的樣本為泥沙。從表7中可以看出，經(jīng)過訓練的自組織競爭神經(jīng)網(wǎng)絡對40個測試樣本進行了分類，分類錯誤率為1/40＝0.025%，將一個沙地樣本錯歸為泥沙樣本，其他各類樣本均沒有出現(xiàn)錯分現(xiàn)象。

表5 訓練樣本的分類結果

表6 測試樣本的分類結果

表7 類矩陣及精度評價

為了驗證本文算法的有效性，采用傳統(tǒng)的最大似然法和基于BP神經(jīng)網(wǎng)絡法與之進行比較，實驗為：本文算法精度為98%，BP神經(jīng)網(wǎng)絡法為89%，最大似然法為90%。由結果可以看出，本文算法精度要明顯高于傳統(tǒng)的分類算法，錯分現(xiàn)象顯著減少，驗證了該算法在海底質(zhì)地分類中的有效性和可行性。最后筆者利用matlab平臺編寫了基于K－L變換的自組織神經(jīng)網(wǎng)絡海底底質(zhì)分類軟件，包括圖像預處理、特征提取、K－L變換和網(wǎng)絡分類4個模塊。

4 結束語

采用基于K－L變換的自組織競爭神經(jīng)網(wǎng)絡分類，利用了該網(wǎng)絡自組織學習特點，K－L變換大大減少了聲圖特征矢量的冗余性，同時又保持了各分量之間的獨立性?；叶裙采仃嚭突叶染党浞址从沉瞬煌５椎踪|(zhì)的特性，由于其存在交叉，采用單一的特征矢量分類不可行，將兩者相結合，可獲得正確的分類結果。

［1］Chakraborty B，Kodagali V，Baracho J.Sea－Floor Classification Using Multibeam Echo－Sounding Angular Backscatter Data：A Real－Time Approach Employing Hybrid Neural Network Architecture［J］.IEEE Journal of Oceanic Engineering，2003，28（1）：121－128.

［2］Preston J M，Christney A C，Bloomer S F，et al.Seabed Classification of Mulibeam Sonar Image［J］.Oceans，2001（4）：2616－2623.

［3］Michalopoulou Z H，Alexandrou D，Moustier D C.Application of neural and statistical classifiers to the problem of seafloor characterization［J］.IEEE Journal of Oceanic Engineering，1995，20（3）：190－197.

［4］Tamseet D.Seabed Characterization and Classification from Power Spectura of Side－scan Sonar Data［J］.Geophys，1993，15：43－64.

［5］陽凡林.多波束和側掃聲納數(shù)據(jù)融合及其在海底底質(zhì)分類中的應用［D］.武漢：武漢大學，2003.

［6］唐秋華，周興華，劉寶華，等.學習向量量化神經(jīng)網(wǎng)絡在多波束底質(zhì)分類中的應用［J］.武漢大學學報：信息科學版，2006，31（3）：44－47.

［7］田艷琴，郭 平，盧漢清.基于灰度共生矩陣的多波段遙感圖像紋理特征的提取［J］.計算機科學，2004（12）：164－165.

［8］田景文，高美娟.人工神經(jīng)網(wǎng)絡算法研究及應用［M］.北京：北京理工大學出版社，2006.