王 璠 楊璐萍

(安徽新華學(xué)院 城市建設(shè)學(xué)院,安徽 合肥 230088)

0 引言

隨著城市化進程和社會經(jīng)濟的發(fā)展,基于遙感影像對城市擴張監(jiān)測顯得尤為重要。因此對區(qū)域尺度上的城市擴張進行有效的監(jiān)測可為土地規(guī)劃和生態(tài)環(huán)境保護等提供科學(xué)參考,具有重要的現(xiàn)實意義[1]。目前已有國內(nèi)許多學(xué)者對城市擴張進行研究,常變?nèi)?、李仁東在分析武漢市建設(shè)用地擴張的時空特征時,利用面向?qū)ο蠓诸惙椒ㄌ崛〗ㄔO(shè)用地,采用空間數(shù)據(jù)庫引擎(spatial database engine, SDE)空間分析方法,得出前期擴張向東北—西南方向演變,后期則不存在方向性的結(jié)論,為后期各種政府機關(guān)決策提供了依據(jù)[2];陳征等人以浙江省為例,利用遙感影像數(shù)據(jù)和夜間燈光數(shù)據(jù)提取城鎮(zhèn)信息,區(qū)分出城鎮(zhèn)用地和非城鎮(zhèn)用地,進行了城鎮(zhèn)時空擴張的特征分析[3];孫善磊等人以環(huán)杭州灣地區(qū)為例,使用決策樹和形狀指數(shù)相結(jié)合提取城鎮(zhèn)用地信息,證明了此方法的可行性,且杭州市向西北和北方向向外擴張[4]。而此研究則是意在探究精度,因此使用支持向量機分類算法,融合不同數(shù)據(jù)源波段,提高分類精度。詳細(xì)流程如圖1所示。

圖1 研究路線

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

安徽省坐落于中國東部,與山東、江蘇、河南、湖北等多省接壤,有省會合肥和蕪湖、馬鞍山等16個地級市。安徽省是中部經(jīng)濟地帶,接北聯(lián)南,承東啟西。安徽省位于我國東部季風(fēng)區(qū),由于淮河從中穿過,因此淮河以北屬于溫帶季風(fēng)氣候,淮河以南屬于亞熱帶季風(fēng)氣候[5]。安徽地形多樣,平原占44.7%、山地占27.8%、丘陵占27.5%。安徽省土地總面積13.96萬km2,土地資源類型多樣,土地構(gòu)成較為合理,土地利用潛力很大。

1.2 數(shù)據(jù)與預(yù)處理

研究使用了2018年安徽省的哨兵2號(Sentinel-2),美國陸地衛(wèi)星(Landsat8)和珞珈一號夜間燈光數(shù)據(jù)。Sentinel-2數(shù)據(jù)是從歐洲航天局官網(wǎng)(https://scihub.copernicus.eu/dhus/#/home)下載獲得,本研究覆蓋范圍為2018年,選取第2、3、4、8波段,空間分辨率為10 m[6-7]。Landsat8從美國地質(zhì)勘探局(United States Geological Survey,USGS)官網(wǎng)(https://earthexplorer.usgs.gov/)獲得,時間與Sentinel-2數(shù)據(jù)一致,選擇了綠光波段、紅光波段、近紅外波段和中紅外波段,空間分辨率是30 m[8]?！扮箸煲惶枴笔侨蚴最w專業(yè)夜光遙感衛(wèi)星,由武漢大學(xué)團隊與相關(guān)機構(gòu)共同研發(fā)制作,在2018年6月2日12時13分,高分六號及其搭載的珞珈一號科學(xué)試驗衛(wèi)星,在中國酒泉衛(wèi)星發(fā)射中心發(fā)射升空,衛(wèi)星進入預(yù)定軌道[9]?！扮箸煲惶枴贝钶d了高靈敏度夜光相機,其精度將達到地面分辨率100 m。數(shù)據(jù)的預(yù)處理基于GEE(Google Earth Engine)平臺,對數(shù)據(jù)進行預(yù)處理工作,包括對不同數(shù)據(jù)源進行去云、定標(biāo)、幾何校正、拼接、裁剪等[10]。

1.3 基于支持向量積分類方法

1.3.1訓(xùn)練樣區(qū)選擇

訓(xùn)練樣本選擇應(yīng)該具有代表性,并且要滿幅均勻選取。本研究選取了森林、農(nóng)田、水體、不透水面四類地物類型[11]。詳細(xì)統(tǒng)計如表1所示。

表1 訓(xùn)練樣區(qū)統(tǒng)計

1.3.2波段計算與分類

由表2可知,研究選取了Landsat8的綠光波段、紅光波段、近紅外、中紅外波段進行指數(shù)運算,計算的地物特征指數(shù)有歸一化植被指數(shù)(normalized vegetation index,NDVI),歸一化水體指數(shù)(normalized water index,NDWI),歸一化建筑物指數(shù)(normalized building index,NDBI),并選取了紅光波段和近紅外波段[12]。同時又選取了Sentinel-2數(shù)據(jù)空間分辨率為10 m的波段,分別為第3(GREEN)、4(RED)、8(NIR)波段,由于哨兵數(shù)據(jù)第11波段SWIR波段對建筑物提取效果較好,因此選取分辨率為20 m的11(SWIR)波段,計算了NDVI、NDBI、NDWI[13]。

表2 指數(shù)運算表

支持向量機分類法(support vector machine,SVM)是監(jiān)督分類的一種模型[14]。可以分析數(shù)據(jù)、識別模式,用于分類和回歸分析。支持向量機方法是建立在統(tǒng)計學(xué)習(xí)理論的VC維理論和結(jié)構(gòu)風(fēng)險最小原理基礎(chǔ)上的,根據(jù)有限的樣本信息在模型的復(fù)雜性(即對特定訓(xùn)練樣本的學(xué)習(xí)精度)和學(xué)習(xí)能力(即無錯誤地識別任意樣本的能力)之間尋求最佳折中,以期獲得最好的推廣能力。研究所使用的SVM不僅能分類出像元屬性,而且也能解決間隔值的問題[15]。

將Landsat8的RED、NIR、NDVI、NDWI、NDBI與Sentinel-2數(shù)據(jù)的RED、NIR、NDVI、NDWI、NDBI和珞珈一號數(shù)據(jù)加入SVM算法進行分類,同時將Landsat8的RED、NIR、NDVI、NDWI、NDBI與珞珈一號數(shù)據(jù)和Sentinel-2數(shù)據(jù)的RED、NIR、NDVI、NDWI、NDBI與珞珈一號數(shù)據(jù)也各做一次分類作為本研究的對比結(jié)果。

2 結(jié)果與分析

2.1 土地利用分類

在進行分類時,將Landsat8的RED、NIR、NDVI、NDWI、NDBI與Sentinel-2數(shù)據(jù)的RED、NIR、NDVI、NDWI、NDBI結(jié)合,生成一個10 m的安徽省土地利用分類圖。Landsat8的RED、NIR、NDVI、NDWI、NDBI與珞珈一號數(shù)據(jù)和Sentinel-2數(shù)據(jù)的RED、NIR、NDVI、NDWI、NDBI與珞珈一號(Luojia-1)數(shù)據(jù)各自生成10 m的安徽省土地利用分類圖結(jié)果,如圖2所示。

(a)Landsat & Sentinel-2 & Luojia-1 (b)Sentinel-2 & Luojia-1 (c)Landsat & Luojia-1

基于SVM分類結(jié)果,不透水面類型提取結(jié)果如圖3所示。

(a)Landsat & Sentinel-2 & Luojia-1 (b)Sentinel-2 & Luojia-1 (c)Landsat & Luojia-1

2.2 分類結(jié)果驗證

通過在谷歌地圖上選取樣本點,并用混淆矩陣對分類結(jié)果進行精度驗證,谷歌選取的精度驗證點如圖4所示,驗證結(jié)果如表3所示。

圖4 精度驗證點審圖號：皖S(2019)13號

表3 多源數(shù)據(jù)源精度驗證結(jié)果與不透水面面積統(tǒng)計

通過對圖2、圖3、表3分析可知,針對2018年安徽省不透水面的提取,用Landsat、Sentinel-2和珞珈一號三種不同來源不同分辨率的數(shù)據(jù),精度達到0.97,說明精度較高。不透水面的面積統(tǒng)計是19194.283434 km2,面積也較準(zhǔn)確。而Landsat 8 &珞珈一號和Sentinel-2 &珞珈一號分類精度比Landsat & Sentinel-2 &珞珈一號三種結(jié)合分類精度要低,且面積統(tǒng)計相差較大。

通過圖5可以看出,使用三種不同數(shù)據(jù)源分類,在提高精度的基礎(chǔ)上,如圖5(a),Landsat & Sentinel-2 &珞珈一號結(jié)合的分類能夠?qū)⒆溆谔锏鼗蛏酱謇锏拇迓浜偷缆窚?zhǔn)確提取。如圖5(b)中,Sentinel-2 &珞珈一號則不能夠準(zhǔn)確提取出此類型的地物。圖5(c)中Landsat 8 &珞珈一號則會出現(xiàn)錯分,易將村落分為林地類。

(a)Landsat & Sentinel-2 & Luojia-1 (b)Sentinel-2 & Luojia-1

3 結(jié)論

通過利用多源高分辨率遙感影像對安徽省不透水面的提取,并對提取的結(jié)果進行了精度驗證與分析,得到了以下結(jié)論：

(1)通過對安徽省不透水面的提取,可以看出，Sentinel-2數(shù)據(jù)、Landsat數(shù)據(jù)和珞珈一號數(shù)據(jù)所組合的有11個波段對不透水面面積提取的精度最高。

(2)通過實驗發(fā)現(xiàn),使用Landsat、Sentinel-2和珞珈一號三種數(shù)據(jù)源結(jié)合提取不透水面,對小村落提取效果較好。相比之下,單純只用Landsat或者Sentinel-2數(shù)據(jù),小村落很容易識別為林地或農(nóng)田,誤差影響較大。

(3)通過選取特征點進行驗證,研究中生成的安徽省不透水面專題圖達到的精度較高,分辨率也較高。安徽省全省不透水面面積接近于19 194.283 434 km2,因此,研究方法切實可行。

(4)本研究所采用的影像分辨率不高,在進行重采樣后造成單個像元存在多種類型的地物,對分類的結(jié)果產(chǎn)生了一定影響。因此,在后續(xù)研究中,可嘗試使用高分辨率的遙感影像。

4 結(jié)束語

隨著三調(diào)的推進,遙感影像解譯的準(zhǔn)確性和效率性日益重要,基于SVM方法的多源高分辨率遙感影像分類及提取能夠滿足準(zhǔn)確性和效率性的社會特征,能夠很好地協(xié)調(diào)節(jié)省人力物力財力。