Data-driven MRSI spectral localization via low-rank component analysis

Improved phosphorus MRSI acquisition through compressed sensing acceleration combined with low-rank reconstruction

OBJECTIVES

Phosphorus-31 magnetic resonance spectroscopic imaging (31P-MRSI) is a non-invasive tool for assessing cellular high-energy metabolism in-vivo. However, its acquisition suffers from a low sensitivity, which necessitates large voxel sizes or multiple averages to achieve an acceptable signal-to-noise ratio (SNR), resulting in long scan times.

MATERIALS AND METHODS

To overcome these limitations, we propose an acquisition and reconstruction scheme for FID-MRSI sequences. Specifically, we employed Compressed Sensing (CS) and Low-Rank (LR) with Total Generalized Variation (TGV) regularization in a combined CS-LR framework. Additionally, we used a novel approach to k-space undersampling that utilizes distinct pseudo-random patterns for each average. To evaluate the proposed method's performance, we performed a retrospective analysis on healthy volunteers' brains and ex-vivo perfused kidneys.

RESULTS

The presented method effectively improves the SNR two-to-threefold while preserving spectral and spatial quality even with threefold acceleration. We were able to recover signal attenuation of anatomical information, and the SNR improvement was obtained while maintaining the metabolites peaks linewidth.

CONCLUSIONS

We presented a novel combined CS-LR acceleration and reconstruction method for FID-MRSI sequences, utilizing a unique approach to k-space undersampling. Our proposed method has demonstrated promising results in enhancing the SNR making it applicable for reducing acquisition time.

Joint learning of nonlinear representation and projection for fast constrained MRSI reconstruction

PURPOSE

To develop and evaluate a novel method for computationally efficient reconstruction from noisy MR spectroscopic imaging (MRSI) data.

METHODS

The proposed method features (a) a novel strategy that jointly learns a nonlinear low-dimensional representation of high-dimensional spectroscopic signals and a neural-network-based projector to recover the low-dimensional embeddings from noisy/limited data; (b) a formulation that integrates the forward encoding model, a regularizer exploiting the learned representation, and a complementary spatial constraint; and (c) a highly efficient algorithm enabled by the learned projector within an alternating direction method of multipliers (ADMM) framework, circumventing the computationally expensive network inversion subproblem.

RESULTS

The proposed method has been evaluated using simulations as well as in vivo  1 $$ {}^1 $$ H and  31 $$ {}^{31} $$ P MRSI data, demonstrating improved performance over state-of-the-art methods, with about 6× $$ \times $$ fewer averages needed than standard Fourier reconstruction for similar metabolite estimation variances and up to 100× $$ \times $$ reduction in processing time compared to a prior neural network constrained reconstruction method. Computational and theoretical analyses were performed to offer further insights into the effectiveness of the proposed method.

CONCLUSION

A novel method was developed for fast, high-SNR spatiospectral reconstruction from noisy MRSI data. We expect our method to be useful for enhancing the quality of MRSI or other high-dimensional spatiospectral imaging data.

Noise-reduction techniques for 1H-FID-MRSI at 14.1 T: Monte Carlo validation and in vivo application

Proton magnetic resonance spectroscopic imaging (1H-MRSI) is a powerful tool that enables the multidimensional non-invasive mapping of the neurochemical profile at high resolution over the entire brain. The constant demand for higher spatial resolution in 1H-MRSI has led to increased interest in post-processing-based denoising methods aimed at reducing noise variance. The aim of the present study was to implement two noise-reduction techniques, Marchenko-Pastur principal component analysis (MP-PCA) based denoising and low-rank total generalized variation (LR-TGV) reconstruction, and to test their potential with and impact on preclinical 14.1 T fast in vivo 1H-FID-MRSI datasets. Since there is no known ground truth for in vivo metabolite maps, additional evaluations of the performance of both noise-reduction strategies were conducted using Monte Carlo simulations. Results showed that both denoising techniques increased the apparent signal-to-noise ratio (SNR) while preserving noise properties in each spectrum for both in vivo and Monte Carlo datasets. Relative metabolite concentrations were not significantly altered by either method and brain regional differences were preserved in both synthetic and in vivo datasets. Increased precision of metabolite estimates was observed for the two methods, with inconsistencies noted for lower-concentration metabolites. Our study provided a framework for how to evaluate the performance of MP-PCA and LR-TGV methods for preclinical 1H-FID MRSI data at 14.1 T. While gains in apparent SNR and precision were observed, concentration estimations ought to be treated with care, especially for low-concentration metabolites.

SNR Enhancement for Multi-TE MRSI Using Joint Low-Dimensional Model and Spatial Constraints

We present a novel method to enhance the SNR for multi-TE MR spectroscopic imaging (MRSI) data by integrating learned nonlinear low-dimensional model and spatial constraints. A deep complex convolutional autoen-coder (DCCAE) was developed to learn a nonlinear low-dimensional representation of the high-dimensional multi-TE 1 H spectroscopy signals. The learned model significantly reduces the data dimension thus serving as an effective constraint for noise reduction. A reconstruction formulation was proposed to integrate the spatiospectral encoding model, the learned model, and a spatial constraint for an SNR-enhancing reconstruction from multi-TE data. The proposed method has been evaluated using both numerical simulations and in vivo brain MRSI experiments. The superior denoising performance of the proposed over alternative methods was demonstrated, both qualitatively and quantitatively. In vivo multi-TE data was used to assess the improved metabolite quantification reproducibility and accuracy achieved by the proposed method. We expect the proposed SNR-enhancing reconstruction to enable faster and/or higher-resolution multi-TE 1 H-MRSI of the brain, potentially useful for various clinical applications.

Whole-brain high-resolution metabolite mapping with 3D compressed-sensing SENSE low-rank 1 H FID-MRSI

There is a growing interest in the neuroscience community to map the distribution of brain metabolites in vivo. Magnetic resonance spectroscopic imaging (MRSI) is often limited by either a poor spatial resolution and/or a long acquisition time, which severely restricts its applications for clinical and research purposes. Building on a recently developed technique of acquisition-reconstruction for 2D MRSI, we combined a fast Cartesian ¹ H-FID-MRSI acquisition sequence, compressed-sensing acceleration, and low-rank total-generalized-variation constrained reconstruction to produce 3D high-resolution whole-brain MRSI with a significant acquisition time reduction. We first evaluated the acceleration performance using retrospective undersampling of a fully sampled dataset. Second, a 20 min accelerated MRSI acquisition was performed on three healthy volunteers, resulting in metabolite maps with 5 mm isotropic resolution. The metabolite maps exhibited the detailed neurochemical composition of all brain regions and revealed parts of the underlying brain anatomy. The latter assessment used previous reported knowledge and a atlas-based analysis to show consistency of the concentration contrasts and ratio across all brain regions. These results acquired on a clinical 3 T MRI scanner successfully combined 3D ¹ H-FID-MRSI with a constrained reconstruction to produce detailed mapping of metabolite concentrations at high resolution over the whole brain, with an acquisition time suitable for clinical or research settings.

Achieving high-resolution 1H-MRSI of the human brain with compressed-sensing and low-rank reconstruction at 7 Tesla

Low sensitivity MR techniques such as magnetic resonance spectroscopic imaging (MRSI) greatly benefit from the gain in signal-to-noise provided by ultra-high field MR. High-resolution and whole-slab brain MRSI remains however very challenging due to lengthy acquisition, low signal, lipid contamination and field inhomogeneity. In this study, we propose an acquisition-reconstruction scheme that combines 1H free-induction-decay (FID)-MRSI sequence, short TR acquisition, compressed sensing acceleration and low-rank modeling with total-generalized-variation constraint to achieve metabolite imaging in two and three dimensions at 7 Tesla. The resulting images and volumes reveal highly detailed distributions that are specific to each metabolite and follow the underlying brain anatomy. The MRSI method was validated in a high-resolution phantom containing fine metabolite structures, and in five healthy volunteers. This new application of compressed sensing acceleration paves the way for high-resolution MRSI in clinical setting with acquisition times of 5 min for 2D MRSI at 2.5 mm and of 20 min for 3D MRSI at 3.3 mm isotropic.

Accelerated radial echo-planar spectroscopic imaging using golden angle view-ordering and compressed-sensing reconstruction with total variation regularization

PURPOSE

To implement a novel, accelerated, 2D radial echo-planar spectroscopic imaging (REPSI) sequence using undersampled radial k-space trajectories and compressed-sensing reconstruction, and to compare results with those from an undersampled Cartesian spectroscopic sequence.

METHODS

The REPSI sequence was implemented using golden-angle view-ordering on a 3T MRI scanner. Radial and Cartesian echo-planar spectroscopic imaging (EPSI) data were acquired at six acceleration factors, each with time-equivalent scan durations, and reconstructed using compressed sensing with total variation regularization. Results from prospectively and retrospectively undersampled phantom and in vivo brain data were compared over estimated concentrations and Cramer-Rao lower-bound values, normalized RMS errors of reconstructed metabolite maps, and percent absolute differences between fully sampled and reconstructed spectroscopic images.

RESULTS

The REPSI method with compressed sensing is able to tolerate greater reductions in scan time compared with EPSI. The reconstruction and quantitation metrics (i.e., spectral normalized RMS error maps, metabolite map normalized RMS error values [e.g., for total N-acetyl asparate, REPSI = 9.4% vs EPSI = 16.3%; acceleration factor = 2.5], percent absolute difference maps, and concentration and Cramer-Rao lower-bound estimates) showed that accelerated REPSI can reduce the scan time by a factor of 2.5 while retaining image and quantitation quality.

CONCLUSION

Accelerated MRSI using undersampled radial echo-planar acquisitions provides greater reconstruction accuracy and more reliable quantitation for a range of acceleration factors compared with time-equivalent compressed-sensing reconstructions of undersampled Cartesian EPSI. Compared to the Cartesian approach, radial undersampling with compressed sensing could help reduce 2D spectroscopic imaging acquisition time, and offers a better trade-off between imaging speed and quality.

Spectral Wavelet-feature Analysis and Classification Assisted Denoising for enhancing magnetic resonance spectroscopy

Magnetic resonance spectroscopy (MRS) is capable of revealing important biochemical and metabolic information of tissues noninvasively. However, the low concentrations of metabolites often lead to poor signal-to-noise ratio (SNR) and a long acquisition time. Therefore, the applications of MRS in detection and quantitative measurements of metabolites in vivo remain limited. Reducing or even eliminating noise can improve SNR sufficiently to obtain high quality spectra in addition to increasing the number of signal averaging (NSA) or the field strength, both of which are limited in clinical applications. We present a Spectral Wavelet-feature ANalysis and Classification Assisted Denoising (SWANCAD) approach to differentiate signal and noise peaks in magnetic resonance spectra based on their respective wavelet features, followed by removing the identified noise components to improve SNR. The performance of this new denoising approach was evaluated by measuring and comparing SNRs and quantified metabolite levels of low NSA spectra (e.g. NSA = 8) before and after denoising using the SWANCAD approach or by conventional spectral fitting and denoising methods, such as LCModel and wavelet threshold methods, as well as the high NSA spectra (e.g. NSA = 192) recorded in the same sampling volumes. The results demonstrated that SWANCAD offers a more effective way to detect the signals and improve SNR by removing noise from the noisy spectra collected with low NSA or in the subminute scan time (e.g. NSA = 8 or 16 s). The potential applications of SWANCAD include using low NSA to accelerate MRS acquisition while maintaining adequate spectroscopic information for detection and quantification of the metabolites of interest when a limited time is available for an MRS examination in the clinical setting.

Constrained Magnetic Resonance Spectroscopic Imaging by Learning Nonlinear Low-Dimensional Models

Magnetic resonance spectroscopic imaging (MRSI) is a powerful molecular imaging modality but has very limited speed, resolution, and SNR tradeoffs. Construction of a low-dimensional model to effectively reduce the dimensionality of the imaging problem has recently shown great promise in improving these tradeoffs. This paper presents a new approach to model and reconstruct the spectroscopic signals by learning a nonlinear low-dimensional representation of the general MR spectra. Specifically, we trained a deep neural network to capture the low-dimensional manifold, where the high-dimensional spectroscopic signals reside. A regularization formulation is proposed to effectively integrate the learned model and physics-based data acquisition model for MRSI reconstruction with the capability to incorporate additional spatiospectral constraints. An efficient numerical algorithm was developed to solve the associated optimization problem involving back-propagating the trained network. Simulation and experimental results were obtained to demonstrate the representation power of the learned model and the ability of the proposed formulation in producing SNR-enhancing reconstruction from the practical MRSI data.

Dynamic hyperpolarized 13 C MR spectroscopic imaging using SPICE in mouse kidney at 9.4 T

This study aims to investigate the feasibility of dynamic hyperpolarized ¹³ C MR spectroscopic imaging (MRSI) using the SPectroscopic Imaging by exploiting spatiospectral CorrElation (SPICE) technique and an estimation of the spatially resolved conversion constant rate (k_pl ). An acquisition scheme comprising a single training dataset and several imaging datasets was proposed considering hyperpolarized ¹³ C circumstances. The feasibility and advantage of the scheme were investigated in two parts: (a) consistency of spectral basis over time and (b) accuracy of the estimated k_pl . The simulations and in vivo experiments support accurate k_pl estimation with consistent spectral bases. The proposed method was implemented in an enzyme phantom and via in vivo experiments. In the enzyme phantom experiments, spatially resolved homogeneous k_pl maps were observed. In the in vivo experiments, normal diet (ND) mice and high-fat diet (HFD) mice had k_pl (s^-1 ) values of medullar (ND: 0.0119 ± 0.0022, HFD: 0.0195 ± 0.0005) and cortical (ND: 0.0148 ±0.0023, HFD: 0.0224 ±0.0054) regions which were higher than vascular (ND: 0.0087 ±0.0013, HFD: 0.0132 ±0.0050) regions. In particular, the k_pl value in the medullar region exhibited a significant difference between the two diet groups. In summary, the feasibility of using modified SPICE for dynamic hyperpolarized ¹³ C MRSI was demonstrated via simulations and in vivo experiments. The consistency of spectral bases over time and the accuracy of the estimated k_pl values validate the proposed acquisition scheme, which comprises only a single training dataset. The proposed method improved the spatial resolution of dynamic hyperpolarized ¹³ C MRSI, which could be used for k_pl estimation using high signal-to-noise ratio spectral bases.

Non-sedated functional imaging based on deep synchronization of PROPELLER MRI and NIRS

BACKGROUND AND OBJECTIVE

Periodically rotated overlapping parallel lines with enhanced reconstruction-echo planar imaging (PROPELLER-EPI) is a promising technique for non-sedated functional imaging due to its unique advantage of motion correction. However, its multiple-blades sampling blood-oxygen-level dependent (BOLD) signal leads to low sampling rate and aliasing of higher frequency physiological signal components such as the cardiac pulsation.

METHODS

In this study, we use near infrared spectroscopy (NIRS) synchronized with pulse sequences of PROPELLER-EPI, utilizing the fact that the optical sensing speed is inherently high. NIRS measures changes of oxyhemoglobin and deoxyhemoglobin to identify the transient states of on-BOLD and off-BOLD, and then labels each blade by temporal co-registration. The labeled blades from multiple epochs of a functional experiment are then used for the k-space data combination and subsequent image reconstruction. An eigenfunction model is proposed for temporal co-registration and to quantify the temporal resolution of the hemodynamic response.

RESULT

The experiment of NIRS labeled PROPELLER-EPI was carried out with the optical sampling rate of 10 Hz and the magnetic pulses repetition time of 1000 ms, and the temporal resolution is 20 times better than that of the state-of-the-art sliding-window PROPELLER-EPI. We compared the functional imaging results against the conventional magnetic resonance echo planar imaging-measured activity and achieved an accuracy of 0.9.

CONCLUSIONS

Using the synchronization of NIRS, the proposed imaging scheme provides an effective way to implement PROPELLER-EPI, which features motion free, high SNR, and enhanced spatial-temporal resolution.

Fast high-resolution brain metabolite mapping on a clinical 3T MRI by accelerated 1 H-FID-MRSI and low-rank constrained reconstruction

PURPOSE

Epitomizing the advantages of ultra short echo time and no chemical shift displacement error, high-resolution-free induction decay magnetic resonance spectroscopic imaging (FID-MRSI) sequences have proven to be highly effective in providing unbiased characterizations of metabolite distributions. However, its merits are often overshadowed in high-resolution settings by reduced signal-to-noise ratios resulting from the smaller voxel volumes procured by extensive phase encoding and the related acquisition times.

METHODS

To address these limitations, we here propose an acquisition and reconstruction scheme that offers both implicit dataset denoising and acquisition acceleration. Specifically, a slice selective high-resolution FID-MRSI sequence was implemented. Spectroscopic datasets were processed to remove fat contamination, and then reconstructed using a total generalized variation (TGV) regularized low-rank model. We further measured reconstruction performance for random undersampled data to assess feasibility of a compressed-sensing SENSE acceleration scheme. Performance of the lipid suppression was assessed using an ad hoc phantom, while that of the low-rank TGV reconstruction model was benchmarked using simulated MRSI data. To assess real-world performance, 2D FID-MRSI acquisitions of the brain in healthy volunteers were reconstructed using the proposed framework.

RESULTS

Results from the phantom and simulated data demonstrate that skull lipid contamination is effectively removed and that data reconstruction quality is improved with the low-rank TGV model. Also, we demonstrated that the presented acquisition and reconstruction methods are compatible with a compressed-sensing SENSE acceleration scheme.

CONCLUSIONS

An original reconstruction pipeline for 2D 1 H-FID-MRSI datasets was presented that places high-resolution metabolite mapping on 3T MR scanners within clinically feasible limits.

Compressed sensing for high-resolution nonlipid suppressed 1 H FID MRSI of the human brain at 9.4T

PURPOSE

The aim of this study was to apply compressed sensing to accelerate the acquisition of high resolution metabolite maps of the human brain using a nonlipid suppressed ultra-short TR and TE ¹ H FID MRSI sequence at 9.4T.

METHODS

X-t sparse compressed sensing reconstruction was optimized for nonlipid suppressed ¹ H FID MRSI data. Coil-by-coil x-t sparse reconstruction was compared with SENSE x-t sparse and low rank reconstruction. The effect of matrix size and spatial resolution on the achievable acceleration factor was studied. Finally, in vivo metabolite maps with different acceleration factors of 2, 4, 5, and 10 were acquired and compared.

RESULTS

Coil-by-coil x-t sparse compressed sensing reconstruction was not able to reliably recover the nonlipid suppressed data, rather a combination of parallel and sparse reconstruction was necessary (SENSE x-t sparse). For acceleration factors of up to 5, both the low-rank and the compressed sensing methods were able to reconstruct the data comparably well (root mean squared errors [RMSEs] ≤ 10.5% for Cre). However, the reconstruction time of the low rank algorithm was drastically longer than compressed sensing. Using the optimized compressed sensing reconstruction, acceleration factors of 4 or 5 could be reached for the MRSI data with a matrix size of 64 × 64. For lower spatial resolutions, an acceleration factor of up to R∼4 was successfully achieved.

CONCLUSION

By tailoring the reconstruction scheme to the nonlipid suppressed data through parameter optimization and performance evaluation, we present high resolution (97 µL voxel size) accelerated in vivo metabolite maps of the human brain acquired at 9.4T within scan times of 3 to 3.75 min.

High-resolution dynamic 31 P-MRSI using a low-rank tensor model

PURPOSE

To develop a rapid ³¹ P-MRSI method with high spatiospectral resolution using low-rank tensor-based data acquisition and image reconstruction.

METHODS

The multidimensional image function of ³¹ P-MRSI is represented by a low-rank tensor to capture the spatial-spectral-temporal correlations of data. A hybrid data acquisition scheme is used for sparse sampling, which consists of a set of "training" data with limited k-space coverage to capture the subspace structure of the image function, and a set of sparsely sampled "imaging" data for high-resolution image reconstruction. An explicit subspace pursuit approach is used for image reconstruction, which estimates the bases of the subspace from the "training" data and then reconstructs a high-resolution image function from the "imaging" data.

RESULTS

We have validated the feasibility of the proposed method using phantom and in vivo studies on a 3T whole-body scanner and a 9.4T preclinical scanner. The proposed method produced high-resolution static ³¹ P-MRSI images (i.e., 6.9 × 6.9 × 10 mm³ nominal resolution in a 15-min acquisition at 3T) and high-resolution, high-frame-rate dynamic ³¹ P-MRSI images (i.e., 1.5 × 1.5 × 1.6 mm³ nominal resolution, 30 s/frame at 9.4T).

CONCLUSIONS

Dynamic spatiospectral variations of ³¹ P-MRSI signals can be efficiently represented by a low-rank tensor. Exploiting this mathematical structure for data acquisition and image reconstruction can lead to fast ³¹ P-MRSI with high resolution, frame-rate, and SNR. Magn Reson Med 78:419-428, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Compartmentalized low-rank recovery for high-resolution lipid unsuppressed MRSI

PURPOSE

To introduce a novel algorithm for the recovery of high-resolution magnetic resonance spectroscopic imaging (MRSI) data with minimal lipid leakage artifacts, from dual-density spiral acquisition.

METHODS

The reconstruction of MRSI data from dual-density spiral data is formulated as a compartmental low-rank recovery problem. The MRSI dataset is modeled as the sum of metabolite and lipid signals, each of which is support limited to the brain and extracranial regions, respectively, in addition to being orthogonal to each other. The reconstruction problem is formulated as an optimization problem, which is solved using iterative reweighted nuclear norm minimization.

RESULTS

The comparisons of the scheme against dual-resolution reconstruction algorithm on numerical phantom and in vivo datasets demonstrate the ability of the scheme to provide higher spatial resolution and lower lipid leakage artifacts. The experiments demonstrate the ability of the scheme to recover the metabolite maps, from lipid unsuppressed datasets with echo time (TE) = 55 ms.

CONCLUSION

The proposed reconstruction method and data acquisition strategy provide an efficient way to achieve high-resolution metabolite maps without lipid suppression. This algorithm would be beneficial for fast metabolic mapping and extension to multislice acquisitions. Magn Reson Med 78:1267-1280, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Denoising of MR spectroscopic imaging data using statistical selection of principal components

OBJECTIVES

To evaluate a new denoising method for MR spectroscopic imaging (MRSI) data based on selection of signal-related principal components (SSPCs) from principal components analysis (PCA).

MATERIALS AND METHODS

A PCA-based method was implemented for selection of signal-related PCs and denoising achieved by reconstructing the original data set utilizing only these PCs. Performance was evaluated using simulated MRSI data and two volumetric in vivo MRSIs of human brain, from a normal subject and a patient with a brain tumor, using variable signal-to-noise ratios (SNRs), metabolite peak areas, Cramer-Rao bounds (CRBs) of fitted metabolite peak areas and metabolite linewidth.

RESULTS

In simulated data, SSPC determined the correct number of signal-related PCs. For in vivo studies, the SSPC denoising resulted in improved SNRs and reduced metabolite quantification uncertainty compared to the original data and two other methods for denoising. The method also performed very well in preserving the spectral linewidth and peak areas. However, this method performs better for regions that have larger numbers of similar spectra.

CONCLUSION

The proposed SSPC denoising improved the SNR and metabolite quantification uncertainty in MRSI, with minimal compromise of the spectral information, and can result in increased accuracy.

COMPARTMENTALIZED LOW-RANK REGULARIZATION WITH ORTHOGONALITY CONSTRAINTS FOR HIGH-RESOLUTION MRSI

We introduce a novel compartmental low rank algorithm for high resolution MR spectroscopic imaging. We model the field inhomogeneity compensated MRSI dataset as the sum of a lipid dataset and a metabolite dataset using the spatial compartmental information obtained from water reference data. Both these datasets are modeled as low-rank subspaces, and are assumed to be orthogonal to each other. We formulate the recovery of the dataset from spiral measurements as a low-rank recovery problem. Experiments using numerical phantom and in-vivo data demonstrates the ability of the algorithm to provide improved spatial resolution and nuisance signal free spectra.

Speeding up dynamic spiral chemical shift imaging with incoherent sampling and low-rank matrix completion

PURPOSE

To improve the temporal and spatial resolution of dynamic ¹³ C spiral chemical shift imaging via incoherent sampling and low-rank matrix completion (LRMC).

METHODS

Spiral CSI data were both simulated and acquired in rats, and undersampling was implemented retrospectively and prospectively by pseudorandomly omitting a fraction of the spiral interleaves. Undersampled data were reconstructed with both LRMC and a conventional inverse nonuniform fast Fourier transform (iNUFFT) and compared with fully sampled data.

RESULTS

Two-fold undersampling with LRMC reconstruction enabled a two-fold improvement in temporal or spatial resolution without significant artifacts or spatiotemporal distortion. Conversely, undersampling with iNUFFT reconstruction created strong artifacts that obscured the image. LRMC performed better at time points with strong metabolite signal.

CONCLUSION

Incoherent undersampling and LRMC provides a way to increase the spatiotemporal resolution of spiral CSI without degrading data integrity. Magn Reson Med 77:951-960, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Magnetic resonance spectroscopic imaging at superresolution: Overview and perspectives

The notion of non-invasive, high-resolution spatial mapping of metabolite concentrations has long enticed the medical community. While magnetic resonance spectroscopic imaging (MRSI) is capable of achieving the requisite spatio-spectral localization, it has traditionally been encumbered by significant resolution constraints that have thus far undermined its clinical utility. To surpass these obstacles, research efforts have primarily focused on hardware enhancements or the development of accelerated acquisition strategies to improve the experimental sensitivity per unit time. Concomitantly, a number of innovative reconstruction techniques have emerged as alternatives to the standard inverse discrete Fourier transform (DFT). While perhaps lesser known, these latter methods strive to effect commensurate resolution gains by exploiting known properties of the underlying MRSI signal in concert with advanced image and signal processing techniques. This review article aims to aggregate and provide an overview of the past few decades of so-called "superresolution" MRSI reconstruction methodologies, and to introduce readers to current state-of-the-art approaches. A number of perspectives are then offered as to the future of high-resolution MRSI, with a particular focus on translation into clinical settings.

High-resolution (1) H-MRSI of the brain using SPICE: Data acquisition and image reconstruction

PURPOSE

To develop data acquisition and image reconstruction methods to enable high-resolution (1) H MR spectroscopic imaging (MRSI) of the brain, using the recently proposed subspace-based spectroscopic imaging framework called SPICE (SPectroscopic Imaging by exploiting spatiospectral CorrElation).

THEORY AND METHODS

SPICE is characterized by the use of a subspace model for both data acquisition and image reconstruction. For data acquisition, we propose a novel spatiospectral encoding scheme that provides hybrid data sets for determining the subspace structure and for image reconstruction using the subspace model. More specifically, we use a hybrid chemical shift imaging /echo-planar spectroscopic imaging sequence for two-dimensional (2D) MRSI and a dual-density, dual-speed echo-planar spectroscopic imaging sequence for three-dimensional (3D) MRSI. For image reconstruction, we propose a method that can determine the subspace structure and the high-resolution spatiospectral reconstruction from the hybrid data sets generated by the proposed sequences, incorporating field inhomogeneity correction and edge-preserving regularization.

RESULTS

Phantom and in vivo brain experiments were performed to evaluate the performance of the proposed method. For 2D MRSI experiments, SPICE is able to produce high-SNR spatiospectral distributions with an approximately 3 mm nominal in-plane resolution from a 10-min acquisition. For 3D MRSI experiments, SPICE is able to achieve an approximately 3 mm in-plane and 4 mm through-plane resolution in about 25 min.

CONCLUSION

Special data acquisition and reconstruction methods have been developed for high-resolution (1) H-MRSI of the brain using SPICE. Using these methods, SPICE is able to produce spatiospectral distributions of (1) H metabolites in the brain with high spatial resolution, while maintaining a good SNR. These capabilities should prove useful for practical applications of SPICE. Magn Reson Med 76:1059-1070, 2016. © 2015 Wiley Periodicals, Inc.

Improved image reconstruction for subspace-based spectroscopic imaging using non-quadratic regularization

A new MR spectroscopic imaging method, called SPICE (SPectroscopic Imaging by exploiting spatiospectral CorrElation), has been recently proposed to enable high-resolution metabolic imaging with good SNR. A key problem within the SPICE framework is image reconstruction from a very noisy and sparsely sampled dataset. This paper addresses this problem by integrating the low-rank model used in SPICE reconstruction with a non-quadratic regularization. An efficient primal-dual based algorithm is described to solve the associated optimization problem. The proposed method has been validated using both simulation and phantom studies and is expected to enhance the unprecedented capability of SPICE for high-resolution metabolic imaging.

Improved Low-Rank Filtering of Magnetic Resonance Spectroscopic Imaging Data Corrupted by Noise and B₀ Field Inhomogeneity

GOAL

To improve the signal-to-noise ratio (SNR) of magnetic resonance spectroscopic imaging (MRSI) data.

METHODS

A low-rank filtering method recently proposed for denoising MRSI data is extended by: 1) incorporating tissue boundary constraints to enable local low-rank filtering, and 2) integrating B0 field inhomogeneity correction by rank-minimization to make the low-rank model more effective.

RESULTS

The proposed method was validated using both simulated and in vivo MRSI data. Its denoising performance is also compared with an upper bound based on the constrained Cramér-Rao lower bound for low-rank filtering.

CONCLUSION

Low-rank filtering can effectively improve the SNR of MRSI data corrupted by both noise and B0 field inhomogeneity.

SIGNIFICANCE

The proposed low-rank filtering method will enhance the practical utility of high-resolution MRSI, where SNR has been a limiting factor.

3D-printed Shepp-Logan phantom as a real-world benchmark for MRI

PURPOSE

As prescribed and reliable geometrical entities, phantoms have served as indispensable validation tools in a variety of MR-related topics. Though a number of phantoms have been conceived, certain applications may warrant highly customized geometries. The purpose of this study was to demonstrate the expediency of rapid prototyping for generating a flexible class of MR-compatible phantom designs.

METHODS

An incarnation of the three-dimensional Shepp-Logan numerical phantom, amended for use in magnetic resonance spectroscopic imaging, was actualized using rapid prototyping. Each of the comprising compartments was filled with a solution containing prepared concentrations of common (1)H brain metabolites. Analytical Fourier expressions for the phantom class were established in order to generate a set of simulated measurements, which were then contrasted with acquired data.

RESULTS

Experimental results for both structural and spectroscopic imaging substantiate the suitability of rapid prototyping for MR phantom applications. The analytically simulated measurements show excellent agreement with the measured data, but also highlight the various consequences effectuated when certain aspects of the acquisition model are disregarded or misrepresented.

CONCLUSION

Rapid prototyping offers a novel and versatile framework for MR phantom-based validation studies. Furthermore, the growing accessibility and open-source compatibility may provide an important link between the often disparate numerical and haptic testing.

A subspace approach to high-resolution spectroscopic imaging

PURPOSE

To accelerate spectroscopic imaging using sparse sampling of (k,t)-space and subspace (or low-rank) modeling to enable high-resolution metabolic imaging with good signal-to-noise ratio.

METHODS

The proposed method, called SPectroscopic Imaging by exploiting spatiospectral CorrElation, exploits a unique property known as partial separability of spectroscopic signals. This property indicates that high-dimensional spectroscopic signals reside in a very low-dimensional subspace and enables special data acquisition and image reconstruction strategies to be used to obtain high-resolution spatiospectral distributions with good signal-to-noise ratio. More specifically, a hybrid chemical shift imaging/echo-planar spectroscopic imaging pulse sequence is proposed for sparse sampling of (k,t)-space, and a low-rank model-based algorithm is proposed for subspace estimation and image reconstruction from sparse data with the capability to incorporate prior information and field inhomogeneity correction.

RESULTS

The performance of the proposed method has been evaluated using both computer simulations and phantom studies, which produced very encouraging results. For two-dimensional spectroscopic imaging experiments on a metabolite phantom, a factor of 10 acceleration was achieved with a minimal loss in signal-to-noise ratio compared to the long chemical shift imaging experiments and with a significant gain in signal-to-noise ratio compared to the accelerated echo-planar spectroscopic imaging experiments.

CONCLUSION

The proposed method, SPectroscopic Imaging by exploiting spatiospectral CorrElation, is able to significantly accelerate spectroscopic imaging experiments, making high-resolution metabolic imaging possible.