A convolutional neural network to filter artifacts in spectroscopic MRI

Purpose Proton MRSI is a noninvasive modality capable of generating volumetric maps of in vivo tissue metabolism without the need for ionizing radiation or injected contrast agent. Magnetic resonance spectroscopic imaging has been shown to be a viable imaging modality for studying several neuropathologies. However, a key hurdle in the routine clinical adoption of MRSI is the presence of spectral artifacts that can arise from a number of sources, possibly leading to false information. Methods A deep learning model was developed that was capable of identifying and filtering out poor quality spectra. The core of the model used a tiled convolutional neural network that analyzed frequency‐domain spectra to detect artifacts. Results When compared with a panel of MRS experts, our convolutional neural network achieved high sensitivity and specificity with an area under the curve of 0.95. A visualization scheme was implemented to better understand how the convolutional neural network made its judgement on single‐voxel or multivoxel MRSI, and the convolutional neural network was embedded into a pipeline capable of producing whole‐brain spectroscopic MRI volumes in real time. Conclusion The fully automated method for assessment of spectral quality provides a valuable tool to support clinical MRSI or spectroscopic MRI studies for use in fields such as adaptive radiation therapy planning

In vivo magnetic resonance spectroscopy: basic methodology and clinical applications

The clinical use of in vivo magnetic resonance spectroscopy (MRS) has been limited for a long time, mainly due to its low sensitivity. However, with the advent of clinical MR systems with higher magnetic field strengths such as 3 Tesla, the development of better coils, and the design of optimized radio-frequency pulses, sensitivity has been considerably improved. Therefore, in vivo MRS has become a technique that is routinely used more and more in the clinic. In this review, the basic methodology of in vivo MRS is described—mainly focused on 1H MRS of the brain—with attention to hardware requirements, patient safety, acquisition methods, data post-processing, and quantification. Furthermore, examples of clinical applications of in vivo brain MRS in two interesting fields are described. First, together with a description of the major resonances present in brain MR spectra, several examples are presented of deviations from the normal spectral pattern associated with inborn errors of metabolism. Second, through examples of MR spectra of brain tumors, it is shown that MRS can play an important role in oncology

Across‐vendor standardization of semi‐LASER for single‐voxel MRS at 3T

The semi‐adiabatic localization by adiabatic selective refocusing (sLASER) sequence provides single‐shot full intensity signal with clean localization and minimal chemical shift displacement error and was recommended by the international MRS Consensus Group as the preferred localization sequence at high‐ and ultra‐high fields. Across‐vendor standardization of the sLASER sequence at 3 tesla has been challenging due to the B1 requirements of the adiabatic inversion pulses and maximum B1 limitations on some platforms. The aims of this study were to design a short‐echo sLASER sequence that can be executed within a B1 limit of 15 μT by taking advantage of gradient‐modulated RF pulses, to implement it on three major platforms and to evaluate the between‐vendor reproducibility of its perfomance with phantoms and in vivo. In addition, voxel‐based first and second order B0 shimming and voxel‐based B1 adjustments of RF pulses were implemented on all platforms. Amongst the gradient‐modulated pulses considered (GOIA, FOCI and BASSI), GOIA‐WURST was identified as the optimal refocusing pulse that provides good voxel selection within a maximum B1 of 15 μT based on localization efficiency, contamination error and ripple artifacts of the inversion profile. An sLASER sequence (30 ms echo time) that incorporates VAPOR water suppression and 3D outer volume suppression was implemented with identical parameters (RF pulse type and duration, spoiler gradients and inter‐pulse delays) on GE, Philips and Siemens and generated identical spectra on the GE ‘Braino’ phantom between vendors. High‐quality spectra were consistently obtained in multiple regions (cerebellar white matter, hippocampus, pons, posterior cingulate cortex and putamen) in the human brain across vendors (5 subjects scanned per vendor per region; mean signal‐to‐noise ratio [less than] 33; mean water linewidth between 6.5 Hz to 11.4 Hz). The harmonized sLASER protocol is expected to produce high reproducibility of MRS across sites thereby allowing large multi‐site studies with clinical cohorts

A novel method to measure T-1-relaxation times of macromolecules and quantification of the macromolecular resonances

Purpose: Macromolecular peaks underlying metabolite spectra influence the quantification of metabolites. Therefore, it is important to understand the extent of contribution from macromolecules (MMs) in metabolite quantification. However, to model MMs more accurately in spectral fitting, differences in T1 relaxation times among individual MM peaks must be considered. Characterization of T1 -relaxation times for all individual MM peaks using a single inversion recovery technique is difficult due to eventual contributions from metabolites. On the contrary, a double inversion recovery (DIR) technique provided flexibility to acquire MM spectra spanning a range of longitudinal magnetizations with minimal metabolite influence. Thus, a novel method to determine T1 -relaxation times of individual MM peaks is reported in this work. Methods: Extensive Bloch simulations were performed to determine inversion time combinations for a DIR technique that yielded adequate MM signal with varying longitudinal magnetizations while minimizing metabolite contributions. MM spectra were acquired using DIR-metabolite-cycled semi-LASER sequence. LCModel concentrations were fitted to the DIR signal equation to calculate T1 -relaxation times. Results: T1 -relaxation times of MMs range from 204 to 510 ms and 253 to 564 ms in gray- and white-matter rich voxels respectively at 9.4T. Additionally, concentrations of 13 MM peaks are reported. Conclusion: A novel DIR method is reported in this work to calculate T1 -relaxation times of MMs in the human brain. T1 -relaxation times and relaxation time corrected concentrations of individual MMs are reported in gray- and white-matter rich voxels for the first time at 9.4T

A fully automated method for tissue segmentation and CSF-correction of proton MRSI metabolites corroborates abnormal hippocampal NAA in schizophrenia

In this report, we describe the implementation and application of a fully automated segmentation routine using SPM99 algorithms and MATLAB for clinical Magnetic Resonance Spectroscopic Imaging (MRSI) studies. By segmenting high-resolution 3-D image data and coregistering the results to the spatial localizer slices of a spectroscopy examination, the program offers the possibility to easily calculate segmentation maps for a large variety of MRSI experiments. The segmented data are corrected for the individual point-spread function, slice and VOI profiles for measurement sequences with selective pulses as well as for the chemical shifts of different metabolites. The new method was applied to investigate discrete hippocampal metabolite abnormalities in a small sample of schizophrenic patients in comparison to healthy controls (15 patients, 15 controls). Only after correction was the N-acetyl-aspartate (NAA) signal significantly lower in patients compared to controls. No differences were found for the corrected signals from the creatine/phosphocreatine (Cr) or choline-containing compounds (Ch). These results are in good agreement with neuropathological and previous MR spectroscopy studies of the hippocampus in schizophrenic patients

Numerical simulation of PRESS localized MR spectroscopy

Numerical simulations of NMR spectra can provide a rapid and convenient method for optimizing acquisition sequence parameters and generating prior spectral information required for parametric spectral analysis. For spatially resolved spectroscopy, spatially dependent variables affect the resultant spectral amplitudes and phases, which must therefore be taken into account in any spectral simulation model. In this study, methods for numerical simulation of spectra obtained using the PRESS localization pulse sequence are examined. A comparison is made between three different simulation models that include different levels of detail regarding the spatial distributions of the excitation functions, and spin evolution during application of the pulses. These methods were evaluated for measurement of spectra from J-coupled spin systems that are of interest for in vivo proton spectroscopy and results compared with experimental data. It is demonstrated that for optimized refocusing pulses it is sufficient to account for chemical shift effects only, although there is some advantage to implementing a more general numerical simulation approach that includes information on RF pulse excitation profiles, which provides sufficient accuracy while maintaining moderate computational requirements and flexibility to handle different spin systems