基于GEE可视化和下载Sentinel2 L2A数据

最新推荐文章于 2025-05-01 08:00:00 发布

燕南路GISer

最新推荐文章于 2025-05-01 08:00:00 发布

阅读量538

点赞数 3

CC 4.0 BY-SA版权

文章标签： python

原文链接：https://blog.csdn.net/m0_56729804/article/details/135970429

本文介绍了如何使用GoogleEarthEngine处理Landsat8卫星数据，包括去云、数据预处理、多光谱图像合成，并将处理后的图像导出到GoogleDrive。作者展示了如何使用maskL8sr函数和Map组件进行地理空间分析。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//作者：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,8) //将地图中心设置为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('2018-04-01', '2018-12-30')
                  .filter(ee.Filter.lte('CLOUD_COVER',10))
                  .map(maskL8sr)
                  .median()
                  .clip(roi)//裁剪
                  
print(dataset);
 
var visualization = {
  bands: ['SR_B4', 'SR_B3', 'SR_B2'],
  min: 0.0,
  max: 0.3,
};
Map.addLayer(dataset, visualization, 'Image');
// 彩色合成显示
 
 
var mergedImage = dataset.select("SR_B1","SR_B2","SR_B3","SR_B4","SR_B5","SR_B6","SR_B7");// "ST_B10"
// 将所选波段合并为一个多波段图像
print(mergedImage)
 
function exportdata(data1) {
  var data = data1.toFloat()
  Export.image.toDrive({
  image: data,//要下载的影像，类型为image对象
  folder: 'Landsat8-L2A',//云盘文件夹
  description: 'AH',
  fileNamePrefix: 'Landsat8-OLI',
  region: roi,
  scale: 30,//分辨率
  crs: "EPSG:4326",//投影坐标系
  maxPixels: 1e13//最大像元数
  });
}
exportdata(mergedImage);
//导出数据