next up previous contents
Next: MATLAB Implementation Up: Mathematical Analysis Previous: Triple-Weighted Integration Method

   
Delay and Dispersion Correction

Since the blood samples used to estimate Ca(t) are taken from the subject's arm, the activity data gathered from them must be corrected for delay and dispersion within the body. In order to perform the delay correction, we used the technique of Iida, by performing a least squares fitting of the equation:

 \begin{displaymath}A^{*}(t) = \frac{\int\limits_{T_1}^{T_2} \left( \alpha \left[...
...
e^{-\beta t} \right] + \gamma C_a(t) \right) dt} {T_2 - T_1}
\end{displaymath} (11)

where A*(t) in this case is the average activity over grey matter, $\int_{T_1}^{T_2}dt$ represents integration over frames, $\alpha$, $\beta$, and $\gamma$ are the fitting parameters, and Ca(t) is the blood data. The delay itself enters this equation by generating Ca(t) from:

 \begin{displaymath}C_a(t) = \bar{g}(t+\delta)
\end{displaymath} (12)

The function $\bar{g}$ is the result of performing dispersion correction on the blood sample data. This is represented by the equation:

 \begin{displaymath}\bar{g}(t) = g(t) + \tau \frac{dg}{dt}
\end{displaymath} (13)

which is an implicit deconvolution of the equation:

\begin{displaymath}g(t) = \bar{g}(t) \otimes \left[ \frac{1}{\tau} e^{\frac{-t}{\tau}}
\right]
\end{displaymath} (14)


next up previous contents
Next: MATLAB Implementation Up: Mathematical Analysis Previous: Triple-Weighted Integration Method
Mark Wolforth <wolforth@bic.mni.mcgill.ca>
Greg Ward <greg@bic.mni.mcgill.ca>
Sean Marrett <sean@bic.mni.mcgill.ca>