基于Terra%2fMODIS时间序列数据洞庭湖地区水域面积动态监测研究.pdfVIP

基于Terra/MODIS 时间序列数据的洞庭湖地区水域面积动态监测研究 Dynamic monitoring of the water area variations in Dongting Lake region using Terra/MODIS time series data 李景刚，李纪人，黄诗峰 LI Jing-gang, LI Ji-ren, HUANG Shi-feng 中国水利水电科学研究院 遥感技术应用中心，北京 100044 博士生 Email: sharp818@163.com 1、绪言 为了有效探讨250m MODIS 数据在湖泊水域面积变化动态监测中的可行性，进而为深入 探讨三峡工程建成运营后江湖关系奠定基础。本文以覆盖东、西、南3 个主要洞庭湖区，北 纬28°39′20″～29°32′30″和东经111°55′10″～113°39′45″间，总面积约 17300km2 的一个矩形区域为研究对象；利用2005 年总共47 景Terra/MODIS 卫星8 天合成的 250m 地表反射率数据（MOD09Q1）构建的时间序列数据集，采用多源信息提取的方法对洞庭 湖地区2005 年水域面积变化进行了动态监测。另外，通过全年最大淹没频率指数的构建， 完成了对洞庭湖地区重点水域淹没风险评价。 2 、研究方法 2.1、基于多源信息的研究区水面提取 通过反复试验，我们确定了如下两条规则： 1）如果由MOD09Q1 数据计算得到的NDVI 指数满足NDVI = 0.1 ，则设定其为水体， 以区分水体和植被； 2 ）如果MOD09Q1 第二波段的DN 值满足DNband2 1000 ，即反射率小于0.1 时，则也 设定其为水体，从而来区分水体与裸地和城镇。 另外，分别对MOD09Q1 第一和第二波段的DN 值设定阈值来消除云的影响；并利用由 SRTM 数字高程数据生成的研究区地形坡度图来消除山体阴影的影响。最后，对研究区2005 年最终采用的35 景数据的水面提取结果，分别进行分组求并运算，得到各月份的最大水域 面积分布图。 2.2 基于多源信息的研究区水面提取结果精度检验 首先，对完成几何纠正的 ASAR 后向散射系数数据，通过阈值法来提取研究区内的水 体面积。然后，同MOD09Q1 数据提取结果进行叠加分析，进而对MOD09Q1 数据提取的 水域面积精度进行检验。 2.3 研究区全年最大淹没频率指数构建 对于全年淹没频率指数，我们定义其为某一区域内被水体淹没天数占全年总天数的比 重。受研究区天气状况以及MODIS 数据本身质量的影响，要获取研究区每天的水体淹没范 围非常困难。因此，用前面MODIS 数据提取的研究区各月最大水体分布范围来代替该月的 水体淹没范围，计算得到的全年淹没频率指数称之为全年最大淹没频率指数。 3 、结果分析 3.1 研究区2005 年水面提取结果精度分析 Terra/MODIS 数据水面提取的制图精度为79.80%，漏提误差为21.20% ；而水面提取的 用户精度则达到了89.33%，错提误差为11.67% （表1）。考虑到MODIS 数据的空间分辨率 以及数据质量等因素的影响，本文认为该提取精度是可以接受的。 表 1 研究区水面提取精度分析矩阵 MODIS 提取结果 非水面/% 水面/% 总和/% ASAR 提取结果 非水面/% 88.26 1.02 89.28 水面/% 2.17 8.55 10.72 总和/% 90.43 9.57 100.00 水面制图精度：8.55 / 10.72 = 79.80%；漏提误差：21.20% 水面用户精度：8.55 / 9.57 = 89.33%；错提误差：11.67% 注：