Research on Extended Kalman Filter in Extracting Cerebral Blood Flow Signals by Near Infrared Spectroscopy
LIU Song-yang1, LIU Guang-da1, LIU Zhuo-ya2, QIU Ji-qing3*, CAI Jing1, ZHU Zhan-peng3, ZHANG Cheng3, QI Yuan3, ZHANG Shang1
1. Jilin University，Changchun 130061, China
2. Zhuhai College of Jilin University, Zhuhai 519041, China
3. The First Hospital of Jilin University，Changchun 130061, China
Abstract：There are two types of hemoglobin in the cerebral blood stream: oxygenated hemoglobin (HbO2) and reduced hemoglobin (HbR). The changes in the concentration of these two hemoglobins in the cerebral blood flow can reflect the neural activity in the brain. Extracting the signals of concentration changes can provide basis and reference for the diagnosis and treatment of related diseases such as epilepsy focus localization and depression. At present, algorithms for extracting cerebral blood flow signals using near-infrared spectroscopy include the EEMD-ICA method principal component analysis (PCA) , independent component analysis (ICA), the coherent averaging method, Adaptive filtering, etc. The above algorithms have their own characteristics and advantages in the extraction of near-infrared brain neural activity signals. However, the above methods all pay attention to various physiological interferences such as respiration and eye movement and ignore measurement interferences that conform to Gaussian distribution during measurements, such as instrument precision and crosstalk in signal transmission. In order to extract signals of changes in the concentration of oxygenated hemoglobin (HbO2) and reduced hemoglobin (HbR) in cerebral blood flow, a functional near infrared spectroscopy (fNIRS) cerebral blood flow parameter acq uisition device is designed in this article. In the device, a light source Diode near-infrared light sources with wavelengths of 750 and 830 nm were selected to collect brain blood flow changes. The extended Klaman Filter (EKF) algorithm was used to establish a corresponding mathematical model of physiological interference and measurement interference. Perform recursive calculation with the minimum principle, and combine the initial state estimation of the system at the next moment with the measured feedback to obtain a state estimate of infinitely close to the real value at that moment.), The change of the optical density signal is converted into a signal of change in oxygenated hemoglobin (HbO2) and reduced hemoglobin (HbR) concentration. The results show that the method proposed in this paper can effectively remove the measurement interference that conforms to the Gaussian distribution. In the Valsava experiment and the visual evoked experiment, the curve of changes in the concentration of oxygenated hemoglobin (HbO2) and reduced hemoglobin (HbR) in the cerebral blood flow can be extracted. Compared with the mainstream EEMD algorithm for extracting brain signals, its RMSE value is increased by 0.96%, and r value is increased by 0.6%, which indicates that the proposed method has certain advantages. The method proposed in this paper provides an effective method for detecting neural activity in related brain diseases.
刘颂阳，刘光达，刘卓娅，邱吉庆，蔡 靖，朱展鹏，张 程，齐 远，张 尚. 扩展的卡尔曼滤波在近红外光谱提取脑血流信号中的研究[J]. 光谱学与光谱分析, 2020, 40(07): 2048-2053.
LIU Song-yang, LIU Guang-da, LIU Zhuo-ya, QIU Ji-qing, CAI Jing, ZHU Zhan-peng, ZHANG Cheng, QI Yuan, ZHANG Shang. Research on Extended Kalman Filter in Extracting Cerebral Blood Flow Signals by Near Infrared Spectroscopy. SPECTROSCOPY AND SPECTRAL ANALYSIS, 2020, 40(07): 2048-2053.
[1] Yi M H, Qiu T, Okubo M, et al. Vibrational Spectroscopy, 2018, 95: 23.
[2] Jobsis F F. Science, 1977, 198(4323): 1264.
[3] ZHA Yu-tong, LIU Guang-da(查雨彤, 刘光达). Spectroscopy and Spectral Analysis(光谱学与光谱分析)，2015，35(10)：2746.
[4] Gosselin R, Rodrigue D, Duchesne C. Computers & Chemical Engineering, 2011, 35(2): 296.
[5] Liu G D, Wei X, Zhang S, et al. Optik, 2018, 160: 168.
[6] Zhan Y H，Brooks D H，Franceschini M A, et al. Journal of Biomedical Optics, 2005, 10(1)：011014.
[7] Virtanen J，Noponen T. Journal of Biomedical Optics, 2009, 14: 054032.
[8] Gratton E, Toronov V, Wolf U, et al. Journal of Biomedical Optics, 2005, 10(1): 011008.
[9] Brunner S, Fomin P, Kargel Ch. Waste Management, 2015, 38: 49.
[10] Obata T, Liu T T, Miller K L, et al. Neuroimage, 2003, 21: 144.
[11] Fritson K J，Josephs O, Rees G. Magnetic Resonance in Medicine 1998, 39: 41.
[12] Delpy D T, Cope M, Zee P V D, et al. Physics in Medicine & Biology, 1988, 33(12): 1433.
[13] Patterson M S, Wilson B C, Wyman D R. Lasers in Medical Science, 1991, 6(2): 155.
[14] Zhang Tao, Shao Changyan, Wang Xinnian. Atmospheric Scattering-Based Multiple Images Fog Removal, 4th International Congress on Image & Signal Processing, 2011, 10.1109/CISP.2011.6100007.
[15] LIU Guo, LÜ Qun-bo, LIU Yang-yang(刘 国，吕群波，刘扬阳). Acta Photonica Sinica(光子学报)，2018，47(2)：210001.
