IGARSS 2011, Vancouver, Canada HYPERSPECTRAL UNMIXING USING A NOVEL CONVERSION MODEL

1. IGARSS 2011, Vancouver, CanadaHYPERSPECTRAL UNMIXING USING A NOVEL CONVERSION MODEL Fereidoun A. Mianji, Member, IEEE, Shuang Zhou, Member, IEEE, Ye Zhang, Member, IEEE Presentation by: Shuang Zhou School of Electronics and Information Technology Harbin Institute of Technology, Harbin, China

2. Hyperspectral Unmixing Using a Novel Conversion Model • Introduction 2. Linear Mixture Model (LMM) and Mean-Based Algorithms for Unmixing 3. Structure of the Proposed Approach 4. Data and Experimental Design 5. Experiments with an Simulated Image made of Real Hyperspectral Data 6. Experiments with Real Hyperspectral Image with Implanted Mixed Pixels 7. Experiments with Real Hyperspectral image with Many Endmembers 8. Discussion and Conclusion on the Proposed Approach

3. 1. Introduction The precise information obtained by unmixing is very helpful in the other applications such as edge sharpening, area estimation, subpixel level target detection, etc. Fig. 1.1. Spectral variability (a) and mixed-pixel interference (b).

4. 1. IntroductionContinue Unmixing generally refers to a process which includes two steps: • an endmember extraction algorithm (EEA) to define the pure (or purest) pixels in the scene as the endmembers; • an endmember abundance quantification algorithm (inversion step) to define the contribution of different endmembers in each pixel.

5. 2. Linear Mixture Model (LMM) and Mean-Based Algorithms for Unmixing There are two main models to deal with an unmixing problem: • the macroscopic mixture model which models a mixed pixel as a linear combination of materials with relative concentrations (LMM), 2) the intimate spectral mixture model which proposes a nonlinear mixing model of materials. Two main constraints which are often imposed in constrained Linear Spectral Mixture Analysis(LSMA) are a) abundance sum-to-one constraint (ASC) b) abundance nonnegativity constraint (ANC).

6. 2. Linear Mixture Model (LMM) and Mean-Based Algorithms for Unmixing Mean-Based Algorithms for Unmixing Fully constrained least squares (FCLS) algorithm can efficiently meet both abundance constraints and is optimal in terms of least squares error (LSE). FCLS assumes that the spectra for all endmembers in the scene are available and unique. Therefore, it is optimized for ideal cases where the information of the pure pixels are a priori known, and where there is no spectral variability in the scene.

7. 2. Linear Mixture Model (LMM) and Mean-Based Algorithms for UnmixingLimitations of Mean-Based Algorithms Two main obstacles for a successful LSMA-based method include: 1)thenecessary a priori knowledge of the pure signatures of materials present in the image scene is usually not available. 2) these signatures rarely represent pure spectral signatures of all materials in the image scene in practical situations due to the endmember spectral variability. The main resulting drawback with the mean-based unmixing algorithms is being limited to the mean spectra of each endmemeber. They ignore the statistics associated with the training samples.

8. 3. Structure of the Proposed Approach The rational behind the proposed method is exploiting the statistical variations of the training samples extracted from the scene rather than being only limited to their mean vectors. The central idea of the proposed method is converting the abundance quantification problem (inversion step) into a classification problem through assigning a set of artificial class labels to different possible abundances (0%-100%) of an endmember in a pixel. To realize the proposed method, an SVM approach is applied. Thus, the proposed method is referred to as unmixing-to-classification conversion model with SVM and is noted as uccm-SVM. The well known one-against-all (OAA) strategy with parallel architecture is adopted in uccm-SVM for the multiendmember unmixing tasks.

9. 3. Structure of the Proposed Approachuccm-SVM

10. p endmembers: p 1 - Designating sample set i training sample as “one” (initial ize with i=1) SVM training set s including 2 - Making synthetic classes extracted pure using remaining p - 1 training sample sets (“rest”) pixel vectors Image Pixel s SVM1: trained for SVM2: trained for SVM101: t rained for … 0% of “one” 1% of “one” 100% of “one” Quantification result for “one” in all pixels (fractional image) Next i (endmember) No Yes Endmember fractions i> p ? rescaling to unity 3. Structure of the Proposed Approachuccm-SVM Fig. 3.1. The layout and process flow of uccm-SVM designed for unmixing with a resolution of 1%.

11. 4. Data and Experimental Design The hyperspectral data used to validate the proposed unmixing method include: • San Diego data for making a simulated image • Indian Pines data as real image with implanted compositions (mixed pixels) • The University of Pavia data set as a multiendmemebr real image.

12. 5. Experiments with a Simulated Image Made of Real Hyperspectral Data Using 30 randomly selected pixels from three classes of San Diego image, sand, concrete, and roof, we simulated a hyperspectral image. We construct the background with sand and concrete and we use roof to implant the targets. A: sand B: concrete C: roof The simulated image is a 1530 pixel vectors (a 51 by 30 hyperspectral image) Fig. 5.1. The simulated image constructed using samples from 3 classes of San Diego.

13. 5. Experiments with a Simulated Image Made of Real Hyperspectral Data Results We expect to see only 3 peaks, with the amplitue&width of 0.1&15, 0.2&10, and 0.3&5 starting in pixel locations in 300, 750, and 1200, respectively. As can be seen, our method presents a much better result. Fig. 5.2 Unmixing result for roof by FCLS. Fig. 5.3. Unmixing result for roof by uccm-SVM.

14. 5. Experiments with a Simulated Image Made of Real Hyperspectral DataResults From the table below it is clear that the proposed method is superior to the FCLS algorithm. Table 5.1 Comparison of uccm_SVM with FCLS in terms of average square error and computational time.

15. 6. Experiments with Real Hyperspectral Image with Implanted Mixed Pixels Regions of Interest covering corn-notill & soybean-notill (ROI1) and wheat & corn-min (ROI2) were selected. Fig. 6.1 Two chosen ROIs on Indian Pine hyperspectral image for experiments.

16. 6. Experiments with Real Hyperspectral Image with Implanted Mixed PixelsImplanting the Border with Mixture of Neighboring Pixels ROI1 is a 21 by 8 pixel area covering two landcovers, corn-notill and soybean-notill. We implanted the pixels of the border line with a varying mixture of the two endmembers using the neighboring pixels. To unmix the implanted pixel vectors, 21 different pixel vectors from the other areas of each landcover class were randomly chosen as the training samples. Fig. 6.2. Implanting a varying mixture of neighboring pixel vectors in the border line between two adjacent landcover classes

17. 6. Experiments with Real Hyperspectral Image with Implanted Mixed PixelsUnmixing the Mixed Pixels Implanted in ROI1 Table 6.1. Comparison of uccm-SVM with FCLS in terms of average square error and computational time for Indian Pine ROI1. Fig. 6.3. (from left to right) True abundance of corn-notill in the border line of ROI1, unmixing results by FCLS, and unmixing result by uccm-SVM.

18. 7. Experiments with Real Hyperspectral image with Many Endmembers 1- Using all the labeled training and test samples, and a maximum likelihood classifier, we classified the image. 2- With a visual and spectral comparison to the result of some other techniques (e.g. SVM), we insured of the high accuracy of our classification result. 3- Then, we downsampled the image by a 3x3 mean filter to obtain an image (203x113) with a priori known compositions of endmembers. 4- We unmixed the downsampled image by FCLS and uccm-SVM with very low, low, and medium number of training samples (from pure samples of the downsampled image). Fig. 7.1. The University of Pavia data set.

19. 7. Experiments with Real Hyperspectral Image with Many EndmembersResults ASE: uccm-SVM performs better than FCLS on majority of single classes and in ASE average over all sizes of training sets. Computationally: uccm-SVM is faster than FCLS for low number of training samples and slower for higher numbers. Table. 7.1. Average square error (ASE) for the obtained fractional images using FCLS and uccm-SVM for unmixing of The University of Pavia data set (downsampled).

20. 8. Discussion and Conclusion on the Proposed Approach A new learning-based unmixing technique through SVM is proposed in this paper. An unmixing-to-classification conversion model which transfers the unmixing task into a classification problem is the main premise of the method. Abundance quantification for each endmember is performed in a one-against-rest framework using a number of artificially made classes. Uccm-SVM considerably outperformed the mean-based technique, FCLS algorithm. The next step would be improving the performance of the proposed method while reducing its computational burden for big sizes of training sets.

21. Thanks For Your Attention Shuang Zhou School of Electronics and Information Technology Harbin Institute of Technology, Harbin, China