Scale considerations for integrating forest inventory plot data and satellite image data for regional forest mapping

To the best of our knowledge, one or more authors of this paper were federal employees when contributing to this work. This is the publisher’s final pdf. The published article is copyrighted by Elsevier and can be found at: http://www.journals.elsevier.com/remote-sensing-of-environment. Supplementary Data can be found at: http://www.sciencedirect.com/science/article/pii/S003442571300343XThe integration of satellite image data with forest inventory plot data is a popular approach for mapping forest\ud vegetation over large regions. Several methodological choices regarding spatial scale, mostly related to spatial\ud resolution or grain, can profoundly influence forest maps developed from plot and imagery data. Yet often the\ud consequences of scaling choices are not explicitly addressed. Our objective was to quantify the effects of several\ud scale-related methods on map accuracy for multiple forest attributes, using a variety of diagnostics that address\ud different map characteristics, to help guide map developers and users. We conducted nearest-neighbor imputation\ud over a large region in the Pacific Northwest, USA, to investigate effects of imputation grain (single pixel or\ud kernel); inclusion of heterogeneous plots; accuracy assessment grain and extent; and value of k (k = 1 and\ud k = 5), where k is the number of nearest-neighbor plots. Spatial predictors were from rasters describing climate\ud and topography and a time-series of Landsat imagery. Reference data were from regional forest inventory plots\ud measured over two decades. All analyses were conducted at a spatial resolution of 30 m × 30 m. Effects of imputation\ud grain and heterogeneous plots on map accuracy were small. Excluding heterogeneous plots slightly improved\ud map accuracy and did not lessen the systematic agreement (AC[subscript SYS]) between our maps and observed\ud plot data. Accuracy assessment grain strongly influenced map accuracy: maps assessed with a multi-pixel\ud block were much more accurate than when assessed with a single pixel for almost all map diagnostics, but this\ud was an artifact of methods rather than reflecting real differences among maps. Unsystematic agreement\ud (AC[subscript UNS]) between our maps and plots, or random error, improved notably with increasing accuracy assessment\ud extent for all scaling methods, indicating that reliability of most map applications can be improved through coarsening\ud the map grain. Value of k strongly influenced map diagnostics. The k = 5 maps were better than k = 1\ud maps for local-scale accuracy, but at the cost of reduced AC[subscript SYS], and loss of variability and poor areal representation\ud of forest conditions over the study region. The k = 1 maps produced notably better predictions of the least\ud abundant forest conditions (early-successional, late-successional, and broadleaf). None of the scaling methods\ud were optimal for all map diagnostics. Nevertheless, given a variety of diagnostics associated with a range of\ud scaling options, map developers and map users can make informed choices about methods and resulting maps\ud that best meet their particular objectives, and we present some general guidelines in this regard. Most of our\ud findings are applicable to mapping with Landsat data in other forested regions with similar forest inventory\ud data, and to other methods for spatial prediction