ICESat-2/ATLAS06卫星数据处理2(附代码)-冰川物质平衡监测

本文主要是介绍ICESat-2/ATLAS06卫星数据处理2(附代码)-冰川物质平衡监测,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

    冰川是冰冻圈系统的重要组成部分,对气候和环境变化具有高度敏感性,是监测气候变化的良好指示器。首先冰川的变化意味着冰川的物质和形式的时空变化。而由于气候变化导致的冰川融化或积聚导致冰川的物质平衡发生变化,最终增加了冰川的面积或融化了冰川的末端。而冰川就像一个坚固的水库,冰川的退缩会导致某些山区缺水,冰川融化的水将形成冰川湖。由湖泊和冰川融化形成的大多数冰川湖是无法控制的堰塞湖,如果冰川湖破裂,下游可能会发生洪水。同样的冰川高程变化直接反映了冰川厚度的变化,并影响着气候、降水等一系列生态问题。因此监测冰川表面高程变化是非常必要的任务,成为许多研究者研究的热点问题。冰川表面高程变化反映了冰川物质平衡变化,是冰川物质平衡研究的重要内容。

 

冰川表面厚度变化的监测常用的方法如花杆测量、雪坑测量、多期DEM比较、摄影测量、冰雷达、航空/航天测高,包括雷达高度计和激光雷达。传统的方法是花杆测量,在冰川表面安插花杆进行检测,但是开展野外工作相当不容易,工作难度大,而且这种方法局限于易于观测的,范围较小,海拔较低,人可以涉足的冰川,而且精度很难保证,这就与测量的时间间隔,测量的方法等有着很大的关系。雪坑测量是通过在研究区域挖掘雪坑,钻取冰芯或雪芯,不同时间序列的冰芯或雪芯也是不同时期冰雪的沉降。在钻取冰芯或雪芯得同时测定不同深度的密度,就可以计算出某一时期内的积累量, 但由于可实施区域范围小,且作业困难,通常不是很好的监测方法。随着GPS精度的不断提高GPS技术逐渐应用到冰川的变化监测中,但是由于冰川环境恶劣,大范围开展GPS测量仍有一定的难度。冰雷达常用于监测极地冰盖厚度、内部构造以及冰下地貌,来重建过去的冰雪累积率与消融率。但冰雷达容易信号衰减形成雷达反射空白区域。卫星测高技术在获取地形特征和高程信息方面取得了较大的进步,为冰川变化监测提供了便利。利用激光测高卫星可以大面积进行观测,监测到冰川的变化。利用激光测高卫星数据,能更好的反映出高程变化的信息,在精度上也有一定的保证。随着激光测高系统的发展,一定会将冰川监测带上新的高度。

 

ICESat-2的传感器ATLAS更加先进,覆盖区域也更大,在冰川区域得到的脚点数据也更加密集。ICESat-2于2018年9月15日上午6:02(太平洋标准时间)在联合发射联盟的最终Delta II火箭上从范登堡空军基地的太空发射中心2号起飞。从那时起,在“冰桥行动”(Operation IceBridge)期间,为了填补ICESat-2发射前的数据缺口,研究人员一直在北极和南极地区飞行特别装备的P-3B和DC-8飞机。ICESat-2卫星上配有ATLAS(高级地形激光测高仪系统),其发射的光子需要3.3毫秒就能到达地球并返回。根据美国国家航空航天局的说法,ATLAS的测量误差可缩小到十亿分之一秒。利用物理定律,科学家们就可以计算出卫星到地球表面的距离,准确度精确到4毫米。ATL06号高程数据产品,它在沿轨和垂轨方向有计算的成对激光束的表面高,表1列举了ICESat-1卫星参数与ICESat-2卫星参数对比。

 

表1 ICESat-1卫星参数与ICESat-2卫星参数对比

Table 1 Comparisonof icesat-1 and icesat-2 satellite parameters

ICESat-1项目

参数

ICESat-2项目

参数

卫星标识

ICESAT-1

卫星标识      

ICESAT-2

运载火箭

Deltall2运载

运载火箭

Deltall2运载

发射地点

加利福尼亚州

发射地点

加利福尼亚州

设计寿命/a

3

设计寿命/a

3

轨道高度/Km

590

轨道高度/Km

500

轨道倾角/度

94

轨道倾角/度

92

轨道类型

非太阳同步

轨道类型

非太阳同步

轨道周期/d

183

轨道周期/d

91

开始工作时间

2003年1月

开始工作时间

2018年9月

波长/nm

1064;532

波长/nm

532

波束数量

1

波束数量

6

脉冲频率/HZ

40

脉冲频率/HZ

10K

光斑直径/m

60~70

光斑直径/m

17.5

脉冲能量/mJ

1064;75

532;32

脉冲能量/mJ

0.04~0.12

光斑间距/km

航向 0.17

旁向 2.5~15

光斑间距/km

航向 0.09

旁向 2.5~3.3

重复周期/d

183

重复周期/d

91

测高精度/m

0.15

测高精度/m

0.1

数据采集方式

模数转换

数据采集方式

光子计数

数据来源:https://nsidc.org/data/ICESat/(美国冰雪数据中心) 

 

ICESat卫星测高原理:

  激光高度计的测高原理是联合多个子系统进行定位,由GLAS传感器发射激光脉冲,通过穿过大气层的激光脉冲在地面产生的后向散射回波,分析激光脉冲回波的传播时间,因此可以计算出激光高度计与地面光斑之间的距离。联合星载GPS与地基卫星激光测距进行精密定轨,以确定传感器的空间位置,利用恒星相机与激光发射参考系统得到卫星的姿态信息和激光束的传播方向,联合这些观测值即可解算地面激光光斑的三维坐标。

 

 

ATLAS的工作模式如图1所示:我们可以看出ATLAS一共发射6束激光脉冲,分3组平行排列,每组之间地表距离约3km,组内2束激光脉冲间隔为90m,且为一强一弱,这样激光脉冲呈强弱间隔排列。由于可见绿光对光子的敏感性,激光波长不再像ICESat-1那样采用1064nm和532nm两个波长,而是使用532nm。

图1 多光束激光测量的示意图

Figure Schematic diagram of multi beamlaser measurement

数据处理代码实现:(详见公众号)

公众号二维码.jpg

 

结果展示:

图2 ICESat-2脚点轨迹以及时间

 

Figure2 trajectory and time of ICESat-2 footpoint

展望:

星载激光高度计因其搭载平台的特殊性,和地面及机载激光测距设备相比,具有许多优势。激光测高技术已有20多年的历史。他的技术已经相当成熟,但随着激光技术的发展和应用需求的扩大,激光测高技术仍有很大的发展空间。激光系统是高频的,随着窄束多光束扫帚的发展,单光子探测技术的发展也可以实现微脉冲距离的测量。未来,随着技术的发展,各种技术问题都将被破坏,激光高度计将有很好的发展和应用前景。未来需要利用时空分辨率更高的遥感影像和以及大量的地面控制点,结合激光高度卫星,以便在空间和时间上覆盖整个研究区域,来获取更精确的计算结果。

 


 

#ICESat-2/ATLAS06号陆地表面高程数据的处理与应用就先介绍到这里了,期待下次相见...

 


 

这篇关于ICESat-2/ATLAS06卫星数据处理2(附代码)-冰川物质平衡监测的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/526223

相关文章

LangChain转换链:让数据处理更精准

1. 转换链的概念 在开发AI Agent(智能体)时,我们经常需要对输入数据进行预处理,这样可以更好地利用LLM。LangChain提供了一个强大的工具——转换链(TransformChain),它可以帮我们轻松实现这一任务。 转换链(TransformChain)主要是将 给定的数据 按照某个函数进行转换,再将 转换后的结果 输出给LLM。 所以转换链的核心是:根据业务逻辑编写合适的转换函

uniapp接入微信小程序原生代码配置方案(优化版)

uniapp项目需要把微信小程序原生语法的功能代码嵌套过来,无需把原生代码转换为uniapp,可以配置拷贝的方式集成过来 1、拷贝代码包到src目录 2、vue.config.js中配置原生代码包直接拷贝到编译目录中 3、pages.json中配置分包目录,原生入口组件的路径 4、manifest.json中配置分包,使用原生组件 5、需要把原生代码包里的页面修改成组件的方

公共筛选组件(二次封装antd)支持代码提示

如果项目是基于antd组件库为基础搭建,可使用此公共筛选组件 使用到的库 npm i antdnpm i lodash-esnpm i @types/lodash-es -D /components/CommonSearch index.tsx import React from 'react';import { Button, Card, Form } from 'antd'

17.用300行代码手写初体验Spring V1.0版本

1.1.课程目标 1、了解看源码最有效的方式,先猜测后验证,不要一开始就去调试代码。 2、浓缩就是精华,用 300行最简洁的代码 提炼Spring的基本设计思想。 3、掌握Spring框架的基本脉络。 1.2.内容定位 1、 具有1年以上的SpringMVC使用经验。 2、 希望深入了解Spring源码的人群,对 Spring有一个整体的宏观感受。 3、 全程手写实现SpringM

代码随想录算法训练营:12/60

非科班学习算法day12 | LeetCode150:逆波兰表达式 ,Leetcode239: 滑动窗口最大值  目录 介绍 一、基础概念补充: 1.c++字符串转为数字 1. std::stoi, std::stol, std::stoll, std::stoul, std::stoull(最常用) 2. std::stringstream 3. std::atoi, std

记录AS混淆代码模板

开启混淆得先在build.gradle文件中把 minifyEnabled false改成true,以及shrinkResources true//去除无用的resource文件 这些是写在proguard-rules.pro文件内的 指定代码的压缩级别 -optimizationpasses 5 包明不混合大小写 -dontusemixedcaseclassnames 不去忽略非公共

基于Springboot + vue 的抗疫物质管理系统的设计与实现

目录 📚 前言 📑摘要 📑系统流程 📚 系统架构设计 📚 数据库设计 📚 系统功能的具体实现    💬 系统登录注册 系统登录 登录界面   用户添加  💬 抗疫列表展示模块     区域信息管理 添加物资详情 抗疫物资列表展示 抗疫物资申请 抗疫物资审核 ✒️ 源码实现 💖 源码获取 😁 联系方式 📚 前言 📑博客主页:

麻了!一觉醒来,代码全挂了。。

作为⼀名程序员,相信大家平时都有代码托管的需求。 相信有不少同学或者团队都习惯把自己的代码托管到GitHub平台上。 但是GitHub大家知道,经常在访问速度这方面并不是很快,有时候因为网络问题甚至根本连网站都打不开了,所以导致使用体验并不友好。 经常一觉醒来,居然发现我竟然看不到我自己上传的代码了。。 那在国内,除了GitHub,另外还有一个比较常用的Gitee平台也可以用于

剑指offer(C++)--平衡二叉树

题目 输入一棵二叉树,判断该二叉树是否是平衡二叉树。 class Solution {public:bool IsBalanced_Solution(TreeNode* pRoot) {if(pRoot==NULL)return true;int left_depth = getdepth(pRoot->left);int right_depth = getdepth(pRoot->rig

众所周知,配置即代码≠基础设置即代码

​前段时间翻到几条留言,问: “配置即代码和基础设施即代码一样吗?” “配置即代码是什么?怎么都是基础设施即代码?” 我们都是知道,DevOp的快速发展,让服务器管理与配置的时间大大减少,配置即代码和基础设施即代码作为DevOps的重要实践,在其中起到了关键性作用。 不少人将二者看作是一件事,配置即大代码是关于管理特定的应用程序配置设置本身,而基础设施即代码更关注的是部署支持应用程序环境所需的