Scripts for index images

arun · November 7, 2018, 10:53pm

I’m not familiar with JS, so want to clarify that NDVI and NDWI in the “Python scripts template” use the standard equations? I think I’ve also seen:
return [B8A - B05 / B8A - B05]

NDVI: (NIR - Red) / (NIR + Red)
Sentinel 2 - Band 8a (NIR) & Band 4 (R )
Landsat 8 - Band 5 (NIR) & Band 4 (R )

NDWI: (SWIR - NIR) / (SWIR + NIR)
Sentinel 2 - Band 8a (NIR) & Band 11 (SWIR)
Landsat 8 – Band 5 (NIR) & Band 6 (SWIR)

Or is the below script (S2 NDVI) just ratios of the bands; values look correct for the downloaded S2 image? Can I similarly put the S2 band names for NDWI? Do similar scripts for L8?

let viz = new Identity();

function evaluatePixel(samples) {
let val = index(samples[0].B8A, samples[0].B04);
return viz.process(val);
}

function setup(ds) {
setInputComponents([ds.B04, ds.B8A]);
setOutputComponentCount(1);
}

gmilcinski · November 7, 2018, 11:09pm

I am not sure, what your question is.
let val = index(samples[0].B8A, samples[0].B04);
is the same as
let val = (B8A-B04)/(B8A+B04);

You can find many configurations here:

arun · November 7, 2018, 11:15pm

Thanks Grega, you answered my question. If they are both same, why use the big script in the template? Or does the big script offer any advantages?

gmilcinski · November 8, 2018, 5:56am

In this specific case, for simple indices, there is no reason to use the “big script”.
There are however more complicated cases, where the structure is required.
E.g. multi-temporal processing (cloudless mosaics):

github.com

sentinel-hub/custom-scripts/blob/master/sentinel-2/cloudless_mosaic/L2A-first_quartille.js

/*
Script works on Sentinel-2 L2A data and requires scene classification (SCL) band. 
It takes one year of data, which is quite compute and time intensive, which is why it is recommended to run it on small area (e.g. 256x256 px).
An example of the results is New Zealand's cloudless mosaic, available here: https://data.linz.govt.nz/layer/93652-nz-10m-satellite-imagery-2017/

For the output value for each pixel it uses the first quartile value of valid values, each band separately. If there are none it uses invalid values instead. 
When using SCL its very important to use nearest neighbor resampling with a resolution of about 20m/px or more. 
*/

function    setup (dss) {
 setInputComponents([dss.B04,dss.B03,dss.B02,dss.SCL]);
 setOutputComponentCount(3);
}
function filterScenes (scenes, inputMetadata) {
   return scenes.filter(function (scene) {
     return scene.date.getTime()>=(inputMetadata.to.getTime()-12*31*24*3600*1000) ;
   });
}

function getValue(values) {

This file has been truncated. show original

or Leaf area index:

github.com

sentinel-hub/custom-scripts/blob/master/sentinel-2/lai/script.js

// Leaf Area Index script as implemented in SNAP but without input and output validation!
// Input/output values which are suspect are not reported or changed. Most values, however, do not fall under this category.
// Visualized as an interval from 0-3. This can be adjusted in the evaluatePixel method.

var degToRad = Math.PI / 180;

function evaluatePixel(samples) {
  var sample = samples[0];
  var b03_norm = normalize(sample.B03, 0, 0.253061520471542);
  var b04_norm = normalize(sample.B04, 0, 0.290393577911328);
  var b05_norm = normalize(sample.B05, 0, 0.305398915248555);
  var b06_norm = normalize(sample.B06, 0.006637972542253, 0.608900395797889);
  var b07_norm = normalize(sample.B07, 0.013972727018939, 0.753827384322927);
  var b8a_norm = normalize(sample.B8A, 0.026690138082061, 0.782011770669178);
  var b11_norm = normalize(sample.B11, 0.016388074192258, 0.493761397883092);
  var b12_norm = normalize(sample.B12, 0, 0.493025984460231);
  var viewZen_norm = normalize(Math.cos(sample.viewZenithMean * degToRad), 0.918595400582046, 1);
  var sunZen_norm  = normalize(Math.cos(sample.sunZenithAngles * degToRad), 0.342022871159208, 0.936206429175402);
  var relAzim_norm = Math.cos((sample.sunAzimuthAngles - sample.viewAzimuthMean) * degToRad)

This file has been truncated. show original