【GEE】使用GEE批量查询下载Landsat8数据

        刚发了一篇Landsat8地表温度反演的博文,顺便分享一下如何使用GEE批量查询、下载Landsat8数据集。代码比较简单就是查询函数和导出函数,然后还有一个显示函数。网上的教程一大堆,都差不多的代码,在这里要感谢一些前辈们的无私奉献。

踩坑1:2021年之前的C01数据集的显示函数不能直接用到C02数据集上,不然影像显示为全白。

一、代码部分

        直接导入自己的研究区即可,无需修改其他内容。然后到Tasks里面就可以看到所有符合条件的影像集合,自己点击下载即可。

//作者:RS迷途小书童
//博客:https://blog.csdn.net/m0_56729804?type=blog
var roi = table //感兴趣的区域信息
var style_set = {color:"red",fillColor:"00000000"}; //设置地图中要素的颜色和填充颜色
Map.addLayer(roi.style(style_set),{},"shape") //使用之前定义的样式集将roi添加到地图中。该地图层默认使用几何形状(例如多边形)来表示区域
Map.centerObject(roi,10) //将地图中心设置为roi对象,并设置缩放级别为10

function maskL8sr(image) {
  // 去云函数
  var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
  //从输入图像中选择QA_PIXEL波段,使用位运算和掩码来识别填充、云、云影等像素
  var saturationMask = image.select('QA_RADSAT').eq(0); //从输入图像中选择QA_RADSAT波段,并识别未饱和的像素。
  // 将缩放因子应用于适当的频带
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  // 从输入图像中选择光学波段,并应用归一化处理。
  var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);
  // 从输入图像中选择热红外波段,并应用归一化处理。
  return image.addBands(opticalBands, null, true)
      .addBands(thermalBands, null, true)
      .updateMask(qaMask)
      .updateMask(saturationMask);}
var dataset = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
                  .filterBounds(roi)
                  .filterDate('2023-01-01', '2023-03-31')
                  .map(maskL8sr)
                  .map(function(image){
                      return image.clip(roi)//裁剪
                  })
                  
print("dataset",dataset)
var visualization = {
  bands: ['SR_B4', 'SR_B3', 'SR_B2'],
  min: 0.0,
  max: 0.3,
};
Map.addLayer(dataset, visualization, 'Image');
// 彩色合成显示

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.toFloat();
      // 20240724修改
      Export.image.toDrive({
        image: image,
        description: indexs[i],
        fileNamePrefix: indexs[i],
        folder: 'Landsat8',
        region: roi,
        scale: 30,
        crs: "EPSG:4326",
        maxPixels: 1e13
      });
    }
  });
}
exportImageCollection(dataset);

二、效果显示

        这个专栏开了之后也不知道怎么分享博文,感觉没什么好写的,不如直接放代码。但是有感觉光放代码有比较单调,很烦。等后期看看如何改进,大家有什么问题可以随时留言交流!

  • 12
    点赞
  • 54
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 43
    评论
GEE是Google Earth Engine的简称,它是一个基于云计算的地理信息处理平台。GEE提供了很多遥感数据集,包括Landsat 8卫星的数据。 要下载Landsat 8的数据,首先需要在GEE上注册一个账号。注册完成后,可以进入GEE数据目录,搜索Landsat 8数据集。在搜索结果中,可以找到Landsat 8 OLI (Operational Land Imager) 和 TIRS (Thermal Infrared Sensor)的数据。 找到Landsat 8数据集后,可以通过GEE的JavaScript API来编写代码进行数据下载。首先,需要定义感兴趣区域(ROI)并设置时间范围。然后,使用Landsat 8数据集的ID创建一个图像集合,通过filter()方法进行筛选,选取所需的波段。 代码示例: ```javascript var roi = ee.Geometry.Rectangle([lon1, lat1, lon2, lat2]); // 设置感兴趣区域的经纬度范围 var startDate = '2021-01-01'; // 设置开始时间 var endDate = '2021-12-31'; // 设置结束时间 var landsat8Collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') // 获取Landsat 8数据集 .filterBounds(roi) // 筛选感兴趣区域的数据 .filterDate(startDate, endDate); // 筛选时间范围的数据 print(landsat8Collection); // 查看筛选后的数据集 // 下载Landsat 8的数据 var landsat8 = ee.Image(landsat8Collection.first()); // 获取第一幅图像 // 设置下载参数 var downloadParams = { name: 'landsat8_data', // 下载文件的名称 scale: 30, // 数据的分辨率 crs: 'EPSG:4326' // 数据的投影 }; // 开始下载 Export.image.toDrive({ image: landsat8, description: 'landsat8_download', fileNamePrefix: 'landsat8_data', scale: 30, region: roi }); ``` 以上代码解释了如何使用GEE下载Landsat 8的数据。根据实际的需求,可以调整感兴趣区域、时间范围和下载参数等。GEE提供了丰富的功能和数据,可以满足用户在地理信息处理和遥感数据分析方面的需求。
评论 43
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

RS迷途小书童

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值