基于高光譜影像湖泊葉綠素a濃度反演分析

樊志強(qiáng), 王志國,*, 黃平平, 董亦凡, 高志奇

基于高光譜影像湖泊葉綠素a濃度反演分析

樊志強(qiáng)1,2, 王志國1,2,*, 黃平平1,2, 董亦凡1,2, 高志奇1,2

1. 內(nèi)蒙古工業(yè)大學(xué)信息工程學(xué)院, 內(nèi)蒙古 010051 2. 內(nèi)蒙古自治區(qū)雷達(dá)技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 內(nèi)蒙古 010051

針對湖泊葉綠素a濃度的快速變化及水體富營養(yǎng)化問題, 以太湖水體為例, 利用太湖區(qū)域2019年10月、11月和12月的珠海一號(hào)高光譜影像, 采用波譜特征分析、主成分分析、皮爾遜相關(guān)系數(shù)、建立反演模型分析江蘇省太湖區(qū)域水體葉綠素a濃度的分布狀態(tài)以及時(shí)空變化狀況。通過主成分分析劃分與水體信息相關(guān)的成分, 提出采用珠海一號(hào)高光譜影像的中紅外波段建立耦合可見光和近紅外輻射信號(hào)的WCI模型。結(jié)果表明: (1)水體光譜曲線特征與主成分分析結(jié)果一致, 分別為光合作用波段、細(xì)胞散射波段和水體增溫波段。(2)改進(jìn)的耦合可見光與中紅外輻射信號(hào)的建立的WCI模型的反演精度最高, 均方根誤差RMSE為(1.4008 mg·L-1), 平均相對誤差(MRE)為10.69%。(3)2019年10月一12月間, 湖泊葉綠素濃度整體變化不大, 太湖葉綠素a濃度的遙感反演結(jié)果總體上呈現(xiàn)為湖西區(qū)域濃度較低, 湖東區(qū)域濃度較高的特征。

高光譜影像; 葉綠素a; 水體富營養(yǎng);模型; 遙感反演

0 前言

葉綠素a(Chlorophyll-a, Chl-a)是影響內(nèi)陸湖泊水色的重要物質(zhì)之一, 其濃度的變化反映了湖泊水質(zhì)的污染狀況, 是湖泊環(huán)境監(jiān)測的重要指標(biāo)[1]。葉綠素a濃度與藻類的種類、數(shù)量等密切相關(guān), 是評價(jià)水體的富營養(yǎng)化程度的重要參數(shù), 因此監(jiān)測葉綠素a濃度具有重要的現(xiàn)實(shí)意義。傳統(tǒng)的葉綠素a測定方法, 從野外采集水樣, 將其帶回實(shí)驗(yàn)室對其進(jìn)行過濾、萃取、測定, 從而確定葉綠素a濃度[2]此類方法的缺點(diǎn)是需要耗費(fèi)大量的人力、物力且無法實(shí)現(xiàn)大面積水域的同步測量[3]。

隨著遙感技術(shù)的發(fā)展, 傳感器空間分辨率和光譜分辨率的性能指標(biāo)進(jìn)一步提高, 通過遙感影像監(jiān)測可實(shí)現(xiàn)高效, 實(shí)時(shí), 大范圍的水質(zhì)參數(shù)定量反演。特別是星載高光譜遙感為內(nèi)陸水體監(jiān)測提供了有利條件, 國內(nèi)外學(xué)者開展了一系列的相關(guān)研究。Purandara等[4]建立了印度文伯納德湖Landsat-TM數(shù)據(jù)與懸浮泥沙間的關(guān)系, 并用現(xiàn)場數(shù)據(jù)進(jìn)行了驗(yàn)證; 韋安娜等[5]采用窮舉法建立了高分五號(hào)數(shù)據(jù)的三波段模型, 結(jié)果表明該模型反演葉綠素a濃度效果較好; 蓋穎穎等[6]采用機(jī)載海洋高光譜數(shù)據(jù)儀光譜數(shù)據(jù)對葉綠素a和懸浮物進(jìn)行反演, 獲得了研究區(qū)的葉綠素a和總懸浮物濃度的空間分布狀況; 孫小涵等[7]利用榮成灣Chl-a濃度實(shí)測值和HJ-1衛(wèi)星數(shù)據(jù), 對3種適用于黃海二類水體的Chl-a濃度反演算法的精度進(jìn)行了驗(yàn)證, 并分析了其濃度隨時(shí)空變化的情況; 郭宇龍等[8]提出了一種新的水體葉綠素a指數(shù)CSI(葉綠素a光譜指數(shù)), 結(jié)果表明CSI對水體葉綠素a濃度大小具有指示作用, 反演精度與三波段算法相近。由于不同研究區(qū)域的光學(xué)特性和所采用數(shù)據(jù)源的不同, 使得最優(yōu)模型結(jié)構(gòu)及參數(shù)存在較大差異, 文獻(xiàn)[9][10]分別對多種模型進(jìn)行比較, 以得出最優(yōu)模型來反演葉綠素a濃度。高光譜影像的光譜分辨率高, 在相同的波長區(qū)間中較多光譜影像有更多的波段可供選擇, 因此更適合用于反演葉綠素a濃度。文獻(xiàn)[11][12]分別基于高光譜影像建立了三波段反演模型, 結(jié)果表明反演效果較好。由上述可知, 國內(nèi)外學(xué)者基于不同遙感數(shù)據(jù)源對于不同自然環(huán)境下的水體進(jìn)行了葉綠素a濃度反演, 但是對于遙感影像特定的波段參與葉綠素a濃度反演的物理機(jī)理研究還不夠深入, 缺少必要的實(shí)證研究。

本研究以太湖水域?yàn)檠芯繀^(qū), 通過2019年11月19日的珠海一號(hào)高光譜影像和同步水質(zhì)化學(xué)實(shí)測數(shù)據(jù), 由主成分分析篩選出與湖泊藻類葉綠素濃度強(qiáng)相關(guān)的波段成分, 并采用對湖泊水體有增溫作用的中紅外波段建立遙感反演模型, 以實(shí)測數(shù)據(jù)為依據(jù)對模型反演結(jié)果精度進(jìn)行驗(yàn)證分析。最后對模型反演所得的葉綠素濃度分布狀況進(jìn)行整體趨勢分析, 以便對內(nèi)陸二類水體葉綠素a濃度的高光譜反演研究提供有價(jià)值的物理模型方法參考。

1 研究區(qū)概況與數(shù)據(jù)源

太湖位于江蘇省南部, 長江中下游地區(qū), 地處于東經(jīng)119°52′32′′—120°36′10′′、北緯30°55′40′′—31°32′58′′之間。太湖水域面積約為2388平方公里, 最大水深不足3 m, 湖岸線393.2公里。太湖常年受到季風(fēng)和臺(tái)風(fēng)的影響, 在風(fēng)浪的作用下, 湖中的沉積物容易發(fā)生再懸浮并使得水體要素組成發(fā)生改變, 進(jìn)而改變了水體主要的光學(xué)響應(yīng)特征, 增加了水體光學(xué)作用過程的復(fù)雜度, 導(dǎo)致一定程度降低了遙感反演模型的精度和普適性。因此, 對于光學(xué)特性復(fù)雜的內(nèi)陸水體, 目前并沒有較好的滿足時(shí)空變化的通用模型, 并且太湖是我國重點(diǎn)治理的富營養(yǎng)化水域之一, 因此研究設(shè)計(jì)反演太湖的葉綠素a濃度揭示其時(shí)空分布規(guī)律具有重要意義。

1.1 數(shù)據(jù)采集

研究檢索衛(wèi)星過境時(shí)間, 同步獲取2019年11月19日均勻分布于太湖區(qū)域的37個(gè)站點(diǎn)的葉綠素a濃度實(shí)測數(shù)據(jù), 各站點(diǎn)數(shù)據(jù)包含監(jiān)測站點(diǎn)的經(jīng)緯度、葉綠素a濃度/(mg·L-1)等。采樣期間選擇未發(fā)生水華且藻類分布比較均勻的水體, 對水樣的各項(xiàng)指標(biāo)進(jìn)行測定。采樣點(diǎn)分布位置如圖1所示, 采樣當(dāng)日天氣晴朗, 風(fēng)速基本保持在5 m·s-1左右。

統(tǒng)計(jì)采樣數(shù)據(jù)得出, 葉綠素a濃度最高值出現(xiàn)在太湖的錫東西南站點(diǎn), 為22.3 mg·L-1。最低值出現(xiàn)在胥湖心(HF)站點(diǎn), 為1.3 mg·L-1。葉綠素a濃度的均值是7.5 mg·L-1。37個(gè)采樣點(diǎn)包含了從低到高的各種葉綠素a濃度的水體, 能較好地代表整個(gè)湖泊的葉綠素a濃度分布變化數(shù)值范圍。

1.2 遙感數(shù)據(jù)與預(yù)處理

1.2.1 遙感數(shù)據(jù)

珠海一號(hào)高光譜數(shù)據(jù)的空間分辨率為10 m, 光譜平均分辨率為2.5 nm, 共有32個(gè)波段, 波譜范圍為400—1000 nm。本文選用的遙感數(shù)據(jù)為2019年11月19日過境太湖的兩景珠海一號(hào)高光譜數(shù)據(jù)。葉綠素a濃度反演方法所使用的波段均能在珠海一號(hào)數(shù)據(jù)中獲得, 從而為空間尺度的水體葉綠素a定量反演提供了較為可靠的影像數(shù)據(jù)。由于遙感影像受到水體與傳感器間大氣輻射傳輸過程影響, 因此需要對影像進(jìn)行相應(yīng)的預(yù)處理。

圖1 采樣點(diǎn)位分布圖

Figure 1 Sampling point distribution map

1.2.2 影像預(yù)處理

由于遙感影像成像過程中的復(fù)雜性, 傳感器接收到的能量與地物本身輻射的能量是有差異的。傳感器記錄輸出的能量信號(hào)包含了由于太陽位置和角度條件、大氣條件、地形影響以及傳感器本身的性能等所引起的各種失真, 這些失真會(huì)對影像的使用和理解造成影響, 必須加以校正或消除。本研究影像的預(yù)處理過程參考?xì)W比特高光譜數(shù)據(jù)使用說明, 主要包括輻射定標(biāo)、大氣校正、正射校正、影像鑲嵌、水域提取等過程。其中, 輻射定標(biāo)將珠海一號(hào)高光譜遙感影像的值轉(zhuǎn)化為輻射亮度值, 大氣校正用于剔除大氣及其他因素的影響獲得表觀反射率, 正射校正對影像成像過程中產(chǎn)生的幾何畸變進(jìn)行校正。由于一景遙感影像無法覆蓋整個(gè)太湖區(qū)域, 將兩景影像拼接并用水體邊界矢量裁剪得到研究區(qū)太湖水體的高光譜遙感影像。

2 模型方法

2.1 單波段模型

單波段模型通常選用特定波長處反射率值作為自變量, 進(jìn)而建立的線性或非線性反演模型[13]。葉綠素a在440 nm、685 nm、700 nm均與葉綠素a濃度相關(guān)性較高, 但是440 nm附近受非色素顆粒物和黃色物質(zhì)的影響較大, 所以685 nm和700 nm附近反射率值較為常用。

2.2 波段比值模型

波段比值模型通常使用近紅外波段較大反射率值與紅波段較小反射率值進(jìn)行組合來估算葉綠素a濃度[13], 即R(NIR)/R(RED)。

2.3 三波段模型

三波段模型的具體形式為:

2.4 NDCI模型

Maisha等早在1997年提出了一種葉綠素估算指數(shù)——?dú)w一化葉綠素指數(shù)(Normalized Difference Chlorophyll Index,),模型的表達(dá)式可以表示為:

2.5 WCI模型

王金梁等[15]通過借鑒植被葉綠素監(jiān)測指數(shù)(Meris Terrestrial Chlorophyll Index)的研究思路, 構(gòu)建了基于水體葉綠素a的監(jiān)測指數(shù)(Water Chlorophyll-a Index)。模型可用式(3)表示:

2.6 四波段模型

Le等[16]在三波段模型的基礎(chǔ)上進(jìn)行改進(jìn)得到了四波段模型。四波段模型采用高光譜數(shù)據(jù)的四個(gè)波段進(jìn)行組合, 近紅外波段的加入減弱了懸浮顆粒的吸收和后向散射影響。四波段模型的表達(dá)式如下:

3 結(jié)果與分析

3.1 光譜特征分析

圖2為2019年11月19日珠海一號(hào)高光譜影像對應(yīng)于采樣點(diǎn)提取出的光譜曲線。

從圖2中可以看出, 太湖的珠海衛(wèi)星遙感反射率曲線呈現(xiàn)典型的二類水體的特征, 各采樣點(diǎn)的光譜曲線具有相似的趨勢。各采樣點(diǎn)光譜曲線可大致分為三部分, 在446—500 nm范圍內(nèi)水體中的藻類進(jìn)行光合作用,

利用光能將二氧化碳轉(zhuǎn)化為儲(chǔ)存能量的有機(jī)物, 由反射率光譜可反映出水體中藻類的葉綠素含量; 在500—700 nm范圍內(nèi)水中藻類表現(xiàn)為弱吸收作用以及細(xì)胞的散射作用, 反射率的高低決定了葉綠素濃度的大小; 700 nm以后由于水對光的強(qiáng)吸收作用使得反射率較低, 同時(shí)700 nm以后的中紅外波段主要為水體增溫進(jìn)而提供給藻類生長繁殖的能量, 對湖泊水體的葉綠素a濃度具有指示作用。

3.2 波段分析

高光譜影像相較于多光譜波段具備更多波段集合, 波段輻射特征表達(dá)的冗余信息較多且波段間相關(guān)系數(shù)較高, 因此本文通過主成分分析將原始高光譜影像的32個(gè)波段精簡為幾個(gè)相關(guān)性低、方差解釋信息充分的成分。本研究通過主成分分析得出三個(gè)主成分。第一主成分主要反映藻類光合作用強(qiáng)度, 波段區(qū)間為B1—B5。第二主成分主要反映藻類的弱吸收和細(xì)胞的散射作用, 波段區(qū)間為B9—B16。第三主成分主要反映水體吸收光使水體增溫為藻類生長提供能量的作用, 波段區(qū)間為B19—B32。通過對波段因子進(jìn)行主成分分析解決了波段間的多重線性問題, 同時(shí)實(shí)現(xiàn)了數(shù)據(jù)降維, 將原本復(fù)雜的32個(gè)波段轉(zhuǎn)化為方差解釋主成分, 降低了建模的復(fù)雜性。

圖2 太湖采樣點(diǎn)的遙感反射率曲線

Figure 2 Spectrum Reflectance Curve of Taihu Sampling Site

皮爾遜相關(guān)系數(shù)用于描述兩個(gè)變量間的線性關(guān)聯(lián)度, 一般用表示,的取值在[-1,1]。若>0, 表明兩個(gè)變量是正相關(guān); 若<0, 表明兩個(gè)變量是負(fù)相關(guān)。計(jì)算公式為:

3.3 不同模型估算葉綠素a濃度分析

將上述確定的波段因子作為自變量, 葉綠素a濃度作為因變量, 兩者進(jìn)行回歸分析, 由各種算法估算的葉綠素a濃度反演模型如表1所示。由表中可以看出, 中紅外波段參與建模的模型、模型和波段比值模型決定系數(shù)較高, 分別為0.3967、0.4169和0.397。而單波段、三波段和四波段模型沒有中紅外波段的參與, 模型的決定系數(shù)較低, 分別為0.1212、0.0144和0.2471。由此可得, 中紅外波段參與反演葉綠素a濃度時(shí)具有更好的效果。

為了驗(yàn)證上述分析, 用剩余7個(gè)實(shí)測數(shù)據(jù)驗(yàn)證上述模型, 分別計(jì)算模型反演值與實(shí)測葉綠素a濃度值的均方根誤差()和相對誤差(), 計(jì)算公式如下:

表1 不同葉綠素a模型對比

Figure 3 Comparison of estimated and measured values of chlorophyll a in different models

3.3 葉綠素a濃度反演

將所建立的模型應(yīng)用于太湖2019年10月29日、2019年11月19日和2019年12月13日的珠海一號(hào)高光譜影像, 得到葉綠素a濃度的空間分布結(jié)果如圖4所示。反演結(jié)果顯示, 太湖葉綠素a濃度在10月份大部分湖區(qū)基本保持高濃度狀態(tài), 呈現(xiàn)出東高西低的分布狀態(tài)。由于湖西沿岸區(qū)域受人類活動(dòng)影響較小且水體含有的水生植物較少, 湖東區(qū)域水生藻類較多。因此反演結(jié)果與實(shí)際情況較為吻合。11月份和12月份由于氣溫的降低使得水中藻類生長緩慢且大部分藻類趨于老化甚至死亡。湖泊藻類受氣溫的影響葉綠素a濃度高值區(qū)域面積明顯減少。11月與12月葉綠素濃度分布狀況變化不大?？傮w上, 10月、11月和12月葉綠素a濃度的分布均呈現(xiàn)東高西低的狀態(tài)。

注: (a)10月份; (b)11月份; (c)12月份。

Figure 4 The results of inversion of chlorophyll a concentration from the hyperspectral data of Zhuhai-1

4 討論

通過前文反演結(jié)果可以看出, 江蘇省太湖區(qū)域葉綠素濃度高值出現(xiàn)在東部湖區(qū), 低值區(qū)域主要出現(xiàn)在西部湖區(qū), 根據(jù)江蘇省人民代表大會(huì)常務(wù)委員會(huì)在2010年修訂的《江蘇省太湖水污染防治條例》, 江蘇省對太湖流域?qū)嵭蟹旨壉Ｗo(hù), 即太湖湖體、沿湖岸五公里區(qū)域、入湖河道上溯十公里以及沿岸兩側(cè)各一公里范圍為一級保護(hù)區(qū); 主要入湖河道上溯10—50公里以及沿岸兩側(cè)各一公里范圍為二級保護(hù)區(qū); 其他地區(qū)為三級保護(hù)區(qū)。湖泊水體的污染主要由于人類活動(dòng)的影響和外源污染物質(zhì)的輸入所導(dǎo)致, 因此湖泊水體的污染狀況和防治措施應(yīng)得到重視。本文的葉綠素分布圖與朱云芳等人[19]采用神經(jīng)網(wǎng)絡(luò)得到的反演結(jié)果相一致, 在梅梁灣、竺山灣和湖心區(qū)北部區(qū)域葉綠素濃度較高。因此本文對江蘇省太湖水體葉綠素濃度的反演分析具有一定的科學(xué)性和合理性, 進(jìn)一步驗(yàn)證了WCI模型應(yīng)用于太湖的適用性。

江蘇省太湖區(qū)域是我國最大的淡水湖之一, 對太湖進(jìn)行污染防治是生態(tài)文明建設(shè)的標(biāo)志性工程。太湖水體屬于典型的內(nèi)陸Ⅱ類水體, 包含除葉綠素濃度外多種影響湖泊水體的其他物質(zhì), 如懸浮顆粒物、黃色物質(zhì)等。受自然環(huán)境的改變、人類圍湖造田以及破壞性的漁業(yè)活動(dòng)使得太湖魚類資源減少。因此, 經(jīng)濟(jì)發(fā)展與生態(tài)保護(hù)并行發(fā)展是十分必要的, 在湖泊葉綠素濃度高值區(qū)域應(yīng)減少污染物排放, 從而降低水體的營養(yǎng)物質(zhì)輸入, 可有效避免藍(lán)藻水華的大面積爆發(fā)。葉綠素濃度低值區(qū)域, 應(yīng)結(jié)合生態(tài)環(huán)境優(yōu)勢, 合理利用湖泊水資源, 保障生態(tài)系統(tǒng)的良好運(yùn)作, 才能最大程度實(shí)現(xiàn)綠水青山就是金山銀山。

本文采用珠海一號(hào)高光譜影像反演內(nèi)陸湖泊水體的葉綠素a濃度從而監(jiān)測湖泊水體的富營養(yǎng)化狀況。與文獻(xiàn)[20]和[21]重點(diǎn)分析可見光波段原始反射率和微分值與葉綠素濃度的關(guān)系有所不同, 本文采用對水體有增溫效果的中紅外波段參與建立反演模型并驗(yàn)證了模型的合理性和適用性, 為湖泊水體的大面積實(shí)時(shí)監(jiān)測、探究葉綠素濃度空間分布格局、制定生態(tài)系統(tǒng)服務(wù)保護(hù)措施提供了參考, 為實(shí)現(xiàn)湖泊生態(tài)環(huán)境保護(hù)和可持續(xù)發(fā)展提供了一定的科學(xué)依據(jù)。但是同時(shí)也存在一些不足之處, 珠海一號(hào)數(shù)據(jù)水色遙感波段的光譜響應(yīng)機(jī)理認(rèn)識(shí)還不夠深刻, 需進(jìn)一步探索相關(guān)光譜響應(yīng)特征, 如峰谷間距、峰谷位置、峰高等。研究這些特征與內(nèi)陸湖泊水體的關(guān)系, 從而使模型更好地解釋水中各個(gè)組分的光學(xué)機(jī)理。

5 結(jié)論

通過主成分分析實(shí)現(xiàn)波段降維, 篩選得出珠海一號(hào)高光譜影像與葉綠素濃度強(qiáng)相關(guān)且包含主要信息的3個(gè)波段成分, 分別是光合作用成分、細(xì)胞散射成分和水體增溫成分。其中, 光合作用成分主要是可見光波段, 水體增溫成分與中紅外波段有緊密關(guān)系。經(jīng)研究得出, 建立耦合可見光與中紅外波段的WCI反演模型可有效實(shí)現(xiàn)湖泊水體葉綠素濃度的反演, 反演結(jié)果滿足湖泊葉綠素a濃度監(jiān)測的需求。將所建立葉綠素濃度反演模型用于不同時(shí)間的影像均得到較好的結(jié)果, 因此本文研究所建模型的時(shí)空穩(wěn)定性較好。與此同時(shí), 本文研究的內(nèi)容與方法對于其他湖泊水體反演葉綠素a濃度同樣適用。

[1] 徐祎凡, 李云梅, 王橋, 等. 基于環(huán)境一號(hào)衛(wèi)星多光譜影像數(shù)據(jù)的三湖一庫富營養(yǎng)化狀態(tài)評價(jià)[J]. 環(huán)境科學(xué)學(xué)報(bào), 2011, 31(1): 81–93.

[2] 吳廷寬, 賀中華, 梁虹, 等. 貴州百花湖葉綠素a高光譜監(jiān)測模型研究[J]. 人民長江, 2017, 48(3): 23–27.

[3] ZHOU Qu, TIAN Liqiao, WAI O, et al. High Frequency Monitoring of Suspended Sediment Variations for Water Quality Evaluation at Deep Bay, Pearl River Estuary, China: Influence Factors and Implications for Sampling Strategy[J]. Water, 2018, 10(3): 323.

[4] PURANDARA B K, JAMADAR B S, CHANDRAMOHAN T, et al. Water Quality Assessment of a Lentic Water Body Using Remote Sensing: A Case Study[J]. Environmental Pollution, 2018: 371–380.

[5] 韋安娜, 田禮喬, 陳曉玲, 等. 基于窮舉法的鄱陽湖葉綠素a濃度高光譜反演模型與應(yīng)用研究——以GF-5衛(wèi)星AHSI數(shù)據(jù)為例[J]. 華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 191(3): 119–125.

[6] 蓋穎穎, 王章軍, 楊雷, 等. 金沙灘近岸水體葉綠素a和懸浮物遙感反演研究[J]. 國土資源遙感, 2020, 32(3): 129–135.

[7] 孫小涵, 胡連波, 馮永亮, 等. 基于HJ-1衛(wèi)星數(shù)據(jù)的榮成灣葉綠素a濃度時(shí)空變化特征分析[J]. 海洋湖沼通報(bào), 2018(5): 72-79.

[8] 郭宇龍, 李云梅, 王橋, 等. 一種新的水體葉綠素a指數(shù)及其應(yīng)用潛力分析[J]. 光學(xué)學(xué)報(bào), 2015, 35(11): 18–30.

[9] 黃啟會(huì), 賀中華, 梁虹, 等. 基于高光譜數(shù)據(jù)的百花湖葉綠素a濃度估算[J]. 環(huán)境科學(xué)與技術(shù), 2019, 42(1): 134– 141.

[10] 宋挺, 周文鱗, 劉軍志, 等. 利用高光譜反演模型評估太湖水體葉綠素a濃度分布[J]. 環(huán)境科學(xué)學(xué)報(bào), 2017, 37(3): 888–899.

[11] 楊煜, 李云梅, 王橋, 等. 基于環(huán)境一號(hào)衛(wèi)星高光譜遙感數(shù)據(jù)的巢湖水體葉綠素a濃度反演[J]. 湖泊科學(xué), 2010(4): 33–41.

[12] 杜聰, 王世新, 周藝, 等. 利用Hyperion高光譜數(shù)據(jù)的三波段法反演太湖葉綠素a濃度[J]. 環(huán)境科學(xué), 2009, 30(10): 2904–2910.

[13] 林劍遠(yuǎn), 張長興. 航空高光譜遙感反演城市河網(wǎng)水質(zhì)參數(shù)[J]. 遙感信息, 2019, 34(2): 23–29.

[14] EI-AIEM A, CHOKMANI K, LAURION I, et al. Comparative Analysis of Four Models to Estimate Chlorophyll-a Concentration in Case-2 Waters Using MODerate Resolution Imaging Spectroradiometer (MODIS) Imagery[J]. Remote Sensing, 2012, 4(8): 2373–2400.

[15] 王金梁, 秦其明, 李軍, 等. 基于高光譜數(shù)據(jù)的水體葉綠素a指數(shù)反演模型的建立[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(3): 128–134.

[16] LE Chengfeng, LI Yunmei, ZHA Yong, et al. A four-band semi-analytical model for estimating chlorophyll a in highly turbid lakes: The case of Taihu Lake, China[J]. Remote Sensing of Environment, 2009, 113(6): 1175–1182.

[17] 劉大召, 唐世林, 付東洋. 珠江口水體葉綠素a濃度高光譜反演研究[J]. 廣東海洋大學(xué)學(xué)報(bào), 2008, 28(1): 49–52.

[18] 蔣城飛, 廖珊, 付東洋, 等. 湛江港海域葉綠素a濃度的高光譜遙感反演[J]. 廣東海洋大學(xué)學(xué)報(bào), 2016, 36(6): 107–113.

[19] 朱云芳, 朱利, 李家國, 等. 基于GF-1 WFV影像和BP神經(jīng)網(wǎng)絡(luò)的太湖葉綠素a反演[J]. 環(huán)境科學(xué)學(xué)報(bào), 2017, 37(1): 130–137.

[20] 鄧實(shí)權(quán), 田禮喬, 李建, 等. 面向GF-5衛(wèi)星高光譜傳感器的渾濁水體葉綠素a濃度反演算法研究一以鄱陽湖為例[J]. 華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 179(3): 123–129.

[21] 孫昊, 周林飛. 石佛寺水庫葉綠素a濃度高光譜遙感反演[J]. 節(jié)水灌溉, 2019, 283(3): 67–70.

Inversion analysis of lake chlorophyll-a concentration based on hyperspectral images

FAN Zhiqiang1,2, WANG Zhiguo1,2,*, HUANG Pingping1,2, DONG Yifan1,2, GAO Zhiqi1,2

1. School of Information Engineering, Mongolia Industrial University, Hohhot 010051, Inner Mongolia, China 2. InnerMongolia Key Laboratory of Radar Technology and Application, Hohhot 010051, Inner Mongolia, China

In order to address the rapid change of chlorophyll-a concentration and eutrophication of lake water, Taihu Lake was taken as an example. The hyperspectral images of Zhuhai-1 in October, November and December 2019, spectral characteristic analysis, principal component analysis and Pearson correlation coefficient were used to establish an inverse model to analyze the distribution and temporal and spatial variation of chlorophyll-a concentration in Taihu Lake of Jiangsu Province. Through principal component analysis, the components related to water information were divided, and the WCI model coupling visible and near-infrared radiation signals was established by using the mid-infrared band of Zhuhai-1 hyperspectral image. The results are as follows. (1) The characteristics of water spectral curve were consistent with the results of principal component analysis, which were photosynthesis band, cell scattering band and water warming band. (2) The improved WCI model coupled with visible light, near-infrared and mid-infrared radiation signals had the highest inversion accuracy, with the root mean square error (RMSE) of 1.4008(mg·L-1) and the mean relative error (MRE) of 10.69%. (3) From October to December 2019, the overall chlorophyll concentration of the lake changed little. The remote sensing inversion results of chlorophyll-a concentration in Taihu Lake showed a low concentration in the west of the lake and a high concentration in the east of the lake.

hyperspectral image;chlorophyll-a; eutrophication; model; remote sensing inversion

樊志強(qiáng), 王志國, 黃平平, 等. 基于高光譜影像湖泊葉綠素a濃度反演分析[J]. 生態(tài)科學(xué), 2023, 42(1): 121–128.

FAN Zhiqiang, WANG Zhiguo, HUANG Pingping, et al. Inversion analysis of lake chlorophyll-a concentration based on hyperspectral images[J]. Ecological Science, 2023, 42(1): 121–128.

10.14108/j.cnki.1008-8873.2023.01.014

TP79

A

1008-8873(2023)01-121-08

2020-11-05;

2020-12-30

國家自然科學(xué)基金項(xiàng)目(61631011); 內(nèi)蒙古自治區(qū)直屬高校基本科研項(xiàng)目(JY20220209); 內(nèi)蒙古自治區(qū)科技計(jì)劃項(xiàng)目(2019GG138)

樊志強(qiáng)(1995—), 男, 烏蘭察布人, 碩士, 主要從事高光譜影像定量反演研究, E-mail: fanzhiqiang_95@163.com

王志國, 男, 博士, 講師, 主要從事環(huán)境信息化監(jiān)測技術(shù)研究, E-mail:jdywzg@163.com