許 瑩，張友靜,2，張 琴

（1.河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京210098；2.河海大學(xué) 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，江蘇 南京210098； 3.南京物聯(lián)網(wǎng)應(yīng)用研究院有限公司，江蘇 南京210013）

基于高分影像的水利空間要素提取規(guī)則集構(gòu)建

許 瑩1，張友靜1,2，張 琴3

（1.河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京210098；2.河海大學(xué) 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，江蘇 南京210098； 3.南京物聯(lián)網(wǎng)應(yīng)用研究院有限公司，江蘇 南京210013）

針對高分遙感水利要素自動(dòng)提取問題，提出了一種在融合處理后的高分一號(hào)（GF-1）遙感影像下表征和提取水利基礎(chǔ)空間要素的方法。依據(jù)水利要素與水體的空間關(guān)系，采取先進(jìn)行土地利用/覆蓋分類，再提取水利設(shè)施要素的技術(shù)思路，構(gòu)建了較為完整的水利空間要素提取規(guī)則集。利用混淆矩陣和矢量化圖像對提取精度和疊合精度進(jìn)行評(píng)價(jià)。結(jié)果表明，該規(guī)則集用于高分辨率影像水利基礎(chǔ)空間要素提取的效果較好。其中，河流提取精度和疊合精度最高，分別為97.1%和94.1%；河流阻斷物的提取精度和疊合精度最低，分別為86.0%和81.4%；其他水利要素總體精度較高。

GF-1影像；水利基礎(chǔ)空間要素；特征提??；面向?qū)ο蠓椒?；?guī)則集構(gòu)建

水利基礎(chǔ)空間要素（簡稱水利要素）是自然與人工構(gòu)造物的有機(jī)集成，主要包括兩類：流域自然特征信息要素，包含流域河湖水系分布、匯水單元及相互連接關(guān)系、流域地表覆蓋類型等；流域水利工程分布要素，包含水庫、堤防、水閘、渠道等。水利要素類型多，空間分布廣。目前利用遙感影像進(jìn)行水利要素提取的研究大多集中于河、湖和水系等流域自然要素的提取，主要方法有差值法[1]、水體指數(shù)法[2]、多波段譜間關(guān)系法[3-5]以及決策樹法[6]。對于山區(qū)細(xì)小水體和半干旱區(qū)枯季山間河道，楊樹文[7]等利用TM影像多波段譜間關(guān)系改進(jìn)法，實(shí)現(xiàn)了較高精度的提取。閆霈[8]等利用增強(qiáng)型水體指數(shù)Modified NDWI提取半干旱地區(qū)水體。駱劍承[9]等提出了“全域-局部”迭代轉(zhuǎn)換的信息提取方法，通過全域分割分類和局部分割分類的迭代算法，實(shí)現(xiàn)了高精度水體信息的自動(dòng)提取。然而對于水利要素中的水利工程要素提取，目前研究較少。由于我國高分辨率影像數(shù)據(jù)的可獲得性得到極大改善，因此研究基于高分辨率影像數(shù)據(jù)的水利要素提取方法，對于深入探討高分辨率影像的地學(xué)響應(yīng)，推動(dòng)高分影像數(shù)據(jù)水利應(yīng)用具有重要意義和實(shí)用價(jià)值。本文通過分析各類水利要素的光譜特征、紋理特征及其在高分辨率遙感影像上的特征，研究水利要素空間分布和空間關(guān)系特點(diǎn)，構(gòu)建基于高分辨率影像的水利要素提取方法和規(guī)則集，為利用高分遙感影像提取和更新水利要素提供技術(shù)支持。

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

選取湖北省荊州市城鄉(xiāng)結(jié)合部為研究區(qū)，區(qū)域內(nèi)土地覆被類型多樣，包括城鎮(zhèn)建設(shè)用地、耕地、道路和水體等；大小河流、渠系縱橫，水利設(shè)施眾多；對于水利要素提取具有較強(qiáng)的代表性。

研究數(shù)據(jù)為2013-08-06的GF-1影像，全色波段空間分辨率為2 m，多光譜波段空間分辨率為8 m，有紅、綠、藍(lán)、近紅外4個(gè)波段。GF-1衛(wèi)星于2013-04-26成功發(fā)射，是我國高分辨率對地觀測系統(tǒng)的首發(fā)星[10-11]，搭載6臺(tái)多光譜相機(jī)（2臺(tái)2 m全色/8 m多光譜相機(jī)，4臺(tái)16 m 多光譜寬幅相機(jī)），重訪周期為4 d，能夠提供區(qū)域內(nèi)同時(shí)相的高分辨率影像。研究所使用的驗(yàn)證數(shù)據(jù)為經(jīng)過矢量化的同范圍影像，空間分辨率為1 m。

為了消除大氣對地物反射的影響，對遙感圖像進(jìn)行輻射校正和圖像融合，將研究區(qū)同一時(shí)刻空間分辨率為8 m的多光譜數(shù)據(jù)和2 m的全色數(shù)據(jù)融合。融合影像既具有全色影像的高分辨率，又具有多光譜特征，有利于目視解譯和地物信息提取。

2 研究方法

2.1 研究思路

利用水利要素空間分布和空間關(guān)系等特點(diǎn)，在進(jìn)行土地利用/覆蓋初步分類的基礎(chǔ)上，提取水利要素。水利要素中流域自然要素和人工建筑要素的分布均與基礎(chǔ)土地覆被類別存在一定的空間關(guān)系，因此首先進(jìn)行土地利用/覆蓋分類，可在一定程度上簡化水利要素的規(guī)則集提取過程，提高提取精度。

2.2 土地利用/覆蓋分類

本文主要采用面向?qū)ο蠓诸惙?。首先將影像分割成對象特征，然后再進(jìn)行分類研究；能夠充分利用影像的光譜信息、紋理特征和拓?fù)潢P(guān)系，并綜合考慮光滑度、緊致度等空間特征，生成同質(zhì)的影像對象，從而大大提高分類精度；并能減少面向像元分類方法可能產(chǎn)生的“椒鹽”現(xiàn)象。

通過對實(shí)際區(qū)域的了解和對照土地利用圖進(jìn)行了初步的目視判讀，將實(shí)驗(yàn)區(qū)大致分為5個(gè)土地覆被類型：水體、耕地、植被、建筑和道路。確定土地覆被最適宜的分割尺度；多尺度分割后，采用模糊分類的方法實(shí)現(xiàn)地物類別的提??；特征選擇中選取影像對象在不同波段光譜特征集或選取空間特征集、紋理特征集。在本文中，根據(jù)各地物類別在高分辨率影像上的表征，構(gòu)建各類別的初始特征空間，主要包括光譜特征、形狀特征、紋理特征及自定義特征。光譜特征包括4個(gè)波段的平均值、最大差分、亮度值、比率和方差等17個(gè)特征；紋理特征采用灰度共生矩陣計(jì)算的平均值、對比度、方差等8 個(gè)特征；形狀特征為面積、形狀指數(shù)、密度、長度/寬度等8個(gè)特征；自定義特征為歸一化植被指數(shù)（NDVI）、歸一化水體指數(shù)等4個(gè)特征，初步確定的特征空間包括37個(gè)特征。通過選取樣本，統(tǒng)計(jì)各類樣本的特征值，建立樣本信息表并離散化得到?jīng)Q策表，利用粗糙集和蟻群算法相結(jié)合的屬性約簡算法對決策表進(jìn)行屬性約簡，最后得到的各要素類別的約簡特征為：NDVI、B1波段比率（Ratio B1）、亮度（Brightness）、長寬比（Length/ Width）、長度（Length）、灰度共生矩陣的熵（GLCM entropy）等6個(gè)特征，通過決策樹C4.5算法，對特征選擇后的樣本特征數(shù)據(jù)建立決策規(guī)則。

2.3 水利要素提取

水利要素主要包括河湖等的自然流域要素以及堤防、水閘等人工建筑要素，與基礎(chǔ)土地覆被類別存在位置等空間關(guān)系。例如，河流、湖泊、水庫和渠道是以水體為基礎(chǔ)類別的二次分類；堤防、水閘一般是以道路、建筑為基礎(chǔ)類別的二次分類；堤防是沿河、湖等兩岸的人工建筑，一般以道路或植被為基礎(chǔ)類別進(jìn)行二次分類。因此，利用土地利用/覆蓋分類的結(jié)果進(jìn)行水利要素的識(shí)別，確定各水利要素適宜的多尺度影像分割尺度，再通過對各水利要素的形狀、大小、紋理、位置關(guān)系、鄰接關(guān)系等特征綜合分析進(jìn)行水利要素的提取。

對水利要素進(jìn)行提取時(shí)，需著重考慮水利要素的形狀、大小、空間位置和鄰接關(guān)系等。本實(shí)驗(yàn)中采用的相關(guān)特征主要有：①形狀特征，包括面積、形狀指數(shù)和長寬比等；②鄰接關(guān)系特征，如堤防與河流、湖泊等存在鄰接關(guān)系，渠道與耕地等存在鄰接關(guān)系。

2.3.1 流域自然要素提取

在土地利用/覆蓋分類獲取的水體類別中進(jìn)行二次分類得到流域自然要素。本文提取的流域自然要素包括河流、湖泊和坑塘。人工建筑要素中的渠道因其形狀特征，也一并進(jìn)行提取。河流面積較大，且在空間上呈不規(guī)則曲線流動(dòng)；渠道面積小，且多為沿某一方向直線延伸；坑塘形狀多呈面積較小的多邊形。據(jù)此，分別建立提取規(guī)則。

2.3.2 水利工程分布要素提取

人工建筑要素的提取主要包括渠道、堤防、水閘和橋梁，其中渠道已提取出來。堤防是指在江、河、湖、海沿岸或水庫區(qū)、分蓄洪區(qū)周邊修建的土堤或防洪墻等，通常與河流相鄰，隨河流呈線狀延伸。實(shí)驗(yàn)區(qū)中堤防主要存在于河流兩岸，與防范對象相鄰接。根據(jù)堤防這一空間特征，以河流為主體對象，尋找與河流相鄰的地物，將其定義為堤防，并進(jìn)行平滑處理，求得主要河流兩側(cè)的堤防。從遙感影像上看，水閘和橋梁為切斷河流的主要地物，本文將橋梁和水閘統(tǒng)稱為河流阻斷物。根據(jù)河流阻斷物的特征，采用距離函數(shù)和方向函數(shù)來實(shí)現(xiàn)河流阻斷物的提取。首先在土地利用/覆蓋分類的基礎(chǔ)上，通過距離函數(shù)，設(shè)定一定閾值，搜尋這一閾值范圍內(nèi)除水體以外的類別，但僅用距離函數(shù)，河流兩旁的道路等地物可能會(huì)被一并提取，因此，需要同時(shí)考慮到河流阻斷物的方向性。河流阻斷物通常與河流相交，而河流兩旁的道路隨河流同向延伸，因此利用方向特征設(shè)定一定的閾值，將河流兩岸距離范圍內(nèi)的地物剔除，實(shí)現(xiàn)河流阻斷物的提取。渠道寬度較小，跨越渠道的河流阻斷物以小路居多，方向以相切為主。本文主要對跨越河流的主要道路進(jìn)行判別。

3 結(jié)果與分析

3.1 規(guī)則集建立

利用面向?qū)ο蠓诸愃枷耄罁?jù)水利要素與水體的空間關(guān)系，構(gòu)建了較為完整的水利要素提取規(guī)則集，如表1。

3.2 土地利用/覆蓋分類結(jié)果評(píng)價(jià)

經(jīng)過初始分類后，得到土地覆被分類結(jié)果如圖1所示。對分類結(jié)果作精度評(píng)價(jià)，實(shí)驗(yàn)中用于驗(yàn)證的樣本數(shù)據(jù)為在ENVI軟件中多邊形隨機(jī)采樣獲得，生成初始分類混淆矩陣見表2。

由表2可知，土地覆被總體分類精度較高，為85.8%；Kappa系數(shù)為0.84。其中水體分類精度最高，達(dá)到98.9%，因其在高分辨率影像上表征明顯，特點(diǎn)突出，易于提取。農(nóng)田與植被光譜特征差異小，雖綜合考慮了光譜特征、紋理特征等因素，但混分仍較嚴(yán)重；城區(qū)中道路和建筑的光譜差異不明顯，道路形狀細(xì)小破碎，建筑和道路也有一定程度的混分。

根據(jù)本文研究思路，水體類別是水利要素提取的基礎(chǔ)要素，與水利要素的提取關(guān)聯(lián)較大。因此認(rèn)為分類結(jié)果較為理想，可作為下一步提取水利要素的基礎(chǔ)。

表1 水利要素提取規(guī)則集

圖1 土地覆被分類

表2 土地覆被分類精度評(píng)價(jià)

3.3 水利要素提取結(jié)果評(píng)價(jià)

水利要素提取的結(jié)果如圖2所示。幾種典型的水利要素提取結(jié)果如圖3所示。本文對水利要素提取結(jié)果的評(píng)價(jià)方式是將實(shí)驗(yàn)中提取的結(jié)果與經(jīng)過矢量化得到的校驗(yàn)數(shù)據(jù)疊加，進(jìn)行空間一致性檢驗(yàn)。該方法能在一定程度上避免采用混淆矩陣等方法帶來的抽樣誤差問題，可靠性和準(zhǔn)確度更高。各水利要素的提取精度見表3。

由分類結(jié)果可知，總體提取精度和疊合精度較好。其中，河流提取精度最高，達(dá)到97.1%；河流阻斷物提取精度最低，為86.0%。河流的疊合精度最高，達(dá)到94.1%；河流阻斷物的疊合精度最低，為81.4%。影響提取精度的可能原因有：

1）河流面積較大，在高分辨率影像上表征明顯，特征突出，提取精度最高。部分細(xì)小河流較為破碎，易與城區(qū)中的建筑物、植被及其陰影等混分，提取不夠完整。檢驗(yàn)數(shù)據(jù)雖與原始影像不是同一時(shí)期，但研究區(qū)為城鄉(xiāng)結(jié)合部，防洪設(shè)施較為完備，水體面積變化不大。

2）坑塘因其面積較小，且形狀各異，易與建筑等地物混分，導(dǎo)致提取的疊合精度低；同時(shí)，提取精度也與分割尺度的選擇有關(guān)。

3）河流阻斷物的類別眾多，大小不一，與水體等地物的空間位置關(guān)系各異，容易導(dǎo)致河流阻斷物的漏提或錯(cuò)提。本文中的河流阻斷物包含水閘但不限于水閘，水閘與小型橋梁在高分辨率影像上的表征差異較小，較難區(qū)分，難以單獨(dú)提取。

4）堤防的提取精度較高，到達(dá)90.1%，但其疊合精度相對較低，僅為83.8%。原因可能是提取時(shí)設(shè)定河流兩側(cè)固定距離的地物為堤防，而實(shí)際堤防與河邊距離并不完全一致，因此少部分堤防提取結(jié)果與實(shí)際堤防不能疊合。

圖2 水利要素提取

圖3 典型水利要素提取結(jié)果

表3 水利要素提取精度評(píng)價(jià)

4 結(jié) 語

1）實(shí)驗(yàn)中通過多尺度分割選擇各地物最適宜的分割尺度，利用粗糙集和蟻群算法相結(jié)合的特征選擇模型，簡化了繁瑣的特征集創(chuàng)建，有效提高了分類的精度和時(shí)間效益。影像土地覆被總體分類精度為85.8%，kappa系數(shù)為0.84。

2）深入挖掘水利要素的各種特征，包括光譜、形狀、紋理、空間及上下文關(guān)系特征等，建立各類水利要素的提取規(guī)則，能夠有效實(shí)現(xiàn)水利要素的識(shí)別。各種水利要素的提取精度均能達(dá)到85% 以上，疊合精度達(dá)到80%以上。

3）研究結(jié)果可為構(gòu)建基于高分辨率影像的水利要素提取方法和規(guī)則集、利用高分遙感影像提取和更新水利要素等提供技術(shù)支持。從提取結(jié)果來看，河流和湖泊等面積較大、表征較為明顯的水體提取精度較高；坑塘和河流阻斷物等地物因其面積小、光譜特征不明顯，提取精度相對較低。河流阻斷物是在研究中針對水閘提出來的包含水閘但不限于水閘的橫跨于河流、渠道上的人工建筑物，雖可通過空間位置關(guān)系進(jìn)行提取，但難以區(qū)分其中的水閘、橋梁等要素。因此，更深入地挖掘光譜、位置關(guān)系等信息，建立更復(fù)雜、有效的規(guī)則集，對各種小型水利要素實(shí)現(xiàn)更細(xì)致、更高精度的提取，是下一步研究需要關(guān)注的重點(diǎn)。

[1] 陸家駒,李士鴻.TM資料水體識(shí)別技術(shù)的改進(jìn)[J].環(huán)境遙感,1992,7(1)∶17-23

[2] 肖乾廣,陳維英,王葳.氣象衛(wèi)星影像用于松花江洪水監(jiān)測[J].遙感信息,1987(4)∶26-27

[3] 楊忠恩,駱劍承,徐鵬煒.利用NAOA-AHVRR資料提取水體信息的初步研究[J].遙感技術(shù)與應(yīng)用,1995,10(1)∶25-29

[4] 都金康,黃永勝,馮學(xué)智,等.SPOT衛(wèi)星影像的水體提取方法及分類研究[J].遙感學(xué)報(bào),2001,5(3)∶214-219

[5] 張明華.用改進(jìn)的譜間關(guān)系模型提取極高山地區(qū)水體信息[J].地理與地理信息科學(xué),2008,24(2)∶14-16

[6] 鄧勁松,王坷,李君,等.決策樹方法從SPOT5衛(wèi)星影像中自動(dòng)提取水體信息研究[J].浙江大學(xué)學(xué)報(bào)(農(nóng)業(yè)與生命科學(xué)版),2005,31(2)∶171-174

[7] 楊樹文,薛重生,劉濤,等.一種利用TM影像自動(dòng)提取細(xì)小水體的方法[J].測繪學(xué)報(bào),2010,39(6)∶611-617

[8] 閆霈,張友靜,張?jiān)?利用增強(qiáng)型水體指數(shù)(EWI)和GIS去噪音技術(shù)提取半干旱地區(qū)水體信息的研究[J].遙感信息,2007(6)∶62-67

[9] 駱劍承,盛永偉,沈占鋒,等.分步迭代的多光譜遙感水體信息高精度自動(dòng)提取[J].遙感學(xué)報(bào),2009,13(4)∶610-615

[10] 路云閣,劉采,王姣.基于國產(chǎn)衛(wèi)星數(shù)據(jù)的礦山遙感監(jiān)測一體化解決方案∶以西藏自治區(qū)為例[J].國土資源遙感,2014,26(4)∶85-90

[11] 白照廣.高分一號(hào)衛(wèi)星的技術(shù)特點(diǎn)[J].中國航天,2013,36(8)∶5-9

P237

B

1672-4623（2016）05-0071-04

10.3969/j.issn.1672-4623.2016.05.023

許瑩，碩士，主要從事遙感技術(shù)機(jī)理與應(yīng)用研究。

2015-04-29。

項(xiàng)目來源：高分辨率對地觀測系統(tǒng)重大專項(xiàng)資助項(xiàng)目（08-Y30B07-9001-13/15）。