基于GIS的河道可視化及斷面提取方法研究

耿　敬，鄭天駒，李明偉，張清明

(1.哈爾濱工程大學船舶工程學院，哈爾濱 150001；2.黑龍江省松花江航運建設管理中心，哈爾濱 150001)

?

科技成果

基于GIS的河道可視化及斷面提取方法研究

耿敬1，鄭天駒1，李明偉1，張清明2

(1.哈爾濱工程大學船舶工程學院，哈爾濱 150001；2.黑龍江省松花江航運建設管理中心，哈爾濱 150001)

摘要：地理信息系統(GIS)是地理學與計算機信息科學的有機結合，具有強大的可視化功能，在空間分析方面優勢尤為顯著。文章在對以往研究學習的基礎上，對利用GIS和CAD軟件進行河道三維可視化及斷面提取的技術流程進行了系統的分析，并深入研究了流程的關鍵技術，同時，對流程在操作過程中可能遇到的問題及解決辦法進行了討論。以松花江依蘭段河道地形測量數據為例，利用GIS和CAD軟件建立了河道三維可視化數字模型，并提取河道斷面信息，驗證了技術流程的可行性。

關鍵詞：地理信息系統；河道可視化；斷面提取；DEM數據

0引言

河道斷面數據是研究河道演變、洪水淹沒及沖淤變化的基礎，是航道整治以及沿河流域規劃等工程設計必不可少的底層數據。傳統方法利用CAD格式的數據進行河道斷面信息提取存在精度低、連續性差且工作量大的不足。地理信息系統(GIS)是在地理學、地圖學、測量學和計算機科學等學科基礎上發展建立起來的一門學科，具有空間數據的獲取、存儲、顯示、編輯、處理、分析、輸出和應用等功能，自身具有一定結構，是一個完整的系統[1]。GIS擁有強大的空間分析能力，能夠將河道的幾何形態數據(包括CAD數據、遙感數據等)處理整合，建立三維可視化數字模型，由一系列地面點(x，y)位置及其相聯系的高程Z組成[2]，利用可視化數字河道模型可以直觀、便捷的提取河道指定位置的斷面數據信息，同時為跨河大橋選址、航電樞紐建設、河道的綜合治理、規劃和開發提供決策依據。

目前，地理信息化系統在河道可視化及斷面提取方面應用的研究已取得一定成績。劉暢等人提出一種在ArcGIS平臺上將二維散點地形圖轉換成DEM(數字高程模型)并提取河道斷面地形數據的方法[3]。王偉等人利用GIS建立河道可視化模型，對河道槽蓄量、沖淤量、沖淤厚度進行計算，使河道演變分析由傳統手工或簡單的計算統計發展到計算機系統性處理與智能分析，實現了質的飛躍[4]。張獻志等人基于GIS組件的內置工具和模塊，采用將二維斷面線進行TIN的三維表面內插法實現河道任意斷面數據的半自動提取，并探討了轉換過程中存在的問題及相應的解決措施[5]。王鐵軍等人分析了CAD到GIS的數據轉換過程中容易出現的問題，并總結了為克服這些問題而在數據生產、轉換的各個步驟上應采取的措施，提出CAD與GIS生產建設中目前存在的問題及解決的思路[6]。賀巧寧基于ArcGIS Engine通過二次開發實現了河道斷面生成與數據獲取，并對DEM、TIN與等高線提取的河道斷面高程數據的差異進行了討論，解決了傳統河床演變斷面分析必須立足固定斷面的局限[7]。楊婷利用AutoCAD.NET API二次開發，提出一套快速提取高程信息的技術及其實現流程[8]。

文章研究了基于GIS和CAD軟件對河道可視化模型建立與斷面提取的方法與技術流程，并以松花江依蘭河段為研究區，對該地區數字地形生成過程中參數的選取及可能出現的相關問題進行了討論，實現了河道的可視化以及任意斷面信息快速提取，為河道水動力數值模擬計算及工程可視化管理平臺建設提供基礎數據。

1研究區概況

松花江流域位于我國東北地區，地理位置為E119°52′～132°31′、N41°42′～51°38′。松花江為黑龍江右岸的一大支流，流域東西長920 km，南北寬1070 km，流域面積56.59×104km2，占整個黑龍江流域面積的30%，占我國境內黑龍江流域面積的63%。行政區劃屬黑龍江省、吉林省、內蒙古自治區和遼寧省，其中內蒙古自治區占27.4%，吉林省占24.1%，遼寧省占0.1%，其余均屬黑龍江省，占全流域面積的48.4%，主要為嫩江左岸及松花江干流部分[9]。

松花江干流哈爾濱以上為平原區，坡度較緩，河道比降為0.05‰，灘地較寬闊，有濕地和牛軛湖，河道主槽寬一般為400m～600m，平槽時平均水深4m～7m；哈爾濱至佳木斯段為高平原和丘陵山區，河谷狹窄；在牡丹江口上游為三姓淺灘，長約41 km，該處河寬2000m左右，島嶼眾多、水流分散、河床比降較陡，河床底質為堅硬的板巖、砂巖和花崗巖，水深較其它河段為淺，一般只有1.2m～1.3m，枯水期僅有0.8m左右；佳木斯以下為低平原區，屬三江平原一部分，河道寬闊，主槽寬達1500m～2000m，平槽平均水深2m～3m，河道中有沙洲，河槽滯蓄能力較大[10]。

研究區域位于依蘭縣境內，松花江干流牡丹江匯合口上游1000m處，下距依蘭縣城2.5 km，上距哈爾濱市330 km(航道距離)。上游約81.5 km處為通河水文站，通河水文站至依蘭壩址之間無大支流匯入，下游2.8 km處為依蘭水文站。地理位置見圖1。

圖1　松花江依蘭段遙感影像圖

2數據預處理與研究方法

2.1數據的預處理

河道的數據通常以紙質地形圖或CAD地形圖保存。紙質地形圖需要進行空間校正、矢量化處理，ArcGIS自帶的ArcScan可直接對紙質地圖的掃描文件進行矢量化處理，但對地形圖本身的標注規范及圖片質量要求較高，很多地形圖無法處理或處理后出現大量錯誤，因此現在仍多采用人工矢量化。考慮到目前大多數河道實測信息多以CAD電子文件格式存儲，本文僅以松花江依蘭段河道CAD數據實現河道可視化過程為例，對紙質地形圖的矢量化不再進行討論。該段河道10年內未出現顯著變化，故使用該段河道2006年1:2000地形圖數據進行研究。

數據處理的目的是從原始數據中提取出可用于在ArcGIS中建立河道數字高程模型及可視化所需的數據信息。處理過程中，若一個河道的信息保存在多份圖紙中，需首先在CAD中對地圖形進行拼接處理，隨后刪除無關圖層、填充符號等冗余數據。完成后使用ArcGIS、南方CASS9.0提取等高線與高程點數據。保存數據時選用高斯克呂格投影坐標系中的Beijing 54坐標系，轉換后的高程點和等高線的屬性表中必須得包含高程信息字段且有具體值，否則需結合注記，對高程信息進行補充。下面對等高線與高程點的提取進行介紹。

2.1.1等高線的處理

經過初步處理的CAD數據仍可能存在問題，對于容易混淆的各種線狀符號，如航標符號、與等高線相交的線等應將之刪除。等高線指的是地形圖上高程相等的相鄰各點所連成的線，應該為閉合曲線[11]。在本研究中，等高線可能由于標注信息、丁壩存在的原因而間斷。對于因標注信息而間斷的等高線，參照臨近的等高線形狀將其重新連接；對于因研究區丁壩間斷的等高線，因為在1∶2 000比例尺的地形圖中丁壩的外形特征已經可以較為明顯的表達，同時丁壩的位置與高程對航電樞紐工程有著重要的影響，故予以保留，將間斷的等高線與丁壩相連。完成以上處理后，參照地形圖對ArcMap中的等高線數據進行推算并賦予高程值。

當遇到地形圖的部分等高線高程值無法確定時，可依據高程點數據進行參考。隨后將等高線屬性表中多余字段刪除，只保留高程字段，將長度過小與無意義的線條刪除。最后檢查每個線要素，將同一條等高線連接起來，確保每條等高線都是一條連續的一條閉合線。至此，獲得帶有高程信息的等高線數據。

2.1.2高程點的處理

為了獲得高精度的數字地形，在生成數字地形時添加高程點作為對等高線的補充。大多數CAD地形圖并未給高程點賦值，而是采用在點的旁邊設置文本標注的方式，為了能夠在ArcGIS中建立三維地形，需進行賦值操作。將注記的內容賦給高程點的方法有很多，本文選擇使用地形圖專業處理軟件南方CASS9.0生成的高程點數據，將高程值保存在Z字段中，而非通過ArcToolBox中的數據轉換功能得到高程點數據，主要原因是實踐發現在CASS下可以更好、更快、更準確的實現將注記文字的高程信息賦給高程點。將賦值后的高程點數據添加到ArcMap中，對高程信息的準確性進行檢查，并將生成的高程點數據導出為.shp格式，從而獲得帶有高程信息的點數據。

2.2研究方法

文章以河道CAD地形圖為原始數據，提出的河道可視化及斷面提取方法主要研究了河道數字高程模型生成及斷面提取實現等關鍵問題，給出了包含數據預處理、河道地形圖等高線與高程點的提取、TIN數據的生成、DEM數據的生成、河道可視化及斷面提取的方法與技術流程。技術路線見圖2。

3數字地形數據生成

3.1河道TIN數據生成

TIN是由不規則三角網構成的一種矢量數據，是三維空間中的分塊線性模型，在整個區域中是連續的，但不可導[12]。我們可以根據這些三角網來顯示三維立體地形。TIN相當于是對采樣點的一個平面三角剖分，因而它是一個連通平面圖。TIN模型的優點是能以不同層次的分辨率來描述地形表面，可有效避免平坦處的數據冗余。利用提取獲得的等高線數據和高程點數據，生成河道TIN數據，等高線的高程使用高程字段值，高程點的高程使用Z字段值。有時候可能由于等高線或高程點賦值錯誤導致生成的TIN面異常，這時候可以在ArcScene中加載生成TIN面，調整夸張系數為20或更大，以方便觀察，加載等高線數據，并在Base Height中選擇use a constant value or expression，將高程顯示設置為elevation，使兩者擬合，通過疊加觀察找出影響精度的等高線或高程點并刪除。獲得的TIN數據見圖3。

圖2技術路線圖

圖3河道TIN地形數據

3.2河道DEM數據生成

TIN數據能夠較好地表示復雜地形，以該格式存儲的三維地形信息具有較好的顯示效果。但由于其存在數據結構復雜，不便于規范化管理，難以與矢量和柵格數據進行聯合分析，在后續的應用中存在諸多不便，需將TIN經3D Analyst模塊，轉換生成河道DEM柵格數據，CELLSIZE是柵格像元的大小采樣距離，直接影響轉換得到DEM的分辨率及文件大小[13]，經過試驗，在本研究中大小設置為5，既滿足了地形數據的精度要求，又保證了加載速度。隨后對獲得的河道DEM柵格數據進行垂直高程夸張、高程重新配色、渲染與光照調整等操作，完成可視化河道生成。

4河道斷面信息提取

4.1斷面信息提取

4.1.1提取指定位置斷面

3D Analyst工具可提取河道指定位置的斷面信息。斷面可沿在橫穿河道的任意線要素生成，在ArcMap中，使用3D Analyst模塊功能創建目標位置斷面圖的表面并進行線插值，在河道數據表面點擊創建斷面多段線，完成后創建斷面。在屬性欄中，可更改斷面圖的布局。

4.1.2特征點生成斷面

按照特定位置高程點生成斷面圖需要提取沿線特征點或沿某一走向特定位置點的高程。當提取的特征點為線要素的折點時，直接調用Features-Feature Vertices to Points命令；當提取的特征點需滿足某種特定距離要求時，如：沿水流動力軸線，每隔5 m提取1個特征點，可先創建一個點圖層，在開啟編輯的狀態下，選中生成點的參考線，然后選擇Construct Points設置相關參數，即可沿線提取滿足某種特定條件的點要素。將提取的特征點和DEM數據加載到ArcMap中，從DEM中提取高程值賦予特征點，提取出的高程存儲在Raster Value字段中。生成河道斷面圖使用ArcGIS創建圖表功能中的創建斷面線功能，使DEM斷面線的橫軸為距離第1個點之間的距離，縱軸為高程值。

以上2種斷面提取方式可基本滿足河道分析工作的所有斷面數據提取需要。將河道DEM數據與遙感影像進行地理配準后，即可根據需求提取任意坐標處的河道斷面信息，斷面數據不再受到實測固定斷面資料的局限。在后續研究中亦可將多年的數據全部進行數字化處理，獲得多組不同時間的河道斷面數據。

4.2河道斷面拓撲關系圖

拓撲是描述地理空間關系的模型，一種維護地理空間實體間空間幾何關系的機制。河道斷面拓撲關系主要用來描述河道斷面的上下游關系，擬在通過河道斷面對象的空間位置，結合河道中心線或深槽線，將河道斷面對象之間的上下游關系建立起來(圖4)。并利用獲得的斷面及其拓補關系，建立河道數值模型。

5結論

河道斷面信息用途廣泛，傳統的獲取方法多存在耗時長、不準確、存在重復性工作等不足。本文深入研究了基于GIS的實現河道可視化及斷面地形數據提取的方法，通過數據處理手段，提取CAD地形圖中的等高線與高程點并賦值，將二維數據轉化為GIS中的三維數據，并聯合兩者生成河道可視化數字高程模型，信息更為可靠。以松花江依蘭段河道的CAD地形圖原始數據為例，使用該方法完成了地形信息數字化與河道三維地形數據生成，利用生成的地形提取了河道斷面數據，將獲得的河道斷面數據以上下游為序建立聯系，為壩址料場位置的選擇及工程可視化管理平臺建設提供基礎數據。該方法通過可視化手段直觀地表達河道地形信息并對所需信息進行提取，可作為洪水演進模擬、河道演變預測分析、航道整治沖淤分析等提供技術支撐。

圖4　河道斷面拓撲關系圖

參考文獻：

[1]王向東，湯立群，丁賢榮.應用GIS計算分析長江鎮揚河段的河床演變[J].泥沙研究，2008(06): 68-73.

[2]胡紅兵，胡光道.基于DEM的河道演變可視化表達和定量分析[J].長江科學院院報，2008(02): 26-28.

[3]劉暢，辛小康，尹煒.基于DEM的河道斷面數據提取方法[J].人民長江，2015(S1): 9-11.

[4]王偉，白亮，李圣偉.基于GIS的河道演變與可視化分析方法研究[J].人民長江，2006(12): 4-7.

[5]張獻志，勇施，陳煉鋼.基于ArcGIS的河道斷面數據提取方法[J].水電能源科學，2012，30(4):139-141.

[6]王鐵軍，麻璽琳.從CAD到GIS的數據轉換[J].東北測繪，1999，22(3): 12-13.

[7]賀巧寧.基于GIS的河床演變斷面分析方法體系研究[D].南京：河海大學，2007.

[8]楊婷.CAD 地形圖高程信息快速提取的技術與實現[J].地理空間信息，2015(01)：162-164.

[9]袁杰，李弘義，李勝軍.松花江流域徑流變化特征及趨勢分析[J].海河水利，2013(05)：33-36.

[10]高爽，劉文斌.基于洪水預報模型的哈爾濱站預報[J].黑龍江水利科技，2009(04)：10-11.

[11]張大偉，權錦，何曉燕，等.基于地理信息系統的一維潰壩洪水分析系統研發及應用[J].水利學報，2013(12)：1475-1487.

[12]李曉印，郭達志，張佃國.DEM與TIN的數據精度與應用領域的對比分析[J].濟南大學學報(自然科學版)，2009(01): 76-79.

[13]程軍蕊，王侃，馮秀麗，等.基于GIS的流域水質模擬及可視化應用研究[J].水利學報，2014(11)：1352-1360.

Research of River Course Visualization and Section Extraction Method based on GIS

GENG Jing1;ZHENG Tian-ju1;LI Ming-wei1and ZHANG Qing-ming2

(1.Ship Engineering Institute Harbin Engineering University，Harbin 150001，China；2.Heilongjiang Provincial Songhua River Shipping Construction Management Center，Harbin 150001，China)

Abstract：Geographic information system is the organic combination of geography and computer information science，having a powerful visual function，particularly advantageous in spatial analysis.Based on the learning to previous studies，GIS and CAD software were used to analyze systematically the technical process and section extraction of three-dimensional visualization of river course，and the key technique of process was researched deeply，meanwhile，the problems that may be met possibly during the process of operation and solved method were discussed in the paper.Taking the topographic survey data of river course at Yilan river section along the Songhua River as the example，GIS and CAD software were used to set up the three-dimensional visualization digital model and extract the section information of river course，verifying the feasibility of technical process.

Key words：geographic information system；visualization of river course；section extraction；DEM data

中圖分類號：TV147

文獻標識碼：A

[作者簡介]耿敬(1968-)，女，江蘇宜興人，教授，碩士生導師，研究方向為河海及水利工程設計理論與建造技術、智能航道技術；鄭天駒(1990-)，男，河南鄭州人，碩士研究生；李明偉(1984-)，男，內蒙古通遼人，博士后，研究方向為河海與水利工程；張清明(1965-)，男，黑龍江哈爾濱人，研究員級高工，主要從事航電樞紐建設管理工作。

[基金項目]國家自然科學基金(51509056);交通運輸部信息化技術研究項目(2014364554050);中國博士后基金面上資助項目(2015M571394);黑龍江省博士后資助經費(LBH-Z14059);中央高校基本科研業務費專向資金資助(HEUCF150108)。

[收稿日期]2016-01-16

文章編號：1007-7596(2016)03-0001-05