 Research
 Open Access
 Published:
Nonlinear estimation of BOLD signals with the aid of cerebral blood volume imaging
BioMedical Engineering OnLine volume 15, Article number: 22 (2016)
Abstract
Background
The hemodynamic balloon model describes the change in coupling from underlying neural activity to observed blood oxygen level dependent (BOLD) response. It plays an increasing important role in brain research using magnetic resonance imaging (MRI) techniques. However, changes in the BOLD signal are sensitive to the resting blood volume fraction (i.e., \(V_0\)) associated with the regional vasculature. In previous studies the value was arbitrarily set to a physiologically plausible value to circumvent the illposedness of the inverse problem. These approaches fail to explore actual \(V_0\) value and could yield inaccurate model estimation.
Methods
The present study represents the first empiric attempt to derive the actual \(V_0\) from data obtained using cerebral blood volume imaging, with the aim of augmenting the existing estimation schemes. Bimanual finger tapping experiments were performed to determine how \(V_0\) influences the model estimation of BOLD signals within a singleregion and multipleregions (i.e., dynamic causal modeling). In order to show the significance of applying the true \(V_0\), we have presented the different results obtained when using the real \(V_0\) and assumed \(V_0\) in terms of singleregion model estimation and dynamic causal modeling.
Results
The results show that \(V_0\) significantly influences the estimation results within a singleregion and multipleregions. Using the actual \(V_0\) might yield more realistic and physiologically meaningful model estimation results.
Conclusion
Incorporating regional venous information in the analysis of the hemodynamic model can provide more reliable and accurate parameter estimations and model predictions, and improve the inference about brain connectivity based on fMRI data.
Background
Functional magnetic resonance imaging (fMRI) offers a noninvasive technology to examine hemodynamic signals in the cerebrovascular system. The hemodynamic balloon model was introduced in 1998 to reveal the coupling dynamics between neural activity and blood oxygen level dependent (BOLD) responses by Buxton et al. [1]. The balloon model describes the causal mechanisms within a hemodynamic process in a certain region of interest (ROI) during brain activation. BOLD responses can be observed via the dynamic changes in cerebral blood volume (CBV) v, cerebral blood flow f, and vein deoxyhemoglobin (dHb) content q. This model is especially helpful in understand the potential consequences of interactions between physiological mechanisms. Since the inception of this model, there has been growing interest in using it to interpret observed fMRI data. The model can be used to infer biologically meaningful parameters that could be employed to investigate the changes in underlying physiological variables during brain activation [2–5], restrict the activation detection process with classic statistical techniques [6, 7], and deduce similar systems or different driving conditions [8–11].
The primary causes of unreliability in model estimation is that the BOLD fMRI technique is sensitive to changes in the signal from venous blood. The change in the signal intensity of a particular voxel is strongly dependent on what fraction of the voxel the vessel occupies. Moreove, changes in BOLD signal intensities during task activation are related not only to multiple physiological states but also regional vessel occupancy, including capillaries and large veins. Indeed, the evaluation of model structure also indicates that the blood volume fraction (BVF) greatly influences the uncertainty of model output [12]. However, this problem has been ignored in all previous studies. Most studies performed to data have avoided the illconditioning problem simply by employing a physiological plausible value of \(V_0=0.02\) instead of investigating the actual value in a particular ROI [2–5, 7, 13, 14] or throughout the brain [6, 15].
Given the importance of the true BVF, efforts are needed to incorporate actual vascular information of the voxel in the hemodynamic model estimation. Firstly, when a voxel includes only brain tissue, the assumption of \(V_0=0.02\) is reasonable [2, 16]. However, when a voxel is mostly or totally occupied by a vessel or vessels, the value might typically be above 0.6 [17]. Secondly, voxels associated with a larger amount of blood are always more likely to show significant BOLD activation due to the inherent nature of the fMRI technique. In this situation, employing an unrealistic \(V_0\) value might yield an unreliable model that does not reflect the physiological reality. This illustrates the importance of taking into account the actual BVF during the estimation procedure.
Several methods have been applied in attempts to obtain the true BVF. We recently showed that magnetic resonance angiography (MRA) might provide a method for roughly estimating the BVF value [18]. The results inferred that the \(V_0\) value in a voxel consists of two derivative components: (1) a constant tissue blood volume component \(V_s=0.02\), which is the smallvessel blood volume that includes capillaries and small postscapulaes, and (2) a variable large blood vessels component \(V_l\), which is the blood volume of large blood vessels. However, this method has not been used to obtain the actual \(V_0\) directly. Indeed, the regional CBV can be measured by another imaging modality, called the dynamic susceptibility contrast (DSC) materialenhanced gradientecho (GE) MR technique [19]. The present study therefore augmented the true BVF acquired from CBV imaging in order to focuses on the influence of \(V_0\) on hemodynamic model estimation and the importance of using the true BVF in the analysis.
This paper is organized as follows. Firstly, we briefly review the hemodynamic Balloon model which constitutes the fundamental component of hemodynamic model estimation. Secondly, we explain the important influence of \(V_0\) with the adoption of a realistic value obtained from the CBV imaging technique. Lastly, the influence of \(V_0\) on model estimation within a singleregion and multipleregions according to the results of a classic bimanual finger tapping experiment is discussed in terms of the impacts of the actual \(V_0\) on parameter estimates and statespace reconstruction.
Hemodynamic balloon model
The hemodynamic balloon model describes the dynamic interrelationship between the blood flow f (neural activity to changes in flow), the regional blood volume information v (changes in flow to changes in blood volume and venous outflow), and the vein dHb content q (changes in flow, volume and oxygen extraction fraction to changes in dHb). The hemodynamic process can be described as the follows:
where \(\tau _s\) reflects signal decay, \(\tau _f\) is the feedback autoregulation time constant, \(\tau _0\) is the transit time, \(\alpha\) is a stiffness parameter, \(\epsilon\) is the neuronal efficacy, u(t) is the neuronal input, and \(E_0\) represents the resting oxygen extraction fraction. The variables f, v, and q are expressed in normalized form, relative to resting values. The balloon model can account for the hemodynamic responses in sparse, noisy fMRI measurements [12, 15]. However, since the above describing equations contain a secondorder time derivative, we can introduce a new variable \(s=\dot{f}\) to express this hemodynamic system as a set of four firstorder ordinary differential equations. Then the observed response BOLD signal could be expressed as follows:
This equation is appropriate when using an fMRI machine 1.5T magnet. The observed y is normalized relative to the value at rest, and \(V_0\) is the resting BVF [2]. Equations 1 and 2 consist of the architecture of hemodynamic inputandoutput system. The model architecture is depicted in Fig. 1.
The BOLD response is associated with all of these parameters, but, we know that parameter \(V_0\) can not be identified along with other parameters simultaneously, instead only their product. Most previous efforts have imposed a physiologically plausible value of \(V_0= 0.02\) to handle the illconditioned nature of the problem [2–10]. Changes in the BOLD signal are strongly affected by \(V_0\), and so an unrealistic \(V_0\) may lead to unreliable model parameter estimation.
Experiment
Two human subjects participated in this study. The experiment was approved by the Health Sciences Research Ethics Committee of Zhejiang University, and written informed consent was obtained from both subjects. Functional images were acquired on a 1.5T scanner using a standard fMRI echo planar imaging protocol (resolution: \(64\times 64\) matrix; repetition time \(\mathrm{TR}=2\) s). In total, 110 acquisitions were made in a blockdesigned finger tapping experiment, giving 11 20s blocks. The conditions for successive blocks alternated between rest and task performance, starting with rest. Furthermore, the CBV imaging sequence consisted of 30 T\(2^*\)weighted images that were collected with a GE sequence (resolution: \(128\times 128\) matrix; 0.1 mmol/kg GdDTPA administered using a powered injector). In order to achieve a sufficient signaltonoise ratio and complete coverage of the brain, TR was increased to 3.1 s, since a typical value is 1 s. The other sequence parameters remained unchanged.
All CBV images were downsampled to make their spatial resolution identical to that of the fMRI image, and thereby allow voxelbyvoxel curve analysis. Concentration–time curves were created for each voxel [20–23]. The calculated \(V_0\) was then used in an existing data estimation procedure [24]. Figure 2 shows an example of an axial CBV image and the observed S and fitted concentration–time curves from one voxel. Data preprocessing and statistical analysis were performed using the SPM5 program (Wellcome Department of Cognitive Neurology, http://www.fil.ion.ucl.ac.uk/spm). The activation map was obtained by applying ttests between all bimanual motor conditions and resting baselines with a cutoff for statistical significance of \(P < 0.001\).
Results
Impact of BVF on singleregion model estimation
We now compare and evaluate the respective impact of the realistic and assumed BVFs on hemodynamic model estimation within a singleregion. Firstly, we chose the maximally activated voxel in the left primary motor cortex (LPM) on the basis of the analyzed fMRI data from SPM5 as the ROI (Fig. 2) and then defined the cluster based on faces and edges excluding corners in order for this voxel to have six neighbors. We extracted the ultimate time series to be analyzed by averaging over the time series of seven voxels. This procedure allowed the model parameters and statespace functions for each of the two subjects to be estimated. Furthermore, for the sake of simplicity, we assumed that the neural parameter had the same value throughout all trails: \(\epsilon _1=\epsilon _2=\cdots =\epsilon _n\), where n denotes the number of trials (i.e., \(n=5\) here). A control random search algorithm was applied in the parameter estimation procedure [25].
Figures 4 and 5 show the BOLD signal and underlying physiological variables of the two subjects for the real \(V_0\) derived from CBV imaging in the maximally activated voxel. The estimated BOLD signal and state variables for an assumed value of \(V_0=0.02\) are also drawn in Figs. 3 and 4 (as dashed lines). The comparison indicates that the assumed and true \(V_0\) could produce similar BOLD estimates in terms of magnitude and shape, with only a slight distinction in the plateau period. This result is consistent with those of previous studies involving the balloon model. However, we also found a large difference between the assumed and actual \(V_0\) values in terms of the reconstructed physiological states. We can conclude that the intensity of changes in the underlying state variables with the assumed \(V_0\) were double those with the true \(V_0\); that is, underestimating \(V_0\) produced an overestimation of the physiological state variables. Moreover, Figs. 4 and 5 indicate that a larger difference between actual \(V_0\) and hypothetical \(V_0\), resulted in a greater difference between estimated physiological state. This means that attention should be paid to ensuring that a realistic \(V_0\) is used in model estimation. The presence of a larger amount of blood in an activated voxel magnifies the effects induced by neuronal activity, lead to an excessive signal for that voxel and unrealistic activity predictions. Similar BOLD changes in a voxel associated with larger veins will change f, v, and q less than for a voxel with a smaller blood fraction. Most activation detection techniques are only capable of indicating the neural activity from changes in BOLD signal or activity maps, and they do not direct infer whether the underlying physiological variation is closely related to \(V_0\) and actually reflects neural activity. Under this circumstance, the use of an arbitrary value of \(V_0\) will influence the spatial specificity of fMRI signals in statistical testing. However, we can assume that functional activated regions induced by an experimental event rather than large regional amounts of blood and the employment of an unrealistic \(V_0\) are suitable when fMRI signal estimation and activation detection are exclusively needed.
Table 1 indicates that the uncertainty of \(V_0\) induces changes in other parameters, with \(V_0\) exerting a complicated, nonlinear, and inconsistent influence on the entire hemodynamic process. Table 1 also indicates that \(V_0\) has a greater influence on the estimated neuronal efficacy parameter \(\epsilon\) than on the other parameters (\(\epsilon\) is 0.3910 with the true \(V_0\), and 0.9089 with the hypothetical \(V_0\)). A previous study found that the uncertainty of model output was more sensitive to variation of \(\epsilon\) than those of other parameters, except \(V_0\) [12]. The defined \(\epsilon\) represents the efficacy with which neural activity causes an increased BOLD signal. As a consequence, if we could use the true \(V_0\), the estimated \(\epsilon\) could offer a better and more intuitive reflection of the activation level, enhancing the functional specificity of fMRI.
Impact of BVF on dynamic causal models
As for balloon model research, dynamic causal modeling (DCM) has been introduced to explore effective connectivity based on hemodynamic observations [8, 9]. DCM extends the balloon model from a single region to multiple regions by utilizing a multipleinput, multipleoutput system. Singleregion model estimation supposes that the extrinsic experimental input consistently accesses all brain regions and that a certain brain area only receives input in this way (\(\epsilon u\) in Eq. 1), whereas DCM assumes that responses (\(x_i\) in Eq. 3) are elicited by two distinct inputs sources: the extrinsic influence of the sensory input (\(\epsilon u\) in Eq. 3) and the intrinsic influence of the interaction regions (\(a_{ij} x_k\) in Eq. 3). In other words, DCM uses estimated neural activities (internal and external) to evaluate the causal correlation among brain areas. While the uncertain \(V_0\) has an important influence on parameter \(\epsilon\) in the hemodynamic model, it is interesting to know how the \(V_0\) influences DCM. In this study we therefore also investigated the effect of \(V_0\) on DCM.
We constructed the simplest tworegion hierarchical system in order to demonstrate the significant effect of BVF on the DCM system. From the two brain areas that interact with and influence each other, we could measure the observed BOLD signals that each of the two regions produced, the relationship can be expressed as follows:
where \(x_1\) and \(x_2\) are the neuronal dynamics in two regions, \(u_1\) and \(u_2\) represent external inputs to the system, \(a_{11}\) and \(a_{22}\) represent the internal connectivity within a region without input, \(a_{12}\) and \(a_{21}\) encode the fixed interregion connectivity without input, and \(c_{11}\) and \(c_{22}\) embody the extrinsic influences of input on neuronal activity. One can augmented the state vector consisting of the model parameters at two regions by concatenating them into a single higher dimensional state space and the measurement vector was also expanded to include two observations in two areas [8]. In the experiment, we adopted a 0–1 squarewave function as two inputs, and the system output was two time series from two regions, \(x_1\) and \(x_2\). While attempting to determine the dimension of the parameters, a more efficient filtering strategy was used to deal with the model estimation problem [26, 27]. The estimation scheme employed for DCM is formally identical to that reported previously [5, 15]. The results of this analysis are presented in Fig. 5, in which the effective connections are presented as directed black arrows along with coupling parameters calculated with the real \(V_0\) and assumed \(V_0\). In order to construct the model system, we chose two regions in the left primary (LPM) and the right primary motor cortex (RPM) containing the two maxima of the activation map. The output regionspecific time series comprised all adjacent (based on faces and edges but not corners) voxels of each maximum (a total of seven voxels), the location is shown in Fig. 2. The conflicts between the motor preparation were interpreted as inhibitory connections between the LPM and RPM [28, 29]. The fixed connectivity from the RPM to the LPM is actually slightly weaker than that from the LPM to the RPM. This indicates that backward influences (RPM to LPM) are stronger than forward connections (LPM to RPM). Furthermore, the fixed connectivity in the RPM is stronger than that in the LPM, indicating that the right pathway is used more frequently than the pathway on the left side. From Fig.6 we conclude that the two different \(V_0\) have different impacts, with the largest deviation being about \(40\,\%\) for the strength of the visual input to the LPM or RPM.
Discussions
This study focused on the important but long ignored issue of how the resting cerebral BVF (i.e., \(V_0\)) impacts hemodynamic models. Previous studies have used a physiologically plausible value of \(V_0= 0.02\) instead of exploring the actual \(V_0\) in the model estimation procedure. However, the intensity of any hemodynamic signal change is greatly affected by the regional BVF, since the active domains subject to model estimation often overlap with those areas characterized by a large BVF [30]. Under such circumstances, an inaccurate \(V_0\) may give rise to inaccurate estimates of the parameters and the reconstructed physiological state. This study used CBV imaging to augment the true \(V_0\) calculated in the hemodynamic model. In order to show the significance of applying the true \(V_0\), we have presented the different results obtained when using the real \(V_0\) and assumed \(V_0\) in terms of singleregion model estimation and DCM. It was found that using the actual \(V_0\), yielded more realistic and physiologically meaningful model estimation results.
The results obtained in this study indicate that \(V_0\) has a rather complicated impact on estimated model parameters. Despite the BOLD responses being similar when using the assumed and real \(V_0\), there was a huge difference in the estimated parameters and the derived physiological state in the ROI. Because the balloon model describes the causal mechanism of a hemodynamic system, its order is higher than the externally observable system, which results in poorly identifiable model parameters due to the nature of nonlinear optimization and temporally sparse sampling. These model parameters have clear physiological meanings, and they should be justified and interpreted with caution [13, 31]. If the actual \(V_0\) is adopted, \(\epsilon\) can be more reliably observe via fMRI measurements. Therefore, \(V_0\) significantly influences the evaluations of brain connectivity. There have recently been extensive discussions on DCM and Granger causal modeling (GCM), with an emphasis on the connectivity among distributed brain systems [32–34]. In order to obtain a more robust understanding of brain causality, we used a biophysical model to eliminate signal bias in imaging procedure and variations of the hemodynamic response in diverse brain domains. However, an unrealistic \(V_0\) might degraded such efforts.
A potential limitation of the present study is to the extent that \(V_0\) as measured by CBV imaging is affected by the amount of blood associated with BOLD signals. We consider that both CBV imaging and the BOLD contrast have tiny difference in terms of the \(V_0\). The former contains the volume of blood across arteries, capillaries, and veins, whereas the latter is relevant to capillaries and veins [35]. Although the arterial fraction of CBV is much less than the venous BVF [36, 37], CBV imaging also partly removes the effect of overestimates about BVF. This is therefore a suitable method for approximating the value of \(V_0\). In addition, this study concentrated on explaining the influence of BVF on hemodynamic model estimation, and the results demonstrated the importance of taking advantage of actual BVF information in the estimation procedure. The argument about the origin of the two modalities were beyond the scope of this paper.
Conclusion
The present study presented the first empiric attempt to derive the actual \(V_0\) from data obtained using CBV imaging, with the aim of augmenting the existing estimation schemes. The results show that \(V_0\) significantly influences the estimation results within a singleregion model estimation and DCM. Using the actual \(V_0\) can provide more reliable and accurate parameterizations and model predictions, and improve brain connectivity estimation based on fMRI data.
Abbreviations
 fMRI:

functional magnetic resonance imaging
 BOLD:

blood oxygen level dependent
 BVF:

blood volume fraction
 CBV:

cerebral blood volume
 DCM:

dynamic causal modeling
 ROI:

region of interest
 dHb:

deoxyhemoglobin
 DSC:

dynamic susceptibility contrast
 LPM:

left primary motor
 RPM:

right primary motor
References
 1.
Buxton RB, Frank LR. A model for the coupling between cerebral blood flow and oxygen metabolism during neural stimulation. J Cerebral Blood Flow Metab. 1997;17:64–72.
 2.
Friston KJ, Mechelli A, Turner R, Price CJ. Nonlinear responses in fMRI: the balloon model, Volterra kernels, and other hemodynamics. NeuroImage. 2000;12:466–77.
 3.
Riera JJ, Watanabe J, Kazuki I, Naoki M, Aubert E, Ozaki T, Kawashima R. A statespace model of the hemodynamic approach: nonlinear filtering of BOLD signals. NeuroImage. 2004;21:547–67.
 4.
Johnston LA, Duff E, Egan GF. Particle filtering for nonlinear BOLD signal analysis. In: 9th international conference on medical image computing and computer assisted intervention (MICCAI), Copenhagen, Denmark. 2006. p. 292–9.
 5.
Hu ZH, Zhao XH, Liu HF, Shi PC. Nonlinear analysis of the BOLD signal. EURASIP J Adv Signal Process. 2009;2009:1–13.
 6.
Deneux T, Faugeras O. Using nonlinear models in fMRI data analysis: model selection and activation detection. NeuroImage. 2006;32:1669–89.
 7.
Hu ZH, Zhang HY, Wang LW, Song XL, Shi PC. Joint estimation for nonlinear dynamic system from fMRI time series. In: 10th international conference on medical image computing and computer assisted intervention (MICCAI), Brisbane, Australia. 2007. p. 734–41.
 8.
Friston KJ, Harrison L, Penny W. Dynamic causal modelling. NeuroImage. 2003;19:1273–302.
 9.
Stephan KE, Kasper L, Harrison LM, Daunizeau J, Ouden HEM, Breakspear M, Friston KJ. Nonlinear dynamic causal models for fMRI. NeuroImage. 2008;42:649–62.
 10.
Li XF, Marrelec G, Hess RF, Benali H. A nonlinear identification method to study effective connectivity in functional MRI. Med Image Anal. 2010;14:30–8.
 11.
Li XF, Coyle D, Maguire L, McGinnity TM, Benali H. A model selection method for nonlinear system identification based fMRI effective connectivity analysis. IEEE Trans Med Imaging. 2011;30(7):1365–80.
 12.
Hu ZH, Shi PC. Sensitivity analysis for biomedical models. IEEE Trans Med Imaging. 2010;29(11):1870–81.
 13.
Johnston LA, Duff E, Mareels I, Egan GF. Nonlinear estimation of the BOLD signal. NeuroImage. 2008;40:504–14.
 14.
Hettiarachchi IT, Pathirana PN, Brotchie P. A state space based approach in nonlinear hemodynamic response modeling with fMRI data. In: 32nd annual international conference of the IEEE EMBS, Buenos Aires, Argentina. 2010. p. 2391–4.
 15.
Hu ZH, Shi PC. Nonlinear analysis of BOLD signal: biophysical modeling, physiological states, and functional activation. In: 2007 IEEE international conference on image processing (ICIP), San Antonio, Texas, USA. 2007. p. 145–8.
 16.
Jezzard P, Matt PM, Smith SM. Functional MRI: an introduction to methods. New York: Oxford University Press; 2001.
 17.
Lu HZ, Law M, Johnson G, Ge Y, van Zijl PCM, Helpern JA. Novel approach to the measurement of absolute cerebral blood volume using vascularspaceoccupancy magnetic resonance imaging. Magn Reson Med. 2005;54:1403–11.
 18.
Hu ZH, Liu C, Liu PS, Liu HF. Exploiting magnetic resonance angiography imaging improves model estimation of BOLD signal. PLoS One. 2012;7(2):31612.
 19.
Rempp KA, Brix G, Wenz F, Becker CR, Lorenz FGWJ. Quantification of regional cerebral blood flow and volume with dynamic susceptibility contrastenhanced MR imaging. Radiology. 1994;193:637–41.
 20.
Rosen BR, Belliveau JW, Buchbinder BR, McKinstry RC, Porkka LM, Kennedy DN, Neuder MS, Fisel CR, Aronen HJ, Kwong KK, Weisskoff RM, Cohen MS, Brady TJ. Contrast agents and cerebral hemodynamics. Magn Reson Med. 1991;19:285–92.
 21.
Norman D, Axel L, Berninger WH, Edwards MS, Cann CE, Redington RW, Cox L. Dynamic computed tomography of the brain: techniques, data analysis, and applications. Am J Roentgenol. 1981;136(4):1–12.
 22.
Madsen MT. A simplified formulation of the gamma variate function. Phys Med Biol. 1992;37(7):1597–600.
 23.
Chan AA, Nelson SJ. Simplified gammavariate fitting of perfusion curves. In: 2th IEEE international symposium on biomedical imaging (ISBI), Arlington, VA, USA. 2004. p. 1067–70.
 24.
Hu ZH, Peng JL, Kong DX, Chen YM, Zhang HY, Lu MH, Liu HF. A novel statistical optimization strategy for estimating intravascular indicator dynamics using susceptibility contrastenhanced MRI. IEEE Trans Med Imaging (submitted)
 25.
Hu ZH, Ni PY, Liu C, Zhao XH, Liu HF, Shi PC. Quantitative evaluation of activation state in functional brain imaging. Brain Topogr. 2012;25:362–73.
 26.
Julier SJ, Uhlmann JK. Unscented filtering and nonlinear estimation. Proc IEEE. 2004;92(3):401–22.
 27.
Merwe R, Wan EA. The squareroot unscented Kalman filter for state and parameterestimation. In: 2001 IEEE international conference on acoustics, speech and signal processing, Salt Lake City, Utah, USA. 2001. p. 3461–4.
 28.
Immisch I, Waldvogel D, VanGelderen P, Hallett M. The role of the medial wall and its anatomical variations for bimanual antiphase and inphase movements. NeuroImage. 2001;14:674–84.
 29.
Weerd PD, Reinke K, Ryan L, McIsaac T, Perschler P, Schnyer D, Trouard T, Gmitrof A. Cortical mechanisms for acquisition and performance of bimanual motor sequences. NeuroImage. 2003;19:1405–16.
 30.
Kim DS, Duong TQ, Kim SG. Highresolution mapping of isoorientation columns by fMRI. Nat Neurosci. 2000;3:164–9.
 31.
David O, Guillemain I, Saillet S, Reyt S, Deransart C, Segebarth C, Depaulis A. Identifying neural drivers with functional MRI: an electrophysiological validation. PLoS Biol. 2008;6(12):e315.
 32.
Roebroeck A, Formisano E, Goebel R. The identification of interacting networks in the brain using fMRI: model selection, causality and deconvolution. NeuroImage. 2011;58:296–302.
 33.
Lohmann G, Erfurth K, Muller K, Turner R. Critical comments on dynamic causal modelling. NeuroImage. 2011;59(3):2322–9.
 34.
Friston KJ, Li BJ, Daunizeau J, Stephan KE. Network discovery with DCM. NeuroImage. 2011;56(2):1202–21.
 35.
Uǧurbil K, Adriany G, Andersen P, Chen W, Gruetter R, Hu XP, Merkle H, Kim DS, Kim SG, Strupp J, Zhu XH, Ogawa S. Magnetic resonance studies of brain function and neurochemistry. Ann Rev Biomed Eng. 2000;2:633–60.
 36.
Ito H, Kanno I, Lida H, Hatazawa J, Shimosegawa E, Tamura H, Okudera T. Arterial fraction of cerebral blood volume in humans measured by positron emission tomography. Ann Nucl Med. 2001;15(2):111–6.
 37.
An HY, Lin WL. Cerebral oxygen extraction fraction and cerebral venous blood volume measurements using MRI: effects of magnetic field variation. Magn Reson Med. 2002;47:958–66.
Author’s contributions
YZ lead data collection, performed the data analysis and drafted the manuscript. ZLW assisted with data collection, data analysis and the drafting of the manuscript. ZZC performed data collection. QL supported partly this study and drafted the manuscript. ZHH conceived of the study, guided its design and coordination, participated in data collection, performed the statistical analysis and drafted the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank the editor and two anonymous referees for their insightful suggestions and valuable comments, which helped to improve the quality of our presented work. This work is supported in part by the National Basic Research Program of China under Grant 2013CB329501, in part by the National High Technology Research and Development Program of China under Grant 2012AA011600, in part by the National Natural Science Foundation of China under Grant 81271645, in part by the Public Projects of Science Technology Department of Zhejiang Province under Grant 2013C33162, and in part by the Zhejiang Provincial Natural Science Foundation of China under Grant LY12H18004.
Competing interests
The authors declare that they have no competing interests.
Author information
Affiliations
Corresponding author
Additional information
Yan Zhang and Zhenghui Hu contributed equally to this work
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhang, Y., Wang, Z., Cai, Z. et al. Nonlinear estimation of BOLD signals with the aid of cerebral blood volume imaging. BioMed Eng OnLine 15, 22 (2016). https://doi.org/10.1186/s1293801601376
Received:
Accepted:
Published:
Keywords
 Blood volume fraction
 Cerebral blood volume imaging
 Dynamic causal modeling