Skip to Main Content
U.S. Forest Service
Caring for the land and serving people

United States Department of Agriculture

Home > Search > Publication Information

  1. Share via EmailShare on FacebookShare on LinkedInShare on Twitter
    Dislike this pubLike this pub
    Author(s): Mercedes Berterretche; Andrew T. HudakWarren B. Cohen; Thomas K. Maiersperger; Stith T. Gower; Jennifer Dungan
    Date: 2005
    Source: Remote Sensing of Environment. 96: 49-­61.
    Publication Series: Scientific Journal (JRNL)
    Station: Rocky Mountain Research Station
    PDF: Download Publication  (1006.07 KB)


    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.

    Publication Notes

    • You may send email to 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.


    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.


    conditional simulation, kriging, reduced major axis regression

    Related Search

    XML: View XML
Show More
Show Fewer
Jump to Top of Page