基于gf-1与landsat-8影像的土地覆盖分类比较-宋军伟.pdf
《基于gf-1与landsat-8影像的土地覆盖分类比较-宋军伟.pdf》由会员分享,可在线阅读,更多相关《基于gf-1与landsat-8影像的土地覆盖分类比较-宋军伟.pdf(9页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第35卷第2期2016年2月地理科学进展Progress in Geography、,0135No2Feb2016基于GF1与Landsat8影像的土地覆盖分类比较宋军伟,张友静+,李鑫川,杨文治(河海大学地球科学与工程学院,南京210098)摘要:高分一号(GF1)卫星具有多种分辨率与大幅宽结合、重访周期短等优势,而Landsat-8卫星具有多波段、高辐射分辨率等优势。针对不同传感器参数特点,利用支持向量机分类器(sVM)对同区域同期两种数据进行土地覆盖分类对比研究。结果表明:两种传感器对应波段决定系数均大于092;典型样本的光谱趋势一致性良好,但在农田与林地、不透水面与裸土的典型样本可分离
2、性方面,Landsat-8优于GF1;GF1与Landsat一8的分类总精度分别为9038和9007,但不同地物类型的分类精度存在差异,波谱响应函数的差异可能是导致Landsat_8对林地的分类精度高于GF1的原因;此外,GF一1对零碎分布地物类型的分类精度高于Landsat8,主要原因是GFl具有更高的空间分辨率。关键词:土地覆盖分类比较;地物波谱;地物可分性;GF1影像;Landsat一8影像1引言土地覆盖主要描述地球表面的自然属性。土地覆盖变化是全球变化研究的热点问题之一,也是区域水文、气候和旱情监测的基础数据(赫英明等,2013;左玉珊等,2014)。遥感数据因其具有大范围、时效性强和
3、成本低的优势成为土地覆盖数据获取的主要数据源(刘爱霞等,2006;杨莹等,2012;陈趁新等,2014)。20世纪以来,许多搭载新型卫星的传感器相继发射,针对不同问题,不同传感器有不同的参数特点。在中低尺度研究中,Landsat系列卫星应用最为广泛,但Landsat数据幅宽相对窄、重访周期长,在较大区域或流域研究中需要不同时相数据的拼接,因而难以保证数据的一致性;此外相对较低的空间分辨率对地物破碎的区域,难以获得较高的分类精度。高空间分辨率影像具有地物内部分异明显,纹理、细节丰富等特点(黄昕等,2007)。现有的高空间分辨率卫星如SPOT、IKONOS、Ouicl(Bird等在区域研究中应用较
4、多(张锦水等,2006),但波段少、幅宽小,在研究中仍需要不同时相数据的拼接,且价格高昂。高分一号(gf_1)卫星具有重访周期短、多种分辨率与大幅宽结合的特点,但也存在波段数少的问题。目前GF1数据已应用于各类作物面积提取(王利民等,2015)、水体提取(李艳华等,2015)以及植被参数提取(刘吉凯等,2014)等方向;同时还与2013年2月发射的同期Landsat8数据进行对比分析,如玉米LAI反演(贾玉秋等,2015)、海冰探测能力(张晰等,2015)等,而GF1与Landsat8影像的土地覆盖分类能力比较研究还未见报道。因此,探讨和比较GF一1数据对区域主要土地覆盖的分类能力和分类精度,
5、对于有效和全面利用GF1数据具有实际的应用价值,也可为国产遥感卫星传感器设计提供参考。为了对比不同传感器参数对土地覆盖分类能力与精度的影响,评价GF一1在土地覆盖分类中的应用价值,选取湖北钟祥市同日Landsat一8和GF1影像数据,分析对应波段相关性,采用SvM分类器收稿日期:201508;修订日期:2叭5一ll。基金项目:国家科技重大专项(08一Y3080790011315)Foundation:National Science and Technology Major Project ofchina,No08-Y30807-900l1315】。作者简介:宋军伟(1991一),男,陕四咸阳
6、人,硕士研究生,主要从事遥感技术机理与应用研究,Email:695276282qqcom。通讯作者:张友静(1955一),男,江苏南京人,教授,主要从事遥感、GIs教学与科研工作,Email:zhangyjhhueducn。瓢用万方数据256 地理科学进展 第35卷进行土地覆盖分类,对比GF1与Landsat一8数据的分类结果差异并分析造成差异的主要原因。2研究区域、数据与方法21研究区概况研究区位于湖北省钟祥市(30042N3 l 036N,112。07E113。00E)(图1)。该市位于汉江中游、汉江平原北端,地处亚热带季风性气候的过渡带上,光、热、水同步,具有多宜性的生态环境。区域内主要
7、土地覆盖类型包括林地、裸土、水体、农田和不透水面等。研究区西部海拔较高地区土地覆盖类型为林地;东部丘陵地带以林地和坡耕地为主;汉江两岸的平原地区以农田为主,农田与分散的居民点呈混合分布状态,对比较不同分辨率影像中零碎分布地物的分类效果较为典型。城区位于研究区中部汉江东侧的平原地区。22数据与方法221高分一号(GF1)数据GF一1卫星搭载有两台2 m分辨率全色和8 m分辨率多光谱高分相机(PMS),四台16 m分辨率多光谱宽幅相机(MV)(白照广,2013;翁进等,2015)。本文选取了一景2013年8月6日获取的高分一号wFV3影像数据,轨道号:210(Row),4(Path),太阳高度角6
8、8440,太阳方位角12591。,波谱范围为045089m,分为4个波段,整个区域内无云,可视性良好。具体参数见表1。222 Landsat一8数据Landsat一8携带有0LI(Operational Land Image)和TIRs(The册alInfrared sensor)j隹扫式成像仪,其中oLI包括9个波段,除全色波段为15 m分辨率外,其余波段的空间分辨率为30 m(徐涵秋等,2013)。TIRS包括2个热红外波段,空间分辨率为100 m。选取同一时间、同一区域的Landsat8卫星oLI传感器数据,与GF一116m分辨率多光谱数据进行土地覆盖分类对比分析。轨道号:38(Row)
9、,124(Path),太阳高度角6435。,太阳方位角12096。,波长范围为043138 um,分为9个波段,本文采用L1级产品。具体参数见表1。同时选择同期的25 m分辨率的sPOT5多光谱图l研究区位置Fig1 Location of the study area表1 Landsat8与GF1卫星参数对比Iab1 Comparison of satellite parameters between Landsat-8 and GF一1万方数据第2期 宋军伟等:基于GF1与Landsat8影像的土地覆盖分类比较 257与全色波段融合影像(数据来自Go091e eanh),主要用于分类训练样
10、本的选择和分类结果的验证。223数据预处理对两景影像进行辐射定标,得出大气上行辐亮度数据,通过FLAASH进行大气校正,将其转换为地面真实反射率数据。利用矢量数据对影像进行裁剪得出研究区。投影选择UTMzone49N投影,wGS84地理坐标系。以GF一1数据为基准,选取控制点采用二次多项式法对Landsat一8数据进行相对几何校正,刚ES为05456。3分类方案及方法31分类方案分类体系建立参考国土资源部土地利用现状分类一级分类系统,并考虑所用数据及研究区土地覆盖类型的特点,最终确定将研究区的土地覆盖类型划分为林地、水体、农田、不透水面和裸土五类。32分类方法目前利用遥感影像进行分类的主要方法
11、有监督分类和非监督分类、面向对象的影像信息提取以及基于专家知识的决策树分类等,选择支持向量机分类器(Support Vector Machine,SVM)对研究区进行土地覆盖分类(骆建承等,2002;张锦水等,2006)。支持向量机(sVM)是20世纪90年代中期发展起来的基于统计学习理论的一种机器学习方法,通过寻求结构化风险最小来提高学习机泛化能力,实现经验风险和置信范围的最小化,从而达到在统计样本量较少的情况下,亦能获得良好统计规律的目的(骆建承等,2002;张友静等,2006)。其独有的特性和高效性使得该分类器在遥感分类中得到广泛的应用。svM利用核函数将低维空间向量机映射到高维空间,从
12、而在高维空间寻找线性分隔超平面使线性分隔最大。选用径向基函数(Radial BasisFullction,简称RBF)作为本次分类的核函数。最常用的径向基函数是高斯核函数,公式为: 1 -口、K(x,xc)=expl一坐掣I(1) Z盯 式中:Kfx,xC)为空间中任一点x到某一中心xC之间欧氏距离的单调函数,xC为核函数中心,盯为函数的宽度参数,控制了函数的径向作用范围。4结果与分析41两种传感器波段对比分析为了对比两种传感器在红、绿、蓝和近红外波段反射率差异,选取各类典型地物(林地、农田、水体、不透水面、裸土)共7293个像元进行相关性分析(Hao et a1,2014;zhong et
13、al,2014),分析结果见图2。j噻0 O05 O1 O15 O2 025 03 O 35 0 O05 O1 0 15 02 025 O_3 O_35 04Landsat一8蓝波段反射率 Landsat_8绿波段反射率O 01 O2 O3 O4 O 5 0 02 O 4 O6 O 8Landsac一8红波段眨射率 Landsat-8近红外波段反射率图2 Landsat8与GF1影像对应波段相关性分析F培2 Co仃elation analysis of corresponding charulels be铆een Landsat-8 and GFl邪”如笛如m晒0万方数据258 地理科学进展
14、第35卷通过相关性分析得到Landsat8和GF1各对应波段的线性回归方程的决定系数分别为094、092、093和097,均达到极显著性水平(P0001),说明两种传感器的反射率存在显著的相关关系。蓝、绿和近红外波段的斜率略小于1,红光波段斜率略大于1,原因是GF1红光传感器光敏范围较大,进人了o7075 um范围,近红外波段传感器也进入了该范围,而Landsat8没有进人这一区域(图3),该范围处于植被“红边”区域。当样本中林地较多时,Landsat8在红光波段的反射率值会偏低,同时“红边”区域只是植被从吸收转向反射的快速过渡区域,植被在这一区域的反射并未达到最强程度,因此植被对进入这一区域
15、的GF1近红外波段反射会有所降低(徐涵秋等,2011;贾玉秋等,2015)。这一差异将会导致GF1对植被的接受信号减弱,在土地覆盖分类时易造成林地与农田的混分。42典型样本分析为了获取高质量且具有代表性的样本数据,利用25 m分辨率的SPOT5多光谱与全色波段融合影像选择分类样本数据(黄婷等,20 1 3),同时保证样本数据在整个研究区均衡分布,采用分层随机抽样方法进行钡0试样本选择(刘旭拢等,2006),训练样本个数及特征见表2。分析在不同传感器中各类地物的反射率差异,计算各类地物在不同波段的反射率平均值,得到典型地物类型样本反射率光谱曲线(图4)。从图4可以看出,两种传感器中对应地物类型的
16、反射率曲线趋势一致性良好,在近红外波段的反射率差异最小。在蓝光波段、绿光波段和红光波段3)() 45f】 )Ll 65(J j(】 H)【) u)() l(1)【l l 5Il波长向,l图3 GF一1与Landsat-8在可见光及近红外波段波谱响应函数Fig3 Spectral response functions(SRFs)of the Visible andNIR channe】s of GF1 and Landsat8Landsat一8各地物的反射率略大于GF1,其中不透水面最为明显,裸土和农田次之,水体的反射率最为接近。为了对比两种数据各类地物的光谱可分离性,计算不同类别样本问的统计距
17、离。类别间的距离基于Je倚iesMatusita距离和转换分离度(Transfornled Dive唱ence)(刘新杰等,2013;许伟等,2015),衡量各类别的可分离性,确定各类别间的差异程度。表3列出了高分一号数据和Ladsat一8数据6组训练样本可分离性的计算结果。从表3看出,农田和林地,不透水面和裸土两组的类别可分离性Landsat8高于GF1,其他各组的差异不大,这说明在光谱可分离性方面,Landsat一8优于GF一1。43分类结果分析选取相同的训练样本(表2),利用支持向量机对表2典型样本选择11ab2 SeIection of typical samples林地水体不适水c|
18、i j?j蓼呈深红(!:!,纹理特征粗l糙大蓝色到深牖包之|1】,纹理均譬红色,确HJJ湿的纹 2993州特征亮度高,里蓝色夹淡 2029绿色图4典型地物反射率光谱曲线Fig4 Spectral curves ofrenectance oftypical surface types一、墨、3。、圈OOOO0O0OO主喜装裂万方数据第2期 宋军伟等:基于GF1与Landsat_8影像的土地覆盖分类比较 259表3 GF1与Landsat8数据可分离性比较Tab3 Comparison of separabil姆be铆een GF一1and Landsat一8 data两景影像进行分类,其中GF一
19、1利用WFv传感器的4个波段、Landsat8利用OLI传感器除去卷云波段和热红外波段的其余8个波段进行分类。对分类结果进行聚类和筛选处理,去除孤立点并对细小斑块进行合并,分类结果见图56。从图56看出,研究区中农田的面积最大,林地主要分布在中西部和东部海拔较高地区。两景影像中水体的分类结果较为一致,林地与不透水面的分类差异较大,分别分布在研究区中西侧和南部汉江两倾l|o图5 GF1土地覆盖分类结果Fig5 LandcoVer classification result of GF一144分类精度评价采用分层随机抽样点作为精度评价方法,利用25 m分辨率sPOT5多光谱与全色波段融合影像目视判
20、读,对分类结果进行精度评价,测试样本个数见表2,表45为分类结果的混淆矩阵。从表4。5混淆矩阵可看出,GF1和Landsat一8的分类总精度较高,分别为9038、9007,分类结果理想。Landsat一8对林地的分类精度高于GF1,主要因为波谱响应函数的差异造成GF一1将更多的林地误分为农田。GF一1对水体与不透水面的分类精度高于Landsat8,因为研究区分布较多水塘、细小河流、村庄和小规模的城镇,地物类型破碎,在更高分辨率GF1影像上混合像元少,纹理与细节信息明显,从而分类效果更好。GF一1将部分收割农田和未利用耕地分为裸土,是造成裸土分类精度低于Landsat8的主要原因。45分类差异比
21、较对两景影像的分类结果图进行叠加,分析分类结果差异,结果图如图7,差异面积比例见表6。图6 Landsat8土地覆盖分类结果Fig6 Land-cover classification result of Landsat一8表4 GF1分类混淆矩阵1ab4 Confusion matrix of GF1 classification分类类别实际类别林地 水体 农田 不透水面 裸土 总体 用户精度林地水体农田不透水面裸土总体生产精度6127261266143686988993984682O总精度:9038 Kappa系数:O875o:,。如矾泓oo。挖:巧刚4l5l78船M蛙。M。”刚30O,5
22、556趴万方数据260 地理科学进展 第35卷总精度:9007 K印pa系数:0853从图7、表6可看出,Landsat一8与GF1的分类结果有1619不一致。农田与林地的混分比例最高达877,其中有790的面积Landsat8将其分为林地而GF一1分为农田,差异较大的原因应为两种传感器波谱响应函数在红波段和近红外波段响应范围差异造成。有301的面积Landsat一8分类结果为不透水面而GF1分类结果为农田和裸土,264的面积Landsat一8分类结果为裸土GF1分类结果为农田,结合表4、表5混淆矩阵可看出,GF一1对不透水面的分类精度高于Landsat8,而对裸土的分类精度低于Landsat
23、8。造成差异的原因一方面由于不透水面与裸土光谱特征相似,易产生混分,部分收割农田和未利用耕地光谱特征与裸土相似,也易混分其中;另一方面高分辨率影像对零碎分布地物类型具有更好的识别度。选择分类差异较大区域,利用25 m分辨率的sPOT5多光谱与全色波段融合影像进行分类结果实地验证。验证图见图89。整体来看,在图8验证区Landsat8的林地面积明显大于GFl。选取差异明显区域A、B、C,与25m分辨率sPOT5影像对比验证。在A区域海拔较盏鋈翟馘ZZO 10 kmI_JZ工林地一农田鬲不透水面一裸土一农|j一不透水面蚕 农【H一裸土。图7 GF1与Landsat8分类差异比较F培7 Di廊ren
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 gf landsat 影像 土地 覆盖 分类 比较 宋军伟
限制150内