林 聰，李二珠，杜培軍

(1. 南京大學(xué)衛(wèi)星測繪技術(shù)與應(yīng)用國家測繪地理信息局重點實驗室，江蘇 南京 210023； 2. 南京大學(xué)江蘇省地理信息技術(shù)重點實驗室， 江蘇 南京 210023)

土地利用/覆蓋變化對陸地生態(tài)系統(tǒng)的生物多樣性，水、碳和養(yǎng)分循環(huán)，能量平衡，以及溫室氣體釋放增加等其他環(huán)境問題具有重要的影響[1]。近年來，隨著遙感技術(shù)的發(fā)展，基于遙感數(shù)據(jù)的解譯方法成為獲取土地利用/覆蓋信息的主要方式之一。遙感影像分類一直是遙感影像處理的重要環(huán)節(jié)[2]。隨著計算機(jī)技術(shù)的不斷發(fā)展，神經(jīng)網(wǎng)絡(luò)分類、專家系統(tǒng)分類、支持向量機(jī)分類等機(jī)器學(xué)習(xí)方法在土地覆蓋遙感分類中取得了廣泛的應(yīng)用[3]。分類方法的發(fā)展極大地滿足了日益龐大的遙感信息提取所需要的數(shù)量與速度，但是樣本的選取成為制約遙感圖像自動解譯發(fā)展的一個重要因素。由于受大氣吸收與散射、物候、傳感器定標(biāo)等因素的影響，同一地物的光譜在兩期影像上不能服從同樣的概率統(tǒng)計分布[4]。而針對大范圍的分類制圖任務(wù)，重新標(biāo)注樣本需要耗費(fèi)大量的時間與人力，制約了遙感數(shù)據(jù)的自動化處理及長時間序列土地利用/覆蓋變化的研究。

近年來，遷移學(xué)習(xí)已經(jīng)引起了廣泛關(guān)注和研究。遷移學(xué)習(xí)是運(yùn)用已有知識對不同但相關(guān)領(lǐng)域問題進(jìn)行求解的一種機(jī)器學(xué)習(xí)方法，它放寬了傳統(tǒng)機(jī)器學(xué)習(xí)中的兩個基本假設(shè)：①用于學(xué)習(xí)的訓(xùn)練樣本與新的測試樣本滿足獨立同分布的條件；②必須用足夠的訓(xùn)練樣本才能學(xué)習(xí)得到一個好的分類模型[5-6]。遷移學(xué)習(xí)是為了解決已有信息、知識的再利用問題而提出的一種機(jī)器學(xué)習(xí)新理論，這一理論為解決上述遙感分類中的樣本選取問題提供了新思路。遷移學(xué)習(xí)在遙感圖像的目標(biāo)識別[7-8]及土地覆蓋分類[4，9-13]中已經(jīng)有了一些研究，但研究多針對高分?jǐn)?shù)據(jù)或高光譜數(shù)據(jù)，而對于在土地利用/覆蓋研究與應(yīng)用中最常用的中等分辨率數(shù)據(jù)如Landsat衛(wèi)星數(shù)據(jù)的研究較匱乏。本文基于Landsat數(shù)據(jù)，從像元級分類發(fā)出，充分挖掘其豐富的光譜信息，并從知識遷移的角度，一定程度上解決了訓(xùn)練樣本在目標(biāo)影像中的優(yōu)化問題，從而構(gòu)建一種基于樣本標(biāo)簽遷移的遙感數(shù)據(jù)自動分類方法。

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

巫山縣位于重慶市東部，處三峽庫區(qū)腹心，素有“渝東北門戶”之稱。巫山縣位于北緯30°45′—31°29′，東經(jīng)109°33′—110°12′之間，轄區(qū)面積2958 km2，地處四川盆地東端，屬四川盆地東部山地地貌，相對高差懸殊。研究區(qū)屬中亞熱帶濕潤季風(fēng)氣候區(qū)，四季分明，無霜期長，雨量充沛，日照時間長[14]。本文選取研究區(qū)位置如圖1所示，位于巫山縣內(nèi)，主要包括縣政府所在的巫峽鎮(zhèn)，以及龍鄉(xiāng)鎮(zhèn)、兩坪鄉(xiāng)、建平鄉(xiāng)大部分地區(qū)和其他鄉(xiāng)鎮(zhèn)的部分地區(qū)，面積約為441 km2。

圖1 研究區(qū)示意圖

研究用到的數(shù)據(jù)主要包括全球30 m地表覆蓋數(shù)據(jù)(GlobeLand30)、研究區(qū)1∶5萬DEM及Landsat TM影像。GlobeLand30數(shù)據(jù)集[15-16]是在863計劃支持下，由國家基礎(chǔ)地理信息中心牽頭完成的30 m分辨率的全球地表覆蓋遙感制圖數(shù)據(jù)產(chǎn)品。該數(shù)據(jù)覆蓋全球南北緯80°內(nèi)的陸地范圍，包括10種地表覆蓋類型，其中研究區(qū)范圍內(nèi)共有耕地、森林、草地、水體、人造地表5類覆蓋類型。如圖2、圖3所示。

圖2 研究區(qū)DEM與Globeland30產(chǎn)品

研究所使用的TM數(shù)據(jù)共兩期(見表1)，作為參考土地覆蓋分類專題圖的GlobeLand30數(shù)據(jù)雖然以2010年為基準(zhǔn)年，但實際用于解譯的影像大多獲取時間為2009—2011年，因此為了充分檢測變化信息，提高最終結(jié)果的精度，選取2011年而非2010年的數(shù)據(jù)作為源領(lǐng)域影像。

表1 TM數(shù)據(jù)信息

2 研究方法

遷移學(xué)習(xí)是基于給定的源領(lǐng)域數(shù)據(jù)集Ds與學(xué)習(xí)任務(wù)Ts，充分挖掘其中的知識以幫助完成或完善Dt中的學(xué)習(xí)任務(wù)Tt的機(jī)器學(xué)習(xí)方法[5]。近年來，隨著時間序列遙感研究的發(fā)展，遷移學(xué)習(xí)逐漸被引入遙感自動分類方法中，如面向?qū)ο蠓诸愔械膶ο蠹墭颖具w移[4]與場景分類中的特征遷移[12]。以Landsat數(shù)據(jù)為代表的中等分辨率的遙感數(shù)據(jù)，可以對全球大多數(shù)區(qū)域長時間連續(xù)觀測，并且有可以提供樣本標(biāo)簽的相關(guān)土地利用/覆蓋產(chǎn)品，其樣本遷移的問題有更大的應(yīng)用研究價值。本文運(yùn)用變化檢測技術(shù)與光譜形狀編碼的方法，構(gòu)建了地物類別標(biāo)簽(關(guān)聯(lián)知識)遷移的自動分類框架，如圖4所示。

圖3 研究區(qū)遙感影像

圖4 基于樣本遷移的自動分類算法框架

2.1 遷移學(xué)習(xí)

衛(wèi)星對同一地區(qū)的重復(fù)觀測可以獲得時間序列遙感數(shù)據(jù)，假設(shè)X={X1,X2,…,XM}為M期同一地區(qū)但獲取時間不同的多光譜影像，其中第m期為有K個波段的多光譜數(shù)據(jù)Xm={xm,1,xm,2,…,xm,N}(其中N為研究區(qū)影像的像元總數(shù))，其中含有大量訓(xùn)練樣本，而其他的影像僅有少量不足以分類甚至沒有任何樣本。通過遷移學(xué)習(xí)的方法可以獲得時間序列內(nèi)其他任意影像的土地覆蓋分類圖。當(dāng)目標(biāo)影像中有少量標(biāo)注樣本時，這種遷移的方法被稱為歸納遷移學(xué)習(xí)(inductive transfer learning)，當(dāng)目標(biāo)影像中沒有任何標(biāo)簽樣本時，則為直推式遷移學(xué)習(xí)(transductive transfer learning )。

2.2 變化檢測

變化檢測是上述算法流程中的一個關(guān)鍵技術(shù)，其主要作用是確定X1和X2兩個時相影像之間的不變信息。相對輻射定標(biāo)可以提升變化檢測效果，在變化檢測之前可以考慮采用暗像元法或者直方圖匹配法來完成相對輻射定標(biāo)。采用變化矢量分析法(CVA)，變化矢量是描述從時相1到時相2某像元光譜變化大小和方向的矢量[17]。變化矢量分析法首先需要對兩個不同時相的影像作差值運(yùn)算，生成多光譜的差值圖像。差值圖像中的第i個像元的變化矢量可以表示為[16]

(1)

(2)

(3)

如上文所述，ρ表示第i個像元的變化強(qiáng)度，在本文的變化檢測中，僅考慮運(yùn)用CVA來確定像元的變化強(qiáng)度從而選擇不變區(qū)域，不考慮變化類型的問題。

對于CVA計算得到的變化檢測結(jié)果，采用大津法(Otsu’s)對結(jié)果進(jìn)行閾值分割。Otsu’s是一種圖像灰度自適應(yīng)的閾值分割算法，該方法按照圖像上的灰度值分布，將圖像分為背景和前景兩部分，遍歷不同的閾值，計算不同閾值下對應(yīng)的背景與前景之間的類內(nèi)方差取得最大值時，對應(yīng)的閾值即為所求的閾值。

通過上述過程，確定不變地物像元對應(yīng)的空間位置，并將相應(yīng)的地類標(biāo)簽標(biāo)注到目標(biāo)影像上。

2.3 光譜形狀計算

對于上文得到的結(jié)果，可以采用隨機(jī)選取的方法抽取樣本對目標(biāo)影像進(jìn)行分類，但是考慮到土地覆蓋分類產(chǎn)品C1本身的精度及分類產(chǎn)品，在直接分類的結(jié)果上進(jìn)行了一系列后處理，且分類產(chǎn)品使用影像與源影像存在一些差異，導(dǎo)致實際上遷移至目標(biāo)影像中的地類標(biāo)簽可能不完全準(zhǔn)確。為了選擇可靠的樣本進(jìn)行訓(xùn)練，采用光譜形狀計算的方法構(gòu)建一種較優(yōu)樣本的選擇策略。

對于中等分辨率的多光譜數(shù)據(jù)，很難獲取到較好的形狀或紋理特征，尤其是像研究區(qū)這樣地形起伏較大、景觀格局破碎度較高的區(qū)域。同時Landsat數(shù)據(jù)有較豐富的光譜信息，重復(fù)挖掘其光譜信息有助于優(yōu)化分類結(jié)果。通過對影像的觀察統(tǒng)計發(fā)現(xiàn)，同類地物如林地由于太陽高度角或生長程度的不同，表現(xiàn)在影像上的各個波段的反射率可能有差異。但是光譜曲線的形狀總是一致或相似，可以通過一種光譜曲線形狀計算的方法獲得[20]。

(4)

(5)

一般認(rèn)為，訓(xùn)練樣本與分類結(jié)果在特征空間上分布越趨于一致，取得的分類效果就越理想，因此直方圖統(tǒng)計常用樣本的選取。利用不變地物的土地覆蓋分類結(jié)果將地類標(biāo)記在目標(biāo)影像上。對于其中的每一類樣本，采用上述公式可以獲得其相應(yīng)的光譜形狀編碼結(jié)果。將計算得到相同的光譜形狀的樣本進(jìn)行聚類，并統(tǒng)計聚類后的直方圖，去掉直方圖內(nèi)值較小的類，并按直方圖內(nèi)的比例選取訓(xùn)練樣本。

3 試驗與結(jié)果分析

研究區(qū)影像尺寸為700×700像素，以2011年獲取的Landsat影像為源領(lǐng)域影像，對應(yīng)的土地覆蓋分類產(chǎn)品為GlobeLand30，待分類的目標(biāo)領(lǐng)域影像為2007年獲取的Landsat影像，采用本文提出的樣本遷移方法，對目標(biāo)數(shù)據(jù)進(jìn)行分類。為了更好地說明光譜曲線形狀在樣本選擇中的作用，測試了不使用光譜曲線形狀下的分類效果；同時考慮研究區(qū)地形起伏較大的特點，加入數(shù)字高程模型與坡度作為特征，發(fā)現(xiàn)對分類結(jié)果有積極的影響。選取支持向量機(jī)(SVM)作為目標(biāo)影像分類的分類器。

分類結(jié)果如圖5所示，其中圖5(a)為未經(jīng)過光譜形狀計算直接分類結(jié)果(SVM)，圖5(b)為光譜形狀計算后的分類結(jié)果(SVM_scs)，圖5(c)(SVM(T))與圖5(d)(SVM(T)_scs)為圖5(a)與圖5(b)基礎(chǔ)上加上地形特征后的分類結(jié)果。結(jié)合2006年的研究區(qū)二調(diào)數(shù)據(jù)及Google Earth高清影像，隨機(jī)選取了2109個有效的驗證樣本點，對分類結(jié)果進(jìn)行精度驗證，同時計算了4次分類結(jié)果的各類地物的生產(chǎn)者精度、總體精度與Kappa系數(shù)，見表2。

從精度驗證的結(jié)果來看，加入地形特征對研究區(qū)的分類結(jié)果的精度有顯著提高，主要表現(xiàn)在對人工地表、耕地、草地的分類精度有明顯提升。對比SVM_scs與SVM分類結(jié)果、SVM(T)_scs與SVM(T)分類結(jié)果，以及SVM(T)_scs與SVM(T)分類結(jié)果發(fā)現(xiàn)，光譜形狀計算對總體精度和Kappa系數(shù)有明顯的提高，其中主要表現(xiàn)在對耕地和草地兩類的精度有明顯的提高。考慮到耕地和草地因為植被的稀疏程度或同一地區(qū)不同作物物候特征存在差異因素，這兩類地類類內(nèi)的光譜可能就存在較大差異。光譜形狀方法對樣本的選取進(jìn)行了一定的優(yōu)化，使得不同光譜曲線形狀的耕地或草地按照各自的比例選取樣本，這樣的訓(xùn)練樣本更能代表地物的真實分布，從而對分類的結(jié)果有一定的提高。

SVM(T)_scs的分類總體精度達(dá)到了93.65%，對應(yīng)的Kappa系數(shù)為0.914 1。由此可見，本文提出的自動分類方法可以得到較好的分類結(jié)果。尤其對人工地表、水體、林地和草地的提取精度較高，均在85%以上。耕地的提取精度相對較低，為84.18%，考慮到研究區(qū)內(nèi)的耕地類型較多樣，且部分區(qū)域分布十分破碎。因此耕地的最終提取結(jié)果還算較為理想，可以在一定程度上滿足應(yīng)用需求。

4 結(jié) 語

本文結(jié)合土地覆蓋遙感分類產(chǎn)品提出了一種基于遷移學(xué)習(xí)的樣本自動選取的方法，并構(gòu)建了面向Landsat的像元級自動分類的算法框架。基于GlobeLand30產(chǎn)品，運(yùn)用基于變化矢量分析的變化檢測與光譜形狀編碼的方法，充分挖掘Landsat數(shù)據(jù)豐富的多光譜信息，通過知識遷移的方法將源領(lǐng)域內(nèi)的地類標(biāo)簽遷移至目標(biāo)影像中，并通過支持向量機(jī)算法完成分類，最終分類精度達(dá)到93.65%。綜上所述，本文所提出的樣本遷移方法大大降低了目視采集的煩瑣性，可以自動選取可靠的樣本，并且能夠得到較好的分類結(jié)果，從而提高了遙感影像分類的自動化程度，在一定程度上滿足了時間序列遙感數(shù)據(jù)的智能解譯需求，具有實際應(yīng)用價值。

參考文獻(xiàn)：

[1] 陳廣生,田漢勤.土地利用/覆蓋變化對陸地生態(tài)系統(tǒng)碳循環(huán)的影響[J].植物生態(tài)學(xué)報,2007，31(2):189-204.

[2] 沈泉飛，曹敏，史照良，等．基于布谷鳥算法的遙感影像智能分類[J]．測繪通報，2017(1): 65-68.

[3] 張增詳,汪瀟,溫慶可,等.土地資源遙感應(yīng)用研究進(jìn)展[J].遙感學(xué)報,2016,20(5):1243-1258.

[4] 吳田軍,駱劍承,夏列鋼,等.遷移學(xué)習(xí)支持下的遙感影像對象級分類樣本自動選擇方法[J].測繪學(xué)報,2014,43(9):908-916.

[5] PAN J S，YANG Q.A Survey on Transfer Learning[J].IEEE Transactions on Knowledge & Data Engineering,2010,22(10):1345-1359.

[6] 莊福振,羅平,何清,等.遷移學(xué)習(xí)研究進(jìn)展[J].軟件學(xué)報,2015,26(1):26-39.

[7] DAN Z,SANG N,CHEN Y,et al. Remote Sensing Object Recognition Based on Transfer Learning[C]∥International Conference on Fuzzy Systems and Knowledge Discovery (FSKD).Shenyang：IEEE,2013:930-934.

[8] 魯恒,付蕭,賀一楠,等.基于遷移學(xué)習(xí)的無人機(jī)影像耕地信息提取方法[J].農(nóng)業(yè)機(jī)械學(xué)報,2015,46(12):274-284.

[9] RAJAN S,GHOSH J,CRAWFORD M M.Exploiting Class Hierarchies for Knowledge Transfer in Hyperspectral Data[J].IEEE Transactions on Geoscience & Remote Sensing,2006,44(11):3408-3417.

[10] PERSELLO C,BRUZZONEL.Active Learning for Domain Adaptation in the Supervised Classification of Remote Sensing Images[J].IEEE Transactions on Geoscience & Remote Sensing,2012,50(11):4468-4483.

[11] DEMIR B,BOVOLO F，BRUZZONE L.Updating Land-Cover Maps by Classification of Image Time Series: A Novel Change-Detection-Driven Transfer Learning Approach[J].IEEE Transactions on Geoscience & Remote Sensing,2013,51(1):300-312.

[12] SUN H,LIU S,ZHOU H,et al.Unsupervised Cross-View Semantic Transfer for Remote Sensing Image Classification[J].Geoscience & Remote Sensing Letters,2016,13(1):13-17.

[13] ZHOU Y,LIAN J,HAN M.Remote Sensing Image Transfer Classification Based on Weighted Extreme Learning Machine[J].Geoscience & Remote Sensing Letters,2016,13(10):1405-1409.

[14] 羅潔瓊.基于GIS的三峽庫區(qū)山地人居環(huán)境自然適宜性動態(tài)評價——以奉節(jié)縣和巫山縣為例[D].重慶:西南大學(xué),2013.

[15] 廖安平，彭舒，武昊，等． 30 m 全球地表覆蓋遙感制圖生產(chǎn)體系與實踐[J]．測繪通報，2015(10): 4-8．

[16] 陳軍,陳晉,廖安平,等.全球30m地表覆蓋遙感制圖的總體技術(shù)[J].測繪學(xué)報,2014,43(6):551-557.

[17] 陳晉,何春陽，卓莉.基于變化向量分析(CVA)的土地利用/覆蓋變化動態(tài)監(jiān)測(Ⅱ)——變化類型的確定方法[J].遙感學(xué)報,2001,5(5):346-353.

[18] 常方正,趙銀娣,劉善磊.遙感影像CVA 變化檢測的CUDA 并行算法設(shè)計[J].遙感學(xué)報，2016,20(1):114-128.

[19] LIU S,BRUZZONE L,BOVOLO F,et al.Sequential Spectral Change Vector Analysis for Iteratively Discovering and Detecting Multiple Changes in Hyperspectral Images[J].IEEE Transactions on Geoscience Remote Sensing,2015,53(8):4363-4378.

[20] LI E,DU P,SAMAT A,et al.An Automatic Approach for Urban Land-cover Classification from Landsat-8 OLI Data[J].International Journal of Remote Sensing,2015,36(24):5983-6007.