Functional MRI and Vision Science

Publications

          

 

    

 

         

 

 

Projects

·       High resolution fMRI

 

 

 

 

·       Real-time fMRI reconstruction and analysis

 

·       Optimizing fMRI contrast

Functional contrast is established by the variable amounts of T2* decay manifested by neural tissue containing  deoxygenated blood. We therefore use a long echo time, TE, to enhance this contrast. If the functional image variations were entirely at low spatial frequencies, best contrast would be obtained for TE = T2*. Because fMRI is often performed using long acquisitions, the T2* contrast variably affects different spatial frequencies as they are acquired. Given the spatial-frequency content of the functional perturbation, it is possible to calculate the optimal TE. Click here for the detailed calculation.

·       Hemodynamic modeling

The response to a brief (e.g., <2 s) interval of neural activity creates a lengthy response when measured using BOLD-contrast fMRI. We have found that this response, the BOLD hemodynamic impulse response (HIR), can be well modeled as the combination of two linear waveforms. One of these corresponds to the flow fluctuations induced by a brief vascular modulation in a hemodynamic network modeled by a 4-element windkessel, a lumped circuit consisting of two resistances, a compliant element, and an inertive element. This initially positive-going response is a damped sinusoid or sinh.

The second component is initially negative going. An early model attributed this component to a second vascular event within the parenchyma that would cause deoxygenated blood to stagnate. This model simply used a separate, unique 4-element windkessel as a model. This 5-parameter model was remarkably successful in describing a very large ensemble of measurements in several subjects; click here for details. A recent, more detailed model postulates that the negative-going component is the consequence of local diffusion from parenchymal capillaries. The corresponding response can be roughly simulated as a gamma function. This form of the model can also describe the same large dataset with similar efficacy. Details will appear soon.