多年冻土活动层厚度估算方法

技术特征：
1.一种多年冻土活动层厚度估算方法，其特征在于，包括：获取估算区域的多个不同轨道的sar影像数据集；分别对多个所述sar影像数据集进行时序insar处理，得到每个所述sar影像数据集的los向地表时序形变；根据每个所述sar影像数据集的los向地表时序形变，获取所述估算区域的地表真实二维形变结果；获取所述估算区域的遥感数据和土壤数据，并根据所述遥感数据和所述土壤数据对所述估算区域的未冻水含量进行预报；根据所述地表真实二维形变结果和所述未冻水含量，对所述估算区域的活动层厚度进行估算。2.根据权利要求1所述的方法，其特征在于，所述根据每个所述sar影像数据集的los向地表时序形变，获取所述估算区域的地表真实二维形变结果，包括：分别将每个所述sar影像数据集的los向地表时序形变分解为垂直、东西和南北三个方向的形变分量，得到每个轨道的观测方程；利用多个轨道的观测方程获取所述估算区域的地表真实二维形变结果。3.根据权利要求2所述的方法，其特征在于，第i个轨道的观测方程为：其中，表示第i个轨道的sar影像数据集的los向地表时序形变，d
u
表示在垂直方向的形变分量，d
e
表示在东西方向的形变分量，d
n
表示在南北方向的形变分量，θ
i
表示第i个轨道的sar影像数据集的入射角，α
i
表示第i个轨道的sar影像数据集的方位角，i＝1,
…
,j，j表示轨道的总数。4.根据权利要求3所述的方法，其特征在于，所述利用多个轨道的观测方程获取所述估算区域的地表真实二维形变结果，包括：利用最小二乘法对进行求解，得到所述估算区域的地表真实二维形变结果；所述地表真实二维形变结果包括d
u
和d
e
。5.根据权利要求1所述的方法，其特征在于，所述遥感数据包括多层土壤的温度数据，所述土壤数据包括每层所述土壤的土质类型以及每种土质类型对应的未冻水含量常数；所述根据所述遥感数据和所述土壤数据对所述估算区域的未冻水含量进行预报，包括：利用公式w
u
＝aθ-b
计算所述估算区域的未冻水含量；其中，w
u
表示所述估算区域的未冻水含量，θ表示所述温度数据的绝对值，a和b均表示所
述未冻水含量常数。6.根据权利要求4所述的方法，其特征在于，在所述根据每个所述sar影像数据集的los向地表时序形变，获取所述估算区域的地表真实二维形变结果之后，所述方法还包括：利用傅里叶拟合方法对所述地表真实二维形变结果中的d
u
进行拟合；从拟合结果中提取融化季时段对应的垂直方向的形变分量7.根据权利要求6所述的方法，其特征在于，所述根据所述地表真实二维形变结果和所述未冻水含量，对所述估算区域的活动层厚度进行估算，包括：通过公式计算所述估算区域从活动层融化开始t0至结束时刻t
k
的活动层厚度h；其中，δh
r
表示所述估算区域在t
r-1
至t
r
时刻融化的活动层厚度，时刻融化的活动层厚度，时刻融化的活动层厚度，表示在t
r
时刻的值，表示在t
r-1
时刻的值，w
ice
表示所述估算区域在融化季时段的t
r-1
至t
r
时刻参与水冰相变过程的冰含量，ρ
water
表示水的密度，ρ
ice
表示冰的密度，w
r
表示所述估算区域在t
r
时刻的土壤水分含量，表示所述估算区域在t
r-1
时刻的未冻水含量。

技术总结
本申请适用于多年冻土监测技术领域，提供了一种多年冻土活动层厚度估算方法，其中该方法包括：获取估算区域的多个不同轨道的SAR影像数据集；分别对多个SAR影像数据集进行时序InSAR处理，得到每个SAR影像数据集的LOS向地表时序形变；根据每个SAR影像数据集的LOS向地表时序形变，获取估算区域的地表真实二维形变结果；获取估算区域的遥感数据和土壤数据，并根据遥感数据和土壤数据对估算区域的未冻水含量进行预报；根据地表真实二维形变结果和未冻水含量，对估算区域的活动层厚度进行估算。本申请能提高活动层厚度估算的精度和可靠性。本申请能提高活动层厚度估算的精度和可靠性。本申请能提高活动层厚度估算的精度和可靠性。


技术研发人员：王琪洁 李佳晨 张亚 杨沙
受保护的技术使用者：中南大学
技术研发日：2022.11.08
技术公布日：2022/12/30