使用GEE对Sentinel-2A数据进行重采样、波段合成、导出并下载

//选择数据集
var geometry = ee.FeatureCollection('users/lijingyu990820/wudinghe')
Map.centerObject(geometry,14);
var sentinel2 = ee.ImageCollection('COPERNICUS/S2_SR')//2A
//'COPERNICUS/S2'为1C
                  .filterDate('2020-03-01', '2020-12-31')
                  // 预先过滤以获得较少的云量。小于10%
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE',10))
                  .filterBounds(geometry);
var imageMedian = sentinel2.median().clip(geometry);

// 定义重采样函数
var resampleTo10m = function(image) {
  // 重采样到 10 米分辨率,使用双线性插值
  var resampledImage = image.resample('bilinear').reproject({
    crs: 'EPSG:4326', //指定目标投影
    scale: 10
  });
  return resampledImage;
};

// 使用重采样函数对单个图像进行重采样
var resampledImage = resampleTo10m(imageMedian);


// 合并所有波段
// 选择要导出的波段
var bandB2 = resampledImage.select('B2');
var bandB3 = resampledImage.select('B3');
var bandB4 = resampledImage.select('B4');
var bandB5 = resampledImage.select('B5');
var bandB6 = resampledImage.select('B6');
var bandB7 = resampledImage.select('B7');
var bandB8 = resampledImage.select('B8');
var bandB9 = resampledImage.select('B9');
var bandB11 = resampledImage.select('B11');
// 将所选波段合并为一个多波段图像
var mergedImage = bandB2.addBands([bandB3, bandB4,bandB5, bandB6,bandB7, bandB8,bandB9, bandB11]);
//可视化 
Map.addLayer(mergedImage, {bands: ['B4', 'B3', 'B2'], min: 0,max: 3000,gamma: 1.4,}, "432");
Map.addLayer(geometry);
//Map.addLayer(dataset.median().clip(geometry), visualization, 'RGB');


// 导出图像
Export.image.toDrive({
  image: mergedImage,
  description: 'sentinel2_merged',
  crs: 'EPSG:4326',
  scale: 10,
  region: geometry,
  maxPixels: 1e13,
  folder: 'sentinel-2'
});

  • 7
    点赞
  • 26
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
以下是使用GEE(Google Earth Engine)进行Sentinel-1数据反演土壤湿度的代码示例: 1.导入Sentinel-1数据: ``` var sentinel1 = ee.ImageCollection('COPERNICUS/S1_GRD') .filterBounds(geometry) .filterDate('2018-01-01', '2021-12-31') .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')) .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')) .filter(ee.Filter.eq('instrumentMode', 'IW')) .filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING')) .select(['VV', 'VH']) ``` 2.计算土壤湿度: ``` var s1 = ee.Image(sentinel1.first()).clip(geometry) var vh = s1.select('VH') var vv = s1.select('VV') var b = 0.7 //地表粗糙度系数 var alpha = 1 //极化反照率系数 var eps = 0.015 //介电常数 //计算表面散射系数 var gamma_naught_vh = ee.Image.constant(10).multiply(vh).log10().multiply(10).subtract(83) var gamma_naught_vv = ee.Image.constant(10).multiply(vv).log10().multiply(10).subtract(83) //计算比例系数 var rho = gamma_naught_vh.divide(gamma_naught_vv) //计算土壤湿度 var soil_moisture = ee.Image.constant(1.5).multiply(rho.divide(Math.sin(Math.toRadians(40))).add(1)).log10().multiply(-10/b).divide(alpha*eps) ``` 3.可视化土壤湿度: ``` var vizParams = { min: -1, max: 1, palette: ['blue', 'white', 'green'] } Map.addLayer(soil_moisture, vizParams, 'soil moisture') ``` 注意:以上代码仅为示例,具体使用时需要根据实际情况进行修改。同时,需要注意在使用GEE时要了解相关的API和限制条件。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值