【ICESat-2处理】用ATL03提取地面点和冠层点（python）

【备注：这一篇真是耗费我太多时间了，我真是提了好久，所以要设置收费，请见谅】

我最近整理了一点ICESat-2的处理代码，我看现在网上的分享内容还挺少的，所以就想把自己整理的发出来吧，方便大家用。我是一个粗人，我不怎么会写代码，不喜欢只分享处理的关键函数，所以我的代码都是从数据读取到最后的可视化都有，直接拿去运行的。请读者们根据自己的需求挑选吧。

前言

这一篇主要是简单还原ATL08的处理过程，ATL08的处理文档是说用的迭代中值滤波法：atl08-v005-userguide_1_0.pdf (nsidc.org)；

原文翻译过来就是：“一条平滑的线将位于冠层顶部和地面之间。一旦这条线建立起来，算法就会迭代地将线上方的点标记为冠层。这个过程重复五次，以消除位于估计表面之上的点以及位于地面以下的噪声点”。

说实在我不懂原理，我是先去噪（看我之前的文：【ICESat-2处理】ATL03的去噪（python）-CSDN博客），然后迭代了几次寻找中值，然后把中值上的点标记为冠层，中值下的点标记为地面。反正我感觉我的方法也是挺糙的&