EVI calculation with Sentinel 2A data

Question

I'm using ArcGIS Pro 2.8.2. I've downloaded Sentinel 2A data, and would like to use raster calculator to calculate EVI. However, my results are outside of the -1 to 1 range, and the image it produces is very strange.

My code is:

2.5*(("B08"-"B04")/(("B08"+6*"B04"-7.5*"B02")+1))

where B02 is the blue band, B04 is the red band, and B08 is the NIR band. I based this code off of the general formula given here , which is:

However, when I enter this, I get a range of -18405 to 18220, and my image looks like:

Any advice on why this might be would be greatly appreciated.

JayantaPoddar · Accepted Answer

The formula considers the bands to be atmospherically corrected (L2A bands).

First you need to do an atmospheric correction using Sen2Cor (Check page No. 63, for the download option).

1. After download, extract the zip file in desired folder (no space in folderpath).

2. From Command Prompt, change directory to the extracted folder. Then execute a similar command as follows.

L2A_Process.bat --resolution 10 <Input Raster Folder path till *.SAFE folder>

This process may take a while.

Next, Add the desired bands from the newly created *L2A* folder (R10m subfolder) to ArcGIS Pro. Use these bands for the calculation as follows.

(2.5 * ("L2A_B08" - "L2A_B04") / (("L2A_B08" + 6.0 * "L2A_B04" - 7.5 * "L2A_B02") + 1.0))/10000

*My Area of Interest (here) is a desert with sparse vegetation.