GEE6：获取每年水体数据

1.JRC Monthly Water History, v1.2 [deprecated]

数据描述：

  该数据集包含1984年至2019年地表水的位置和时间分布地图，并提供了关于这些水面的范围和变化的统计数据。这些数据使用1984年3月16日至2019年12月31日期间从陆地卫星5号、7号和8号获得的4185439个场景生成的。使用专家系统将每个像素分别分类为水/非水，并将结果整理成整个时间段和两个时代(1984-1999,2000-2019)的月度历史，以进行变化检测。“月历史” 集合包含逐月的水检测的整个历史。该图集包含430张图片，从1984年3月到2019年12月，每个月一张。

2. GEE code

//获取哪一年的数据
for(var i=0;i<9;i++){
  if (i<20){  var year='200'+i;}
  if (i>20||i==20){  var year='20'+i;}
  var yearn = parseInt(JSON.parse(year));
  get_yearly_water(yearn);
}

function get_yearly_water(year) {
  
    //设置需要提取的区域，由于是上传的shp文件，需要转为geometry的格式
    var roi = table.geometry();

    //设置需要提取的年份
    var startDate = ee.Date.fromYMD(year, 1, 1);
    var endDate = ee.Date.fromYMD(year, 12, 31);
    
    //筛选JRC水体数据
    var myjrc = jrc.filterBounds(roi).filterDate(startDate, endDate);
    
    //在每个月份的影像中添加一个obs属性的波段，一个像素如果有数据，则为1，没有数据则为0
    myjrc = myjrc.map(function(img){
      var obs = img.gt(0);
      return img.addBands(obs.rename('obs').set('system:time_start', img.get('system:time_start')));
    });

    //在每个月份的影像中添加一个onlywater属性的波段，一个像素如果有水则为1，没有水则为0
    myjrc = myjrc.map(function(img){
      var water = img.select('water').eq(2);
      return img.addBands(water.rename('onlywater').set('system:time_start', img.get('system:time_start')));
    });
    
    //计算每个像素点在一年12景影像中， 有数据的次数
    var totalObs = ee.ImageCollection(myjrc.select('obs')).sum().toFloat();
    
    //计算每个像素点在一年12景影像中， 有水的次数
    var totalWater = ee.ImageCollection(myjrc.select('onlywater')).sum().toFloat();
    
    //统计每个像素点在一年中有水的比例
    var floodfreq = totalWater.divide(totalObs).multiply(100);
    
    //删除没有值的像素
    var myMask = floodfreq.eq(0).not();
    floodfreq = floodfreq.updateMask(myMask);
    
    var viz = {min:0, max:50, palette: ['ffffff', 'fffcb8', '0905ff']};
    var floodfreq1=floodfreq.clip(roi);
    var year_folder=year+"folder_gte";
    
    //如果某个像素一年有7个月有水，则为水体
    var gte60=floodfreq1.gte(60)
    //加载范围
    Map.addLayer(roi)
    //加载影像
    Map.addLayer(floodfreq.clip(roi),viz,year_folder)
    
    //导出影像
    Export.image.toDrive({
      image: gte60,
      region: roi,
      crs: "EPSG:4326",
      scale: 30,
      maxPixels : 1e13,
      folder:year_folder,
      description:year_folder});

     
    //计算计算水体面积
    var stats2 = gte60.reduceRegion({
      reducer: ee.Reducer.sum(),
      geometry: roi,
      scale: 250,
    });

    print(year_folder);
    print(stats2);
}

结果展示：