如果你要绘制研究区概况或者是做一个产品的比较
需要整个研究区的Landsat8影像的话
你可以使用下面代码
// 定义研究区域
var table = ee.FeatureCollection("projects/ee-1261423515/assets/EHEimage/EHE");
function rmCloud (image){
var qa = image.select('pixel_qa')
var cloudMask = qa.bitwiseAnd(1 << 3).eq(0)
var cloudShadowMask = qa.bitwiseAnd(1 << 5).eq(0)
var mask_all = cloudMask.and(cloudShadowMask)
return image.updateMask(mask_all)
}
//去云算法
// 设置日期范围
var startDate = '2020-06-01';
var endDate = '2020-09-01';
// 加载Landsat 8数据集
var landsat = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR')
.filterDate(startDate, endDate)
.filterBounds(table)
// .filter(ee.Filter.lte('CLOUD_COVER',5))
.map(rmCloud);
// 对图像进行融合
var meanImage = landsat.select('B4', 'B3', 'B2').median().clip(table);
// 可视化融合影像
Map.centerObject(table, 5);
Map.addLayer(table, {}, 'Study Area');
var RGB_show = {
min: 0.0,
max: 3000,
bands: ['B4', 'B3', 'B2'],
};
Map.addLayer(meanImage, RGB_show, 'Landsat 8 Merged Image');
// 导出融合影像到Google Drive
Export.image.toDrive({
image: meanImage,
description: 'Landsat8_Merged_2020',
folder: 'Landsat8',
region: table.geometry(),
scale: 30,
crs: 'EPSG:4326',
maxPixels: 1e13
});
整体效果还不错
注意
// 对图像进行融合 var meanImage = landsat.select(‘B4’, ‘B3’,
‘B2’).median().clip(table);
这里我选用的是RGB,如果你需要其他波段的话,请在这里添加
如果你仅仅是使用影像来可视化的话,追求没有云的效果
你可以适当的调整时间,来看看哪一张那个融合后的影像效果最好
// 设置日期范围 var startDate = ‘2020-06-01’; var endDate = ‘2020-09-01’;