王風帆,馬永
(國家海洋信息中心,天津 300171)
海底地形是海洋科學研究、海洋工程和軍事應用等方面的重要基礎數據,開展海底地形測量是海洋測繪的重要內容[1-2]。多波束測深是當下海底地形測量的主要手段之一,可以對海底實現全覆蓋測量[3]。目前,海上多波束測量的測線布設一般借助于ArcGIS、Global Mapper 等GIS 平臺[4-5]。由于這些平臺缺少針對性功能開發,導致操作效率不高,設計人員需要同時掌握兩個領域的專業知識與技術,提高了熟練進行測線布設與修改的學習成本。
另外,在南北極等特殊區域,因常年漂浮海冰等障礙物[6-8],在此區域開展多波束測深存在一些特殊注意事項。例如在測量作業期間,因調查船躲避浮冰會導致測線出現大量彎曲;而船舶破冰產生的噪聲會造成波束腳印縮短,這些都會影響全覆蓋測量的實現。為了填補測量間隙,研究人員需要現場根據多波束數據預處理結果進行測線設計與修改。因此,開發專門的多波束測量測線布設系統,提高測線布設效率具有重要意義。
本文首先根據極地海域多波束水深測量作業需求,開展了系統總體設計、架構、功能設計與開發,然后闡述了功能實現的關鍵技術,最后選取區域進行應用試驗及討論。
多波束測量測線布設系統針對作業前實施方案編制和現場補測工作設計了兩種場景,主要功能包括測線、條幅的生成以及布設結果評價三個方面。測線生成主要分為根據預設參數自動生成和手動設計兩種方式,分別對應前面兩種場景;為了使生成條幅盡可能接近實際測量結果,提高預測準確性,引入海底地形網格數據等相關資料作為參考;對測量結果的評價需要依據相關標準指標[9],并綜合考慮實際場景以及功能實現的難易度。
根據極地海域多波束測量測線布設實際需求,系統技術路線(圖1)設計如下:
圖1 技術路線
(1)輸入測區邊界點地理坐標、多波束儀掃描開角、條幅重疊率指標等初始參數;確定測線方向,或者利用測區范圍內國際公開地形網格數據(先驗地形數據)計算獲得建議測線方向,后者主要用于主測線設計。
(2)針對主測線設計,根據步驟(1)參數,結合先驗地形數據,自動生成測線;針對檢查測線設計,通過輸入測線間隔自動生成測線;針對現場補漏測線設計,通過手動設計方式生成測線。
(3)結合先驗地形數據,沿測線生成多波束測量條幅多邊形。
(4)利用主測線、檢查測線和補漏測線的條幅多邊形以及測區坐標參數,統計總覆蓋率和主測線條幅重疊率等指標,并通過可視化技術展示。
系統使用Python 語言以及numpy、matplotlib、scipy、shapely、gdal等第三方函數庫,以PyCharm為開發平臺,采用分層模式進行架構[10](圖2),包括數據層、應用層和表示層。
圖2 系統架構
(1)數據層:包括支撐系統功能實現的基礎數據(先驗地形數據等)和系統運行所產生的產品數據,例如測線、條幅等。
(2)應用層:提供系統各項功能,實現各種操作,包括參數推薦、測線生成、條幅制作、結果評價等。
(3)表示層:為用戶提供交互界面,可輸入各項參數、顯示設計過程與結果、手動點擊生成測線,以及導入測線數據等。
(1)建議測線方向計算
利用測區坐標參數,制作測區多邊形(可以是非矩形)。選用IBCAO 等國際公開地形網格[11-12]數據源作為先驗地形數據,利用測區多邊形進行裁剪,獲得測區內地形網格數據。生成等深線,計算等深線走向,以獲得的所有等深線走向平均數作為測線走向建議值。
(2)測線生成
測線自動生成包括自動測線間隔與固定測線間隔兩種方式。前者以滿足條幅重疊率指標為目標,通過腳印估計的方式保證測線間條幅重疊率大致達到指標要求,這里需要用到先驗地形數據、多波束開角等;后者則以滿足測線間距要求為目標,固定測線之間距離,生成均勻排布的測線。
手動設計方式則是通過人機互動方式,手動點擊測區內位置生成測點,依次連接測點形成測線。
(3)條幅制作
輸入測線數據,包括自動生成測線和手動生成測線。利用先驗地形網格獲得沿測線水深,結合多波束開角等信息計算多波束腳印,利用腳印端點坐標生成條幅多邊形。
另外,為了進一步滿足現場設計需求,系統還需能夠導入現場多波束數據處理結果,并形成覆蓋區域多邊形,用于現場補漏設計參考及結果評價。
(4)結果評價
統計并檢查測線總長度,計算設計測線形成的總覆蓋率和條幅重疊率[9]。利用測區參數和條幅多邊形,計算測量覆蓋面積與測區面積之比,作為總覆蓋率。條幅重疊率主要針對主測線,統計相鄰測線重疊區域寬度與測線間距之比的平均值,即平均條幅重疊率γ:
式中:Li為采樣點處測線間條幅重疊的寬度,Di為采樣點處測線間距,n為采樣點數。
測線間條幅重疊率可通過可視化手段盡可能展示相關信息,低于預設指標的重疊區域可在圖中突出顯示。
(5)測線導出
完成設計后,將形成測線導出為shp 格式或者文本文件等多種形式,前者主要用于同其他制圖軟件配合,方便于方案圖件制作,后者是為了作業導航,用于與調查船導航系統配合。
系統測線自動生成功能有兩種方式:測線自動排布與固定間隔排布。二者主要區別在于測線間隔的確定方法,主測線排布需要滿足條幅重疊率指標的要求[9],特別是在地形變化比較劇烈的區域,如果采用等間隔測線容易造成成本與標準之間的矛盾,間隔太大容易造成重疊率不達標,甚至出現呈帶狀的測量間隙,間隔太小會導致工作量太大,成本太高,手工設計比較耗時,因此需要進行測線間隔的自動判斷,以提供最為合適的測線布設方案。而檢查測線沒有條帶覆蓋率的要求,因此可以采用效率較高的固定間隔方式生成。
(1)測線自動排布
測線走向確定后,再確定測線通過的一點即可生成一條直線,然后利用測區多邊形對直線切割即可得到測線(圖3)。
圖3 測線自動排布原理圖
這里為了便于實現,首先將測區多邊形簡化為矩形,然后在一條通過矩形頂點并與測線垂直的線段上確定測線點,這條線段稱為參考線。
通過插值獲取參考線上的水深,然后利用式(2)計算各水深點的單側腳印寬度W。
式中:θ為多波束開角,H為水深。
最后在測線上依次確定滿足條幅重疊率[9]的最遠的點,作為測線點。
(2)固定間隔排布
與測線自動排布方法類似,區別在于參考線改為簡化矩形的對角線。在輸入固定的測線間隔后,以一定間隔在參考線上進行測線點的排布。
最后,根據測線走向和錨定的點生成直線,再利用測區邊界多邊形進行切割,獲得測線。
輸入測線,統一轉換為LineString 格式,利用先驗地形網格數據,插值獲取測線上的水深信息,結合多波束儀開角,利用式(1)計算垂直測線方向的腳印長度,進而獲得腳印端點坐標,順時針依次連接測線兩側腳印端點,獲取條幅多邊形。
輸入的多波束預處理數據一般是geotiff 格式,利用gdal庫解析,獲取坐標、投影、水深等信息,進而獲取測量空白的位置,并生成覆蓋區域多邊形,用于補漏測線設計及結果評價。
結果評價功能中,比較關鍵的是根據測區、條幅和多波束覆蓋等形成的多邊形對總覆蓋率和條幅重疊率進行計算。
根據《海洋多波束水深測量規程》(DZT 0292-2016),多波束測深需滿足全覆蓋測量及達到一定的條幅重疊率。融合條幅獲取總覆蓋多邊形,該多邊形面積與測區面積之比達到100%即為全覆蓋測量。在前期設計階段,計算總覆蓋率所用到的主要是主測線和檢查測線條幅多邊形的融合結果,現場設計階段用到的主要是現場處理結果所形成的多邊形與檢查測線條幅多邊形的融合結果。
條幅重疊率按照規程為重疊寬度與測線間距的百分比,應用于對主測線設計的評估。系統采用可視化手段,對測線間條幅重疊率進行展示,方便設計者直觀評價設計結果。受地形變化影響,條幅一般不是規則矩形,測線間的條幅重疊率是變化的,因此采用重疊率平均值以及標準差作為評價依據。
(1)前期方案設計
選取某一海區進行試驗(圖4),測區形狀根據地形大致走向設計為平行四邊形,截取區域內水深網格數據,計算得到建議主測線布設方向為218.7°。設置條幅重疊率指標為10%。
圖4 方案設計主測線與檢查測線
通過自動測線排布生成主測線(圖4 中藍色線),通過固定間隔排布(24 km)生成檢查測線(圖4 中黑色線)。利用條幅制作功能生成主測線條幅(圖4中黃色條帶)與檢查測線條幅(圖4中綠色條帶),其中,黃色條帶間加深部分為條幅重疊區域。最后利用測線評估得到總覆蓋率(98%)和測線間條幅重疊率(圖5)。
圖5 方案設計主測線間條幅重疊率
可以看出,總覆蓋率基本能夠滿足全覆蓋要求,但也存在一些小的空白區域。而在實際測量中,調查船改換測線時這些區域順便就被覆蓋掉,因此影響不大。相鄰測線間平均條幅重疊率存在大于20%的情況,這是由于部分測線上水深變化較大引起的,這一點可以從圖5 中的標準差看出,但考慮到實際測量中調查船航跡存在不可避免的彎曲,相較而言對結果的影響不會很顯著。
(2)補漏測線設計
利用歷史多波束測量處理數據,生成多波束覆蓋區域圖(圖6),現場作業時為了躲避海冰等障礙物,測量產生了許多空白區域。通過點擊方式設計4 條補漏測線(圖7 中藍色線),繼而生成條幅多邊形(圖7 中黃色條帶),可以直觀看到測線條幅已基本將空白區域覆蓋掉,最后利用評價功能計算總覆蓋率為99.7%,基本實現全覆蓋。
圖6 現場處理多波束數據(白色區域為漏測區域)
圖7 補漏測線設計
試驗證明,多波束測量測線布設系統可以方便、高效地設計出多波束測線,并對設計結果進行快速評價,可應用于前期與現場設計等場景。系統還將繼續在諸多地方進行改進,以期為設計者提供一個功能齊全、使用方便的多波束測量測線設計工具。
本文針對極地海域多波束測量測線布設進行了系統設計,并對其中關鍵技術進行了詳細的分析探討,最后利用實際資料進行試驗與討論。研究結果表明:
(1)系統設計明確,架構簡捷有效,可方便用于系統功能開發,與用戶交互友好,設計人員上手難度小,適用于多波束測量測線布設。
(2)測線自動生成技術可快速、有效生成測線,制作條幅與可視化展示便于結果評價、信息統計對標相關規程,并針對現場補測場景開展了實用設計,方便高效。
(3)為進一步滿足多波束測量測線布設實際需求,系統會進行界面化設計,并進一步提升兼容性與擴展性,考慮繼續優化測線生成與腳印預測功能,增加外部測線輸入等功能,不斷提高系統的可用性。