看到身边不少同学,都喜欢用GEE来下载遥感影像数据。一是因为非常便捷,二是因为GEE平台上的遥感影像多数已经过地形校正或大气校正,使用起来非常喜人。Hence,本期介绍如何在GEE上下载哨兵2号影像数据

//////define your study area here!
var studyarea = area1;
Map.addLayer(studyarea)

////////////////////////////prepare sentinel imagecollection/////////////////////////////////////
/**
 * Function to mask clouds using the Sentinel-2 QA band
 * @param {ee.Image} image Sentinel-2 image
 * @return {ee.Image} cloud masked Sentinel-2 image
 */
function maskS2clouds(image) {
  var qa = image.select('QA60');
  // Bits 10 and 11 are clouds and cirrus, respectively.
  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;
  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));
  return image.updateMask(mask).divide(10000);
}  //去云处理


var dataset = ee.ImageCollection('COPERNICUS/S2_SR')
                  .filterDate('2020-01-01', '2020-12-31')   //影像时间
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE',20))  //云量
                  .map(maskS2clouds)  //去云
                  .map(function(image) {
                   var ndvi = image.normalizedDifference(['B8', 'B4']).rename('B5');
                   return image.addBands(ndvi);
                   })   //添加指数
                   .map(function(image){
                   var B1 = image.select("B2").rename("B1");//blue
                   var B2 = image.select("B3").rename("B2");//green
                   var B3 = image.select("B4").rename("B3");//red
                   var B4 = image.select("B8").rename("B4");//nir
                   return image.addBands(B1)
                               .addBands(B2)
                               .addBands(B3)
                               .addBands(B4)
                   })  //重命名
                   .select(["B1","B2","B3","B4","B5"/*,"B6","B7","B8","B9"*/]);  //影像集
                   
var study_img = dataset.median().clip(studyarea);  //中值影像,裁剪出研究区影像
Map.addLayer(study_img, {bands:["B4", "B3", "B2"], min:0, max:0.3}, 'RGB');

//导出哨兵2号影像数据
Export.image.toDrive({
  image:study_img.select(["B1","B2","B3","B4","B5"]),
  description:"sentinel2_AREA",
  region:studyarea,
  scale:10,
  maxPixels:1e13
})  
Logo

DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。

更多推荐