劉建強(qiáng), 葉小敏, 陳 鋆

(1.國家衛(wèi)星海洋應(yīng)用中心, 北京 100081； 2.自然資源部空間海洋遙感與應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 北京 100081；3.江西水利規(guī)劃設(shè)計(jì)研究院, 南昌 330029)

湖泊水體是人類寶貴的可利用水資源，是全球水循環(huán)和生態(tài)系統(tǒng)極其重要的組成部分，在調(diào)節(jié)區(qū)域氣候和控制干旱、洪澇等災(zāi)害方面發(fā)揮著不可替代的作用[1].在全球變暖和人類活動(dòng)加劇的情況下，近幾十年來內(nèi)陸湖泊的數(shù)量、水位、面積都發(fā)生了劇烈的變化[2].因此，及時(shí)準(zhǔn)確地監(jiān)測(cè)湖泊的特征變化對(duì)有效管理水資源、評(píng)估人類活動(dòng)和氣候變化對(duì)湖泊的影響至關(guān)重要.水體范圍作為表征湖泊流域水文狀態(tài)的關(guān)鍵信息，在水資源調(diào)查、生態(tài)環(huán)境監(jiān)測(cè)、旱洪災(zāi)害評(píng)估等領(lǐng)域具有廣泛的應(yīng)用[3-5].傳統(tǒng)湖泊水體范圍主要依靠人工測(cè)量湖岸線進(jìn)行監(jiān)測(cè)，需要消耗大量人力物力且數(shù)據(jù)更新周期長(zhǎng)，難以滿足旱洪災(zāi)害的動(dòng)態(tài)跟蹤需求，相比人工監(jiān)測(cè)，遙感技術(shù)具有低成本、大范圍、高時(shí)空分辨率等優(yōu)勢(shì)，目前已在全球湖泊水體范圍監(jiān)測(cè)上得到廣泛應(yīng)用.現(xiàn)階段，已有大量研究基于Landsat系列衛(wèi)星TM數(shù)據(jù)、Terra/Aqua衛(wèi)星MODIS數(shù)據(jù)和環(huán)境減災(zāi)衛(wèi)星CCD數(shù)據(jù)等衛(wèi)星遙感多光譜數(shù)據(jù)，結(jié)合歸一化水體指數(shù)法、歸一化植被法等方法，實(shí)現(xiàn)水體范圍的識(shí)別與提取，并通過人工目視解譯、混淆矩陣等方式進(jìn)行精度評(píng)估，也開展了部分洪澇監(jiān)測(cè)應(yīng)用探索[4,6-15].

2018年9月7日發(fā)射成功的海洋一號(hào)C星(HY-1C)是我國第一顆業(yè)務(wù)化海洋衛(wèi)星，與2020年6月11日發(fā)射的海洋一號(hào)D星(HY-1D)可形成上、下午雙星組網(wǎng).HY-1C/D衛(wèi)星搭載的海岸帶成像儀(coastal zone imager，CZI)為四波段CCD相機(jī)，其星下點(diǎn)像元空間分辨率≤50 m，幅寬≥950 km，重訪周期3 d，具有較高分辨率，大幅寬等觀測(cè)優(yōu)勢(shì)[16]，已經(jīng)在海岸帶生態(tài)環(huán)境、應(yīng)急減災(zāi)、水環(huán)境監(jiān)測(cè)和農(nóng)業(yè)遙感等方面得到初步應(yīng)用，如，Cai等[17]基于HY-1C衛(wèi)星數(shù)據(jù)分析了島嶼對(duì)舟山近海懸浮泥沙濃度的影響，結(jié)果表明島嶼附近懸浮泥沙濃度變化明顯且有一定的規(guī)律性；梁超等[18]基于HY-1C CZI多光譜影像構(gòu)建的多個(gè)指數(shù)實(shí)現(xiàn)了海岸帶紅樹林的高精度自動(dòng)化提?。簧騺喎宓萚19]探索了HY-1C光學(xué)載荷對(duì)海面溢油的識(shí)別能力，并指出CZI能夠進(jìn)一步實(shí)現(xiàn)對(duì)溢油污染類型的區(qū)別；周屈等[20]利用HY-1C衛(wèi)星CZI數(shù)據(jù)對(duì)疫情期間武漢火神山、雷神山醫(yī)院附近知音湖、黃家湖水體濁度變化進(jìn)行監(jiān)測(cè)，初步證明了應(yīng)急醫(yī)院建設(shè)與運(yùn)營對(duì)鄰近湖泊水環(huán)境無擾動(dòng)影響；王利民等[21]基于HY-1C衛(wèi)星CZI數(shù)據(jù)實(shí)現(xiàn)了對(duì)冬小麥分布區(qū)的提取.目前還未有利用HY-1C/D衛(wèi)星組網(wǎng)的CZI多光譜影像開展內(nèi)陸湖泊水體面積較長(zhǎng)時(shí)序分析并服務(wù)于洪澇災(zāi)害的相關(guān)研究.

本文以中國最大的淡水湖——鄱陽湖洪澇監(jiān)測(cè)為例，利用HY-1C/D衛(wèi)星組網(wǎng)的CZI多光譜影像基于歸一化水體指數(shù)(normalized difference water index，NDWI)-歸一化植被指數(shù)(normalized difference vegetation index，NDVI)聯(lián)合方法開展水體范圍提取與面積時(shí)間序列分析，并結(jié)合水文站水位現(xiàn)場(chǎng)數(shù)據(jù)開展水面積-水文定量關(guān)系研究與洪澇風(fēng)險(xiǎn)分析.

1研究區(qū)域數(shù)據(jù)與方法

1.1研究區(qū)概況

鄱陽湖位于江西省北部、長(zhǎng)江中下游以南(28°22′N～29°45′N、115°47′E～116°45′E)，是長(zhǎng)江流域最大的通江湖泊(見圖1).鄱陽湖承納贛江、撫河、信江、饒河、修河等重要河流的來水，經(jīng)調(diào)蓄后由湖口注入長(zhǎng)江，同時(shí)其水位還受長(zhǎng)江徑流水位的影響，是一個(gè)過水性、吞吐型、季節(jié)性湖泊，對(duì)調(diào)節(jié)長(zhǎng)江徑流流量、維護(hù)區(qū)域生態(tài)平衡至關(guān)重要[22-24].鄱陽湖流域地處亞熱帶季風(fēng)氣候區(qū)，一半以上的降水集中在3月至6月，因此具有極大的季節(jié)性水體面積波動(dòng)(約2 300 km2).同時(shí)，鄱陽湖流域還是長(zhǎng)江中下游的重要水源地、重要的生物物種遺傳基因庫、我國最大的淡水魚產(chǎn)區(qū)、國際重要濕地和世界著名的候鳥越冬棲息地[3,25].鄱陽湖旱洪災(zāi)害頻發(fā)，對(duì)利用星地聯(lián)合多源數(shù)據(jù)進(jìn)行水體范圍動(dòng)態(tài)跟蹤提出了更高需求.

1.2衛(wèi)星影像與現(xiàn)場(chǎng)水位數(shù)據(jù)

1.2.1衛(wèi)星影像

1) HY-1C/D衛(wèi)星CZI數(shù)據(jù).本文使用的是2019年3月-2020年8月43景晴空無云經(jīng)過大氣瑞利校正后的L2A級(jí)HY-1C/D衛(wèi)星CZI反射率數(shù)據(jù).各月份的有效數(shù)據(jù)量基本相等.HY-1C/D衛(wèi)星CZI波段設(shè)置和性能指標(biāo)信息如表1所示.

2) GF-3號(hào)衛(wèi)星SAR數(shù)據(jù).本文采用2020年8月17日高分三號(hào)(GF-3)雙極化標(biāo)準(zhǔn)條帶SAR影像數(shù)據(jù)來評(píng)估近似同步的HY-1C衛(wèi)星CZI提取水體范圍結(jié)果的可靠性，其空間分辨率為25 m，幅寬130 km.

圖1 研究區(qū)域與星子水文站站位圖Fig.1 Study area and location of the Xingzi hydrological station

表1 HY-1C/D衛(wèi)星CZI波段設(shè)置與性能指標(biāo)Tab.1 Bands and specifications of CZI onboard HY-1C/D

1.2.2水位數(shù)據(jù)

為了構(gòu)建水體面積-水位之間的關(guān)系模型，從2019年3月10日至2020年8月25日的星子水文站水位數(shù)據(jù)(基于黃海基準(zhǔn)面的水位高度值)中遴選出同步的日水位觀測(cè)值.

1.3CZI水體識(shí)別方法

歸一化差分水體指數(shù)(normalized difference water index, NDWI)作為主要水體識(shí)別方法[30]，通過水體反射率最大的綠色波段和水體反射率最小的近紅外波段之間的比值，可有效增強(qiáng)陸地和水體像元之間的光譜差異，計(jì)算方法如公式(1)所示.通常情況下，NDWI提取水體像素選擇0作為分割閾值即可，但由于湖泊地區(qū)地物類型復(fù)雜多變，為了提高水體識(shí)別準(zhǔn)確率，同時(shí)采用歸一化差分植被指數(shù)(normalized difference vegetation index，NDVI)作為輔助信息，NDVI計(jì)算方法如公式(2)所示.

(1)

(2)

式中，ρG、ρR和分別為HY-1C/D衛(wèi)星CZI傳感器綠光(中心波長(zhǎng)為560 nm的第2波段)、紅光(中心波長(zhǎng)為650 nm的第3波段)和近紅外波段(中心波長(zhǎng)為825 nm的第4波段)的反射率.為了進(jìn)一步保證閾值的穩(wěn)定性，研究采用經(jīng)過大氣瑞利散射校正后的HY-1C/D衛(wèi)星CZI傳感器L2A級(jí)反射率數(shù)據(jù)計(jì)算NDWI和NDVI.結(jié)合目視解譯分析結(jié)果確定鄱陽湖HY-1C/D衛(wèi)星CZI影像水體像元識(shí)別標(biāo)準(zhǔn)為：NDWI≥0或者-0.2

1.4水體面積-水位定量關(guān)系模型

假定較短時(shí)間內(nèi)湖泊地形變化不大，則可通過分析水體面積和水位數(shù)據(jù)獲得其定量關(guān)系.為了獲得一般情況下鄱陽湖水體面積與水位的關(guān)系，可將少數(shù)異常數(shù)據(jù)剔除.采用最小二乘法對(duì)其水位和水體面積數(shù)值進(jìn)行線性擬合，得到的鄱陽湖地區(qū)水體面積-水位定量關(guān)系模型：

S=a·H+b，

(3)

式中，H為水位，單位：m；S為水體面積，單位：km2，a、b為最小二乘擬合系數(shù).

2結(jié)果與分析

2.1CZI鄱陽湖水體識(shí)別結(jié)果與精度評(píng)價(jià)

利用1.3節(jié)描述的數(shù)據(jù)處理方法和水體識(shí)別判定標(biāo)準(zhǔn)，進(jìn)行數(shù)據(jù)處理實(shí)現(xiàn)了HY-1C/D衛(wèi)星CZI影像鄱陽湖水體識(shí)別.圖2為2020年2月18日HY-1C衛(wèi)星與2020年7月17日HY-1D衛(wèi)星CZI傳感器獲取的鄱陽湖遙感影像(R、G、B=3、2、1)及提取的水體范圍疊加圖，由圖2可見，無論枯水期(圖2(a))還是豐水期(圖2(b))，HY-1C/D衛(wèi)星CZI影像均能得到較好的水體范圍識(shí)別結(jié)果.

圖2 (a) 2020年2月18日(北京時(shí)間：10∶53) HY-1C衛(wèi)星與(b)2020年7月17日(北京時(shí)間： 13∶45)HY-1D衛(wèi)星CZI傳感器獲取的鄱陽湖遙感影 像(R、G、B=3、2、1)及提取的水體范圍疊加圖Fig.2 February 18, 2020 (Beijing time: 10∶53) HY-1C satellite (a) and July 17, 2020 (Beijing time: 13∶45) HY-1D satellite (b) CZI sensor acquired Poyang Lake remote sensing Image (R, G, B=3, 2, 1) and the superimposed map of the extracted water bodyrange

為進(jìn)一步驗(yàn)證HY-1C/D衛(wèi)星CZI影像水體范圍識(shí)別結(jié)果的精度，選取同步的GF-3衛(wèi)星SAR影像數(shù)據(jù)的識(shí)別結(jié)果來進(jìn)行對(duì)比驗(yàn)證.星載合成孔徑雷達(dá)(synthetic aperture radar，SAR)是主動(dòng)探測(cè)微波成像雷達(dá)，對(duì)云雨不敏感.由于以鏡面反射為主的水體和以散射為主的粗糙地表對(duì)微波反射信號(hào)存在較大差異，SAR影像上水陸像元存在明顯的亮暗差異，能夠顯著區(qū)分.已有研究表明，SAR影像能夠作為識(shí)別提取地表水體的可靠數(shù)據(jù)源和光學(xué)多光譜遙感影像的有效補(bǔ)充和驗(yàn)證[26-29].圖3(a)為2020年8月16日10∶51 HY-1C衛(wèi)星CZI在28.75°N～29.15°N，116.48°E～116.85°E范圍內(nèi)水體識(shí)別結(jié)果局部放大圖.圖3(b)為2020年8月17日06:19(北京時(shí)間)成像的GF-3衛(wèi)星標(biāo)準(zhǔn)條帶成像的C波段HH極化SAR影像及水體提取結(jié)果疊加圖.可以看出，利用CZI提取的水體范圍信息與SAR影像顯示的水體信息基本一致.結(jié)合前序目視解譯結(jié)果分析，CZI提取的面積為5 824.58 km2，SAR提取的面積為5 591.32 km2兩者相對(duì)差異為4.0 %.部分地區(qū)存在邊界差異可能與兩個(gè)衛(wèi)星傳感器的空間分辨率、獲取時(shí)間等方面差異有關(guān)(HY-1C/D CZI分辨率50 m，GF-3 SAR分辨率25 m，時(shí)間相差約1 d).

圖3 (a) 2020年8月16日10∶52(北京時(shí)間)HY-1C 衛(wèi)星CZI影像識(shí)別的水體范圍結(jié)果與 (b)2020年8月17日06∶19(北京時(shí)間)GF-3衛(wèi)星SAR影像水體范圍提取結(jié)果的比較Fig.3 Comparison of land and water recognition results from CZIdata with SAR image at 10∶52 (Beijing Time) on August16, 2020 by HY-1C Satellite CZI (a) and at 6∶19 (Beijing Time) on August17, 2020 by GF-3 Satellite SAR (b)

2.2鄱陽湖水體面積時(shí)序變化與洪澇風(fēng)險(xiǎn)分析

為了分析鄱陽湖水體面積變化與洪澇風(fēng)險(xiǎn)之間的關(guān)系，對(duì)自2019年3月10日至2020年8月25日共43個(gè)時(shí)相的HY-1C/D衛(wèi)星CZI數(shù)據(jù)進(jìn)行鄱陽湖地區(qū)的水體識(shí)別；經(jīng)過幾何配準(zhǔn)后疊加所有影像，并統(tǒng)計(jì)每個(gè)像元被識(shí)別為水體的次數(shù)，獲得了此間約18個(gè)月時(shí)間跨度的鄱陽湖水體淹沒范圍頻次分布圖，如圖4(a).在所有的43景影像中，通過頻率分布直方圖及目視解譯易知，鄱陽湖大面積水體的覆蓋頻次均在19次以上，這可能與季節(jié)性因素有關(guān).因此本文將被水淹沒頻次小于19次的區(qū)域判定為洪澇風(fēng)險(xiǎn)區(qū).如圖4(a)所示，洪澇風(fēng)險(xiǎn)區(qū)主要在鄱陽湖東南區(qū)域和西部區(qū)域.例如2020年7月汛期江西省鄱陽縣發(fā)生了洪災(zāi)，蓮湖鄉(xiāng)和昌州鄉(xiāng)大面積土地被洪水大面積淹沒.對(duì)比2019年8月16號(hào)和2020年7月17號(hào)兩景CZI影像的局部放大圖(圖4(b))，可以明顯看出洪災(zāi)的影響范圍，通過對(duì)比計(jì)算得到該次洪災(zāi)的洪水淹沒面積約為1 320 km2.

2.3鄱陽湖水體面積-水位關(guān)系分析

通過對(duì)比水體面積和水位波動(dòng)的結(jié)果，可以發(fā)現(xiàn)鄱陽湖水體面積和水位變化呈現(xiàn)極為一致的增減規(guī)律，即水位升高，則同步的水體面積增大；水位降低，同步水體面積減小(見圖5).在個(gè)別水位日期升高或降低的情況中，水體面積變化相對(duì)于滯后，可能與個(gè)別低洼地區(qū)或汛期洪水淹沒地區(qū)有關(guān).低洼的內(nèi)澇區(qū)及洪災(zāi)淹沒區(qū)，在水位降低時(shí)部分積水的地區(qū)排洪不及時(shí)；或者水位升高時(shí)，由于水體漫過防洪堤壩，低洼地區(qū)被突然淹沒所致.例如，2020年7月9日鄱陽縣昌洲鄉(xiāng)汛期期間決堤后，洪水無法在水位降低時(shí)及時(shí)排出，從而導(dǎo)致2020年7月鄱陽湖水位較快降低時(shí)，水體面積卻緩慢減小.此外，從提取結(jié)果的時(shí)序變化來看，鄱陽湖水體面積具有明顯的季節(jié)特征.9月中旬至1月，水體面積最小，約1.6×103km2；2月至6月，水體面積增大至2×103～3×103km2；7月和8月，水體面積最大，最大面積達(dá)4.7×103km2.鄱陽湖作為一個(gè)通江湖泊，提取結(jié)果明顯表征其水體面積季節(jié)性變化劇烈.

采用最小二乘方法對(duì)水體面積和水位數(shù)據(jù)進(jìn)行線性擬合，擬合系數(shù)計(jì)算為：a=293.10,b=-787.49，進(jìn)而得到鄱陽湖水體面積-水位定量關(guān)系模型:

S=293.10H-787.49，

(4)

鄱陽湖水體面積和水位相關(guān)性較高，其相關(guān)系數(shù)R高達(dá)0.97.該結(jié)果表明鄱陽湖水體面積與水位可以根據(jù)兩者擬合關(guān)系互相補(bǔ)充.有實(shí)測(cè)水位時(shí)可以在一定程度上估算同期鄱陽湖水體面積，為云雨霧環(huán)境條件下湖泊水體面積監(jiān)測(cè)提供參考.在極端情況下，甚至可以根據(jù)水文站水位數(shù)據(jù)計(jì)算出水體面積，也可以基于不同區(qū)域的淹沒頻率粗略判斷出潛在的淹沒區(qū)域.同時(shí)，也可由HY-1C/D遙感數(shù)據(jù)所獲得的水體面積估算鄱陽湖的水位，即可通過遙感手段間接估算其水位值.當(dāng)然，由于受人類活動(dòng)或自然災(zāi)害(如洪澇災(zāi)害引起決堤等)等影響，地形可能發(fā)生改變，則其水位與水體面積關(guān)系也可能發(fā)生相應(yīng)改變.

圖4 (a)2019年3月-2020年9月HY-1C/D衛(wèi)星CZI傳感器獲取的鄱陽湖水體淹沒范圍頻次分布圖及 (b)淹沒范圍局部放大圖(以2019年8月16日和2020年7月17日為例)Fig.4 (a) Frequency distribution map of the submerged area of Poyang Lake obtained by the HY-1C/D satellite CZI sensor from March 2019 to September 2020 and (b) Partial enlarged view of inundation range (take August 16, 2019 and July 17, 2020 as examples)

圖5 HY-1C/D衛(wèi)星CZI數(shù)據(jù)監(jiān)測(cè)的鄱陽湖地區(qū)水體面積與水位實(shí)測(cè)值隨時(shí)間變化曲線Fig.5 Time series of water area and water level in Poyang Lake area monitored by CZI data of HY-1C and HY-1D Satellites

3結(jié)語

本文通過HY-1C/D衛(wèi)星CZI影像數(shù)據(jù)計(jì)算歸一化水體指數(shù)(NDWI)和歸一化植被指數(shù)(NDVI)，并建立像元識(shí)別規(guī)則獲取了2019年3月至2020年8月鄱陽湖水體面積.結(jié)合星子站的水位數(shù)據(jù)分析了鄱陽湖水體面積-水位定量關(guān)系，識(shí)別出鄱陽湖洪澇風(fēng)險(xiǎn)區(qū)域.可得到如下結(jié)論：1) HY1-C/D衛(wèi)星數(shù)據(jù)重訪周期短、幅寬大、信噪比優(yōu)良，可快速獲取圖像細(xì)節(jié)特征清晰的光學(xué)遙感數(shù)據(jù)，具有沿海信息監(jiān)測(cè)、內(nèi)陸水體提取的潛力和價(jià)值，拓展了國產(chǎn)衛(wèi)星的應(yīng)用領(lǐng)域，為洪澇災(zāi)害監(jiān)測(cè)提供了新的備用數(shù)據(jù)源；2) 鄱陽湖水體面積季節(jié)性變化較大，易發(fā)生洪澇區(qū)域主要位于湖東南區(qū)域和西部區(qū)域；3) 鄱陽湖地區(qū)的水體面積與水位具有高相關(guān)性，該探索對(duì)于鄱陽湖水體面積的星地聯(lián)合動(dòng)態(tài)監(jiān)測(cè)具有一定的積極意義.