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

相关文章

Python实例题之pygame开发打飞机游戏实例代码

《Python实例题之pygame开发打飞机游戏实例代码》对于python的学习者,能够写出一个飞机大战的程序代码,是不是感觉到非常的开心,:本文主要介绍Python实例题之pygame开发打飞机... 目录题目pygame-aircraft-game使用 Pygame 开发的打飞机游戏脚本代码解释初始化部

Java中Map.Entry()含义及方法使用代码

《Java中Map.Entry()含义及方法使用代码》:本文主要介绍Java中Map.Entry()含义及方法使用的相关资料,Map.Entry是Java中Map的静态内部接口,用于表示键值对,其... 目录前言 Map.Entry作用核心方法常见使用场景1. 遍历 Map 的所有键值对2. 直接修改 Ma

从基础到进阶详解Pandas时间数据处理指南

《从基础到进阶详解Pandas时间数据处理指南》Pandas构建了完整的时间数据处理生态,核心由四个基础类构成,Timestamp,DatetimeIndex,Period和Timedelta,下面我... 目录1. 时间数据类型与基础操作1.1 核心时间对象体系1.2 时间数据生成技巧2. 时间索引与数据

深入解析 Java Future 类及代码示例

《深入解析JavaFuture类及代码示例》JavaFuture是java.util.concurrent包中用于表示异步计算结果的核心接口,下面给大家介绍JavaFuture类及实例代码,感兴... 目录一、Future 类概述二、核心工作机制代码示例执行流程2. 状态机模型3. 核心方法解析行为总结:三

python获取cmd环境变量值的实现代码

《python获取cmd环境变量值的实现代码》:本文主要介绍在Python中获取命令行(cmd)环境变量的值,可以使用标准库中的os模块,需要的朋友可以参考下... 前言全局说明在执行py过程中,总要使用到系统环境变量一、说明1.1 环境:Windows 11 家庭版 24H2 26100.4061

pandas实现数据concat拼接的示例代码

《pandas实现数据concat拼接的示例代码》pandas.concat用于合并DataFrame或Series,本文主要介绍了pandas实现数据concat拼接的示例代码,具有一定的参考价值,... 目录语法示例:使用pandas.concat合并数据默认的concat:参数axis=0,join=

C#代码实现解析WTGPS和BD数据

《C#代码实现解析WTGPS和BD数据》在现代的导航与定位应用中,准确解析GPS和北斗(BD)等卫星定位数据至关重要,本文将使用C#语言实现解析WTGPS和BD数据,需要的可以了解下... 目录一、代码结构概览1. 核心解析方法2. 位置信息解析3. 经纬度转换方法4. 日期和时间戳解析5. 辅助方法二、L

Python使用Code2flow将代码转化为流程图的操作教程

《Python使用Code2flow将代码转化为流程图的操作教程》Code2flow是一款开源工具,能够将代码自动转换为流程图,该工具对于代码审查、调试和理解大型代码库非常有用,在这篇博客中,我们将深... 目录引言1nVflRA、为什么选择 Code2flow?2、安装 Code2flow3、基本功能演示

IIS 7.0 及更高版本中的 FTP 状态代码

《IIS7.0及更高版本中的FTP状态代码》本文介绍IIS7.0中的FTP状态代码,方便大家在使用iis中发现ftp的问题... 简介尝试使用 FTP 访问运行 Internet Information Services (IIS) 7.0 或更高版本的服务器上的内容时,IIS 将返回指示响应状态的数字代

MySQL 添加索引5种方式示例详解(实用sql代码)

《MySQL添加索引5种方式示例详解(实用sql代码)》在MySQL数据库中添加索引可以帮助提高查询性能,尤其是在数据量大的表中,下面给大家分享MySQL添加索引5种方式示例详解(实用sql代码),... 在mysql数据库中添加索引可以帮助提高查询性能,尤其是在数据量大的表中。索引可以在创建表时定义,也可