modis数据的处理方法.docVIP

MODIS数据的处理方法(ENVI) 美国RSI公司（Research Systems Inc.）的产品ENVI能很好地支持HDF数据格式。ENVI（The Environment for Visualizing Images）遥感影像处理软件，是分析、处理并显示多光谱数据、高光谱数据和雷达数据的高级工具。ENVI能接受大量的传感器数据，是世界目前唯一 17 5.MODIS数据的处理 5.1植被指数的选择 植被指数定量地表明了植被活力，它们比用单波段来探测绿色植被更具有灵敏 性，已经被广泛应用于遥感影像的解译中，如土地利用变化监测、植被密度覆盖评价、 森林、农作物测定和农作物预报等方面。 归一化植被指数NDVI是普遍应用的40几种植被指数中应用最为广泛的，归一 化植被指数NDVI（Normalization Difference Vegetation Index）被定义为近红外波段与 可见光红波段数值之差和这两个波段数值之和的比值。即 NIRRED NIRRED ρρ ρρ + ?N DVI =5－1 其中 RED ρ和 NIR ρ分别是红光（第1波段）和近红外（第2波段）波段的反射率。 NDVI可以用于季节性植被变化的研究，也可用于年际间植被变化的分析、干旱 地区植被变化规律等等，本文研究区——内蒙古自治区为干旱半干旱地区，且要分析 内蒙古地区植被变化中的月际、年际变化，因此，本文选用NDVI作为多年植被变化 分析的植被指数。 5.2MODIS数据的预处理 论文所用的MODIS数据来源于内蒙古自治区遥感与地理信息系统重点实验室 MODIS接收处理系统存档的2002－2006年5－8月每日的TERRA卫星接收的250m 分辨率的MODIS 1B数据。 MODIS 1B数据预处理主要包括反射率数据的生成、几何纠正这两步。这些工作 都由遥感图像处理软件ENVI 4.2来实现。 5.2.1几何纠正及Bowtie处理 由于MODIS采用“多元归并”的对地遥感方式，即多个探测器同时扫描，导致 通道探测器之间的光谱响应不可能做到完全相同，尽管在发射前做了相对定标，但随 着时间推移或仪器本身的问题，探测器之间还是会产生差异，这就会在图像上会产生 “条纹”，这就是所谓的Bowtie效应。这个问题会不同程度的对MODIS数据的使用 造成影响，在越接近星下点的区域，Bowtie效应越小，在越接近轨道边缘的区域， Bowtie效应影响越大。所以一定要对MODIS数据的Bowtie现象进行处理。18 由于MODIS 1B数据（HDF文件）中已经包含有经纬度数据，我们可以利用这 些经纬度数据进行几何纠正。利用ENVI 4.2软件，将HDF文件中经纬度数据导出为 地面控制点（GCP）文件，然后采用ENVI 4.2软件的MAP功能模块下的几何纠正功 能进行影像纠正。利用ENVI 4.2的这一功能，可以达到0.3个像元的精度。同时ENVI 4.2在几何纠正的同时对Bowtie也进行了处理。 在几何纠正的过程中采用的投影为Albers等面积圆锥投影（采用克拉索夫斯基 椭球体坐标系）具体参数如下： 中央子午线：105°00′00″E 南基准纬线：25°00′00″N 北基准纬线：47°00′00″N 中央纬线：00°00′00″N 5.2.2 1、2波段反射率数据的计算 MODIS 1B数据对于反射太阳波段提供反射率数据产品，为了减小浮点型数据占 用的存储空间，MODIS 1B采用16-bit整数和尺度转换的方法，通过偏移量（offset） 和尺度因子（scale）两个参数转换为16-bit的整数型数据。为了得到具有物理意义的 数据产品，需要再次通过尺度转换的方法把整形数据转换成具有物理意义的浮点数 据。 太阳反射波段反射率计算公式为： Rreflactance_scales(DNreflatance_offsets) iiii =×?5－2 其中 i R是i波段对应的反射率值， i DN为MODIS 1B数据中i波段的整数型 数据， i reflactance_scales为波段反射系数、为 i reflatance_offsets波段偏移量。 5.3植被指数的计算 5.3.1 NDVI的计算 用MODIS数据1、2波段的反射率数据反演出归一化植被指数NDVI。 NIRd NIRd NDVI Re Re ρρ ρρ + ?= 5－3 其中 Red ρ为红光波段反射率（对应于MODIS数据的第1波段）， NIR ρ为近红外 波段反射率（对应于MODIS数据的第2波段）。5.3.2 MVC法合成月NDVI图 用最大值合成法MVC合成月NDVI图（NDVI-MON）。 max() i,ji,k NDVI=NDVI5－4 其中 NIR ρ为近红外波段反射率（对应于MO