动态时间弯曲技术支持下时序NDVI数据的耕地分布信息提取.PDFVIP

Resources Science 第36卷第9期 2014年9月 2014，36（9）：1977-1984 Vol.36，No.9 Sept.，2014 文章编号：1007-7588（2014）09-1977-08 动态时间弯曲技术支持下时序NDVI 数据的 耕地分布信息提取 ——以江苏省苏南地区为例 1 1，2 1，2 1，2 孙越凡 ，钟礼山 ，程 亮 ，李满春 （1.南京大学地理信息科学系，南京210023；2.江苏省地理信息技术重点实验室，南京210023 ） 摘 要：我国耕地保护形势十分严峻，及时获取耕地现状信息对耕地保护工作的开展具有重要参考意义。文 章探讨了动态时间弯曲技术在时序NDVI数据的耕地分布信息提取中的可行性，并以江苏省苏南地区为实验区，用 2010年MODISNDVI时间序列数据进行耕地提取实验。首先使用频域低通滤波对时序NDVI数据进行重构，抑制 大气、光照等噪声的影响，接着选取耕地训练样本，构建耕地参考时间序列，然后根据动态时间弯曲（DTW）方法计 算待分类像元时间序列与耕地参考时间序列的相似性程度，对计算结果使用阈值分割，得到研究区耕地信息。实 验中耕地提取结果正确率达84.08%，完整率达85.68%，证实了该方法的适用性。研究表明，基于时序NDVI数据相 似性分析可以快速有效提取大面积耕地信息，为农业政策的制定提供参考依据。 关键词：耕地分布；NDVI；时间序列；相似性分析；DTW 1 引言 且经常因为大气环境因素导致影像不可用。Terra、 耕地是地表覆盖的主要类型之一，是人类赖以 Aqua等中低分辨率EOS系列卫星扫描幅宽大，重访 生存的基本资源，耕地的变化对粮食安全和生态环 周期短，通过2颗卫星搭载的中分辨率成像光谱仪 境稳定性有着非常重要的影响。近年来高速的经 （MODerate-resolution Imaging Spectroradiometer， 济发展与城市化进程占用了大量耕地，加上生态退 MODIS）传感器协同工作，每天可以得到2次覆盖全 耕还林、农业结构调整等原因，我国耕地总量持续 球的观测数据，经过最大值合成处理等处理去除云 下降，人均耕地面积不足世界平均水平的40%[1]。 雾影响，能够形成稳定的遥感时间序列产品。 为保障国家粮食安全，国务院提出坚守18亿亩耕地 全球耕地类型多样，根据所生长作物的种类和 红线，实现耕地资源的可持续利用，而及时获取耕 生长阶段不同，耕地的影像光谱特征、纹理特征也 地现状信息对开展耕地保护工作有重要参考意义。 复杂多变。然而某一基准年度全球范围内多季相 卫星遥感可以从不同尺度、以不同分辨率和光 TM/ETM+影像通常难以获取，全球耕地提取只能基 谱进行对地观测，能够快速获取大范围耕地种植信 于单时相影像进行。仅依靠单时相TM/ETM+影像 息，在资源调查与变化监测等方面已经得到广泛应 提取耕地，部分区域光谱纹理特征不明显，容易造 用。不同种类耕地在不同季相和物候期有不同的 成耕地与其他类型的混分[3]。为此，在现有方法的 光谱响应，利用多时相遥感数据进行耕地提取相对 基础上，引入多时相的归一化植被指数（Normalized 单一时相而言更具优势[2]。中高分辨率卫星如 DifferenceVegetationIndex，NDVI）用于耕地提取， Landsat、SPOT、IKNOS等卫星扫描幅宽小，重访周期 提高耕地提