王 祥,程先富,夏 超
(1.安徽師范大學(xué) 國土資源與旅游學(xué)院,安徽 蕪湖 241003;2.安徽自然災(zāi)害過程與防控研究省級實(shí)驗(yàn)室,安徽 蕪湖 241003;3.南京農(nóng)業(yè)大學(xué) 資源與環(huán)境科學(xué)學(xué)院,江蘇 南京 210095)
基于高分一號(hào)影像的干旱半干旱區(qū)地物信息提取方法研究
王 祥1,2,程先富1,2,夏 超3
(1.安徽師范大學(xué) 國土資源與旅游學(xué)院,安徽 蕪湖 241003;2.安徽自然災(zāi)害過程與防控研究省級實(shí)驗(yàn)室,安徽 蕪湖 241003;3.南京農(nóng)業(yè)大學(xué) 資源與環(huán)境科學(xué)學(xué)院,江蘇 南京 210095)
高分一號(hào)影像是由我國自主研發(fā)的高分辨率對地觀測系統(tǒng)的首發(fā)星所獲取的影像?;诟叻忠惶?hào)影像,采用基于像元和面向?qū)ο蟮姆诸惙椒ú⒔Y(jié)合野外調(diào)查和目視解譯對干旱半干旱區(qū)高分一號(hào)影像的分類方法進(jìn)行研究。以鄂爾多斯市烏審旗段為例,對區(qū)域內(nèi)地物信息進(jìn)行提取,經(jīng)檢驗(yàn)表明將面向?qū)ο蟮姆诸惙椒ㄍ恳暯庾g相結(jié)合,總體分類精度能到達(dá)81%,基本滿足應(yīng)用要求。
高分一號(hào)影像;干旱半干旱區(qū);面向?qū)ο?;像元;分類方?/p>
鄂爾多斯盆地蘊(yùn)藏著豐富的能源資源和礦產(chǎn)資源,但該地地處干旱半干旱區(qū),水資源短缺,生態(tài)環(huán)境十分脆弱。利用國產(chǎn)高分一號(hào)影像在室內(nèi)進(jìn)行地物信息提取,能及時(shí)了解土地利用現(xiàn)狀和發(fā)展變化趨勢,尤其是對該地區(qū)植被分布狀況的掌握,對于客觀評價(jià)該地區(qū)的生態(tài)現(xiàn)狀,采取相對應(yīng)的治理措施,從而為實(shí)現(xiàn)基于生態(tài)環(huán)境保護(hù)的地區(qū)開發(fā)理念提供可靠的技術(shù)支撐。
以往的研究中國內(nèi)常采用國外的MODIS、Landsat、SPOT等遙感影像進(jìn)行地物信息提取[1],左玉珊等通過對MODIS數(shù)據(jù)波段光譜知識(shí)的充分挖掘[2],建立CART分類決策樹,實(shí)現(xiàn)了區(qū)域范圍內(nèi)土地覆被的自動(dòng)分類;胡榮明等以SPORT5高空間分辨率遙感影像為研究數(shù)據(jù)[3],采用面向?qū)ο蟮倪b感影像分類方法,將基于像元的分類精度與基于面向?qū)ο蟮姆诸惥冗M(jìn)行對比分析,結(jié)果表明面向?qū)ο蟮姆诸惙椒ū然谙裨姆诸惙椒ǚ诸惥雀?。國外遙感影像獲取周期長,高分辨率影像成本高等因素制約著研究中影像的獲取和應(yīng)用[4]。近些年,隨著國產(chǎn)高分辨率衛(wèi)星傳感器資源一號(hào)02C和高分一號(hào)等陸續(xù)發(fā)射,國產(chǎn)遙感影像在地物信息提取中得到越來越多的關(guān)注和應(yīng)用,牛增懿等利用高分一號(hào)影像[5],采用面向?qū)ο蟮姆诸惙椒▽ρ芯繀^(qū)的土壤鹽漬化信息進(jìn)行了提取,結(jié)果表明高分一號(hào)影像在土壤鹽漬化信息提取的應(yīng)用中表現(xiàn)出較好的效果。
在干旱半干旱區(qū)利用遙感影像進(jìn)行地物信息提取中,早期張冬梅等研究了一種輔助草原遙感影像目視解譯的半透明疊加方法對TM影像大尺度內(nèi)的地物進(jìn)行了目視解譯[6],提高了目視解譯的精度。隨著遙感解譯技術(shù)的發(fā)展和影像數(shù)據(jù)的不斷豐富,出現(xiàn)了多種計(jì)算機(jī)解譯算法和研究針對某一種影像數(shù)據(jù)在某一特定區(qū)域中應(yīng)用的解譯方法,張興余等利用ALOS、SPOT5高分辨率遙感影像[7],結(jié)合野外調(diào)查,完成了烏蘭布和沙漠研究區(qū)梭梭林的目視解譯及植被制圖,解譯精度達(dá)到93.3%;齊紅超等采用機(jī)器學(xué)習(xí)C5.0決策樹算法[8],綜合利用地物波譜、NDVI、TC、紋理等信息,對研究區(qū)進(jìn)行了地物信息提取,提取精度可以滿足西北干旱區(qū)大面積的土地覆被變化制圖的需要;朱海濤等構(gòu)建了基于多特征的面向?qū)ο鬀Q策樹分類方法[9],對半干旱區(qū)地物信息進(jìn)行了提取,提取精度達(dá)86.9%,滿足研究的需要。
以上研究中在數(shù)據(jù)層面上主要以國外影像為數(shù)據(jù)源;從區(qū)域研究尺度層面上在干旱半干旱區(qū)以大尺度研究為主;從研究方法層面上方法種類多,差別大;然而總體上采用國產(chǎn)高分一號(hào)影像為數(shù)據(jù)源,在中小尺度區(qū)域上對干旱半干旱區(qū)地物信息進(jìn)行提取的研究很少。
1.1 研究區(qū)概況
研究區(qū)選在內(nèi)蒙古自治區(qū)鄂爾多斯市烏審旗境內(nèi),屬于海流兔河流域。該流域位于毛烏素沙漠與陜北黃土高原的過渡帶,地理上介于108°35′~109°20′E和38°00′~39°00′N之間,海拔高度在982 m到1 480 m之間。氣候上屬溫帶大陸性季風(fēng)氣候,全年多風(fēng)少雨,形成了典型的干旱、半干旱草原氣候環(huán)境。
圖1 研究區(qū)位置示意圖
研究區(qū)位于海流兔河流域?yàn)鯇徠於?,面積約421 km2,地貌和土壤類型主要分為湖濱灘地,土壤類型為草甸土;沙地,土壤類型為風(fēng)沙土;沙梁,土壤類型為栗鈣土。海拔在1 300 m左右。地物信息主要有林地(旱柳、沙柳和檸條)、草地(沙蒿、牧草和低矮沙植)、耕地、水域、居民地、道路和其他用地(沙地和鹽堿地)共12個(gè)地物類型。
1.2 數(shù)據(jù)來源
高分一號(hào)衛(wèi)星是中國高分辨率對地觀測系統(tǒng)的首發(fā)星,其遙感影像全色波段的空間分辨率為2 m,多光譜波段的最大空間分辨率達(dá)到8 m。多光譜波段包括近紅外、紅色、綠色和藍(lán)色4個(gè)波段信息。本文所使用的高分一號(hào)遙感影像是由西安地質(zhì)調(diào)查中心提供的,成像于2014年6月26日,影像為多景全色影像和多景多光譜影像。對兩種影像進(jìn)行輻射定標(biāo)、大氣校正和正射校正后,將全色影像與多光譜影像融合,進(jìn)行影像拼接,最后裁剪出研究區(qū)的影像。
1.3 研究方法
圖2 技術(shù)流程圖
1.3.1 野外調(diào)查及影像特征判讀
通過查看文獻(xiàn),發(fā)現(xiàn)以往對此區(qū)域研究中多采用中低分辨率遙感影像作為研究的數(shù)據(jù)源,土地利用信息分類也沒有本研究精細(xì)。且此高分一號(hào)遙感影像顏色上非真彩色,主要有兩種顏色,淺黑色和淺黃色。通過以往的歷史資料和對于新入手此研究的新手來說現(xiàn)有的作業(yè)經(jīng)驗(yàn)無法對高分一號(hào)遙感影像上的土地利用信息做初步的判讀。進(jìn)行野外調(diào)查獲得研究區(qū)實(shí)地的地物信息,對研究區(qū)的地物類型,長勢,分布等有大致的了解,不僅可以提高遙感影像土地利用信息提取人的知識(shí)儲(chǔ)備,建立解譯標(biāo)志庫,且對計(jì)算機(jī)解譯后的手動(dòng)修改和人工目視解譯有很大幫助,還可以提供后期分類精度評價(jià)的實(shí)測數(shù)據(jù)[10-11]。
1.3.2 基于像元分類
長期以來,傳統(tǒng)的基于像元的分類方法是進(jìn)行遙感影像信息提取的基本方法,其方法相對于面向?qū)ο蠓诸惙椒ú僮骱啽?,對分類人員的技術(shù)要求較低,適合初學(xué)者和對分類精度要求一般的情況下使用[12]。在對研究區(qū)高分一號(hào)遙感影像土地利用信息進(jìn)行提取時(shí),首先考慮采用操作簡便的基于像元的分類方法。本研究中采用基于像元分類中的監(jiān)督分類方法,運(yùn)用常規(guī)的最大似然分類算法進(jìn)行分類。最大似然分類算法是遙感影像監(jiān)督分類的經(jīng)典算法,在土地利用信息提取中有廣泛的應(yīng)用,很多相關(guān)研究及文獻(xiàn)都把它作為標(biāo)準(zhǔn)來衡量其他分類方法的優(yōu)劣[13]。因?yàn)樵谝巴庹{(diào)查之前已經(jīng)將道路提取出來,因此這里分類類別為11個(gè)地物,利用ENVI5.2軟件進(jìn)行第一次分類。通過野外調(diào)查和對影像的判讀知道,旱柳、沙柳、耕地和居民地分布面積小,布局零散,去除這4種地物,進(jìn)行第二次分類。
表1 研究區(qū)地物影像特征及預(yù)采取的分類方法
1.3.3 面向?qū)ο蠓诸?/p>
面向?qū)ο蠓诸惙椒ㄏ鄬τ诨谙裨姆诸惙椒ㄔ诓僮魃蠈Ψ诸惾藛T的技術(shù)要求較高,其構(gòu)建分類規(guī)則庫是個(gè)不斷試驗(yàn)和摸索的過程,花費(fèi)時(shí)間較長[14]。隨著高分辨率遙感影像在生產(chǎn)和科學(xué)研究中得到廣泛的運(yùn)用,新的遙感影像分類方法面向?qū)ο蠓诸惙椒ㄒ蚱溥m合高精度遙感影像分類且分類精度較高隨之新興起來,越來越多的人開始學(xué)習(xí)和使用面向?qū)ο蠓诸惣夹g(shù)[15]。在采用傳統(tǒng)的基于像元的分類方法中分類結(jié)果不夠理想,旱柳、沙柳、耕地和住宅用地不能提取,各地物分類精度差別很大,總體分類精度不夠高,需要采用面向?qū)ο蠓诸惙椒ㄟM(jìn)行分類。
面向?qū)ο蠓诸惔笾掳ǘ喑叨确指?,特征信息提取和選用一種分類方法進(jìn)行分類這三個(gè)步驟[16]。本研究中利用eCognition軟件在參考前人的分類方法和相應(yīng)地物的分類參數(shù)后根據(jù)本研究區(qū)地物的光譜特征、形狀特征、紋理特征和鄰里關(guān)系等影像特性,反復(fù)試驗(yàn)摸索出適合本研究區(qū)地物信息提取的具體分類參數(shù),建立適合本研究區(qū)土地利用信息提取的分類規(guī)則庫,采用面向?qū)ο蠓诸惖某R?guī)分類方法最鄰近距離法進(jìn)行分類。在分類中,反復(fù)試驗(yàn)發(fā)現(xiàn)去除旱柳和沙柳后能有效提高其他地物的分類精度,因此分類類別為9個(gè)地物,分別為檸條、沙蒿、牧草、低矮沙植、耕地、水域、住宅用地、沙地和鹽堿地。
1.3.4 目視解譯
目視解譯是許多遙感應(yīng)用項(xiàng)目的先遣工作,是遙感應(yīng)用的基礎(chǔ)。對已提取的地物信息進(jìn)行修改與調(diào)整,以提高最終結(jié)果的準(zhǔn)確度并根據(jù)要求進(jìn)行制圖。在采用以上兩種方法提取地物信息后,剩下旱柳和沙柳未提取。根據(jù)先前對遙感影像的熟悉,已經(jīng)能很好地知道旱柳和沙柳分布的一些基本規(guī)律。旱柳多分布在居民點(diǎn)、耕地和部分水域附近,分布相對于沙柳在小區(qū)域上較集中,外形上呈塊狀,顏色相對較深,紋理粗糙。沙柳絕大部分分布在大片沙地區(qū),很少與其他植被混雜到一塊,區(qū)域分布廣,但面積小,外形上呈條帶狀由東北朝西南走向,顏色相對較淺。根據(jù)這些影像特性解譯出旱柳和沙柳。
如圖3所示在基于像元分類中當(dāng)分類類別越多,類別越細(xì)小時(shí),分類的結(jié)果中出現(xiàn)的椒鹽現(xiàn)象越嚴(yán)重。而面向?qū)ο蟮姆诸惙椒梢酝ㄟ^多尺度分割將多個(gè)像元看做一個(gè)整體來克服圖斑零碎的問題,使圖文效果更好。對于圖中分布面積不大的旱柳、居民地和耕地,可以根據(jù)其分布規(guī)律采用人工目視解譯的方法進(jìn)行提取,這樣在實(shí)際解譯工作中比其他分類方法要更實(shí)用。
圖3 不同分類方法分類結(jié)果圖
次數(shù)各類別分類精度/%旱柳沙柳檸條沙蒿牧草低矮沙植耕地水域居民地沙地鹽堿地總體精度/%一673741182411738556827二//73696835/81/837162
表2中/表示在該分類方法中沒有對此種地物類型進(jìn)行分類。從表2中可以看出第一次分類總體精度只有27%,遠(yuǎn)低于第二次分類的精度。其原因可能有兩方面,一方面,分類類別過多增加了具有相似光譜特征的地物區(qū)分度;類別中耕地、低矮沙植同沙地光譜相似度很高,部分牧草、鹽堿地、水域同居民地的相似度很高。另一方面,基于像元的分類方法主要是基于光譜信息特征,該方法大多用于區(qū)域內(nèi)較大地物類型區(qū)分類;類別中旱柳、沙柳、耕地和居民地4種地物在高分一號(hào)遙感影像中分布零散,總體面積小,屬于較小零碎地物,難以與其他地物區(qū)分。當(dāng)去除這4種地物后進(jìn)行第二次分類后各種地物的分類精度明顯提高。第一次分類和第二次分類中的水域和鹽堿地分類精度差別不大,兩次分類精度都較高,可能因?yàn)閮煞N地物光譜顏色具有顯著的特征,能很好地同其他地物區(qū)分出來。
表3 基于像元分類和面向?qū)ο蠓诸惥仍u價(jià)
表3是對基于像元分類(7種地物)和面向?qū)ο蠓诸?9種地物)兩種分類方法的分類精度評價(jià)統(tǒng)計(jì),表中/表示在此種分類方法中沒有對此種地物類型進(jìn)行分類。從表3中可以看出基于像元分類方法對水域和沙地的分類精度較高,其中沙地達(dá)到83%,此分類精度達(dá)到實(shí)際應(yīng)用要求。面向?qū)ο蠓诸惙椒ㄔ诒然谙裨诸惙椒ǘ喑鰞煞N地物的情況下總體精度比基于像元的分類精度還高出了19個(gè)百分點(diǎn),且其植被部分分類精度較高,除了居民地分類精度較差外,其余地物類型經(jīng)過手動(dòng)修改后基本能達(dá)到應(yīng)用要求。面向?qū)ο蟮姆诸惙椒ㄔ诜诸惥壬厦黠@高于基于像元的分類方法。
本文采用了基于像元和面向?qū)ο蟮姆诸惙椒ú⒔Y(jié)合野外調(diào)查和目視解譯對干旱半干旱區(qū)高分一號(hào)影像的分類方法進(jìn)行了研究,且將基于像元的分類方法同面向?qū)ο蟮姆诸惙椒ㄏ啾容^,得出以下結(jié)論:
(1)干旱半干旱區(qū)地物中植被的種類較少,分布的區(qū)域和面積差別較大,在高分一號(hào)影像上分布面積較小的旱柳和沙柳不易與其他地物區(qū)分。
(2)面向?qū)ο蟮姆诸惙椒ㄏ啾容^基于像元的分類方法在形狀因子和紋理因子上更能區(qū)分出地物,分類精度更高,并能克服圖像的椒鹽現(xiàn)象。
(3)高分一號(hào)影像在中小尺度上對干旱半干旱區(qū)地物信息進(jìn)行詳細(xì)提取的應(yīng)用中表現(xiàn)出較好的效果,結(jié)合目視解譯,其分類精度能達(dá)到81%,基本滿足應(yīng)用要求。
[1]初慶偉,張洪群,吳業(yè)煒,等.Landsat-8衛(wèi)星數(shù)據(jù)應(yīng)用探討[J].遙感信息,2013,28(4):110-114.
[2]左玉珊,王衛(wèi),郝彥莉,等.基于MODIS影像的土地覆被分類研究——以京津冀地區(qū)為例[J].地理科學(xué)進(jìn)展,2014,33(11):1556-1565.
[3]胡榮明,魏曼,楊成斌,等.以SPORT5遙感數(shù)據(jù)為例比較基于像素與面向?qū)ο蟮姆诸惙椒╗J].遙感技術(shù)與應(yīng)用,2012,27(3):366-371.
[4]郭燕,武喜紅,程永政,等.用高分一號(hào)數(shù)據(jù)提取玉米面積及精度分析[J].遙感信息,2015,30(6):31-36.
[5]牛增懿,丁建麗,李艷華,等.基于高分一號(hào)影像的土壤鹽漬化信息提取方法[J].干旱區(qū)地理,2016,39(1):171-181.
[6]張冬梅,袁清,王加亭.一種輔助草原遙感影像目視解譯方法的探討——影像與基礎(chǔ)圖件半透明疊加法[J].遙感技術(shù)與應(yīng)用,2006,21(6):560-564.
[7]張興余,劉勇,許寶榮,等.烏蘭布和沙漠高分辨率遙感影像梭梭林解譯方法探討[J].遙感技術(shù)與應(yīng)用,2010,25(6):828-835.
[8]齊紅超,祁元,徐瑱.基于C5.0決策樹算法的西北干旱區(qū)土地覆蓋分類研究——以甘肅省武威市為例[J].遙感技術(shù)與應(yīng)用,2009,24(5):648-653.
[9]朱海濤,張霞,王樹東,等.基于面向?qū)ο鬀Q策樹算法的半干旱地區(qū)遙感影像分類[J].遙感信息,2013,28(4):50-56.
[10]杜佳虹.遙感影像目視解譯矢量化分析[J].測繪與空間地理信息,2011,34(2):87-91.
[11]楊博,劉義.農(nóng)業(yè)遙感影像目視解譯技術(shù)要點(diǎn)[J].現(xiàn)代化農(nóng)業(yè),2008,(4):37-39.
[12]錢茹茹.遙感影像分類方法比較研究[D].西安:長安大學(xué)(碩士學(xué)位論文),2007.
[13]李全,王海燕,李霖.基于最大似然分類算法的土地覆蓋分類精度控制研究[J].國土資源科技管理2005,(22):42-45.
[14]Trias-Sanz R,Stamon G,Louchet J.Using Colour,Texture,and Hierarchial Segmentation for High-resolution Remote Sensing[J].ISPRS Journal of Photogrammetry and Remote Sensing,2008,63(2):156-168.
[15] Kurtz C,Passat N,Gancarski P,et al.Extractionof Complex Patterns from Multi-resolution Remote Sensing Images:A Hierarchical Top-down Methodology[J].Pattern Recognition,2012,45(2):685-706.
[16] Whiteside T G,Boggs G S,Maier S W. Comparing Objectbased and Pixel-based Classifications for Mapping Savannas[J].International Journal of Applied Earth Observation and Geoinformation,2011,13(6):884-893.
Research on the Methods of Feature Information Extraction in Aridand Semi-arid Areas from the GF-1 Images
WANG Xiang1,2, CHENG Xianfu1,2, XIA Chao3
(1,SchoolofTerritorialResourcesandTourism,AnhuiNormalUniversity,Wuhu241003,China;2.AnhuiKeyLaboratoryofNaturalDisasterProcessandPrevention,Wuhu241003,China;3.Resources&EnvironmentScienceDepartment,NanjingAgriculturalUniversity,Nanjing210095,China)
GF-1 images by our independent research and development of high resolution earth observation system of star image can be obtained from the start. Based on the GF-1 images, pixel-based and object-oriented classification method and combined with field investigation and the visual interpretation, to arid and Semi-arid areas happened in GF-1 images classification method were studied. Our study area is Wushen County in Ordos city,to extract the mainland area information. The test shows that the object-oriented classification method combined with visual interpretation, the overall classification accuracy can reach 81%, basic meet application requirements.
GF-1 images; arid and Semi-arid areas; object-oriented; pixel; method of classification
2016-12-01
王祥(1990-),男,安徽蕪湖人,碩士研究生,研究方向:區(qū)域環(huán)境與遙感應(yīng)用。
P237
A
1009-9735(2017)02-0104-06