基于國(guó)產(chǎn)高分辨率光學(xué)遙感影像的水體提取

鄧富亮，章欣欣，花利忠，李宗梅

(廈門理工學(xué)院空間信息技術(shù)研究所，福建 廈門 361024)

基于國(guó)產(chǎn)高分辨率光學(xué)遙感影像的水體提取

鄧富亮，章欣欣，花利忠，李宗梅

(廈門理工學(xué)院空間信息技術(shù)研究所，福建 廈門 361024)

遙感圖像中地表水體同山體、建筑物等地物產(chǎn)生的陰影在光譜特征上存在較高的類間相似性，導(dǎo)致提取過程中容易出現(xiàn)混淆和錯(cuò)分的情況。針對(duì)此問題，提出一種基于面向?qū)ο蠛腿斯し淙旱牡乇硭w提取方法。該方法首先對(duì)遙感圖像進(jìn)行分割以獲取分割對(duì)象的光譜、比率、幾何形狀等統(tǒng)計(jì)特征，以彌補(bǔ)高分遙感圖像波段數(shù)目少，信息量不足的缺陷；并借助人工蜂群算法在解決復(fù)雜問題最優(yōu)化方面的優(yōu)勢(shì)，選取水體同陰影二值分類的幾何平均正確率作為算法的適應(yīng)度函數(shù)，最終獲取地表水體的最優(yōu)化提取規(guī)則。選取廈門市大嶝島和湖南省資興市部分區(qū)域，基于國(guó)產(chǎn)高分一號(hào)、二號(hào)遙感數(shù)據(jù)進(jìn)行水體提取，并與傳統(tǒng)SVM分類結(jié)果進(jìn)行比較。實(shí)驗(yàn)結(jié)果表明本算法提取水體的總體精度和Kappa系數(shù)均優(yōu)于傳統(tǒng)SVM分類器，表明該方法可應(yīng)用于高分遙感圖像的地表水體提取。

人工蜂群；地表水體提??；高分；面向?qū)ο?；遙感圖像分類

水體的空間分布信息是水文模型分析的數(shù)據(jù)基礎(chǔ)，如何快速、準(zhǔn)確地獲取水體信息對(duì)地表水資源調(diào)查[1]、水利設(shè)施規(guī)劃、洪澇[2]等自然災(zāi)害監(jiān)測(cè)等領(lǐng)域具有重要意義。隨著我國(guó)高分辨率對(duì)地觀測(cè)系統(tǒng)的可持續(xù)發(fā)展，以及國(guó)產(chǎn)高空間分辨率衛(wèi)星，如高分、資源遙感衛(wèi)星遙感數(shù)據(jù)的普及，針對(duì)國(guó)產(chǎn)高分遙感影像的水體信息的智能提取和識(shí)別成為當(dāng)前國(guó)內(nèi)研究的熱點(diǎn)。

1 水體提取模型研究進(jìn)展

目前國(guó)內(nèi)外學(xué)者針對(duì)不同分辨率影像的水體提取提出諸多方法和模型，這些方法從構(gòu)建原理大致可分為四種類型[3]：?jiǎn)尾ǘ伍撝捣指睢⒍嗖ǘ伪戎的Ｐ?、基于分類算法的提取模型以及基于混合像元線性分解(Linear unmixing)模型。

1.1 基于單/多波段閾值分割

單一波段的閾值分割模型最為簡(jiǎn)單，易于實(shí)現(xiàn)，但受限于單個(gè)波段信息豐度的制約，其提取精度較低；多波段比值關(guān)系方法綜合利用了地物在不同波段光譜的特征，通過設(shè)計(jì)比值計(jì)算方法進(jìn)一步削弱了傳感器成像條件、大氣輻射等影響導(dǎo)致的誤差，提高了分類的準(zhǔn)確性[4]。例如，在中低分辨率影像的水體提取方面，目前已被廣泛使用的有NDWI[5]、MNDWI[6]、SWI[7]、AWEI[3]等各種指數(shù)模型，均取得很好的效果[8]。但此類方法的分割閾值的判定需要結(jié)合專家知識(shí)和具體實(shí)際情況，經(jīng)過多次嘗試從而獲得最佳的結(jié)果，在不同研究區(qū)中存在較大差異，無法保證其方法的普適性[9]。此外，上述大部分方法需要借助中紅外等光譜波段特征，無法應(yīng)用于僅有四個(gè)波段的高分遙感影像。

1.2 基于分類器模型

基于分類算法的城市地物提取，包括水體，目前也取得一系列成果。但由于高分圖像中“同譜異物”現(xiàn)象明顯，如水體和建筑物陰影兩種地物類別四個(gè)波段的灰度值近似，很容易造成混淆，因此借助各種分類器的水體特征算法更為有效，常見的分類器模型如決策樹[10]、最大似然(Maximum Likelihood，ML)[11]、支持向量機(jī)(Support Vector Machine, SVM)[12]等，但此類算法在水體提取應(yīng)用方面也同樣存在一定限制，如最大似然方法需要預(yù)先設(shè)定實(shí)際地物類別，對(duì)分類樣本的選取要符合統(tǒng)計(jì)上的獨(dú)立分布規(guī)律，其分類結(jié)果和提取的特征參數(shù)基于全局?jǐn)M合精度，因而無法直接使用于單個(gè)指定的地物類別；SVM算法將特征映射到高緯空間，利用核函數(shù)尋找最優(yōu)支持向量進(jìn)行分類，其本質(zhì)上屬于“黑箱”算法，難以對(duì)分類規(guī)則進(jìn)行直觀解釋。

1.3 基于混合像元分解

混合像元分解模型基于光譜成像機(jī)理，模擬光譜的混合過程，并通過改進(jìn)分解模型以提高分類精度。但在實(shí)際應(yīng)用中，由于高分辨率光譜單個(gè)像元包含的信息量較少、光譜信息不穩(wěn)定，與端元樣本的信息對(duì)比困難，往往無法同分解模型較好地吻合。此外，高分辨率的多光譜數(shù)據(jù)對(duì)于純凈端元的提取需要基于研究區(qū)實(shí)際地物情況和具體的成像參數(shù)，針對(duì)水體邊界的水分含量的估算需要做進(jìn)一步的優(yōu)化，實(shí)現(xiàn)精細(xì)提取和識(shí)別仍需深入研究。

此外，傳統(tǒng)基于像素的分類方法應(yīng)用于高分圖像上容易引起“椒鹽效應(yīng)”，一些學(xué)者通過引入基于面向?qū)ο蟮膱D像分割技術(shù)，將圖像劃分成統(tǒng)計(jì)特性相近的對(duì)象塊，一方面可對(duì)特征象元特征進(jìn)行聚類分割，使得更高層影像的分析和理解成為可能；另一方面，可針對(duì)對(duì)象幾何、光譜、紋理的特征進(jìn)行二次提取，增加分割對(duì)象屬性的特征維度，可有效提高提取結(jié)果的精度[13]，成為提高高分遙感影像分類精度的關(guān)鍵。

基于高空間分辨率遙感數(shù)據(jù)的水體信息提取模型，其優(yōu)勢(shì)通常體現(xiàn)在以下兩個(gè)方面：一是借助空間分辨率的提高，模型除辨識(shí)出一般的大面積水體之外，還能有效地識(shí)別圖像中存在的細(xì)小水體。細(xì)小水體的空間特征一般表現(xiàn)為細(xì)小狹長(zhǎng)狀，且易受周邊地物的光譜影響，在水體邊緣表現(xiàn)為混合像元；二是改進(jìn)后的提取或分類模型能更為有效地去除暗區(qū)的干擾信息，如建筑物陰影、山區(qū)陰影等。針對(duì)以上問題，本文提出一種基于面向?qū)ο蠓指詈腿斯し淙核惴ǖ某鞘兴w提取方法。該方法首先基于面向?qū)ο笏枷雽?duì)高分遙感圖像進(jìn)行多尺度的圖像分割，并統(tǒng)計(jì)不同地物的特征信息，包括單波段灰度值、多波段比值、幾何特征等。隨后利用人工蜂群算法在解決復(fù)雜系統(tǒng)問題最優(yōu)化的優(yōu)勢(shì)，選取合理的適應(yīng)度函數(shù)，進(jìn)行最優(yōu)規(guī)則的提取。最后以廈門大嶝島和湖南省資興市部分區(qū)域的高分衛(wèi)星為研究區(qū)驗(yàn)證模型的分類精度。

2 面向?qū)ο蟮娜斯し淙核w提取方法

本文提出的水體信息提取方法包括兩部分內(nèi)容：(1)采用基于面向?qū)ο蟮膱D像分割方法對(duì)高分?jǐn)?shù)據(jù)進(jìn)行多尺度的對(duì)象分割，根據(jù)像元特征的相似性分割成不同的對(duì)象。同時(shí)對(duì)各個(gè)對(duì)象的光譜、比率、幾何形狀等特征進(jìn)行二次統(tǒng)計(jì)，提高特征信息的維度；(2)提出一種基于人工蜂群算法的水體特征提取算法，針對(duì)水體提取研究中，水體同山區(qū)、建筑物陰影容易混淆的現(xiàn)象，以水體、陰影二者分類的幾何平均正確率作為算法的適應(yīng)度函數(shù)，并采用監(jiān)督分類方式構(gòu)建提取模型，提取規(guī)則基于模糊規(guī)則系統(tǒng)(Fuzzy rule-based system)，可有效辨識(shí)不同特征對(duì)地物的影響程度。模型流程圖如圖 1所示。

圖1 模型流程圖Fig.1 Flowchart of the proposed model

2.1 面向?qū)ο蟮母叻诌b感圖像分割

目前，基于面向?qū)ο蟮膱D像分割算法很多，但大多算法未能整合光譜信息和空間信息，且不同算法在效率和結(jié)果上存在較大差異，缺乏統(tǒng)一的評(píng)價(jià)體系；另一方面，不同區(qū)域空間特征復(fù)雜多樣，理論上圖像分割算法的最優(yōu)尺度無法得到，常見做法是通過專家知識(shí)進(jìn)行人工判定。本文采用的圖像分割算法步驟如下：首先采用快速掃描法對(duì)原始遙感圖像進(jìn)行初始分割，形成分割初始對(duì)象，并構(gòu)建區(qū)域鄰接圖；隨后采用基于異質(zhì)性最小區(qū)域合并算法進(jìn)行區(qū)域合并；最后根據(jù)某種層次輸出規(guī)則，在區(qū)域合并過程中輸出滿足層次輸出條件的中間分割結(jié)果，直至滿足分割終止條件。

結(jié)合高分辨率遙感影像和水體特征，除原始四個(gè)波段外，借助ArcGIS波段計(jì)算工具，生成第五(Blue/Green)、第六(Green/Red)、第七(Red/NIR)、第八(NDVI)四個(gè)比率波段，以提高圖像分割算法和后續(xù)分類器的分類精度。

針對(duì)分割對(duì)象進(jìn)行二次統(tǒng)計(jì)，可獲取每個(gè)對(duì)象的光譜、比率、幾何形狀等特征信息。并通過相關(guān)性檢驗(yàn)法對(duì)因子進(jìn)行篩選。

2.2 基于人工蜂群的水體提取模型

人工蜂群算法(Artificial Bee Colony， ABC)是群智能算法中的一種新興技術(shù)[14]，其原理同螞蟻、遺傳算法類似，通過模擬個(gè)體間的相互協(xié)作和競(jìng)爭(zhēng)關(guān)系，進(jìn)而反映出種群的群體性智能行為，因此適合于解決一些非線性的復(fù)雜問題[15]。ABC蜂群算法的優(yōu)勢(shì)在于其在迭代過程可同時(shí)進(jìn)行全局和局部搜索，具有收斂快、適用性強(qiáng)的優(yōu)點(diǎn)。且算法本身靈活易于實(shí)現(xiàn)，可同其它技術(shù)結(jié)合以改進(jìn)原算法。已被廣泛應(yīng)用到如復(fù)雜問題優(yōu)化[16]、圖像變化檢測(cè)[17]、圖像分割和分類[18]等多方面，取得較好的成果。

泛用的ABC模型應(yīng)至少包括四方面內(nèi)容，即優(yōu)化規(guī)則定義、適應(yīng)度函數(shù)、全局/周邊鄰域搜索步驟以及模型預(yù)測(cè)。針對(duì)遙感影像水體提取的步驟，本文對(duì)ABC模型進(jìn)行如下改進(jìn)：

水體提取規(guī)則的定義：采用了基于模糊規(guī)則(Fuzzy rule-based)的IF…THEN形式。對(duì)于分割閾值選取可采用最大最小值的閾值區(qū)間表示。相比傳統(tǒng)的“硬”分類器，如最大似然、SVM方法，基于模糊規(guī)則的分類體系無需依賴具體的數(shù)學(xué)公式，所提取的規(guī)則更為直觀、易于解釋。

適應(yīng)度函數(shù)選?。哼m應(yīng)度函數(shù)的選取是ABC算法的核心，決定了ABC算法優(yōu)化的最終目標(biāo)。對(duì)于水體提取而言，影響水體分類精度的主要原因是水體同陰影二者在圖像上均呈現(xiàn)暗色調(diào)，即相同的灰度值從而引起的混淆。因此，針對(duì)水體提取精度的優(yōu)化可以轉(zhuǎn)換為如何最大化地區(qū)分水體同城市建筑物陰影兩種地物類別。實(shí)際情況中，由于水體同建筑物陰影在數(shù)量級(jí)存在較大差異，建筑物陰影對(duì)象相對(duì)較為破碎，數(shù)量較多。因此，本文引入幾何平均正確率PG-mean指數(shù)作為ABC算法的適應(yīng)度函數(shù)，其公式如下：

(1)

表1 二值分類結(jié)果統(tǒng)計(jì)指標(biāo)Table 1 Statistic values of the binary classification

全局/周邊鄰域查找：在蜂群算法的初始化構(gòu)建過程采用隨機(jī)初始化，全局搜索一般采用隨機(jī)搜索策略，便于跳出局部最優(yōu)解。而周邊鄰域查找則可使算法逐步收斂，是提高適應(yīng)度評(píng)價(jià)的關(guān)鍵步驟。周邊鄰域查找可采用如下公式表達(dá)：

(2)

2.3 提取結(jié)果和精度驗(yàn)證

根據(jù)ABC算法得到的最終分類規(guī)則可用于水體的類別提取。由于基于模糊規(guī)則得到的提取規(guī)則普遍規(guī)則類別重疊的現(xiàn)象，即一個(gè)樣本數(shù)據(jù)可能滿足多個(gè)分類規(guī)則。因此，在本文算法中額外引入一個(gè)指標(biāo)，即規(guī)則的樣本覆蓋度Pcover，該指標(biāo)反映了當(dāng)前規(guī)則分類正確的樣本數(shù)占所屬類別樣本總數(shù)的比例。當(dāng)出現(xiàn)樣本重疊時(shí)，通過綜合比較適應(yīng)度和覆蓋度確定最終預(yù)測(cè)結(jié)果，計(jì)算公式為：

P=α×PG-mean+(1-α)×Pcover

(3)

其中，α為權(quán)重閾值，N為樣本總數(shù)。α的取值在[0,1]之間，其值反映了兩個(gè)指標(biāo)對(duì)預(yù)測(cè)結(jié)果的影響關(guān)系。為了得到最優(yōu)結(jié)果，可以對(duì)α進(jìn)行迭代計(jì)算統(tǒng)計(jì)，本文以0.1為步長(zhǎng)，利用公式對(duì)樣本數(shù)據(jù)進(jìn)行循環(huán)計(jì)算，統(tǒng)計(jì)P值以及閾值α，最終得到提取結(jié)果。

為進(jìn)一步對(duì)模型提取結(jié)果的精度進(jìn)行驗(yàn)證和評(píng)價(jià)，實(shí)驗(yàn)采用監(jiān)督分類方法的傳統(tǒng)評(píng)價(jià)方法，即分類精度和Kappa系數(shù)，對(duì)分類結(jié)果進(jìn)行驗(yàn)證。同時(shí)結(jié)合可視化分類結(jié)果進(jìn)行目視判讀分析，并對(duì)提取的規(guī)則進(jìn)行分析，綜合評(píng)價(jià)ABC算法應(yīng)用于水體提取的有效性。

3 實(shí)驗(yàn)

3.1 數(shù)據(jù)預(yù)處理

本文選取廈門市大嶝島的GF-1號(hào)和湖南省資興市的GF-2號(hào)部分?jǐn)?shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù)。GF-1作為中國(guó)首個(gè)高空間分辨率對(duì)地觀測(cè)衛(wèi)星，于2013年4月26號(hào)發(fā)射，可提供8 m分辨率多光譜和2 m分辨率全色波段影像。而GF-2的星下點(diǎn)像元分辨率則進(jìn)一步提高到0.8 m全色和3.2 m多光譜波段。借助ENVI遙感處理軟件，原始數(shù)據(jù)經(jīng)過正射校正、幾何校正、輻射定標(biāo)、Flaash大氣校正、全色波段融合等預(yù)處理，獲得四個(gè)多光譜波段，包括藍(lán)(0.45～0.52 μm)、綠(0.52～0.59 μm)、紅(0.63～0.69 μm)以及近紅外(0.77～0.89 μm)。

隨后根據(jù)圖像中水系的空間分布特征進(jìn)行裁剪，最終獲得研究區(qū)圖像(圖2)，其中大嶝島數(shù)據(jù)行列大小為1 276×1 103像素，資興市為1 866×1 592像素。兩個(gè)實(shí)驗(yàn)區(qū)水體的空間典型特征主要體現(xiàn)在：廈門大嶝島數(shù)據(jù)中水體主要包括沿海海域、養(yǎng)殖、農(nóng)田水塘，不同水體面積差別較大；而資興市水體則包括河流、池塘、農(nóng)田灌溉引導(dǎo)等細(xì)小水體。

圖2 實(shí)驗(yàn)區(qū)假彩色合成影像(RGB波段：432)Fig.2 False color composite images of the study areas

由于本文主要研究水體信息提取，因此圖像裁剪完畢后，首先統(tǒng)計(jì)各個(gè)象元的亮度值，隨后采用閾值分割法進(jìn)行暗區(qū)掩膜處理，排除如建設(shè)用地、農(nóng)田、未利用地等地類像元。隨后借助ENVI的ROI樣本提取工具，采用人工判讀法對(duì)掩膜后的圖像進(jìn)行隨機(jī)采樣，所選擇土地利用類型包括水體、建設(shè)用地、農(nóng)田、未利用地、道路以及陰影，其中，非水體地物類型在分類完成后通過重采樣合為一類。兩幅遙感影像共隨機(jī)提取約450個(gè)矩形框樣本，其中大嶝島提取1 213/58個(gè)樣本點(diǎn)/矩形用于訓(xùn)練，11 614/98個(gè)樣本點(diǎn)/矩形用于驗(yàn)證；資興市則提取了2 743/83個(gè)樣本點(diǎn)/矩形用于訓(xùn)練，22 176/223個(gè)樣本點(diǎn)/矩形用于驗(yàn)證。

3.2 多尺度圖像分割結(jié)果

將遙感影像的四個(gè)原始波段、四個(gè)比率波段，共八個(gè)波段作為輸入特征，利用多尺度圖像分割算法獲得分割結(jié)果，尺度選取方面則設(shè)置起始尺度70，終止尺度250，總計(jì)輸出8個(gè)尺度的分割結(jié)果。以廈門市為例，不同尺度的對(duì)象分割數(shù)目如表2所示。

表2 多尺度分割圖像的統(tǒng)計(jì)結(jié)果Table 2 Statistic results of the multi-scale segmentation of imagery

根據(jù)實(shí)驗(yàn)區(qū)的地物空間分布特征，結(jié)合專家知識(shí)判斷可發(fā)現(xiàn)，在70尺度下，如圖3(a)及其局部放大圖3(c)所示，分割對(duì)象最多，對(duì)應(yīng)的小型建筑物、小面積水體(如農(nóng)田周邊的池塘)、建筑物周邊綠化設(shè)施的輪廓分割的較為清晰，但道路、林地、草地、大面積水體(海域)等則較為破碎；而在100尺度下，如圖3(b)及其局部放大圖3(d)所示，道路、大面積水體的輪廓較為符合實(shí)際，但細(xì)小水體容易同周邊地物合并；后續(xù)的尺度則普遍存在不同地物類別被合并為同一對(duì)象的情況，如用于分類會(huì)降低最后的分類精度。因此確定尺度70的分割圖像用于識(shí)別。

3.3 圖像對(duì)象特征統(tǒng)計(jì)及分類

針對(duì)分割結(jié)果選取不同的統(tǒng)計(jì)指標(biāo)，可有效地增加分割對(duì)象的信息維度，因此，實(shí)驗(yàn)對(duì)前四個(gè)單一光譜波段分別計(jì)算圖像對(duì)象的光譜亮度、標(biāo)準(zhǔn)差、最大/小像素值、最大最小像素值比率四個(gè)光譜特征；后四個(gè)比值波段則計(jì)算器圖像對(duì)象的光譜均值、光譜值標(biāo)準(zhǔn)差兩個(gè)特征；以及對(duì)象的幾何特征，包括光滑度、粗糙度、正圓比、長(zhǎng)寬比四個(gè)形狀特征因子。所得的對(duì)象樣本特征值，利用SPSS統(tǒng)計(jì)軟件進(jìn)行相關(guān)性檢驗(yàn)，排除相關(guān)性系數(shù)大于0.7的特征因子。通過特征篩選，最后共選取了包括像素值、波段均值、標(biāo)準(zhǔn)差、多光譜比值、形狀等共13個(gè)特征，具體描述見表3。

表3 分割對(duì)象的特征提取信息Table 3 Information of the segmentation object

以上特征因子作為ABC算法的輸入特征要素。ABC算法參數(shù)設(shè)置如下：蜂群規(guī)模200，采蜜蜂100，搜尋次數(shù)200，循環(huán)次數(shù)600。輸入樣本數(shù)據(jù)執(zhí)行算法訓(xùn)練模型并提取分類規(guī)則，隨后利用驗(yàn)證樣本進(jìn)行模型精度評(píng)價(jià)，驗(yàn)證結(jié)果以混淆矩陣形式描述，并同傳統(tǒng)SVM分類器分類結(jié)果進(jìn)行比較，具體結(jié)果見表 4。

圖3 多尺度分割結(jié)果Fig.3 Results of the multi-scale segmentation

SVMABC地物類別水體陰影未分類產(chǎn)品精度/%水體陰影未分類產(chǎn)品精度/%水體4975139028778253327942678340陰影1304545546870518248521879293用戶精度/%9745765896708594總體精度/%81978769Kappa系數(shù)0658807633

由表 4可知，利用ABC算法構(gòu)建的水體提取模型，其驗(yàn)證樣本的分類精度為87.69%，Kappa系數(shù)為0.7633，其中水體提取生產(chǎn)精度為83.40%，非水體(主要為建筑物陰影)精度為92.93%；而用戶精度則分別為96.70%和85.94%。而采用傳統(tǒng)基于像素的SVM分類器，其總體分類精度為81.97%，Kappa系數(shù)為0.6588，其中水體地類的產(chǎn)品精度為77.82%，陰影為87.05%，對(duì)應(yīng)的用戶精度則為97.45%和76.58%。兩種分類器精度對(duì)比可知，基于SVM相比ABC可識(shí)別更多的水體樣本，但其中有1 390個(gè)水體樣本被錯(cuò)分為陰影，多于ABC的794樣本。同時(shí)，ABC分類器可識(shí)別出4 852個(gè)陰影樣本，同樣高于SVM分類器的4 545個(gè)樣本?？傮w而言，基于ABC的分類器在高分影像水體提取的精度方面，要優(yōu)于傳統(tǒng)SVM分類器，表明本文所提出的ABC分類器能更好地滿足實(shí)際應(yīng)用需求。

3.4 水體提取結(jié)果

基于ABC分類模型進(jìn)一步提取大嶝島的地表水體分布圖，結(jié)果如圖 4所示。采用目視判讀法可知，基于此模型由于預(yù)先進(jìn)行了圖像分割，因此水體對(duì)象的幾何形狀更為完整，此外，二次統(tǒng)計(jì)引入的對(duì)象幾何特征進(jìn)一步彌補(bǔ)了原始的光譜波段信息，為分類器提供了更多的特征屬性，提高了分類器對(duì)地表水體和建筑物陰影的辨識(shí)能力，進(jìn)一步提高了分類精度。

圖4 大嶝島分類結(jié)果比較Fig.4 Comparison of the extraction results of the Dadeng Island

采用相同模型輸入?yún)?shù)，對(duì)資興市的地表水提進(jìn)行提取，提取結(jié)果如圖5所示。通過目視判讀比較可知，基于面向圖像分割和ABC優(yōu)化的分類算法所提取的細(xì)小水體，形狀更為完整，而采用SVM分類器所提取的水體則表現(xiàn)出破碎、不連續(xù)現(xiàn)象，表明采用圖像分割后的地物對(duì)象有助于提取狹窄的細(xì)小水體。

圖5 ABC和SVM分類器的分類結(jié)果比較Fig.5 Comparison of the extraction results of Zixing

4 結(jié)論

針對(duì)高分遙感影像中水體提取面臨的問題，如細(xì)小水體、易同陰影混淆情況，提出一種基于面向?qū)ο蠓指詈腿斯し淙核惴ǖ乃w提取方法。該方法首先基于面向?qū)ο髮?duì)高分遙感圖像進(jìn)行多尺度圖像分割，引入額外四個(gè)光譜比值波段作為輸入特征，提高圖像的分割精度。并統(tǒng)計(jì)分割圖像的特征信息，包括單波段灰度值、多波段比值、幾何形狀特征等。隨后利用人工蜂群算法在解決復(fù)雜系統(tǒng)問題最優(yōu)化的優(yōu)勢(shì)，選取水體和建筑物陰影二值分類的幾何平均正確率作為適應(yīng)度函數(shù)，進(jìn)行最優(yōu)規(guī)則的提取。實(shí)驗(yàn)基于國(guó)產(chǎn)GF-1、GF-2衛(wèi)星數(shù)據(jù)。提取結(jié)果表明本算法可有效地區(qū)分水體與建筑物陰影。根據(jù)所提取的規(guī)則進(jìn)一步分析可發(fā)現(xiàn)，雖然二者地物在大部分特征屬性上相似，但在Red/NIR的波段均值、分割對(duì)象的亮度和標(biāo)準(zhǔn)差3個(gè)特征表現(xiàn)出一定程度的差異性，可作為辨識(shí)水體同陰影的參考依據(jù)。后續(xù)研究可進(jìn)一步結(jié)合其它指數(shù)模型，如陰影水體指數(shù)SWI，進(jìn)一步擴(kuò)展特征向量的信息維度；同時(shí)可進(jìn)一步優(yōu)化目前的分割算法，如引入GLCM等幾何紋理特征，進(jìn)一步提高圖像分割對(duì)象的精度。

[1] 陶明剛. Landsat-TM遙感影像岸線變遷解譯研究——以九龍江河口地區(qū)為例[J]. 水文地質(zhì)工程地質(zhì), 2006, 33(1):107-110.[TANG M G. Interpretation of coastline changes based on Landset-TM remote sensing images: an example of the mouth area of Jiulongjiang River[J]. Hydrogeology & Engineering Geology, 2006, 33(1): 107-110. (in Chinese)]

[2] 李政國(guó), 張茂省, 李林. 面向?qū)ο蟮倪b感影像承災(zāi)體邊緣提取方法[J]. 水文地質(zhì)工程地質(zhì), 2011, 38(1):128-132.[LI Z G, ZHANG M S, LI L. A study of object-oriented extracting edge of hazard-affected body in remote sensing image[J]. Hydrogeology & Engineering Geology, 2011, 38(1):128-132. (in Chinese)]

[3] FEYISA G L, MEILBY H, FENSHOLT R,etal. Automated water extraction index: a new technique for surface water mapping using landsat imagery[J]. Remote Sensing of Environment, 2014, 140(1):23-35.

[4] 陳文倩, 丁建麗, 李艷華, 等. 基于國(guó)產(chǎn)GF-1遙感影像的水體提取方法[J]. 資源科學(xué), 2015, 37(6):1166-1172.[CHEN W Q, DING J L, LI Y H，etal. Extraction of water information based on China-made GF-1 remote sense image[J]. Resources Science, 2015, 37(6):1166-1172.(in Chinese)]

[5] 周藝, 謝光磊, 王世新, 等. 利用偽歸一化差異水體指數(shù)提取城鎮(zhèn)周邊細(xì)小河流信息[J]. 地球信息科學(xué)學(xué)報(bào), 2014, 16(1):102-107.[ZHOU Y, XIE G L, WANG S X,etal. Information extraction of thin rivers around built-up lands with false NDWI[J]. Journal of Geo-information Science, 2014, 16(1):102-107.(in Chinese)]

[6] 徐涵秋. 利用改進(jìn)的歸一化差異水體指數(shù)(MNDWI)提取水體信息的研究[J]. 遙感學(xué)報(bào), 2005, 9(5):589-595. [XU H Q. A Study on information extraction of water body with the Modified Normalized Difference Water Index(MNDWI)[J]. Journal of Remote Sensing, 2005, 9(5): 589-595.(in Chinese)]

[7] 曹凱, 江南, 李新國(guó), 等. 基于SPOT-5圖像的城市水體自動(dòng)提取模型研究[J]. 國(guó)土資源遙感, 2005, 17(4):24-27.[CAO K, JIANG N, LI X G,etal. A model for automatic extraction of water information form the city zone based on SPOT-5 Image[J]. Remote Sensing for Land & Resources, 2005, 17(4): 24-27.(in Chinese)]

[8] 王志輝, 易善楨. 不同指數(shù)模型法在水體遙感提取中的比較研究[J]. 科學(xué)技術(shù)與工程, 2007, 7(4):534-537.[WANG Z H, YI S Z. Comparison and research on the different index models used in water extraction by remote sensing[J]. Science Technology and Engineering, 2007, 7(4):534-537.(in Chinese)]

[9] LI B, ZHANG H, XU F. Water extraction in high resolution remote sensing image based on hierarchical spectrum and shape features[C]//Proceedings of the IOP Conference Series: Earth and Environmental Science, F, 2014.

[10] 陳靜波, 劉順喜, 汪承義, 等. 基于知識(shí)決策樹的城市水體提取方法研究[J]. 遙感信息, 2013, 28(1):29-33.[CHEN J B, LIU S X, WANG C Y,etal. Research on urban water body extraction using knowledge-based decision tree[J]. Remote Sensing Information, 2013, 28(1):29-33.(in Chinese)]

[11] 賈科利, 常慶瑞. 利用主成分分析法提取水體信息[J]. 淮陰師范學(xué)院學(xué)報(bào)(自然科學(xué)版), 2006, 5(3):242-245.[JIA K L, CHANG Q R. Extracting water body information using principal components analysis[J]. Journal of HuaiYin Teachers College(Natural Science Edition), 2006, 5(3):242-245.(in Chinese)]

[12] 朱長(zhǎng)明, 張新, 駱劍承, 等. 基于樣本自動(dòng)選擇與SVM結(jié)合的海岸線遙感自動(dòng)提取[J]. 國(guó)土資源遙感, 2013, 25(2):69-74.[ZHU C M, ZHANG X, LUO J C,etal. Automatic extraction of coastline by remote sensing technology based on SVM and auto-selection of training samples[J]. Remote Sensing for Land and Resources, 2013, 25(2):69-74.(in Chinese)]

[13] 徐濤, 譚宗坤, 閆小平. 面向?qū)ο蟮某鞘兴w信息提取方法[J]. 地理空間信息, 2010, 8(3):64-66.[XU T, TAN Z K, YAN X P. Extraction techniques of urban water bodies based on object-oriented[J]. Geospatial Information, 2010, 8(3):64-66.(in Chinese)]

[14] KARABOGA D. An idea based on honey bee swarm for numerical optimization[J/OL]. ResearchGATE,2005.[http://www.researchgate.net/publication/][2005-01]

[15] KARABOGA D, AKAY B. A survey: algorithms simulating bee swarm intelligence[J]. Artificial Intelligence Review, 2009, 31(1-4):61-85.

[16] KARABOGA D, AKAY B. A modified artificial bee colony (ABC) algorithm for constrained optimization problems[J]. Applied Soft Computing, 2011, 11(3): 3021-3031.

[17] 肖永豪, 余衛(wèi)宇. 基于蜂群算法的圖像邊緣檢測(cè)[J]. 計(jì)算機(jī)應(yīng)用研究, 2010, 27(7):2 748-2 750.[XIAO Y H, YU W Y. Bee colony algorithm for image edge detection[J]. Application Research of Computers, 2010, 27(7):2748-2750.(in Chinese)]

[18] MARINAKI M, MARINAKIS Y, ZOPOUNIDIS C. Honey bees mating optimization algorithm for financial classification problems[J]. Applied Soft Computing, 2010, 10(3):806-812.

責(zé)任編輯：汪美華

A surface water body extraction method based on domesticremote sensing imagery of high resolution

DENG Fuliang, ZHANG Xinxin, HUA Lizhong, LI Zongmei

(CollegeofComputerandInformationEngineering,XiamenUniversityofTechnology,Xiamen,Fujian361024，China)

Due to the high spectral similarity existing in water and shadow, extraction of remote sensing imagery is easily confused and misclassified. To address this problem, we propose a method combined with the object-oriented image segmentation and the artificial bee colony algorithm (ABC) to extract surface water body from remote sensing imagery. Firstly, a series of statistic factors, such as spectrum, ratio and sharp features, are calculated during image segmentation. We used these factors to make up the defect of insufficient information existing in high-resolution imagery. Then, with the strength of solving complicate problem by ABC algorithms, we chose the geometric mean of accuracies between surface water bodies and shadows as the fitness function of classifier to generate the optimal extraction rules. The experiments are carried out in the Dadeng island of Xiamen in Fujian and part of the city of Zixing in Hunan, which are based on the domestic GF-1 and GF-2 remote sensing imageries. The results are compared with the SVM classifier and show that the proposed method can achieve better overall accuracy and Kappa coefficient, indicating that the proposed method is suitable for extraction of surface water bodies from remote sensing imagery of high spatial resolution.

artificial bee colony; water extraction; GF; object-oriented; remote sensing classification

2016-07-19;

2017-01-16

國(guó)家自然科學(xué)基金資助項(xiàng)目(41401475、41471366、41501448)；福建省測(cè)繪地理信息局(2015JX04)資助

鄧富亮(1982-)，男，博士后，講師，主要從事遙感圖像處理及應(yīng)用研究。E-mail：fldeng8266@xmut.edu.cn

章欣欣(1983-)，男，博士后，副教授，主要從事遙感圖像處理及應(yīng)用研究。E-mail：zhangxinxin@xmut.edu.cn

10.16030/j.cnki.issn.1000-3665.2017.03.21

P231.1

A

1000-3665(2017)03-0143-08