GEE学习笔记 七十八：干涸的洪泽湖

今天看了一篇报道直击60年一遇气象干旱：洪泽湖缩小近一半，鱼蟹受灾严重！_新华报业网（直击60年一遇气象干旱：洪泽湖缩小近一半，鱼蟹受灾严重！），既然玩GEE那就要玩出点花样来，我们看看洪泽湖那边究竟变成什么样子了～～。

第一部分统计结果

（1）洪泽湖的矢量边界，这个直接使用的现有的矢量边界（自己没去重新做?），然后统计的面积是1700多平方千米。（百度上说是1700多平方千米，维基百科说是1500多平方千米，都有道理）

（2）统计去年2018年洪泽湖的水面面积是1400多平方千米，这个误差主要是矢量边界会把湖边界中的湿地也算进去。

（3）今年2019年7月份的数据如下，水面面积只有1000多平方千米，和新闻报道中的数据（900多平方千米）差不多，而且主要消失的水体集中在成子湖（最上面这部分）。

相关代码截图

    具体代码就不讲了，这里说一个我在代码中常用到的方法，生成年度合成数据，具体来讲就是将一个时间序列数据按照时间年来合并，比如我这个取的是中值。使用循环来做只是其中一种方式，如果大家熟悉join，还可以使用join来做。

1. function getYearCol(sDate, eDate, lxCol, region) {

3.   var yearList = ee.List.sequence(ee.Date(sDate).get("year"), ee.Number(ee.Date(eDate).get("year")).subtract(1));

5.   var yearImgList = yearList.map(function(year) {

7.     year = ee.Number(year);

9.     var _sdate = ee.Date.fromYMD(year, 1, 1);

11.     var _edate = ee.Date.fromYMD(year.add(1), 1, 1);

13.     

15.     var tempCol = lxCol.filterDate(_sdate, _edate);

17.     var img = tempCol.median().clip(region);

19.     img = img.set("year", year);

21.     img = img.set("system:index", ee.String(year.toInt()));

23.     return img;

25.   });

27.   

29.   var yearImgCol = ee.ImageCollection.fromImages(yearImgList);

31.   print("yearImgCol",yearImgCol);

33.   return yearImgCol;

35. }

第二部分原因分析

    事先声明，这些都只是我从现有的遥感数据以及气象数据分析，仅代表个人非常不成熟的观点。造成湖面消失的原因无非是温度过高水体蒸发速度过快；降水量比较小，水源无法及时补充；再有就是人们用水量突增等。我这里只分析一下温度以及降水（这两个数据最好找～～）。

（1）温度，我这里对比了去年的状况，6、7月份比去年温度稍高一些，由于这是计算的平均值，所以可以想到整体区域温度应该会比往常要高。而且我们从其他途径也可以了解到今年全球气候反常，问题已经达到历史之最。

看看今年6、7月份气象影像，确实是一片红（热！）

（2）降水量

大致可以看到相比去年的每月的累积降水量，今年确实少一些。如果图看的不是非常清晰，我们统计一下具体数据，如下图，可以看到累计降水量少了将近200毫米。

相关代码截图

第三部分结论分析

    高温少雨，同时造成周围用水量增加，直接很大程度上影响了洪泽湖的蓄水量，如果不采取更多措施，那么只会加剧洪泽湖目前的窘状。

来源请引用：地理遥感生态网科学数据注册与出版系统。