今天看了一篇報道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