Title: Comparison of regression and geostatistical methods for mapping Leaf Area Index (LAI) with Landsat ETM+ data over a boreal forest.
Author: Berterretche, Mercedes; Hudak, Andrew T.; Cohen, Warren B.; Maiersperger, Thomas K.; Gower, Stith T.; Dungan, Jennifer
Source: Remote Sensing of Environment. 96: 49-61.
Description: This study compared aspatial and spatial methods of using remote sensing and field data to predict maximum growing season leaf area index (LAI) maps in a boreal forest in Manitoba, Canada. The methods tested were orthogonal regression analysis (reduced major axis, RMA) and two geostatistical techniques: kriging with an external drift (KED) and sequential Gaussian conditional simulation (SGCS). Deterministic methods such as RMA and KED provide a single predicted map with either aspatial (e.g., standard error, in regression techniques) or limited spatial (e.g., KED variance) assessments of errors, respectively. In contrast, SGCS takes a probabilistic approach, where simulated values are conditional on the sample values and preserve the sample statistics. In this application, canonical indices were used to maximize the ability of Landsat ETM+ spectral data to account for LAI variability measured in the field through a spatially nested sampling design. As expected based on theory, SGCS did the best job preserving the distribution of measured LAI values. In terms of spatial pattern, SGCS preserved the anisotropy observed in semivariograms of measured LAI, while KED reduced anisotropy and lowered global variance (i.e., lower sill), also consistent with theory. The conditional variance of multiple SGCS realizations provided a useful visual and quantitative measure of spatial uncertainty. For applications requiring spatial prediction methods, we concluded KED is more useful if local accuracy is important, but SGCS is better for indicating global pattern. Predicting LAI from satellite data using geostatistical methods requires a distribution and density of primary, reference LAI measurements that are impractical to obtain. For regional NPP modeling with coarse resolution inputs, the aspatial RMA regression method is the most practical option.
Keywords: conditional simulation, kriging, reduced major axis regression
View or Print this Publication (1,006 KB)
- We recommend that you also print this page and attach it to the printout of the article, to retain the full citation information.
- This article was written and prepared by U.S. Government employees on official time, and is therefore in the public domain.
Berterretche, Mercedes; Hudak, Andrew T.; Cohen, Warren B.; Maiersperger, Thomas K.; Gower, Stith T.; Dungan, Jennifer. 2005. Comparison of regression and geostatistical methods for mapping Leaf Area Index (LAI) with Landsat ETM+ data over a boreal forest. Remote Sensing of Environment. 96: 49-61.
Get the latest version of the Adobe Acrobat reader or Acrobat Reader for Windows with Search and Accessibility