HANTS时间序列谐波分析法.docx

)——时间序列HANT(theHarmonicanalysisoftimeseries谐波分析法时间序列谐波分析法(HarmonicAnalysisofTimeSeries,HANTS是平滑和滤波两种方法的综合,它能够充分利用遥感图像存在时间性和空间性的特点,将其空间上的分布规律和时间上的变化规律联系起来。 时间序列谐波分解法进行影像重构时充分考虑了植被生长周期性和数据本身的双重特点, 能够用代表不同生长周期的植被频率曲线重新构建时序 NDVI影像,真实反映植被的周期性变化规律。时间序列谐波分析法是对快速傅立叶变换的改进, 它不仅可以去除云污染点,而且对时序图像的要求不象快速傅立叶变换( FFT那么严格,它可以是不等时间间隔的影像。因此同快速傅立叶变换相比, HANTSS频率和时间系列长度的选择上具有更大的灵活性。时间序列谐波分析法进行时序影像的重构也是基于云对NDVI的负值影响,但是它与最大值去除云污染的影响是两个完全不同的方法。它是首先通过傅立叶变换得到非零频率的振幅和相位,然后将所有的点进行最小二次方拟合。通过观测资料与拟合曲线的比较,对于那些明显低于拟合曲线的点被作为云污染点通过把它们的权重赋为零而拒绝参与曲线的拟合。 建立在剩余点上进行新的曲线拟合,通过这种反复进行的迭代过程实现图像的重构。HANTS勺核心算法是最小二乘法和傅立叶变换,通过最小二乘法的迭代拟合去除时序 NDVI值中受云污染影响较大的点,借助于傅立叶在时间域和频率域的正反变换实现曲线的分解和重构,从而达到时序遥感影像去云重构的目的。采用时间序列谐波分析法(HANTS可以对时间谱数据进行平滑。其核心算法是傅立叶变换和最小二乘法拟合,即将时间谱数据分解成有限个谐波(正弦波或余弦波) ,从中选取若干个能反映影像时序特征的谐波进行叠加,达到重构时序数据的目的。(左丽君MODIS/NDVI和MODIS/EVI在耕地信息提取中的应用及对比分析 )原始NDVI曲线与HANTS滤波后曲线Hants软件操作步骤1:projict设置输入输出、图像大小、波段等。data_list为时间序列数据的文件列表,需和数据文件放同一个目录下(数据为裸数据格式,envi标准格式即可),格式如下图所示(IMG只是加的一个标识,用于其他处理,并不是文件的格式)。口叵)区I口叵)区IHelp|Hefp需要处理的原始图像存放的路径需要处理的放的路径处理后的原始图塚存放的路径衣:::二1||需要处理的原始图像名称(X文件)处理后的原始图塚存放的路径卩需妾处理的原始图像名称(皿文件)壽妾处理的原始图像样本数卩要处理的原始图像样Cdata_list-txt-记事本文件(E)绸辑(E)格式(Q)査帮助迫)_kn_16_daps_HDUI号s)DUI[krCKdayjDUIZkin2l6Zdays2NDUI2kn_16_days_NDUI[\h23u02*<^kn^16_days_HDUIM0D13A2•A2£!01GF7・h23u©2・Ot6・2B06359031508・h(ll:・FllaCril6[(l2ys[l1DUIM0D13A2-A2S01033・h23u02・>*ZdayCHDUIMOD13A2-A2D01049・・tl(iF・4一1<皿一16一(1包9$一时0111M0D13A2,-,,lZkP)Zl6^,^kn^ldZdaps^HDUIM0D13A2-A2001129・h23u024OSS«~一一~(