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.