因为需要对比查看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
}
});
*/