ICESat-2 ATL08链接到ATL03光子

ICESat-2 ATL08链接到ATL03光子

ICESat-2作为目前最先进的机载激光雷达系统，备受国内科研民工关注，作为一个开源数据，是广大童鞋赖以毕业的重要资源。但ICESat-2文档冗长，劝退了大量初次接触该数据的童鞋。同时，由于ICESat-2应用时间尚短，国内外相关的数据解析代码也寥寥无几。本着造福广大科研同胞，避免重复踩坑的宗旨，我决定把目前掌握的数据解析方法分享出来。

ICESat-2数据介绍

目前大家用的比较多的其实也就ATL03和ATL08两种数据。其中ATL03主要记录了原始的光子数据，主要包括所有光子的经纬度坐标。而ATL08为应用层数据，包括分类后的信号光子信息，提取的分段地面高程、树高等信息。具体包括的内容大家最好还是下个数据看一下，本文重点不在此，故不再过多赘述。

ICESat-2 ATL03和ATL08数据的对应关系

ATL08中的信号光子是没有记录光子坐标的，仅记录了分类属性。那么想要获取这些属性对应的光子坐标，就需要将ATL08链接到ATL03数据上。但这其中的链接逻辑是比较复杂的，我也是在痛苦摸索几天之后才搞明白。
此处仅介绍如何链接分类属性到光子坐标：
首先需要搞清楚的是，ATL08中的分类属性仅包括信号光子的属性，而不包括太阳背景噪声等造成的噪声光子，也就是做过去噪之后的数据。
那么此处我们先讲要用到的ATL08和ATL03中的字段，
ATL08： ph_segment_id (光子对应的20m分段的id)、d_flag(光子分类属性)、classed_pc_indx(光子在当前段中的索引)
ATL03：h_ph(光子高程)、lat_ph(光子纬度)、lon_ph(光子经度)、segment_id(20m分段的段号)、segment_ph_cont(每个分段中的光子数)、ph_index_beg(每个分段中第一个光子在所有光子中的索引)

看到用到的数据其实大概就已经能猜到两个文件的对应关系了，它是通过20m分段来进行对应的。那么我们要实现每个分段中光子的对应就需要按以下3步来进行：

1. 找到ATL08中数据在ATL03中的分段。
2. 按照ATL08分段中每个光子的索引对应查找其在ATL03所有光子中的索引。
3. 按照查找到的ATL03索引找到对应的经纬度和高程坐标，然后将其与ATL08中的属性组合起来。
   

根据ATL08数据找到对应的光子坐标数据

具体而言：
① 拿到ATL08一个光子。
② 用该光子的20m分段id减去ATL03数据中第一个20m分段id，得到这个光子对应于所有20m分段的索引idex。
③ 根据该index得到这个分段中第一个光子在所有光子数据中的索引 idx_all。
④ 根据idx_all加上ATL08光子在当前段内是第几个光子的索引，得到该光子在所有光子中的索引 idx。
⑤ 根据idx找到对应的经纬度坐标、高程等信息。
⑥ 如果需要的话将经纬度信息投影到投影坐标系。

以上即可得到ATL08分类结果所对应的光子数据：