[发明专利]一种MODIS遥感积雪信息的处理方法和装置

权利要求书

1.一种MODIS遥感积雪信息的处理方法，其特征在于，包括：

基于数字高程模型数据，将预报流域划分为若干个自然流域单元，以及进行高程带划分，将自然流域单元和高程带叠置，得到预报流域的计算单元；

获取MODIS反演雪深信息，将MODIS反演雪深信息与所述计算单元在空间上叠置，得到预报流域内各计算单元对应网格的雪深信息；

对所述雪深信息的有效性进行判别，对具有无效雪深信息的网格比例超过设定阈值的计算单元中的雪深信息进行插补处理，计算处理后的所有计算单元内积雪区域的平均积雪深度；

对计算单元内积雪区域的平均积雪深度进行时间尺度的插值计算，得到计算单元积雪区域逐日的雪深信息。

2.根据权利要求1所述的一种MODIS遥感积雪信息的处理方法，其特征是：所述对所述雪深信息的有效性进行判别，对具有无效雪深信息的网格比例超过设定阈值的计算单元中的雪深信息进行插补处理，包括：

逐一遍历每个计算单元内以网格数据形式呈现的积雪深度信息，判别并标识计算单元内每个网格积雪深度数据的有效性，统计无效雪深信息的网格数在计算单元内的占比，若无效雪深信息的网格数大于总网格*设定的阈值，则将该计算单元定义为待处理计算单元，否则为正常计算单元；

结合地形和历史同期积雪覆盖信息对待处理单元内无效积雪数据进行插补。

3.根据权利要求2所述的一种MODIS遥感积雪信息的处理方法，其特征是：所述结合地形和历史同期积雪覆盖信息对待处理单元内无效积雪数据进行插补，包括：

依据待处理单元各网格地形和历史积雪信息，建立单元内高程、坡度与积雪分布的函数关系：H＝a₁×(Z-Z_b)+a₂×(S-S_b)+b，其中，H为网格雪深；Z为网格高程，Z_b为基准高程；S为网格坡度，S_b为基准坡度，a₁，a₂，b分别为函数中的常数，

将待处理单元内具有无效雪深信息的网格的坡度、高程代入函数关系式，计算得到待处理单元各网格的积雪深度。

4.根据权利要求1所述的一种MODIS遥感积雪信息的处理方法，其特征是：计算处理后的计算单元内积雪区域的平均积雪深度，包括：

计算计算单元内积雪深度大于0的网格的总积雪深度，并除以雪深大于0的网格总数，得到计算单元内积雪区域内的平均雪深。

5.根据权利要求1所述的一种MODIS遥感积雪信息的处理方法，其特征是：所述计算单元内积雪区域逐日的雪深信息，包括：

所述计算单元中某网格的第n日积雪深度H_n计算公式为：H_n＝H_n-1-H_m，n+H_s，n，其中，H_n-1为有积雪信息日第n-1日的该网格积雪深度，H_m，n为第n日的融雪深度，H_s，n为第n日的降雪深度，若第n日的气温超过临界温度，H_s，n为0，若第n日的气温低于临界温度，H_m，n为0；

其中，C_s为融雪径流系数；α为度日因子；ΔT是计算单元与基准气象站因高程不同而产生的气温调整值；T_n为第n日实测气温；S_n为积雪覆盖率，ρ_w、ρ_s分别为水、雪密度；

计算所述计算单元中所有网格的第n日的积雪深度，进而计算得到所述计算单元内第n日积雪区域雪深信息，进而得到所述计算单元内积雪区域逐日的雪深信息。

6.根据权利要求5所述的一种MODIS遥感积雪信息的处理方法，其特征是：若前一日的积雪深度H_n-1小于第n日的融雪深度H_m，n，则该网格积雪深度H_n为0。