X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 13 May 2015 10:13 AM by  anon
Vegetation Indices - output doesn't make sense.
 2 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
13 May 2015 10:13 AM
    Hello, Apologies if these are simple questions, I'm new to ENVI and struggling to understand much of the help online! I'd appreciate any advice on any of these issues. I am trying to calculate indices to look at vegetation damage/stress over time using landsat 8 images with the bands stacked. I've used both the band math and spectral index tools to try and calculate Enhanced Vegetation Index. In the band math tool I entered the expression (2.5*((B2-B1)/(B2+6*B1-7.5*B3+1))) where B1= red B2=NIR and B3=blue. In both cases the result is very strange, with values ranging from -10000 to +50000, even though the range should be -1 to +1. I'm having a similar problem with the Green Vegetation Index, getting values of several thousand. However my NDVI and NBR indexes, calculated using the NDVI tool, are fine. I'm at a loss as to how to work out what I've done wrong. In addition two quicker questions: When manually entering wavelengths using edit ENVI header do you just enter the lowest value? i.e. 0.43 for OLI band 1 (0.43-0.45)? I couldn't find a way to enter a range. Is it possible to calculate indexes like the Photochemical Reflectance Index (=(531nm-570nm)/(531+570nm)) with this sort of imagery? Many thanks for your time.

    Deleted User



    Basic Member


    Posts:228
    Basic Member


    --
    13 May 2015 11:54 AM
    Hi Rachael, It sounds like you are running into an issue with data type in your Band Math expression. If your input data and values in the expression are integer data type (versus floating point values), then your expression will be calculated using integer math, which won't allow the output to be anything other than integers. Which is going to end up looking like nonsense. In this case, you really want to do the calculation using all floating point math, so that you can get floating point (decimal) values in the output. The easiest way to indicate that in your Band Math expression is to change it so that the integers in your equation are all floating point values, by adding .0 to the end of them. So it will look like this: (2.5*((B2-B1)/(B2+6.0*B1-7.5*B3+1.0))) Does that give you the values you expect? When manually entering wavelengths in the ENVI header, I think it is best to enter the center of the wavelength range for the band. ENVI doesn't really care what you use. It will interpret and report the pixel values as being for whatever specific wavelength you supply. Mainly this is used when viewing pixel spectra, because you need to know which wavelength to assign each pixel value to. With Landsat data, you don't have measurements at 531 nm and then again at 570 nm. Instead you have one band (the visible green band) that covers that whole range, from about 520 nm to about 600 nm. So the equation for the PRI isn't meaningful in that situation. Regards, Peg Exelis VIS

    Deleted User



    New Member


    Posts:59
    New Member


    --
    13 May 2015 12:18 PM
    Hi, I believe the help is incorrect where it says, "Values range from -1 to 1." I would expect that range with the normalized difference indices (x-y)/(x+y), but not with the EVI equation. The ranges that you see should be correct. We may need to add some details to the help on how to interpret the data values of EVI.
    You are not authorized to post a reply.