Google Earth Engine下载Sentinel2数据

var district = table;
var dsize = district.size();
var district_geometry = district.geometry();
print(dsize);
Map.centerObject(district_geometry,7);

//去云的范围内的范围内,
function maskS2clouds(image) {
var qa = image.select('QA60');

// Bits 10 and 11 分别是云和卷云
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")
.filterBounds(district_geometry)
.filterDate('2020-10-28', '2020-10-30')
// 选择云量真彩色真彩色,0到到100的范围内
.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
.map(maskS2clouds)
.select(['B2', 'B3', 'B4','B5','B6','B7','B8A','B11','B12']);

//可视化参数,按2,3,4波段合成真彩色
var rgbVis = {
min: 0.0,
max: 0.3,
bands: ['B2', 'B3', 'B4'],
};

//按中值合成指定指定时间段的影像
Map.addLayer(dataset.median(), rgbVis, 'RGB');

//设置显示样式:color代表边界颜色;fillcolor代表填充颜色
var styling={color:'red',fillColor:'00000000'}

//FH就是矢量边界,这里可以不用管
Map.addLayer(district.style(styling),{},"FH");

var IMAGE = dataset.mosaic();

//下载数据
Export.image.toDrive({
image:IMAGE,
description:'Sentinel2',
scale:30,
maxPixels: 1e13,
region:district_geometry,
crs: "EPSG:32649",
fileFormat: 'GeoTIFF',
formatOptions: {
cloudOptimized: true
}
});

posted @   mylifeforaiur  阅读(25)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· 没有Manus邀请码?试试免邀请码的MGX或者开源的OpenManus吧
· 【自荐】一款简洁、开源的在线白板工具 Drawnix
· 园子的第一款AI主题卫衣上架——"HELLO! HOW CAN I ASSIST YOU TODAY
点击右上角即可分享
微信分享提示