利用Modis―NDVI时序数据分析滇中地区绿色覆被变化特征.docVIP

利用Modis―NDVI时序数据分析滇中地区绿色覆被变化特征 　　摘要： 本文通过Modis 1B数据计算了2009～2013年滇中地区内NDVI值，并用最大值合成法（MVC）获得了NDVI月时间序列数据，利用相关性分析和差值分析定量描述了农田、森林、灌丛和草地4种绿色覆被NDVI年内变化曲线在年际时间尺度上的协同性和差异性。 　　Abstract： This paper calculates the NDVI values within the 2009-2013 in the central Yunnan Province by Modis 1 B data， and NDVI monthly time series data is obtained with maximum synthesis（MVC）. By using correlation and difference analysis， the synergies and differences of NDVI curve during the year on interannual time scales of the 4 types of green cover： farmland， forest， scrub and grassland are described quantitatively. 　　关键词： MODIS；NDVI；相关系数；协同性；滇中地区 　　Key words： MODIS；NDVI；correlation coefficient；collaborative；central Yunnan Province 　　中图分类号：Q948 文献标识码：A 文章编号：1006-4311（2014）22-0320-02 　　1 研究区概况 　　滇中地区位于我国“两纵三横”城市化战略格局中包昆通道纵轴的南端，目前是我国主体功能区规划的重点开发地区。滇中地区位于云南省中部，介于100°00′27.11″～103°4′20.41″E，22°9′55.19″～27°00′24.71″N之间，地处红河、南盘江、金沙江、澜沧江四大水系的分水岭地区，地形属于高原，地形小波状起伏，平均海拔为2000m，多为起伏和缓的低山和丘陵，发育着类型多样的岩溶地形。其西北地区处于横断山脉的边缘，相对高差大，地势险峻。而东南地区的地势较为缓和，有着非常开阔的河谷，海拔一般在800m～1000m之间。整个区域的总面积为17.4万km2，占整个云南省面积的45.29%，主要包括昆明、曲靖、玉溪、楚雄、丽江、大理、红河七个市、地区的49个县。其主要的绿色覆被有林地、草地、灌丛和农田四大类别。 　　2 数据资料与研究方法 　　2.1 数据资料及其预处理 本文选用了美国的对地观测计划ESO的ERRA卫星携带的MODIS传感器的轨道数据，综合考虑了影像获取的难度、影响的质量、所面临的工作量以及研究目标等因素。在2009～2013年每天接收的资料中，选取分辨率为250m的1B级数据作为其基础性研究资料（表1），共计为1261个时期的影像。在处理数据时，研究人员选用ENVI软件进行影像批处理。同时，以第四代交互式计算机编程语言IDL为基础，开发出了影像批预处理程序，主要在ENVI软件平台上运行。该程序采用了Delaunay三角网坐标转换法以及双线性插值法，对获取到的影像进行了几何校正，并在校正过程中采用后向映射算法有效地消除了“双眼皮”现象，最终校正完成后的坐标系统为UTM投影，WGS84坐标系统。研究人员在此基础上采用了辐射定标和去云处理的方式裁剪出了滇中研究区的影像。 　　此外，为了研究不同植被类型NDVI时间序列数据季相及年际间的变化特征问题，在Modis MOD12（Modis土地利用/覆被数据，分辨率为500m）数据、云南省1：25万地形图、云南省1：50万植被盖度图等资料的辅助下，对2005年Modis MOD02（Modis经过辐射定标后的数据，分辨率为250m）数据进行了决策树分类。参照国际上IGBP的土地覆被分类系统，分出林地、灌丛、草地、农田、居民地、水体等六类地物类别。 　　2.2 NDVI年际变化特征的协同性分析 植被NDVI值在时间尺度上的变化反映了植被年内生长的NDVI变化过程，不同植被类型的NDVI时间曲线具有不同的变化特征，通过对不同植被类型的NDVI时间曲线的分析可以了解该植被的年内生长状况。本文采用NDVI特征值的分析方法计算滇中地区NDVI值的年际变化特征。选取相关系数r分析NDVI 曲线年际间的相关性，指示植被变化年际间的协同性。相关系数r的计算公式： 　　rxy=■（1） 　　式（1）中，rxy为第x年和第y年每一象元ND