Generating a raster map presentation of a forest resource by solving a transportation problem

Necessary tools for raster map generation, for the approach based on the calibration estimator, were developed and implemented. The allocation of the area weight of each pixel to sample plots was formulated as a transportation problem, using a spectral distance measure as a transportation cost, and solved using the transportation simplex algorithm. Pixel level accuracy was calculated for the methods based on the calibration estimator so that the results could be compared with the results of the nearest neighbour estimation, the reference sample plot method (RSP) at pixel level. Local averaging in a 3 × 3 window was performed for each generated raster map as a postprocessing phase to smooth the map. Test plot results were calculated both for the unfiltered raster map and the filtered raster map. RSP produced the smallest RMSE in the pooled test data. Local averaging with a 3 × 3 filter decreased the pixel level error – and the bias – and the differences between the methods are smaller. Without local averaging, the pixel level errors of the methods based on solving the transportation problem were high. Raster map generation using the methods of this study forms an optional part – followed possibly by the classification of the pixel level results – of the whole computation task, when the area weight computation is based on the calibration estimation. For larger areas than in the present study, such as municipalities, the efficiency of the method based on the transportation model must be improved before it is a usable tool, in practice, for raster map generation. For nearest neighbour methods, the area size is not such a problem, because the inventory area is processed pixel by pixel.

[1]  Mats Nilsson,et al.  Simultaneous use of Landsat-TM and IRS-1C WiFS data in estimating large area tree stem volume and aboveground biomass , 2002 .

[2]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[3]  Erkki Tomppo,et al.  Calibration of small-area estimates for map errors in multisource forest inventory. , 2000 .

[4]  Arto Haara Kasvuennusteiden luotettavuuden selvittäminen knn-menetelmällä ja monitavoiteoptimoinnilla , 1970 .

[5]  Timo Tokola,et al.  The Influence of Field Sample Data Location on Growing Stock Volume Estimation in Landsat TM-Based Forest Inventory in Eastern Finland , 2000 .

[6]  C. Särndal,et al.  Calibration Estimators in Survey Sampling , 1992 .

[7]  Albert R. Stage,et al.  Most Similar Neighbor: An Improved Sampling Inference Procedure for Natural Resource Planning , 1995, Forest Science.

[8]  Erkki Tomppo,et al.  Stratification by ancillary data in multisource forest inventories employing k-nearest-neighbour estimation , 2002 .

[9]  Annika Kangas,et al.  Calibrating predicted diameter distribution with additional information. , 2000 .

[10]  E. Tomppo Multi-source national forest inventory of Finland. , 1994 .

[11]  Ronald E. McRoberts,et al.  Stratified estimation of forest area using satellite imagery, inventory data, and the k-Nearest Neighbors technique , 2002 .

[12]  N. Altman An Introduction to Kernel and Nearest-Neighbor Nonparametric Regression , 1992 .

[13]  Juho Pitkänen,et al.  Point accuracy of a non-parametric method in estimation of forest characteristics with different satellite materials , 1996 .

[14]  E. Tomppo,et al.  Selecting estimation parameters for the Finnish multisource National Forest Inventory , 2001 .

[15]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[16]  Erkki Tomppo,et al.  Kunnittaiset metsävaratiedot 1990–94 , 1970 .

[17]  J. Lappi Forest inventory of small areas combining the calibration estimator and a spatial model , 2001 .

[18]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[19]  Hamdy A. Taha,et al.  Operations research: an introduction / Hamdy A. Taha , 1982 .

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .