Sentinel-2 with cloudy parts replaced by Sentinel-1

Show script or download it.

      // VERSION=3
// SAR data substitution for clouded Sentinel-2 images V1.0.0
// Twitter: Pierre Markuse (@pierre_markuse)
// CC BY 4.0 International - https://creativecommons.org/licenses/by/4.0/

var setup = () => ({
  input: [
    {datasource: "s2l1c", bands:["B02", "B03", "B04","CLM", "CLP"], units: "REFLECTANCE", mosaicking: "ORBIT"},
    {datasource: "s1grd", bands:["VV", "VH"], units: "REFLECTANCE", mosaicking: "ORBIT"},
  ],
  output: [
    { id: "default", bands: 3, sampleType: SampleType.AUTO }
  ],
})

function evaluatePixel(samples, inputData, inputMetadata, customData, outputMetadata) {
 var S2L1C = samples.s2l1c[0];
 var S1 = samples.s1grd[0];
 let SI = 54272; // Yes, I'm that old... 
 let WAT = 25; // Water Threshold for SAR
 let CLP = S2L1C.CLP / 2.55; // Cloud Propability
 let CLPT = 70; // Cloud Propabilty Threshold in percent
 
 
  
if ( (CLP > CLPT && S1.VV/S1.VH <= WAT)) 
 {
    return {
      default: [S1.VV*3.0,S1.VV*1.1+S1.VH*8.75,S1.VH*1.75]
    }}
  
  if ( (CLP > CLPT && S1.VV/S1.VH > WAT)) 
 {
    return {
      default: [S1.VV*1,S1.VV*8,0.5+S1.VV*3+S1.VH*2000]
    }}

  
    let val = [3*S2L1C.B04, 3*S2L1C.B03, 3*S2L1C.B02];
  return {
    default: val
  }
}

Evaluate and visualize

• EO Browser

General description

The script uses cloud masks to identify cloudy sentinel-2 areas and replace the data with terrain visualisation based on Sentinel-1 data.

Description of representative images

Sevastopol with cloudy areas in Sentinel-1, and non-cloudy areas in Sentinel-2.

Contributors:

• Pierre Markuse, (Twitter: @pierre_markuse)

License

• CC BY 4.0 International