庫安邦，楊 鯤

(1.交通運輸部天津水運工程科學研究所，天津 300456；2.天津水運工程勘察設(shè)計院 天津市水運工程測繪技術(shù)重點實驗室，天津 300456)

海底環(huán)境錯綜復雜，受海洋水動力、海流的沖刷等影響，海底管道下方沉積物容易被掏蝕從而出現(xiàn)裸露、懸空等現(xiàn)象，懸跨長度過長就會使管道發(fā)生斷裂[1]，出現(xiàn)油氣泄漏，經(jīng)濟蒙受巨大損失的同時對海洋生態(tài)環(huán)境造成惡劣的影響。因此人們越來越重視海底管道的檢測和維護，側(cè)掃聲吶系統(tǒng)以其低成本、高效率、高分辨率等優(yōu)勢在海底管道檢測方面應(yīng)用越來越廣泛[2-6]。來向華[7-8]等針對單波束和側(cè)掃聲吶探測海底管道進行了探究，驗證了單波束和側(cè)掃聲吶系統(tǒng)檢測海底管道的實用性；安永寧[9]等提出了針對管道懸空具體的治理方法；荊少東[10]考慮了不同海底地貌對海底管道穩(wěn)定性的影響；王雷[11]等分析了海底管道存在的不同狀態(tài)，根據(jù)側(cè)掃聲吶系統(tǒng)探測結(jié)果，針對不同的管道狀態(tài)提出了對應(yīng)的計算方法；張彥昌[12]、張永明[13]、馮百全[14]等人針對綜合側(cè)掃聲吶、淺地層剖面儀、多波束系統(tǒng)等多種手段探測海底管線進行了介紹，對比分析了各個聲學探測方法的優(yōu)劣性，最后得出了綜合探測效果最佳的結(jié)論。

雖然前人對側(cè)掃聲吶探測海底管道進行了大量的總結(jié)和論述，但是傳統(tǒng)的分析和研究多基于商業(yè)軟件，各廠商軟件雖然可以實現(xiàn)簡單的瀑布圖顯示和處理，但鑒于各個公司軟件側(cè)重功能不同，無法滿足針對性的需求。對感興趣的區(qū)域研究或直接在圖像顯示軟件上截取灰度圖，針對截取的灰度圖做進一步的分析，顯示軟件上的目標灰度圖像經(jīng)過了系統(tǒng)的壓縮、轉(zhuǎn)換處理等，在一定程度上存在失真和畸變，為保證后續(xù)數(shù)據(jù)二次開發(fā)的準確性，本文以XTF格式側(cè)掃聲吶數(shù)據(jù)為例，解析原始數(shù)據(jù)并對其進行一系列的改正處理，最終得到完整的、較準確的管道聲吶圖像，并將其集成開發(fā)了聲吶管道數(shù)據(jù)可視化系統(tǒng)。

1 XTF格式數(shù)據(jù)結(jié)構(gòu)

XTF文件由一個頭文件開始，包含文件的標識和數(shù)據(jù)通道信息。一個頭文件至少是1 KB，文件信息占256個字節(jié)，包含聲吶類型、通道數(shù)量等信息；一個頭文件結(jié)構(gòu)體包含六個通道空間，主要記錄通道的類型、采樣精度等信息。頭文件之后就是各個數(shù)據(jù)包，每個數(shù)據(jù)包包含各自的頭文件，用于識別數(shù)據(jù)包的類型和大小等，各數(shù)據(jù)包相互獨立，可以按需讀取。XTF文件中主要的數(shù)據(jù)包類型有測深、聲吶等，數(shù)據(jù)包頭文件后邊就是通道數(shù)據(jù)，通道數(shù)據(jù)包含通道頭文件和采樣數(shù)據(jù)，通道頭文件中包含聲吶系統(tǒng)采樣方式、采樣率等信息，對于側(cè)掃聲吶采樣數(shù)據(jù)就是回波強度。其他幾個數(shù)據(jù)包的結(jié)構(gòu)分布和聲吶數(shù)據(jù)包的結(jié)構(gòu)是一致的，本次研究僅針對聲納數(shù)據(jù)開展研究，其他就不再贅述。

2 關(guān)鍵技術(shù)

從采集的原始側(cè)掃聲吶數(shù)據(jù)入手，分析XTF格式側(cè)掃聲吶數(shù)據(jù)成圖關(guān)鍵步驟，著重介紹了數(shù)據(jù)解碼、海底線檢測和航跡處理等部分，最后集成各關(guān)鍵步驟算法，集成開發(fā)了海底管道可視化系統(tǒng)，高效的管理和使用側(cè)掃聲吶海底管道數(shù)據(jù)。

2.1 數(shù)據(jù)解碼

根據(jù)XTF文件的組成結(jié)構(gòu)，按照圖1順序讀取XTF文件，XTF頭文件主要包括字段、字節(jié)偏移量和注釋三部分，字段表示各部分名稱，字節(jié)偏移量對應(yīng)字段開始的位置，根據(jù)這一特點，讀取數(shù)據(jù)時可合理利用指針工具，根據(jù)字節(jié)偏移量針對性的讀取聲吶數(shù)據(jù)包內(nèi)的所有內(nèi)容。利用指針偏移讀取數(shù)據(jù)包頭文件后就進入到通道記錄的數(shù)據(jù)，經(jīng)過通道數(shù)據(jù)頭文件的“放行”后，就見到了重要的聲強采樣數(shù)據(jù)，根據(jù)高低頻，左右舷的先后順序，依次讀取每個通道中的采樣數(shù)據(jù)，再依次按順序組合。

圖1 XTF數(shù)據(jù)流結(jié)構(gòu)圖Fig.1 XTF structure diagram of data flow

2.2 閾值法檢測海底線

2-a 閾值法檢測海底線 2-b 海底線跟蹤圖2 閾值法檢測海底線Fig.2 Threshold method for detection of submarine line

海底線由海底點連接形成，每個海底點是單Ping掃描序列的第一個海底強回波。傳播至海底的回波信號會發(fā)生階躍性變化，可通過設(shè)定合適的閾值I，按照發(fā)射接收的先后順序，提取每個回波序列中首個大于I的回波并認為該回波為到達海底的第一個回波(圖2)，也就是拖魚至海底的高度，連接每一個海底點就得到了海底線，本文設(shè)定閾值檢測首個回波采用左-右綜合提取法，當一側(cè)“水柱區(qū)域”存在懸浮物時，多個異常點會對首個回波檢測產(chǎn)生干擾，此時利用海底線左右對稱原則，檢測另一側(cè)海底點對稱該對應(yīng)區(qū)域，此法對于復雜海況聲吶數(shù)據(jù)海底線跟蹤適用性更強。

表1 聲吶圖像部分Ping對應(yīng)的時間和位置信息Tab.1 The time and location with part of sonar image

2.3 航跡處理

表1為解碼出的部分Ping對應(yīng)的時間和位置信息，可以看出側(cè)掃聲吶采樣間隔約為0.12 s；GPS記錄通常以秒為單位，前5Ping的位置信息是一樣的，7、8Ping位置信息也是一樣的，因此要想得到每一Ping的真實位置，保證地理編碼的準確性，需要對其位置進行重采樣。

如圖3可知，原始航跡點非常稀疏，經(jīng)等距插值后每一Ping均具有位置信息，但當航向不穩(wěn)定或有轉(zhuǎn)角的地方，插值后仍存在扇形裂縫無法消除，將會影響圖像的細節(jié)特征，因此對插值后的航跡采取三次B樣條函數(shù)平滑，本次實驗數(shù)據(jù)平滑后的航跡有效的刪除了扇形裂縫和重疊，更加符合真實效果。當測量船急轉(zhuǎn)彎或避讓測線上障礙物而導致的航跡劇烈變化時，此時需在對原始航跡點插值前進行預(yù)處理編輯，如剔除異常點之后再進行插值擬合，或直接將航跡以突變區(qū)為節(jié)點截斷分開處理，最后在對應(yīng)時間位置進行地理編碼。

2.4 地理編碼

像素地理編碼是指根據(jù)處理后的航跡數(shù)據(jù)，將每一個采樣點的位置納入指定的坐標系，賦予聲吶圖像每一個像素點位置信息，主要步驟包括：(1)建立地理編碼模型；(2)定義像素分辨率；(3)計算地理編碼矩形；(4)像素點坐標變換?？紤]篇幅大小，在此不再贅述具體過程。

3-a 等距離插值航跡圖3-b 三次B樣條函數(shù)平滑航跡圖圖3 等距離插值和三次B樣條函數(shù)擬合航跡圖Fig.3 Equidistance interpolated and three B spline fitted track map

圖4 側(cè)掃聲吶海底管道數(shù)據(jù)可視化系統(tǒng)Fig.4 Visualization system of side scan sonar data of submarine pipeline

3 海底管道可視化系統(tǒng)

前面介紹了從數(shù)據(jù)解碼一直到地理編碼完成的關(guān)鍵步驟。為更方便快捷的處理側(cè)掃聲吶海底管道數(shù)據(jù)，集成關(guān)鍵算法，編寫了一套海底管道可視化系統(tǒng)，界面簡潔明了，實現(xiàn)了側(cè)掃聲吶海底管道數(shù)據(jù)預(yù)處理的友好的人機交互界面。系統(tǒng)的主要框架如圖4所示，主要包括數(shù)據(jù)顯示、處理和評價三個方面，具體包括原始側(cè)掃聲吶管道數(shù)據(jù)的可視化顯示，可以直觀的展示懸空裸露段的海底管道，然后根據(jù)需要進行針對性的摘取分析、處理，提高了數(shù)據(jù)分析的效率，另外還可以對原始文件里提取的管道數(shù)據(jù)和輔助信息進行二次開發(fā)利用，并對其數(shù)據(jù)質(zhì)量進行客觀的評價。

檢測到海底線后，可以對聲吶圖像進行畸變校正等操作，主要包括斜距改正和輻射校正，圖5為改正前后的對比圖，斜距改正后的聲吶圖像消除了距離向的畸變，經(jīng)過輻射校正消除了縱向灰度不均勻產(chǎn)生的圖像畸變。斜距改正的好壞在于海底線跟蹤的準確性，本文采用左-右綜合閾值檢測法跟蹤海底線，可靠性和準確性更強；輻射校正先根據(jù)式(1)計算出聲吶圖像每列像素灰度值相對整個聲吶圖像灰度的改正因子，然后根據(jù)(2)式改正每列聲吶圖像像素灰度的不均衡，從而達到消除輻射畸變的目的。

(1)

式中：i是行號；j是列號；n為每列像素個數(shù)；m為聲吶圖像總像素數(shù)；xij表示i行j列像素的灰度值。

(2)

圖5 聲吶圖像畸變改正前后對比圖Fig.5 Comparison of sonar images before and after distortion correction

另外該軟件還可以提取姿態(tài)、時間等其他輔助信息，并將它們以文本的形式保存，如圖6為提取的部分輔助信息并將其以文本形式存儲以備二次利用，根據(jù)不同的轉(zhuǎn)換參數(shù)，在航跡處理環(huán)節(jié)，可以提取原始數(shù)據(jù)中記錄的坐標信息，可實現(xiàn)大地坐標和直角坐標等之間的相互轉(zhuǎn)換，改正后的航跡應(yīng)用到地理編碼中，最終得到具有地理位置信息的無畸變聲吶圖像。

圖6 輔助信息存儲文本Fig.6 Auxiliary information stores by text

4 結(jié)論

本文從實測側(cè)掃聲吶海底管道數(shù)據(jù)入手，介紹了數(shù)據(jù)處理過程的關(guān)鍵部分，最后對各個環(huán)節(jié)進行了整合和封裝，開發(fā)了一套海底管道可視化系統(tǒng)，可以快速的對實測聲吶管道數(shù)據(jù)進行處理和相關(guān)信息提取，為數(shù)據(jù)的二次開發(fā)利用提供了很好的條件，本系統(tǒng)的創(chuàng)新性主要在于以聲吶管道數(shù)據(jù)處理為切入點，綜合各個關(guān)鍵步驟算法，可以對海底顯著特征物進行針對性的分析，軟件相對簡單實用，另外可以根據(jù)需要針對性的添加對應(yīng)的處理算法；但仍存在需要完善的地方，由于開發(fā)語言存在一定的低效性，所以對于單個數(shù)據(jù)量很大的聲吶條帶存在一定的局限性，后續(xù)筆者會繼續(xù)深入研究并優(yōu)化各個算法，增加條帶鑲嵌等功能，使整個系統(tǒng)更加完善實用。