Adaptive denoising of event-related functional magnetic resonance imaging data using spectral subtraction

A new adaptive signal-preserving technique for noise suppression in event-related functional magnetic resonance imaging (fMRI) data is proposed based on spectral subtraction. The proposed technique estimates a parametric model for the power spectrum of random noise from the acquired data based on the characteristics of the Rician statistical model. This model is subsequently used to estimate a noise-suppressed power spectrum for any given pixel time course by simple subtraction of power spectra. The new technique is tested using computer simulations and real data from event-related fMRI experiments. The results show the potential of the new technique in suppressing noise while preserving the other deterministic components in the signal. Moreover, we demonstrate that further analysis using principal component analysis and independent component analysis shows a significant improvement in both convergence and clarity of results when the new technique is used. Given its simple form, the new method does not change the statistical characteristics of the signal or cause correlated noise to be present in the processed signal. This suggests the value of the new technique as a useful preprocessing step for fMRI data analysis.

[1]  S. Sidney Burrus.  Computer-Based Exercises for Signal Processing Using MATLAB , 1993 .

[2]  Guillermo Sapiro,et al.  Anisotropic 2-D and 3-D averaging of fMRI signals , 2001, IEEE Transactions on Medical Imaging.

[3]  Erkki Oja,et al.  Independent Component Analysis , 2001 .

[4]  S E Petersen,et al.  Detection of cortical activation during averaged single trials of a cognitive task using functional magnetic resonance imaging. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[5]  X Hu,et al.  Wavelet transform‐based Wiener filtering of event‐related fMRI data , 2000, Magnetic resonance in medicine.

[6]  Aapo Hyvärinen,et al.  Fast and robust fixed-point algorithms for independent component analysis , 1999, IEEE Trans. Neural Networks.

[7]  G. Lohmann,et al.  Using replicator dynamics for analyzing fMRI data of the human brain , 2002, IEEE Transactions on Medical Imaging.

[8]  B. P. Lathi,et al.  Modern Digital and Analog Communication Systems , 1983 .

[9]  H. Gudbjartsson,et al.  The rician distribution of noisy mri data , 1995, Magnetic resonance in medicine.

[10]  Markus Svensén,et al.  Probabilistic modeling of single-trial fMRI data , 2000, IEEE Transactions on Medical Imaging.

[11]  Ravi S. Menon,et al.  Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[12]  X Hu,et al.  Retrospective estimation and correction of physiological artifacts in fMRI by direct extraction of physiological activity from MR data , 1996, Magnetic resonance in medicine.

[13]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[14]  Iwao Kanno,et al.  Activation detection in functional MRI using subspace modeling and maximum likelihood estimation , 1999, IEEE Transactions on Medical Imaging.

[15]  Saeed Vaseghi,et al.  Advanced Signal Processing and Digital Noise Reduction , 1996 .

[16]  S Makeig,et al.  Analysis of fMRI data by blind separation into independent spatial components , 1998, Human brain mapping.

[17]  T. Sejnowski,et al.  Functionally Independent Components of the Late Positive Event-Related Potential during Visual Spatial Attention , 1999, The Journal of Neuroscience.

[18]  S Makeig,et al.  Spatially independent activity patterns in functional MRI data during the stroop color-naming task. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Terrence J. Sejnowski,et al.  Independent Component Analysis of Simulated ERP Data , 2000 .

[20]  Lars Kai Hansen,et al.  Modeling the hemodynamic response in fMRI using smooth FIR filters , 2000, IEEE Transactions on Medical Imaging.

[21]  Jean-Francois Mangin,et al.  What is the best similarity measure for motion correction in fMRI time series? , 2002, IEEE Transactions on Medical Imaging.

[22]  V D Calhoun,et al.  Independent component analysis of fMRI data in the complex domain , 2002, Magnetic resonance in medicine.

[23]  W. Orrison,et al.  Functional Brain Imaging , 1995 .