我想在谷歌地球引擎(GEE(中对哨兵-2图像应用大气校正。我看到了用Python编写的Sammurphy代码,不幸的是它对我不起作用。我使用下面的代码(Java(尝试了暗像素减法,但它导致我感兴趣的区域上出现完全黑暗的图像。
我是地球引擎和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']);
}
如果您不需要特定的大气校正,则可以使用GEE中已有的2A级哨兵-2数据。这是数据集信息的链接。该数据集的大气校正由sen2cor执行。请注意,数据可用的时间段,因为 2A 级数据不适用于整个数据存档。