許墅，戴文伯，周雨淼，季明

（中交疏浚技術(shù)裝備國家工程研究中心有限公司，上海 200082）

0 引言

國外于20 世紀(jì)80 年代提出數(shù)字化的疏浚監(jiān)視系統(tǒng)，提高了疏浚船舶的施工效率[1]。國外以IHC 公司為代表，其疏浚施工監(jiān)控系統(tǒng)可大大提高船舶航行與作業(yè)的安全性、可靠性和經(jīng)濟(jì)性，國內(nèi)的中交疏浚研究中心等單位都自主研發(fā)了疏浚施工監(jiān)視系統(tǒng)，包含船舶定位、剖面顯示、施工狀態(tài)監(jiān)控以及歷史回放等功能，在國內(nèi)眾多挖泥船上得到廣泛應(yīng)用[2-5]。

目前的疏浚施工監(jiān)視系統(tǒng)大多是二維系統(tǒng)，缺少三維施工監(jiān)視功能，導(dǎo)致系統(tǒng)監(jiān)視效果不夠直觀形象，而且系統(tǒng)都是基于C/S（Client/Server）架構(gòu)，導(dǎo)致每臺客戶機(jī)都需安裝客戶端軟件，并進(jìn)行調(diào)試，當(dāng)有多臺客戶機(jī)時(shí)，安裝部署和維護(hù)較為麻煩，當(dāng)有跨平臺需求時(shí)，對不同操作系統(tǒng)需要單獨(dú)開發(fā)客戶端軟件[6]。而基于B/S 架構(gòu)的系統(tǒng)在處理多終端需求時(shí)，只需要在服務(wù)器上部署系統(tǒng)，更新較為方便，客戶端只需普通網(wǎng)頁瀏覽器即可訪問系統(tǒng)[7]。隨著網(wǎng)絡(luò)通信技術(shù)的發(fā)展，瀏覽器業(yè)務(wù)處理能力逐漸增強(qiáng)，基于B/S 架構(gòu)的監(jiān)控和管理系統(tǒng)已在多個(gè)領(lǐng)域得到應(yīng)用[7-8]，而隨著WebGL 的逐漸成熟，也讓三維Web 開發(fā)成為未來的發(fā)展趨勢[8]。

針對目前耙吸船施工監(jiān)視系統(tǒng)存在的問題，本文設(shè)計(jì)了一套基于B/S 的耙吸船施工三維監(jiān)視系統(tǒng)，利用WebGL 技術(shù)將水下地形和耙吸挖泥船模型在網(wǎng)頁中三維顯示，并根據(jù)實(shí)時(shí)獲取的傳感器數(shù)據(jù)驅(qū)動船體和耙臂模型運(yùn)動，同時(shí)支持通過實(shí)時(shí)計(jì)算耙頭與水下地形接觸時(shí)的耙頭位置，對地形進(jìn)行重構(gòu)并顯示。該系統(tǒng)為疏浚作業(yè)提供可交互的三維視角，為高精度施工提供新的工具，讓疏浚人員直觀了解施工區(qū)域水深情況和疏浚軌跡。

1 系統(tǒng)概況

1.1 系統(tǒng)架構(gòu)

系統(tǒng)采用B/S 架構(gòu)，分為數(shù)據(jù)庫、服務(wù)器和瀏覽器3 個(gè)部分（如圖1 所示）。

數(shù)據(jù)庫存放著耙吸船的船體和耙臂3D 模型、實(shí)時(shí)數(shù)據(jù)、水深數(shù)據(jù)以及由水深數(shù)據(jù)在服務(wù)器端生成的地形數(shù)據(jù)。服務(wù)器負(fù)責(zé)響應(yīng)前端瀏覽器的請求并做處理，返回3D 模型數(shù)據(jù)、地形數(shù)據(jù)和施工數(shù)據(jù)，其中地形數(shù)據(jù)在服務(wù)端生成，利用地形構(gòu)建算法，由水深數(shù)據(jù)構(gòu)建生成，并保存到數(shù)據(jù)庫中，此外，前端瀏覽器會將地形重構(gòu)的結(jié)果發(fā)到服務(wù)器端，由服務(wù)器負(fù)責(zé)將地形重構(gòu)結(jié)果進(jìn)行保存。前端瀏覽器根據(jù)從服務(wù)器獲取的船舶3D模型實(shí)現(xiàn)船舶3D 顯示，并添加水面模擬效果；根據(jù)從服務(wù)器獲取到的施工數(shù)據(jù)驅(qū)動船體模型運(yùn)動，實(shí)時(shí)顯示船體方位、耙臂角度、耙頭姿態(tài)、吊架姿態(tài)等信息；根據(jù)從服務(wù)器獲取的地形數(shù)據(jù)構(gòu)建水下地形并渲染顯示，同時(shí)根據(jù)耙頭位置實(shí)時(shí)重構(gòu)地形，模擬挖掘效果，同時(shí)向服務(wù)器返回地形重構(gòu)結(jié)果。

1.2 系統(tǒng)功能

本系統(tǒng)主要包含場景交互、施工可視化和場景模擬三大功能。場景交互實(shí)現(xiàn)用戶對虛擬場景的交互操作，包含視角控制、視點(diǎn)切換和船體透明等子功能，其中視角控制可實(shí)現(xiàn)移動觀察、靜止觀察和圍繞船體3 種漫游模式的自由切換，以適應(yīng)不同場景下漫游查看的需要。移動觀察模式下用戶可自由調(diào)整視角，且視角跟隨船體移動，從而便于觀察船舶姿態(tài)變化，防止船舶移出畫面；靜止觀察模式下用戶可自由調(diào)整視角，但視角不跟隨船體移動，從而便于觀察船舶運(yùn)動狀態(tài)和周圍地形；圍繞船體模式下，視角以船為中心繞船移動，且跟隨船體移動，從而方便觀察船體四周情況。視點(diǎn)切換功能以畫中畫形式將小窗口視景快速定位到指定的視角位置，如耙頭、泥艙等位置，從而方便用戶觀察特定區(qū)域的狀態(tài)。施工可視化實(shí)現(xiàn)耙吸船疏浚過程的可視化監(jiān)控，基于實(shí)船數(shù)據(jù)驅(qū)動場景仿真模擬，包含航行模擬、吃水模擬、耙臂姿態(tài)模擬、地形挖掘等功能。場景模擬實(shí)現(xiàn)耙吸船施工場景的再現(xiàn)，包含船舶3D 顯示、水面模擬、地形可視化和光照模擬等。整個(gè)系統(tǒng)功能如圖2 所示。

圖2 系統(tǒng)功能圖Fig.2 System function diagram

2 主要技術(shù)問題及解決方案

2.1 數(shù)字地形模型構(gòu)建

通過多波束等聲吶設(shè)備采集到的現(xiàn)場施工區(qū)域水深數(shù)據(jù)是離散點(diǎn)形式，數(shù)據(jù)點(diǎn)間距一般各不相同，適合以不規(guī)則三角網(wǎng)（TIN）構(gòu)成面片，因此，本文采用Delaunay 方法構(gòu)建地形三角網(wǎng)。先將三維空間的水深離散點(diǎn)集投影到XY平面，對3D 離散點(diǎn)集進(jìn)行降維操作，變?yōu)?D 平面離散點(diǎn)集，將3D Delaunay 三角化變?yōu)?D Delaunay 三角化，降低了三角化的復(fù)雜性。然而Delaunay 方法具有凸包特性，針對凹包形態(tài)離散點(diǎn)會生成多余三角網(wǎng)[9]，圖3—圖5 說明凹包形態(tài)離散點(diǎn)采用Delaunay 方法構(gòu)網(wǎng)存在的問題，圖3 是凹包形態(tài)的原始離散點(diǎn)，圖4 是期望的正確三角構(gòu)網(wǎng)結(jié)果，圖5 是采用Delaunay 方法三角構(gòu)網(wǎng)的錯(cuò)誤結(jié)果，其中虛線部分三角形是生成的多余三角網(wǎng)，從而使得三角網(wǎng)邊界形態(tài)與原始離散點(diǎn)邊界形態(tài)不一致，導(dǎo)致地形構(gòu)建的區(qū)域超過原始離散點(diǎn)覆蓋的區(qū)域，影響地形構(gòu)建的準(zhǔn)確性和后期渲染的美觀性。施工現(xiàn)場采集到的水深離散點(diǎn)邊界不規(guī)則，形態(tài)各異，因此采用Delaunay 方法必須首先要解決生成錯(cuò)誤網(wǎng)格的問題。

圖3 離散點(diǎn)Fig.3 Discrete point

圖4 正確構(gòu)網(wǎng)結(jié)果Fig.4 Correct triangulation results

圖5 錯(cuò)誤構(gòu)網(wǎng)結(jié)果Fig.5 Error triangulation results

針對錯(cuò)誤網(wǎng)格問題，本文通過在Delaunay 三角化中引入三角網(wǎng)索引和拓?fù)鋵傩裕谏傻娜蔷W(wǎng)基礎(chǔ)上借助三角形鏈表和鄰邊關(guān)系，通過遍歷三角網(wǎng)最外側(cè)三角形，即鄰邊三角形個(gè)數(shù)小于3 的三角形，從而剔除邊長超過閾值的三角形，同時(shí)更新三角形鏈表和鄰邊關(guān)系，然后重新循環(huán)此過程，直至最外側(cè)三角形邊長都小于閾值才結(jié)束，進(jìn)而優(yōu)化三角網(wǎng)形態(tài)。圖6 是數(shù)字地形模型的完整構(gòu)建過程。

圖6 數(shù)字地形模型構(gòu)建過程Fig.6 Construction process of digital terrain model

本文邊長閾值根據(jù)計(jì)算的離散水深點(diǎn)平均間距而確定，即將平均間距乘以1 個(gè)放大系數(shù)后作為閾值，經(jīng)過測試，系數(shù)取1.5~2.0 比較合適。圖7 是未優(yōu)化的錯(cuò)誤三角網(wǎng)效果，圖中方框區(qū)域可明顯看出存在的錯(cuò)誤網(wǎng)格，圖8 是優(yōu)化后的正確三角網(wǎng)效果。

圖7 錯(cuò)誤三角網(wǎng)效果Fig.7 Error triangulation effect

圖8 正確三角網(wǎng)效果Fig.8 Correct triangulation effect

2.2 地形重構(gòu)

耙吸挖泥船在疏浚施工過程中耙頭與地形相互作用，使得地形實(shí)時(shí)發(fā)生變化，耙頭作用區(qū)域地形會發(fā)生塌陷，形成挖槽，因此，耙吸船施工三維監(jiān)視需要研究實(shí)時(shí)的地形重構(gòu)技術(shù)，實(shí)現(xiàn)耙頭挖掘地形的模擬。針對地形的動態(tài)修改，本質(zhì)上是對地形幾何頂點(diǎn)空間位置的修改，因此首先需要獲得耙頭挖掘影響區(qū)域的頂點(diǎn)，如圖9 所示，然后對這些頂點(diǎn)位置進(jìn)行修改。由于需要實(shí)時(shí)展現(xiàn)耙頭挖掘地形的效果，但地形空間范圍廣，頂點(diǎn)數(shù)量多，因此需要一個(gè)高效的算法快速查找到耙頭挖掘過程中經(jīng)過區(qū)域的地形頂點(diǎn)。

圖9 耙頭影響區(qū)域內(nèi)的頂點(diǎn)Fig.9 Vertices within the influence area of the rake head

本文采用BVH（Bounding Volume Hierarchies，層次包圍盒）方法實(shí)現(xiàn)高效快速獲得耙頭影響區(qū)域的頂點(diǎn)。BVH 方法的核心思想是用體積略大而幾何特征簡單的包圍盒來近似地描述復(fù)雜的幾何對象，所有包圍盒分層逐次地再次包圍，獲得一個(gè)更大的包圍盒，直到包圍住所有物體，實(shí)際上，它是一個(gè)樹形結(jié)構(gòu)，因此可以仿照樹的結(jié)構(gòu)，將2 個(gè)或3 個(gè)小的包圍體包圍成1 個(gè)更大的包圍體，以此類推。圖10 是BVH 方法構(gòu)建層次包圍盒的示意圖。

圖10 層次包圍盒示意圖Fig.10 BVH schematic

本文在地形構(gòu)建過程中采用BVH 方法同步構(gòu)建地形層次包圍盒，將地形劃分為若干包圍盒包圍的細(xì)小區(qū)域，從而對耙頭與相交的區(qū)塊做進(jìn)一步的相交測試，大大加速耙頭與地形相交測試的過程，更快獲得相交區(qū)域的地形頂點(diǎn)。圖11 是建立的地形層次包圍盒效果圖。

圖11 地形層次包圍盒Fig.11 Terrain BVH

2.3 數(shù)據(jù)壓縮

B/S 架構(gòu)下，前端WebGL 程序需要向服務(wù)端請求大量三維場景數(shù)據(jù)，而在瀏覽器中顯示高精度三維場景易受到網(wǎng)絡(luò)帶寬的限制，為將加載時(shí)間控制在可接受的范圍內(nèi)，從而縮短響應(yīng)時(shí)間，提升用戶體驗(yàn)感，本文使用Draco 壓縮算法對服務(wù)端存放的耙吸船模型和地形等數(shù)據(jù)進(jìn)行壓縮，后傳輸給前端解壓并使用，經(jīng)過測試，耙吸船模型數(shù)據(jù)壓縮率達(dá)到50%左右，地形數(shù)據(jù)平均壓縮率為70%左右。

3 實(shí)現(xiàn)效果

系統(tǒng)后端針對上傳的水深文件進(jìn)行地形構(gòu)建處理，生成地形模型，前端三維窗口初始化時(shí)加載船體模型和地形模型，并根據(jù)從后臺獲取耙吸船的船位、航向、潮位、耙臂姿態(tài)角度等數(shù)據(jù)在三維場景初始化船體狀態(tài)，最終利用WebGL 技術(shù)渲染顯示船體、水面和地形，其中地形根據(jù)設(shè)定的水深顏色表進(jìn)行分級顯示，同時(shí)在地形上疊加上光照和陰影，使得地形渲染效果更有立體感。三維場景初始化完成后，前端從后臺定時(shí)獲取施工數(shù)據(jù)，從而在三維場景中動態(tài)更新船舶施工運(yùn)行狀態(tài)，包括船舶航行軌跡、耙臂和耙頭姿態(tài)等，同時(shí)利用地形重構(gòu)技術(shù)，結(jié)合耙頭位置信息，實(shí)現(xiàn)耙吸船疏浚過程水下泥面的實(shí)時(shí)動態(tài)更新。

圖12 是前端系統(tǒng)界面圖，在頁面上方的菜單欄中可以顯示挖泥船和周圍環(huán)境的一些基本狀態(tài)信息，如航速、船體吃水、航向、GPS 位置、潮位等，并有能控制船體透明開啟或關(guān)閉的切換按鈕，以及視角控制模式切換的按鈕，視圖窗口右側(cè)是2 個(gè)小窗口畫面，方便用戶觀察特定區(qū)域的場景畫面。

圖12 前端系統(tǒng)界面圖Fig.12 System interface diagram

圖13 是利用系統(tǒng)實(shí)現(xiàn)的耙吸船疏浚過程三維監(jiān)視效果，用戶可直觀了解施工區(qū)域水深情況和疏浚軌跡，實(shí)時(shí)查看水下地形動態(tài)重構(gòu)效果，并能直觀查看耙臂實(shí)時(shí)姿態(tài)。

圖13 耙吸船疏浚過程三維監(jiān)視效果Fig.13 3D monitoring effect of TSHD during dredging process

經(jīng)測試，系統(tǒng)在處理器I7、顯卡NVIDIA GTX1080、內(nèi)存16G 的電腦上運(yùn)行，渲染幀速達(dá)到30FPS 以上，可流暢操作，而且利用相同的施工數(shù)據(jù)同步驅(qū)動二維施工監(jiān)視系統(tǒng)和三維監(jiān)視系統(tǒng)，測試結(jié)果表明，三維監(jiān)視系統(tǒng)的船舶位置、耙臂姿態(tài)和施工軌跡等與二維監(jiān)視系統(tǒng)保持一致，滿足施工控制要求。圖14 是二維施工監(jiān)視系統(tǒng)和三維施工監(jiān)視系統(tǒng)同步對比的效果。

圖14 二維和三維施工監(jiān)視系統(tǒng)同步對比Fig.14 Synchronous comparison of 2D and 3D construction monitoring systems

4 結(jié)語

針對基于C/S 架構(gòu)的二維施工監(jiān)視系統(tǒng)存在的不足，本文設(shè)計(jì)并開發(fā)了一套基于B/S 架構(gòu)的耙吸船三維施工監(jiān)視系統(tǒng)，解決了數(shù)字地形模型構(gòu)建、地形重構(gòu)和數(shù)據(jù)壓縮等關(guān)鍵技術(shù)難題。系統(tǒng)采用WebGL 技術(shù)建立了三維可視化系統(tǒng)，并接入施工數(shù)據(jù)，從而對當(dāng)前船體運(yùn)行狀態(tài)、耙臂姿態(tài)等進(jìn)行三維動態(tài)展示，并實(shí)現(xiàn)耙吸船疏浚過程水下泥面實(shí)時(shí)三維動態(tài)更新。系統(tǒng)開發(fā)成果為耙吸船疏浚作業(yè)提供全新和可交互的三維視角，將助力疏浚作業(yè)高質(zhì)量施工。