Google Earth Engine 中 Sentinel-2 图像的大气校正答案

【问题标题】：Atmospheric correction for Sentinel-2 imagery in Google Earth EngineGoogle Earth Engine 中 Sentinel-2 图像的大气校正
【发布时间】：2020-07-03 22:26:44
【问题描述】：

我想在 Google 地球引擎 (GEE) 中对 Sentinel-2 图像应用大气校正。我看到了用 Python 编写的 Sammurphy 代码，不幸的是它对我不起作用。我使用下面的代码（Java）尝试了暗像素减法方法，但它导致我感兴趣的区域上的图像完全变暗。

我是 Earth Engine 和 JavaScript 的新手。有没有人尝试过在 GEE 中对 Sentinel 2 图像使用暗像素减法或任何其他大气校正（最好是用 Java 编写的代码）？

var toa=maskedComposite1;
var thresh= 0.5;

var dark=findDarkPixels(toa, thresh)
print(dark);


//Function to find dark pixels from threshold on sum of NIR, SWIR1, & SWIR2 bands
//Returns classified image with binary [0,1] 'dark' band
// toa: Sentinel 2 image converted to surface radiance
// thresh: threshold (0.2 - 0.5) value for sum of NIR, SWIR1 & SWIR2 bands
function findDarkPixels(toa, thresh) {
  var darkPixels1 = toa.select(['B8','B11','B12']);
  var darkPixels = darkPixels1.reduce(ee.Reducer.sum()).lt(thresh);
  var filtered = darkPixels.focal_mode(0.1, 'square', 'pixels');
  Map.addLayer(filtered,{},'darkPixel');
  return filtered.rename(['dark']);
}

【问题讨论】：

标签： google-earth sentinel2

【解决方案1】：

如果您不需要特定的大气校正，则可以使用 GEE 中已有的 Level-2A Sentinel-2 数据。 Here is the link to dataset info. 此数据集的大气校正由 sen2cor 执行。请注意数据可用的时间段，因为 2A 级数据不适用于整个数据存档。

【讨论】：

感谢 Agris 的评论。目前，我正在尝试对不同空间和时间尺度的 S2 图像（使用 1C 级 -TOA）实施大气校正。因此，我认为最好自己应用大气校正，而不是直接使用表面反射率（2A 级-BOA），因为您指出的感兴趣区域和时间段的图像不可用。