http://www.iisc.ernet.in/
Assimilation of Endmember Variability in Spectral Mixture Analysis for Urban Land Cover Extraction
http://wgbis.ces.iisc.ernet.in/energy/
Uttam Kumar1,2,3          S. Kumar Raja4          Chiranjit Mukhopadhyay2           T.V. Ramachandra1,5,6,*
1Energy and Wetlands Research Group, Centre for Ecological Sciences [CES], 2Department of Management Studies, 5Centre for Sustainable Technologies (astra),
6Centre for infrastructure, Sustainable Transportation and Urban Planning [CiSTUP], Indian Institute of Science, Bangalore – 560012, India.
3International Institute of Information Technology (IIITB), Bangalore-560100, India.
4EADS Innovation Works, Airbus Engineering Centre India, Xylem No 4, Mahadevapura Post, Whitefield Road, Bangalore - 560 048, India.
*Corresponding author:
cestvr@ces.iisc.ernet.in

Results

Simulated data

Accuracy assessment of the abundance maps for 3 classes with small, medium and large variability obtained from VECLS algorithm with reference abundance maps is depicted in Table 5, 6 and 7 respectively. Fig.4, 5 and 6 shows the impact of intra-class variability on the accuracy of unmixing. Note that there are no pixels with abundances from 20% up to 70% in the abundance maps. The scatter in the relationship between predicted and actual class composition increases with an increase in the degree of intra-class spectral variation (from Fig.4 to 6 for each class), reducing the accuracy of the sub-pixel class composition estimates derived in Table 5. However, these accuracies are far better when compared to abundance estimation without endmember variability, also evident from Table 5.

The correlation and RMSE between predicted and actual/real/reference proportions with spatially degraded (down sampled) images are shown in Table 6 (at 50% of the original image size i.e. 50 x 50) and Table 7 (at 10% of the original image size i.e. 10 x 10). It is evident from these Tables that VECLS has an improved accuracy over CLS method. That is, the variability in endmembers has been accounted by the proposed algorithm than assuming single endmember per class. By degrading the images, it was observed that the accuracy comes down marginally and there is no drastic reduction. However, one important conclusion that can be drawn is that it is difficult to capture more variability with higher accuracy. Controlled variability and the effect of changing resolution resulted in two important findings: (i) As variability increases, the accuracy of classification decreases. In fact, the accuracy of sub-pixel fraction estimates linearly decreases with intra-class variability (Barducci & Mecocci, 2005; Settle, 2006).This is rational because higher the intra-class variability, more is the likelihood that the actual spectral characteristics of the endmembers in a pixel will deviate from the fixed endmembers. (ii) Down sampling of image pixels (from high to low resolution) affects the overall accuracy marginally with linear unmixing when a) the endmembers have been identified properly at each original image resolution with respect to the ground conditions and,b) the data fit well into the linear model.

Table 5 : Accuracy of prediction derived from low, medium and large variability data compared to absence of endmember variability (single spectra for each endmember) for 100 x 100 dimension images.

Small variability
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
Class Correlation coefficient (r) (p < 2.2e-16) RMSE
1 0.9998 0.7332 0.0029 0.2312
2 0.9999 0.7403 0.0005 0.2077
3 0.9999 0.7406 0.0004 0.2088
Medium variability
1 0.9948 0.7155 0.0386 0.2506
2 0.9989 0.7380 0.0153 0.2104
3 0.9954 0.7321 0.0308 0.2181
Large variability
1 0.9851 0.7897 0.0563 0.2114
2 0.9968 0.7344 0.0259 0.2712
3 0.9850 0.8076 0.0561 0.1399

Table 6 : Accuracy of prediction derived from low, medium and large variability data compared to absence of endmember variability at half the original image dimension.

Small variability
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
Class Correlation coefficient (r) (p < 2.2e-16) RMSE
1 0.9944 0.7105 0.0039 0.1732
2 0.9988 0.7211 0.0015 0.1223
3 0.9988 0.7213 0.0015 0.1221
Medium variability
1 0.9921 0.7341 0.0089 0.1022
2 0.9973 0.7241 0.0011 0.1210
3 0.9921 0.7222 0.0088 0.1289
Large variability
1 0.9821 0.7879 0.0673 0.0779
2 0.9933 0.7339 0.0079 0.1052
3 0.9841 0.8069 0.0601 0.0649

Table 7 : Accuracy of prediction derived from low, medium and large variability data compared to absence of endmember variability at one-tenth of the original image dimension.

Small variability
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
Class Correlation coefficient (r) (p < 2.2e-16) RMSE
1 0.9911 0.7101 0.0029 0.2111
2 0.9985 0.7199 0.0009 0.2101
3 0.9983 0.7201 0.0011 0.2091
Medium variability
1 0.9889 0.7977 0.0005 0.2061
2 0.9890 0.7810 0.0004 0.2405
3 0.9900 0.7001 0.0414 0.2813
Large variability
1 0.9237 0.7297 0.1841 0.2417
2 0.9122 0.7544 0.1977 0.2106
3 0.9114 0.8076 0.1985 0.1917


Fig. 4. The impact of small intra-class variability on the accuracy of unmixing using VECLS (a-c) and using CLS without variability (d-f).


Fig. 5. The impact of medium intra-class variability on the accuracy of unmixing using VECLS (a-c) and using CLS without variability (d-f).


Fig. 6. The impact of large intra-class variability on the accuracy of unmixing using VECLS (a-c) and using CLS without variability (d-f).

Remote sensing data

The IKONOS, Landsat ETM+ and MODIS images were unmixed using VECLS method taking into account covariance due to variability in endmembers and later, subsequently assuming endmembers with zero covariance for comparison (Fig.7). The MODIS unmixed images are not shown in Fig. 7 due to space constraint. The obtained abundances were validated using correlation and RMSE as shown in Tables 8, 9 and 10.

Table 8 : Accuracy of prediction for IKONOS images for 2 and 3 classes with and without endmember variability

IKONOS (4 m)
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
2 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
vegetation 0.9951 0.8122 0.0527 0.1745
non-vegetation 0.9751 0.8222 0.0751 0.1688
 
3 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
Urban 0.9307 0.7998 0.0715 0.2173
vegetation 0.9461 0.7843 0.0602 0.2207
Water 0.9747 0.7678 0.0321 0.0302


Fig. 7. Classified output for 2 classes (vegetation, non-vegetation), 3 classes (urban, vegetation and water) from IKONOS, and 4 classes (urban, vegetation, water and open area) from Landsat ETM+

Fig. 8, 9 and 10 indicates that, in general, abundances obtained from VECLS technique are much better than CLS method. For 2 and 3 classes, with decrease in the spatial resolution, the accuracy also decreases. However, this may not be always true when using Landsat and MODIS datain case of vegetation (Fig. 8) and in case of vegetation, water and open area for 4 classes (Fig. 10). The accuracy depends on the actual number of endmembers present in the image correlated to ground conditions and also the capability of the spectral bands to distinguish those objects.

Table 9 : Accuracy of prediction for Landsat images for 2, 3 and 4 classes with and without endmember variability

Landsat (25 m)
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
2 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
vegetation 0.9307 0.7998 0.0744 0.2109
non-vegetation 0.9461 0.7843 0.0609 0.2346
 
3 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
urban 0.8905 0.7290 0.1165 0.2846
vegetation 0.9104 0.7967 0.0935 0.2195
        water       0.9214 0.7519 0.0827 0.2516
 
4 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
urban 0.8413 0.6714 0.1424 0.3375
vegetation 0.8525 0.6385 0.1589 0.3787
water 0.8825 0.6192 0.1205 0.3926
open area 0.8161 0.6439 0.1749 0.3596

Table 10 : Accuracy of prediction for MODIS images for 2, 3 and 4 classes with and without endmember variability

MODIS (250 m)
  VECLS CLS (0 covariance) VECLS CLS (0 covariance)
2 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
vegetation 0.9273 0.7761 0.0911 0.2389
non-vegetation 0.9137 0.7561 0.0855 0.2522
 
3 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
urban 0.8621 0.7036 0.1490 0.3054
vegetation 0.8882 0.7055 0.1229 0.3003
water 0.8414 0.7165 0.1620 0.2946
 
4 Class Correlation coefficient (r) (p < 2.2e-16) RMSE
urban 0.8231 0.6606 0.1777 0.4470
vegetation 0.8731 0.6804 0.12093 0.3206
water 0.8930 0.6842 0.1142 0.3292
open area 0.7971 0.6492 0. 2105 0.3629

Also, with a very high spatial resolution data such as IKONOS with 4 m spatial resolution, LMM would not be always very useful. Ideally, in this case, we may get abundances as 1 or 0 since there are high chances that the intrinsic scale of objects on the ground may be of the order of 1 IKONOS MS pixel. On the other hand, LMM might work best for low spatial resolution and high spectral resolution (such as MODIS) as evident from Fig.10 for vegetation, water and open area classes. Since class separability depends on the number of spectral bands, so the separation between endmembers of different classes is important for obtaining higher accuracy. As also evident from simulated data analysis, increase in variability may reduce classification accuracy, but degrading the spatial resolution may not always decrease the accuracies of individual classes compared to accuracies obtained from classification of high spatial resolution data, provided a minimum residual error linear model with proper endmember selection is achieved. Urban class was identified with higher accuracy in IKONOS and Landsat bands compared to MODIS bands. Usually, urban areas have more heterogeneity with contrasting features and variability compared to other LC classes. Moreover, spatial extent of the urban objects are typically of high to medium spatial resolution pixel size, so Landsat data would be appropriate for urban land use / land cover types of studies while IKONOS data would be best for micro level urban studies.


Fig. 8. Correlation between actual and predicted proportion for 2 classes with changing spatial resolutions


Fig. 9. Correlation between actual and predicted proportion for 3 classes with different sensor data


Fig. 10. Correlation between actual and predicted proportion for 4 classes with different sensor data.

BACK  «  TOP   »   NEXT
Citation : Uttam Kumar, S. Kumar Raja, Chiranjit Mukhopadhyay and T.V. Ramachandra., 2013, Assimilation of endmember variability in spectral mixture analysis for urban land cover extraction., Advances in Space Research, Volume 52, Issue 11, 1 December 2013, Pages 2015-2033.
* Corresponding Author :
Dr. T.V. Ramachandra
Energy & Wetlands Research Group, Centre for Ecological Sciences, Indian Institute of Science, Bangalore – 560 012, India.
Tel : +91-80-2293 3099/2293 3503-extn 107,      Fax : 91-80-23601428 / 23600085 / 23600683 [CES-TVR]
E-mail : cestvr@ces.iisc.ernet.in, energy@ces.iisc.ernet.in,     Web : http://wgbis.ces.iisc.ernet.in/energy, http://ces.iisc.ernet.in/grass
E-mail    |    Sahyadri    |    ENVIS    |    GRASS    |    Energy    |      CES      |      CST      |    CiSTUP    |      IISc      |    E-mail