一、背景

这是一篇做冰湖溃决的文章，作者主要使用了random forest来检测喜马拉雅地区的冰湖溃决现象，这项成果发表在了Remote Sensing of Environment上。

文献连接：https://doi.org/10.1016/j.rse.2017.12.025

文献引用：Georg Veha, Oliver Korupa, Sigrid Roessnerb, Ariane Walz. Detecting Himalayan glacial lake outburst floods from Landsat time series. Remote Sensing of Environment, vol. 207, pp. 87-94, 2018.

二、具体工作

首先，作者在摘要里提到，使用的数据是Landsat从1988年到2016年，用于寻找由于GLOF导致的冰湖面积变化的现象。一共使用了2491张图，用随机森林来生成模糊地表覆盖图，总体的精度为91%。使用基于似然的变化点检测方法，来估计GLOF发生的时间，总计找到了10/11文件记录到的冰湖溃决事件。总体来说，自动冰湖愧疚检测算法仍然是一个具有挑战的任务，需要解决的问题不局限于图像配准（co-registration），误分类，图像噪声等。

1. 背景

文章的开始，作者介绍了很多冰湖溃决监测的意义，并总结出了近期一些文献关于冰湖溃决事件的整理，具体如下：

这里作者能中RSE，很重要的一个原因是：

To our knowledge, no study has systematically explored the Landsat archive for retrospective GLOF detection in the HKH, although it offers a largely continuous, nearly 30-year time series with regional coverage every 16 days.

没有研究是基于Landsat来系统性的检测冰湖事件的。

其原因是两方面的：一方面是数据的问题，高原地区常常受到季风的影响，导致云层较为密集，且山体阴影、图像噪声对冰湖也有影响；另一方面，冰湖制图的方法多是半自动化的，而且需要大量的后处理工作。

因此本研究考虑采用了一种基于机器学习的方法，即随机森林来做冰湖溃决的检测。简单来说，这篇文章的工作主要是：

(1) a Random-Forest based land cover classification.

(2) a novel, likelihood based change-point algorithm to approximate the time stamp of GLOFs.

在过去70年中，一共文献记录了36起GLOF，本文找到了11起，这11起如上表1所示，他们的实际地理位置则如下图所示：

2. 数据

下面是作者所使用到的数据以及来源：

具体来说：

①Landsat数据。下载地址可从USGS上下载到，这11起GLOF涉及到19景影像。作者收集了9月-11月之间，云量小于60%的，1988-2016年间的所有影像，一共下载了2491景，并进行了TOA校正，所有的数据情况如下图所示：

其中，喜马拉雅地区中部拍摄最为密集的地区，每个像素达到429景影像。

②CFMask数据。即用C波段来检测云层，去除Landsat影像中的云层影响。下载的话还是在USGS上下载。它是将一景影像分为云，阴影，冰雪，水，空白。

③DEM数据。包括SRTM DEM和ALOS world 3D DEM，这两个数据分别在USGS上和JAXA上下载，之后再用9*9的高斯滤波对其平滑。

④冰川数据。冰川数据使用的ICIMOD和RGI，ICIMOD确实中国地区的冰川，因此使用RGI进行补齐。

3. 方法

整个流程结构如下图所示：

整体上来说分为三部分：

①数据预处理。

试验区跨了6个UTM投影带，因此需要先做重投影。然后将试验区分成30km*30km的小块进行处理。同时在冰川周围3km处做了缓冲区，并去除掉了非常陡峭的区域。

②随机森林分类。

首先定义了6种与冰湖溃决有关的土地覆盖类型。即为水，沉积物，冰雪，云，阴影，土地（不属于其他5类的全部属于这一类）。然后这11景包含GLOF的影像作为训练数据，其他的57景影像用作测试数据，用于检测未记录的GLOF。之后即使用随机森林对其进行分类，并将最终分类结果的沉积物和土地合并为一类。这里需要注意的是，作者的训练数据选择的是样本点。随机森林设置了1000棵数。评估随机森林的分类，采用了对数损失函数：