Skip to content
Snippets Groups Projects
Commit 453a44b4 authored by TGermain's avatar TGermain
Browse files

evi change for baseline processing

parent e4cfe004
No related branches found
No related tags found
No related merge requests found
......@@ -622,11 +622,11 @@ def create_raw_evi(
) # ignore warnings when dividing by zero
if version.parse(baseline) >= version.parse("4.0"):
logger.info('Baseline >= 4.0 detected, changing DN calculation')
evi = (2.5 * ((nir-1000)/10000 - (red-1000)/10000) / ((nir-1000)/10000 + 6.0 * (red-1000)/10000 - 7.5 * (blue-1000)/10000) + 1).astype(np.int32)
evi = (2.5 * (((nir-1000)/10000) - ((red-1000)/10000)) / (((nir-1000)/10000) + 6.0 * ((red-1000)/10000) - 7.5 * ((blue-1000)/10000) + 1)*10000).astype(np.int32)
evi_masked = np.where(red != 0, evi, 32767)
else:
logger.info('Baseline <= 4.0 detected, sticking with original DN calculation')
evi = (2.5 * ((nir/10000 - red/10000)) / ((nir/10000 + 6.0 * red/10000 - 7.5 * blue/10000) + 1)).astype(np.int32)
evi = (2.5 * ((nir/10000 - red/10000)) / ((nir/10000 + 6.0 * red/10000 - 7.5 * blue/10000) + 1)*10000).astype(np.int32)
evi_masked = np.where(red != 0, evi, 32767)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment