成 芳，胡迺成

（中國人民解放軍91439部隊，遼寧 大連 116041）

基于三次樣條函數(shù)的側(cè)掃聲納定位數(shù)據(jù)平滑方法研究

成 芳，胡迺成

（中國人民解放軍91439部隊，遼寧 大連 116041）

針對側(cè)掃聲納定位采樣率與回波采樣率不同步的問題，提出基于三次樣條函數(shù)的定位數(shù)據(jù)平滑方法，將定位數(shù)據(jù)與回波數(shù)據(jù)一一對應(yīng)，避免聲納圖像產(chǎn)生空白區(qū)或像素重疊。實驗數(shù)據(jù)結(jié)果表明：該方法可對定位數(shù)據(jù)進行較好的平滑，有效消除拖魚因航向改變或航跡不穩(wěn)而產(chǎn)生的航跡線折點，從而消除了聲納圖像中的條形、扇形裂縫，解決了由采樣率不同引起的聲納圖像不連續(xù)現(xiàn)象，為后續(xù)地理編碼圖的處理提供了前提。

側(cè)掃聲納；圖像；Bezier函數(shù)；回波數(shù)據(jù)；GPS數(shù)據(jù)；采樣率；Ping

側(cè)掃聲納圖像是側(cè)掃聲納系統(tǒng)對海底的直觀反映[1-4]。如何根據(jù)各時刻回波數(shù)據(jù)形成聲納圖像是精細描繪海底地形和準確提取水下目標空間屬性信息的基礎(chǔ)。而目前，由于GPS采樣率與側(cè)掃聲納回波采樣率的不同[5-18]，造成在定位數(shù)據(jù)采樣時間間隔內(nèi)，不同Ping對應(yīng)同一位置，經(jīng)過地理編碼后側(cè)掃聲納圖像存在重疊，這將嚴重影響海底地形的精細描繪和海底目標物的準確提取，因此，側(cè)掃聲納成圖過程中定位數(shù)據(jù)的平滑便成為聲納圖像（地理編碼后）是否準確的制約因素。為此，本文重點研究GPS采樣率與回波采樣率不同步問題，提出基于三次樣條函數(shù)的定位數(shù)據(jù)平滑方法，實現(xiàn)定位數(shù)據(jù)與回波數(shù)據(jù)的完美融合，為后續(xù)聲圖的地理編碼提供了良好的前提條件。

1 基于三次樣條函數(shù)的定位數(shù)據(jù)平滑方法

1.1 方法提出

設(shè)Ti，Ti+1為定位系統(tǒng)中兩個相鄰點位Pi(yi, xi)，Pi+1(yi+1,xi+1)處的記錄時間，通過將Ti，Ti+1與Ping更新時間比較，可確定Ti～Ti+1時間內(nèi)，拖魚接收到的Ping線數(shù)量n，如圖1所示。

圖1 Ti～Ti+1時間內(nèi)，拖魚接收的Ping數(shù)量

設(shè)L為Ti～Ti+1時間內(nèi)拖魚航行距離。一般認為相鄰時間內(nèi)船速恒定，因此，在Ti～Ti+1時間段內(nèi)對航跡進行等距插值，使Ping線在Pi(yi,xi)，Pi+1(yi+1,xi+1)內(nèi)航跡均勻分布，則第k條Ping線對應(yīng)的航跡點Pk(yk,xk)距Pi(yi,xi)的距離為：

設(shè)過Pi(yi,xi)，Pi+1(yi+1,xi+1)的航跡函數(shù)為x= f(y)，則Ti～Ti+1時間內(nèi)拖魚航行距離為：

事實上，不同側(cè)掃聲納的采樣率及不同GPS的采樣率均不盡相同。但普遍地，側(cè)掃聲納采樣率要大于GPS采樣率，即fc=kf（G其中，fc為側(cè)掃聲納采樣率，fG為GPS采樣率，k＞1）。經(jīng)上述方法等距插值后，可消除由采樣率不同造成的航跡分布不均而引起的條帶裂縫，但由于航向改變或航跡不穩(wěn)造成的Ping線交叉，使聲納圖像依然存在扇形裂縫，見圖2。

圖2 等距插值后航跡轉(zhuǎn)彎處可能出現(xiàn)的裂縫

事實上，高次插值函數(shù)的計算量大，有劇烈振蕩，且數(shù)值穩(wěn)定性差；分段插值中，分段線性插值在分段點上僅連續(xù)而不可導，分段三次埃爾米特插值有連續(xù)的一階導數(shù)，如此光滑程度常不能滿足需要，樣條函數(shù)可以同時解決這兩個問題，使插值函數(shù)既是低階分段函數(shù)，又是光滑的函數(shù)，并且只需在區(qū)間端點提供某些導數(shù)信息。因此利用三次樣條函數(shù)（三次Bezier曲線）對航跡點進行平滑處理，來消除扇形裂縫。即由Bezier曲線方程x=fBezier(y)替換x=f(y)，由LBezier替換L，由FBezier(y)替換F(y)，于是式（1）和式（2）可化為式（4）和式（5）：

即可消除聲納圖像中由航跡不穩(wěn)造成的扇形裂縫，其效果見圖3。

圖3 航跡三次樣條函數(shù)平滑后的效果

1.2 算例分析

選取某海域原始側(cè)掃聲納數(shù)據(jù)作為實驗數(shù)據(jù)，該數(shù)據(jù)由EdgeTech4200-FS型側(cè)掃聲納[7]測量獲得。在進行該測線測量時，發(fā)射脈沖中心頻率為403 kHz，發(fā)射脈沖寬度為3 ms，量程為75 m，數(shù)據(jù)采樣間隔為0.015 6 ms。

按該型號側(cè)掃聲納的數(shù)據(jù)格式對原始數(shù)據(jù)進行解析[19]，分離出每條掃描線的序號和拖魚航跡坐標經(jīng)緯度。從分離出的數(shù)據(jù)可看出，定位數(shù)據(jù)約5Ping更新一次，實例說明了定位采樣率與回波采樣不同步的現(xiàn)象。若不進行定位數(shù)據(jù)的平滑處理，在像素重采樣時，易造成相鄰多條掃描線相互重合，而在航跡特征點之間形成條帶裂縫，為后續(xù)地理編碼帶來不便。

采用上節(jié)提出的方法對定位數(shù)據(jù)進行處理，處理結(jié)果如圖4～圖6。

圖5 定位數(shù)據(jù)經(jīng)三次樣條函數(shù)平滑后的聲納圖像

圖4 為定位數(shù)據(jù)經(jīng)等距插值后的聲納圖像，圖5為對等距插值后的定位數(shù)據(jù)，利用三次樣條函數(shù)（Bezier曲線）平滑形成聲納圖像。圖4中雖然消除了條形裂縫現(xiàn)象，但仍存在扇形裂縫，其主要分布在航跡拐點處，且航跡拐角越大，扇形裂縫也就越大。圖5中海底成像效果較好，聲納圖像平滑，消除了扇形裂縫。

通過以上實例，可以得到以下結(jié)論：

（1）由于定位采樣率相對回波采樣率低，造成聲納圖像中灰度重合或產(chǎn)生裂縫現(xiàn)象，其中灰度重合表現(xiàn)為相鄰多Ping回波歸算點位重合，裂縫現(xiàn)象表現(xiàn)為相鄰Ping間隔大；

（2）定位數(shù)據(jù)的等距內(nèi)插可有效消除聲納圖像中條形裂縫現(xiàn)象，但扇形裂縫仍存在；

（3）對等距插值后的定位數(shù)據(jù)進行三次樣條函數(shù)擬合，定位數(shù)據(jù)得到較好的平滑，消除了拖魚因航向改變或航跡不穩(wěn)而產(chǎn)生的航跡線折點，從而消除了扇形裂縫。

2 結(jié)論

本文對側(cè)掃聲納定位采樣率與回波采樣率不同步的問題進行了研究，旨在解決如何根據(jù)側(cè)掃聲納系統(tǒng)對海底或測量環(huán)境的離散采樣數(shù)據(jù)，獲取平滑的海底聲學圖像。通過本文研究，得到以下結(jié)論：

（1）由于定位采樣率與回波采樣率不同，在定位數(shù)據(jù)采樣時間間隔內(nèi)，不同Ping對應(yīng)同一個位置，而航跡特征點之間又無ping線，經(jīng)過地理編碼后側(cè)掃聲納圖像存在重疊或裂縫；

（2）利用等距插值法對定位數(shù)據(jù)進行平滑，可有效消除聲納圖像中的條形裂縫，但當拖魚航向發(fā)生改變或航跡不穩(wěn)時，無法消除折點處的扇形裂縫；

（3）采用本文提出的基于三次樣條函數(shù)的定位數(shù)據(jù)平滑方法對等距插值后的定位數(shù)據(jù)進行平滑，消除了拖魚航跡線上的折點，從而消除了扇形裂縫。

[1]來向華,潘國富,茍諍慷,等.側(cè)掃聲納系統(tǒng)在海底管道檢測中應(yīng)用研究[J].海洋工程,2011,29(3):117-121.

[2]劉雁春,暴景陽,李明叁.我國海洋測繪技術(shù)的新進展[J].測繪通報,2007(3):1-7.

[3]劉雁春,翟京生,翟國君,等.海洋測繪學科發(fā)展現(xiàn)狀[R].中國測繪學科發(fā)展藍皮書(2004卷),北京:測繪出版社,2005.

[4]FlemmingBW.Side scan sonar:a practical guide[J].The International Hydrographic Review,1976,53(1):65-92.

[5]Discover 4200-MP Side Scan Processor Software User’s Manual[R].edgetech,2008.

[6]4200 Series High Definition Dual FrequencySide Scan Sonar System[R].edgetech,2008.

[7]4200-FSHigh Definition,Multi-pulse,Dual FrequencySide Scan SystemUser’s Hardware Manual[R].edgetech,2008.

[8]4200 Side Scan Sonar SystemUser Hardware Manual[R].edgetech,2010.

[9]J-Star Full SpectumDiagnostic Software Technical&User’s Manual[R].edgetech,2010.

[10]2000 Series Combined Side Scan Sonar and Sub-BottomProfilingSystemUser’s Manual[R].edgetech,2010.

[11]3100-P Portable Sub-BottomUser Hardware Manual[R].edgetech,2014.

[12]3100-P Portable Sub-BottomProcessor Sotfware[R].edgetech,2007.

[13]3200-XSSub-BottomSystemUser Hardware Manual[R].edgetech,2015.

[14]4125 Side Scan Sonar SystemUser Hardware Manual[R].edgetech,2010.

[15]4125 Magnetometer Interface User Hardware Manual[R].edgetech,2015.

[16]Discover 4125 Side Scan Processor Software User’s Manual[R].edgetech,2010.

[17]6205 Bathymetry&Side Scan Sonar User Hardware Manual[R].edgetech,2015.

[18]劉江,蔡伯根,唐濤，等.基于GPS的列控軌道地圖數(shù)據(jù)生成方法研究[J].測繪學報,2011,40(1):111-117.

[19]Description ofthe Edgetech(.jsf)File Format[R].edgetech,2011.

Research on the Smoothing Method of Side Scan Sonar Positioning Data Based on Triple Spline Function

CHENG Fang,HU Nai-cheng
Troop 91439 of the People's Liberation Army,Dalian 116041,Liaoning Province,China

Aimed at solving the problem of inconformity of positioning and echo sampling rates,this paper puts forward the smoothing method of side scan sonar positioning data based on the Triple Spline Function.Through the one-to-one correspondence between the positioning data and echo data,blank area and pixel overlap of sonar images are eliminated.The results indicate the feasibility of this method,which can effectively eliminate the line inflection and thus bars or sector gaps,so as to solve the problem of image discontinuities caused by different sampling rates.It offers a premise for processing subsequent geocoding.

side scan sonar;image;Bezier function;echo data;GPS data;sampling rate;Ping

P229.1

A

1003-2029（2017）03-0023-04

10.3969/j.issn.1003-2029.2017.03.005

2017-01-08

國家自然科學基金資助項目（41374108，41476087）

成芳（1981-），女，博士，工程師，主要從事系統(tǒng)總體與海洋環(huán)境測量研究。E-mail：chengfangtl@163.com