From DavinciWiki
Jump to: navigation, search


Generates a Modified Soil-Adjusted Vegetation Index map

Arguments and Return Values

Arguments: 1- or 2-dimensional greyscale images representing near-infrared and red spectral bands (both are required)

Return Value: A float greyscale image


Syntax: msavi(nir, red)

This function is part of the larger class of vegetation-index methods used in satellite image-processing. See ndi or this discussion of related functions for more detailed information.

MSAVI is of the subtype which will return both positive and negative values. Positive values generally indicate the presence of vegetation (with greater values indicating healthier vegetation). Negative values generally indicate a lack of vegetation (water, rock, soil, etc.).

The formula used is:

MSAVI = NIR + 0.5 - (0.5 * sqrt((2 * NIR + 1)^2 - 8 * (NIR - (2 * Red)))

This formula was obtained from "Above-Ground Biomass Estimation of Successional and Mature Forests Using TM Images in the Amazon Basin" (Dengsheng Lu, Paul Mausel, Eduardo Brondizio, and Emilio Moran, 2002).


dv> nir
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> r
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> result = msavi(nir, r)
Generating a Modified Soil-Adjusted Vegetation Index image.
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> max(result)

dv> min(result)


Full set of LandSat 7 imagery:

NYC-Band Comparison.PNG

Output of msavi():

Vegetation Index-NYC-MSAVI.PNG

DavinciWiki Mini-Nav Bar


Contact Developers

  • davinci-dev [AT]

All other topics

  • See navigation on the left

Recent Library Changes

Created On: 07-16-2012
Modified On: 05-02-2013

Personal tools