李紅凌，劉 斌，巴艷君

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，烏魯木齊 830046；2.新疆維吾爾自治區(qū)測繪科學(xué)研究院，烏魯木齊 830046)

0 引 言

我國干旱/半干旱區(qū)生態(tài)環(huán)境十分脆弱，該區(qū)域濕地的變化非常活躍，是環(huán)境變化的敏感指示器[1]，其轉(zhuǎn)移，收縮以及消失都能夠反映出該區(qū)域內(nèi)生態(tài)環(huán)境的質(zhì)量和容量。當(dāng)下，遙感在濕地方面的研究重點體現(xiàn)在濕地資源檢測、濕地信息分類與濕地的變化監(jiān)測等方面。就濕地遙感信息提取方法來說，自動化智能化一直是研究者追求的目標(biāo)[2]。針對高分辨率的遙感影像，為了改善地物目標(biāo)的分類精度，學(xué)者通常選擇面向?qū)ο蟮挠跋穹指罴夹g(shù)，并根據(jù)分割結(jié)果從對象中提取出光譜特征以及空間關(guān)系特征來改善對象間光譜混合問題[3]。

在面向?qū)ο蠓指钸^程中，當(dāng)分割閾值選擇不合理時會引起“欠分割”、“過分割”、“邊緣不匹配”等問題，為了滿足不同地物的最佳分割效果，可以利用地物特征異質(zhì)性來尋找各類別最適宜的分割尺度，最終獲得不同類型的影像對象[4]。由Baatz等人[5]于2000年提出的分形網(wǎng)絡(luò)演化方法(Fractal Net Evolution Approach，F(xiàn)NEA)是一種非常適用于高分辨率遙感影像的圖像分割算法，已被應(yīng)用與多個國際項目中，實踐證明這種基于統(tǒng)計的影像分割技術(shù)在影像分析中頗具潛力[6,7]。 FNEA技術(shù)的依據(jù)是局部最佳適應(yīng)法,該方法確保像元在統(tǒng)計對象異質(zhì)性最小的結(jié)果后進(jìn)行合并,即單個像元或像元集合與其相鄰的像元同時進(jìn)行計算與合并,兩個相鄰對象之間的差異通過啟發(fā)式處理過程進(jìn)行比較[8]。

在濕地信息遙感提取方面，范登科[9]等發(fā)現(xiàn)歸一化差異綜合水體指數(shù)在水體提取中水陸區(qū)分度高。朱長明[10]和姜浩[11]等根據(jù)水體的線性特征提出解決影細(xì)小河流的提取的方法。鄭麗娟等人提出面向?qū)ο蠼Y(jié)合DEM數(shù)據(jù)對研究區(qū)濕地進(jìn)行提取[12];謝靜、李曉娟等人采用多時相HJ_1影像利用面向?qū)ο蠓椒ㄍ瓿闪送赀_(dá)山以北的濕地遙感分類[13]。上述方法在水體提取上取得了較好的效果，對濕地信息提取有重要的參考價值。

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

塔里木河流域位于新疆天山山脈以南，昆侖山山脈以北，由葉爾羌河、和田河、阿克蘇河等匯合而成。塔里木河從西向東蜿蜒于塔里木盆地北部，上游地區(qū)地勢較為起伏，多為不平緩的沙漠地區(qū)，流域范圍從葉爾羌河流域起，到塔里木河尾閭的臺特瑪湖，全長2 421 km，是世界上第一最大內(nèi)陸河流域。該流域深處歐亞大陸腹地，四周環(huán)山，地理位置偏離海洋，日照充足，導(dǎo)致該區(qū)域內(nèi)降水量小，蒸發(fā)量大，表現(xiàn)為干旱大陸性氣候。根據(jù)《中華人民共和國濕地分類標(biāo)準(zhǔn)( GB/T 24708-2009)》采用二級分類系統(tǒng)，結(jié)合塔里木盆地水文地貌條件和濕地發(fā)育特征，濕地一級類型包括天然濕地與人工濕地，二級類型包括湖泊濕地、河流濕地、沼澤濕地、水庫、坑塘和人工養(yǎng)殖場六個類型。

試驗采用的遙感影像數(shù)據(jù)為2014年10月葉爾羌河中游區(qū)域RapidEye高分遙感影像，地理坐標(biāo)為39°4′58.58″N，77°36′33.01E～39°0′47.84″N，77°42′40.90″E，空間分辨率為5 m，波段范圍分別為藍(lán)波段(440～510 nm)，綠波段(520-590 nm)， 紅波段(630～685 nm)，紅邊波段(690～730 nm)和近紅外波段(760～850 nm)；1∶50 000 DOM數(shù)據(jù)用于研究區(qū)影像正射校正；土地利用現(xiàn)狀數(shù)據(jù)為研究區(qū)第二次土地調(diào)查項目成果，該數(shù)據(jù)為地表覆蓋分類輔助判讀的參考數(shù)據(jù)及分類結(jié)果精度驗證的主要數(shù)據(jù)源。研究區(qū)內(nèi)濕地類型有河流、水庫、沼澤、人工養(yǎng)殖場和坑塘。原始影像經(jīng)過輻射校正、正射校正、融合、鑲嵌等處理。

2 研究方法

2.1 主成分分析

主成分分析(Principal Component Analysis, PCA)由Pearson[14]于1901年發(fā)明，其原理是通過對協(xié)方差矩陣進(jìn)行特征分析，從而達(dá)到減少數(shù)據(jù)維數(shù)的目的，并且該過程要遵循數(shù)據(jù)集對方差的貢獻(xiàn)最大原則[15]。影像經(jīng)過主成分分析后得到的5個分量，前3個分量PC1、PC2和PC3累計方差百分比為99.65%,將前3個分量分別放入紅、綠、藍(lán)通道合成假彩色影像，與原真彩色影像相比色彩鮮明，各地物的區(qū)分度更明顯，其中深色區(qū)域反映出可能為水體的信息。PC1分量的灰度圖像中可以看出水庫、沼澤等地物顯示高亮度值，限定閾值可提取相應(yīng)目標(biāo)地物。根據(jù)研究區(qū)內(nèi)典型地物在PCA各分量中的波譜曲線顯示，水體在PC2波段上的值VALUEwater<-800，裸地、植被、建筑物等地物值VALUE(not water)>-800。從PC2分量的灰度圖像中可以看出水體呈現(xiàn)暗色調(diào)，其他地物為亮色調(diào).由此可知，在面向?qū)ο蟮臐竦靥崛≈?，選擇PC2分量作為特征參數(shù)，調(diào)整閾值可以有效提取水體信息。

圖1 主成分分析圖Fig.1 The spectrum of typical features in PCA plot

圖2 典型地物在 PCA 的波譜曲線圖Fig.2 Principal components analysis images

2.2 多尺度分割

FNEA的多尺度分割(Multiresolution Segmentation)[16]是一種基于像元自下而上的分割思維，其原理是從任一種子像元開始,將小的像元合并為大的對象，分割過程受尺度參數(shù)、形狀參數(shù)和緊致度的影響。尺度參數(shù)在軟件應(yīng)用中是不是一個具體的術(shù)語，其作用是限定影像對象的異質(zhì)性閾值，通過調(diào)整參數(shù)改變分割對象的規(guī)模(顏色同質(zhì)性是基于光譜顏色的標(biāo)準(zhǔn)差)；形狀參數(shù)用于限定對象的形狀以及優(yōu)化對象的均質(zhì)度[17]；緊致度參數(shù)可以進(jìn)一步表示對象的均質(zhì)度[18]；合并規(guī)則為合并后對象的異質(zhì)性不得大于事先給定的閾值。

由于研究區(qū)存在多種濕地類型，為解決不同地物目標(biāo)適宜的分割尺度，本研究通過ESP(Estimation of Scale Parameter)尺度評價工具[19]統(tǒng)計影像同質(zhì)性的局部變化LV(Local Variance)以及LV變化率值ROC(Rates Of Change)。ROC計算公式為：

式中：Li表示目標(biāo)層第i層對象層的平均標(biāo)準(zhǔn)差，而Li-1表示目標(biāo)層第i-1層中對象層的平均標(biāo)準(zhǔn)差。當(dāng)ROC最大時變化率曲線呈現(xiàn)峰值，峰值點有可能為目標(biāo)地物對應(yīng)最優(yōu)尺度。

參照ESP評價結(jié)果結(jié)合目視判別法，最終判斷尺度在240時適合水庫、沼澤等目標(biāo)較大的濕地類型提??；尺度為80時人工養(yǎng)殖場的邊緣匹配狀況較好；尺度為40時能準(zhǔn)確對河流、坑塘進(jìn)行分割。

圖3 ESP尺度分割評價圖Fig.3 ESP scale segmentation evaluation image

2.3 濕地分層提取

濕地包含的信息較為復(fù)雜，利用單一波段的光譜閾值難以區(qū)分。本研究在針對濕地分類所提取得特種光譜、紋理、形狀特征主要包括：

(1)主成分分析的PC1、PC2分量。

(2)歸一化植被指數(shù)NDVI(Normalized Difference Vegetation Index),通過植被在紅光和近紅外波段輻射差異，計算兩個波段的反射率和差、比值來判別地表植被覆蓋和長勢，其表達(dá)式為：

式中：Bnir表示近紅外波段；Br表示紅光波段。

(3)歸一化差異綜合水體指數(shù)CIWI(Combined Index of NDVI and NIR for Water- Identification)[20]。水體隨波長增長，其光譜反射率降低，而植被和土壤則相反，將NDVI和近紅外結(jié)合，利用乘性因子可以增強水體光譜特征值，抑制植被、河漫灘等地物信息，從而區(qū)分水體和其他地物。表達(dá)式為：

(4)紋理特征主要有灰度共生同質(zhì)性指數(shù)(GLCM Homogeneity)。

(5)形狀特征中形狀指數(shù)(Shape Index)和對象的像元長度(Length)。

分割過程中設(shè)置參數(shù)詳情：緊致度為0.5，形狀指數(shù)為0.4。分別在240、80和40三個尺度的分割結(jié)果上建立分類規(guī)則集，濕地提取流程如圖4所示。

圖4 濕地信息分層提取流程圖Fig.4 Flowchart of wetland information extraction

首先，在分割尺度為240層，通過PC2閾值對影像進(jìn)行水體粗提取，提取出可能是水體的區(qū)域，并根據(jù)歸一化差異綜合水體指數(shù)CIWI與形狀指數(shù)ShapeIndex進(jìn)一步從水體中提取水庫信息,同時，通過植被覆蓋指數(shù)NDVI與PC1兩個特征從非水體中提取出沼澤濕地。其次，在尺度為80時，通過調(diào)整CIWI指數(shù)和紋理指數(shù)GLCM的閾值來確立人工養(yǎng)殖場。然后，在尺度為40時通過PC2,CIWI提取較為細(xì)小的水體信息，將提取的細(xì)小水體圖斑合并，通過限定對象長度Length值區(qū)分河流和庫塘信息。同時，在該尺度上對可能為陸地的區(qū)域重新檢測水體信息，再次利用CIWI指數(shù)來提取坑塘信息。最后對分類結(jié)果進(jìn)行合并、平滑以及去除細(xì)小圖斑等處理。

3 結(jié)果與分析

將按照試驗步驟所得到的研究區(qū)濕地分類結(jié)果，可以得出結(jié)論基于FNEA的多尺度分層提取可以實現(xiàn)研究區(qū)內(nèi)，水庫、沼澤、河流、人工養(yǎng)殖場和坑塘等濕地類型的自動解譯(見圖5)，基本滿足生產(chǎn)應(yīng)用的需求，也對國內(nèi)外其他干旱區(qū)、內(nèi)陸地區(qū)的濕地遙感解譯提供參考。

為了更為客觀的驗證該提取結(jié)果的精度，本文參考第二次土地調(diào)查專題數(shù)據(jù)和野外實測數(shù)據(jù)，在研究區(qū)范圍內(nèi)各濕地類型隨機選取30～50個檢驗樣本點，統(tǒng)計樣本點在分類影像中的混淆矩陣，如表1所示?；煜仃囷@示，基于FNEA分割尺度下的多層次濕地信息提取方法在塔里木河流域內(nèi)濕地遙感解譯中總精度達(dá)到了87.7%，Kappa系數(shù)為0.842,各濕地類型的用戶精度與制圖精度明顯高于80%,取得了較好的結(jié)果。其中，水庫與沼澤的制圖均達(dá)到了90%，人工養(yǎng)殖場達(dá)到93.33%。另外，人工養(yǎng)殖場與河流的用戶精度甚至高于95%，達(dá)到了96.55%與96.15%。由于受到植被、河漫灘、陰影等地物的干擾，導(dǎo)致分類結(jié)果中出現(xiàn)個別錯分、漏分等現(xiàn)象。但是，整體而言，基于FNEA多尺度分割結(jié)合分層規(guī)則集的建立在研究區(qū)內(nèi)濕地解譯和提取中達(dá)到了較為理想的效果。

圖5 濕地分類圖Fig.5 The classifications of wetland

表1 濕地自動解譯混淆矩陣Table.1 The confusion Matrix of wetland Automatic interpretation

4 結(jié) 語

干旱/半干旱地區(qū)的稀疏植被、河漫灘、建筑物陰影等地物的光譜亮度值與濕地中沼澤、河流等地物非常相似，傳統(tǒng)的解譯方法很少對空間特征信息加以應(yīng)用，為解決單一光譜信息中同譜異物、同物異譜造成的地物界限不明確和椒鹽現(xiàn)象，該研究提出了基于面向?qū)ο缶W(wǎng)絡(luò)分形演化法并融入空間關(guān)系，多尺度多層次的提取濕地信息。結(jié)果表明，面向?qū)ο笏季S在針對高分辨率的遙感影像的濕地提取中體現(xiàn)出一定優(yōu)勢：

(1)FNEA分割結(jié)合ESP分割評價工具能夠快速、準(zhǔn)確的確立與地物相適應(yīng)的分割尺度并有效地將像元分割成對象目標(biāo)，避免花費大量的人力和時間去確立合理的分割尺度，同時，改善了后續(xù)分類中的地物混淆和椒鹽現(xiàn)象。

(2)主成分分析法不僅能夠去除影像中的冗余信息達(dá)到減維降噪的目的，研究顯示，主成分分析中的第一分量與第二分量是分離水體與其他信息的有效特征，在水體解譯中非常實用。

(3)在FNEA最優(yōu)分割尺度下的分類規(guī)則集能夠快速、自動提取研究區(qū)內(nèi)濕地信息，其架構(gòu)靈活可控，適于被移植和擴(kuò)展，經(jīng)換區(qū)試驗證明具有一定的普適性。同時，對其他區(qū)域的濕地信息識別與提取提供參考。

該研究針對多層次、多特征的干旱區(qū)/半干旱區(qū)濕地信息提取還不夠完善，有待針對干旱區(qū)/半干旱區(qū)的區(qū)域特點發(fā)現(xiàn)和拓展更多特征，探索更為高效準(zhǔn)確的規(guī)則集模型。

□

[1] 黃 方,張樹文,Jan de Leeuw,等.20世紀(jì)70年代以來松嫩平原烏裕爾河中下游流域的濕地演變[J].自然資源學(xué)報,2006,21(6):949-956.

[2] 朱長明，李均力，張 新，等.面向?qū)ο蟮母叻直媛蔬b感影像濕地信息分層提取[J].測繪通報,2014,(10):23-28.

[3] 汪 閩,萬其明,張大騫,等.光譜、形狀特征結(jié)合的多精度圖像分割算法與應(yīng)用[J].地球信息科學(xué)學(xué)報,2010,12(2):261-268.

[4] 黃慧萍，吳炳方．地物大小、對象尺度、影像分辨率的關(guān)系分析[J]．遙感技術(shù)與應(yīng)用,2006,21(3):243-248.

[5] Baatz M, Sch?pe A. Multiresolution segmentation: An opti-mization approach for high quality multi-scale image seg-mentation[C]∥ Angewandte Geographische Informations-verarbeitung XII, Strobl J, Blaschke T, Griesebner G(Eds),Wichmann, Heidelberg,Germany,2000.

[6] Schiewe J.Segmentation of high-resolution remotely senseddata-concepts, applications and problems[J]. InternationalArchives of Photogrammetry Remote Sensing And Spatial Information Sciences, 2002,34(4):380-385.

[7] Blaschke T, Burnett C, Pekkarinen A. Image segmentation methods for object-based analysis and classification[C]∥ Remote sensing image analysis: Including the spatial do-main.Springer Netherlands, 2006.

[8] Schiewe J, Tufte L, Ehlers M. Potential and problems of multi-scale segmentation methods in remote sensing[J].GeoBIT/GIS, 2001,6(1):34-39.

[9] 范登科,李 明,賀少帥.基于環(huán)境小衛(wèi)星CCD影像的水體提取指數(shù)法比較[J].地理與地理信息科學(xué),2012,28(2):14-19.

[10] 朱長明,駱劍承,沈占鋒,等.DEM 輔助下的河道細(xì)小線性水體自適應(yīng)迭代提取[J].測繪學(xué)報,2013,42(2):277-283.

[11] 姜 浩,馮 敏,肖 桐,王昌佐.基于線狀特征增強的 TM 遙感影像細(xì)小河流提取方法[J].測繪學(xué)報,2014,43(7):705-710.

[12] 鄭麗娟,李小娟,胡德勇,等.基于對象和輔助DEM的濕地信息提取 以洪澤湖濕地為例[J].遙感技術(shù)與應(yīng)用，2009,24(3):346-351．

[13] 謝 靜,王宗明,毛德華,等.基于面向?qū)ο蠓椒ê投鄷r相HJ-1影像的濕地遙感分類——以完達(dá)山北三江平原為例[J].濕地科學(xué),2012,10(4) :429-438．

[14] 王 鶯,王 靜,姚玉璧,等.基于主成分分析的中國南方干旱脆弱性評價[J].生態(tài)環(huán)境學(xué)報,2014,23(12):1 897-1 904.

[15] Pearson K. Principal components analysis[J]. The London Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1901,6(2):559.

[16] Benz UC, Hofmann P, Willhauck G, Lingenfelder M. Multi-Resolution, Object Oriented Fuzzy Analysis of Remote Sensing Data for GIS-Ready Information[J]. SPRSJournal of Photogrammetry & Remote Sensing, 2004,(58):239-258.

[17] Beijing Space Eye Innovation Technology Co, LT D. eCognition Training Files[EB/OL]. [2008 29224].http:∥www.bsei.com.cn/xz ecognition rjzl.html〗.

[18] Definiens AG. Definiens Developer 7 Reference Book[M]. München:Trappentreustr Press,2007:23-24.

[19] Dragut L, Tiede D, Levick SR.ESP: a tool to estimate scale parameter for multiresolution image segmentation of remotely sensed data[J]．International Journal of Geographical Information Science,2010,24(6):859-871.

[20] 莫偉華,孫涵.基于EOS/MODIS衛(wèi)星數(shù)據(jù)的洪澇災(zāi)害遙感監(jiān)測應(yīng)用技術(shù)研究[D]．南京：南京信息工程大學(xué)，2006.

·信息·

歡迎訂閱2016年《中國農(nóng)村水利水電》

《中國農(nóng)村水利水電》是由中國灌溉排水發(fā)展中心、水利部農(nóng)村水電及電氣化發(fā)展局、武漢大學(xué)、中國國家灌排委員會主辦的水利水電專業(yè)性技術(shù)期刊，內(nèi)容包括水利水電研究、設(shè)計、建設(shè)和管理的各個方面，信息量大，讀者面廣，且注重應(yīng)用技術(shù)的研究和推廣，理論聯(lián)系實際，具有很強的實用性。是全國(省、部)優(yōu)秀科技期刊，并榮獲“首屆國家期刊獎”，是“中國期刊方陣”的雙高期刊。是全國中文核心期刊、中國科技論文統(tǒng)計源期刊、中國學(xué)術(shù)期刊綜合評價數(shù)據(jù)庫來源期刊。

欄目設(shè)置：農(nóng)田水利、灌區(qū)建設(shè)與管理、節(jié)水灌溉、供水工程、水工建筑、泵站技術(shù)、水文水資源、水環(huán)境與水生態(tài)、水利經(jīng)濟(jì)、水電建設(shè)、防汛與抗旱、國外技術(shù)等。

月刊，大16開，每月15日出版，定價10.00元/冊，全年120元。

國內(nèi)訂閱：全國各地郵局，郵發(fā)代號38-49

國外訂閱：中國國際圖書貿(mào)易總公司(北京399信箱)，國外代號M4211，也可直接從編輯部訂閱

地 址：武漢大學(xué)(二區(qū))《中國農(nóng)村水利水電》編輯部 郵 編：430072

電 話：(027)68776133 68776880 傳 真：(027)68776133

聯(lián) 系 人：關(guān)良寶E-mail：xsdbjb@188.com

單位名稱：《中國農(nóng)村水利水電》讀者服務(wù)部

開戶銀行：中國銀行武漢市東湖開發(fā)區(qū)支行轉(zhuǎn)水大支行

賬 號：556057528310。