Sentinel-2数据目前在GEE中共有两种产品:一种是普通的没有经过大气校正的数据;另一种是经过大气校正后的数据(不全)。这两种数据去云的操作方法都是通过QA波段来完成的。
代码
var roi = /* color: #0b4a8b */ee.Geometry.Polygon(
[[[115.22447911987308, 38.97207826950874],
[117.64147130737308, 39.07450027191289],
[117.55358068237308, 40.67633196985795],
[114.91686193237308, 40.576272152256934]]]);
Map.centerObject(roi, 7);
//利用Sentinel-2自带的质量波段QA60通过按位运算去除云或阴影
function rmCloud(image) {
var qa = image.select('QA60');
var cloudBitMask = 1 << 10;
var cirrusBitMask = 1 << 11;
var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
.and(qa.bitwiseAnd(cirrusBitMask).eq(0));
return image.updateMask(mask);
}
var rawImage = ee.Image("COPERNICUS/S2/20180118T031039_20180118T031037_T50SMJ");
var visParams = {
bands: ['B4', 'B3', 'B2'],
min: 0,
max: 3000
};
print("rawImage", rawImage);
Map.addLayer(rawImage, visParams, "rawImage");
var cleanImage = rmCloud(rawImage);
print("cleanImage", cleanImage);
Map.addLayer(cleanImage, visParams, "cleanImage");
结果
Aikox66: 比例因子0.01
hyh354: 博主你好,想问一下有没有该产品去云的的相关代码
p13870656002: 请问有下载土壤水分数据图像的代码吗?只有土壤水分指数下载的代码吗?
TR0727: 您好!请问您的问题解决了吗?我也想知道
m0_70538913: 你好,想问一下era5官网中再分析数据中地表感热通量小时数据是什么?怎么处理呀,因为我计算的月均感热通量是负数,最终得到的应该是正数,不知道怎么处理数据?