一种时序NDVI数据序列的优化方法技术

本发明专利技术提供一种时序NDVI数据序列的优化方法，该方法包括如下步骤：获取初始的第一时序NDVI数据序列；在限制条件下，采用最佳坡度系数截取法BISE对所述第一时序NDVI数据序列进行优化，得到第二时序NDVI数据序列；采用小波变换法对所述第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列。采用本发明专利技术所述的方法优化初始的时序NDVI数据序列，能有效的去除初始的时序NDVI数据序列中的噪声数据，同时保持了时序NDVI数据的真实性，使得时序NDVI数据能够更为准确的体现被监测目标的真实变化信息。

An optimization method for sequential NDVI data sequences

The present invention provides an optimization method of time series NDVI data sequence, which includes the following steps: acquiring initial first sequence NDVI data sequence; optimizing the first sequence NDVI data sequence by optimum slope coefficient interception method BISE under limited conditions, and obtaining the second sequence NDVI data sequence by wavelet transform; The conversion method optimizes the second time sequence NDVI data sequence and obtains the target time sequence NDVI data sequence. The method of the invention can optimize the initial sequence NDVI data sequence, effectively remove the noise data in the initial sequence NDVI data sequence, and maintain the authenticity of the sequence NDVI data, so that the sequence NDVI data can more accurately reflect the real change information of the monitored target.

【技术实现步骤摘要】
一种时序NDVI数据序列的优化方法
本专利技术涉及地理信息与遥感
，尤其涉及一种时序的NDVI数据序列的优化方法。
技术介绍
NDVI(NormalizedDifferenceVegetationIndex)归一化差分植被指数，也称为生物量指标变化。NDVI是植被生长状态及植被覆盖度的最佳指示因子。NDVI能够准确的反映出地表植被的活力状态和植被季相变化特征。基于这一特性，NDVI现已被广泛的应用于各种变化的监测中，如城市扩张监测、森林扰动监测、植被长势监测、气候变化监测、开采扰动监测等。原始的时序NDVI数据是通过对卫星遥感器获取的遥感数据进行辐射定标、大气校正等计算得到的。遥感器获取遥感数据的过程容易受到云、气溶胶、太阳高度角和地物双向性反射等一系列不良因素的影响，导致遥感器获取的遥感数据中包含了一定量的噪声数据。这也导致在此基础上计算得出的时序NDVI数据中存在一定量的噪声数据。当上述时序NDVI数据用于变化检测时，噪声数据会在一定程度上影响检测结果的真实性，导致检测结果不能体现被检测对象的真实变化，使人们根据检测结果做出错误的判断。常见的时序NDVI数据序列的优化方法，如采用最佳坡度系数截取法BISE优化时序NDVI数据序列，虽能很大程度地去除时序NDVI数据序列中的噪声数据，但优化后的时序NDVI数据序列仍然具有较高的离散性，尤其是在连续噪声处BISE表现不佳；采用Savitzky-Golay滤波法优化时序NDVI数据序列，在去除时序NDVI数据序列中的噪声数据的同时，还可以降低优化后的时序NDVI数据序列中数据的离散性，但优化后的时序NDVI数据序列与真实的数据序列的时相与幅度存在一定误差，目前该误差无法消除。综上所述，如何去除时序NDVI数据序列中的这些噪声数据，并尽可能的保证时序NDVI数据的真实性，是目前我们面临的首要问题。
技术实现思路
为解决上述问题，本专利技术提供了一种时序NDVI数据序列的优化方法，包括如下步骤：获取初始的第一时序NDVI数据序列；在限制条件下，采用最佳坡度系数截取法BISE对第一时序NDVI数据序列进行优化，得到第二时序NDVI数据序列；采用小波变换法对第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列。优选地，获取初始的第一时序NDVI数据序列步骤，包括：获取原始卫星遥感数据；对原始卫星遥感数据进行辐射定标；对定标后的数据进行大气校正；对大气校正后的数据进行配准和几何校正,并对计算得出的NDVI数据按时间先后顺序重新组合，得到第一时序NDVI数据序列。优选地，特定的限制条件的表达式为：其中xi-1、xi、xi+1为在第一时序NDVI数据序列中选取的3个连续的数据,i为大于等于2的整数。优选地，在限制条件下，采用最佳坡度系数截取法BISE对第一时序NDVI数据序列进行优化,得到第二时序NDVI数据序列步骤，包括:数据采集：在第一时序NDVI数据序列中依次选取3个的连续的数据：xi-1、xi、xi+1；数据预处理：根据公式和公式分别计算出xi-1和xi的斜率slop(i-1,i)，以及xi+1和xi的斜率slop(i+1,i)；数据优化：设定一个阈值Ti，当斜率slop(i-1,i)的绝对值|slop(i-1,i)|和斜率slop(i+1,i)的绝对值|slop(i+1,i)|同时大于阈值Ti，且满足时，xi被公式取代，即：优化数据的整理：按时序先后逐一选取第一时序NDVI数据序列中的数据，并对所选取的数据进行优化，优化后的数据按时序由先到后重新组合，得到第二时序NDVI数据序列。优选地，采用小波变换法对第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列步骤，包括：采用预先设定的小波对第二时序NDVI数据序列进行N层分解计算，N为预先设定的分解层数，N≥2；采用预先设定的阈值模式Y，对第1层到第N层的每一层高频系数进行阈值量化处理；根据小波分解的第N层的低频系数和经过量化处理后的第1层到第N层的高频系数，对第二时序NDVI数据序列进行小波重构，得到目标信号波；从目标信号波中提取NDVI数据，将提取到的NDVI数据按时间先后顺序重新组合得到目标时序NDVI数据序列。优选地，预先设定的小波为db7小波；分解层数N具体为：2；特定阈值模式Y设置为：极大极小阈值准则：minimaxi，且设置为软阈值：soft。采用本专利技术的方法优化初始的时序NDVI数据序列，能有效的去除初始的时序NDVI数据序列中的噪声数据，同时保持了时序NDVI数据的真实性，使得NDVI能够更为准确的体现被监测目标的真实变化信息。附图说明图1为本专利技术实施例提供的一种时序NDVI数据序列的优化方法流程图；图2为本专利技术获取初始的第一时序NDVI数据序列步骤的流程图；图3为本专利技术在限制条件下，采用最佳坡度系数截取法BISE对所述第一时序NDVI数据序列进行优化,得到第二时序NDVI数据序列步骤的流程图；图4为采用小波变换法对所述第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列步骤的流程图；图5为检测区域一：中国内蒙宝日希勒镇图6为检测区域二：美国阿拉巴拉契亚马丁镇图7为第一样点的数据优化对比图；图8为第二样点的数据优化对比图；图9为第三样点的数据优化对比图；图10为第四样点的数据优化对比图；图11为采用最佳坡度系数截取法BISE优化后的时序NDVI数据序列得到的变化检测结果；图12为采用Savitzky-Golay滤波法优化后的时序NDVI数据序列得到的变化检测结果；图13为采用本专利技术所记载的新方法优化后的时序NDVI数据序列得到的变化检测结果；图14为采用原始时序NDVI数据序列得到的变化检测结果。具体实施方式下面结合附图和实施例，对本专利技术的技术方案做进一步的详细描述。图1为本专利技术实施例提供的一种时序NDVI数据序列的优化方法流程图。如图1所示，时序NDVI数据序列的优化方法包括以下步骤：S101，获取初始的第一时序NDVI数据序列。具体地，该步骤具体的实施过程，如图2所示，包括：首先，获取原始卫星遥感数据；然后，对该原始卫星遥感数据进行辐射定标；紧接着，对定标后的数据进行大气校正；然后，对大气校正后的数据进行配准和几何校正；最后，根据配准和几何校正后的遥感数据计算遥感数据的NDVI数据并将计算得出的NDVI数据按时间先后顺序重新组合，得到第一时序NDVI数据序列。S102，在限制条件下，采用最佳坡度系数截取法BISE对第一时序NDVI数据序列进行优化,得到第二时序NDVI数据序列。具体地，限制条件的表达式为：其中xi-1、xi、xi+1为在第一时序NDVI数据序列中选取的3个连续的数据，i为大于等于2的整数。该步骤具体的实施过程，如图3所示，包括：首先，进行数据采集，即：在第一时序NDVI数据序列中依次选取3个的连续的数据：xi-1、xi、xi+1；然后，对采集到的数据进行预处理：根据公式和公式分别计算出xi-1和xi的斜率slop(i-1,i)，以及xi+1和xi的斜率slop(i+1,i)；紧接着，对采集到的数据进行优化：设定一个阈值Ti，当斜率slop(i-1,i)的绝对值|slop(i-1,i)|和斜率slop(i+1,i)的绝对值|slop(i+本文档来自技高网...

【技术保护点】
1.一种时序NDVI数据序列的优化方法，其特征在于，包括如下步骤：获取初始的第一时序NDVI数据序列；在限制条件下，采用最佳坡度系数截取法BISE对所述第一时序NDVI数据序列进行优化，得到第二时序NDVI数据序列；采用小波变换法对所述第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列。

【技术特征摘要】
1.一种时序NDVI数据序列的优化方法，其特征在于，包括如下步骤：获取初始的第一时序NDVI数据序列；在限制条件下，采用最佳坡度系数截取法BISE对所述第一时序NDVI数据序列进行优化，得到第二时序NDVI数据序列；采用小波变换法对所述第二时序NDVI数据序列进行优化，得到目标时序NDVI数据序列。2.根据权利要求1所述的方法，其特征在于，所述获取初始的第一时序NDVI数据序列步骤，包括：获取原始卫星遥感数据；对所述原始卫星遥感数据进行辐射定标；对定标后的数据进行大气校正；对大气校正后的数据进行配准和几何校正；对配准和几何校正后的数据进行NDVI计算，并对计算得出的NDVI数据按时间先后顺序重新组合，得到所述第一时序NDVI数据序列。3.根据权利要求1所述的方法，其特征在于，所述限制条件的表达式为：其中xi-1、xi、xi+1为在所述第一时序NDVI数据序列中选取的3个连续的数据，i为大于等于2的整数。4.根据权利要求3所述的方法，其特征在于，所述在限制条件下，采用最佳坡度系数截取法BISE对所述第一时序NDVI数据序列进行优化，得到第二时序NDVI数据序列步骤，包括:数据采集：在所述第一时序NDVI数据序列中依次选取3个的连续的数据：xi-1、xi、xi+1；数据预处理：根据公式和公式分别计算出xi-1和xi的斜率slop(i-1,i)，以及xi+1...

【专利技术属性】
技术研发人员：李晶，杨震，闫萧萧，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京,11