張　毅，張彥峰

（武漢大學(xué)遙感信息工程學(xué)院，湖北　武漢 430079）

基于光譜-空間特征的遙感影像水體提取方法

張毅，張彥峰

（武漢大學(xué)遙感信息工程學(xué)院，湖北武漢 430079）

地表水體是人類生存和社會(huì)發(fā)展不可取代的基本要素，使用遙感影像提取地表水體相比傳統(tǒng)實(shí)地調(diào)研的方式，具有高效、省時(shí)、省力的特點(diǎn)。傳統(tǒng)的水體提取方法一般僅使用一種或少量類型特征的結(jié)合，在處理高分辨率遙感影像時(shí)容易出現(xiàn)誤提取。針對(duì)以上問題，提出一種結(jié)合多種光譜和空間特征的自動(dòng)化水體提取方法。通過兩幅高分一號(hào)影像上的定性和定量實(shí)驗(yàn)結(jié)果表明，該方法能夠比傳統(tǒng)方法更穩(wěn)定更精確的提取影像中的水體信息。

遙感；水體提取；影像分類；光譜特征；空間特征

地表水體是人類生存和社會(huì)發(fā)展不可取代的基本要素。地表水體的研究對(duì)氣候模型，農(nóng)業(yè)穩(wěn)定性，河流變化，濕地監(jiān)管，地表水體調(diào)研和管理，環(huán)境監(jiān)測(cè)等領(lǐng)域起著重要的作用［1，2］。準(zhǔn)確地繪制出地表水體的空間分布是水文相關(guān)領(lǐng)域研究的前提，也是制定相關(guān)政策非常關(guān)鍵的一步［3］。傳統(tǒng)的以實(shí)地調(diào)研為主的水體提取方式需要耗費(fèi)大量的人力物力資源，且耗時(shí)較長(zhǎng)難以獲得及時(shí)的信息。遙感是一種利用非接觸手段獲取信息的方式，通過從遙感影像中解讀水體信息，大大降低了獲得地表水體信息的成本。

目前的國內(nèi)外研究中，遙感影像提取水體的主要研究主要集中在特征提取和水體提取算法兩個(gè)方面。合適的特征對(duì)正確的從遙感影像提取水體至關(guān)重要，目前在水體提取中研究的側(cè)重點(diǎn)為水體特征的研究主要包括：

1）光譜特征：遙感影像獲取的原始信息進(jìn)行輻射定標(biāo)、大氣校正后的結(jié)果，是水體提取過程中最基本的特征，大部分水體提取過程中均會(huì)使用光譜特征?？傮w上來說，水體在藍(lán)綠兩波段上呈現(xiàn)反射特性，在近紅外波段上呈現(xiàn)吸收的特性。根據(jù)水體和背景像素的光譜反射率差異，得到基于規(guī)則的方法得到影像中的水體，或使用光譜反射率作為特征訓(xùn)練分類器如支持向量機(jī)（SVM）［4］，對(duì)影像進(jìn)行二值分類得到水體和非水體信息。

2）水體指數(shù)：通過分析不同地物對(duì)遙感影像波段的響應(yīng)情況對(duì)波段進(jìn)行組合、代數(shù)運(yùn)算，從而增強(qiáng)水體信息，并抑制其他地物信息。常見的水體指數(shù)有：NDWI［5］，MNDWI［6］，AEWI［7］，HRWI［4］等。

3）紋理特征：使用紋理特征或紋理特征與光譜的結(jié)合，提取水體信息。文獻(xiàn)［8］在水體提取問題上提出了方向方差算子和紋理表相結(jié)合的多紋理特征結(jié)合的水體提取算法，在高分辨率灰度影像的實(shí)驗(yàn)上得到了很好的提取結(jié)果。文獻(xiàn)［9］比較了僅使用顏色特征和顏色特征加灰度共生矩陣角二階矩特征的提取效果，實(shí)驗(yàn)表明，增加紋理特征后，提取結(jié)果有了很大的提高。

4）顏色特征：對(duì)影像的紅綠藍(lán)波段反射率進(jìn)行顏色空間變化如HIS變換、CIELAB變換等，利用變化后的顏色空間特征作為水體提取的特征。文獻(xiàn)［10］研究了HIS變換后的色調(diào)（H）、飽和度（S）和亮度（I），對(duì)高分辨率遙感影像進(jìn)行水體提取。文獻(xiàn)［9］通過CIELAB變換，使用亮度分量（L）進(jìn)行分割，并使用顏色分量（a和b）和亮度分量的紋理特征構(gòu)成特征向量提取航空影像上的水體。

5）形狀特征：通過分割的方式得到對(duì)象，提取對(duì)象的形狀特征，利用水體和非水體對(duì)象在形狀特征上的差異，提取影像中的水體對(duì)象。文獻(xiàn)［11］使用緊致度（Compactness）、主要點(diǎn)（Criticalpoints）和對(duì)稱性（Symmetry）三個(gè)形狀特征成功區(qū)分了影像中的陰影和水體。

目前研究中，一般水體提取方法僅使用以上一類或兩類特征，然而，由于地表水體的復(fù)雜性，以上一種或兩種特征的結(jié)合并不能很好的從影像中提取水體信息。首先，“同物異譜”導(dǎo)致同樣的水體可能具有不同的光譜、水體指數(shù)、顏色特征，而“同譜異物”現(xiàn)象使得水體和背景信息在這些特征上存在混淆。其次，一般的紋理特征在計(jì)算時(shí)需要進(jìn)行取窗口的操作，在水體邊緣上的紋理信息會(huì)產(chǎn)出模糊的現(xiàn)象，無法完全反映水體的特點(diǎn)。最后，自然界中的水體往往具有不確定的形狀，因此，很多情況下形狀特征無法準(zhǔn)確區(qū)分水體和背景地物。針對(duì)以上問題，本文提出一種基于光譜-空間特征的遙感影像水體提取算法，結(jié)合了光譜、水體指數(shù)、像元長(zhǎng)度指數(shù)和面積特征，多種類型特征的結(jié)合更好的描述了遙感影像上水體的特征，魯棒地提取遙感影像中的水體信息。

1　基于光譜-空間特征的水體提取算法

本文方法使用了光譜、水體指數(shù)、像元長(zhǎng)度特征和面積特征，通過結(jié)合光譜和空間特征，一方面減輕了“同物異譜”和“同譜異物”問題對(duì)水體和非水體像素的影響，另一方面引入的像元長(zhǎng)度特征綜合影像中像素的局部結(jié)構(gòu)特征，有效的反映了像元鄰域內(nèi)的相似性，且避免了求取紋理特征時(shí)滑動(dòng)窗口操作可能產(chǎn)生的問題。本文方法總體上可以分為兩個(gè)部分，第一部分使用雙閾值分割去除影像中的大部分非水體像素，第二部分將第一步得到的水體待選像素聚集成水體待選對(duì)象，并使用面積特征和水體指數(shù)去除非水體對(duì)象，最后對(duì)得到的水體對(duì)象進(jìn)行區(qū)域生長(zhǎng)，得到精確的水體邊界信息。方法的整體流程如圖1所示。

圖1　水體提取方法整體流程圖

1.1輻射定標(biāo)

首先對(duì)遙感影像進(jìn)行輻射定標(biāo)預(yù)處理，根據(jù)公式（1）將衛(wèi)星各載荷的通道觀測(cè)技術(shù)值DN轉(zhuǎn)換為衛(wèi)星載荷入瞳處等效表觀輻射亮度：

式中：Gain為定標(biāo)斜距；DN為衛(wèi)星載荷觀測(cè)值；Bias為定標(biāo)截距，單位為；以上參數(shù)均由實(shí)際衛(wèi)星數(shù)據(jù)給出；

再根據(jù)如下公式，得到大氣頂端的表觀輻射量：

其中，L為上一步得到的入瞳處的表觀輻射量，d為年積日即影像獲取時(shí)間為一年中的第幾天，θ為太陽高度角，以上數(shù)據(jù)可從影像的元數(shù)據(jù)獲得，ESUN為太陽表觀輻射率均值，不同的衛(wèi)星傳感器會(huì)給定相應(yīng)的取值。

1.2基于光譜-空間特征的雙閾值分割

本步驟的目標(biāo)為結(jié)合光譜和像元的局部結(jié)構(gòu)特征通過兩步閾值分割處理，去除不符合要求的非水體像素。水體指數(shù)從光譜的角度反映了水體和非水體像素的差異，其取值范圍為［-1，1］，水體的NDWI取值一般為正，且取值越大該像素為水體的可能性也越大，水體指數(shù)的計(jì)算公式如下所示：

其中，green為遙感影像上的綠波段，NIR為近紅外波段的輻射量；

其次，本文引入像元長(zhǎng)度指數(shù) （PixelLength Index，PLI），衡量水體內(nèi)部的相似性。像元長(zhǎng)度指數(shù)是像元結(jié)構(gòu)特征集合［13］（StructuralFeatureSet，SFS）中的一種特征，該特征考慮了相鄰像元的相似性，綜合了像素間的上下文信息，具有減少同質(zhì)性區(qū)域的噪聲和光譜變化的特性。像元長(zhǎng)度特征反映了像元在其鄰域內(nèi)相似性的最大值。由于同一水體對(duì)象內(nèi)部的光譜差異較小，且不同水體對(duì)象之間的形態(tài)差異較大，使用像元長(zhǎng)度指數(shù)既可以反映影像內(nèi)部水體之間的相似性，也減小了水體對(duì)象形態(tài)差異造成的像元結(jié)構(gòu)差異。

在獲得影像的像元長(zhǎng)度指數(shù)時(shí)，為了加大水體之間的相似性，以及水體和非水體之間的差異性，本文使用拉伸后的水體指數(shù) SNDWI（Stretched NDWI）特征計(jì)算像元長(zhǎng)度指數(shù)。由公式得到水體指數(shù)特征（NDWI），并根據(jù)公式（4）進(jìn)行拉伸。

使用SNDWI得到的光譜相似性測(cè)度為下式：

其中，phi表示第i個(gè)方向線上中心像素和鄰域像素之間的異質(zhì)性，Psncen為中心像素的SNDWI特征，Psncen為鄰域像素的SNDWI特征。

對(duì)預(yù)處理后的影像數(shù)據(jù)提取水體指數(shù)和像元長(zhǎng)度指數(shù)，通過雙閾值分割，得到水體待選像素，具體步驟如下所示：首先，對(duì)影像的SNDWI特征計(jì)算得到的像元長(zhǎng)度特征使用給定閾值進(jìn)行二值化，從像元的局部相似性，排除了影像中大量的非水體像素。獲得上步中通過閾值分割的像素的NDWI特征值，使用大津法對(duì)其進(jìn)行閾值分割，利用水體和非水體像素在水體指數(shù)上的差異進(jìn)一步排除影像中的非水體像素。

1.3基于雙閾值分割結(jié)果的水體對(duì)象提取

經(jīng)過上一步的處理，得到了通過光譜和空間特征約束的像素的二值結(jié)果。首先，根據(jù)二值圖上像素的連通性，將8鄰域內(nèi)相鄰的連接像素聚集為對(duì)象，生成水體待選對(duì)象。

根據(jù)聚集的水體對(duì)象的面積進(jìn)行閾值分割，給定面積閾值，若對(duì)象的面積大于給定的閾值則認(rèn)為該對(duì)象為大面積水體待選對(duì)象，否則將其作為小面積水體待選對(duì)象。給定一個(gè)較大的面積閾值，如10000，則得到的大面積水體待選對(duì)象為面積很大的平滑區(qū)域，且具有較高的NDWI值。由于遙感影像上存在面積大且水體指數(shù)高的弱紋理非水體地物概率非常小，故接受大面積水體待選對(duì)象為大面積水體對(duì)象。而由于小面積的對(duì)象中含有背景地物的概率較高，故對(duì)小面積水體待選對(duì)象將進(jìn)行進(jìn)一步的處理以排除背景地物。

獲得小面積水體待選對(duì)象覆蓋范圍內(nèi)所有像素的NDWI特征，根據(jù)其直方圖形狀使用峰谷的方法（peak-and-valleymethod）獲得閾值再一次進(jìn)行閾值分割。為了減小分割時(shí)的漏檢率，因此選擇最左邊的谷點(diǎn)對(duì)應(yīng)的NDWI值作為分割的閾值。判斷直方圖上峰谷的方式為：

1）如果pn≥pn+1＆pn≥pn-1，則為直方圖上的峰值；

2）如果pn≥pn+1＆pn≥pn-1，則為直方圖上的谷值；

其中，pn表示直方圖上n位置對(duì)應(yīng)的頻率。

然而，由上述條件得到的結(jié)果往往會(huì)包含較多的偽峰和偽谷，因此，通過對(duì)直方圖進(jìn)行核函數(shù)估計(jì)（Kerneldensityestimation，KDE）得到平滑后的概率密度曲線（直方圖），從平滑后的結(jié)果中提取出峰值和谷值。核函數(shù)密度估計(jì)是非參數(shù)的估計(jì)隨機(jī)變量概率密度函數(shù)的一種方法，常被從來實(shí)現(xiàn)數(shù)據(jù)平滑的目的。其計(jì)算公式如下所示：

其中，Kh（·）為核函數(shù)，本文使用高斯核函數(shù)。圖2給出了經(jīng)過第一步處理后剩余像素的NDWI直方圖，以及對(duì)該直方圖進(jìn)行核函數(shù)估計(jì)的結(jié)果。從圖中可以看出，藍(lán)色部分的直方圖中存在偽峰值，而平滑后的概率密度曲線上不存在偽峰值，因此，使用平滑后的結(jié)果可以得到正確的閾值。

圖2　影像NDWI直方圖核函數(shù)估計(jì)結(jié)果

如果概率密度估計(jì)得到的結(jié)果中不包含谷值點(diǎn)，則說明經(jīng)過第一步的處理后，剩下的小面積水體待選對(duì)象中不包含非水體像素，則不對(duì)其進(jìn)行進(jìn)一步的閾值分割。

由以上兩步可以看出，本文方法旨在對(duì)大面積的水體使用較小的NDWI閾值進(jìn)行分割，而對(duì)面積較小的水體使用較高的NDWI閾值分割，并且在使用PLI進(jìn)行閾值分割時(shí)去除了大量的與水體具有類似NDWI取值的背景像素。

由于水體邊緣像素的PLI和NDWI特征比水體中心像素小，故在上述處理中可能被錯(cuò)誤的去除，因此，在得到水體對(duì)象后，通過區(qū)域生長(zhǎng)的方式得到更精確的水體邊界。將所有水體邊緣像素加入種子點(diǎn)堆棧，若種子點(diǎn)邊緣像素與其鄰域像素之間異質(zhì)性測(cè)度小于給定閾值，則認(rèn)為該鄰域像素為水體，并將其加入種子點(diǎn)堆棧，重復(fù)進(jìn)行上述操作直到所有種子點(diǎn)均被處理完成。這里使用的異質(zhì)性測(cè)度為光譜角，因?yàn)樵撎卣骶C合考慮了影像的所有光譜特征，計(jì)算公式如（7）所示，給定的閾值為thsa。

其中，SA（x，y）表示光譜向量之間的光譜角，xi和yi分別表示種子點(diǎn)像素和鄰域像素在i波段上的光譜反射率，n表示影像的波段個(gè)數(shù)。對(duì)于具有可見光和近紅外四個(gè)波段的高分辨率遙感影像來說，x和y均為四波段波譜反射率組成的四維向量。

2　實(shí)驗(yàn)與結(jié)果

本文使用兩幅高分一號(hào)衛(wèi)星的8米分辨率多光譜影像進(jìn)行水體提取實(shí)驗(yàn)。高分一號(hào)衛(wèi)星是我國高分辨率對(duì)地觀測(cè)系統(tǒng)重大專項(xiàng)的第一個(gè)衛(wèi)星，具有重要的戰(zhàn)略意義和地位。其8m分辨率的多光譜衛(wèi)星有近紅外和可見光四個(gè)波段，重訪周期為4d，具有較高空間分辨率和較短的回訪周期的特點(diǎn)，因此，具有較高的應(yīng)用價(jià)值。本文使用兩幅高分一號(hào)數(shù)據(jù)進(jìn)行實(shí)驗(yàn)，通過對(duì)比人工勾選的水體提取結(jié)果定性和定量的比較了本文提取方法與使用光譜特征的支持向量機(jī)（典型的水體提取方法）的提取結(jié)果。

本文使用如下方法進(jìn)行定量的精度評(píng)點(diǎn)：對(duì)比參考影像使用混淆矩陣統(tǒng)計(jì)，得到正確提取的水體像素（TP），錯(cuò)誤提取的水體像素（FP），正確提取的非水體像素（TN），錯(cuò)誤提取的非水體像素（FN），四類的像素個(gè)數(shù)。使用統(tǒng)計(jì)結(jié)果得到水體的生產(chǎn)者精度PA=TP/（TP+FN），生產(chǎn)者精度反映了水體的漏檢率，使用者精度UA=TP/（TP+FP），使用者精度反映了水體的誤檢率，總體精度OA=（TP+TN）/N，其中N為影像中的像素總體個(gè)數(shù)和kappa系數(shù)（如式（8）所示）。

其中，

圖3（a）顯示了本文實(shí)驗(yàn)中所用的第一幅影像（影像一），圖3（b）顯示了相應(yīng)的參考水體掩膜。影像一相幅大小為4548*4500，經(jīng)度范圍為120.0874度到 120.5829度，緯度范圍為 40.0612度到40.4397度。從影像中可以看出大面積水體的光譜特征并不統(tǒng)一，其大部分面積呈藍(lán)色，水體左下角顏色變深，并且在水體內(nèi)部可以看到條紋狀的差異。影像中的背景地物主要由黃色的裸地、深色山體和亮建筑物組成。

圖3　實(shí)驗(yàn)影像一（a為偽彩色顯示結(jié)果，b為參考水體掩膜）

圖4　影像一的提取結(jié)果（a為SVM結(jié)果，b為本文方法結(jié)果）

表1　影像一精度評(píng)定結(jié)果

圖4給出了使用SVM和本文方法提取的結(jié)果，表1給出了影像一提取結(jié)果的定量評(píng)價(jià)。從提取結(jié)果中可以看出，兩種方法均較好的提取出了遙感影像中的水體，本文方法得到的大面積水體更為完整，小面積的誤提取更少。從定量的結(jié)果可以看出，本文方法在四個(gè)指標(biāo)上均優(yōu)于SVM的提取結(jié)果。其中，使用者精度提高了約3%，說明本文方法極大的減小了誤檢情況的發(fā)生。

圖5　實(shí)驗(yàn)影像二（a為偽彩色顯示結(jié)果，b為參考水體掩膜）

圖5（a）顯示了本文實(shí)驗(yàn)中所用的第二幅影像（影像二），圖5（b）為其相應(yīng)的參考水體掩膜。影像二相幅大小為4548*4500，經(jīng)度范圍為120.1701度到120.6267度，緯度范圍為39.1023度到30.4952度。從影像中可以看出，大面積的水體本身的光譜發(fā)生了較大的變化，左下方主要呈現(xiàn)藍(lán)色，而右上方呈現(xiàn)綠色，說明右上方的水體在紅波段的反射率較高。不僅如此，影像中存在小面積的暗水體。影像二的背景信息復(fù)雜，既包括亮建筑物，也包括暗建筑物，還存在較多的植被、裸地以及陰影信息。

圖6給出了SVM和本文方法對(duì)影像二的提取結(jié)果，表2給出了定量的精度評(píng)價(jià)結(jié)果。與人工勾選的參考影像相比可以看出，SVM和本文方法均有效提取出了影像中的主要水體信息，且在小面積水體上均存在一定的誤檢問題，但整體上本文方法得到的結(jié)果更為完整，且誤檢的像素更少。從定量評(píng)價(jià)的角度可以看出，本文的生產(chǎn)者精度比SVM高約10%，使用者精度高約13.5%，整體精度高2.3%，因此，使用本文方法較好的提高了水體的提取精度。

圖6　影像二的提取結(jié)果（a為SVM結(jié)果，b為本文方法結(jié)果）

表2　影像二精度評(píng)定結(jié)果

從總體上了來說，本文方法的提取精度均高于SVM的提取精度，多種特征的結(jié)合以及像元長(zhǎng)度指數(shù)的引入，有效的提高了本文方法在大面積水體提取結(jié)果上的完整度，且有效的去除了小面積的噪聲。不僅如此，SVM為監(jiān)督分類方法，使用該方法進(jìn)行水體提取時(shí)，需要人工選擇訓(xùn)練樣本，并通過多次實(shí)驗(yàn)選取合適的參數(shù)。而本文方法為自動(dòng)提取方法，計(jì)算速度快。

3　結(jié)束語

針對(duì)目前遙感影像水體提取中使用的特征存在的多種問題，本文提出了一種基于光譜-空間特征的遙感影像水體提取算法。本文方法結(jié)合像元長(zhǎng)度指數(shù)、光譜指數(shù)以及對(duì)象的面積特征，從光譜和空間特征的角度區(qū)分影像中的水體和非水體像素。通過兩幅高分一號(hào)數(shù)據(jù)，定性和定量的比較了本文方法的有效性。從總體上來說，本文方法能夠高效自動(dòng)化的提取遙感影像上的水體。

［1］ Du Z，Linghu B，Ling F，et al.Estimating surface water area changes using time-series Land sat data in the Qingjiang River Basin，China［J］.Journal of Applied Remote Sensing，2012，6 （1）：063609-063609.

［2］ Sun F，Sun W，Chen J，et al.Comparison and improvement of methods for identifying water bodies in remotely sensed imagery［J］.International journal of remote sensing，2012，33 （21）：6854-6875.

［3］ Gessner M O，Hinkelmann R，Nützmann G，et al.Urban water interfaces［J］.Journal of Hydrology，2014，514：226-232.

［4］ Yao F，Wang C，Dong D，et al.High-Resolution Mapping of Urban Surface Water Using ZY-3 Multi-Spectral Imagery［J］. Remote Sensing，2015，7（9）：12336-12355.

［5］ McFeeters S K.The use of the Normalized Difference Water Index（NDWI）in the delineation of open water features［J］.International journal of remote sensing，1996，17（7）：1425-1432.

［6］ Xu H.Modification of normalised difference water index （NDWI）to enhance open water features in remotely sensed imagery［J］.International Journal of Remote Sensing，2006，27 （14）：3025-3033.

［7］ Feyisa G L，Meilby H，F(xiàn)ensholt R，et al.Automated Water Extraction Index：A new technique for surface water mapping using Landsat imagery［J］.Remote Sensing of Environment，2014，140：23-35.

［8］ Wang H，Pan L，Zheng H.Multi-texture-model for water extraction based on remote sensing image［C］//Image and Signal Processing，2008.CISP'08.Congress on.IEEE，3：710-714.

［9］ ZhaoM，ZhangY，ShangH，et al.Water Area Extraction from RGB Aerophotograph Based on Chromatic and Textural Analysis［C］//GEO Processing 2011，The Third International Conference on Advanced Geographic Information Systems，Applications，and Services：46-52.

［10］ 柳稼航.基于視覺特征的高分辨率光學(xué)遙感影像目標(biāo)識(shí)別與提取技術(shù)研究［D］.上海交通大學(xué)，2011.

［11］ Li B，Zhang H，Xu F.Water extraction in high resolution remote sensing image based on hierarchical spectrum and shape features［C］//IOP Conference Series：Earth and Environmental Science.IOP Publishing，2014，17（1）：012123.

［12］ 黃昕.高分辨率遙感影像多尺度紋理，形狀特征提取與面向?qū)ο蠓诸愌芯浚跠］.武漢：武漢大學(xué)，2009.

P407.8