基于全球航行的海洋氣象演示系統設計?

郭 鹍 趙甲文 陶冠時 蔣曉龍

（1.中船航海科技有限責任公司 北京 100070）（2.天津航海儀器研究所 天津 300131）

1 引言

利用傳統氣象環境系統進行輔助決策需要非常專業知識，隨著航運業的發展，岸海一體及智能航行技術呈現跨越性的進步，基于大數據氣象信息輔助決策已經由一種趨勢轉變為現實。通過衛星通信方式，傳輸基于grib2格式的氣象要素［1］，實時掌握世界范圍內的天氣動態，對航海人員的日程規劃、航海路線制定、貨物裝卸、室外活動起著非常大的輔助作用［2］。為了保證航海工作人員能以形象化的手段獲取所制定航線范圍內的天氣狀態，本文采用質點運動動畫算法以及標量渲染線性插值算法，形象地展示了全球范圍內海洋氣象環境各種氣象要素的質點運動或渲染［3］，最后形成海洋氣象演示系統，船舶駕駛人員可以通過該演示系統觀察周邊海域以及所關心海域的氣象信息，幫助其及時作出決策。

2 海洋氣象演示系統架構設計

海洋氣象演示系統采用B/S結構進行設計［2］，客戶端采用使用SWT瀏覽器組件進行設計。SWT瀏覽器組件提供了在Java應用程序中內嵌瀏覽器應用，并可通過Java API與瀏覽器自身進行交互，通過在客戶端程序內嵌一個內置瀏覽器完成軟件渲染和演示［7，15］。

海洋氣象演示系統所解析的氣象數據源來自于美國國家環境預測中心全球預報系統，其數據格式為國際氣象組織規定的標準文件格式grib2，對grib2編碼進行分析、解析其中的相關氣象數據，然后分類處理并轉換為Json數據格式。Json數據比較方便數據庫應用，B/S客戶端將Json數據渲染到瀏覽器頁面。

全球海洋氣象演示系統在啟動之后首先會到指定網址上下載grib2數據，存放到指定的目錄下，也可以直接讀取從外部磁盤導入的grib2文件。當解析文件時，則調用Grib2json系統進行解析，解析完成之后，全球海洋氣象演示系統會把json文件存放到自己的publicdata路徑下，系統按照類別（大氣或海洋、時間、分辨率）去讀取合適的json文件來進行渲染，最終在高清屏幕上以可視化的方式進行展示，包括風速、渦旋、垂直速率、溫度、相對濕度、云混合率、可降水量、云水量、洋流等，其中，風速矢量、洋流矢量需要以動畫形式展示出來。并通過按鈕對氣象要素按照海拔、分辨率、時間來進行自由切換，完成各類氣象數據的瀏覽、漫游、查詢、顯示功能，并根據不同類型的氣象數據提供不同渲染模式。

海洋氣象演示系統中通信和信息傳輸主要依靠移動通信網絡和移動衛星網絡進行保障，移動衛星網絡系統采用VSAT通信或者Inmarsat通信系統［9～10］。Inmarsat通信系統是由國際海事衛星組織（Inmarsat）自1970年代開始發展的衛星移動通信系統，可覆蓋全球除兩極區之外的全部地區。最高速率可達 492 kbit/s［13～14］，其成本相對較高。考慮成本問題，在離岸較近的情況下，氣象數據下載會選用陸地的3G/4G網絡，在遠洋時優先使用VSAT通信以及Inmarsat通信系統提供的衛星網絡［12］。

全球海洋氣象演示系統的工作流程如圖1所示。

圖1 全球海洋氣象演示系統工作流程

3 數據下載及格點數據解析

海洋氣象演示系統通過FTP方式連接到美國國家環境預測中心獲取全球預報系統數據，該數據提供了4個分辨率：2.5度分辨率、1度分辨率、0.5度分辨率、0.25度分辨率，海洋氣象演示系統可根據網絡狀態以及實際網絡下行速率選用不同數據分辨率，因數據較大需要采用斷點續傳，在條件允許的情況下一般選用0.25度分辨率的數據，數據大約197MB，在遠洋航行時，由于通信成本可采用2.5度分辨率數據，數據約3.7MB。

本文中的海洋氣象演示系統的數據為grib2格式，grib2是國際氣象組織的最新標準格式［1］，grib2具備更好的自描述性、可擴展性和靈活性，與以往的氣象標準數據格式相比如GRID或者GARF相比是一種更為方便的傳輸大容量數據格式，我國氣象預報分析工作也采用該格式數據作為參考，歐洲格點報、地面繪圖報、正點實況報均采用該格式。研究掌握grib2的文件數據分析使用可進一步拓展全球海洋環境信息演示系統應用范圍。

grib2的文件內容是由十六進制字符（alphanumeric character）組成的。比如：FF表示十六進制的FF，數值為255，并且是按照字節依次存放的。每2個字節之間會有空格作為分割符號。grib2文件是有若干的區段組成的文件，每個區段表示一定的含義，存放一定的數據。

grib2的區段描述如表2所示：

表2 grib2區段描述

grib2具有多種優勢，尤其在多維數據描述方面，因此美國NOAA的數據包中存儲了多個格點數據、多個產品，具備相當的模塊特性，而且為擴展和增加提供方便的途徑［11］。

4 頁面渲染設計

海洋氣象演示系統頁面渲染分別為氣象要素標量渲染和矢量渲染，根據不同要素渲染選取色彩，通過不同顏色色條的深淺來表達值域大小，非常直觀。

4.1 標量渲染

海洋氣象演示系統的溫度、濕度均為標量數據，在渲染圖上的每個網格都對應著一個渲染點，根據數據值的范圍渲染出對應的顏色。由于數據采樣點是離散的，因此溫度和濕度等信息需要采用算法進行插值。

標量渲染采用了單線性線性插值［4～5］，基本思想是根據點（x，y）的2個鄰點灰度值，經過線性插值計算求出（xc，yc），用數學方法來描述輸入、輸出圖像點之間的空間關系時，其幾何運算的一般定義為

這種算法是逐象素、逐行地產出輸出圖像。每個象素的灰度級由2個象素參與的插值所唯一確定。這種方法對于氣象標量要素渲染邏輯執行更為有效。

如圖2所示，A區域渲染點可以根據從json中提取的數據進行顏色渲染，渲染為天藍色，B區域渲染點無值，C區域渲染點同樣進行了顏色渲染，渲染為綠色，那么，此時B區域處的每個渲染點就要利用上述雙線性插值的方法進行顏色的過渡渲染，使得A區域渲染點和B區域渲染點之間的交匯時的顏色柔和自然。

圖2 氣象標量要素渲染示意圖

4.2 矢量渲染方式

矢量渲染用于風速和洋流的渲染，而矢量渲染則會有兩個值，由這兩個值確定一條方形的對角線。

矢量的渲染與其他氣象要素的渲染有所不同，其他分量用標量渲染時，每個點都有一個對應的渲染值，而矢量渲染則會有U（U-Component）V（V-Component）兩個值，由這兩個值確定一條方形對角線。實際情況中，如果存在的每個矢量都被描繪出來，則會有非常多的細微的線條，從效果上來看，就會混沌不清，不易辨識。此時需要從平行的對角線中抽象出一條，沿著此條對角線的方向繼續抽象渲染，用來清晰準確地表達矢量走向。并在此線上進行粒子描繪，形成加速動畫。

矢量的渲染采用了雙線性插值［4～5］進行線段抽象，基本思想是根據點（x，y）的4個鄰點灰度值，經過線性插值計算求出（x0，y0），用數學方法來描述輸入、輸出圖像點之間的空間關系時，其幾何運算的一般定義為

雙線性插值是逐象素、逐行地產出輸出圖像。每個象素的灰度級由四個象素參與的插值所唯一確定。這種方法對于氣象標量要素渲染邏輯執行更為有效。如圖3所示。

圖3 雙線性插值示意圖

具體f（x，y）在點P（x，y）的插值算法如下：假設已知函數 f在 P11=（x1，y1）、P12=（x1，y2），P21=（x2，y1）以及 P22=（x2，y2）四個點的值。首先在 x 方向進行線性插值，得到

然后在y方向進行線性插值，得到

最后綜合起來就是雙線性插值最后的結果：

通過線性插值得到抽象線段F。

以圖4下抽象線段F為例。

圖4 矢量動畫運行圖示意

此線段抽象出了幾個主要的點，并在此線段上進行粒子描繪，亮度通過不同灰度值的白色定義，數值越大則越接近白色，步驟如下：

第一步：P1（x1，y1）點高亮顯示，假設灰度值為10。

第二步：P2（x2，y2）點亮度值變為 10，P1（x1，y1）點變為稍暗的灰度值顯示（灰度值值為5）。

第三步：P3（x3，y3）點灰度值值為 10，P2（x2，y2）點灰度值值為5，P1（x1，y1）點灰度值值為3。

第四步：P4（x4，y4）點灰度值值為 10，P3（x3，y3）點灰度值值為 5，P2（x2，y2）點灰度值值為 3，P1（x1，y1）點灰度值為1。

由上幾步，可以顯示一個高亮的粒子拖曳著尾巴（顯示劃過的痕跡）從抽象線段上劃過，通過亮度的變化并無限循環，就形成了風速的質點動畫。實際處理時可采用整數進行運算，本文采用了整數計算代替浮點數運算，比如左移11位變11位小數為整數，頁面渲染實例如圖5所示。

圖5 頁面渲染示意圖

5 軟件整體運行

海洋氣象演示系統其他要素的處理方法也都與風、溫度等要素處理辦法類似，這里不再描述，在本系統完成數據下載、解析后，通過渲染后的工作界面如圖6所示。

6 結語

本文首先對基于全球航行的海洋氣象演示系統設計架構進行了說明，根據海洋氣象演示系統分析了grib2數據。在完成grib2數據解析的基礎上，針對溫度、濕度數據采用線性插值渲染，針對風速、風向矢量運動的特點，設計質點動畫算法。該演示系統綜合運用了計算機通訊技術、圖像處理技術，顯示全球海域的氣象狀態清晰、簡明并易于學習，一方面可為船舶系統電子配套設計人員作為開展智能船舶設計參考，另一方面可輔助航海部門人員知悉航行過程中海洋環境各個氣象要素的變化情況，幫助船舶進行航路規劃，規避航路上危險天氣。