Iterative blind separation of Gaussian data of unknown order

A method for blind separation of noisy jointly Gaussian multivariate signals X is presented, where X = AP + N, X is an observed set of vectors, A is the mixing matrix, P is the unknown signal matrix, and N is white noise. The objective is to estimate all matrices on the right-hand side when even their dimensions (the system order) are unknown. The algorithms developed are extensions of the Iterative Order and Noise (ION) estimation algorithm [10]. Improvements made within the iterative structure of ION to better estimate the order and noise yield ION'. The addition of a second-order blind identification algorithm (SOBI, [4]) subsequently yields ONA, which fully characterizes a data set by estimating the (O)rder, (N)oise, and mixing matrix (A). Metrics are developed to evaluate the performance of these algorithms, and their applicability is discussed. Optimum algorithm constants for ION' and ONA are derived, and their range of applicability is outlined. The algorithms are evaluated through application to three types of data: (1) simulated Gaussian data which spans the problem space, (2) a set of non-Gaussian factory data with 577 variables, and (3) a hyperspectral image with 224 channels. The ONA algorithm is extended to 2D (spatial) hyperspectral problems by exploiting spatial rather than time correlation. ONA produces a full characterization of the data with high signal-to-noise ratios for most unknown parameters in the Gaussian case, though the jointly Gaussian P is shown to be most difficult to retrieve. In all three cases, ONA reduces the noise in the data, identifies small sets of highly correlated variables, and unmixes latent signals. The spatial ONA identifies surface features in the hyperspectral image and retrieves sources significantly more independent than those retrieved by PCA. Further exploration of the applicability of these algorithms to other types of data and further algorithmic improvement is recommended. Thesis Supervisor: David H. Staelin Title: Professor of Electrical Engineering

[1]  J. Cardoso,et al.  Maximum Likelihood Source SeparationBy the Expectation-Maximization Technique : Deterministic and Stochastic Implementation , 1995 .

[2]  S. Macenka,et al.  Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) , 1988 .

[3]  D. Rowe A Bayesian approach to blind source separation , 2002 .

[4]  R. Okafor Maximum likelihood estimation from incomplete data , 1987 .

[5]  Eric Moulines,et al.  A blind source separation technique using second-order statistics , 1997, IEEE Trans. Signal Process..

[6]  D. H. Staelin,et al.  Iterative signal-order and noise estimation for multivariate data , 2001 .

[7]  Pierre Comon Independent component analysis - a new concept? signal processing , 1994 .

[8]  Robert A. Schowengerdt,et al.  Remote sensing, models, and methods for image processing , 1997 .

[9]  Hagai Attias,et al.  Independent Factor Analysis , 1999, Neural Computation.

[10]  J. Cardoso,et al.  On optimal source separation based on second and fourth order cumulants , 1996, Proceedings of 8th Workshop on Statistical Signal and Array Processing.

[11]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[12]  Terrence J. Sejnowski,et al.  An Information-Maximization Approach to Blind Separation and Blind Deconvolution , 1995, Neural Computation.

[13]  H. Hartley Maximum Likelihood Estimation from Incomplete Data , 1958 .

[14]  Junehee Lee Blind noise estimation and compensation for improved characterization of multivariable processes , 2000 .

[15]  John A. Richards,et al.  Remote Sensing Digital Image Analysis: An Introduction , 1999 .

[16]  W. Michael Conklin,et al.  Multivariate Bayesian Statistics: Models for Source Separation and Signal Unmixing , 2005, Technometrics.

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