基于时谱秦岭中段农作物提取研究.docVIP

基于时谱秦岭中段农作物提取研究 　　摘要以秦岭中段为研究区域,采用MODIS影像数据,利用时谱分析方法进行山区农作物提取研究。结果表明,利用时谱分析方法可以快速提取山区农作物,精度为87.9?。 　　关键词时谱;秦岭中段;农作物提取 　　中图分类号 S127 文献标识码A文章编号 1007-5739(2010)05-0014-01 　　 　　所谓时谱,主要指利用单时期反射率影像所获取的目标地物瞬时状态的光谱特征,以及利用多时次的遥感影像(尤其是植被指数数据)所获取的目标地物的季节性变化特征。植被是土地覆被的构成主体,而正是由于植被的主导地位,土地覆被的变化主要表现为植被之间以及植被与其他覆被类型间的变化,变化的结果可能造成其形态属性永久性改变,暂时性改变周期性变化。对于绝大多数土地覆被类型来说,几乎都具有特定的周期性变化规律,特别是在气候等驱动力作用下表现出明显的季节性变化。因此,季节性变化是土地覆被最本质的特性之一,也是土地覆被的光谱特征随时间变化的一种重要表现形式。随着国家财政对粮食生产补贴力度的不断加大,为政府部门提供准确的农作物种植信息十分重要。 　　1研究区概况 　　秦岭中段是秦岭山脉位于陕西境内的部分,位于北纬32°50′～34°30′,东经105°30′～110°0′;年均温在6～13℃;降水的季节变化明显,主要集中于夏季;平均海拔在2 000m以上,70%以上的山地海拔在500～2 500m[1];研究区内主要有两大河流,汉江源于宁强县,流经汉中市和安康地区;渭河由宝鸡市进入陕西省境内,横贯八百里秦川[2]。根据查阅资料和实地野外考察,秦岭中段农作物基本上分布在这两大流域。 　　2主要数据源及数据预处理 　　2.1数据预处理 　　采用的遥感数据主要包括空间分辨率为250m、时间分辨率为16d的MODIS反射率影像以及MODIS/NDVI和MO DIS/EVI植被指数数据。由于所下载的MODIS反射率数据和植被指数数据都是横跨h26和h27这2个条带的数据集,因此对2008年23个时相的反射率数据和46个时相植被指数数据,先利用ENVI软件对原始数据进行裁剪、拼接、重投影以及格式转换等方面的处理,最后得到初始时序的MODIS/反射率影像、MODIS/NDVI和MODIS/EVI影像数据。 　　连续3次分别于2007年8月、2008年8月、2009年5月前往研究区进行野外实地考察。考察的主要内容是利用GPS工具,对秦岭中段南北坡典型土地覆被类型进行准确定位,并详细记录样本点的属性;3次总共确定了1 000多个样点。 　　2.2时序MODIS植被指数重构 　　对比分析生成的初始时序的MODIS影像数据发现,有些时相的遥感数据上云污染较重,因此对初始时序的MODIS遥感数据作进一步的去云处理。采用White等人提出的利用BISE算法来消除云影响,具体的方法是:对于任一像元,如果其在某一时刻(ti)的NDVI(或EVI)值同时满足:(Pixel■- Pixel■)/Pixel■≥20%,并且 (Pixel■- Pixel■)/Pixel■≥20%(式中Pixel■、Pixel■、Pixel■分别为ti时刻、ti-1时刻和ti+1时刻的NDVI(或EVI)值),则表明该像元受到了云的影响。为此,可将ti-1和ti+1时刻的平均值作为该地物在ti时刻的NDVI(或EVI)值,即: 　　NDVI■′=(NDVI■- NDVI■)/2,或 　　EVI■′=(EVI■- EVI■)/2 　　据此,可消除原始的时序MODIS/NAVI和MODIS/EVI植被指数数据上残留云的影响。 　　3秦岭中段农作物提取结果及精度评价 　　3.1农作物及主要土地覆被类型时谱分析 　　土地覆被类型丰富,主要包括常绿覆被、季节性绿色覆被(例如落叶林和作物)、灰色覆被和蓝色覆被等[3]。在研究区冬季MODIS影像上,落叶林全株枯死或地上部分枯死,此时的色调最浅;而作物和常绿覆被呈现绿色调等。据此分别采集反射率及MODIS植被指数时谱数据。 　　从图1、图2可以看出,主要土地覆被类型MODIS反射率数据在近红外波段达到峰值;但是常绿覆被、作物由近红外波段到中红外波段呈现出骤然下降的趋势,而其他土地覆被类型下降缓慢或趋平。常绿覆被和作物的时谱曲线差异显著,常绿覆被呈现单峰型波动曲线,并且峰值出现时间位于每年的5月初至8月底,其中MODIS/NDVI最大值约为0.9,出现于每年的7月底至8月初,而MODIS/EVI最大值约为0.7;谷值在每年的11月下旬至翌年2月底出现,并且NDVI最小值约为0.45,EVI最小值则约为0.1;作物在1年周期内,即n1(MODIS/NDVI第1期)至n23(MODIS/NDVI第23