基于HYDRO1K、SRTM3和ASTER GDEM的韓江流域水文地形信息對比

劉 遠，周買春

(華南農業大學水利與土木工程學院，廣州 510642)

數字高程模型(Digital Elevation Model, DEM)是地表曲面或粗糙度的數字化表達，是計算機地形分析基礎數據。像過去地形圖一樣，世界上許多國家都把DEM數據作為一種國家基礎設施來建設，就像道路交通、通訊設施一樣，通過政府立項，由具有技術力量的科技機構競標，按一定標準生成，委派專門機構管理和負責更新，并鼓勵廣泛地推廣應用。隨著DEM數據的豐富和社會需求迫切，一種趨勢是開放DEM數據供公眾免費下載，如幾種精度不一、覆蓋全球可免費下載的DEM有HYDRO1K[1]、SRTM3[2]和ASTER GDEM[3]。不同的DEM數據產品由不同的機構發布，其數據基礎和生成方式的不同隱藏了數據質量的差異。

DEM數據在流域水文研究中有廣泛的應用。在流域地表水、地下水運動中，地形代表重力對水的作用，使得流域徑流隨地形起伏從高處流向低處，許多分布式水文模型基于DEM數據，通過提取的數字河網及相關地形信息，將流域產流匯出流域。除流域匯流外，TOPMODEL模型還通過流域土壤—地形指數的分布來模擬流域土壤蓄滿狀況及徑流形成區域，DEM數據為土壤—地形指數計算提供了極大的方便和效率。此外，在利用站網氣象數據(如降雨、大氣壓強、氣溫等)通過插值形成流域空間分布數據時，也常用DEM代表的地面高程對插值結果進行修正，如在降雪融雪模型、土壤凍融模型和冰川水文模擬中，近地氣溫需考慮隨高程的直減率。DEM數據質量對流域水文研究成果的可靠性和精度有極大的影響，如筆者根據HYDRO1K、SRTM3和ASTER GDEM 3種DEM數據提取的韓江流域數字河網則存在明顯差異[4]。

在DEM易于獲得的情況下，有關DEM提取地形信息的研究主要集中在信息提取方法[5-9]和水平分辨率對提取信息精度的影響[10-14]2大方面。然而，地形信息的可靠性和精度歸根結底取決于DEM是否精確地反映了地面高程，基于不同DEM產品提取流域水文地形信息進行比較研究，國內外還鮮有報道。為此，本文繼續深入比較HYDRO1K、SRTM3(version 4)和ASTER GDEM 提取的流域水文地形信息，分析不同DEM數據源對流域地形特征描述的差異，確定最真實反映流域水文地形狀況的DEM數據，為流域水文研究選擇恰當的DEM數據提供依據。

1 基于DEM提取流域水文地形信息的方法

流域水文地形信息采用分布式流域水文模型BTOPMC[15-17]的地形子模型來提取。BTOPMC的地形子模型由DEM填洼運算、生成數字河網、劃分子流域及其他提取地形相關信息(如水流方向、坡度、地形指數等)的算法等構成。

(1)填洼。DEM實際上是通過某種數學插值方法得到的柵格高程，是連續地面的一種離散表達，生成DEM的高程數據源密度和數學插值方法都會造成DEM中或多或少存在某些虛假洼地，即一些柵格上的水流流不出流域出口(對于一個閉合的外流流域，不存在真實洼地)，這些洼地可認為是DEM數據的一種誤差。要提取一個完整的流域數字河網，首先必須對原DEM進行洼地處理，即通過修改洼地DEM柵格的高程值，使它高于與之相鄰的8個柵格其中之一，使得洼地柵格的水流能夠順暢地流出。由于一個(或一群)洼地柵格可能套合在一個更大的洼地中，填洼通常是一個反復的過程，使全部的洼地得以消除，它是地形子模型中時間消耗最長的步驟。BTOPMC采用逐步填洼，直到全部洼地填平，循環終止，每次填洼采用隨地貌傾斜度變化的小高程增量：

(1)

式中：i、j分別是洼地的行坐標和列坐標；hc是用戶定義的一個小高程增量閾值(避免造成新的洼地)，取hc=0.1 m(考慮到HYDRO1K、SRTM3和ASTER GDEM的高程精度是米)；α、β是x、y向地形傾角的權重系數，考慮韓江流域的實際情況取α=β=1；Nr、Nc是填洼前DEM的行、列數；填洼后，洼地高程h(i,j)=hmin(i,j)+dh(i,j)；hmin(i,j)是與洼地柵格(i,j)相鄰的最低柵格的高程。

(3)上游集水面積和流域邊界。上游集水面積是指水流流入當前柵格的所有單元格面積(包括當前柵格)。計算任一個單元格的上游集水面積，只需沿著水流方向反向追蹤，直至到達流域邊界(上游集水面積為1的柵格)，當整個水流方向矩陣搜索完畢后，位于追蹤線路上的全部柵格數目即為上游集水面積(以柵格的個數表示)。確定流域邊界首先要確定流域的出口，即上游集水面積最大的柵格。水流最終流入到流域出口柵格的所有單元格都處于該流域范圍，這些柵格所占區域的邊界即為該流域邊界。

2 研究區域和數據來源

2.1 研究區域

韓江流域是廣東除珠江以外的第2大流域，流域面積30 112 km2，位于粵東、閩西南，經緯度為115°13′～117°09′E、23°17′～26°05′N(見圖1)。韓江的主源是梅江，發源于紫金縣烏突山七星崠，在大埔縣三河壩與發源于寧化縣武夷山南段木馬山北坡的汀江匯合(三河壩以下始稱韓江)，最后經韓江三角洲，分北、東、西溪在汕頭市出南海，全長470 km。韓江流域以山地為主，約占流域總面積的70%。此外，丘陵和平原大約為25%和5%，丘陵主要分別分布在梅江流域和其他干支流谷地，平原主要分布在韓江三角洲。

圖1 韓江流域地理位置和水系Fig.1 The geographic location and drainage network of Hanjiang River basin

2.2 DEM數據來源及填洼處理

2.2.1 3種常用DEM數據的介紹

HYDRO1K[1]是美國地質勘測局地球資源觀測系統數據中心(US Geological Survey’s EROS Data Center)與聯合國環境計劃全球資源信息中心(United Nations Environmental Program/Global Resources Information Database)1996年發布的1 km分辨率的全球DEM(不包括格陵蘭島和極地)。HYDRO1K是在另一個基于多源數據的全球DEM GTOPO30[18](分辨率30″，約1 km)的基礎上，結合DCW(Digital Chart of the World)的流域河網和邊界，增強了水文特性。

SRTM[2](Shuttle Radar Topography Mission)是在由美國國家航空航天局(National Aeronautics and Space Administration，NASA)和國家空間信息情報局(National Geospatial-Intelligence Agency，NGA)合作完成。2000年2月，通過裝載干涉雷達于“奮進”號航天飛機在太空飛行11 d，獲得地球北緯60°至南緯56°之間的雷達影像數據，覆蓋地球80%以上的陸地表面。SRTM于2003年6月發布，按精度可以分為SRTM1和SRTM3，分辨率分別為1″(約30 m)和3″(約90 m)，前者僅有北美地區數據，后者覆蓋全球。經過多版的改進，SRTM已更新至第4版。

ASTER GDEM[3](Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model)是由美國國家航空航天局和日本經濟產業省于2009年6月發布的30m分辨率的全球DEM。ASTER GDEM覆蓋范圍占全球陸地面積的99%(北緯83°至南緯83°)，是目前覆蓋最廣的地形數據。該數據是根據NASA的新一代對地觀測衛星TERRA的詳盡觀測結果制作而成。從2000年發射至今，ASTER還在不斷更新數據，算法也不斷改進，目前已發布了第2版ASTER GDEM數據。

2.2.2 DEM數據處理

用韓江流域所在的矩形邊界分別裁剪3種DEM數據，見圖2。3種DEM反映的地形變化趨勢基本一致，流域地勢高低分明，高程較低的山谷、河流清晰可見；具有高分辨率的ASTER GDEM和SRTM3更能反映地形的精細度。統計它們的特征值，結果見表1。同時，統計各個高程值的柵格出現的頻率，得到累積分布曲線(見圖3)。表1顯示，HYDRO1K、SRTM3和ASTER GDEM在矩形邊界內的平均值很接近，但DEM分辨率越低，最大高程越小，最小高程越大，即DEM的變化范圍越小，反映地形起伏的柵格高程的標準差也越小。由圖3可以看到，3條高程累積分布曲線很接近，表明3種DEM反映的研究區域內的地貌變化是一致的，特別是SRTM3和ASTER GDEM；HYDRO1K中，大高程柵格的比例較SRTM3和ASTER GDEM小，說明較低的分辨率使得DEM反映的地形趨于平坦。

圖2 研究區域3種DEM的空間分布Fig.2 Spatial distributions of 3 DEMs in the study region

圖3 研究區域的高程累積分布曲線Fig.3 Elevation frequency curve of the study region

表1 研究區域3種DEM的基本特征 m

利用BTOPMC地形子模型對裁剪的研究區域內3種DEM的原始數據進行填洼處理，生成各自的無洼地DEM，用于計算流域的水文地形信息。

3 結果分析

3.1 水流方向對比

在填洼后的DEM基礎上，采用D8法確定研究區域各柵格的水流方向[見圖4(d)]，并裁剪韓江流域邊界(潮安水文站以上)內的部分，所用的流域邊界是由各自DEM提取得到的，結果見圖4(a)～(c)。3種DEM提取的水流方向柵格在空間上分布基本一致，分辨率較高的SRTM3和ASTER GDEM吻合程度相對較高；HYDRO1K在汀江中、上游平原地區和ASTER GDEM在梅江的中游出現了較大面積的同流向區域，在這些地方提取的河網出現與實際不相符的平行河道或錯誤河道[4]。這與采用的D8法有關，根據D8法對流向的定義，其計算得到的流向往往會在 45°或其倍數的方向上產生平行徑流。

圖4 3種DEM提取的韓江流域水流方向的空間分布Fig.4 Spatial distributions of drainage directions in Hanjiang river basin extracted from 3 DEMs

由圖5的對比可以看到，3組數據8個方向的比例基本一致，SRTM3和ASTER GDEM更為接近；以水平和垂直方向居多，HYDRO1K最為明顯，它在該4個方向的比例達到66%，又以北向最少；ASTER GDEM雖也以水平和垂直方向居多，但8個方向的比例較其他2組數據均勻。這同樣與采用的D8法有關，因為沿水平或垂直方向的路徑比沿對角方向短，在相鄰8個柵格的高程相差不大的情況下，根據最陡坡度法則確定的流向多為水平或垂直方向；另一方面，DEM分辨率越低，柵格尺寸越大，沿水平或垂直方向的路徑與沿對角方向的路徑就相差越大，但是低分辨率DEM使地勢趨于平坦，相鄰柵格間的高程變化變小，所以DEM分辨率越低就越容易產生水平或垂直的流向。

圖5 3種DEM提取的韓江流域水流方向對比Fig.5 Comparisons of drainage directions in Hanjiang river basin extracted from 3 DEMs

3.2 坡度對比

在填洼后的DEM基礎上，沿水流方向計算流域內各柵格的坡度，結果見圖6。3種DEM提取的坡度柵格在空間上分布基本一致，高山、平地清晰能辨；分辨率較高的SRTM3和ASTER GDEM的吻合程度相對較高，反映的地面坡度變化更有層次，地形更加精微；分辨率較低的HYDRO1K反映的地形較為

單調，呈現大面積的平坦區域，與實際地形有所不符。3組數據的統計特征見表2，HYDRO1K的平均值和標準差都比SRTM3和ASTER GDEM小很多，說明低分辨率DEM使得韓江流域的地形趨于平坦、均勻。

圖7是根據3組坡度數據統計得到的各級坡度柵格出現頻率的累積分布曲線。HYDRO1K的曲線明顯偏離SRTM3和ASTER GDEM的曲線，其坡角大于3°的柵格百分比還不到50%，嚴重坦化了流域的地形，將影響由坡度計算的其他地形因子(如地形指數)的精度，同時也影響由HYDRO1K提取的數字河網的精度，其數字河網在坡度較小的成片區域中生成了與實際不相符的平行河道[4]。SRTM3和ASTER GDEM的坡度累積頻率曲線吻合較好，說明它們反映的韓江流域地勢具有較好一致性；ASTER GDEM的曲線呈躍階分布，說明ASTER GDEM雖具有較高的水平分辨率，但垂直精度不如SRTM3(前者垂直精度±20 m，后者±16 m)，導致坡度的躍階。

圖6 3種DEM提取的韓江流域坡度的空間分布Fig.6 Spatial distributions of surface slope in Hanjiang river basin extracted from 3 DEMs

圖7 3種DEM提取的韓江流域坡度的累積分布曲線Fig.7 Surface slope distribution curve of Hanjiang river basin extracted from 3 DEMs

3.3 地形指數對比

采用單流向算法計算韓江流域的地形指數，結果見圖8。3種DEM提取的地形指數柵格在空間上分布基本一致，地形指數較大區域，即土壤缺水容量小或土壤達到了飽和狀態，能呈現出流域的河流、溝谷形狀；分辨率較高的SRTM3和ASTER GDEM吻合程度相對較高，其地形指數的空間變化更有層次，即更能反映流域土壤水分空間變化，符合山地、丘陵流域的實際情況，這與高分辨率DEM更能反映流域地形變化的結果是一致的。3組數據的統計特征見表2，由于HYDRO1K分辨率低使得地形區域均坦化，流域的坡度變小，所以其地形指數在全流域的平均值和標準差都明顯較其他2者大，據此進行流域水文模擬時，流域土壤水分易趨于飽和。

圖9是根據3組地形指數數據統計得到的各級地形指數柵格出現頻率的累積分布曲線，HYDRO1K的分布曲線與另外2條曲線相差較大，其地形指數大小分布相對均勻，而SRTM3和ASTER GDEM的地形指數集中分布在10以下，分別達到了82.3%和75.0%，說明流域大部分區域的坡度較大，土壤水分不易于飽和，與流域的山地、丘陵地形相一致；SRTM3與ASTER GDEM的分布曲線吻合較好，在地形指數較小處，幾乎不存在差別，雖然ASTER GDEM的坡度分布曲線呈明顯的階梯分布，但其地形指數已幾乎不存在這一現象，說明由DEM垂直精度引起的誤差，在地形因子的計算過程中得到一定程度的衰減。

圖8 3種DEM提取的韓江流域地形指數的空間分布Fig.8 Spatial distributions of topographic index in Hanjiang river basin extracted from 3 DEMs

DEM數據坡度平均值最大值最小值標準差地形指數平均值最大值最小值標準差HYDRO1K0．0720．6108．600×10-80．07312．60433．4097．4025．115SRTM30．2422．5692．400×10-80．1838．07037．4883．5473．639ASTERGDEM0．2384．3841．330×10-70．1908．66937．0072．6163．750

圖9 3種DEM提取的韓江流域地形指數的累積分布曲線Fig.9 Topographic index distribution curve of Hanjiang river basin extracted from 3 DEMs

4 討論與結論

筆者[17]對HYDRO1K、SRTM3和ASTER GDEM提取的韓江流域數字河網的對比分析表明，DEM的垂直精度對提取的河網精度起控制作用，所以ASTER GDEM提取的河網精度不如SRTM3；低分辨率的HYDRO1K在大尺度流域上能提取到一定精度的河網，但在流域地勢平坦區域的河網精度較低。進一步對比分析這3種DEM數據提取的地形信息，得到與之相輔相成的結論：①在大尺度的韓江流域上，低分辨率的HYDRO1K提取的地形信息在一定程度上能反映流域的地形變化，但坦化了流域地形，使得基于地形指數的流域土壤水分分布趨于均勻，在產流表現上易趨于飽和，從水文意義上達不到刻畫像韓江這一尺度流域的山地、丘陵地形。②DEM的水平分辨率和垂直精度決定DEM反映地表變化的精細度，DEM水平分辨率和垂直精度越高，反映的地形越精微。SRTM3和ASTER GDEM提取的韓江流域地形信息能精微地反映流域的地形變化和土壤水分空間格局變化。③當水平分辨率達到一定精度時，與水平分辨率相比，DEM的垂直精度對提取地形信息的真實性起控制作用。ASTER GDEM的垂直精度 不如SRTM3，所以提取的地形信息不如SRTM3精確。④地形信息提取精度與DEM的填洼算法和流向算法有關。BTOPMC采用隨地貌傾斜度變化的小高程增量逐步填洼，填洼過程對區域地形信息的擾動較小；在采用D8單流向算法確定柵格水流方向時，沿最陡坡度流向比沿最大落差流向更容易產生水平和垂直流向，這在低分辨率的HYDRO1K中尤為明顯；但由于HYDRO1KDEM生成時融合了DCW河網信息，2種流向算法對數字河網和地形指數的影響只限于那些不具有DCW河網信息的源頭區域。

不同DEM數據源提取的流域地形信息存在差別，其反映的流域地形變化、土壤水分空間格局也不同，將對分布式水文模型的產、匯流模擬產生影響。所以，流域水文模擬應當選擇最能反映流域真實情況的DEM數據。根據HYDRO1K、SRTM3和ASTER GDEM 3種常用DEM提取的流域地形信息和數字河網的對比分析，就韓江流域而言，選擇SRTM3進行水文模擬是最合適的。雖HYDRO1K的精度較低，但由于融合了DCW河網信息使得水文特性得到加強，所以提取的韓江流域地形信息和數字河網仍具有相當的精度，可用于流域的水文模擬。可以設想，若在SRTM3中融入河網信息(DCW，或較之精度更高的河網)，必將進一步提高其水文應用的精度。

□

[1] U. S. Department of the Interior, U. S. Geological Survey. Hydro1k[EB/OL]. (2015-01-30)[2015-07-16]. https:∥ lta.cr.usgs.gov/HYDRO1K.

[2] Cgiar-Csi. SRTM 90 m digital elevation data[EB/OL]. (2008-08-19)[2015-07-16]. http:∥srtm.csi.cgiar.org.

[3] Systems Japan Space. Aster gdem[EB/OL]. (2012-12-30)[2015-07-16]. http:∥www.jspacesystems.or.jp/ersdac/ GDEM/E/index.html.

[4] 劉 遠, 周買春, 陳芷菁, 等. 基于不同DEM數據源的數字河網提取對比分析——以韓江流域為例[J]. 地理科學, 2012,32(9):1 112-1 118.

[5] Quinn P, Bever K, Chevalier V, et al. The prediction of hillslope flow paths for distributed hydrological modeling using digital terrain models[J]. Hydrological Processes, 1991,5(1):59-79.

[6] 解河海, 黃國如. 地形指數若干計算方法探討[J]. 河海大學學報(自然科學版), 2006,34(1):46-50.

[7] 鄔 倫, 汪大明, 張 毅. 基于DEM的水流方向算法研究[J]. 中國圖象圖形學報, 2006,11(7):998-1 003.

[8] 雍 斌, 張萬昌, 陳艷華. TOPMODEL中地形指數ln(a/tanβ)的新算法[J]. 地理研究, 2007,26(1):37-45.

[9] O'Callaghan J F, Mark D M. The extraction of drainage networks from digital elevation data[J]. Computer Vision, Graphics, and Image Processing, 1984,28(3):323-344.

[10] 楊 光, 李慶和, 孫保平. 基于GIS的不同精度DEM提取地形指標誤差分析——以鹽池縣南部山區為例[J]. 水土保持研究, 2008,15(5):57-60.

[11] Wolock D M, Price C V. Effects of digital elevation model map scale and data resolution on a topography-based watershed model[J]. Water Resource Research, 1994,30(11):3 041-3 052.

[12] 秦福來, 王曉燕, 竇培謙, 等. 不同精度DEM流域地形特征提取分析及坡度誤差控制研究[J]. 國土資源遙感, 2005,(4):56-59.

[13] 吳險峰, 劉昌明, 王中根. 柵格DEM的水平分辨率對流域特征的影響分析[J]. 自然資源學報, 2003,18(2):148-154.

[14] 張微微, 武 偉, 劉洪斌. 不同比例尺DEM提取地形信息的比較分析——以重慶市為例[J]. 西南大學學報(自然科學版), 2007,29(7):153-157.

[15] Takeuchi K, Ao T. Introduction of block-wise use of TOPMODEL and Muskingum-Cunge method for the hydro-environment simulation of a large ungauged basin[J]. Hydrological Sciences Journal, 1999,44(4):633-646.

[16] Ao T. Development of a Distributed Hydrological Model for large river basins and its application to Southeast Asian Rivers[D]. Kofu, Japan: University of Yamanashi, 2001.

[17] Takeuchi K, Hapuarachchi P, Zhou M, et al. A BTOP model to extend TOPMODEL for distributed hydrological simulation of large basins[J]. Hydrological Processes, 2008,22(17):3 236-3 251.

[18] U. S. Department of The Interior, Survey U S G. Global 30 arc-second elevation (GTOPO30)[EB/OL]. (2015-01-01) [2015-07-16]. https:∥lta.cr.usgs.gov/GTOPO30.