基于Landsat的粵港澳大灣區水面率演變分析

姜欣彤，廖小龍，全栩劍，薛 嬌，鐘逸軒，于百順

(中水珠江規劃勘測設計有限公司，廣東 廣州 510610)

關鍵字：水面率；演變分析；Landsat；粵港澳大灣區

水面通常指由河流(江、河、渠等)，湖泊(天然或人工湖泊)，水庫，濕地(天然或人工濕地)等形成的水體表面[1]，水面率是指承載水域功能的區域面積占區域總面積的比率。適宜的城市水面率對改善城市地區生態質量和提高城市生活空間質量具有至關重要的作用[2]。浙江省、江蘇省、上海浦東新區、廣州市等地均通過河道調查及全國水利普查結果開展過本區域水面率的相關研究[3]，隨著3S技術的迅速發展，RS與GIS相結合的方法，逐漸應用在水體信息監測中[4]。利用遙感數據進行水體信息提取的數據源主要包含雷達遙感數據、光學遙感數據以及兩者結合數據[5]。其中，基于雷達遙感SAR數據作為數據源的應用范圍主要為洪澇監測[6-7]、城市水體和開放水體研究[8-9]、河流及海岸帶形態變化監測[10]、各種提取算法研究[11]等；基于光學遙感數據的研究中，應用最多的是免費開放的Landsat數據，主要集中在大區域尺度的研究以及對較大河流、湖泊等水體的提取和分析上[12-13]。MODIS、AVHRR和CBERS等中、低分辨率衛星影像主要應用在海洋、大型河流和湖泊等水體的識別和信息提取方面[14-15]。

粵港澳大灣區是繼紐約灣區、舊金山灣區、東京灣區之后的世界第四大灣區，包括9市和港澳2個特別行政區，隨著“9+2”粵港澳城市群互聯互通格局的構建，迅速的城市化發展也為粵港澳大灣區的生態環境建設和維護提出了更高的要求[16-17]，但粵港澳大灣區位于珠江三角洲，該地區河網密布交錯，徑流潮流相互作用，水沙變化頻繁，既存在洪(潮)澇災害頻發的老問題，又在經濟社會發展過程中暴露出水資源供需矛盾、水生態損害、水環境污染等新問題，因此基于遙感影像對近20年來粵港澳大灣區城市群水面率進行演變分析，為粵港澳大灣區城市群水資源合理利用和水生態建設工作提供理論和方法支持，同時對推進粵港澳大灣區生態文明建設和加強環境保護具有重要的參考意義。城市水體提取受到高層建筑導致的陰影影響，相比山區、丘陵地區等區域的水體提取較為困難，因此，選取空間分辨率為30 m、免費開放的Landsat作為數據源。

1 研究區域概況

粵港澳大灣區位于廣東省(21°25′～24°30′N，111°12′～115°35′E)，其城市群由廣州、深圳、珠海、佛山、東莞、中山、惠州、江門、肇慶9個地級市和香港、澳門2個特別行政區組成，地處珠江下游，區域面積55 080 km2。區域氣候屬于南亞熱帶的海洋性季風氣候，年均溫22℃左右，夏季多雨且臺風頻繁，冬季冷暖變化無常。2017年末常住人口約6 956.93萬人，地區生產總值為10.18萬億元，折合美元1.6萬億元，位列四大灣區第二位，僅次于東京灣區。以粵港澳大灣區九市二區為研究范圍，由于澳門面積甚小，所以將澳門與珠海市做合并分析，研究區域概況見圖1。

2 數據資料與技術路線

2.1 數據來源

從美國地質調查局(https://earthexplorer.usgs.gov/)下載Landsat 7、Landsat 8系列衛星的相應圖像。2013年發射的Landsat 8是一種新型的Landsat衛星。該衛星攜帶2個傳感器，分別稱為操作陸地成像儀(OLI)和熱紅外傳感器(TIRS)。Landsat 8數據具有16位輻射計分辨率，每幅圖像覆蓋距離地面185 km×185 km。Landsat 7于1999年4月15日由德爾塔11導彈發射，其內置了ETM+傳感器。覆蓋粵港澳大灣區核心城市群的影像軌道號為[122，44]、[121，45]、[122，45]，共三景。篩選云量覆蓋低于20%的影像，選取2000、2005、2010、2015、2018年共5年數據，各年份數據具體日期見表1。城市群行政邊界采用2015年國家地市行政邊界，數據來源于中科院地理科學與資源研究所(http://www.resdc.cn/Default.aspx)。

表1 各年份數據源具體日期

2.2 技術路線

水體提取技術路線見圖2。

圖2 水體邊界提取技術路線

2.3 數據預處理

2003年以后的Landsat 7影像由于傳感器發生故障造成部分條帶受損，故首先在ENVI中對Landsat 7進行條帶修復。LandsatL1T數據文件已經經過了帶DEM的地形校正，坐標精度基本能滿足中小比例尺的要求，但還未做輻射定標，大氣校正是多光譜遙感數據進行地表參數定量分析的前提，能夠消除或減少大氣分子和氣溶膠的散射和吸收對地物反射率的影響，故要在ENVI5.3中對Landsat 8和修復好的Landsat 7數據完成輻射定標和大氣校正數據預處理工作。

3 研究方法

3.1 水體提取方法

目前，國內外學者基于Landsat遙感影像的水體光譜曲線特征提出了多種水體信息自動識別的遙感模型，主要可歸納為單波段閾值法、多波段法和水體指數法3種方法[18]。單波段閾值法主要利用水體在近紅外/中紅外波段強烈吸收的原理來識別水體信息，難以忽視陰影造成的影響，會遺漏部分水體信息，存在一定局限性。多波段法主要是綜合利用遙感影像特征波段的光譜特征來識別水體信息，其中，譜間關系法通過分析地物在遙感影像原始波段及圖像轉換獲取的特征波段的光譜特征曲線來獲取識別水體信息的光譜信息，并利用邏輯判斷規則提取水體信息。水體指數法則基于譜間特征關系，通過特征波段間的比值運算，增強水體與其他地物之間的反差度來識別水體信息。最常見的為基于綠光/近紅外波段構成的歸一化水體指數NDWI(Normalized Difference Water Index)[19]，該指標可以突出影像中的水體信息，并能最大程度抑制植被信息，水體具有正值，而植被和土壤為0或負值，計算公式為:

(1)

式中 Green——綠光波段；NIR——近紅外波段。

然而NDWI指數忽略了土壤/建筑物因素的干擾，土壤/建筑物在綠光和近紅外波段的波譜特征與水體幾乎一致。因此考慮引入另一指標作為水體判別條件：歸一化植被指數NDVI(Normalized Difference Vegetation Index)[20]。NDVI是植被監測中最重要的指標之一，可以通過該指標區分植被、土壤和水體等不同覆蓋物，大于0.2表示植被，小于0表示水體，其計算公式為：

(2)

式中 Red——紅光波段；NIR——近紅外波段。

采用決策樹模型作為自動水體判別模型，以此對圖像進行分類，按照該規則分類的值被分類到一個分支中，其他值被分類到另一個類別中。該方法可以較好地提取復雜背景中的水體信息，有效地區分水體與山體陰影，提取結果比單閾值法精度要高[21]。在城市水體提取中，NDWI指標難以區分高樓建筑物，因此給NDWI設定一個閾值，以作為判別建筑物和水體的基準，經反復多次試驗(0，0.1，0.2，0.3)后，選定NDWI的閾值為0.2，故使用圖3所示判別規則進行水體提取。

圖3 水體提取判別法則

使用NDWI和NDVI作為判別指標進行水體提取后，圖像中存在許多噪點，故采用眾數濾波方式對圖像進行平滑處理，將處理后的圖像轉化為二值圖，即除水體外的像元設置為空值。在Google Earth中使用目視解譯法對提取后的水體進行校驗，對誤提的陰影、草地等部分進行刪除，以修正水體提取結果。最后將提取結果轉化為矢量圖，以統計水體面積。

3.2 土地利用轉移矩陣

土地利用轉移矩陣能反映在研究時段內初期和末期各地類相互轉化的動態過程，包括各地類相互轉化的面積及轉化方向[22]。其數學表達式為:

(3)

式中A——面積；n——土地利用類型數；i、j——研究初期和末期的土地利用類型；A12——研究初期至末期第1種土地類型轉為第2種土地類型的面積；Ann——研究初期至末期某種土地類型轉化為自身的面積。

4 分析與討論

4.1 水域面積變化

按照上述方法，計算分析得到各基準年粵港澳大灣區水域空間分布(圖4)，水域總面積變化(圖5)，大灣區城市群水面率變化(圖6)，2000、2018年粵港澳大灣區城市群水域面積對比(表2)。根據計算結果分析，可得出以下結論。

a)粵港澳大灣區近20年水域面積總體呈先減少后增加的趨勢，水域總面積從2000年2 841.94 km2下降到2015年2 629.83 km2，水面率由5.17%減少到4.79%。隨后水域面積略有增加趨勢，到2018年，水域總面積為2 664.89 km2，水面率為4.85%。

b)對比2000、2018年粵港澳大灣區城市群水域面積，發現總體上來說，各城市近20年來均呈減少趨勢，其中江門市水域面積減少最多，減少了44.20 km2，其次是佛山市，水域面積減少32.34 km2，珠海、東莞、深圳、中山、廣州水域面積依次減少了27.21、19.57、16.84、13.57、12.07 km2，其余城市水域面積減少不超過10 km2。

c)從2000—2018年各城市水面率變化分析來看，珠海(含澳門)、佛山兩市的水面率持續減少，珠海市水面率減少最為嚴重，從2000年19.02%減少至2018年17.11%，減少率為-0.11%/a。佛山市水面率從2000年13.17%減少至2018年12.32%，減少率為-0.05%/a。其余各城市水面率均在2015年后略有上漲。

a)2000年

b)2005年

圖5 粵港澳大灣區水體面積和水面率變化趨勢

圖6 粵港澳大灣區城市群水體面積變化趨勢

表2 2000、2018年粵港澳大灣區城市群水域面積對比 單位：km2

4.2 土地利用轉移變化

選用300 m精度的土地利用(LUCC)遙感監測數據，時間選取2000、2005、2010、2015、2018年共5期，數據來源于歐洲航天局(European Space Agency)。土地利用轉移分析結果見表3。粵港澳大灣區各市對應水域轉移方向分析結果見表4。根據計算成果分析，得出以下幾點結論。

a)粵港澳大灣區土地利用主要類型分為耕地、林地、草地、灌木、濕地、水域、建筑地和裸地，其中占比最高的類型為耕地，占總土地面積的42%，其次是林地，占總土地面積的36%。

b)近20年來耕地與草地面積均有大幅度減少，減少面積占國土總面積比例分別為6.35%、1.30%；建筑地面積大幅度增加，增加面積占國土總面積比例7.91%。

c)不同土地利用類型內部之間有相互轉化。各類型用地主要的轉換方向均為建筑用地。耕地和草地轉向建筑用地面積分別為3 472、721 km2，濕地、水域向建筑用地轉化面積分別為39、88 km2。

d)2000年來粵港澳大灣區各市水域面積減少主要原因是由于城市擴張，建筑用地侵占了水域范圍，導致城市水面率下降。江門市水域面積減少主要是轉移為濕地，而珠海市水域面積減少原因主要是轉移為耕地和濕地，建筑用地侵占面積較小。

表3 粵港澳大灣區2000—2018年土地利用轉移矩陣 單位：km2

表4 粵港澳大灣區各市水域轉移方向 單位：km2

5 結論與展望

在考慮到數據源可免費獲取，且精度相對較高的條件下，選擇Landsat7和8作為水面提取分析的數據源。考慮到城市水體提取中存在高層建筑、植被土壤等較多影響因素，同時選取水體指數和植被指數閾值作為地物分類基準，對粵港澳大灣區城市群主要河流、湖泊和較大型水庫進行了有效識別，得到主要結論如下：①粵港澳大灣區近20年水域面積總體呈先減少后增加的趨勢，水面率由2000年4.48%減少至2015年的4.15%，2018年增加至4.22%;②粵港澳大灣區各城市水域面積近20年來均有不同程度的減少趨勢，江門和佛山市水域面積減少最多，依次減少了44.20、32.34 km2，東莞、深圳、廣州水域面積減少不超過20 km2，其余城市水域面積減少不超過10 km2；③近20年來澳門、佛山兩市的水面率持續減少，澳門市水面率減少最為嚴重，減少率為-0.24%/a，佛山市減少率為-0.04%/a，其余各城市水面率均呈現先減后增的趨勢，在2015年后有所上漲；④造成粵港澳大灣區水域面積減少的主要原因是由于城市擴張，建筑用地侵占了水域范圍，導致城市水面率下降，江門水域面積減少是由于水域轉移為濕地，珠海水域面積減少則是由于水轉移為耕地和濕地。

值得注意的是，由遙感影像提取得到的水域面積數據，與土地利用數據中的水域面積數據有所偏差，這可能是由于圖像精度、云體、建筑覆蓋等不確定性因素。為得到更為清晰的水域邊界，可進一步考慮更為精細化尺度的遙感數據源對細小尺度的水體進行提取，將多源遙感數據提取結果進行融合，以計算更為精確的水面率。