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

相关文章

C#实现千万数据秒级导入的代码

《C#实现千万数据秒级导入的代码》在实际开发中excel导入很常见,现代社会中很容易遇到大数据处理业务,所以本文我就给大家分享一下千万数据秒级导入怎么实现,文中有详细的代码示例供大家参考,需要的朋友可... 目录前言一、数据存储二、处理逻辑优化前代码处理逻辑优化后的代码总结前言在实际开发中excel导入很

SpringBoot+RustFS 实现文件切片极速上传的实例代码

《SpringBoot+RustFS实现文件切片极速上传的实例代码》本文介绍利用SpringBoot和RustFS构建高性能文件切片上传系统,实现大文件秒传、断点续传和分片上传等功能,具有一定的参考... 目录一、为什么选择 RustFS + SpringBoot?二、环境准备与部署2.1 安装 RustF

Python实现Excel批量样式修改器(附完整代码)

《Python实现Excel批量样式修改器(附完整代码)》这篇文章主要为大家详细介绍了如何使用Python实现一个Excel批量样式修改器,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一... 目录前言功能特性核心功能界面特性系统要求安装说明使用指南基本操作流程高级功能技术实现核心技术栈关键函

Redis实现高效内存管理的示例代码

《Redis实现高效内存管理的示例代码》Redis内存管理是其核心功能之一,为了高效地利用内存,Redis采用了多种技术和策略,如优化的数据结构、内存分配策略、内存回收、数据压缩等,下面就来详细的介绍... 目录1. 内存分配策略jemalloc 的使用2. 数据压缩和编码ziplist示例代码3. 优化的

Python 基于http.server模块实现简单http服务的代码举例

《Python基于http.server模块实现简单http服务的代码举例》Pythonhttp.server模块通过继承BaseHTTPRequestHandler处理HTTP请求,使用Threa... 目录测试环境代码实现相关介绍模块简介类及相关函数简介参考链接测试环境win11专业版python

Python从Word文档中提取图片并生成PPT的操作代码

《Python从Word文档中提取图片并生成PPT的操作代码》在日常办公场景中,我们经常需要从Word文档中提取图片,并将这些图片整理到PowerPoint幻灯片中,手动完成这一任务既耗时又容易出错,... 目录引言背景与需求解决方案概述代码解析代码核心逻辑说明总结引言在日常办公场景中,我们经常需要从 W

使用Spring Cache本地缓存示例代码

《使用SpringCache本地缓存示例代码》缓存是提高应用程序性能的重要手段,通过将频繁访问的数据存储在内存中,可以减少数据库访问次数,从而加速数据读取,:本文主要介绍使用SpringCac... 目录一、Spring Cache简介核心特点:二、基础配置1. 添加依赖2. 启用缓存3. 缓存配置方案方案

MySQL的配置文件详解及实例代码

《MySQL的配置文件详解及实例代码》MySQL的配置文件是服务器运行的重要组成部分,用于设置服务器操作的各种参数,下面:本文主要介绍MySQL配置文件的相关资料,文中通过代码介绍的非常详细,需要... 目录前言一、配置文件结构1.[mysqld]2.[client]3.[mysql]4.[mysqldum

Python多线程实现大文件快速下载的代码实现

《Python多线程实现大文件快速下载的代码实现》在互联网时代,文件下载是日常操作之一,尤其是大文件,然而,网络条件不稳定或带宽有限时,下载速度会变得很慢,本文将介绍如何使用Python实现多线程下载... 目录引言一、多线程下载原理二、python实现多线程下载代码说明:三、实战案例四、注意事项五、总结引

IDEA与MyEclipse代码量统计方式

《IDEA与MyEclipse代码量统计方式》文章介绍在项目中不安装第三方工具统计代码行数的方法,分别说明MyEclipse通过正则搜索(排除空行和注释)及IDEA使用Statistic插件或调整搜索... 目录项目场景MyEclipse代码量统计IDEA代码量统计总结项目场景在项目中,有时候我们需要统计