曹萌萌，青 松，杜雨春子，袁瑞強(qiáng)，順布日

(內(nèi)蒙古師范大學(xué) 地理科學(xué)學(xué)院，內(nèi)蒙古 呼和浩特 010022)

1 研究背景

隨著社會(huì)發(fā)展和人類活動(dòng)強(qiáng)度的增加，我國(guó)乃至全世界的生態(tài)環(huán)境遭到破壞，其中水體富營(yíng)養(yǎng)化引發(fā)的藻華現(xiàn)象已成為社會(huì)最關(guān)注的問(wèn)題之一[1]。藻華不僅破壞湖泊景觀，導(dǎo)致生物量減少，同時(shí)也會(huì)造成水體自凈能力失衡，嚴(yán)重影響湖泊水生態(tài)環(huán)境。因此，對(duì)湖泊藻華進(jìn)行實(shí)時(shí)動(dòng)態(tài)監(jiān)測(cè)是預(yù)防與治理藻華的關(guān)鍵[2-4]。

與傳統(tǒng)方法相比，遙感技術(shù)不僅可以有效地繪制藻華的空間分布，也可以實(shí)時(shí)監(jiān)測(cè)其生長(zhǎng)情況，所以國(guó)內(nèi)外學(xué)者采用了多種方法對(duì)藻華進(jìn)行提取分析。其中，光譜指數(shù)法被廣泛應(yīng)用于藻華識(shí)別研究，例如：NDVI指數(shù)[5-6]、FAI指數(shù)[7-9]、AFAI指數(shù)[10]、EVI指數(shù)[11]、NDICB指數(shù)[12]等。同時(shí)，也有學(xué)者進(jìn)一步結(jié)合多個(gè)光譜指數(shù)建立了分類決策樹(shù)，并用來(lái)提取復(fù)雜水體環(huán)境中的藻華[13-15]。此外，除了常見(jiàn)的監(jiān)督或非監(jiān)督分類方法[16-18]，還有少部分研究結(jié)合了野外實(shí)測(cè)光譜數(shù)據(jù)[19]、目視解譯方法、水質(zhì)參數(shù)反演法[20]和專家系統(tǒng)[21]對(duì)藻華進(jìn)行識(shí)別。綜合對(duì)比以上幾種分類方法，發(fā)現(xiàn)它們各自都具有缺陷。結(jié)合光譜指數(shù)的決策樹(shù)分類算法缺乏伸縮性，難以針對(duì)大量訓(xùn)練數(shù)據(jù)集中確定唯一的有效閾值，對(duì)缺失值敏感，容易過(guò)擬合，從而加大誤差。監(jiān)督分類方法的缺點(diǎn)在于人為主觀性強(qiáng)、選擇樣本的成本較高，而且僅限于識(shí)別已定義的類別。非監(jiān)督分類方法的分類精度無(wú)法保證，往往需要后期大量的分析處理。專家系統(tǒng)則在知識(shí)獲取和量化等方面頗為困難。此外，由于獲取野外實(shí)測(cè)數(shù)據(jù)與高分辨率影像的成本高、操作難度大，所以難以對(duì)藻華進(jìn)行高效率的識(shí)別。

自1995年粒子群優(yōu)化算法被提出后[22]，已成功應(yīng)用于洪水監(jiān)測(cè)[23]、水質(zhì)參數(shù)反演[24]、生物量估算[25]和土地利用變化[26]等多個(gè)方面。其中，基于離散粒子群優(yōu)化的光譜匹配(spectrum matching based on discrete particle swarm optimization, SMDPSO)算法是由Jia等[27]提出的一種用于提取地表水的方法，他們運(yùn)用該算法較精確地提取了全球8個(gè)地區(qū)不同類型的水體。與上文中提到的其他方法相比，SMDPSO算法具有精度高、成本低、參數(shù)少、先驗(yàn)知識(shí)少和自動(dòng)化潛力大的優(yōu)點(diǎn)。但目前基于SMDPSO算法的藻華監(jiān)測(cè)研究鮮有報(bào)道，而且在以往的藻華研究中，很少涉及到中國(guó)東北部冰封期較長(zhǎng)的湖泊。因此，本文擬運(yùn)用SMDPSO算法對(duì)2009-2018年呼倫湖的藻華進(jìn)行提取，并與光譜指數(shù)FAI進(jìn)行對(duì)比驗(yàn)證。然后，總結(jié)藻華的時(shí)空分布特征，并對(duì)SMDPSO算法的適用性和優(yōu)缺點(diǎn)進(jìn)行評(píng)價(jià)。本研究的成果可以為分析呼倫湖藻華發(fā)生發(fā)展規(guī)律和進(jìn)一步治理和保護(hù)湖泊生態(tài)環(huán)境提供參考。

2 數(shù)據(jù)來(lái)源與研究方法

2.1 研究區(qū)概況

呼倫湖位于內(nèi)蒙古呼倫貝爾市，是內(nèi)蒙古第一大湖泊，其地理位置及水系分布見(jiàn)圖1。湖泊冰封期從11月到次年5月，時(shí)間跨度長(zhǎng)達(dá)6-7個(gè)月[28]。隨著氣候環(huán)境變化和人類活動(dòng)強(qiáng)度的增大，呼倫湖面積逐漸縮減，水體富營(yíng)養(yǎng)化加劇。尤其是自20世紀(jì)80年代以來(lái)，呼倫湖開(kāi)始頻繁出現(xiàn)藻華爆發(fā)的現(xiàn)象，引起了嚴(yán)重的生態(tài)環(huán)境問(wèn)題[4]。呼倫湖藻華的優(yōu)勢(shì)門為藍(lán)藻門，其次是綠藻門[29]。

圖1 呼倫湖地理位置及水系分布

2.2 遙感數(shù)據(jù)及預(yù)處理

基于時(shí)間分辨率、空間分辨率、影像波段數(shù)量和衛(wèi)星發(fā)射時(shí)間等因素的綜合考慮，本文選取2009年至2018年間非冰封期(5月1日-11月1日)的Landsat系列影像對(duì)呼倫湖藻華進(jìn)行遙感識(shí)別。為了避免云的干擾，本文共篩選出53景無(wú)云且質(zhì)量較好的影像進(jìn)行研究。影像預(yù)處理主要包括輻射定標(biāo)和大氣校正。大氣校正通常包括瑞利散射和氣溶膠矯正，但由于呼倫湖水體情況復(fù)雜, 受陸地的影響較大, 氣溶膠變化較為強(qiáng)烈, 因此難以對(duì)呼倫湖進(jìn)行精確的氣溶膠矯正。本研究采用Vanhellemont等[30-31]的算法消除瑞利散射對(duì)遙感影像的影響。另外，為了消除陸地信息對(duì)提取結(jié)果的干擾，進(jìn)一步選用MNDWI(改進(jìn)歸一化差異水體指數(shù))提取呼倫湖水域。

2.3 SMDPSO算法

運(yùn)用基于離散粒子群優(yōu)化的光譜匹配(SMDPSO)算法提取呼倫湖的藻華，需要兩個(gè)步驟：(1)對(duì)遙感影像進(jìn)行光譜匹配獲得藻華概率圖；(2)運(yùn)用DPSO(離散粒子優(yōu)化)算法識(shí)別概率圖中的藻華[27]。

本研究選取呼倫湖中藻華的特征波段作為標(biāo)準(zhǔn)光譜進(jìn)行匹配。標(biāo)準(zhǔn)光譜曲線和目標(biāo)光譜曲線兩者之間的相似程度為藻華的概率，兩者差異越小，則像元被識(shí)別為藻華的可能性就會(huì)越高。計(jì)算方式如下：

PB=cos(B,O)·dist(B,O)

(1)

(2)

(3)

式中：PB為像元是藻華的概率；變量B=(b1,b2,…,bc)和O=(o1,o2,…,oc)分別為藻華的標(biāo)準(zhǔn)光譜矢量和任一像素的光譜矢量；c為波段數(shù)量。

為擴(kuò)大不同波段之間的相對(duì)差異，將藻華的概率值進(jìn)行歸一化處理。概率值越大表明該像素為藻華的可能性就越大，而像元是非藻華的概率值可以用公式(4)表示。

PNB=1-PB

(4)

然后將光譜匹配后得到的概率圖分割為行×列(即rows×cols)的多個(gè)圖塊，將每個(gè)圖塊當(dāng)做自變量代入自適應(yīng)函數(shù)T中，通過(guò)多次迭代，DPSO算法可以求出各圖塊的最大自適應(yīng)值 maxT，而該值所對(duì)應(yīng)的各像素的位置即為分類結(jié)果，最后將各個(gè)圖塊的分類結(jié)果合并即可。具體公式如下：

(5)

(6)

如果將概率圖分為n個(gè)圖塊，則解析過(guò)程重復(fù)n次，而此處只考慮1個(gè)圖塊。本研究設(shè)定1個(gè)圖塊的大小為D=4×4，其適宜解{x1,x2,…,xd}由公式(7)決定。

(7)

(8)

(9)

(10)

(11)

(12)

DPSO算法運(yùn)用于單個(gè)圖塊的具體步驟如下：

(1)隨機(jī)生成I個(gè)D維粒子，即1個(gè)I行D列的矩陣X。矩陣X中的值為0或1。同時(shí)，隨機(jī)生成速度矩陣V，V的取值范圍為0～1。

(2)對(duì)于X的每一行，圖塊中的像元被分為兩類：藻華和非藻華。根據(jù)公式(5)計(jì)算X每行的自適應(yīng)值T，記錄每行的歷史最優(yōu)解，X中最大自適應(yīng)值Tmax所在的那行為全局最優(yōu)解。

(3)運(yùn)用公式(8)、(9)、(11)更新速度矩陣V，運(yùn)用公式(8)～(10)和公式(12)更新位置矩陣X。

(4)再次計(jì)算X每行的自適應(yīng)值T，并比較本次與上次迭代時(shí)X每行的自適應(yīng)值T，選取較大T值所在位置更新歷史最優(yōu)解和全局最優(yōu)解。

(5)當(dāng)k=kmax時(shí)，進(jìn)行下一步，否則返回第(3)步。

(6)全局最優(yōu)解是該圖塊的最終分類結(jié)果。之后依次對(duì)其他圖塊重復(fù)上述步驟。

2.4 光譜指數(shù)法

浮游藻類指數(shù)FAI(floating algae index)是由Hu[32]提出的一種監(jiān)測(cè)海洋與近海岸帶藻類的方法，該指數(shù)對(duì)環(huán)境條件的適應(yīng)能力強(qiáng)，可以有效地提取漂浮藻類，已被廣泛應(yīng)用于藻類的業(yè)務(wù)化監(jiān)測(cè)中。因此，本研究選取FAI對(duì)呼倫湖藻華進(jìn)行遙感識(shí)別，并將結(jié)果作為驗(yàn)證數(shù)據(jù)來(lái)檢驗(yàn)SMDPSO算法的精度與可靠性。

2.5 精度驗(yàn)證方法

以FAI的分類結(jié)果為標(biāo)準(zhǔn)，運(yùn)用格網(wǎng)統(tǒng)計(jì)法對(duì)SMDPSO算法進(jìn)行精度驗(yàn)證。具體原理為：在呼倫湖內(nèi)創(chuàng)建若干個(gè)大小為3 km×3 km的正方形格網(wǎng)，以格網(wǎng)為單位分別統(tǒng)計(jì)SMDPSO算法和FAI提取的藻華面積，之后構(gòu)建散點(diǎn)圖，并計(jì)算兩者之間的決定系數(shù)(R2)和均方根誤差(RMSE)。

3 結(jié)果與分析

3.1 標(biāo)準(zhǔn)光譜曲線與概率圖

首先，對(duì)標(biāo)準(zhǔn)假彩色合成的遙感影像進(jìn)行交互式色彩拉伸，基于先驗(yàn)知識(shí)目視選取藻華與水體的樣點(diǎn)各164個(gè)。然后，計(jì)算樣點(diǎn)在各個(gè)波段的光譜反射率的平均值，得到呼倫湖藻華和水體的標(biāo)準(zhǔn)光譜曲線，如圖2所示。

圖2 呼倫湖藻華和水體的標(biāo)準(zhǔn)光譜曲線

由圖2可以看出，藻華的光譜特征與典型植被相似，且在各個(gè)波段的反射率均大于水體。雖然藻華和水體在綠光波段都具有1個(gè)小的反射峰, 但水體的反射率在綠光波段之后就隨波長(zhǎng)增加而逐漸降低，直至短波紅外波段便趨于0，曲線整體呈現(xiàn)出1個(gè)上下翻轉(zhuǎn)的“對(duì)勾”形狀。與之不同的是，藻華在紅光波段具有明顯的吸收谷，隨后在紅光波段到近紅外波段之間形成一個(gè)“陡坡”，反射率急劇上升,在近紅外波段達(dá)到頂峰。

基于上述特征，本研究以藻華標(biāo)準(zhǔn)光譜曲線中的紅光波段、近紅外波段和短波紅外波段進(jìn)行光譜匹配，并得到藻華概率圖如圖3所示。

圖3 呼倫湖藻華概率圖(2017-08-20)

3.2 分類結(jié)果及驗(yàn)證

以MATLAB R2018軟件為平臺(tái)，運(yùn)用SMDPSO算法提取呼倫湖藻華。其中，輸入對(duì)象為L(zhǎng)andsat影像圖，輸出結(jié)果為二值圖像。之后,采用ArcGIS10.2軟件整理出圖，最終得到2009-2018年呼倫湖藻華分類圖，如圖4所示。由于2012和2016年的大部分影像有云覆蓋，且其余影像也未發(fā)現(xiàn)有藻華爆發(fā)，因此圖4缺少這兩年的分類圖。

圖4 基于SMDPSO算法的2009-2018年呼倫湖藻華分類圖

本文又選取了2015年7月28日和2017年8月2日的遙感影像，對(duì)比SMDPSO算法和FAI的分類結(jié)果，并結(jié)合格網(wǎng)統(tǒng)計(jì)法繪制散點(diǎn)圖(見(jiàn)圖5)，對(duì)結(jié)果進(jìn)行定量評(píng)價(jià)。

由圖5可以看出，SMDPSO算法和FAI的藻華提取結(jié)果十分相似，無(wú)論是面狀、塊狀，還是絮狀的藻華都能被準(zhǔn)確提取。散點(diǎn)圖中絕大多數(shù)點(diǎn)也緊密分布于1∶1線附近，整體擬合度較好。不同分類方法之間的決定系數(shù)R2在2015年7月28日和2017年8月2日分別為0.98和0.97，均方根誤差RMSE均小于0.22 km2，這表明SMDPSO算法可以準(zhǔn)確識(shí)別呼倫湖藻華。

圖5 基于SMDPSO算法和FAI的呼倫湖藻華分類結(jié)果及散點(diǎn)圖(以2015-07-28和2017-08-02遙感影像為例)

3.3 呼倫湖藻華時(shí)空分布特征

通過(guò)計(jì)算各影像中藻華的面積，得到了2009-2018年呼倫湖的藻華面積統(tǒng)計(jì)圖，如圖6所示。

由圖6可知，呼倫湖藻華大多爆發(fā)于每年7和8月，其次是6和9月，藻華面積較大的日期集中在7月下旬至8月中旬(2009年8月12日：125 km2；2011年8月9日：115.5 km2；2015年7月28日: 87.6 km2)。

圖6 2009-2018年呼倫湖藻華面積統(tǒng)計(jì)

為了探究呼倫湖藻華的空間分布特征，本文以像元為單位統(tǒng)計(jì)了2009-2018年藻華出現(xiàn)的次數(shù)，結(jié)果見(jiàn)圖7。

由圖7可知，藻華主要出現(xiàn)在湖泊南部、西南部和北部的邊緣，而湖區(qū)中心卻很少出現(xiàn)藻華，水體情況較好，這可能是周圍土地類型和水流流速所導(dǎo)致的。一方面，湖區(qū)東北部為呼倫湖的旅游景點(diǎn)，人類的頻繁活動(dòng)對(duì)湖水產(chǎn)生了巨大影響；另一方面，由于湖區(qū)西南部距克魯倫河入湖口較近，而此處長(zhǎng)期受到放牧和工廠排放的影響，所以導(dǎo)致河口富營(yíng)養(yǎng)化加劇，很容易出現(xiàn)藻華。

圖7 2009-2018年呼倫湖藻華次數(shù)空間分布

3.4 算法適用性評(píng)價(jià)

為了進(jìn)一步驗(yàn)證SMDPSO算法在其他區(qū)域的適用性，以2016年7月2日黃海區(qū)域的OLI影像為數(shù)據(jù)源對(duì)滸苔(屬綠藻門，呈綠色)進(jìn)行識(shí)別，并與FAI的分類結(jié)果進(jìn)行比較，見(jiàn)圖8。其中圖8(b)和8(c)分別為SMDPSO和FAI方法的分類結(jié)果，其中綠色表示滸苔，藍(lán)色表示海水。

圖8 不同方法對(duì)黃海滸苔的提取結(jié)果及散點(diǎn)圖(2016-07-02)

圖8中的對(duì)比結(jié)果表明，滸苔的空間分布具有高度一致性，SMDPSO算法能夠很好地捕捉和描繪滸苔條紋的細(xì)節(jié)與形狀，即使是很小的斑塊也能清晰地識(shí)別出來(lái)。此外，結(jié)合格網(wǎng)統(tǒng)計(jì)法繪制的SMDPSO和FAI分類結(jié)果的散點(diǎn)圖也表現(xiàn)出非常好的線性關(guān)系，擬合線幾乎與1∶1線重合，R2達(dá)到0.97，RMSE為0.015 km2，這說(shuō)明SMDPSO算法可以有效地識(shí)別黃海的滸苔。同時(shí)，也為算法的進(jìn)一步發(fā)展和應(yīng)用提供了方向。

與光譜指數(shù)法相比，SMDPSO算法既保留了其分類精度高的特點(diǎn)，也在其他方面具有獨(dú)特的優(yōu)勢(shì)。首先，SMDPSO算法的參數(shù)少、成本低，且無(wú)需人工干預(yù)。此算法所涉及的參數(shù)僅包括標(biāo)準(zhǔn)光譜和圖塊大小，而且一旦確定就無(wú)需改變。但是，使用光譜指數(shù)法進(jìn)行遙感分類時(shí)，經(jīng)常需要運(yùn)用目視解譯方法對(duì)每一幅影像的閾值進(jìn)行調(diào)整以保證分類精度，這種方式不僅具有很大的主觀性，而且也會(huì)消耗許多時(shí)間和人力。所以，SMDPSO算法比光譜指數(shù)法更加簡(jiǎn)單方便。其次，由于SMDPSO算法是基于目標(biāo)地物的標(biāo)準(zhǔn)光譜，通過(guò)計(jì)算像元之間的相似性來(lái)識(shí)別地物，所以如果標(biāo)準(zhǔn)光譜發(fā)生變化，SMDPSO算法同樣可以識(shí)別不同種類的藻華或其他類別的地物，這是單一的光譜指數(shù)不能實(shí)現(xiàn)的。然而，SMDPSO算法也存在一些缺點(diǎn)，即在運(yùn)行過(guò)程中消耗的時(shí)間較長(zhǎng)，所以之后的研究重點(diǎn)在于如何在保證精度的前提下提升算法的效率。

4 結(jié) 論

本文利用基于離散粒子群優(yōu)化的光譜匹配(SMDPSO)算法對(duì)2009-2018年呼倫湖的藻華進(jìn)行識(shí)別，并用浮游藻類指數(shù)(FAI)的分類結(jié)果進(jìn)行檢驗(yàn)。主要結(jié)論如下：

(1)SMDPSO算法對(duì)呼倫湖藻華的提取效果較好，與FAI高度一致。

(2)呼倫湖藻華的高發(fā)期為7-8月，且主要爆發(fā)于湖泊邊緣，在湖中心出現(xiàn)的次數(shù)較少。

(3)SMDPSO算法的分類精度與FAI相似，但其成本更低、所需參數(shù)較少，且無(wú)需人工干預(yù)，甚至可以通過(guò)改變標(biāo)準(zhǔn)光譜來(lái)識(shí)別其他地物(如識(shí)別黃海的滸苔)；

(4)在之后的研究中，還需要進(jìn)一步驗(yàn)證SMDPSO算法在多種水生植物并存環(huán)境下的適用性，并將其推廣應(yīng)用于其他具有高分辨率的多光譜遙感影像中。