Skip to Main Content
Comparison of regression and geostatistical methods for mapping Leaf Area Index (LAI) with Landsat ETM+ data over a boreal forest.Author(s): Mercedes Berterretche; Andrew T. Hudak; Warren B. Cohen; Thomas K. Maiersperger; Stith T. Gower; Jennifer Dungan
Source: Remote Sensing of Environment. 96: 49-61.
Publication Series: Scientific Journal (JRNL)
Station: Rocky Mountain Research Station
PDF: Download Publication (1006.07 KB)
DescriptionThis 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.
- You may send email to email@example.com to request a hard copy of this publication.
- (Please specify exactly which publication you are requesting and your mailing address.)
- 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.
CitationBerterretche, 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.
Keywordsconditional simulation, kriging, reduced major axis regression
- Integration of lidar and Landsat ETM+ data for estimating and mapping forest canopy height.
- An improved strategy for regression of biophysical variables and Landsat ETM+ data.
- Estimation of crown biomass of Pinus pinaster stands and shrubland above-ground biomass using forest inventory data, remotely sensed imagery and spatial prediction models
XML: View XML