【GEE】基于GEE批量下载Landsat8 L1C数据(整幅)

 

        之前发过一篇使用GEE下载Landsat8的文章,然后有很多小伙伴私信我各种问题,如L1C、L2数据代码怎么修改,如何镶嵌,如何去云、 如何裁剪等一系列问题。正好快过年了,手头的事也没有多少了,所以这两天整理了一下GEE的相关代码,后续会陆续发出来。

        今天给大家带来的是Landsat8 L1C数据的批量下载,裁剪已经做过了。导出的结果为分幅的研究区影像(即原始影像裁剪后,无镶嵌)。

一、代码部分

//作者:RS迷途小书童
//博客:https://blog.csdn.net/m0_56729804?type=blog
var shape = table;
var style_set = {color:"red",fillColor:"00000000"};
Map.addLayer(shape.style(style_set),{},"shape")
//加载矢量

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 FindData = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR')
                  .filterBounds(shape)
                  .filterDate('2021-01-01', '2021-12-31')
                  .filter(ee.Filter.lte('CLOUD_COVER',5))
                  .map(rmCloud)//执行去云函数
                  .map(function(image){
                      return image.clip(shape)//裁剪
                  })
print("FindData",FindData)
var Count =FindData.size();
print('Count:',Count); 
//查找数据,裁剪影像

var RGB_show = {
  min: 0.0,
  max: 3000,
  bands: ['B4', 'B3', 'B2'],
};
Map.addLayer(FindData, RGB_show, 'FindData');
Map.centerObject(shape,10)//缩放级别
//彩色合成显示

var mergedImage = FindData.select("B1","B2","B3","B4","B5","B6","B7","B10","B11");
// 将所选波段合并为一个多波段图像
print(mergedImage)
function exportImageCollection(imgCol) {
  var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"])
                        .get("list");
  indexList.evaluate(function(indexs) {
    for (var i=0; i<indexs.length; i++) {
      var image = imgCol.filter(ee.Filter.eq("system:index", indexs[i])).first();
      image = image.toInt16();
      Export.image.toDrive({
        image: image,
        description: indexs[i],
        fileNamePrefix: indexs[i],
        folder: 'Landsat8',
        region: shape,
        scale: 30,//分辨率
        crs: "EPSG:4326",//投影坐标系
        maxPixels: 1e13
      });
    }
  });
}
exportImageCollection(mergedImage);
//导出数据

二、效果显示

posted @ 2024-02-01 11:42  RS迷途小书童  阅读(109)  评论(0编辑  收藏  举报