本发明涉及一种转换m天遥感数据到月时间分辨率的平滑方法,属于遥感技术领域。
背景技术:
凡是纪录各种地物电磁波大小的胶片或相片,都称为遥感影像。遥感影像主要是指航空相片和卫星相片。遥感探测器按一定的时间周期重复采集数据,这种重复周期,又称回归周期。它是由飞行器的轨道高度、轨道倾角、运行周期、轨道间隔、偏栘系数等参数所决定。这种重复观测的最小时间间隔称为时间分辨率(temporalresolution)。而时间分辨率不仅是表示遥感影像采样间隔,还可以表示数据产品时间间隔的一项性能指标。比如:中分辨率成像光谱仪(modis)系列影像搭载的上午星(terra卫星)和下午星(aqua卫星)的重放周期为一天四次,而对应出产的很多modis产品时间分辨率则为8天。
实际上,遥感影像数据时间分辨率通常均不是1月,而一年12个月每月的天数也为28-31天不等。因此,月时间分辨率的计算作为遥感影像数据应用中不可或缺的一环,其计算方法直接影响到研究结果的精确度。现有的月时间分辨率采用的是绝对平均值法,这种方法获得的月时间分辨率不够平滑。
技术实现要素:
针对上述不足,本发明提供一种能获得更加平滑的月时间分辨率的转换m天遥感数据到月时间分辨率的平滑方法。
本发明的转换m天遥感数据到月时间分辨率的平滑方法,所述方法包括:
步骤一:采集遥感影像数据,时间分辨率为m天,一年共获得n个时间点,作为样本;m为大于1且小于31的正整数;
步骤二:根据前一个月最后一个样本对本月样本的影响,获得本月的遥感影响数据的月时间分辨率。
优选的是,
式中,
优选的是,所述ν为:
在本月内本月第一个样本之前的无样本天数与m的比值。
本发明还提供一种转换m天遥感数据到月时间分辨率的平滑方法,所述方法包括:
步骤一:采集遥感影像数据,时间分辨率为m天,一年共获得n个时间点,作为样本;m为大于1且小于31的正整数;
步骤二:根据下一个月第一个样本对本月样本的影响,获得本月的遥感影响数据的月时间分辨率。
优选的是,
式中,
优选的是,所述μ为:
本月最后一个样本之后的无样本天数与m的比值。
本发明还提供一种转换m天遥感数据到月时间分辨率的平滑方法,所述方法包括:
步骤一:采集遥感影像数据,时间分辨率为m天,一年共获得n个时间点,作为样本;m为大于1且小于31的正整数;
步骤二:根据前一个月最后一个样本和下一个月第一个样本和对本月样本的影响,获得本月的遥感影响数据的月时间分辨率。
优选的是,
式中,
优选的是,
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天数据所属的当月时,其可用下式表示:
式中,
步骤3:在步骤2的基础上,除了要考虑分布在本月的样本值外,还要考虑前一个月最后一个样本对本月样本的影响,即:采用前接平均值法获得本月的遥感影响数据的月时间分辨率:
因为两个样本日期相差8天,所以本步骤是将前一个月最后一个样本均分为8份,前一个月最后一个样本在本月所占的权重为本月第一个样本之前的无样本天数比上天数8天:
式中,
步骤4:在步骤2的基础上,除了要考虑分布在本月的样本外,还要考虑次月第一个样本值对本月的影响,即:采用后续平均值法获得本月的遥感影响数据的月时间分辨率:
因为两个样本日期相差8天,所以此方法是将次月第一个样本均分为8份,次月第一个样本在本月所占的权重为本月最后一个样本之后的无样本天数比上天数8天,
式中,
步骤5:本步骤不仅考虑到次月月初的样本对本月的影响,同时考虑到本月月末的样本对次月的影响,本实施方式所选数据样本时间分辨率为8天,每一个样本前后都有7天无样本日,所以一个样本所占时间跨度为15天。将每月月初与月末的样本都分为15份,本月月末的样本在本月所占的权重为此样本分布在本月的天数比上所有的跨度天数15天,本月月末的样本在次月所占的权重为此样本分布在次月得天数比上所有的跨度天数15天;同理,次月月初的样本在本月所占的权重为此样本分布在本月的跨度天数比上所有的跨度天数15天,次月月初的样本值在次月所占的权重为此样本分布在次月的天数比上总的跨度天数15天,即:采用综合平均值方法获得本月的遥感影响数据的月时间分辨率:
式中,
x为前一个月最后一个样本分布在本月的天数与2m-1的比值,即:本月第一个样本之前无样本天数与15的比值;
y为本月第一个样本分布在本月的天数与2m-1的比值,即:本月第一个样本之前无样本天数加上8天后与15的比值;
z为本月最后一个样本分布在本月的天数与2m-1的比值,即:本月最后一个样本之后无样本天数加上8天后与15的比值;
w为下一个月第一个样本分布在本月的天数与2m-1的比值,即:本月最后一个样本之后无样本天数与15的比值;
平年和闰年二月天数分别为28天、29天,因此要分别列出其每个月的遥感影像数据求取月时间分辨率。其中,平年从1月到12月的计算公式如下:
对应着,闰年从1月到12月的平滑值计算公式可以分别表达为:
本实施方式是基于时间分辨率为8天的遥感影像数据,提供一种更加平滑的月时间分辨率获取方法,本实施方式的方法同样适用于任何时间分辨率为非1天且小于1月的遥感影像数据。
虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其他所述实施例中。