基于google earth engine的中等分辨率全国水质反演

我写博客的工作不像论文，假大空，我们直接上干货，之所以取一个这么大的名字，当然是我们能做到的。。。

不多说，我们对全国水体进行水质参数反演，不用MODIS，太粗，我们直接用哨兵，这样就可以直接做到大型水库或长河流观测了。当然了。

算法选择很重要，要选择什么算法呢？我这里主要参考这篇论文---《Hendrik V D W , Marcel W . Hue-Angle Product for Low to Medium Spatial Resolution Optical Satellite Sensors[J]. Remote Sensing, 2018, 10(2):180-.192》；

这是2018年发表的一篇水色遥感文章，发表在RS期刊上面，这里面重点讨论了基于Hu距反演水色的可行，尤其是在高分辨率上，并验证了可行性，论文也给了相关校正参数等。具体的算法细节实现如下：

1.水体提取，关于水体提取算法很多，我自己选择了2018年发表的一篇水体指数算法，参考文献-《A robust Multi-Band Water Index (MBWI) for automated extraction of surface water from Landsat 8 OLI imagery》

当然了，这个是用的Landat8传感器，但是我把它改成了哨兵2影像，结果也很好，尤其是山体阴影、冰雪等地区，表现确实比其他水体算法要好，大家可以参考下。

2.光谱角度转换，这个可以参见原文所示，我这里用了原文的60米分辨率参数，用了前5个波段。

3.FUI指数计算，这个就不多说了，但是GEE上面实现javascript代码有点复杂，这也是GEE写复杂算法的短板，它没有复杂的数据结构，确实很难编码。

4.透明度计算，水色级别等等，就不说了。。。

好了，废话不多说，我们看一下结果：

图 原始真彩色影像

图 Muwi水体指数

 图 FUI指数

我这里也在天津地区做了一个GUI demo示例，类似于webgis的界面功能，整体界面如下图所示（访问链接：https://wang749195.users.earthengine.app/view/watermonitor，特别注意，需要翻墙才能访问！）：

图 GUI 界面水体质量监测

我这里是选择2019年夏季天津入海口局部地区作为示例，其他地区大家可以自行测试，总的来说，效果还是不错的，但是也有一定的缺陷，

就是这个精度怎么验证呢？我个人的思路是，这个东西可以当做一个趋势分析产品，不能当做是一个业务化运行产品，用这个东西去分析水质参数，

意义还是非常大的，再结合GEE平台，几乎可以做到及时响应，比本地下载数据，处理，算法编写，要强的多。。。

最后，如果大家需要代码或者产品，请加我qq吧1044625113。。。备注：水质参数反演