【GEE】下载研究区的Landsat8去云清晰影像

Twcat_tree

发布于 2024-04-10 08:44:31

1810

发布于 2024-04-10 08:44:31

文章被收录于专栏：二猫の家二猫の家

如果你要绘制研究区概况或者是做一个产品的比较

需要整个研究区的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
});