高潤喜

(包頭鐵道職業(yè)技術學院，內(nèi)蒙古 包頭014060)

0 引 言

海上原油泄漏近年來頻頻發(fā)生，造成了較為嚴重的環(huán)境污染和經(jīng)濟損失。為了避免此類災難事件影響范圍的持續(xù)擴大，在事故發(fā)生之后，應當立即向有關個人和公眾告知此類事件的嚴重程度和影響范圍，從而能夠及時采取相應的防范措施，使得負面影響能夠降到最低程度。因此，及時監(jiān)視和發(fā)布此類災難事件的地理位置等信息十分重要。由于地理信息系統(tǒng)(GIS)在地理信息分析上具有較為強大的功能，因而GIS 被廣泛應用于記錄和分析災難事件的分布、區(qū)域等地理和特征信息，并在最近一些實際應用中取得了良好效果。

自20 世紀60 年代出現(xiàn)以來，GIS 已經(jīng)被應用到多個領域。目前出現(xiàn)的一些商用GIS 軟件，具備強大功能和高昂價格，能夠滿足多種組織、公司、團體等對于地理信息記錄和分析的需求。然而對于大多數(shù)實際應用來說，此類商用軟件的功能過于復雜，同時擁有較為昂貴的使用成本［1－2］。

在實際應用中，大多數(shù)情況僅僅需要GIS 的部分功能即可以滿足需求，然而采用商用GIS 軟件卻不得不為多數(shù)不需要的功能付出高昂的代價。因此，當前業(yè)界的趨勢為采用開源GIS (OpenGIS)系統(tǒng)為平臺，根據(jù)現(xiàn)實實際需求，開發(fā)適用的定制化GIS系統(tǒng)。

OpenGIS 的意義為開源的地理信息系統(tǒng)，即開發(fā)者能夠按照一定的協(xié)議獲得開源軟件的源代碼，并根據(jù)自己的需要進行下載、修改、應用和發(fā)布。此類地理信息系統(tǒng)的開發(fā)往往依靠于個人開發(fā)者的不斷改進，而非企業(yè)團隊的系統(tǒng)開發(fā)。其最顯著的特征在于免費、自由，但同時開源軟件的功能也受到了較大限制［3］。

常見的OpenGIS 系統(tǒng)有MapWinGIS、DotSpatial等，其本質是一組開源的. NET 開發(fā)組件庫。本文采用DotSpatial 為平臺，設計和開發(fā)了一種海上信息分析和發(fā)布系統(tǒng)，并采用海上原油泄漏監(jiān)視信息分析發(fā)布為應用場景，對該系統(tǒng)中信息分析和發(fā)布的關鍵技術進行研究，同時對GIS 軟件存在的問題提出了相應的解決方案，并且通過實例證明本文提出方法具有較好的可行性和實用性。

1 應用分析

本文應用場景為海上原油泄漏信息的分析和發(fā)布，因而本文的目標在于構建一種有效的海上原油泄漏信息分析發(fā)布系統(tǒng)，采用DotSpatial 作為開發(fā)平臺，并采用GDAL 作為DotSpatial 的補充。

1.1 開發(fā)模式

MapWinGIS 是一組開源的開發(fā)工具，用于中等規(guī)模的 GIS 應用開發(fā)，其核心是一組稱為MapWinGIS 的ActiveX 控件。通過這些工具，開發(fā)者能夠向自己的系統(tǒng)中添加多種GIS 特性。作為MapWinGIS 的升級版本，DotSpatial 擁有更加強大的功能，同時能夠支持更多的開發(fā)環(huán)境，因此，本文采用DotSpatial 作為系統(tǒng)開發(fā)平臺。

地理空間數(shù)據(jù)抽象庫 (Geo Spatial Data Abstraction Library，GSDL)是一個開源庫［4］，能夠實現(xiàn)對于柵格化地理數(shù)據(jù)的讀取、轉換和編輯。當前有許多開源軟件，甚至一些商用軟件均使用該庫。GSDL 擁有自己的數(shù)據(jù)模型和API，并且提供一系列命令行用來編譯和處理數(shù)據(jù)［5］。

1.2 系統(tǒng)功能結構

海上原油泄漏信息監(jiān)視、分析和發(fā)布系統(tǒng)包含地圖顯示、信息展示、數(shù)據(jù)統(tǒng)計、數(shù)據(jù)映射和輸出、公告發(fā)布等功能。其中核心功能為數(shù)據(jù)顯示、編譯、映射和輸出。

柵格化的地理數(shù)據(jù)與傳統(tǒng)的矢量數(shù)據(jù)具有較大不同，當需要在DotSpatial 電子地圖中顯示此類數(shù)據(jù)時，諸如重投影等數(shù)據(jù)處理過程應當事先執(zhí)行完畢，然后才能實現(xiàn)柵格化數(shù)據(jù)和矢量數(shù)據(jù)的分別顯示以及2 種數(shù)據(jù)的相互匹配。在本系統(tǒng)中，柵格化的原油泄漏點跡被轉化為矢量數(shù)據(jù)，然后數(shù)據(jù)被數(shù)值化、顯示，之后每個點跡的長軸和短軸均被確定，從而確定原油泄漏區(qū)域的范圍。這些處理能夠向用戶直觀地展示泄漏事件發(fā)生地點和發(fā)展趨勢。當需要將原油泄漏地點及信息映射到電子地圖并輸出時，用戶使用鼠標拖拽一個矩形，定義需要查看的區(qū)域，然后以HTML 格式發(fā)布相應結果。圖1 為系統(tǒng)的總體結構。

系統(tǒng)的功能總體包含5 個部分，其中地圖加載與顯示、數(shù)據(jù)編譯與分析、數(shù)據(jù)映射與輸出為系統(tǒng)的核心功能。

圖1 系統(tǒng)基本結構圖Fig.1 The structure of system

2 關鍵技術研究

在開發(fā)海上信息分析發(fā)布系統(tǒng)過程中，仍然有若干問題需要克服，包括地理空間數(shù)據(jù)的匹配、數(shù)據(jù)格式轉換、矢量數(shù)據(jù)顯示、映射和輸出，網(wǎng)頁設計和制作等。本文針對以上若干問題進行研究，并給出相應的解決方案。

2.1 數(shù)據(jù)顯示與控制

DotSpatial 在處理和顯示矢量數(shù)據(jù)時具有較好的性能和功能，而在處理柵格化數(shù)據(jù)時功能較弱。例如，DotSpatial 只能處理大小在40 M 以下的數(shù)據(jù)，因此，較大的圖片無法被添加到地圖中去。同時，由于柵格化數(shù)據(jù)與矢量數(shù)據(jù)有較大不同，DotSpatial也無法正常地讀取柵格化數(shù)據(jù)的坐標系信息。

除此之外，DotSpatial 在將圖片加入電子地圖時，無法對圖片進行拉伸，因此當電子地圖的網(wǎng)格數(shù)值不規(guī)則時，地圖圖片就無法正常反應不同網(wǎng)格之間的差異。

本文采用的圖片為ENVISAT SAR 圖片，其大小超過100 M。因此，使用GDAL 進行重采樣、重投影，并在圖片加入系統(tǒng)之前進行拉伸。由于重采樣會降低圖片的分辨率，因此需要將圖片的大小減少到40 M 以下。重投影過程由GDALWARP 實現(xiàn)，能夠實現(xiàn)圖片與地圖層次的匹配。通過以上步驟，實現(xiàn)地圖圖片與電子地圖的匹配，之后系統(tǒng)采用合理的灰度進一步強化地圖的視覺效果，在本系統(tǒng)中采用256 度灰度。

2.2 矢量數(shù)據(jù)提取與顯示

在對原油泄漏點跡進行分析的過程中，有多種算法能夠提取點跡的地理空間信息并進行計算，然后對多時態(tài)數(shù)據(jù)進行顯示。其中對于原油泄漏點跡的處理步驟為:

1)提取原油點跡的空間地理結構。泄漏原油的擴散范圍能夠通過點跡形狀信息反應出來，如幾何形態(tài)等，但是每個原油點跡的形狀總是不規(guī)則的。為了能夠直觀地長時間顯示每個原油點跡移動的變化情況，必須通過一些數(shù)據(jù)來表示其空間形狀，如區(qū)域信息、中點坐標、長軸、短軸、斷點坐標等。通過計算這些數(shù)據(jù)，能夠較為精確地反映出原油點跡的空間地理結構和幾何學形態(tài)，從而為判斷其運動趨勢和擴散范圍提供依據(jù);

2)對數(shù)據(jù)進行轉換。原油點跡信息是由RS 圖像轉換和提取而來的，其數(shù)據(jù)結構為柵格化數(shù)據(jù)，并不適合對空間結構信息進行分析。因此，將柵格化的地理空間數(shù)據(jù)轉換為矢量數(shù)據(jù)，從而更加便于表示每個原油點跡的空間形態(tài)。由于DotSpatial 并不具備該數(shù)據(jù)轉換的功能，因此采用GDAL 將SAR 圖像中的原油點跡圖形轉換為多邊形類型的幾何數(shù)據(jù)集。

3)計算泄漏原油區(qū)域的長軸和短軸。通過步驟2，將泄漏原油的擴散區(qū)域表示為一個多邊形。那么長軸即為多邊形邊界上距離最遠的2 個點跡之間的連線，而短軸則是垂直于長軸，同時被多邊形截取的線段部分。長軸和短軸共同定義了原油泄漏區(qū)域的輪廓，其中對于長軸和短軸的表示方法分別如圖2 和圖3 所示。

圖2 區(qū)域長軸的表示方法Fig.2 The major axis of the area

圖3 區(qū)域短軸的表示方法Fig.3 The minor axis of the area

4)連接各個原油點跡。假設原油泄漏區(qū)域共有n 個點跡，分別為P(0)，P(1)，P(2)，…，P(n)。將每個點跡分別與其他點跡互聯(lián)，則可以得到以下線段:

5)表示出長軸和短軸。通過步驟4 得到若干線段，通過比較線段的長度，假設其中最長的一條連接了P(i)和P(j)，那么L(i，j)即為區(qū)域的長軸。將長軸等分為200 份，在每一等份的邊界點處做長軸的垂線，則垂線被多邊形區(qū)域截得若干條線段。由于多邊形的輪廓可能不規(guī)則，因此截得的線段可能存在多條，在這種情況下，將該垂線上被截得2 個最遠的交點相連，獲得的線段作為最終線段。最后，比較截得的200 條線段，選取其中最長的一條，作為區(qū)域的短軸。

6)多時態(tài)數(shù)據(jù)顯示。多時態(tài)數(shù)據(jù)反映了數(shù)據(jù)在一段時間序列中的變化情況。由于受到海洋中海浪、海風、洋流等影響，泄漏原油的擴散區(qū)域總是在不斷移動和變化，因此，掌握該區(qū)域的移動趨勢和變化軌跡，對于采取進一步的處理和應對措施至關重要。為了實現(xiàn)以上目標，需要對多個時間點下的區(qū)域變化軌跡進行顯示。由于對于原油點跡和泄漏區(qū)域的數(shù)據(jù)，提取自RS 圖像，因而需要采集不同時間點的RS 圖像，并對數(shù)據(jù)進行采集和分析，從而得出一段時間序列下的數(shù)據(jù)變化情況。

7)數(shù)據(jù)可視化。由于DotSpatial 不能夠支持多時態(tài)數(shù)據(jù)顯示，因而本文給出了一種結合色彩填充和圖形填充的方法，來表示一段時間內(nèi)原油點跡的變化軌跡。同時，本文還給出了一種判斷數(shù)據(jù)可信性的方法，采用4 個級別(高、中、低、不可信)來判斷數(shù)據(jù)的可行性。最終，本文在數(shù)字地圖上，使用經(jīng)緯度、海拔、等高線和不同的點線表示出以上數(shù)據(jù)，如圖4 所示。

圖4 多時態(tài)數(shù)據(jù)顯示結果Fig.4 The result of multi-temporal data rendering

8)匹配地圖。DotSpatial 中提供了一系列的控件和API，能夠滿足GIS 系統(tǒng)的地圖繪制和匹配需求。例如，其能夠自動繪制出經(jīng)緯線，同時能夠使用鼠標繪制出多邊形等。

9)輸出數(shù)據(jù)。通過步驟1 ～步驟8，能夠在開源的GIS 系統(tǒng)DotSpatial 中，顯示海上原油泄漏區(qū)域的動態(tài)信息，并能夠采用不同的圖示和色彩表示更加豐富的信息。在Windows GDI +的幫助下，通過鼠標拖拽選中一部分區(qū)域，則DotSpatial 可以自動地將該區(qū)域中的相關信息以HTML 格式進行發(fā)布，通過構建合適的后端Web 發(fā)布系統(tǒng)，則可以發(fā)布經(jīng)過DotSpatial 處理的直觀地理信息。

3 結 語

隨著人類海上活動的日益頻繁，越來越需要一種海上信息分析和發(fā)布系統(tǒng)。本文選取海上原油泄漏事件為應用場景，討論了在該應用中，對于信息監(jiān)視、收集、分析和發(fā)布的各種需求，并使用開源GIS 平臺DotSpatial 構建了海上原油泄漏信息分析與發(fā)布系統(tǒng)，并對該系統(tǒng)的關鍵技術，如地理信息提取、編譯、多時態(tài)顯示、輸出等進行研究，成功開發(fā)了一種開源的信息分析與發(fā)布平臺。同時本文的研究成果能夠應用于廣泛的應用領域，具有較強的通用性。

［1］LI Yin，LAN Guo-xin，LIU Bin-xin，et al．Dynamic analysis on oil spill in Mexico Bay based on remote sensing and GIS［J］．Marine Environmental Science，2012(1):80 －82．

［2］JIAO Jun-chao．Research on the prediction of oil spill based on GIS In Bohai Bay［J］．Marine Environmental Science，2011(5):96 －101．

［3］胡安康．船舶數(shù)字化設計的主脈——CAE 信息鏈［J］．艦船科學技術，2008，30(1):66 －70．HU An-kang．The main vein of digital ship design——CAE information chain［J］．Ship Science and Technology，2008，30(1):66 －70．

［4］LIU Chang-ming，CHEN Luo．GDAL multi-source spatial data access middleware［J］．Geo-spatial Information，2011(5):58 －62．

［5］任席闖，左麗芬，李樹明，等．基于GIS 的超低頻通信服務區(qū)的預測［J］．艦船科學技術，2013，35(4):117 －119．REN Xi-chuang，ZUO Li-fen，LI Shu-ming，et al．The prediction of VLF communication service zone based on GIS［J］．Ship Science and Technology，2013，35(4):117 －119．