Functional Magnetic Resonance Imaging Processing (eBook)

Dr. Xingfeng Li obtained his first degree in automation control, master degree of engineer in power system control, and Ph.D. degree in pattern recognition and machine intelligence in 1996, 2001, and 2004, respectively. Since then, he has been working in various research institutions in different countries on MRI and PET image analysis. He worked as postdoc research fellow from 2004-2009 at McGill University in Canada and INSERM, Paris 6th University in France. From 2009–2013, he has been a research fellow at the University of Ulster, UK. He is currently working on applying nonlinear system identification theory for studying nonlinear dynamic brain system. He conducts extensive research work using fMRI, diffusion weighted imaging, perfusion weighted imaging, structural MRI, and PET methods to investigate human brain system. His research interests include functional medical imaging analysis, numerical analysis, statistical analysis, nonlinear system identification, and optimization algorithms. He has published dozen of papers in the journals NeuroImage, IEEE Transaction on Medical Imaging, and Medical Image Analysis. He is also a member of the editorial board of Journal of Nonlinear Dynamics.

Chapter 1: MRI perfusion weighted imaging analysis 1.1      Perfusion imaging. 21.1.1       Indicator-dilution theory for DSC-MRI 31.1.2       MTT and CBV calculation. 51.1.3       DSC-MRI time series analysis. 71.2      Gamma-variate fitting. 91.2.1       Linear regression method for Gamma-variate fitting. 101.2.2       Nonlinear regression method for Gamma-variate fitting. 111.2.3       Baseline elimination for gamma-variate fitting. 151.2.4       Linear method and nonlinear method for gamma-variate fitting. 181.3      AIF selection. 181.3.1       Robust method for AIF determination. 191.3.2       Deconvolution calculation and residual function estimation. 211.3.3       SVD method for deconvolution. 221.3.4       L2 norm regularization for PWI study. 241.3.5       Piecewise linear method for ridge regression parameter estimation. 251.3.6       CBF, MTT, CBV, arrive time, and T-max maps. 291.4      Dispersion effects in DSC-MRI 321.4.1       Local density random walk for concentration time course. 321.4.2       Convolution method to study disperse effect 331.5      Summary of the PWI algorithm.. 33  Chapter 2: First level fMRI data analysis for activation detection 2.1      fMRI experimental design. 22.1.1       Block design. 22.1.2       Random ER design. 32.1.3       Phase-encoded design. 62.2      fMRI data pre-processing. 92.2.1       fMRI data motion correction. 92.2.2       fMRI time series normalization. 102.3      Activation detection: model free and model based methods. 112.3.1       Model free method: two sample t test for activation detection. 112.3.2       Correlation analysis method. 122.4      Models for hemodynamic response function and drift 122.4.1       HRF models for activation detection. 122.4.2       Drift models for activation detection. 152.5      General linear model (GLM) for activation detect 162.5.1       Generalized linear model (GLM) for activation detection. 162.5.2       Ordinary least square for parameters estimation in GLM.. 172.5.3       FOS to solve the inverse problem.. 182.5.4       Weighted least square estimation. 202.5.5       AR(1) model 202.5.6       AR(q) model 212.6      Hypothesis test and threshold correction. 222.6.1       Hypothesis test for the activation detection. 222.6.2       Bonferroni and FDR/FWE threshold correction. 242.6.3       Number of independent tests. 272.6.4       Permutation/random test 272.7      Summary of algorithm for 1st level fMRI data analysis. 28   Chapter 3: 2nd level fMRI data analysis using mixed model   3.1      Mixed model for fMRI data analysis. 23.1.1       Fixed and random effects in fMRI analysis. 33.1.2       Generalized linear mixed model for fMRI study. 33.1.3       Mixed model and its numerical estimations. 43.2      Numerical analysis for mixed effect models. 53.2.1       Two stage model for 2nd level fMRI analysis. 53.2.2       Maximum likelihood method for variance estimation. 63.2.3       Different runs combination. 63.2.4       Group comparison in the mixed model 83.3      Iterative trust region method for ML estimation. 103.3.1       Levenberg–Marquardt (LM) algorithm.. 103.3.2       LM algorithm implementation. 113.3.3       T and likelihood (LR) tests for the mixed model 123.3.4       Modified EM algorithm for group average. 123.3.5       One simulation example for the numerical processing. 133.3.6       Simulation to combine 2 runs. 173.3.7       Combination of 100 runs. 183.4      Exception trust region algorithm for second level fMRI data analysis. 203.4.1       Average runs within subject 213.4.2       Comparing fMRI response within subject 223.4.3       Compare group of subjects. 253.4.4       Numerical implementation details. 273.4.5       Further numerical improvement: BFGS method. 273.4.6       Potential applications and further development 283.5      Degree of freedom (DF) estimation. 293.5.1       Estimation of DF for T distribution. 293.5.2       ML estimation of mixture of t distributions for mixed model 303.5.3       Hessian matrix calculation for trust region algorithm.. 313.5.4       Trust region and expectation trust region algorithms for df estimation. 323.6      fMRI data analysis future directions. 333.7      Second level fMRI data processing algorithm summary. 33  Chapter 4 : fMRI effective connectivity study 4.1      Nonlinear system identification method for fMRI effective connectivity analysis. 24.1.1       Current methods for fMRI effective connectivity analysis. 24.1.2       Nonlinear system identification theory. 34.1.3       Granger causality (GC) tests. 64.1.4       Directionality indices. 74.1.5       Network structure and regional time series extraction. 74.1.6       Examples to apply NSIM to study effective connectivity. 94.2      Model selections for effective connectivity study. 114.2.1       Nonlinear model for fMRI effective connectivity study. 114.2.2       Model selection for NSIM in effective connectivity study. 124.2.3       AIC and AICc criteria for model selection. 134.2.4       MLARS algorithm for model selection. 134.2.5       Nonlinear interaction terms for the effective connectivity analysis. 164.2.6       Advantages and disadvantages of NSIM.. 164.3      Robust method for second level analysis. 174.3.1       Robust regression and breakdown-point 174.3.2       Least trimmed squares for second level effective connectivity analysis. 184.4      Effective connectivity for resting-state fMRI data. 204.4.1       Resting-state fMRI 204.4.2       Example of applying NSIM to RSN from rfMRI 214.5      Limitations for fMRI effective connectivity in this study. 224.6      Summary of the algorithm for fMRI effective connectivity study. 23  Chapter 5: Diffusion weighted imaging analysis 5.1      Basic principle of diffusion MRI and DTI data analysis. 25.1.1       Physical background of MRI diffusion equation. 25.1.2       Apparent diffusion coefficient (ADC) map and DTI calculation. 35.1.3       Invariant indices for DTI analysis. 55.1.4       High order DTI data analysis. 75.2      Fiber tracking. 85.2.1       Color encoding method to represent fiber 85.2.2       Fiber tracking and 3D representation. 95.3      High angular resolution diffusion imaging (HARDI) analysis. 115.3.1       Q-ball imaging (QBI) 115.3.2       ODF representation. 125.3.3       ODF reconstruction theory. 125.3.4       Spherical harmonics (SH) 135.3.5       Least square method with constraints. 155.3.6       Testing the algorithm on rat data. 165.4      Adaptive Q-ball imaging regularization. 185.4.1       Generalized cross validation (GCV) algorithm for regularization. 185.4.2       Regularization or not regularization?. 195.4.3       GFA and ODF maps from rat data. 215.4.4       GCV method for human QBI ODF regularization. 225.5      Diffusion spectrum imaging. 235.5.1       Difference between QBI and DSI acquisition. 245.5.2       DSI image analysis. 255.5.3       DSI GFA map using fixed  and GCV regularization method. 265.5.4       ODF map for DSI using fixed method and GCV method. 275.6      Summary and future directions. 285.7      Summary of DTI, QBI and DSI image analysis methods. 28   Chapter 6: Voxel based morphometry and its application to Alzheimer’s disease study  6.1      Background for voxel based morphometry analysis. 26.1.1       MRI image segmentation. 26.1.2       MRI image registration. 36.1.3       Statistical methods for VBM analysis. 36.2      Enhanced VBM.. 36.2.1       Histogram match. 46.2.2       Apply to AD study. 56.3      Longitudinal VBM and its application to AD study. 86.3.1       Longitudinal VBM preprocessing steps. 86.3.2       Results of longitudinal VBM for AD study. 96.4      Effective connectivity for longitudinal data analysis. 116.4.1       AR model within subjects for effective connectivity study. 116.4.2       An example from longitudinal AD structural MRI 126.4.3       Advantage and disadvantages of this study. 136.5      Other type of sMRI data analysis. 146.5.1       AD classification. 146.5.2       Structural covariance. 146.6      Summary of (longitudinal) VBM analysis methods. 14  Appendixes  Question answers and hints