ICESat-2点云数据读取与简单处理（ATL03和ATL08）

闪光蛇纹岩

已于 2024-11-10 16:20:37 修改

阅读量3.6k

点赞数 40

文章标签： python matlab

于 2024-11-07 00:45:13 首次发布

本文链接：https://blog.csdn.net/m0_60667301/article/details/143448669

版权

一、数据下载

1、数据介绍

ICESat-2 是 NASA 系列任务中的第三项，ICESat-2通过接收发射的激光光子来获得地表信息（高程），理想情况下，其高程精度可以达到亚米级（0.1m），是研究冰川高程、植被覆盖、地表坡度的重要数据。其中，ATL03是L2A级别数据，ATL08是L3A级别数据。（关于详细的数据介绍可以访问官网：ICESat-2 (nsidc.org)）

2、数据下载

官网（Data | National Snow and Ice Data Center (nsidc.org)）提供了好几种下载方式，在这里，我演示从NASA官网（Earthdata Search | Earthdata Search (nasa.gov)）下载。

首先，①进入网站先登陆账号（右上角Login）；②选择目标数据的时间和空间范围（在这里我筛选从2019-01-01到2019-12-31的河北保定区域的数据）；③搜索栏检索“ATL03”；④可以看到右侧出现了目标数据集，点击该数据集。

可以将鼠标放在感兴趣的数据上，右侧就会出现绿色高光显示该数据的空间位置，也可以观察数据的预览图查看数据的完整性（有时候，虽然数据空间上与研究区重合，但在研究区范围内会因为地形、大气等因素存在数据缺失的情况，表现为预览图上数据的不连续），点击下载数据（ATL08的数据下载同理，只是检索的时候将“ATL03”替换为“ATL08”）。

二、数据读取

ICESat-2的点云数据是存储在HDF5文件中的，Hierarchical Data Format(HDF)是一种针对大量数据进行组织和存储的文件格式，HDF5文件可以组织成类似于文件系统的层次结构，所以在读取数据前需要找到正确的“路径”。

1、利用网站查看文件结构

myHDF5 是欧洲同步辐射光源ESRF（European Synchrotron Radiation Facility）开发和维护的在线HDF5文件查看服务，其界面如下图所示，选择下载的ICESat-2文件上传。打开文件可以看到其中有gt1l、gt1r、gt2l、gt2r、gt3l、gt3r（后面用 gtx 指代），分别对应了ICESat-2 ATLAS探测的6个光子轨道，点击 gtxl\heights\h_ph 就可以看到高程数据，也可以点击 gtxl\geophys_corr\dem_h

关于 h_ph 和 dem_h 两种高程的区别，可以查看官方文档（Documents），里面对每个数据集进行了详细的说明。

2、利用MATLAB读取点云数据

利用MATLAB读取三维点云数据，即包含经纬度坐标的 h_ph 高程信息的光子，其中经度存储在 gtx\heights\lon_ph，纬度存储在 gtx\heights\lat_ph，高程信息则在 gtx\heights\h_ph（gtx为gt1l、gt1r、gt2l、gt2r、gt3l、gt3r），在这里用MATLAB读取轨道gt1l的点云数据，代码如下（其它轨道的类似）：

最低0.47元/天解锁文章