去除大气噪声影响重建遥感监测数据集的方法技术

本发明专利技术涉及一种去除大气噪声影响重建遥感监测数据集的方法，首次利用移动窗口来降低噪声点极值对数字信号波形的影响；通过移动窗口对已经降低了极值影响的时序数据进行多项式拟合；选择迭代运算实现最好拟合效果；整体上，遥感数据重建原理简单，很好的保留了原有数字信号的波形特点，突出了数字信号连续变化地趋势。有效的解决了数字信号处理时平滑过度的问题。首次采用先降低数字信号中噪声点对函数拟合效果的影响；二次多项式函数更好地贴合原数据的变化趋势和曲线走向，函数拟合参数动完成，解决了需要试验确定最优运算参数的问题，现实了全程自动化处理。经试验，结果数据与研究区地表实际NDVI信息的相关性强，具有指示性和代表性。

Method for removing atmospheric noise and reconstructing remote sensing monitoring data set

The invention relates to a method for removing atmospheric noise affecting the reconstruction method of remote sensing monitoring data sets, for the first time by using moving window to reduce the noise impact on the extreme digital signal waveform; fitting the time series data have reduced the extreme influence through the moving window; to achieve the best fitting effect of selection of iterative computation; on the whole, the remote sensing data reconstruction the principle is simple, it is good to retain the original waveform characteristics of the digital signal to digital signal to highlight the trend of continuous change. It effectively solves the problem of over smoothing in digital signal processing. For the first time is used to reduce the noise impact point digital signal on the fitting effect function; the two trend curve and polynomial function to better fit the original data to the fitting parameters, function completed, solved the experiment to determine the optimal operation parameters of the problem, the whole process automation a reality. The results show that the correlation between the data and the actual NDVI information of the study area is strong, which is indicative and representative.

【技术实现步骤摘要】

：本专利技术涉及一种去除大气噪声影响重建遥感监测数据集的方法，尤其是利用移动窗口进行多项式拟合迭代运算滤除噪声来重建时间序列数据集。
技术介绍
：现有常用于地表遥感监测的卫星序列数据产品有几十种，时间序列较长的数字信号数据主要为NOAA/AVHRR、SPOT/VEGETATION、TERRA(AQUA)/MODIS等卫星传感器获得的NDVI数据集。由于受大气状况、太阳高度角和地表的复杂状况等方面的影响，这些数据集本身存在着很多噪声，信号源出现较多、较大的异常值，这破坏了原有地表信号连续性的特点，对地表遥感监测的进一步研究造成了很大困难。时间序列数据的滤波重建目的在于利用多元统计和数值分析方法，参照地表植被的季节/年度变化趋势，来插值补充受噪声影响发生异常的数值，优化数字信号数据。CihlarJ(CihlarJ，LyH，LiZ，1997，MultitemporalAVHRRDataSetsforLandBiosphereStudies，RemoteSensingofEnvironment，60(1):35-57)、KobayashiH(KobayashiH，DyeDG，2005，AtmosphericconditionsformonitoringtheLong-termDynamicsintheAmazonUsingNormalizedDifferenceVegetationIndex，RemoteSensingofEnvironment，97(4):519-525)等提出的数据序列重建方法使用广泛。江东等(江东，付晶莹，黄耀欢，庄大方，2011，地表环境参数时间序列重构的方法与应用分析，地球信息科学学报，13(4):439-446)提出的地表环境参数时间序列重构的方法，为有关研究提供进一步可以分析计算的基础数据。目前，依据核心算法原理，时间序列数据的重建方法可以分为以下几大类：阈值法、滤波法、函数拟合法、综合方法等。耿丽英等(耿丽英，马明国，2014，长时间序列NDVI数据重建方法比较研究进展，遥感技术与应用，29(2):362-368)公开了长时间序列NDVI数据重建方法，并总结归纳了现有数字信号滤波重建方法，指出数字信号滤波效果主要受噪声点的影响，不同的噪声点对重建效果也有着不同影响。所有重建方法的去噪效果是不同的，并指出FFT(傅里叶变换)和HANTS的保真性较弱。梁守真等(梁守真，施平，邢前国，2011，MODISNDVI时间序列数据的去云算法比较，国土资源遥感，88(1):33-36)公开了MODISNDVI时间序列数据的去云方法，以山东省MODISNDVI时间序列数据作为试验数据，从不同角度比较了3种数字信号滤波重建能力，指出HANTS算法和S-G算法会改变几乎所有像元数值，得到的一个比较平滑的时间序列曲线，滤波算法的计算参数没有统一标准，需要多次试验才能确定最优参数。综上，现有重建方法存在问题有以下几个方面：一是数字信号重建过程中，平滑过度，失真明显；二是噪声点降低了最终数字信号滤波处理的效果，噪声点影响不可忽略。
技术实现思路
：本专利技术的目的就在于针对上述现有技术的不足，为解决数字信息滤波去噪过程中出现的诸如计算过程复杂，平滑过度和存在人工设定运算参数等问题而提供一种利用移动窗口进行多项式拟合去除大气噪声影响重建遥感监测数据集的方法。本专利技术的目的是通过以下技术方案实现的：去除大气噪声影响重建遥感监测数据集的方法，包括以下步骤：A、选用SpotVegetation提供的NDVI数据集；B、将影像DN值转化成NDVI数据，根据研究区范围和空间分辨率对图像进行裁剪；C、重复上述A-B步骤，得到NDVI遥感图像数据集；D、在R语言平台上编写算法，设定移动窗口大小为5；5为某待确定点的前2期与后2期为窗口内统计数据的实际范围，统计仅限窗口内数值的平均值与标准差；E、根据窗口内的实际统计值，划定信号源数据的判定范围-S≤Ndvi≤+S；F、在NDVI序列数据集上移动窗口完成整个信号源数据的极值判定与处理；G、选择多项式拟合函数为二次多项式，设定移动窗口大小为5对窗口内极值处理后的信号源数值进行二次多项式拟合；H、设定迭代次数，对上一步拟合后的序列数据进行2次拟合，得到函数拟合的最终结果；I、根据计算初始设定的循环次数，判定是否结束本次计算；J、对结果数据赋投影信息，投影参数与原数据地图参数一致；K、完成时间序列遥感数据集的重建，保存数据。步骤E所述的判定范围为：-S≤Ndvi≤+S；正常值-S≤Ndvi≤+S；偏高+S<Ndvi<+1.5S；显著偏高+1.5S≤Ndvi<+2S；异常偏高+2S≤Ndvi<+∞；偏低-S<Ndvi<-1.5S；显著偏低-1.5S≤Ndvi<-2S；异常偏低-2S≤Ndvi<-∞。设定信号源数据大于+1.5S或小于-1.5S为噪声点异常值，将噪声点异常值以窗口内数字信号点的均值±1.5S代替。有益效果：首次利用移动窗口来消除极大值与极小值对整体波形拟合的影响；对窗口内消除极值影响的时序数据进行多项式拟合；选择迭代运算实现最好拟合效果。有效的解决了有关算法对数字信号处理时平滑过度的问题，首次采用先降低数字信号中噪声点对函数拟合效果的影响，二次多项式函数更好地贴合原数据的变化趋势和曲线走向，函数拟合参数由算法自动完成，解决了需要试验确定最优运算参数的问题，现实了全程自动化运算处理。经试验，结果数据与研究区地表实际NDVI信息的相关性强，具有指示性和代表性。附图说明：图1为去除大气噪声影响重建遥感监测数据集的方法流程图；图2为遥感卫星数据头文件编辑图；图3为不同植被类型的NDVI时序曲线重建图；图4为长时间序列的数字信号重建图；具体实施方式下面结合附图和实施例作进一步的详细说明：(1)地表植被生长有着明显的季节性、连续性特点，进行数字信号滤波处理时，数据的拟合效果受噪声影响明显。由于噪声的存在，这种连续性时常被异常值打断，所以本发明首先对数字信号中的噪声点进行了分析处理。(2)二次多项式函数对有着明显变化趋势的数据序列有着很好的拟合效果。经过降低噪声点异常值影响后的数字信号序列采用二次多项式拟合，通过2-3次迭代算法对原始数据序列集进行拟合。(3)综合以上，构建利用移动窗口进行多项式拟合迭代运算的滤波算法以重建数字信号序列。本模型建设过程主要包括噪声极值处理、多项式拟合、迭代算法等步本文档来自技高网...

【技术保护点】
一种去除大气噪声影响重建遥感监测数据集的方法，其特征在于：包括以下步骤：A、选用Spot Vegetation提供的NDVI数据集；B、将影像DN值转化成NDVI数据，根据研究区范围和空间分辨率对图像进行裁剪；C、重复上述A‑B步骤，得到NDVI遥感图像数据集；D、在R语言平台上编写算法，设定移动窗口大小为5；5为某待确定点的前2期与后2期为窗口内统计数据的实际范围，统计仅限窗口内数值的平均值与标准差；E、根据窗口内的实际统计值，划定信号源数据的判定范围‑S≤Ndvi≤+S；F、在NDVI序列数据集上移动窗口完成整个信号源数据的极值判定与处理；G、选择多项式拟合函数为二次多项式，设定移动窗口大小为5对窗口内极值处理后的信号源数值进行二次多项式拟合；H、设定迭代次数，对上一步拟合后的序列数据进行2次拟合，得到函数拟合的最终结果；I、根据计算初始设定的循环次数，判定是否结束本次计算；J、对结果数据赋投影信息，投影参数与原数据地图参数一致；K、完成时间序列遥感数据集的重建，保存数据。

【技术特征摘要】
1.一种去除大气噪声影响重建遥感监测数据集的方法，其特征在于：包括以下步骤：
A、选用SpotVegetation提供的NDVI数据集；
B、将影像DN值转化成NDVI数据，根据研究区范围和空间分辨率对图像进行裁剪；
C、重复上述A-B步骤，得到NDVI遥感图像数据集；
D、在R语言平台上编写算法，设定移动窗口大小为5；
5为某待确定点的前2期与后2期为窗口内统计数据的实际范围，统计仅限窗口内数值
的平均值与标准差；
E、根据窗口内的实际统计值，划定信号源数据的判定范围-S≤Ndvi≤+S；
F、在NDVI序列数据集上移动窗口完成整个信号源数据的极值判定与处理；
G、选择多项式拟合函数为二次多项式，设定移动窗口大小为5
对窗口内极值处理后的信号源数值进行二次多项式拟合；
H、设定迭代次数，对上一步拟合后的序列数据进行2次拟合，得到函数拟合的最终结
果；

【专利技术属性】
技术研发人员：李晓东，曾发梁，姜琦刚，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22