# Motion artifact cancellation in NIR spectroscopy using discrete Kalman filtering

- Meltem Izzetoglu
^{1}Email author, - Prabhakar Chitrapu
^{2}, - Scott Bunce
^{3}and - Banu Onaral
^{1}

**9**:16

https://doi.org/10.1186/1475-925X-9-16

© Izzetoglu et al; licensee BioMed Central Ltd. 2010

**Received: **3 August 2009

**Accepted: **9 March 2010

**Published: **9 March 2010

## Abstract

### Background

As a continuation of our earlier work, we present in this study a Kalman filtering based algorithm for the elimination of motion artifacts present in Near Infrared spectroscopy (NIR) measurements. Functional NIR measurements suffer from head motion especially in real world applications where movement cannot be restricted such as studies involving pilots, children, etc. Since head movement can cause fluctuations unrelated to metabolic changes in the blood due to the cognitive activity, removal of these artifacts from NIR signal is necessary for reliable assessment of cognitive activity in the brain for real life applications.

### Methods

Previously, we had worked on adaptive and Wiener filtering for the cancellation of motion artifacts in NIR studies. Using the same NIR data set we have collected in our previous work where different speed motion artifacts were induced on the NIR measurements we compared the results of the newly proposed Kalman filtering approach with the results of previously studied adaptive and Wiener filtering methods in terms of gains in signal to noise ratio. Here, comparisons are based on paired t-tests where data from eleven subjects are used.

### Results

The preliminary results in this current study revealed that the proposed Kalman filtering method provides better estimates in terms of the gain in signal to noise ratio than the classical adaptive filtering approach without the need for additional sensor measurements and results comparable to Wiener filtering but better suitable for real-time applications.

### Conclusions

This paper presented a novel approach based on Kalman filtering for motion artifact removal in NIR recordings. The proposed approach provides a suitable solution to the motion artifact removal problem in NIR studies by combining the advantages of the existing adaptive and Wiener filtering methods in one algorithm which allows efficient real time application with no requirement on additional sensor measurements.

## Background

Near infrared spectroscopy is an emerging technology which enables the measurement of changes in the concentration of deoxygenated hemoglobin (deoxy-Hb) and oxygenated hemoglobin (oxy-Hb) noninvasively during functional brain activation in humans [1]. The technology allows the design of portable, safe, affordable, non-invasive and negligibly intrusive monitoring systems which makes it suitable for many operations, including the monitoring of ongoing cognitive activity under routine working conditions and in the field [1–3].

Typically, an optical apparatus consists of a light source by which the tissue is irradiated and a light detector that receives the light after it has interacted with the tissue. In NIR spectroscopy, the range of light used is between 700 to 900 nm since biological tissues are relatively transparent to light in this range [1]. This is mainly due to the fact that within this so called "optical window", the absorbance of the main constituents in the human tissue such as water, oxy- and deoxy-hemoglobin is small allowing the light to penetrate the tissue. Among the main absorbers (chromophores) in the tissue, oxy- and deoxy-Hb are strongly linked to tissue oxygenation and metabolism. Fortunately, in the optical window, the absorption spectra of oxy- and deoxy-Hb remain significantly different allowing spectroscopic separation of these compounds using only a few sample wavelengths. In functional brain imaging studies, since the demand and the consumption of these main absorbers in the brain change during cognitive activity, monitoring the change in their concentrations using NIR spectroscopy provides information about brain function [1–3]. In functional NIR applications, two other variables, namely oxygen index and blood volume, are commonly used to extract information about the cognitive activities performed. They are derived from the change in the concentrations of oxy-Hb and deoxy-Hb extracted from NIR measurements using Beer-Lambert Law. (Detailed information on the calculation of oxygen index and blood volume can be found in [1, 3]).

Due to many attractive atributes, NIR is an ideal candidate for monitoring cortical function in the brain while subjects are engaged in various real life or experimental tasks. However, functional NIR measurements suffer from head motion [4], especially in real world applications where movement cannot be restricted such as in studies involving pilots, children, etc. Head movement can cause fluctuations unrelated to metabolic changes in the blood due to the cognitive activity. These artifacts are often due to the loss of contact of NIR detectors with skin resulting in measurements of either the ambient light or the light emitted directly from the NIR sources. Furthermore, head movement can cause the blood to move towards or away from the measurement area causing the amount of oxygen to increase or decrease in the region of interest. Therefore, removal of motion artifacts from NIR signal is necessary for reliable assessment of cognitive activity in the brain, hence critical to its deployment as a brain monitoring technology suitable for real life applications.

In this article, we propose a new solution for the motion artifact removal from the NIR signal based on Kalman filtering. To our knowledge, adaptive filtering and Wiener filtering are the only techniques used to solve this problem [4]. Both techniques have been widely used for noise reduction in many biomedical, communication, speech processing applications [4–8]. An adaptive filter is usually a finite impulse response (FIR) filter which has an adaptation algorithm that monitors the environment with additional sensors and hardware and varies the filter transfer function according to the changing input signal's characteristics [5–8]. Like adaptive filtering, Wiener filter is an optimal filtering method in the mean square sense, however it uses the statistics of the signals involved to estimate the filter coefficients without the need for additional sensor information [5–8]. Wiener filtering in general demands stationary data and may not be applied in real time efficiently.

In our application, Kalman filtering approach overcomes the problem of using additional sensors and extra wiring requirement of the adaptive filtering. Due to its recursive nature it further allows efficient real time implementation even without stationarity requirements on the data. Results obtained by Kalman filtering achieve better signal to noise ratios (SNR) than the adaptive filtering and are comparable in SNR to Wiener filtering. The performance of the Kalman filter technique combined with the additional benefits of efficient implementation without requiring additional sensors makes the proposed approach a suitable solution for the motion artifact removal problem for NIR studies.

## Methods

### Discrete Kalman Filtering

_{ k }and measurement z

_{ k }dynamics in the form of

as well as assumed statistics on system and measurement noise w
_{
k
}and v
_{
k
}, respectively such as being independent of each other, white and with Gaussian distributions w
_{
k
}~ *N(0,* Q *)*, v
_{
k
}~ *N(0,* R). The Kalman filter is in essence a recursive solution to a least-squares problem.

If all the state space representation matrices; the transition matrix A and the output matrix H; are known, the same system can be easily established and the states and the outputs can be estimated if the initial conditions are known. However, since this will be an open loop system, the estimates will not be robust. Thus, in the Kalman filter, the estimated states are obtained by using a form of feedback control where the error term obtained from the original measurements are fed back to the original system model, whose effect is determined by the Kalman gain matrix. Detailed explanation on the theory and implementation of discrete Kalman filter structure can be found in [9–16].

_{ k-1}) forward in time to obtain the a priori estimates for the next time step ( , ) in other words prediction of the next time step estimates. Note that a priori and a posteriori error , e

_{ k }and error covariance estimates , P

_{ k }respectively are defined as:

The measurement update equations are responsible for the feedback control which incorporates a new measurement (z
_{
k
}) into the a priori estimate (
) through the use of optimal Kalman gain matrix (K
_{
k
}) to obtain an improved a posteriori estimate (
) in other words the correction of the a posteriori estimates. The optimal Kalman gain matrix K
_{
k
}is found such as to minimize the a posteriori error covariance P
_{
k
}in the minimum mean squares sense.

The discrete Kalman filter algorithm starts with initial estimates of a posteriori state and error covariance estimates. Once the time update equations are applied to predict a priori state and error covariance estimates of the next time step, the measurement update equations are applied to these a priori values to find their corrected a posteriori estimates at the same time step using the measurement and the optimal Kalman gain values in the feedback structure that minimizes the a posteriori error covariance matrix in the minimum mean squares sense. Then this procedure is recursively applied using the same time and measurement update pair with the newly generated a posteriori estimates in the place of initial estimates until the final time step is reached. This recursive nature makes the Kalman filter very appealing compared to other techniques (i.e. Wiener filter) since it makes practical implementations much more feasible [12].

### NIR Data Collection Protocol

In this paper, we use the data set we have collected in our previous work [4]. The protocol we had generated was composed of three types of 20 seconds of head movement periods, where the subject was asked to move his/her head up and down continuously and 20 seconds of rest periods in between the head movement periods, where subject was asked to stay still by a prompt on a computer screen. The speed of the head movements was kept constant within each of the three types of head movements, however it was gradually increased from one region to another, starting slow, then medium, then fast in order to capture the effects of different speed head movements on the NIR measurements and to test the performance of all three methods during such conditions. This procedure was repeated two times. A total of eleven subjects participated in this study. All participants signed informed consent statements approved by the Human Subjects Institutional Review Board at Drexel University.

### NIR System Used for Data Collection

### NIR Data Processing for the Application of Kalman, Wiener and Adaptive Filtering

where z
_{
k
}was the motion corrupted NIR measurement, x
_{
k
}was the motion free NIR signal and the measurement noise v
_{
k
}was the motion artifact. The variance of the measurement noise, the motion artifact, σ_{v}
^{2} required in order to be able to perform Kalman filtering was obtained from the data regions where there is head movement. The variance of the system noise, σ_{w}
^{2} is estimated using the AR model parameters and the variance of a prototype motion free NIR signal obtained during the resting period. In any real life situation, the prototypes for the noiseless NIR data and different types of motion artifact can be collected before the protocol starts. This way the variance of the motion artifact and the system model parameters can be estimated before the protocol starts.

By using the estimated system and measurement models we applied the Kalman filter to three different speed head movement data to estimate the noise free NIR signal on eleven subjects. The results were tested on NIR's one channel blood volume data for slow, medium and fast speed head movement regions in comparison with Wiener and adaptive filtering results. Note that for Wiener filtering spectral density estimates were derived from separate motion free and one trial motion corrupted data segments for each of the three motion types. The corresponding Wiener filter of each motion type was applied offline to the remaining trial region with motion artifact for noise suppression. For adaptive filtering we obtained the correlated motion data required for the technique to be applied properly using the measurements simultaneously gathered by an accelerometer attached to the forehead with the NIR sensor. This technique provided real-time application with the drawback of using an extra sensor. Detailed explanation of these practical issues in the application of these previously proposed techniques can be found in [4].

## Results and Discussion

_{x}

^{2}is the variance of motionless NIR data,

*x(n)*, and σ

_{e}

^{2}is the variance of the estimation error,

*e(n)*, which is the difference between the motionless NIR data and motion compensated data after filtering, (

*n*), as

*e(n) = x(n)-*(

*n*). The input SNR was calculated as

_{ν}

^{2}is the variance of motion artifact. Then we obtained ΔSNR = SNR

_{e}-SNR

_{i}[4] for the Kalman, adaptive and Wiener filtering results in order to show the improvements in SNRs on the estimates. A sample result of ΔSNRs for the subject whose data are given in Figure 3 are summarized in Table 1.

ΔSNR (in dBs) for adaptive, Wiener and Kalman filtering for slow, medium and fast head movements

Head Movement Speed | ΔSNR (dB) (Adaptive Filter) | ΔSNR (dB) (Wiener Filter) | ΔSNR (dB) (Kalman Filter) |
---|---|---|---|

Slow | 3.3560 | 5.2526 | 8.5055 |

Medium | 4.1722 | 9.0539 | 7.8306 |

Fast | 2.7906 | 5.7574 | 6.6282 |

The statistical analysis results of ΔSNR (in dBs) for slow, medium and fast head movements.

Head Movement Speed | Statistical analysis for ΔSNR (Kalman vs adaptive filter) | Statistical analysis for ΔSNR (Kalman vs Wiener filter) |
---|---|---|

Slow | S. (t = 2.760, p < 0.020) | N.S. (t = -1.953, p < 0.079) |

Medium | S. (t = 2.783, p < 0.019) | N.S. (t = -0.385, p < 0.708) |

Fast | S. (t = 3.820, p < 0.003) | N.S. (t = -1.274, p < 0.231) |

## Conclusions

In this paper we present a novel approach for motion artifact removal from NIR measurements using Kalman filtering. The proposed approach provides a suitable solution to the motion artifact removal problem in NIR studies by combining the advantages of the existing adaptive and Wiener filtering methods in one algorithm. The results of this preliminary study suggest that the proposed algorithm performs better than the adaptive filtering algorithm providing better SNRs while still holding the real time applicability with the further advantage of no additional sensor requirement. Our results also indicate that the proposed algorithm is comparable in SNR to Wiener filtering, without the constraints on the stationarity and with efficient real time application capability.

## Appendix

## Declarations

### Acknowledgements

Authors would like to thank Mr. Ajit Devaraj for his helps in the collection of the data. This work was sponsored in part by funds from the US ARMY Telemedicine and Advanced Technology Research Center (TATRC), Defense Advanced Research Projects Agency (DARPA) Augmented Cognition Program, the Office of Naval Research (ONR) and Homeland Security, under agreement numbers, W81XWH-08-053, N00014-02-1-0524, N00014-01-1-0986 and N00014-04-1-0119.

## Authors’ Affiliations

## References

- Villringer A, Chance B:
**Non-invasive optical spectroscopy and imaging of human brain function.***Trends in Neuroscience*1997,**20:**435–442. 10.1016/S0166-2236(97)01132-6View ArticleGoogle Scholar - Izzetoglu K, Bunce S, Onaral B, Pourrezaei K, Chance B:
**Functional optical brain imaging using NIR during cognitive tasks.***Int J of Human-Comp Interaction. Special Issue on Augmented Cognition*2004,**17**(2):211–227.View ArticleGoogle Scholar - Izzetoglu M, Izzetoglu K, Bunce S, Onaral B, Pourrezaei K:
**Functional Near Infrared Neuroimaging.***IEEE Trans on Neural Systems and Rehabilitation Engineering*2005,**13**(2):153–159. 10.1109/TNSRE.2005.847377View ArticleGoogle Scholar - Izzetoglu M, Devaraj A, Bunce S, Onaral B:
**Motion artifact cancellation in NIRs using Wiener filtering.***IEEE Trans on BME*2005,**52**(5):934–938. 10.1109/TBME.2005.845243View ArticleGoogle Scholar - Widrow B, Glover JR, McCool JM:
**Adaptive noise canceling: Principles and applications.***Proc IEEE*1975,**63:**1692–1716. 10.1109/PROC.1975.10036View ArticleGoogle Scholar - Haykin S:
*Adaptive filter theory.*Prentice Hall, Upper Saddle River; 2001.Google Scholar - Vasegi S:
*Adaptive signal processing and digital noise reduction.*Wiley & Teubner; 1997.Google Scholar - Widrow B, Stearns S:
*Adaptive signal processing.*Prentice Hall International Inc; 1985.Google Scholar - Kalman RE:
**A New Approach to Linear Filtering and Prediction Problems.***Trans of the ASME--Journal of Basic Engineering*1960, 35–45.Google Scholar - Grewal MS, Andrews AP:
*Kalman filtering: Theory and practice.*Prentice Hall; 1993.Google Scholar - Gelb A:
*Applied optimal estimation.*The MIT Press; 1974.Google Scholar - Welch G, Bishop G:
*An introduction to the Kalman Filter.*Univ. of North Carolina at Chapel Hill, Dept. of Computer Science; 2004.Google Scholar - Rauch HE, Tung F, Striebel CT:
**Maximum likelihood estimates of linear dynamic systems.***AIAA Journal*1965,**3**(8):1445–1450. 10.2514/3.3166MathSciNetView ArticleGoogle Scholar - Crassidis JL, Junkins JL:
*Optimal estimation of dynamic systems.*Boca Raton : Chapman & Hall/CRC; 2004.View ArticleGoogle Scholar - Yazici B, Izzetoglu M, Onaral B, Bilgutay N:
**Kalman Filtering for Self-Similar Processes.***Signal Processing*2006,**86**(4):760–775. 10.1016/j.sigpro.2005.06.009MATHView ArticleGoogle Scholar - Crassidis JL, Markley FL:
**Predictive filtering for nonlinear systems.***J of Guidance, Control and Dynamics*1997,**20**(3):566–572. 10.2514/2.4078MATHView ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.