Optimal control in a magnetization-prepared rapid acquisition gradient-echo sequence.

This article proposes a numerical framework to determine the optimal magnetization preparation in a three-dimensional magnetization-prepared rapid gradient-echo (MP-RAGE) sequence to obtain the best achievable contrast between target tissues based on differences in their relaxation times. The benefit lies in the adaptation of the algorithm of optimal control, GRAdient Ascent Pulse Engineering (GRAPE), to the optimization of magnetization preparation in a cyclic sequence without full recovery between each cycle. This numerical approach optimizes magnetization preparation of an arbitrary number of radio frequency pulses to enhance contrast, taking into account the establishment of a steady state in the longitudinal component of the magnetization. The optimal control preparation offers an optimized mixed T 1 / T 2 contrast in this traditional T 1 -weighted sequence. To show the versatility of the proposed method, numerical and in vitro results are described. Examples of contrasts acquired on brain regions of a healthy volunteer are presented for potential applications at 3 T.

Intraoperative identification of functional brain areas with RGB imaging using statistical parametric mapping: Simulation and clinical studies.

Complementary technique to preoperative fMRI and electrical brain stimulation (EBS) for glioma resection could improve dramatically the surgical procedure and patient care. Intraoperative RGB optical imaging is a technique for localizing functional areas of the human cerebral cortex that can be used during neurosurgical procedures. However, it still lacks robustness to be used with neurosurgical microscopes as a clinical standard. In particular, a robust quantification of biomarkers of brain functionality is needed to assist neurosurgeons. We propose a methodology to evaluate and optimize intraoperative identification of brain functional areas by RGB imaging. This consist in a numerical 3D brain model based on Monte Carlo simulations to evaluate intraoperative optical setups for identifying functional brain areas. We also adapted fMRI Statistical Parametric Mapping technique to identify functional brain areas in RGB videos acquired for 12 patients. Simulation and experimental results were consistent and showed that the intraoperative identification of functional brain areas is possible with RGB imaging using deoxygenated hemoglobin contrast. Optical functional identifications were consistent with those provided by EBS and preoperative fMRI. We also demonstrated that a halogen lighting may be particularity adapted for functional optical imaging. We showed that an RGB camera combined with a quantitative modeling of brain hemodynamics biomarkers can evaluate in a robust way the functional areas during neurosurgery and serve as a tool of choice to complement EBS and fMRI.

Short echo time dual-frequency MR Elastography with Optimal Control RF pulses.

Magnetic Resonance Elastography (MRE) quantifies the mechanical properties of tissues, typically applying motion encoding gradients (MEG). Multifrequency results allow better characterizations of tissues using data usually acquired through sequential monofrequency experiments. High frequencies are difficult to reach due to slew rate limitations and low frequencies induce long TEs, yielding magnitude images with low SNR. We propose a novel strategy to perform simultaneous multifrequency MRE in the absence of MEGs: using RF pulses designed via the Optimal Control (OC) theory. Such pulses control the spatial distribution of the MRI magnetization phase so that the resulting transverse magnetization reproduces the phase pattern of an MRE acquisition. The pulse is applied with a constant gradient during the multifrequency mechanical excitation to simultaneously achieve slice selection and motion encoding. The phase offset sampling strategy can be adapted according to the excitation frequencies to reduce the acquisition time. Phantom experiments were run to compare the classical monofrequency MRE to the OC based dual-frequency MRE method and showed excellent agreement between the reconstructed shear storage modulus G'. Our method could be applied to simultaneously acquire low and high frequency components, which are difficult to encode with the classical MEG MRE strategy.

Constant gradient elastography with optimal control RF pulses.

This article presents a new motion encoding strategy to perform magnetic resonance elastography (MRE). Instead of using standard motion encoding gradients, a tailored RF pulse is designed to simultaneously perform selective excitation and motion encoding in presence of a constant gradient. The RF pulse is designed with a numerical optimal control algorithm, in order to obtain a magnetization phase distribution that depends on the displacement characteristics inside each voxel. As a consequence, no post-excitation encoding gradients are required. This offers numerous advantages, such as reducing eddy current artifacts, and relaxing the constraint on the gradients maximum switch rate. It also allows to perform MRE with ultra-short TE acquisition schemes, which limits T2 decay and optimizes signal-to-noise ratio. The pulse design strategy is developed and analytically analyzed to clarify the encoding mechanism. Finally, simulations, phantom and ex vivo experiments show that phase-to-noise ratios are improved when compared to standard MRE encoding strategies.

Optimal control design of preparation pulses for contrast optimization in MRI.

This work investigates the use of MRI radio-frequency (RF) pulses designed within the framework of optimal control theory for image contrast optimization. The magnetization evolution is modeled with Bloch equations, which defines a dynamic system that can be controlled via the application of the Pontryagin Maximum Principle (PMP). This framework allows the computation of optimal RF pulses that bring the magnetization to a given state to obtain the desired contrast after acquisition. Creating contrast through the optimal manipulation of Bloch equations is a new way of handling contrast in MRI, which can explore the theoretical limits of the system. Simulation experiments carried out on-resonance quantify the contrast improvement when compared to standard T1 or T2 weighting strategies. The use of optimal pulses is also validated for the first time in both in vitro and in vivo experiments on a small-animal 4.7T MR system. Results demonstrate their robustness to static field inhomogeneities as well as the fact that they can be embedded in standard imaging sequences without affecting standard parameters such as slice selection or echo type. In vivo results on rat and mouse brains illustrate the ability of optimal contrast pulses to create non-trivial contrasts on well-studied structures (white matter versus gray matter).

Active control of the spatial MRI phase distribution with optimal control theory.

This paper investigates the use of Optimal Control (OC) theory to design Radio-Frequency (RF) pulses that actively control the spatial distribution of the MRI magnetization phase. The RF pulses are generated through the application of the Pontryagin Maximum Principle and optimized so that the resulting transverse magnetization reproduces various non-trivial and spatial phase patterns. Two different phase patterns are defined and the resulting optimal pulses are tested both numerically with the ODIN MRI simulator and experimentally with an agar gel phantom on a 4.7T small-animal MR scanner. Phase images obtained in simulations and experiments are both consistent with the defined phase patterns. A practical application of phase control with OC-designed pulses is also presented, with the generation of RF pulses adapted for a Magnetic Resonance Elastography experiment. This study demonstrates the possibility to use OC-designed RF pulses to encode information in the magnetization phase and could have applications in MRI sequences using phase images.

Isotropic reconstruction of a 4-D MRI thoracic sequence using super-resolution.

PURPOSE: Four-dimensional (4D) thoracic magnetic resonance imaging (MRI) sequences have been shown to successfully monitor both tumor and lungs anatomy. However, a high temporal resolution is required to avoid motion artifacts, which leads to volumes with poor spatial resolution. This article proposes to reconstruct an isotropic 4D MRI thoracic sequence with minimum modifications to the acquisition protocols. This could be an important step toward the use of 4D MRI for thoracic radiotherapy applications. METHODS: In a postacquisition step, three orthogonal 4D anisotropic acquisitions are combined using super-resolution to reconstruct a series of isotropic volumes. A new phantom that simulates lung tumor motion is developed to evaluate the performance of the algorithm. The proposed framework is also applied to real data of a lung cancer patient. RESULTS: Subjective and objective evaluations show clear resolution enhancement and partial volume effect diminution. The isotropic reconstruction of patient data significantly improves both the visualization and segmentation of thoracic structures. CONCLUSIONS: The results presented here are encouraging and suggest that super-resolution can be regarded as an efficient method to improve the resolution of 4D MRI sequences. It produces an isotropic 4D sequence that would be impossible to acquire in practice. Further investigations will be required to evaluate its reproducibility in various clinical applications.

The effect of sample age and prediction resolution on myocardial infarction risk prediction.

Myocardial infarction (MI) is one of the leading causes of death in many developed countries. Hence, early detection of MI events is critical for effective preventative therapies, potentially reducing avoidable mortality. One approach for early disease prediction is the use of risk prediction models developed using machine learning techniques. One important component of these models is to provide clinicians with the flexibility to customize (e.g., the prediction range) and use the risk prediction model that they deemed most beneficial for their patients. Therefore, in this paper, we develop MI prediction models and investigate the effect of sample age and prediction resolution on the performance of MI risk prediction models. The cardiovascular health study dataset was used in this study. Results indicate that the prediction model developed using SVM algorithm is capable of achieving high sensitivity, specificity, and balanced accuracy of 95.3%, 84.8%, and 90.1%, respectively, over a time span of 6 years. Both sample age and prediction resolution were found not to have a significant impact on the performance of MI risk prediction models developed using subjects aged 65 and above. This implies that risk prediction models developed using different sample age and prediction resolution is a feasible approach. These models can be integrated into a computer aided screening tool which clinicians can use to interpret and predict the MI risk status of the individual patients after performing the necessary clinical assessments (e.g., cognitive function, physical function, electrocardiography, general changes to health/lifestyle, and medications) required by the models. This could offer a means for clinicians to screen the patients at risk of having MI in the near future and prescribe early medical intervention to reduce the risk.

Texture analysis of bone marrow in knee MRI for classification of subjects with bone marrow lesion - data from the Osteoarthritis Initiative.

Visualization of bone marrow lesion (BML) can improve the diagnosis of many bone disorders that are associated with it. A quantitative approach in detecting BML could increase the accuracy and efficiency of diagnosing those bone disorders. In this paper, we investigated the feasibility of using magnetic resonance imaging (MRI)-based texture to (a) identify slices and (b) classify subjects with and without BML. A total of 58 subjects were studied; 29 of them were affected by BML. The ages of subjects ranged from 45 to 74years with a mean age of 59. Texture parameters were calculated for the weight-bearing region of distal femur. The parameters were then analyzed using Mann-Whitney U test and individual feature selection methods to identify potentially discriminantive parameters. Forward feature selection was applied to select features subset for classification. Classification results from eight classifiers were studied. Results show that 98 of the 147 parameters studied are statistically significantly different between the normal and affected marrows: parameters based on co-occurrence matrix are ranked highest in their separability. The classification of subjects achieved an area under the receiver operating characteristic curve (AUC) of 0.914, and the classification of slices achieved an AUC of 0.780. The results show that MRI-texture-based classification can effectively classify subjects/slices with and without BML.