Optimizing brain tissue contrast with EPI: A simulated annealing approach

A new magnetization preparation and image acquisition scheme was developed to obtain high‐resolution brain images with optimal tissue contrast. The pulse sequence was derived from an optimization process using simulated annealing, without prior assumptions with regard to the number of radiofrequency (RF) pulses and flip angles. The resulting scheme combined two inversion pulses with the acquisition of three images with varying contrast. The combination of the three images allowed separation of gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) based on T1, contrast. It also enabled the correction of small errors in the initial T1 estimates in postprocessing. The use of three‐dimensional (3D) sensitivity‐encoded (SENSE) echo‐planar imaging (EPI) for image acquisition made it possible to achieve a 1.153 mm3 isotropic resolution within a scan time of 10 min 21 s. The cortical GM signal‐to‐noise ratio (SNR) in the calculated GM‐only image varied between 30 and 100. The novel technique was evaluated in combination with blood oxygen level‐dependent (BOLD) functional magnetic resonance imaging (fMRI) on human subjects, and provided for excellent coregistration of anatomical and functional data. Magn Reson Med 54:373–385, 2005. Published 2005 Wiley‐Liss, Inc.

[1]  G M Bydder,et al.  MR Imaging: Clinical Use of the Inversion Recovery Sequence , 1985, Journal of computer assisted tomography.

[2]  B. Argall,et al.  Unraveling multisensory integration: patchy organization within human STS multisensory cortex , 2004, Nature Neuroscience.

[3]  Renxin Chu,et al.  Magnetic Resonance in Medicine 51:22–26 (2004) Signal-to-Noise Ratio and Parallel Imaging Performance of a 16-Channel Receive-Only Brain Coil Array at , 2022 .

[4]  J. Mugler,et al.  Optimization of parameter values for complex pulse sequences by simulated annealing: Application to 3D MP‐RAGE imaging of the brain , 1994, Magnetic resonance in medicine.

[5]  Jorge Herbert de Lira,et al.  Two-Dimensional Signal and Image Processing , 1989 .

[6]  Steven G. Louie,et al.  A Monte carlo simulated annealing approach to optimization over continuous variables , 1984 .

[7]  R S Balaban,et al.  MR relaxation times in human brain: measurement at 4 T. , 1996, Radiology.

[8]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[9]  D R Haynor,et al.  Partial volume tissue classification of multichannel magnetic resonance images-a mixel model. , 1991, IEEE transactions on medical imaging.

[10]  S. Holland,et al.  NMR relaxation times in the human brain at 3.0 tesla , 1999, Journal of magnetic resonance imaging : JMRI.

[11]  Keiji Tanaka,et al.  Human Ocular Dominance Columns as Revealed by High-Field Functional Magnetic Resonance Imaging , 2001, Neuron.

[12]  Jeih-San Liow,et al.  Qualitative and Quantitative Evaluation of Six Algorithms for Correcting Intensity Nonuniformity Effects , 2001, NeuroImage.

[13]  A. Mackay,et al.  In vivo measurement of T2 distributions and water contents in normal human brain , 1997, Magnetic resonance in medicine.

[14]  V. Haughton,et al.  T1 and T2 measurements on a 1.5-T commercial MR imager. , 1989, Radiology.

[15]  T W Redpath,et al.  Technical note: use of a double inversion recovery pulse sequence to image selectively grey or white brain matter. , 1994, The British journal of radiology.

[16]  P. Santago,et al.  Quantification of MR brain images by mixture density and partial volume modeling , 1993, IEEE Trans. Medical Imaging.

[17]  Matthew O'Donnell,et al.  Optimization of two-dimensional spatially selective NMR pulses by simulated annealing , 1988 .

[18]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[19]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[20]  J. Mugler,et al.  Three‐dimensional magnetization‐prepared rapid gradient‐echo imaging (3D MP RAGE) , 1990, Magnetic resonance in medicine.

[21]  R. Ordidge,et al.  Frequency offset corrected inversion (FOCI) pulses for use in localized spectroscopy , 1996, Magnetic resonance in medicine.

[22]  N Ohyama,et al.  Parameter estimation for in vivo magnetic resonance spectroscopy (mrs) using simulated annealing , 1990, Magnetic resonance in medicine.

[23]  S G Kim,et al.  Accurate T1 determination from inversion recovery images: Application to human brain at 4 Tesla , 1994, Magnetic resonance in medicine.

[24]  P S Tofts,et al.  Sources of T1 variance in normal human white matter. , 1991, Magnetic resonance imaging.

[25]  Anil K. Jain,et al.  Two-dimensional phase unwrapping using a block least-squares method , 1999, IEEE Trans. Image Process..

[26]  R P Velthuizen,et al.  MRI segmentation: methods and applications. , 1995, Magnetic resonance imaging.

[27]  J. Duyn,et al.  EPI‐BOLD fMRI of human motor cortex at 1.5 T and 3.0 T: Sensitivity dependence on echo time and acquisition bandwidth , 2004, Journal of magnetic resonance imaging : JMRI.

[28]  J. Duyn,et al.  Design of a SENSE‐optimized high‐sensitivity MRI receive coil for brain imaging , 2002, Magnetic resonance in medicine.

[29]  W. Press,et al.  Numerical Recipes in C++: The Art of Scientific Computing (2nd edn)1 Numerical Recipes Example Book (C++) (2nd edn)2 Numerical Recipes Multi-Language Code CD ROM with LINUX or UNIX Single-Screen License Revised Version3 , 2003 .

[30]  R. Turner,et al.  Optimization of 3-D MP-RAGE Sequences for Structural Brain Imaging , 2000, NeuroImage.

[31]  Michael Brady,et al.  Estimating the bias field of MR images , 1997, IEEE Transactions on Medical Imaging.

[32]  Jagath C. Rajapakse,et al.  Statistical approach to segmentation of single-channel cerebral MR images , 1997, IEEE Transactions on Medical Imaging.

[33]  A. Oppenheim,et al.  Nonlinear filtering of multiplied and convolved signals , 1968 .

[34]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.