Problems with NDVI value from PlanetScope

Hi,

You are actually downloading the “visualisation” of NDVI as you can see in your screenshot. To download the actual NDVI value, you need to set your own configuration with a proper layer, i.e., a layer that returns NDVI value as the output.

Here is a webinar that should give you a well introduction of EO Browser, and you can find a step-by-step instruction to create a configuration in the Create your own configurations section.

Hi,
Thanks a lot for your explanation.

Hi 
I tried to create my own configuration as you suggested.
I followed the tutorial that can be found in FAQ.
Then I selected the script for the Planet Scope NDVI that can be found here:

Sentinel-Hub custom scripts

Normalized difference vegetation index, PlanetScope

A repository of custom scripts that can be used with Sentinel-Hub services.

But this is the result:

Can you help me to understand where is my error?
Thanks a lot

Hi,

Sorry for the confusion. PlanetScope does not have B3 and B4 bands. The available bands are listed in our documentation. We will fix the example script in the custom script repository.

Since EO Browser is designed for visualisation, we need to manipulate the NDVI index a little bit so the index can fit into EO Browser. The idea is that we multiply the index by 10000 and add 10000 to make the index range from 0 to 20000. Please check out this blog post for more details.

The following Evalscript returns NDVI as UINT16.

//VERSION=3

function setup() {
  return {
    input: ["Red", "NIR", "dataMask"],
    output: { bands: 1, sampleType: "UINT16" }
  };
}

function evaluatePixel(sample) {
  let ndvi = index(sample.NIR , sample.Red);
  return [10000 * ndvi+ 10000];
}

Once the images are downloaded, you can easily rescale back to the range -1 to 1 by applying the following operation to your image: ndvi_float = (ndvi-10000)/10000.

Thanks a lot. 
Does this process works also for NDWI?

Hi,

Yes, it will work with NDWI too. You would just need to use the Green band instead of the Red band when calculating the index.

//VERSION=3

function setup() {
return {
input: [“Green”, “NIR”, “dataMask”],
output: { bands: 1, sampleType: “UINT16” }
};
}

function evaluatePixel(sample) {
let ndwi = index(sample.NIR , sample.Green);
return [10000 * ndwi+ 10000];

So the layer will be something like this and to calculate the real value of the index on QGIS my formula will be: ndwi_float = (ndwi-10000)/10000 .
Am I right?

Yes, that’s correct.