[发明专利]一种检测区域积雪深度的遥感估算方法

权利要求书

1.一种检测区域积雪深度的遥感估算方法，其特征在于：所述遥感估算方法是按下列步骤进行的：

1）确定检测区域的降雪时间

通过查阅气象部门资料获取区域开始降雪的具体时间；

2）获取检测区域的合成孔径雷达数据

在步骤1）确定区域具体降雪时间的基础上，下载覆盖区域降雪前与积雪深度估算时间之间的所有Sentinel-1A原始数据，利用专业软件对原始数据进行预处理，得到N景Sentinel-1A标准格式数据；

3）检测区域合成孔径雷达数据的D-InSAR处理

对步骤2）中获取的区域连续N景Sentinel-1A标准格式数据进行D-InSAR处理，所述处理的方法是，对第1景与第2景进行D-InSAR处理，对第2景与第3景进行D-InSAR处理，对第3景与第4景进行D-InSAR处理......，依次类推，对第N-1景与第N景进行D-InSAR处理；

4）获取检测区域的积雪连续相位变化

在步骤3）区域合成孔径雷达数据D-InSAR处理的基础上，获取区域的积雪相位变化，具体方法是，得到第1景与第2景之间相位变化T_2-1、第2景与第3景之间相位变化T_3-2、第3景与第4景之间相位变化T_4-3、......、第N-1景与第N景之间相位变化T_N-N-1；

5）计算基于线性模型的检测区域积雪相位变化的积雪深度变化

在步骤4）获取区域积雪连续相位变化的基础上，利用线性模型计算区域积雪相位变化的积雪深度变化，公式如下：

(1)

式中，T为积雪相位变化，H为积雪深度变化（cm），计算得到区域第1景与第2景之间相位变化T_2-1对应的积雪深度变化H_2-1、第2景与第3景之间相位变化T_3-2对应的积雪深度变化H_3-2、第3景与第4景之间相位变化T_4-3对应的积雪深度变化H_4-3、......、第N-1景与第N景之间相位变化T_N-N-1对应的积雪深度变化H_N-N-1；

6）计算基于指数模型的检测区域积雪相位变化的积雪深度变化

在步骤4）获取区域积雪连续相位变化的基础上，利用指数模型计算区域积雪相位变化的积雪深度变化，公式如下：

(2)

式中，T为积雪相位变化，D为积雪深度变化（cm），计算得到区域第1景与第2景之间相位变化T_2-1对应的积雪深度变化D_2-1、第2景与第3景之间相位变化T_3-2对应的积雪深度变化D_3-2、第3景与第4景之间相位变化T_4-3对应的积雪深度变化D_4-3、......、第N-1景与第N景之间相位变化T_N-N-1对应的积雪深度变化D_N-N-1；

7）获取检测区域积雪相位变化的积雪深度变化

在步骤5）及步骤6）计算基于线性、指数模型的区域积雪相位变化的积雪深度变化的基础上，获取区域积雪相位变化的积雪深度变化；

具体方法是求算两者的平均值，得到第1景与第2景之间相位变化T_2-1对应的积雪深度变化（H_2-1+D_2-1）/2、第2景与第3景之间相位变化T_3-2对应的积雪深度变化（H_3-2+D_3-2）/2、第3景与第4景之间相位变化T_4-3对应的积雪深度变化（H_4-3+D_4-3）/2、......、第N-1景与第N景之间相位变化T_N-N-1对应的积雪深度变化（H_N-N-1+D_N-N-1）/2；

8）获取检测区域的积雪深度

在步骤7）获取检测区域积雪相位变化的积雪深度变化基础上，获取检测区域的积雪深度，公式如下：

(3)

式中，A为检测区域的积雪深度cm。