Search tips
Search criteria 


Logo of bmcmiBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Medical Imaging
BMC Med Imaging. 2012; 12: 8.
Published online Apr 23, 2012. doi:  10.1186/1471-2342-12-8
PMCID: PMC3402999
Volumetric BOLD fMRI simulation: from neurovascular coupling to multivoxel imaging
Zikuan Chencorresponding author1 and Vince Calhoun1,2
1The Mind Research Network and LBERI, Albuquerque 87106, NM, USA
2Department of Electrical and Computer Engineering, Albuquerque 87131, NM, USA
corresponding authorCorresponding author.
Zikuan Chen: zchen/at/; Vince Calhoun: vcalhoun/at/
Received December 28, 2011; Accepted April 23, 2012.
The blood oxygenation-level dependent (BOLD) functional magnetic resonance imaging (fMRI) modality has been numerically simulated by calculating single voxel signals. However, the observation on single voxel signals cannot provide information regarding the spatial distribution of the signals. Specifically, a single BOLD voxel signal simulation cannot answer the fundamental question: is the magnetic resonance (MR) image a replica of its underling magnetic susceptibility source? In this paper, we address this problem by proposing a multivoxel volumetric BOLD fMRI simulation model and a susceptibility expression formula for linear neurovascular coupling process, that allow us to examine the BOLD fMRI procedure from neurovascular coupling to MR image formation.
Since MRI technology only senses the magnetism property, we represent a linear neurovascular-coupled BOLD state by a magnetic susceptibility expression formula, which accounts for the parameters of cortical vasculature, intravascular blood oxygenation level, and local neuroactivity. Upon the susceptibility expression of a BOLD state, we carry out volumetric BOLD fMRI simulation by calculating the fieldmap (established by susceptibility magnetization) and the complex multivoxel MR image (by intravoxel dephasing). Given the predefined susceptibility source and the calculated complex MR image, we compare the MR magnitude (phase, respectively) image with the predefined susceptibility source (the calculated fieldmap) by spatial correlation.
The spatial correlation between the MR magnitude image and the magnetic susceptibility source is about 0.90 for the settings of TE = 30 ms, B0 = 3 T, voxel size = 100 micron, vessel radius = 3 micron, and blood volume fraction = 2%. Using these parameters value, the spatial correlation between the MR phase image and the susceptibility-induced fieldmap is close to 1.00.
Our simulation results show that the MR magnitude image is not an exact replica of the magnetic susceptibility source (spatial correlation ≈ 0.90), and that the MR phase image conforms closely with the susceptibility-induced fieldmap (spatial correlation ≈ 1.00).
Keywords: Bold fMRI, Neurovascular coupling, Neuroactive blob (NAB), Blood magnetism, Intravoxel dephasing, Voxelization, Magnetic fieldmap, Spatial correlation
Articles from BMC Medical Imaging are provided here courtesy of
BioMed Central