Pixel purity index-based algorithms for endmember extraction from hyperspectral imagery
Farzeen Chaudhry,Chao-Cheng Wu,Weimin Liu,Chein-I Chang,Antonio Plaza
2006-01-01
Abstract:Pixel purity index (PPI) has been widely used in hyperspectral image analysis for endmember extraction because of its publicity and availability in the ENvironment for Visualizing Images (ENVI) developed by Research Systems. Unfortunately, its detailed implementation has never been made Correspondence/Reprint request: Dr. Chein-I Chang, Remote Sensing Signal and Image Processing Laboratory, Department of Computer Science and Electrical Engineering, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA. E-mail: cchang@umbc.edu Farzeen Chaudhry et al. 30 available in the literature. This chapter investigates the PPI based on limited published sources and reinterprets the PPI from a statistical signal processing view point. So, the first step is to develop a version of the PPI which can be implemented in any general-purpose scientific package without appealing for the ENVI, so that the users with no access to the ENVI can still run their own PPI algorithm (using for instance the MATLAB software). Since the main idea of the PPI is to generate a sufficiently large number of random vectors, referred to as skewers to cover all possible directions, the final set of PPIgenerated endmembers may vary due to the use of different sets of skewers. In order to resolve such an inconsistent final selection set of endmembers, two algorithms are also developed. One is called Automatic PPI (APPI) which takes advantages of random nature of the skewers to implement the PPI as many runs as possible with different sets of skewers, then finds their common endmembers for all runs. The APPI is terminated automatically when the common endmembers remain unchanged between two consecutive runs, in which case the final common endmembers are desired endmembers. In contrast to the APPI, a second algorithm, to be called Fast Iterative PPI (FIPPI) makes use of an initialization algorithm to produce an appropriate set of initial endmembers so that the final set of endmembers converges to desired endmembers. In order to ease computational complexity as well as to compact the information, PPI-based algorithms generally require either the Principal Components Analysis (PCA) or Maximum Noise Fraction (MNF) transform to perform dimensionality reduction. Unfortunately, the PCA and MNF are second order statistics-based transforms, they do not necessarily preserve information of high order statistics. Therefore, this chapter further investigates the ICA-based PPI-based algorithms where the MNF originally used in the PPI-based algorithms is replaced by the ICA to perform dimensionality reduction. Experimental results demonstrate that the ICA-based PPI generally performs more effectively than the MNF-based PPI in endmember extraction.