【技术实现步骤摘要】
多年冻土活动层厚度的遥感估算方法、装置及可读存储介质
本专利技术涉及遥感
，具体而言，涉及一种多年冻土活动层厚度的遥感估算方法、装置及可读存储介质。
技术介绍
目前，多年冻土活动层厚度观测主要依赖于地面测量，主要测量手段有冻融管或钻孔观测，其空间代表性差，费时费力，成本高且时效性不强，且难以获取较大范围空间连续分布的活动层厚度信息。而在空间上估计活动层厚度多以斯蒂芬公式和库德里亚夫采夫经验公式为主，这些经验模型依赖于气象站的气温或土壤表面温度观测，由于多年冻土区气象站极其稀少，这严重限制了对多年冻土活动层厚度的监测。
技术实现思路
为了克服上述现有技术中的不足，本专利技术提供了一种多年冻土活动层厚度的遥感估算方法、装置及可读存储介质。为了实现上述目的，本专利技术实施例所提供的技术方案如下所示：第一方面，本专利技术实施例提供一种多年冻土活动层厚度的遥感估算方法，用于对多年冻土活动层的厚度进行估算，所述遥感估算方法包括：获取被估算区域的遥感数据和辅助数据，所述遥感数据包括：地表温度融化指...
【技术保护点】
1.一种多年冻土活动层厚度的遥感估算方法，其特征在于，所述遥感估算方法包括：/n获取被估算区域的遥感数据和辅助数据，所述遥感数据包括：地表温度融化指数DDT、年平均叶面积指数LAI和积雪日数SnowD，所述辅助数据包括相关土壤数据，所述相关土壤数据的类别包括：土壤沙粒含量Sand、土壤粘粒含量Clay、土壤粉粒含量Silt、土壤有机碳含量Org和土壤密度Bld；/n根据预设的多年冻土活动层厚度的计算公式、所述遥感数据及所述辅助数据估算所述被估算区域的多年冻土活动层厚度；其中，所述多年冻土活动层厚度的计算公式为：/n
【技术特征摘要】
1.一种多年冻土活动层厚度的遥感估算方法，其特征在于，所述遥感估算方法包括：
获取被估算区域的遥感数据和辅助数据，所述遥感数据包括：地表温度融化指数DDT、年平均叶面积指数LAI和积雪日数SnowD，所述辅助数据包括相关土壤数据，所述相关土壤数据的类别包括：土壤沙粒含量Sand、土壤粘粒含量Clay、土壤粉粒含量Silt、土壤有机碳含量Org和土壤密度Bld；
根据预设的多年冻土活动层厚度的计算公式、所述遥感数据及所述辅助数据估算所述被估算区域的多年冻土活动层厚度；其中，所述多年冻土活动层厚度的计算公式为：
其中，a、b、c、d、e、f、g、h和T为预先获取到的估算参数；ALT为所述被估算区域的多年冻土活动层厚度。
2.根据权利要求1所述的估算方法，其特征在于，在根据预设的多年冻土活动层厚度的计算公式、所述遥感数据及所述辅助数据估算所述被估算区域的多年冻土活动层厚度之前，所述估算方法还包括：
预先获取所述被估算区域的所述遥感数据及所述辅助数据的多组数据以及所述多组数据对应的所述被估算区域的多年冻土活动层厚度、以及所述计算公式；
根据所述多组数据、所述多组数据对应的所述多年冻土活动层厚度及所述多年冻土活动层厚度的计算公式，通过参数估计方法得到所述估算参数。
3.根据权利要求1所述的估算方法，其特征在于，获取所述地表温度融化指数DDT的步骤包括：
对获取的一日四次的白天和夜晚地表温度，根据与所述白天和夜晚地表温度一起获取的相应质量控制符QC，选择质量最好的进行白天和夜晚地表温度的配对，对配对成功的像元直接进行平均，得到该像元的日平均地表温度；
提取配对成功像元的日地表温度变幅，利用克里金方法插值到每个像元，得到全年每日全部像元的地表温度变幅D；
对于没有形成配对的像元，根据当天全部像元的所述地表温度变幅D，计算该天没有形成配对像元的日平均地表温度；
根据前三步获取的日平均地表温度，利用时空增强方法估计完全无值像元的日平均地表温度值，得到时空完整的日平均地表温度；
根据时空完整的每日平均地表温度，确定一年中的大于0℃的每日平均地表温度的代数和为地表温度融化指数DDT。
4.根据权利要求1所述的估算方法，其特征在于，获取所述年平均叶面积指数LAI的步骤包括：
判断是否存在时间上缺失的被估算区域的叶面积指数；
根据时间序列重建方法对所述缺失的叶面积指数进行...
【专利技术属性】
技术研发人员：冉有华，
申请(专利权)人：中国科学院寒区旱区环境与工程研究所，
类型：发明
国别省市：甘肃;62
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。