Simo Sarkka1, Aapo Nummenmaa1,2,
1Department of Biomedical
Engineering & Computational Science, Aalto University, Espoo, Finland; 2Athinoula
A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital,
Charlestown, MA, United States; 3Harvard-MIT Division of Health
Sciences & Technology, Harvard University, Cambridge, MA, United States; 4Advanced
Magnetic Imaging Centre, Low Temperature Laboratory, Aalto University, Espoo,
Finland; 5Institute of Biomedical Engineering, National Taiwan
University, Taipei, Taiwan
In this work, we propose a statistical model based method for removal and analysis of physiological noise in fast BOLD fMRI acquisition methods. The proposed stochastic state space model allows for accurate dynamic tracking of time-varying physiological signal frequencies and the estimation method is based on the Interacting Multiple Models (IMM) Kalman filter (KF) algorithm, which is widely used in real time target tracking applications. The method forms statistically the best possible separation of the spatiotemporal BOLD and physiological signals into separate components, which allows for further processing of the de-noised BOLD signal or analysis of the spatial characteristics of the physiological signals. The proposed method was applied to a three-slice EPI data and the results indicate that the method is able to accurately separate the cardiac and respiration signals from the BOLD signal.