hants时间序列谐波分析法.docxVIP

HANTS ( the Harmonic analysis of time series) 时间序列 谐波分析法 时间序列谐波分析法（Harmonic Analysis of Time Series, HANTS）是平 滑和滤波两种方法的综合，它能够充分利用遥感图像存在时间性和空间性的特 点，将其空间上的分布规律和时间上的变化规律联系起来。时间序列谐波分解法 进行影像重构时充分考虑了植被生长周期性和数据本身的双重特点，能够用代表 不同生长周期的植被频率曲线重新构建时序NDVI影像，真实反映植被的周期性 变化规律。时间序列谐波分析法是对快速傅立叶变换的改进，它不仅可以去除云 污染点，而且对时序图像的要求不象快速傅立叶变换（FFT）那么严格，它可以 是不等时间间隔的影像。因此同快速傅立叶变换相比，HANTS在频率和时间系 列长度的选择上具有更大的灵活性。时间序列谐波分析法进行时序影像的重构也 是基于云对NDVI的负值影响，但是它与最大值去除云污染的影响是两个完全不 同的方法。它是首先通过傅立叶变换得到非零频率的振幅和相位，然后将所有的 点进行最小二次方拟合。通过观测资料与拟合曲线的比较，对于那些明显低于拟 合曲线的点被作为云污染点通过把它们的权重赋为零而拒绝参与曲线的拟合。建 立在剩余点上进行新的曲线拟合，通过这种反复进行的迭代过程实现图像的重 构。 HANTS的核心算法是最小二乘法和傅立叶变换，通过最小二乘法的迭代拟合去除时序NDVI值中受 云污染影响较大的点，借助于傅立叶在时间域和频率域的正反变换实现曲线的分解和重构，从而达到时序 遥感影像去云重构的目的。 采用时间序列谐波分析法（HANTS）可以对时间谱数据进行平滑。其核心算法是傅立 叶变换和最小二乘法拟合，即将时间谱数据分解成有限个谐波（正弦波或余弦波），从中选 取若干个能反映影像时序特征的谐波进行叠加，达到重构时序数据的目的。 （左丽 君MODIS/NDVI和MODIS/EVI在耕地信息提取中的应用及对比分析） 原始NDVI曲线与HANTS滤波后曲线Hants软件操作步骤1： projict设置输入输出、图像大小、波段等。6Q海Q2000 原始NDVI曲线与HANTS滤波后曲线 Hants软件操作步骤 1： projict 设置输入输出、图像大小、波段等。 6QQ0 4000 400C 海Q 2000 -2000 -2000 File Edit Options Plct^Function Help Spe-trd Profile File Edit OptiarLS Flct^Functiflii Help Spmgl Pro?— BaSd Nunnber 般a Nurnber 亨林2 Spectral Prof i le : Input dirOutput dinoutput dmprojectUnpackdata_list为时间序列数据的文件列表，需和数据文件放同一个目录下（数据为裸数据格式， envi标准格式即可），格式如下图所示（IMG只是加的一个标识，用于其他处理，并不是文 件的格式）。需要处理的 原始图像存 放的路径需要处理的 .京始图像样 本数m图像的宽和 孙」 Input dir Output din output dm project Unpack data_list为时间序列数据的文件列表，需和数据文件放同一个目录下（数据为裸数据格式， envi标准格式即可），格式如下图所示（IMG只是加的一个标识，用于其他处理，并不是文 件的格式）。 需要处理的 原始图像存 放的路径 需要处理的 .京始图像样 本数m 图像的宽和 孙」 需要处理的 原始图像的 数据类型旧 fWojeSl Ana^is | Sghe湛 | IHS | BaMile Lisi off dat \ n1e$ |data_l$t 处理后的原 始图壕存放 的路径丫 List/t aux data files Number olwnages 需要处理的 原始图像名 称Mxt文件L 处理中产生 的文件名， 可以随便输 a type I Byte dalajB bits/pixel 0 - 255) ory rnput ff data_list. txt -记事本 文件(E)编辑(E)格式(Q)查看Of)帮助(H) M0D13A2.^2901001 .h23u02.905 -2008279062503 ? hdf . 1_kn_16 ?days_NDUI MOD13A2.A20Ol017.^23ua2,O?,2Oa6359a3l5OS.hdf ￡k?36ZdW^NDUI I40D13A2.A2001033.h23u02.005.2008272005854.hdR6 days