[发明专利]地表反照率产品的生成方法及系统

权利要求书

1.一种地表反照率产品的生成方法，其特征在于，包括步骤：

B：根据太阳天顶角、观测天顶角和相对方位角对MODIS地表反射率数据、MODIS表观反射率数据和MODIS地表反照率数据进行网格划分；

C：对POLDER BRDF数据进行优选，按照地表类型对经过优选的POLDER BRDF数据进行分类；

D：根据配套光谱数据模拟得到模拟POLDER BRDF数据和模拟MODIS地表反射率数据，建立所述模拟POLDER BRDF数据和所述模拟MODIS地表反射率数据之间的多元线性转换关系，根据所述多元线性转换关系和经过分类的POLDER BRDF数据，计算得到样本MODIS地表反射率数据；

E：根据经过分类的POLDER BRDF数据，得到POLDER反照率数据，转换所述POLDER反照率数据得到宽波段地表反照率数据；

F：根据所述样本MODIS地表反射率数据和所述宽波段地表反照率数据建立第一转换模型，根据所述第一转换模型和经过网格划分的MODIS地表反射率数据，得到第一地表反照率数据；

G：根据所述POLDER BRDF数据得到地表反射特性参数，根据所述地表反射特性参数和大气状态参数计算得到样本MODIS表观反射率数据；

H：根据所述样本MODIS表观反射率数据和所述宽波段地表反照率数据建立第二转换模型，根据所述第二转换模型和经过网格划分的MODIS表观反射率数据，得到第二地表反照率数据；

I：整合所述第一地表反照率数据、所述第二地表反照率数据和经过网格划分的MODIS地表反照率数据，生成时空连续的地表反照率产品。

2.如权利要求1所述的方法，其特征在于，在所述步骤B之前还包括步骤A：输入POLDER BRDF数据、MODIS地表反射率数据、MODIS表观反射率数据、MODIS地表反照率数据和配套光谱数据。

3.如权利要求1所述的方法，其特征在于，所述步骤I具体包括步骤：

I1：对所述第一地表反照率数据和所述第二地表反照率数据进行插值弥补；

I2：对所述第一地表反照率数据、所述第二地表反照率数据和经过网格划分的MODIS地表反照率数据进行归一化处理；

I3：对归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据进行数据融合，得到融合数据；

I4：对所述融合数据进行平滑处理；

I5：对平滑处理后的融合数据进行反归一化处理，得到所述时空连续的地表反照率产品。

4.如权利要求3所述的方法，其特征在于，所述反归一化处理的公式如下：

ξ=ξ′σ‾-μ‾]]>

其中，ξ是所述平滑处理后的融合数据，ξ是所述时空连续的地表反照率产品，是归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据的均值，是归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据的消除随机误差后的标准差的平均值。

5.如权利要求1所述的方法，其特征在于，所述步骤C具体包括步骤：

C1：从所述POLDER BRDF数据中剔除由于云和气溶胶影响以及地表状态发生变化而不满足二向反射模型的数据；

C2：对经过所述步骤C1处理后的POLDER BRDF数据，按照植被、冰雪和裸地三种地表类型进行分类。