转换M天遥感数据到月时间分辨率的平滑方法技术

本发明专利技术提供一种能获得更加平滑的月时间分辨率的转换M天遥感数据到月时间分辨率的平滑方法，属于遥感技术领域。本发明专利技术在原始的求取月时间分辨率的方法基础上，不仅考虑到分布于所求月的样本对本月的影响，还考虑到分布在本月月初、月末的两个样本对相邻月份的影响，以及相邻前一个月最后一个样本对本月的影响、相邻后一个月份的第一个样本对本月的影响；这些影响以权重的形式赋予样本值，综合考虑后续平均和前接平均值，最终得出更加平滑、最接近原始的月时间分辨率。

【技术实现步骤摘要】
转换M天遥感数据到月时间分辨率的平滑方法
本专利技术涉及一种转换M天遥感数据到月时间分辨率的平滑方法，属于遥感

技术介绍
凡是纪录各种地物电磁波大小的胶片或相片，都称为遥感影像。遥感影像主要是指航空相片和卫星相片。遥感探测器按一定的时间周期重复采集数据，这种重复周期，又称回归周期。它是由飞行器的轨道高度、轨道倾角、运行周期、轨道间隔、偏栘系数等参数所决定。这种重复观测的最小时间间隔称为时间分辨率(Temporalresolution)。而时间分辨率不仅是表示遥感影像采样间隔，还可以表示数据产品时间间隔的一项性能指标。比如：中分辨率成像光谱仪(MODIS)系列影像搭载的上午星(terra卫星)和下午星(aqua卫星)的重放周期为一天四次，而对应出产的很多MODIS产品时间分辨率则为8天。实际上，遥感影像数据时间分辨率通常均不是1月，而一年12个月每月的天数也为28-31天不等。因此，月时间分辨率的计算作为遥感影像数据应用中不可或缺的一环，其计算方法直接影响到研究结果的精确度。现有的月时间分辨率采用的是绝对平均值法，这种方法获得的月时间分辨率不够平滑。
技术实现思路
针对上述不足，本专利技术提供一种能获得更加平滑的月时间分辨率的转换M天遥感数据到月时间分辨率的平滑方法。本专利技术的转换M天遥感数据到月时间分辨率的平滑方法，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据前一个月最后一个样本对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。优选的是，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，ν为本月之前的一个月最后一个样本在本月所占的权重。优选的是，所述ν为：在本月内本月第一个样本之前的无样本天数与M的比值。本专利技术还提供一种转换M天遥感数据到月时间分辨率的平滑方法，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据下一个月第一个样本对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。优选的是，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，μ为下一个月第一个样本在本月所占的权重。优选的是，所述μ为：本月最后一个样本之后的无样本天数与M的比值。本专利技术还提供一种转换M天遥感数据到月时间分辨率的平滑方法，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据前一个月最后一个样本和下一个月第一个样本和对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。优选的是，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，x为前一个月最后一个样本在本月所占的权重，y为本月第一个样本在本月所占的权重，z为本月最后一个样本在本月所占的权重，w为下一个月第一个样本在本月所占的权重。优选的是，x为前一个月最后一个样本分布在本月的天数与2M-1的比值；y为本月第一个样本分布在本月的天数与2M-1的比值；z为本月最后一个样本分布在本月的天数与2M-1的比值；w为下一个月第一个样本分布在本月的天数与2M-1的比值。本专利技术的有益效果在于，本专利技术在原始的求取月时间分辨率的方法基础上，不仅考虑到分布于所求月的样本对本月的影响，还考虑到分布在本月月初、月末的两个样本对相邻月份的影响，以及相邻前一个月最后一个样本对本月的影响、相邻后一个月份的第一个样本对本月的影响；这些影响以权重的形式赋予样本值，综合考虑后续平均和前接平均值，最终得出更加平滑、最接近原始的月时间分辨率。附图说明图1为本专利技术具体实施方式的流程示意图；图2为本专利技术具体实施方式中在平年对叶面积指数LAI平滑后的效果图，其中原点表示平滑前的，线段表示平滑后的；图3为本专利技术具体实施方式中在闰年对叶面积指数LAI平滑后的效果图，其中原点表示平滑前的，线段表示平滑后的。具体实施方式下面结合附图和具体实施例对本专利技术作进一步说明，但不作为本专利技术的限定。结合图1说明本实施方式，本实施方式M＝8，本实施方式的转换8天遥感数据到月时间分辨率的平滑方法，包括：步骤1：采集遥感影像数据，本实施方式采用中国地区叶面积指数LAI数据，时间分辨率为8天，一年共获得46个时间点，作为样本；步骤2：采用绝对平均值法获取遥感影响数据的月时间分辨率：在仅考虑8天数据所属的当月时，其可用下式表示：式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为8天，i＝1+8n(n＝0，1，2，3…，n<47)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数；步骤3：在步骤2的基础上，除了要考虑分布在本月的样本值外，还要考虑前一个月最后一个样本对本月样本的影响，即：采用前接平均值法获得本月的遥感影响数据的月时间分辨率：因为两个样本日期相差8天，所以本步骤是将前一个月最后一个样本均分为8份，前一个月最后一个样本在本月所占的权重为本月第一个样本之前的无样本天数比上天数8天：式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，v为本月之前的一个月最后一个样本在本月所占的权重。步骤4：在步骤2的基础上，除了要考虑分布在本月的样本外，还要考虑次月第一个样本值对本月的影响，即：采用后续平均值法获得本月的遥感影响数据的月时间分辨率：因为两个样本日期相差8天，所以此方法是将次月第一个样本均分为8份，次月第一个样本在本月所占的权重为本月最后一个样本之后的无样本天数比上天数8天，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，μ为下一个月第一个样本在本月所占的权重。步骤5：本步骤不仅考虑到次月月初的样本对本月的影响，同时考虑到本月月末的样本对次月的影响，本实施方式所选数据样本时间分辨率为8天，每一个样本前后都有7天无样本日，所以一个样本所占时间跨度为15天。将每月月初与月末的样本都分为15份，本月月末的样本在本月所占的权重为此样本分布在本月的天数比上所有的跨度天数15天，本月月末的样本在次月所占的权重为此样本分布在次月得天数比上所有的跨度天数15天；同理，次月月初的样本在本月所占的权重为此样本分布在本月的跨度天数比上所有的跨度天数15天，次月月本文档来自技高网...

【技术保护点】
1.一种转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据前一个月最后一个样本对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。

【技术特征摘要】
1.一种转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据前一个月最后一个样本对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。2.根据权利要求1所述的转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月内的样本对应的首尾天数，T为本月样本的总天数，v为本月之前的一个月最后一个样本在本月所占的权重。3.根据权利要求2所述的转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，所述v为：本月第一个样本之前的无样本天数与M的比值。4.一种转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，所述方法包括：步骤一：采集遥感影像数据，时间分辨率为M天，一年共获得N个时间点，作为样本；M为大于1且小于31的正整数；步骤二：根据下一个月第一个样本对本月样本的影响，获得本月的遥感影响数据的月时间分辨率。5.根据权利要求4所述的转换M天遥感数据到月时间分辨率的平滑方法，其特征在于，式中，表示月时间分辨率，DOYi为样本对应的实际日期，两个相邻样本时间间隔为M天，i＝1+M*n(n＝0，1，2，3…，n<N+1)，a和b分别表示本月...

【专利技术属性】
技术研发人员：崔耀平，都义娟，朱连奇，蒋琳，刘素洁，李楠，
申请(专利权)人：河南大学，
类型：发明
国别省市：河南,41