今天看了一篇报道http://news.xhby.net/js/jj/201908/t20190801_6283881.shtml(直击60年一遇气象干旱:洪泽湖缩小近一半,鱼蟹受灾严重!),既然玩GEE那就要玩出点花样来,我们看看洪泽湖那边究竟变成什么样子了~~。
第一部分统计结果
(1)洪泽湖的矢量边界,这个直接使用的现有的矢量边界(自己没去重新做?),然后统计的面积是1700多平方千米。(百度上说是1700多平方千米,维基百科说是1500多平方千米,都有道理)
文章图片
(2)统计去年2018年洪泽湖的水面面积是1400多平方千米,这个误差主要是矢量边界会把湖边界中的湿地也算进去。
文章图片
(3)今年2019年7月份的数据如下,水面面积只有1000多平方千米,和新闻报道中的数据(900多平方千米)差不多,而且主要消失的水体集中在成子湖(最上面这部分)。
文章图片
相关代码截图
文章图片
具体代码就不讲了,这里说一个我在代码中常用到的方法,生成年度合成数据,具体来讲就是将一个时间序列数据按照时间年来合并,比如我这个取的是中值。使用循环来做只是其中一种方式,如果大家熟悉join,还可以使用join来做。
function getYearCol(sDate, eDate, lxCol, region) {var yearList = ee.List.sequence(ee.Date(sDate).get("year"), ee.Number(ee.Date(eDate).get("year")).subtract(1));
var yearImgList = yearList.map(function(year) {year = ee.Number(year);
var _sdate = ee.Date.fromYMD(year, 1, 1);
var _edate = ee.Date.fromYMD(year.add(1), 1, 1);
var tempCol = lxCol.filterDate(_sdate, _edate);
var img = tempCol.median().clip(region);
img = img.set("year", year);
img = img.set("system:index", ee.String(year.toInt()));
return img;
});
var yearImgCol = ee.ImageCollection.fromImages(yearImgList);
print("yearImgCol",yearImgCol);
return yearImgCol;
}
第二部分原因分析
事先声明,这些都只是我从现有的遥感数据以及气象数据分析,仅代表个人非常不成熟的观点。造成湖面消失的原因无非是温度过高水体蒸发速度过快;降水量比较小,水源无法及时补充;再有就是人们用水量突增等。我这里只分析一下温度以及降水(这两个数据最好找~~)。
(1)温度,我这里对比了去年的状况,6、7月份比去年温度稍高一些,由于这是计算的平均值,所以可以想到整体区域温度应该会比往常要高。而且我们从其他途径也可以了解到今年全球气候反常,问题已经达到历史之最。
文章图片
看看今年6、7月份气象影像,确实是一片红(热!)
文章图片
文章图片
(2)降水量
文章图片
大致可以看到相比去年的每月的累积降水量,今年确实少一些。如果图看的不是非常清晰,我们统计一下具体数据,如下图,可以看到累计降水量少了将近200毫米。
文章图片
相关代码截图
文章图片
第三部分结论分析
高温少雨,同时造成周围用水量增加,直接很大程度上影响了洪泽湖的蓄水量,如果不采取更多措施,那么只会加剧洪泽湖目前的窘状。
文章图片
大家如果有问题需要交流或者有项目需要合作,可以微信联系我,加微信好友请留言加上“GEE”。
知乎专栏:https://zhuanlan.zhihu.com/c_123993183
CSDN:https://blog.csdn.net/shi_weihappy
【网络前沿知识|GEE学习笔记 七十八(干涸的洪泽湖)】
推荐阅读
- 网络前沿知识|GEE学习笔记 七十六(【GEE之Python版教程十】字典)
- GEE|GEE学习笔记 五十九(GEE中mosiac、qualityMosiac、max区别)
- python|GEE学习笔记 七十五(【GEE之Python版教程九】数值)
- python|GEE学习笔记 七十(【GEE之Python版教程四】Python基础编程二)
- GEE|GEE学习笔记 六十二(保密是头等大事)
- GEE|GEE学习笔记 五十八(设置数值数据保留N位)
- GEE|GEE学习笔记 五十一(Fusion Table将在2019年12月3日关闭)
- GEE|GEE学习笔记 五十六(GEE中如何把文件导出到Google Drive的子目录)
- GEE|GEE学习笔记 五十三(GEE导出的影像波段名称哪去了())