譚克，鄒麗娜（杭州應用聲學研究所聲納技術(shù)國防科技重點實驗室，浙江杭州310023）

?

正態(tài)分布轉(zhuǎn)換算法在聲納圖像處理中的應用

譚克，鄒麗娜
（杭州應用聲學研究所聲納技術(shù)國防科技重點實驗室，浙江杭州310023）

摘要：根據(jù)聲納圖像匹配處理的實時性需求，研究了正態(tài)分布轉(zhuǎn)換算法。該方法基于對圖像中目標點的概率分布進行匹配，不需要提取圖像的幾何特征或灰度特征，可有效縮短匹配時間。通過對均值濾波法、中值濾波法、小波去噪法、形態(tài)學平滑法4種圖像去噪和迭代閾值法、最大熵法、一維Otsu法和二維Otsu法4種分割算法，進行仿真并對比，選出符合實時性處理的方法，即中值濾波法和一維Otsui法?；诿绹轮輧x器公司的數(shù)字信號處理芯片，將圖像去噪、分割和匹配算法在硬件上實現(xiàn)，硬件實現(xiàn)結(jié)果和仿真結(jié)果一致，實現(xiàn)了聲納圖像匹配的實時處理。

關(guān)鍵詞：聲學；圖像匹配；數(shù)字信號處理器；聲納圖像處理；正態(tài)分布轉(zhuǎn)換算法

0　引言

聲納圖像處理在水下信號研究領(lǐng)域具有重要地位，在經(jīng)濟、軍事上有著重要的研究價值。圖像匹配技術(shù)是圖像處理技術(shù)中的重要組成部分，常常作為其他圖像處理應用的前處理步驟來使用。最近十幾年里，國內(nèi)外的許多學者提出了各種各樣面對許多不同領(lǐng)域、不同要求的匹配算法，力求在各個領(lǐng)域中，加快圖像匹配算法的速度，提高匹配精度，而且許多算法都具有良好的通用性。如David［1］提出基于尺度不變特征的圖像匹配方法，Anuta［2］提出了快速傅里葉變換互相關(guān)圖像匹配方法，高晶等［3］提出新的基于角點檢測的匹配算法，韓冰等［4］提出基于一種快速歸一化積相關(guān)算法，劉佳等［5］提出了一種基于改進尺度不變特征變換算法的圖像匹配方法，等等。

圖像匹配技術(shù)的應用越來越廣泛，如飛行器導航、武器的末端制導、指紋識別等。而這種技術(shù)在聲納圖像處理上的應用還不多。本文選用Biber等［6］提出的正態(tài)分布轉(zhuǎn)換（NDT）算法，結(jié)合匹配步驟前的一些預處理算法，實現(xiàn)聲納圖像去噪、分割和匹配這一過程。選用蛙人圖像，對4種預處理算法進行分析對比，選出了符合實時性要求的預處理算法。使用NDT算法對預處理后的圖像進行匹配處理，并將該處理流程在硬件上實現(xiàn)，可用來實時地匹配連續(xù)兩幅蛙人圖像，為后續(xù)蛙人圖像識別和圖像融合奠定了基礎(chǔ)。

1　NDT算法基本原理

NDT算法已經(jīng)在移動機器人定位中有一定的應用，但在水聲領(lǐng)域中應用較少。NDT算法的基本思想是：在匹配過程中無需查找點或者其他幾何特征之間的對應關(guān)系，把二維圖像中的像素點集合轉(zhuǎn)換成為分段連續(xù)可微的概率分布，對基準圖像和待匹配圖像的概率分布進行匹配［7］，可以快速地完成圖像匹配。

用連續(xù)可微的概率分布對二維的圖像進行表示，其步驟如下：

1）對二維圖像進行區(qū)域劃分，按照一定的大小規(guī)則，分成相同的小格。

2）對含有像素點的小格，以概率密度的形式進行描述，通過正態(tài)分布概率密度函數(shù)來表示柵格中每一個像素點，如（1）式～（3）式所示。

式中：q為每個小柵格中包含像素點樣本的坐標平均值所構(gòu)成的矩陣；∑為協(xié)方差矩陣；r為像素點坐標；ri為第i個像素點的坐標。

成像聲納所拍攝的兩幅不同圖像坐標系的關(guān)系為

式中：（tx，ty）T是兩幅圖像坐標系的相對平移量；φ是圖像間的旋轉(zhuǎn)角度。

基于NDT的匹配算法步驟如下：

1）選一幀圖像作為基準圖像，建立該圖像數(shù)據(jù)的NDT；

2）給坐標變換參數(shù)設(shè)置一個值，也就是對其進行初始化（可以用“0”作為初始坐標參數(shù)）；

3）根據(jù)坐標參數(shù)，將待匹配圖像中的像素點映射到基準圖像中去，用基準圖像的坐標來對待匹配圖像中的像素點進行表示；

4）每一個映射點對應于一個小柵格內(nèi)的概率密度函數(shù)，計算這個點的正態(tài)分布；

5）對每個點的概率密度求和，做為本次匹配過程的分數(shù)值；

6）對本次的分數(shù)值使用Hessian矩陣法進行優(yōu)化，得到新的坐標參數(shù)；

7）回到步驟3，直到滿足收斂條件。

設(shè)坐標變換參數(shù)為p=（tx，ty，φ）T，則評價在該坐標變換參數(shù)下匹配效果的分數(shù)值為

式中：r′i為r根據(jù)變換參數(shù)p映射到基準圖像上的在基準圖像坐標系下的坐標；∑i是r′i對應的協(xié)方差矩陣；qi是 r′i的均值。NDT算法就是尋找使得score（p）最大的坐標變換參數(shù)，對于圖像匹配，該參數(shù)就是在NDT算法下的最佳匹配位置。

因為NDT算法是一種基于概率模型的匹配方法，在匹配過程中不需要對圖像的灰度、邊緣特征或者點特征進行提取和計算，NDT算法僅僅是計算一幅圖在另一幅圖上的映射點概率密度之和，能使概率密度之和最大的點即為匹配點，所以，該方法運算速度快，易于實現(xiàn)。

2　仿真研究

由于海洋環(huán)境中存在隨機噪聲，艦船本身也存在自噪聲，以及受不同海洋環(huán)境（如海底地形、水溫等）的影響，得到的聲納圖像往往不是最適合進行匹配的圖像，所以在進行匹配之前，對待匹配圖像進行圖像預處理是非常重要的，圖像預處理也是圖像匹配過程中不可缺少的環(huán)節(jié)［8］。在圖像匹配前對圖像進行降噪處理，將一些基本的圖像去噪方法應用于聲納圖像，包括均值濾波法、中值濾波法、小波去噪法和形態(tài)學平滑法。本文使用的均值濾波法和中值濾法的遍歷窗大小為3×3.

NDT算法是對圖像中目標像素點的坐標值進行計算，而與灰度值無關(guān)，所以，在對圖像進行匹配之前，先對圖像進行分割處理，將圖像中的目標和背景分開，為下一步的圖像匹配算法做準備。本文對4種分割算法進行了分析和仿真，分別研究了迭代閾值法、最大熵法、一維Otsu法和二維Otsu法［9］。

2.1去噪算法仿真

原始圖像選用蛙人圖像，對其加入高斯白噪聲，噪聲的均值為0，標準差 σ分別為0.1、0.2、0.3、0.4，分別利用均值濾波法、中值濾波法、小波去噪法以及形態(tài)學平滑法，對加入高斯白噪聲的圖像進行噪聲去除。然后計算出峰值信噪比（PSNR）和均方誤差（MSE），以這兩個值來進行比對，評價去噪算法的好壞。圖1是不同去噪方法對σ=0.4的噪聲圖像進行處理的效果對比，表1為σ=0.4時不同算法的運算時間以及峰值信噪比和均方誤差，表2為去噪算法去噪效果對比。

圖1　噪聲圖像去噪方法對比Fig.1 Comparison of several image denoising methods

表1　去噪算法時間對比Tab.1 Time of image denoising

表2　去噪算法去噪效果對比Tab.2 Parameters of image denoising algorithms

由表1和表2可知，從時間上看，利用均值濾波法和中值濾波法對蛙人圖像進行噪聲去除所用的時間要比小波去噪法和形態(tài)學平滑法所用的時間短，根據(jù)峰值信噪比，可以看出均值濾波法和中值濾波法較另外兩種算法要高一些，而均方誤差也是均值濾波法和中值濾波法較小。從對聲納圖像進行噪聲去除的實驗結(jié)果來看，形態(tài)學平滑法和小波去噪法對噪聲的去除效果最好，但同時也模糊了圖像，同均值濾波法相比，中值濾波法能在對圖像進行去噪的同時更好地保留一些邊緣信息。同小波去噪法和形態(tài)學平滑法相比，中值濾波法實現(xiàn)起來較為簡單方便，而且運行時間短，適合實際應用。經(jīng)過上述分析比較，選用中值濾波法作為本文的圖像預處理算法。

2.2分割算法仿真

對原始聲納圖像分別使用迭代閾值法、最大熵法、一維Otsu法和二維Otsu法進行分割處理，結(jié)果如圖2所示。

由表3可以看出，一維Otsu法和迭代閾值法運行所需時間要比一維熵法和二維Otsu法短。從分割閾值選擇上來看，一維熵法選擇的閾值最大；從圖像的分割結(jié)果可以看出，蛙人的一部分邊緣被歸為背景類，存在一定程度的過度分割。選取不同的分割方法會對下一步圖像匹配的效果造成影響，在考慮分割效果的同時，還要考慮算法實現(xiàn)的實時性，也就是運行速度，綜合考量后決定選用一維Otsu法。

2.3NDT算法仿真

NDT算法的仿真結(jié)果如圖3所示。

圖2　不同分割算法的結(jié)果對比Fig.2 Processed images of image segmentation methods

表3　不同分割算法的閾值和時間Tab.3 Time and threshold of image segmentation methods

圖3（b）和圖3（d）分別是基準圖像和待匹配圖像經(jīng)過一維Otsu法分割后的圖像，兩幅圖中用不同顏色對目標進行渲染，圖3（e）中白色的區(qū)域是兩幅圖像匹配重合區(qū)域，不同顏色的其他小區(qū)域分別為兩幅圖中沒有匹配的區(qū)域。從圖3中可以看出，兩幅圖像區(qū)別較大的區(qū)域是圖像上方的一部分，由于NDT算法的基本思想是基于概率密度之和最大化的匹配過程，在匹配過程中上面小區(qū)域占的概率比重小，而且兩幅圖上部區(qū)域包含的像素點位置不同，所以基本上沒有匹配，而蛙人區(qū)域所包括的像素點多，概率比重大，所以大部分都有匹配。

圖3　NDT算法的仿真Fig.3 Matlab results of NDT

3　硬件實現(xiàn)

用于NDT算法實現(xiàn)的數(shù)字信號處理器（DSP）開發(fā)平臺是以美國德州儀器公司產(chǎn)TMS320C6713［10］為核心的處理單元，軟件開發(fā)平臺選用CCS3.3.CCS不能讀取通用格式的圖像文件，利用Matlab軟件來實現(xiàn)對聲納圖像數(shù)據(jù)的讀取，將圖片保存為“.dat”格式的文件。在CCS中將.dat文件像數(shù)據(jù)下載到DSP的外圍存儲器（SDRAM）中，然后再對圖像數(shù)據(jù)做下一步的處理。

對NDT算法進行C語言編寫，實現(xiàn)了整個圖像匹配過程，主要包括對聲納圖像的去噪、圖像分割算法和圖像匹配算法。圖4為DSP實現(xiàn)的軟件流程圖，圖5為開發(fā)平臺的硬件實物圖，圖6～圖8為去噪、分割和匹配算法的硬件實現(xiàn)結(jié)果。

圖4　圖像匹配算法程序流程圖Fig.4 Flow chart of image matching processing

圖5　硬件實物圖Fig.5 Hardware

圖6　3×3中值濾波DSP硬件運行結(jié)果Fig.6 3×3 median filtering on DSP chip

4　結(jié)論

為滿足聲納圖像匹配的實時性需求，本文研究了NDT算法，該算法在匹配過程中無需對灰度特征或其他幾何特征進行運算，而是對基準圖像和待匹配圖像中目標點的概率分布進行匹配，運算速度快，更易于硬件實現(xiàn)。經(jīng)過比對分析，在匹配前選用中值濾波法和一維Otsu法對圖像進行去噪和分割處理，最后將圖像匹配處理全過程在DSP芯片上實現(xiàn)，并對蛙人圖像進行了實時的匹配處理，DSP實時處理結(jié)果與仿真結(jié)果一致。實驗結(jié)果表明NDT算法可實時地對聲納圖像進行匹配處理，具有良好的應用前景。

圖7　分割算法DSP硬件運行結(jié)果Fig.7 Image segmentation on DSP

圖8　NDT算法的DSP硬件運行結(jié)果Fig.8 NDT image matching algorithm on DSP chip

圖像配準計算量大，如何在保證匹配質(zhì)量的同時，提高匹配速度仍然是圖像匹配領(lǐng)域面臨的主要問題，可以從以下兩個方面進行深入研究：1）結(jié)合算法特性，選擇性能更優(yōu)越的DSP處理器，這樣可以大大加快算法運行過程中的運算速度；2）對NDT算法進行優(yōu)化。NDT匹配算法是基于圖像的概率模型匹配，在匹配之前需要給圖像分成小柵格，然后用概率密度函數(shù)對圖像進行表示。可以嘗試根據(jù)圖像長寬比來改變小柵格的形狀，來判斷是否可以得到更好的匹配效果。

參考文獻（References）

［1］ David G L.Distinctive image features from scale-invariant keypoints［J］.International Journal of Computer Vision，2004，60（2）：91-110.

［2］ Anma P E.Spatial registration of multispectral and mutitemporal digital image using fast Fouerier transformation technique［J］. IEEE Transactions on Geoscience Electronics，1970，8（4）：355-368.

［3］ 高晶，吳育峰，吳昆，等.基于角點檢測的圖像匹配算法［J］.儀器儀表學報，2013，34（8）：1717-1725. GAO Jing，WU Yu-feng，WU Kun，et al.Image matching method based on corner detection［J］.Chinese Journal of Scientific Instrument，2013，34（8）：1717-1725.（in Chinese）

［4］ 韓冰，王永明.基于一種快速歸一化積相關(guān)算法的圖像匹配研究［J］.兵工學報，2010，31（2）：160-165. HAN Bing，WANG Yong-ming.Research of image matching based on a fast normalized cross correlation algorithm［J］.Acta Armamentarii，2010，31（2）：160-165.（in Chinese）

［5］ 劉佳，傅衛(wèi)平，王雯，等.基于改進SIFT算法的圖像匹配［J］.儀器儀表學報，2013，34（5）：1107-1112. LIU Jia，F(xiàn)U Wei-ping，WANG Wen，et al.Image matching based on improved SIFT algorithm［J］.Chinese Journal of Scientific Instrument，2013，34（5）：1107-1112.（in Chinese）

［6］ Biber P，Strasser W.The normal distributions transform：a new approach to laser scan matching［C］∥2003 IEEE/RJS International Conference on Intelligent Robots and Systems.Las Vegas，US：IEEE，2003：2743-2748.

［7］ 蔡則蘇，洪炳镕，魏振華.使用NDT激光掃描匹配的移動機器人定位方法［J］.機器人，2005，27（5）：414-419. CAI Ze-su，HONG Bing-rong，WEI Zhen-hua.Localization of mobile robots by NDT laser scan matching algorithm［J］.Robot，2005，27（5）：414-419.（in Chinese）

［8］ 陳明生.圖像配準技術(shù)研究與應用［D］.長沙：國防科學技術(shù)大學，2006. CHEN Ming-sheng.Research on the image registration technology and its application［D］.Changsha：National University of Defense Technology，2006.（in Chinese）

［9］ Otsu A.Threshold selection method from gray-level histogram［J］. IEEE Transactions on Systems，Man，and Cybernetics，1979，9（1）：62-66.

［10］ Texas Instruments Incorporated.TMS320C6713B floating-point digital signal processor［M/OL］.［2005-08-12］.http：∥www. ti.com/lit/ds/symlink/tms320c6713.pdf.

中圖分類號：TB56

文獻標志碼：A

文章編號：1000-1093（2016）06-1052-06

DOI：10.3969/j.issn.1000-1093.2016.06.012

收稿日期：2015-11-16

基金項目：國家“863”計劃項目（2011AA09A107）

作者簡介：譚克（1990—），男，助理工程師，碩士。E-mail：tanke6353064@163.com；鄒麗娜（1976—），女，高級工程師，博士。E-mail：zoulina0717@163.com

The Application of Normal Distribution Transformation Algorithm in Sonar Image Processing

TAN Ke，ZOU Li-na
（Key Laboratory of Sonar Technology for National Defence，Hangzhou Applied Acoustics Research Institute，Hangzhou 310023，Zhejiang，China）

Abstract：For the requirements of the real-time processing of sonar image matching processing，the normal distributions transform（NDT）algorithm and its application in the sonar image processing are studied.The NDT algorithm is based on probability model.The coordinate of target point is calculated instead of the gray value of image，which speeds up the calculation speed of the NDT algorithm.Several methods of image denoising and segmentation are compared to give the results and analysis of each method.Based on TI DSP chip，the NDT image matching algorithm is achieved.

Key words：acoustics；image matching；digital signal processor；sonar image processing；normal distribution transform algorithm