Mary Kociuba1 and Daniel Rowe1,2
1Mathematics, Statistics, and Computer Science, Marquette University, Milwaukee, WI, United States, 2Biophysics, Medical College of Wisconsin, Milwaukee, WI, United States
Synopsis
Discarding the phase
component of the time-series removes relevant biological information from a
complex-valued signal. Although, commonly implemented retrospective image
correction techniques fail to account for physiological artifacts in both the
magnitude and phase components of the time-series. Using the CSF signal, observed during the
data acquisition, as a complex-valued regressor increases the statistical power
of fMRI analysis, through reducing unwanted physiological variability in the
complex-valued signal of interest. The improved performance of implementing the
complex-valued image correction methods is demonstrated with a comparison of
magnitude-only and complex-valued spatial correlations. Purpose & Background
In
an fMRI time-series, CSF signal fluctuations independent of the task-activated
BOLD signal are inherent within spatial frequency measurements [1]. Retrospective
image correction (RETROICOR) regression methods substantially reduce
undesirable physiological variability in the time-series [2]. Although in
RETROICOR methods, only the magnitude component of the time-series is
considered, despite evidence of task related change observed in the phase
time-series [3]. Discarding the phase or neglecting to account for the unwanted
physiological signal within the phase decreases the statistical computing power
of the model. This study demonstrates the advantage of using a complex–valued
CSF signal, estimated from an axial slice inferior to the slice of interest, as
a complex RETROICOR regressor in complex-valued fMRI data. Implementing complex
RETROICOR regressors reduce nuisance signals and more accurately estimates
spatial correlations. The metric for showing the increased performance is a
comparison of magnitude-only (MO) and complex-valued (CV) spatial correlations.
Methods
Experimental
human fMRI data is collected with a single subject on a 3.0 T scanner from a
bilateral finger tapping experiment, performed for sixteen 22-second periods,
with an echo planar pulse sequence (TR/TE =1000/39 ms, BW = 125 kHz, 4 mm thick
axial slices, matrix = 96×96, no. of slices = 10, FOV = 24.0 cm, flip angle =
45°, repetitions = 720). The k-space
readout was Nyquist ghost corrected, IFT reconstructed, and TOAST dynamic B0
corrected [4]. The complex-valued CSF signal is estimated from
isolating the ventricles of an inferior axial slice, and measuring the mean CV
signal. The magnitude and phase regression coefficients,
β and
γ are computed
with a slice corresponding to the task-activated motor cortices, with a
complex-valued regression model [5]. The magnitude and phase
components originating from the CSF signal are removed from the task-activated
slice of interest from a MO and CV time-series separately, and the data is
spatially smoothed with a Gaussian kernel with a full-width-half-max of 3
voxels. The MO and CV spatial correlations in terms of their temporal
frequencies [6] before and after the CSF regression, are compared for two seed
voxels:
a which is located in the
white matter, and
b which is located in the region superior to the ventricles.
Results & Discussion
The
seed voxels,
a and
b, in Fig. 1A show significant global
correlations in the magnitude. Comparing Fig. 1A to Fig. 1B, the MO and CV
spatial correlations before CSF regression both exhibit global correlations,
with noticeable overall phase distortions within the complex-valued data.
Discarding the phase component prematurely eliminates valuable biological
information, which will be appropriately removed in the complex-valued
regression if the phase component consists of extraneous signal. Fig. 1C shows
improvement of the spatial correlation maps in the magnitude after CSF
regression, compared to Fig. 1A. Comparing Fig. 1C to Fig. 1D, demonstrates the
performance of the complex-valued CSF signal as a RETRIOCOR regressor. The CV
spatial correlations show a more significant reduction in physiological
variability in the data from including the CSF phase component in the
regression. The complex regressor is robust enough to significantly remove
unwanted physiological signal, and improve spatial correlation maps in
comparison of the CV to the MO spatial correlations. Only removing the
physiological noise in the magnitude component of the time-series will result
in CV correlation maps in Fig. 1D to be identical to Fig. 1C.
Conclusion
While
the results show improvement implementing the magnitude-only regression, a
striking decrease in variability is notable with the complex regression. Using
the CSF signal as a complex-valued regressor in retrospective
correction decreases physiological noise and increases the accuracy of spatial
correlations in complex-valued fMRI. Failing to account for physiological
artifacts in the magnitude and phase components of the time-series in fMRI data
results in an overestimation of spatial correlations.
Acknowledgements
No acknowledgement found.References
1. Dagli et al.,
NeuroImage 1998. 2. Glover et al.,
MRM 2000. 3. Menon et al., MRM 2002.
4. Hahn et al., HBM 2012[d1] . 5. Rowe, NeuroImage 2005. 6. Cordes et al., J. of Am.
NeuroRadiology 2000. 7. Cordes et
al., J. of Am. NeuroRadiology 2001.