msavi

From DavinciWiki
Jump to: navigation, search

Description

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

Usage

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).

Examples

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)
1.000001431

dv> min(result)
-0.8502268791

Images

Full set of LandSat 7 imagery:

NYC-Band Comparison.PNG

Output of msavi():

Vegetation Index-NYC-MSAVI.PNG

DavinciWiki Mini-Nav Bar

Contents


Contact Developers

  • davinci-dev [AT] mars.asu.edu

All other topics

  • See navigation on the left

Recent Library Changes

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

Personal tools