Retrieval of sub-pixel snow cover information in the Himalayan region using medium and coarse resolution remote sensing data

Numerous constrained and unconstrained algorithms have been used to retrieve sub-pixel snow-cover information quantitatively using medium and coarse spatial resolution multispectral images from the Advanced Wide Field Sensor (AWiFS) and Moderate Resolution Imaging Spectrometer (MODIS) sensors over the Himalayan region. Both the methods give slow convergence rates and inaccurate estimation of sub-pixel components analysed using root mean square (RMS) error and image deviation. Multiplicative iterative algorithms such as the Expectation Maximization Maximum Likelihood Method (EMML) and the Image Space Reconstruction Algorithm (ISRA) based on the minimization of least squares and Kullback–Leibler distances have been attempted to compute the endmembers' abundances in unmixing of satellite data. In this paper we discuss the eigenvalues of minimum noise fraction (MNF) transformation bands, data noise removal using MNF transformation and selection of pure endmembers using satellite images. The normalized difference snow index (NDSI) is also estimated using field spectral reflectance results and satellite images in green and shortwave infrared (SWIR) wavelength regions in order to carry out a comparative analysis for its variations with sub-pixel snow cover fractions. The present analysis shows the advantage of iterative over direct (constrained and unconstrained) methods; constraints are easily handled and allow better regularization of the solution for the ill-conditioned cases. Iterative methods are found to be faster compared to those of direct methods and can be used operationally for all resolution data for accurate estimation of sub-pixel snow cover.

[1]  D S Biggs,et al.  Acceleration of iterative image restoration algorithms. , 1997, Applied optics.

[2]  P. Pilewskie,et al.  Pinatubo and pre‐Pinatubo optical‐depth spectra: Mauna Loa measurements, comparisons, inferred particle size distributions, radiative effects, and relationship to lidar data , 1993 .

[3]  Jeff Dozier,et al.  Snow Mapping and Classification from Landsat Thematic Mapper Data , 1987, Annals of Glaciology.

[4]  Yosio Edemir Shimabukuro,et al.  The least-squares mixing models to generate fraction images derived from remote sensing multispectral data , 1991, IEEE Trans. Geosci. Remote. Sens..

[5]  Jeff Dozier,et al.  Snow Reflectance from LANDSAT-4 Thematic Mapper , 1984, IEEE Transactions on Geoscience and Remote Sensing.

[6]  John A. Antoniades,et al.  Use of filter vectors in hyperspectral data analysis , 1995, Optics & Photonics.

[7]  Toshio Koike,et al.  Retrieval of snow reflectance from Landsat data in rugged terrain , 2002, Annals of Glaciology.

[8]  J. Srinivasulu,et al.  Estimation of spectral reflectance of snow from IRS-1D LISS-III sensor over the Himalayan terrain , 2004 .

[9]  Michel Fily,et al.  Snow grain-size determination from Landsat imagery over Terre Adélie, Antarctica , 1993 .

[10]  S. Sandmeier,et al.  Radiometric corrections of topographically induced effects on Landsat TM data in an alpine environment , 1993 .

[11]  Rudolf Richter,et al.  Correction of atmospheric and topographic effects for high-spatial-resolution satellite imagery , 1997, Defense, Security, and Sensing.

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

[13]  D. Civco.  Topographic normalization of landsat thematic mapper digital imagery , 1989 .

[14]  Victor Haertel,et al.  Spectral linear mixing model in low spatial resolution image data , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[15]  V. V. Salomonsona,et al.  Estimating fractional snow cover from MODIS using the normalized difference snow index , 2004 .

[16]  P. Switzer,et al.  A transformation for ordering multispectral data in terms of image quality with implications for noise removal , 1988 .

[17]  Anil V. Kulkarni,et al.  Algorithm to monitor snow cover using AWiFS data of RESOURCESAT‐1 for the Himalayan region , 2006 .

[18]  Jouni Pulliainen,et al.  Improved linear interpolation method for the estimation of snow-covered area from optical data , 2002 .

[19]  C. Duguay,et al.  Estimating Surface Reflectance and Albedo from Landsat-5 Thematic Mapper over Rugged Terrain , 1992 .

[20]  C. Justice,et al.  The topographic effect on spectral response from nadir-pointing sensors , 1980 .

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

[22]  D. Hall,et al.  Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data , 1995 .

[23]  Chein-I Chang,et al.  Hyperspectral image classification and dimensionality reduction: an orthogonal subspace projection approach , 1994, IEEE Trans. Geosci. Remote. Sens..

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

[25]  Andrew P. Roberts,et al.  Geodynamic implications of paleomagnetic data from Tertiary sediments in Sakhalin, Russia (NW Pacific) , 2003 .

[26]  J. L. Barker,et al.  Landsat MSS and TM post-calibration dynamic ranges , 1986 .

[27]  S.,et al.  Radiometric Calibration of Landsat Thematic Mapper Multispectral Images , 2008 .

[28]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[29]  C. Conese.  Topographic normalization of TM scenes through the use of an atmospheric correction method and digital terrain model , 1993 .

[30]  Charles L. Byrne,et al.  RBI-EMML signal separation for imaging techniques , 2002, SPIE Defense + Commercial Sensing.

[31]  Yuichi Shimamura,et al.  Evaluation of a useful method to identify snow‐covered areas under vegetation – comparisons among a newly proposed snow index, normalized difference snow index, and visible reflectance , 2006 .

[32]  N. DiGirolamo,et al.  MODIS snow-cover products , 2002 .

[33]  R. Tateishi,et al.  Spectral and temporal linear mixing model for vegetation classification , 2004 .

[34]  Joseph W. Boardman,et al.  Inversion of high spectral resolution data , 1990, Other Conferences.

[35]  M. Daube-Witherspoon,et al.  An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT , 1986, IEEE Transactions on Medical Imaging.

[36]  Miguel Velez-Reyes,et al.  Iterative algorithms for unmixing of hyperspectral imagery , 2003, SPIE Defense + Commercial Sensing.

[37]  F. D. van der Meer,et al.  Spectral mixture modelling and spectral stratigraphy in carbonate lithofacies mapping , 1996 .

[38]  Yu. S. Popkov,et al.  Muliplicative algorithms for determining positive solutions of nonlinear equations , 1996 .

[39]  M. S. Moran,et al.  Evaluation of simplified procedures for retrieval of land surface reflectance factors from satellite sensor output , 1992 .

[40]  J. Dozier.  Spectral Signature of Alpine Snow Cover from the Landsat Thematic Mapper , 1989 .

[41]  Dan Tarpley,et al.  Mapping and monitoring of the snow cover fraction over North America , 2003 .

[42]  J. Townshend,et al.  A stepwise regression tree for nonlinear approximation: Applications to estimating subpixel land cover , 2003 .

[43]  Rune Solberg,et al.  Subpixel mapping of snow cover in forests by optical remote sensing , 2003 .

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

[45]  Yu Hen Hu,et al.  Optimal linear spectral unmixing , 1999, IEEE Trans. Geosci. Remote. Sens..