Spectral normalization of SPOT 4 data to adjust for changing leaf phenology within seasonal forests in Cambodia

Abstract As cloud cover exacerbates the application of optical satellite data for forest monitoring in tropical wet and dry regions during the rainy season, data acquisition is mainly restricted to the dry season. When analyzing wide areas, large numbers of single scenes obtained at different times of the dry season are often handled. Such imagery is characterized by changes of spectral reflectance due to vegetation phenology, varying atmospheric effects and solar geometries. In order to allow batch processing with automatic classification techniques, inter-scene comparability is required and data have to be radiometrically normalized. Cambodia is characterized by a mixture of evergreen, semi-evergreen and deciduous forest types, the latter two experiencing at least partial leaf shedding over the course of the dry season. Using spatial medium resolution SPOT 4 data and a manually delineated base map a season adjustment model was developed. The model is adapting the land cover specific spectral signatures of a slave scene (acquired in the middle of the dry season with its seasonal forests defoliated) to an adjacent master scene (from the beginning of the dry season, showing the same forest types with leafs). The relative position of every pixel reflectance was determined in relation to the mean reflectance and its standard deviation for each land cover type and sensor band of the unadjusted slave scene. For seasonality adjustment these pixel reflectance values were transformed (rescaled) to the corresponding position in spectral space defined by the band mean reflectance and standard deviation derived from the corresponding land cover class of the master scene. While the variability of spectral profiles of the pixels in the slave scene is rescaled, the mean reflectance value of the land cover class in the slave scene is conformed to the mean reflectance of the corresponding land cover class in the master scene. The Transformed Divergence (TD) separability index was used to indicate the performance of the adjustment process by characterizing the spectral distance for each land cover type comparing a reference dataset to the uncorrected and to the seasonality corrected scene respectively. While the TD values of all forest types showed a sharp decline, highlighting the good performance of the model, the TD values of the agriculture/urban class remained high, indicating limited normalization of this heterogeneous land cover type. In order to further demonstrate the performance of the model, an object-based land cover classification was applied to the unadjusted as well as to the corresponding adjusted scene. A comparison of the results showed a highly significant improvement of overall accuracy from 32.2% to 75.8% when applying seasonality adjustment.

[1]  J. Cihlar,et al.  Radiometric normalization of multitemporal high-resolution satellite images with quality control for land cover change detection , 2002 .

[2]  T. McMahon,et al.  Updated world map of the Köppen-Geiger climate classification , 2007 .

[3]  A. Gillespie,et al.  Topographic Normalization of Landsat TM Images of Forest Based on Subpixel Sun–Canopy–Sensor Geometry , 1998 .

[4]  Alice G. Laborte,et al.  Spectral Signature Generalization and Expansion Can Improve the Accuracy of Satellite Image Classification , 2010, PloS one.

[5]  D. Jackson,et al.  Trends in Global Cloud Cover in Two Decades of HIRS Observations , 2005 .

[6]  L. Curran,et al.  Utility of Landsat 7 satellite data for continued monitoring of forest cover change in protected areas in Southeast Asia , 2006 .

[7]  Makoto Araki,et al.  Leaf-Shedding Phenology in Lowland Tropical Seasonal Forests of Cambodia as Estimated From NOAA Satellite Images , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[8]  E. Helmer,et al.  A comparison of radiometric normalization methods when filling cloud gaps in Landsat imagery , 2007 .

[9]  Frédéric Achard,et al.  A satellite data set for tropical forest area change assessment , 2011 .

[10]  E. Helmer,et al.  Cloud-Free Satellite Image Mosaics with Regression Trees and Histogram Matching. , 2005 .

[11]  G. Asner,et al.  Cloud cover in Landsat observations of the Brazilian Amazon , 2001 .

[12]  Matthew H. Godfrey,et al.  Regional Management Units for Marine Turtles: A Novel Framework for Prioritizing Conservation and Research across Multiple Scales , 2010, PloS one.

[13]  S. Adler-Golden,et al.  Atmospheric Correction for Short-wave Spectral Imagery Based on MODTRAN 4 , 2000 .

[14]  Hermann Kaufmann,et al.  Comparison of Topographic Correction Methods , 2009, Remote. Sens..

[15]  Philip H. Swain,et al.  Pattern recognition: A basis for remote sensing data analysis , 1972 .

[16]  Chengquan Huang,et al.  Use of a dark object concept and support vector machines to automate forest cover change analysis , 2008 .

[17]  John R. Jensen,et al.  Introductory Digital Image Processing: A Remote Sensing Perspective , 1986 .

[18]  William J. Volchok,et al.  Radiometric scene normalization using pseudoinvariant features , 1988 .

[19]  Gail P. Anderson,et al.  Atmospheric correction for shortwave spectral imagery based on MODTRAN4 , 1999, Optics & Photonics.

[20]  D. Fuller,et al.  Tropical forest monitoring and remote sensing: A new era of transparency in forest governance? , 2006 .

[21]  G. Foody Thematic map comparison: Evaluating the statistical significance of differences in classification accuracy , 2004 .

[22]  P. Chavez Image-Based Atmospheric Corrections - Revisited and Improved , 1996 .

[23]  D. Roy,et al.  Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data , 2008 .

[24]  R. Congalton,et al.  Evaluating seasonal variability as an aid to cover-type mapping from Landsat Thematic Mapper data in the Northeast , 1995 .

[25]  Jeffrey G. Masek,et al.  Large Area Scene Selection Interface (LASSI): Methodology of Selecting Landsat Imagery for the Global Land Survey 2005. , 2009 .

[26]  Ronald J. Hall,et al.  Operational mapping of the land cover of the forested area of Canada with Landsat data: EOSD land cover program , 2003 .

[27]  Frédéric Achard,et al.  An automated approach for segmenting and classifying a large sample of multi-date Landsat imagery for pan-tropical forest monitoring , 2011 .

[28]  Dengsheng Lu,et al.  Assessment of atmospheric correction methods for Landsat TM data applicable to Amazon basin LBA research , 2002 .

[29]  Natsuko Yoshifuji,et al.  A review of evapotranspiration estimates from tropical forests in Thailand and adjacent regions , 2008 .

[30]  Michael A. Wulder,et al.  Landsat continuity: Issues and opportunities for land cover monitoring , 2008 .

[31]  C. Woodcock,et al.  Classification and Change Detection Using Landsat TM Data: When and How to Correct Atmospheric Effects? , 2001 .

[32]  S. Goetz,et al.  Radiometric rectification - Toward a common radiometric response among multidate, multisensor images , 1991 .

[33]  O. Sonnentag,et al.  Climate change, phenology, and phenological control of vegetation feedbacks to the climate system , 2013 .

[34]  Frédéric Achard,et al.  Pre-processing of a sample of multi-scene and multi-date Landsat imagery used to monitor forest cover changes over the tropics , 2011 .

[35]  Nick A. Chappell,et al.  Forest Environments in the Mekong River Basin. , 2007 .

[36]  T. Brandeis,et al.  Mapping tropical dry forest height, foliage height profiles and disturbance type and age with a time series of cloud-cleared Landsat and ALI image mosaics to characterize avian habitat , 2010 .

[37]  W. Cohen,et al.  Modelling forest cover attributes as continuous variables in a regional context with Thematic Mapper data , 2001 .

[38]  Eriko Ito,et al.  Principal Forest Types of Three Regions of Cambodia: Kampong Thom, Kratie, and Mondolkiri , 2007 .