基于MODIS数据洞庭湖水面面积估算方法.pdfVIP

第4()卷第 l4期 人 民 长 江 Vo1．40，No．14 2009年 7月 Yangtze River July， 2009 文章编号：1001—4t79(2009)14—0040—04 基于MODIS数据的洞庭湖水面面积估算方法 龚 伟 杨 大 文 钱 群 (清华大学 水利水电工程系，北京 100084) 摘要：利用美国Terra卫星上MODIS传感器提供的NDVI数据产品，采用两种方法分别在年和月时间尺度上对 2000～2007年间洞庭湖的水面面积进行了估算。在年尺度上，采用 ISODATA方法对NDVI时间序列进行非监 督分类，区分水面及各种湿地类型，并分析了各类湿地的面积变化趋势。在 月尺度上，采用 16d平均NDVI数 据，通过NDVI阈值提取水面，在不同的季节选取不同的闽值 ，并将基于NDV1数据得到的水面面积与基于水 位一面积曲线得到的水面面积进行了对比分析。在此基础上，探讨 了两类不同尺度上分类方法的优点、缺点及 适 用范围。 关 键 词：水面变化；MODIS；NDVI；非监督分类；洞庭湖 中图分类号：X171．1 文献标识码 ：A 本文采用了美国NASA网站上提供的TerraMODIS16d合 1 数据来源 成的植被指数产品MOD13Q1，该数据的空间分辨率为250m。 本文使用的卫星遥感数据为美国Terra卫星的中分辨率成 MOD13Q1产品选择为 16d中无云数据的合成，合成过程中优 像光谱仪(MODIS)传感器提供的归一化植被指数(NDVI)数据。 先选择近星下点无云的像元，并尽可能减小云、阴影、大气中气 Terra(EOS—AM1)是美国地球观测系统(EOS)的先进极地轨道 溶胶的影响，并应用二向反射率(BRDF)模型对观测角度进行 环境遥感卫星，于 1999年 12月 18日发射 ，可提供2000年以后 了订正 J。研究区域范围为北纬 28．5。一29．5。、东经 111．75。 的数据。该卫星通过洞庭湖的时间为当地时间10：30和22：30。 一 ll3．25。，时间范围是 2000年 2月至2o07年 l2月。 其中MODIS是EOS—AM1系列卫星的主要探测器之一 ，它是新 在数据处理过程中，首先采用 MODIS投影变换软件 (MRT 一 代 “图谱合一”的光学遥感仪器，具有36个光谱通道 ，分布在 工具)进行拼接和投影变换，参考坐标系选用WGS84椭球体和 0．4—14 m的电磁波谱范围内，地面分辨率分别为250、500111 UTM投影，重采样后的分辨率分别为250m。尽管 16d合成后 以及 l000nl。多波段数据可以同时提供反映陆地、云边界、云 NDVI数据质量相对于每天数据有了较大提高，但是受云和阴影 特征、海洋水色、浮游植物、生物地理、化学、大气中的水汽、地表 以及辐射校正精度的影响，下载的NDVI数据仍会出现色调不 温度、云顶温度、大气温度 、臭氧、云顶高度等特征信息。 一 致的情况，影像上常出现各种暗影、黑斑、白斑。为保证成果 植被指数是一种通过多波段数据组合来识别地表土壤 、植 质量 ，参考了MODIS陆地工作组(MODLANDST)提供的数据质 被 、水体的综合指标。植被在红光波段 (0．60—0．70 m)有一 量评估结果 】，如果研究区域内超过20％的像元质量不佳，则 个强烈的吸收带，它与叶绿素密度成反比，植物的反射光谱曲线 不采用这幅影像。2000—2007年共有 180幅影像 ，剔除质量不 因叶绿素对这部分光的强烈吸收而呈波谷形态 ；而在近红外波 佳的34幅之后，还剩 146幅质量较佳的影像用于本文的分析。 段(0．75一1．3 m)有一个较高的反射峰，它与植物叶绿素密度