Linear regression modelling for the estimation of oil palm age from Landsat TM

This paper investigates the accuracy with which the age since field planting of oil palm (Elaeis guineensis Jacq.) can be estimated from Landsat Thematic Mapper (TM) radiance at pixel and stand scales. The study site, a commercial plantation 30 km south-east of Kuala Lumpur in Selangor, Malaysia, consisted of even-aged blocks from 4 to 21 years old. Spectral data were the six reflective TM bands and three spectral indices. Nonlinear negative relationships between spectral variables and age are compared to published trends in leaf area, stem height and per cent canopy cover for oil palm and other tree plantations. Correlation coefficients between log age and log radiance are moderate and highly significant (p<0.01) for bands 2-5 and 7 (-0.214 to-0.776) at the pixel scale, and increase at the stand scale (r 2=0.985 for log band 5, p<0.01). Relationships are strongest for the mid-infrared bands, especially band 5 (r 2=0.585, p <0.01) and the infrared index (IRI), a normalized difference index of bands 4 and 5 (r 2= 0.48, p<0.01). Direct and inverse linear regression models for log age with log band and log age with IRI squared (IRIsq) were constructed at both scales. Equivalent age was estimated from the models using independent test sets for differing scales and degrees of aggregation of the age classes. Single age classes cannot be estimated accurately at the pixel or stand scales; the lowest RMS error was obtained from the direct model using all bands (RMS error=3.9 years at pixel scale, 2.7 at stand scale). A posteriori aggregation into generalized age classes (<5, 6-10, 11-15, 16-21 years) improved the RMS error but the results were still unacceptably high (2.2, 2.3, 2.7, 6.0 years respectively for direct model 3 using all bands). Acceptable RMS errors down to 0.58 years were obtained for models using IRIsq with generalized age classes developed and applied at the stand scale when variations in ground cover and other variables were averaged out. The spatial pattern of error in equivalent age deserves investigation for precision crop management.

[1]  J. E. Pinder,et al.  Indications of Relative Drought Stress in Longleaf Pine from Thematic Mapper Data , 1999 .

[2]  Giles M. Foody,et al.  An assessment of radiance in Landsat TM middle and thermal infrared wavebands for the detection of tropical forest regeneration , 1996 .

[3]  Simulation des échanges radiatifs sur maquettes informatiques de Elaeis guineensis , 1994 .

[4]  Robert B. Waide,et al.  Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data , 1989 .

[5]  J. A. Schell,et al.  Monitoring vegetation systems in the great plains with ERTS , 1973 .

[6]  Kevin P. Price,et al.  Empirical relationships between structural and spectral factors of Yellowstone lodgepole pine forests , 1997 .

[7]  Jean-Philippe Gastellu-Etchegorry,et al.  Sensitivity of Texture of High Resolution Images of Forest to Biophysical and Acquisition Parameters , 1998 .

[8]  M. Fiorella,et al.  Determining successional stage of temperate coniferous forests with Landsat satellite data , 1993 .

[9]  Jennifer L. Dungan,et al.  Seasonal LAI in slash pine estimated with landsat TM , 1992 .

[10]  C. Goulding,et al.  Estimation of timber volume in a coniferous plantation forest using Landsat TM , 1997 .

[11]  W. Cohen,et al.  Estimating structural attributes of Douglas-fir/western hemlock forest stands from Landsat and SPOT imagery , 1992 .

[12]  Frank W. Davis,et al.  Thematic mapper analysis of tree cover in semiarid woodlands using a model of canopy shadowing , 1991 .

[13]  Robin L. Graham,et al.  Estimating forest productivity with Thematic Mapper and biogeographical data , 1989 .

[14]  Steven W. Running,et al.  Remote sensing of temperate coniferous forest leaf area index The influence of canopy closure, understory vegetation and background reflectance , 1990 .

[15]  J. McMorrow,et al.  Monitoring tropical rain forest conversion to estate agriculture in Sabah Malaysia using Landsat images and GIS , 1996 .

[16]  Eric R. Ziegel,et al.  Applied Multivariate Data Analysis , 2002, Technometrics.

[17]  H. Larsson,et al.  Linear regressions for canopy cover estimation in Acacia woodlands using Landsat-TM, -MSS and SPOT HRV XS data , 1993 .

[18]  P. Chavez.  Radiometric calibration of Landsat Thematic Mapper multispectral images , 1989 .

[19]  B. Markham,et al.  Radiometric Calibration of Landsat , 1997 .

[20]  A. U. Lubis,et al.  A preliminary assessment of the possibilities of using spatial remote sensing to study developments on an oil palm plantation in north Sumatra. , 1990 .

[21]  S. Ekstrand,et al.  Landsat TM-based forest damage assessment : correction for topographic effects , 1996 .

[22]  P. J. Curran,et al.  The importance of measurement error for certain procedures in remote sensing at optical wavelengths , 1986 .

[23]  F. Gemmell,et al.  Effects of forest cover, terrain, and scale on timber volume estimation with Thematic Mapper data in a rocky mountain site , 1995 .

[24]  P. Chavez.  An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data , 1988 .

[25]  Janet Franklin,et al.  Thematic mapper analysis of coniferous forest structure and composition , 1986 .

[26]  Philip Lewis.  Three-dimensional plant modelling for remote sensing simulation studies using the Botanical Plant Modelling System , 1999 .

[27]  J. Colby.  Land cover classi ® cation using Landsat TM imagery in the tropical highlands : the in ̄ uence of anisotropic re ̄ ectance , 2000 .

[28]  Jerome K. Vanclay,et al.  Utility of Landsat Thematic Mapper data for mapping site productivity in tropical moist forests , 1990 .

[29]  Gautam Badhwar,et al.  Spectral Characterization of Biophysical Characterstics in a Boreal Forest: Relationship between Thematic Mapper Band Reflectance and Leaf Area Index for Aspen , 1986, IEEE Transactions on Geoscience and Remote Sensing.

[30]  Julia Mcmorrow,et al.  Relation of oil palm spectral response to stand age , 1995 .

[31]  B. Markham,et al.  Radiometric properties of U.S. processed Landsat MSS data , 1987 .

[32]  Robert De Wulf,et al.  Extraction of forest stand parameters from panchromatic and multispectral SPOT-1 data , 1990 .

[33]  Alan H. Strahler,et al.  Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: effect of crown shape and mutual shadowing , 1992, IEEE Trans. Geosci. Remote. Sens..