李子彧++甘文霞++蔡曉斌

摘 要：實時獲取水體空間覆蓋范圍，對于災害的防治、水資源的合理調配具有重要的意義。基于數(shù)字高程模型和實測水位數(shù)據(jù)模擬水體覆蓋范圍的方法具有簡單易實現(xiàn)的特點，已有研究通常采用單一水位值作為參數(shù)，該文考慮水位的空間分布，對現(xiàn)有方法進行改進，提出了基于實測水位空間插值的水體淹沒范圍模擬方法。通過在鄱陽湖入江水道這一水體動態(tài)特征明顯的區(qū)域進行水域范圍模擬實驗，對方法的精度進行了綜合評價。結果表明該文方法能夠有效模擬水體范圍，相比現(xiàn)有方法模擬精度提高，當水位空間分布差異較大時，這一優(yōu)勢更為明顯。這一研究成果對水資源綜合管理相關研究與實際工作具有較為重要的實用參考價值。

關鍵詞：水體覆蓋范圍模擬 數(shù)字高程模型 空間分布 鄱陽湖

中圖分類號：P208 文獻標識碼：A 文章編號：1674-098X（2017）06（b）-0095-04

實時獲取水體空間覆蓋范圍，可以更好地了解對象的水文過程，為政府決策提供參考信息，對于水資源的綜合管理、水資源合理調配、洪旱災害的防治具有重要的意義。長期以來，對水文情勢的觀測與描述主要是基于傳統(tǒng)水文測驗手段，但水文站數(shù)量少，觀測范圍小，觀測結果僅能反映水文站所在區(qū)域的水文情況，對于整個區(qū)域的水文情勢的把握有其不足之處。基于衛(wèi)星遙感觀測技術可以大尺度獲取地表信息，但受傳感器時間、空間分辨率限制以及天氣因素的影響，經常無法獲得目標時刻的水體覆蓋范圍。

對于給定水體，水域的空間分布與水位之間有其內在聯(lián)系，都是從其特定角度出發(fā)對水量和水體形態(tài)的反映，從這個思路出發(fā)，出現(xiàn)了采用數(shù)字高程模型（Digital Elevation Model，DEM）與實測水位相結合的水體覆蓋范圍模擬方法，以水位數(shù)據(jù)作為閾值，與DEM柵格對應高程值比較，確定該柵格點是否被淹沒，這種基于柵格的方法數(shù)據(jù)要求少、計算簡單、易實現(xiàn)且能獲得精度較高的水體覆蓋范圍模擬結果。當前研究中通常使用代表水文站點水位或者平均水位作為閾值進行模擬。然而，水位具有空間異質性，而非空間不變，現(xiàn)有方法中缺少對這一實際特點的考慮。

基于此，該文對現(xiàn)有方法進行改進，提出了顧及水位空間特征的水體淹沒范圍模擬方法。以鄱陽湖入江水道這一高動態(tài)水域作為實驗區(qū)域，通過實驗對方法的精度和與現(xiàn)有方法的差異進行了綜合評價。

1 方法

1.1 基于DEM和水位空間插值的水體動態(tài)范圍模擬

該文的方法通過水文站點的水位數(shù)據(jù)進行空間插值得到模擬的水位空間分布，結合水位、淹沒源和DEM地形數(shù)據(jù)，采用有源淹沒算法判斷淹沒區(qū)域。方法可以概括為以下4步。

（1）數(shù)據(jù)獲取。獲取待模擬區(qū)域的數(shù)字高程模型、待模擬時期的實測水位數(shù)據(jù)。

（2）水體淹沒源區(qū)獲取。待判斷點是否與源區(qū)連通是判斷點是否被淹沒的必須條件。因此需要提前獲取淹沒源區(qū)。

（3）實測水位空間插值。采用空間插值方法將空間上離散的點狀水位數(shù)據(jù)空間化，獲取待模擬區(qū)域的水位柵格圖。該文中選擇反距離權重法進行水位的空間插值，它以插值點到已知點之間的距離為權重進行加權平均。

（4）采用有源淹沒算法，以第（3）步獲取的水位數(shù)據(jù)作為閾值，如果柵格的高程值低于給定的水位閾值且與已淹沒點聯(lián)通，則認為該點被淹沒，反之不被淹沒。

1.2 有源淹沒算法實現(xiàn)

當前有源淹沒算法的計算機實現(xiàn)通常采用堆棧法，按照水位值初始化淹沒源堆棧，之后向鄰域8個方向搜索，該文結合柵格DEM的固有特征，提出一種高速有源淹沒算法——數(shù)字柵格外推式掃描法，算法思想為：首先將所有淹沒源區(qū)的點集放入堆棧，然后先沿著x方向掃描已經淹沒區(qū)域柵格點左右兩個方向鄰域點，將高程低于水面的柵格點加入集合；然后沿著y軸方向，掃描集合內所有點上下方向鄰域點，得到新增被淹沒點并將其加入集合；重復上述步驟對集合內所有點進行操作，當集合不再發(fā)生變化時，所有集合的點即為最終淹沒區(qū)域。相比現(xiàn)用的八鄰域掃描法，該文的數(shù)字柵格外推式掃描法采用x、y方向交叉掃描，算法效率可以提高近16倍。

1.3 精度驗證

大量研究表明，基于遙感手段能夠準確的提取水體，獲取水體覆蓋范圍，提取精度能達到90%以上。該文以遙感提取結果為真實值，通過模擬準確比例、模擬錯誤比例值進行統(tǒng)計，對模擬結果進行了綜合精度評價。

1.4 衛(wèi)星遙感水體范圍提取方法

該文采用改進歸一化水體指數(shù)（Modified Normalized Water Index，MNDWI）閾值分割法從衛(wèi)星遙感影像數(shù)據(jù)中提取水體覆蓋范圍[8]。水體的MNDWI值較低而植被、土壤等地表其他地物的MNDWI值較高，設置恰當?shù)拈撝导纯蓞^(qū)分水體和其他地物，提取出水體覆蓋范圍。

2 實驗與結果

2.1 研究區(qū)域

鄱陽湖地處江西省北部、長江中下游南岸，是中國第一大淡水湖，湖區(qū)長期以來飽受水旱災害的困擾[9]。鄱陽湖區(qū)以松門山島為界，可以分為東（南）部主湖區(qū)和西（北）部入江水道。湖底自東向西、由南向北傾斜，高程由最高處的12 m降至湖口1 m左右[10]。流域內水域年內變化顯著，該文選取鄱陽湖入江水道作為研究區(qū)域（圖1），有助于綜合評估方法在不同水文情況下的模擬效果。

2.2 實驗數(shù)據(jù)

（1）實測水位數(shù)據(jù)。該文中所使用的水文數(shù)據(jù)來自江西省水利規(guī)劃設計院，各站刊布水位值轉換至吳淞基準高程。

（2）湖底地形DEM數(shù)據(jù)。該文采用長江水利委員會2000年1∶2.5萬的湖底地形圖數(shù)據(jù)，通過數(shù)字化、TIN模型構建，內插成30 m空間分辨率的柵格DEM數(shù)據(jù)，采用1985年國家高程系統(tǒng)。

（3）衛(wèi)星遙感數(shù)據(jù)，該文選用了Landsat TM/ETM+傳感器所獲取的晴空條件下的遙感影像，共12景（見表1），數(shù)據(jù)來自美國NASA的美國地球資源觀測中心（Earth Resources Observation and Science Center，EROS），幾何校正精度在一個像素以內，采用WGS 84全球投影坐標系統(tǒng)。

2.3 水體覆蓋范圍模擬精度評價

根據(jù)水文站點的空間分布，該文選擇了湖口、星子、都昌這3個水文站的實測水文數(shù)據(jù)，以歷年遙感影像中所提取的最小水體覆蓋范圍作為淹沒源，對表1中所列日期水體覆蓋范圍進行了模擬，以遙感影像水體提取結果作為真實值對模擬結果精度進行了評價（見表2）。

從表中可以看出，模擬結果精度均在80%以上，精度最高達98.09%。圖2直觀地表現(xiàn)了模擬結果以及模擬精度的空間分布狀況。

可以很明顯地看到大多數(shù)淹沒區(qū)域已被成功模擬，整體來說，豐水期的結果優(yōu)于枯水期模擬結果，枯水期入江水道中間洲灘容易被遺漏，靠近湖口的區(qū)域內較易出現(xiàn)多提，可能是因為枯水期湖口水位加速下泄，導致插值出的水位想比實際水位偏高，進而導致已經出露的洲灘模擬時仍被判定為淹沒區(qū)域。

2.4 與現(xiàn)有方法的對比

此節(jié)對該文方法與現(xiàn)有方法（采用單一閾值進行作為參數(shù)）進行了對比實驗，表3是方法在豐水期和枯水期的模擬精度，表3直觀地表現(xiàn)了不同方法提取的結果空間上的對照分析。

總體來說該文方法由于現(xiàn)有方法，從表3可以看出方法的精度有所有提高，從圖3中可以發(fā)現(xiàn)，兩種方法所提取的枯水期的水體覆蓋范圍有較為明顯的差異，主要表現(xiàn)在入江水道中部靠近湖口的區(qū)域。該文方法可以準確捕捉到入江水道靠近湖口區(qū)域的洲灘的出露狀況，而使用單一水位作為閾值則無法模這一區(qū)域的淹沒狀況?？紤]的水位變化后，漏提明顯減少，提高了模擬精度。在豐水期，兩種方法提取結果并沒有特別大的差別，其原因可能是因為豐水期水位較為平穩(wěn)，水位空間插值結果差異較小且與平均水位相差較小。

綜上表明，實測水位空間插值結果能夠更準確地描述水位的分布特征，作為輸入參數(shù)進行水體淹沒范圍模擬可以有效地改善模擬精度，當水位空間差異較大時，這一優(yōu)勢更為明顯。

3 結語

基于實測水位和數(shù)字高程模型進行模擬是獲取水體覆蓋范圍的重要手段之一，對于這一問題的研究對于水文、水環(huán)境、防災減災等問題具有重要的作用。該文認為在這一過程中，需要考慮水位的空間差異，提出了一種改進的基于數(shù)字高程模型和實測水位數(shù)據(jù)的水體淹沒范圍模擬方法，采用由多站點水位空間插值結果作為水域范圍模擬的閾值參數(shù)?；跉v史水位數(shù)據(jù)的水位空間分布特征分析為這一改進提供了理論支撐，在鄱陽湖入江水道這一高動態(tài)水體的實驗表明，對水位進行插值得到水位面，替代單一水位作為水體覆蓋范圍模擬的閾值，可以較大程度地提高模擬精度，尤其是當水位存在大的空間分布差異時，效果更加明顯，這一結論可以為相關研究與實踐工作提供參考。該文在相關的研究內容取得了一定的進展，但仍然有很多需要改進與提高之處。衛(wèi)星測高技術作為一項測量水面高度的新技術，從20世紀70年代逐漸興起，相比傳統(tǒng)水文觀測，具有不受環(huán)境影響，觀測范圍大的特點，使用衛(wèi)星測高數(shù)據(jù)獲取空間密度更高的水位分布狀況，將可以對水位的空間分布提供更加詳盡和準確的描述，進一步提高模擬精度。

參考文獻

[1] 黃淑娥，鐘茂生.基于遙感和流域降水的鄱陽湖水體覆蓋模型的研究[A].2003遙感科技論壇暨中國遙感應用協(xié)會2003年年會[C].2003.

[2] 黃淑娥，王保生.江西98特大洪澇衛(wèi)星遙感監(jiān)測及鄱陽湖水體覆蓋模型的研究[J].江西農業(yè)大學學報， 1999，21（2）：245-248.

[3] 李香顏，陳懷亮，李有.洪水災害衛(wèi)星遙感監(jiān)測與評估研究綜述[J].中國農業(yè)氣象[J].2009，30（1）：102，108.

[4] 鄧雪清，張永生，李波，等.湖泊水下數(shù)字高程模型建模方法及應用模型[J].華東地質學院學報，2003，26（2）：192-196.

[5] 李春紅，任立良，左振魯，等.基于DEM的三峽區(qū)間洪水淹沒范圍模擬[J].水文，2005，25（1）：1-4.