GEE加载哨兵数据以及ESRI、EAS、FROM地表覆被产品

因为需要对比查看ESRI、EAS、FROM数据在水体提取上的效果
所以下面是用GEE加载出哨兵数据和这三种土地覆被数据的代码:
esri2020年地表覆被数据:

// 边界数据

var core_area = ee.FeatureCollection("users/selenebai1998/core_middle");

Map.centerObject(core_area,4);

Map.addLayer(core_area,{},'shp');


/*****************sentinel2***************************/
//去云

function maskS2clouds(image)
{
  var qa = image.select('QA60');
 
  // Bits 10 and 11 are clouds and cirrus, respectively.
  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;
 
  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));
 
  return image.updateMask(mask).divide(10000);
}

//处理数据,生成数据集
var dataset = ee.ImageCollection('COPERNICUS/S2')
                  .filterDate('2020-01-01', '2020-12-31')
                  .filterBounds(core_area)
                  // 选择云量
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 90))
                  .map(maskS2clouds) ;            
                
var image = dataset.median().clip(core_area)                 
//可视化参数,按843波段合成                
var rgbVis = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};

//按中值合成指定指定时间段的影像
Map.addLayer(image, rgbVis, 'Sentinel2-image');



/*****************esri***************************/

// 裁剪函数
function clip(Image){
  return Image.clip(core_area)
}
 
// 加载数据集
var esri_lulc10 = ee.ImageCollection("projects/sat-io/open-datasets/landcover/ESRI_Global-LULC_10m").map(clip)
 
 
// 定义字典和相应的颜色
var dict = {
  "names": [
    "Water",
    "Trees",
    "Grass",
    "Flooded Vegetation",
    "Crops",
    "Scrub/Shrub",
    "Built Area",
    "Bare Ground",
    "Snow/Ice",
    "Clouds"
  ],
  "colors": [
    "#1A5BAB",
    "#358221",
    "#A7D282",
    "#87D19E",
    "#FFDB5C",
    "#EECFA8",
    "#ED022A",
    "#EDE9E4",
    "#F2FAFF",
    "#C8C8C8"
  ]};
 
// 加载影像到地图上,这里的颜色用的就是上面定义的字典

Map.addLayer(esri_lulc10.mosaic(), {min:1, max:10, palette:dict['colors']}, 'ESRI LULC 10m')


/**********************ESA*************************/
// 加载数据集
var ESA = ee.ImageCollection("ESA/WorldCover/v100").map(clip)
 
 
// 定义字典和相应的颜色
var imageVisParam_ESA = {"opacity":1,"bands":["Map"],"min":10,"max":110,"palette":["8dffda","14ff31","686dff","ff9b54","ff712d","ff66b4"]};
 
// 加载影像到
ESA = ESA.mosaic()

Map.addLayer(ESA,imageVisParam_ESA,"ESA_landcover")

//下载数据
/*
Export.image.toDrive({
      image:IMAGE,
      description:'image',
      scale:10,
      maxPixels: 1e13,
      region:district_geometry,
      fileFormat: 'GeoTIFF',
      formatOptions: {
        cloudOptimized: true
      }
  });
  */

  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值