1.一种极地反照率遥感数据的重构方法，其特征在于，包括以下的步骤：

步骤S1：获取南北极区域反照率、云光学厚度和太阳天顶角的遥感数据，统一所述遥感数据的空间分辨率，得到反照率数值序列、云光学厚度数值序列及太阳天顶角数值序列；

步骤S2：根据空缺像元相邻时间和空间的有效像元的反照率数值的均值填补空缺像元的反照率数值，空缺像元的填补函数为：其中，(m,n)代表空间上第(m,n)个像元的二维序号，所述像元具有固定的行列号，i代表时序为i的像元的序号，num表示相邻时空有效像元的数量，yM,N,I表示空缺像元相邻时空的某个有效像元的反照率数值，ym,n,i表示被填补的空缺像元的反照率数值；

步骤S3：将步骤S2中填补的空缺像元分别作为目标像元，选取包括该目标像元的时序上连续的多个像元构成第一像元组，选取与该目标像元时间上相同的有效像元作为对比像元，选取包括该对比像元的时序上连续的多个像元构成第二像元组，所述第一像元组与所述第二像元组的时序相同，计算所述第一像元组与所述第二像元组的反照率数值时序数据的相似性，若所述第一像元组与所述第二像元组的反照率数值时序数据满足相似性条件，则使用所述对比像元的反照率数值替换所述目标像元的反照率数值；

步骤S4：分别选取所述空缺像元所在空间位置上的所有像元的反照率数值时序数据构成填补数据序列[Yi]，利用Whittaker算法对所述填补数据序列Yi进行滤波，生成拟合数据序列[Zi]；Yi表示时序为i的像元的反照率数值，Zi表示对应Yi的拟合数值；拟合数据序列[Zi]与填补数据序列[Yi]的拟合关系式如下：‑1

Zi＝(I+λDz'Dz) Yi

其中，I为单位矩阵，λ为参数，Dz表示zi的二阶差分，Dz＝zi‑3zi‑1+3zi‑2‑zi‑3，Dz′表示Dz对zi的导数；

步骤S5：对于每个像元，计算权重W，权重W的计算公式为：

利用Whittaker算法构建拟合指数F，拟合指数F的计算公式为：

对所述拟合指数F进行迭代计算，当拟合指数F符合条件Fk‑1≥fk≤fk+1时退出迭代，计算得到该像元的重构的反照率数值 其中k表示迭代次数，k≥2，步骤S6：判断空缺像元是否对应被冰雪覆盖的区域，若是则利用所述空缺像元对应的云光学厚度数值和太阳天顶角数值对该空缺像元的所述重构的反照率数值进行修正，得到修正的反照率数值；对该空缺像元的所述重构的反照率数值进行修正的计算公式如下：其中， τ、S分别表示修正的反照率数值、重构的反照率数值、云光学厚度数值、太阳天顶角数值，a、b、c、d为回归系数。

2.如权利要求1所述的重构方法，其特征在于：所述步骤S3中，衡量所述第一像元组与所述第二像元组的反照率数值时序数据的相似性使用皮尔逊相关系数，其计算公式为：其中，X表示第一像元组的反照率数值时序数据，Y表示第二像元组的反照率数值时序数据，ρ(X,Y)表示第一像元组与第二像元组的反照率数值时序数据的皮尔逊相关系数，E表示数学期望，D表示方差，con(X,Y)表示协方差。

3.如权利要求2所述的重构方法，其特征在于：ρ(X,Y)大于0.9时，所述第一像元组与所述第二像元组的反照率数值时序数据满足相似性条件。

4.如权利要求1所述的重构方法，其特征在于：λ设置为2。

5.如权利要求1所述的重构方法，其特征在于：a、b、c、d分别设置为‑0.0491243、

1.06756、0.0217075、0.017950。

6.一种极地反照率遥感数据的重构装置，其特征在于，包括：

数据获取模块，用于获取南北极区域反照率、云光学厚度和太阳天顶角的遥感数据，统一所述遥感数据的空间分辨率，得到反照率数值序列、云光学厚度数值序列及太阳天顶角数值序列；

数值填补模块，用于根据空缺像元相邻时间和空间的有效像元的反照率数值的均值填补空缺像元的反照率数值，空缺像元的填补函数为：其中，(m,n)代表空间上第(m,n)个像元的二维序号，所述像元具有固定的行列号，i代表时序为i的像元的序号，num表示相邻时空有效像元的数量，yM,N,I表示空缺像元相邻时空的某个有效像元的反照率数值，ym,n,i表示被填补的空缺像元的反照率数值；

数值替换模块，用于将数值填补模块中填补的空缺像元分别作为目标像元，选取包括该目标像元的时序上连续的多个像元构成第一像元组，选取与该目标像元时间上相同的有效像元作为对比像元，选取包括该对比像元的时序上连续的多个像元构成第二像元组，所述第一像元组与所述第二像元组的时序相同，计算所述第一像元组与所述第二像元组的反照率数值时序数据的相似性，若所述第一像元组与所述第二像元组的反照率数值时序数据满足相似性条件，则使用所述对比像元的反照率数值替换所述目标像元的反照率数值；

数值拟合模块，用于分别选取所述空缺像元所在空间上的所有像元的反照率数值时序数据构成填补数据序列[Yi]，利用Whittaker算法对所述填补数据序列Yi进行滤波，生成拟合数据序列[Zi]；Yi表示时序为i的像元的反照率数值，Zi表示对应Yi的拟合数值；拟合数据序列[Zi]与填补数据序列[Yi]的拟合关系式如下：‑1

Zi＝(I+hDz'Dz) Yi

其中，I为单位矩阵，λ为参数，Dz表示zi的二阶差分，Dz＝zi‑3zi‑1+3zi‑2‑zi‑3，Dz′表示Dz对zi的导数；

迭代计算模块，用于对于每个像元，计算权重W，权重W的计算公式为：

利用Whittaker算法构建拟合指数F，拟合指数F的计算公式为：

对所述拟合指数F进行迭代计算，当拟合指数F符合条件Fk‑1≥fk≤fk+1时退出迭代，计算得到该像元的重构的反照率数值 其中k表示迭代次数，k≥2，修正模块，用于判断空缺像元是否对应被冰雪覆盖的区域，若是则利用所述空缺像元对应的云光学厚度数值和太阳天顶角数值对该空缺像元的所述重构的反照率数值进行修正，得到修正的反照率数值；对该空缺像元的所述重构的反照率数值进行修正的计算公式如下：其中， τ、S分别表示修正的反照率数值、重构的反照率数值、云光学厚度数值、太阳天顶角数值，a、b、c、d为回归系数。

7.一种计算机设备，其特征在于，包括处理器和存储器，存储在存储器上并可在处理器上运行的极地反照率遥感数据的重构程序，该重构程序被处理器执行时执行如权利要求1‑

5中任一项所述的重构方法的步骤。