楊 波 汪 偉 劉燁瑤 李欣國 梁 濤郭衛(wèi)振 廖佳偉 潘 鋒

①(中國科學(xué)院聲學(xué)研究所 北京 100190)

②(北京海洋聲學(xué)裝備工程技術(shù)研究中心 北京 100190)

③(中國科學(xué)院大學(xué) 北京 100049)

④(中國科學(xué)院信息工程研究所 北京 100093)

1 引言

單波束前視掃描聲吶廣泛應(yīng)用于水下機器人、載人潛水器等運動平臺。由于其機械掃描的工作特性產(chǎn)生了以下兩個問題：為了獲取較快的成像速度，機掃聲吶角度步進較大，導(dǎo)致分辨率較低；水下載體快速航行引起聲學(xué)中心的位置變化，影響成像效果。這些影響將會增加后續(xù)避障[1]、目標(biāo)檢測[2]、特征提取[3]、目標(biāo)識別[4]、跟蹤[5–7]等功能的難度。

針對第一個問題，許多研究者引入了插值的方法獲得高分辨率圖像。最鄰近插值法[8]用距離待插值點最近的整數(shù)坐標(biāo)點像素值替代該點，算法簡單、速度快，但圖像會出現(xiàn)方塊效應(yīng)和鋸齒效應(yīng)。雙線性插值法[9]將待插值點的灰度用與該點周圍4個點的灰度值線性擬合而成，比最鄰近插值平滑、運算速度較快，但平滑操作損失了圖像中許多重要的高頻信息。Keys[10]提出了用待插值點周圍16個相鄰點的灰度值加權(quán)平均來計算插值，相比雙線性插值提高了插值精度，但不能很好地保持圖像的細(xì)節(jié)部分。Trucco等人[11]將線性插值引入前掃聲吶成像中，在極坐標(biāo)中根據(jù)不同距離進行分段波束插值。Zhang等人[12]在Trucco的基礎(chǔ)上進行改進，在波束域做拉格朗日2次插值，提高了圖像質(zhì)量。前面介紹的聲吶圖像都是基于前向映射，將極坐標(biāo)系點映射到直角坐標(biāo)系中得到的是非整數(shù)點，經(jīng)四舍五入后導(dǎo)致極坐標(biāo)系中點與直角坐標(biāo)系中點形成多對一映射?？苜t聰[7]在前視聲吶成像中引入了后向映射，解決了上述問題，圖像更加平滑。Wang等人[13]在采用后向映射的同時，分別在波束域以及每個波束內(nèi)進行了8點2次插值，解決了采樣點稀疏的問題。本文在研究前人方法的基礎(chǔ)上，采用后向映射方法映射出需要繪制的像素點在極坐標(biāo)系的坐標(biāo)，在極坐標(biāo)系作雙線性插值得到準(zhǔn)確位置像素點的值，然后重構(gòu)圖像。

對于運動平臺的前視掃描聲吶成像，載體在水下航行時圖像將會產(chǎn)生徑向誤差。載體運動對前視掃描聲吶成像的影響，前人分析的多是載體的姿態(tài)變化對前視掃描聲吶成像的影響。梁艷[14]介紹了利用DVL(聲多普勒計程儀)測得的速度和姿態(tài)對前視掃描聲吶圖像中目標(biāo)物的位置進行校正，補償了由于載體運動造成的圖像失真。但梁艷直接將前視掃描聲吶測得目標(biāo)數(shù)據(jù)當(dāng)作目標(biāo)在水平面中與聲吶的距離，然后在水平面中對目標(biāo)位置進行校正，忽略了聲吶測得的目標(biāo)距離是在空間中斜距。前視掃描聲吶成像顯示并不是目標(biāo)在水平面中的距離，而是根據(jù)目標(biāo)與聲吶陣的距離映射到水平面中。本文通過建立載體運動模型，分析載體運動時目標(biāo)成像的徑向誤差，推導(dǎo)出圖像域運動補償量，并用泰勒公式對補償量進行化簡。隨后用簡化的補償量處理前視聲吶回波數(shù)據(jù)，并與未經(jīng)運動補償?shù)膱D像進行比較，經(jīng)過運動補償?shù)膱D像能準(zhǔn)確表示目標(biāo)整體信息，驗證了圖像域近似運動補償?shù)挠行浴?/p>

2 單波束前視掃描聲吶成像系統(tǒng)

單波束前視掃描聲吶每次定時發(fā)射一個具有一定脈寬的信號，在某個固定方向形成一個探測波束。接收回波就能獲得該方向上物體信息。對回波信號進行采樣并與本地信號進行相關(guān)，相關(guān)結(jié)果代表了該方向不同距離回波強度。在發(fā)射回波后經(jīng)過一段固定時間，按固定角度步進自動轉(zhuǎn)動聲吶陣，進入第2個方向掃描。經(jīng)過多次轉(zhuǎn)動就可獲得一個扇區(qū)的回波信息。

3 后向映射及插值成像

基于單波束前視掃描聲吶工作特點，為提高成像速度需增大角度步進，從而導(dǎo)致波束間隔變大，生成的圖像間隙較大，容易產(chǎn)生馬賽克現(xiàn)象。成像前需要對波束數(shù)據(jù)作插值處理，再將極坐標(biāo)數(shù)據(jù)轉(zhuǎn)換成笛卡兒坐標(biāo)數(shù)據(jù)。轉(zhuǎn)換后的坐標(biāo)需要進行四舍五入取整，導(dǎo)致近距離的像素點會被多次計算，產(chǎn)生像素模糊，遠(yuǎn)距離的像素點反而會缺失灰度值。坐標(biāo)點經(jīng)四舍五入后位置與真實位置會產(chǎn)生偏差。基于上述因素考慮本文采用后向映射成像方法[9]，即在笛卡兒坐標(biāo)系下對像素進行處理。將需要重繪的像素坐標(biāo)映射至極坐標(biāo)系，在極坐標(biāo)系下查找以及插值。

具體實現(xiàn)時根據(jù)單波束前視掃描聲吶機械掃描面的工作特點，每接收到1ping數(shù)據(jù)，只更新當(dāng)前ping區(qū)域的圖像灰度值。這樣不僅加快了成像速度，也減少了計算量。首先要確定需要在笛卡爾坐標(biāo)系下重繪的區(qū)域，如圖1中兩根紅線之間的小扇形區(qū)域，可以根據(jù)前1ping和當(dāng)前ping的角度來確定。

先計算出小扇形區(qū)域的兩條邊界線，即

圖1 坐標(biāo)轉(zhuǎn)換

其中， H是前視掃描聲吶輸出的轉(zhuǎn)動方向，i 表示ping計數(shù)。然后在笛卡兒坐標(biāo)系下矩形區(qū)域，逐行搜索需要映射的像素點。其中矩形區(qū)的像素為512×256。小扇形區(qū)的每一行的起始像素點可由l1計算出，終止像素點由l2計算出。第m 行的起始像素點為 (256 ?m)/(tan(H(i)/180×π))，終止像素點 為 (256 ?m)/(tan(H(i ?1)/180×π))。實 際搜索的區(qū)域由小扇形區(qū)近似成三角區(qū)。然后在極坐標(biāo)系中對徑向距離超過250的點直接賦值0，從而刪除超出小扇形區(qū)的點。后向映射的坐標(biāo)轉(zhuǎn)換公式為

其中， m, n 分別是像素點的行和列，α , r 分別是極坐標(biāo)系下方向角和徑向距離。由于接收的波束數(shù)據(jù)是以方向角為行、徑向距離為列的矩陣進行存儲，后向映射點的坐標(biāo) α需要轉(zhuǎn)換成矩陣的行。具體轉(zhuǎn)換公式為

其中，前視掃描聲吶的波束數(shù)據(jù)的方向角范圍為±67.2?， 波束角度步進為2.4°, θ, r分別是存儲矩陣的行和列。

本文采用的雙線性插值在波束方向和徑向都做線性插值。具體作法如圖2所示。

圖2 雙線性插值

根據(jù)后向映射求得 M點坐標(biāo)，經(jīng)過向下、向上取整找到相鄰的 P1, P2, Q1和Q24 點。P1和P2作線性插值得到 P12, Q1和Q2作線性插值得到Q12，之后P12和Q12作線性插值得到M 。

4 圖像域運動補償

載體運動會使得聲學(xué)中心與目標(biāo)點的相對位置發(fā)生變化，同一目標(biāo)點在不同ping的掃描角度也會產(chǎn)生差異。由于前視掃描聲吶的工作特性是每更新1ping數(shù)據(jù)，刷新一次圖像。這樣更新一個完整的扇形圖像，從第1ping到最后1ping需要一個較長的時間。更新時間計算公式為

其中， θ是機掃的半個扇區(qū)的角度，? θ是機掃角度步進量， ?t表示機掃步進一次時長?！吧詈Ｓ率刻枴陛d人潛水器在水下最大航速1.25 m/s，扇區(qū)范圍± 67.2?，角度步進2.4°，角度步進時長0.14 s，機掃整個扇區(qū)，潛水器會前進9.8 m。聲吶實際工作距離50 m，對于一個扇形區(qū)內(nèi)，第1ping和最后1ping成圖在運動方向誤差將近似達(dá)到作用距離的19.6%。這將直接影響后續(xù)圖像處理以及操作人員判斷。針對上述問題本文提出圖像域運動補償，即在每次成圖前，考慮航行器的運動對成像影響。由于前視掃描聲吶成像顯示的并不是直接在水平面的投影，Hurtos等人[15]分析了前視聲吶的成像投影模型。由于前視聲吶成像與前視雷達(dá)成像有許多相似之處，本文在Hurtos的成像投影模型的基礎(chǔ)上同時參考了張寅的運動平臺前視雷達(dá)成像[16]，在此基礎(chǔ)上進行成像的徑向誤差分析。載體運動幾何示意圖如圖3所示。

圖3 前視聲吶掃描成像幾何示意圖

載體在距底高度 Ht的 位置以速度vt沿箭頭方向運動，前視聲吶發(fā)射陣按速度 ω進行掃描。獲得1ping載體與目標(biāo)的距離后，將距離信息呈現(xiàn)在扇形圖中。由于載體運動，在繪制當(dāng)前ping數(shù)據(jù)時，上1ping的目標(biāo)在當(dāng)前圖像中位置已發(fā)生變化。為了實時準(zhǔn)確顯示前視圖像，需要知道目標(biāo)在前視成像中不同ping的距離關(guān)系，并在當(dāng)前ping成像時進行糾正。我們以前視成像區(qū)域中單個點目標(biāo)為例，分析不同時刻載體與目標(biāo)的距離如圖4所示。

載體初始時刻 t0，與目標(biāo)S (x,y)的 距離R0，發(fā)射陣的角度為θ 。在t 時刻，載體與目標(biāo)的斜距可表示為

圖4 不同時刻載體與目標(biāo)的距離

徑向誤差表現(xiàn)在圖像域，則體現(xiàn)在目標(biāo)在載體運動方向的移動，如圖5所示。

圖5 目標(biāo)在圖像域徑向誤差

當(dāng) t 與t0的間隔為一個掃描步進時間間隔時，式(10)左邊可以理解為? d，即為一次掃描步進需要補償?shù)牧?。?dāng)θ=0, φ =30, R0=20時，由式(10)可計算出 d =0.1514 m。前視聲吶的頻率f 為675 kHz，對應(yīng)的波長λ 為0.0022 m。由此可見誤差d 遠(yuǎn)大于波長，誤差的補償具有明顯意義?？梢詫κ?10)右邊第2項進行泰勒展開，忽略2次項及高階項得到：

當(dāng)θ=0, φ =30, R0=20時，由式(11)可計算出?d為0.1516，偏差僅為0.13%，能夠滿足補償要求。式(11)可以理解為成像時顯示的是一個與水平面夾角為 φ的平面，其中φ =arcsin(H0/R0)， H0為t0時刻對底高度。水平面的移動需要對應(yīng)拉伸到該平面內(nèi)。每次成像在計算小扇區(qū)像素前，先將上次圖像往聲吶運動方向搬移 ?d。加上運動補償后再計算小扇區(qū)的像素，然后更新整幅圖像。這樣就能實時修正由載體運動帶來的運動方向上誤差。此處的? d 理論上需要根據(jù)目標(biāo)的每個像素點對應(yīng)的φ 計算得到，實際工作時可以根據(jù)聲吶作用范圍對? d 做近似處理，目標(biāo)點取相同的? d 。

5 仿真結(jié)果

對一組實際采集的前視聲吶數(shù)據(jù)分別進行前向映射不做插值、前向映射做2次插值、后向映射2次插值、后向映射雙線性插值。具體的結(jié)果如圖6所示。

圖6 4種插值成像

從圖6(a)可以看出前向映射不做插值的圖像中，各波束數(shù)據(jù)在遠(yuǎn)端有明顯的間隙，經(jīng)過2次兩點插八點后圖像間隙消除了，如圖6(b)所示，但在遠(yuǎn)端形成了明顯的蝶形圖案。改為后向映射2次插值兩點插四點后，可消除蝶形圖案，但還存鋸齒效應(yīng)，如圖6(c)所示。改為后向映射雙線性插值后，可改善鋸齒效應(yīng)的影響，如圖6(d)所示。

為了驗證本文對前視聲吶運動成像的補償效果，選擇一個長40 m的長條目標(biāo)作為前視聲吶成像的仿真目標(biāo)。前視聲吶設(shè)備距底高度10 m，目標(biāo)位于前視聲吶前方水平距離為20 m。分別對前視聲吶靜止、以1.25 m/s前進兩種狀態(tài)進行仿真。具體的仿真結(jié)果如下。

圖7 聲吶靜止、前進、直接補償、本文補償對比

從圖7(b)可以看出，前視聲吶在前進時左側(cè)的目標(biāo)圖像由于未及時更新而產(chǎn)生6.2 m的距離誤差，約為工作距離的12.4%。圖7(c)是直接利用聲吶前進速度補償未及時更新的目標(biāo)的位置，左側(cè)的目標(biāo)與前視聲吶的縱向距離偏差在減小，但補償?shù)牧科?，?dǎo)致左側(cè)目標(biāo)與前視聲吶的縱向距離更近，左側(cè)目標(biāo)圖像距離誤差達(dá)到約1.4 m，約為工作距離的2.8%。圖7(d)采用的是本文的算法對目標(biāo)進行補償，可以看出左右兩側(cè)的目標(biāo)與前視聲吶的縱向距離基本一致，左右兩端的偏差小于0.2 m，約為工作距離的0.4%。從上述仿真結(jié)果來看本文的算法能夠補償前視聲吶前進造成的圖像偏差，并且比直接補償法效果更好。

6 實驗結(jié)果

實驗采用是在南海1000 m深度區(qū)域的前視數(shù)據(jù)。前視聲吶的具體工作參數(shù)如表1所示。

前視聲吶具體的圖像后置處理流程如圖8所示。

先用更新的單ping數(shù)據(jù)作后向映射、雙線性插值得到小扇形區(qū)的圖像。其余區(qū)域的圖像通過之前生成的圖像先根據(jù)艏向角變化進行旋轉(zhuǎn)調(diào)整，然后進行圖像域運動補償來更新。圖像域運動補償采用第4部分推出的式(11)，其中 vt通過傳感器實時獲取， t是單波束前視掃描聲吶的聲吶陣轉(zhuǎn)動的步進時間間隔，可以通過設(shè)備參數(shù)獲取， φ是通過傳感器實時獲取對底高度，并通過公式φ=arcsin(H/R)來獲得，其中 H為對底高度，R 為聲源與目標(biāo)點的距離。利用得到 ?d對上一次整個圖像整體往聲吶運動方向作搬移。然后再用圖8中更新的小扇形區(qū)的數(shù)據(jù)替換補償后圖像的相應(yīng)位置數(shù)據(jù)。然后輸出的就是經(jīng)過實時運動補償?shù)那耙晥D像。數(shù)據(jù)用的是海試實驗對海底觀測的數(shù)據(jù)。

圖9、圖10和圖11分別是潛水器連續(xù)航行時聲吶陣按逆時針轉(zhuǎn)動至33.6°, 0°和–33.6°時的前視成像。每一個位置都顯示了不補償、直接補償和本文補償3種成像結(jié)果，具體的偏差數(shù)值如表2所示。

從表2可以看出本文算法和直接補償算法都能夠補償前視聲吶前進造成的圖像偏差，本文算法比直接補償法效果更好。

表1 前視聲吶工作參數(shù)

圖8 前視圖像后置處理流程

圖9 逆時針轉(zhuǎn)動至33.6°時成像

圖10 逆時針轉(zhuǎn)動至0°時成像

圖11 逆時針轉(zhuǎn)動至-33.6°時成像

表2 前視聲吶成像偏差

7 結(jié)論

本文分析了運動平臺前視掃描聲吶成像的各種插值方法的優(yōu)缺點，在此基礎(chǔ)上提出了一種基于后向映射、雙線性插值和運動補償相結(jié)合的成像方法。通過分析載體運動，推導(dǎo)出圖像域徑向誤差和運動補償量，并用泰勒公式展開進行簡化。經(jīng)過仿真和海試數(shù)據(jù)驗證，經(jīng)過運動補償后的圖像，糾正了不同ping之間的徑向誤差，能準(zhǔn)確地表達(dá)目標(biāo)信息。