Truncation artifact reduction in spectroscopic imaging using a dual-density spiral k-space trajectory

Truncation artifacts arise in magnetic resonance spectroscopic imaging (MRSI) of the human brain due to limited coverage of k-space necessitated by low SNR of metabolite signal and limited scanning time. In proton MRSI of the head, intense extra-cranial lipid signals “bleed” into brain regions, thereby contaminating signals of metabolites therein. This work presents a data acquisition strategy for reducing truncation artifact based on extended k-space coverage achieved with a dual-SNR strategy. Using the fact that the SNR in k-space increases monotonically with sampling density, dual-SNR is achieved in an efficient manner with a dual-density spiral k-space trajectory that permits a smooth transition from high density to low density. The technique is demonstrated to be effective in reducing “bleeding” of extra-cranial lipid signals while preserving the SNR of metabolites in the brain.     

基于同伦l0范数最小化重建的三维动态磁共振成像

高欠采倍数的动态磁共振图像重建具有重要意义,是同时实现高时间分辨率和高空间分辨率动态对比度增强成像的重要环节.本研究提出一种结合黄金角变密度螺旋采样、并行成像和基于同伦l₀范数最小化的压缩感知的图像重建的三维动态磁共振成像方法.黄金角变密度螺旋采样轨迹被用来连续获取k空间数据,具有数据采集效率高、对运动不敏感等优点.在重建算法中,将多线圈稀疏约束应用于时间总变分域,使用基于l₀范数最小化的非线性重建算法代替传统的l₁范数最小化算法,进一步提高了欠采样率.仿真实验和在体实验表明本文所提的方法在保持图像质量的同时,也可以实现较高的空间分辨率和时间分辨率,初步验证了基于同伦l₀范数最小化重建在三维动态磁共振成像上的优势和临床价值.     

An improved 1H magnetic resonance spectroscopic imaging technique for the human breast: preliminary results

The high sensitivity but poor specificity of magnetic resonance imaging for detecting breast cancer has stimulated interest in magnetic resonance spectroscopic imaging (MRSI) as a tool to improve specificity and reduce the number of benign biopsies. The challenge of applying 1H MRSI to the diagnosis of cancer in the human breast is the need for robust lipid suppression and a clinically acceptable acquisition time. We present an improved 1H MRSI technique that uses an independently optimized chemical-shift-selective for lipid suppression and weighted elliptical k-space sampling combined with a Hamming filter for improved sampling efficiency.     

Three dimension double inversion recovery gray matter imaging using compressed sensing

The double inversion recovery (DIR) imaging technique has various applications such as black blood magnetic resonance imaging and gray/white matter imaging. Recent clinical studies show the promise of DIR for high resolution three dimensional (3D) gray matter imaging. One drawback in this case however is the long data acquisition time needed to obtain the fully sampled 3D spatial frequency domain (k-space) data. In this paper, we propose a method to solve this problem using the compressed sensing (CS) algorithm with contourlet transform. The contourlet transform is an effective sparsifying transform especially for images with smooth contours. Therefore, we applied this algorithm to undersampled DIR images and compared with a CS algorithm using wavelet transform by evaluating the reconstruction performance of each algorithm for undersampled k-space data. The results show that the proposed CS algorithm achieves a more accurate reconstruction in terms of the mean structural similarity index and root mean square error than the CS algorithm using wavelet transform.     

Phase correction-based singularity function analysis for partial k-space reconstruction

Partial k-space acquisition is a conventional method in magnetic resonance imaging (MRI) for reducing imaging time while maintaining image quality. In this field, image reconstruction from partial k-space is a key issue. This paper proposes an approach fundamentally different from traditional techniques for reconstructing magnetic resonance (MR) images from partial k-space. It uses a so-called singularity function analysis (SFA) model based on phase correction. With such a reconstruction approach, some nonacquired negative spatial frequencies are first recovered by means of phase correction and Hermitian symmetry property, and then the other nonacquired negative and/or positive spatial frequencies are estimated using the mathematical SFA model. The method is particularly suitable for asymmetrical partial k-space acquisition owing to its ability of overcoming reconstruction limitations due to k-space truncations. The performance of this approach is evaluated using both simulated and real MR brain images, and compared with existing techniques. The results demonstrate that the proposed SFA based on phase correction achieves higher image quality than the initial SFA or the projection-onto-convex sets (POCS) method.     

Fast nosologic imaging of the brain

Magnetic resonance spectroscopic imaging (MRSI) provides information about the spatial metabolic heterogeneity of an organ in the human body. In this way, MRSI can be used to detect tissue regions with abnormal metabolism, e.g. tumor tissue. The main drawback of MRSI in clinical practice is that the analysis of the data requires a lot of expertise from the radiologists. In this article, we present an automatic method that assigns each voxel of a spectroscopic image of the brain to a histopathological class. The method is based on Canonical Correlation Analysis (CCA), which has recently been shown to be a robust technique for tissue typing. In CCA, the spectral as well as the spatial information about the voxel is used to assign it to a class. This has advantages over other methods that only use spectral information since histopathological classes are normally covering neighbouring voxels. In this paper, a new CCA-based method is introduced in which MRSI and MR imaging information is integrated. The performance of tissue typing is compared for CCA applied to the whole MR spectra and to sets of features obtained from the spectra. Tests on simulated and in vivo MRSI data show that the new method is very accurate in terms of classification and segmentation. The results also show the advantage of combining spectroscopic and imaging data.     

Elliptical magnetic resonance spectroscopic imaging with GRAPPA for imaging brain tumors at 3 T

Magnetic Resonance Spectroscopic Imaging (MRSI) is a technique for imaging spatial variation of metabolites and has been very useful in characterizing biochemical changes associated with disease as well as response to therapy in malignant pathologies. This work presents a self-calibrated undersampling to accelerate 3D elliptical MRSI and an extrapolation-reconstruction algorithm based on the GRAPPA method. The accelerated MRSI technique was tested in three volunteers and five brain tumor patients. Acceleration allowed larger spatial coverage and consequently, less lipid contamination in spectra, compared to fully sampled acquisition within the same scantime. Metabolite concentrations measured from the accelerated acquisitions were in good agreement with measurements obtained from fully sampled MRSI scans.     

Phased array 3D MR spectroscopic imaging of the brain at 7 T

Ultra-high-field 7 T magnetic resonance (MR) scanners offer the potential for greatly improved MR spectroscopic imaging due to increased sensitivity and spectral resolution. Prior 7 T human single-voxel MR Spectroscopy (MRS) studies have shown significant increases in signal-to-noise ratio (SNR) and spectral resolution as compared to lower magnetic fields but have not demonstrated the increase in spatial resolution and multivoxel coverage possible with 7 T MR spectroscopic imaging. The goal of this study was to develop specialized radiofrequency (RF) pulses and sequences for three-dimensional (3D) MR spectroscopic imaging (MRSI) at 7 T to address the challenges of increased chemical shift misregistration, B1 power limitations, and increased spectral bandwidth. The new 7 T MRSI sequence was tested in volunteer studies and demonstrated the feasibility of obtaining high-SNR phased-array 3D MRSI from the human brain.     

A correction algorithm for undersampled images using dynamic segmentation and entropy based focus criterion

A post-processing technique is presented for correcting images undersampled in k-space. The method works by taking advantage of the image's background zeros (dynamically segmented through the application of a threshold) to extrapolate the missing k-space samples. The algorithm can produce good quality images from a small set of k-space frequencies with only a few iterations of simple matrix operations, using the image entropy as the focus criterion. It does not require any special patient preparation, extra pulse sequences, complex gradient programming or specialized hardware. This makes it a good candidate for any application that requires short scan times or where only few frequencies can be sampled.     

Exploiting the wavelet structure in compressed sensing MRI

Sparsity has been widely utilized in magnetic resonance imaging (MRI) to reduce k-space sampling. According to structured sparsity theories, fewer measurements are required for tree sparse data than the data only with standard sparsity. Intuitively, more accurate image reconstruction can be achieved with the same number of measurements by exploiting the wavelet tree structure in MRI. A novel algorithm is proposed in this article to reconstruct MR images from undersampled k-space data. In contrast to conventional compressed sensing MRI (CS-MRI) that only relies on the sparsity of MR images in wavelet or gradient domain, we exploit the wavelet tree structure to improve CS-MRI. This tree-based CS-MRI problem is decomposed into three simpler subproblems then each of the subproblems can be efficiently solved by an iterative scheme. Simulations and in vivo experiments demonstrate the significant improvement of the proposed method compared to conventional CS-MRI algorithms, and the feasibleness on MR data compared to existing tree-based imaging algorithms.     

Spatially encoded NMR and the acquisition of 2D magnetic resonance images within a single scan

An approach that enables the acquisition of multidimensional NMR spectra within a single scan has been recently proposed and demonstrated. The present paper explores the applicability of such ultrafast acquisition schemes toward the collection of two-dimensional magnetic resonance imaging (2D MRI) data. It is shown that ideas enabling the application of these spatially encoded schemes within a spectroscopic setting, can be extended in a straightforward manner to pure imaging. Furthermore, the reliance of the original scheme on a spatial encoding and subsequent decoding of the evolution frequencies endows imaging applications with a greater simplicity and flexibility than their spectroscopic counterparts. The new methodology also offers the possibility of implementing the single-scan acquisition of 2D MRI images using sinusoidal gradients, without having to resort to subsequent interpolation procedures or non-linear sampling of the data. Theoretical derivations on the operational principles and imaging characteristics of a number of sequences based on these ideas are derived, and experimentally validated with a series of 2D MRI results collected on a variety of model phantom samples.     

High spatial resolution in vivo 2D (1)H magnetic resonance spectroscopic imaging of human muscles with a band-selective technique.

This report demonstrates a 2D (1)H magnetic resonance spectroscopic imaging (MRSI) technique that can address some technical difficulties often encountered in MRS studies of human muscles. A preliminary application of this whole-slice technique in human skeletal muscles demonstrates clearly noticeable differences in (1)H metabolite spectra between different human muscles. This observation illustrates the importance of multi-voxel and high spatial resolution in a heterogeneous environment. This technique is robust, can be easily implemented on a commercial MR scanner, and should prove useful for investigators in both basic and clinical (1)H MRS studies.     

Magnetic resonance image reconstruction using trained geometric directions in 2D redundant wavelets domain and non-convex optimization

Reducing scanning time is significantly important for MRI. Compressed sensing has shown promising results by undersampling the k-space data to speed up imaging. Sparsity of an image plays an important role in compressed sensing MRI to reduce the image artifacts. Recently, the method of patch-based directional wavelets (PBDW) which trains geometric directions from undersampled data has been proposed. It has better performance in preserving image edges than conventional sparsifying transforms. However, obvious artifacts are presented in the smooth region when the data are highly undersampled. In addition, the original PBDW-based method does not hold obvious improvement for radial and fully 2D random sampling patterns. In this paper, the PBDW-based MRI reconstruction is improved from two aspects: 1) An efficient non-convex minimization algorithm is modified to enhance image quality; 2) PBDW are extended into shift-invariant discrete wavelet domain to enhance the ability of transform on sparsifying piecewise smooth image features. Numerical simulation results on vivo magnetic resonance images demonstrate that the proposed method outperforms the original PBDW in terms of removing artifacts and preserving edges.     

3D sensitivity encoded ellipsoidal MR spectroscopic imaging of gliomas at 3T

Purpose

The goal of this study was to implement time efficient data acquisition and reconstruction methods for 3D magnetic resonance spectroscopic imaging (MRSI) of gliomas at a field strength of 3T using parallel imaging techniques.

Methods

The point spread functions, signal to noise ratio (SNR), spatial resolution, metabolite intensity distributions and Cho:NAA ratio of 3D ellipsoidal, 3D sensitivity encoding (SENSE) and 3D combined ellipsoidal and SENSE (e-SENSE) k-space sampling schemes were compared with conventional k-space data acquisition methods.

Results

The 3D SENSE and e-SENSE methods resulted in similar spectral patterns as the conventional MRSI methods. The Cho:NAA ratios were highly correlated (P<.05 for SENSE and P<.001 for e-SENSE) with the ellipsoidal method and all methods exhibited significantly different spectral patterns in tumor regions compared to normal appearing white matter. The geometry factors ranged between 1.2 and 1.3 for both the SENSE and e-SENSE spectra. When corrected for these factors and for differences in data acquisition times, the empirical SNRs were similar to values expected based upon theoretical grounds. The effective spatial resolution of the SENSE spectra was estimated to be same as the corresponding fully sampled k-space data, while the spectra acquired with ellipsoidal and e-SENSE k-space samplings were estimated to have a 2.36–2.47-fold loss in spatial resolution due to the differences in their point spread functions.

Conclusion

The 3D SENSE method retained the same spatial resolution as full k-space sampling but with a 4-fold reduction in scan time and an acquisition time of 9.28 min. The 3D e-SENSE method had a similar spatial resolution as the corresponding ellipsoidal sampling with a scan time of 4:36 min. Both parallel imaging methods provided clinically interpretable spectra with volumetric coverage and adequate SNR for evaluating Cho, Cr and NAA.     

Visually weighted reconstruction of compressive sensing MRI

Compressive sensing (CS) enables the reconstruction of a magnetic resonance (MR) image from undersampled data in k-space with relatively low-quality distortion when compared to the original image. In addition, CS allows the scan time to be significantly reduced. Along with a reduction in the computational overhead, we investigate an effective way to improve visual quality through the use of a weighted optimization algorithm for reconstruction after variable density random undersampling in the phase encoding direction over k-space. In contrast to conventional magnetic resonance imaging (MRI) reconstruction methods, the visual weight, in particular, the region of interest (ROI), is investigated here for quality improvement. In addition, we employ a wavelet transform to analyze the reconstructed image in the space domain and fully utilize data sparsity over the spatial and frequency domains. The visual weight is constructed by reflecting the perceptual characteristics of the human visual system (HVS), and then applied to ?₁ norm minimization, which gives priority to each coefficient during the reconstruction process. Using objective quality assessment metrics, it was found that an image reconstructed using the visual weight has higher local and global quality than those processed by conventional methods.     

Fast multivoxel two-dimensional spectroscopic imaging at 3 T

The utility of multivoxel two-dimensional chemical shift imaging in the clinical environment will ultimately be determined by the imaging time and the metabolite peaks that can be detected. Different k-space sampling schemes can be characterized by their minimum required imaging time. The use of spiral-based readout gradients effectively reduces the minimum scan time required due to simultaneous data acquisition in three k-space dimensions (k(x), k(y) and k(f(2))). A 3-T spiral-based multivoxel two-dimensional spectroscopic imaging sequence using the PRESS excitation scheme was implemented. Good performance was demonstrated by acquiring preliminary in vivo data for applications, including brain glutamate imaging, metabolite T(2) quantification and high-spatial-resolution prostate spectroscopic imaging. All protocols were designed to acquire data within a 17-min scan time at a field strength of 3 T.     

Perfusion MR imaging and proton MR spectroscopic imaging in differentiating necrotizing cerebritis from glioblastoma multiforme

We describe a lesion with the magnetic resonance imaging (MRI) characteristics of a glioblastoma mutiforme and demonstrate how perfusion MRI and proton MR spectroscopic imaging can be used to differentiate necrotizing cerebritis from what appeared to be a high-grade glioma. A 43-year-old woman presented to her physician complaining of progressive visual disturbance and headache for several weeks. Conventional MRI demonstrated a parietal peripherally enhancing mass with central necrosis and moderate to severe surrounding T2 hyperintensity, suggesting an infiltrating high-grade glioma. However, advanced imaging, including dynamic susceptibility contrast MRI (DSC MRI) and magnetic resonance spectroscopic imaging (MRSI), suggested a nonneoplastic lesion. The DSC MRI data demonstrated no hyperperfusion within the lesion and surrounding T2 signal abnormality, and the MRSI data showed overall decrease in metabolites in this region, except for lactate. Because of the aggressive appearance to the lesion and the patients' worsening symptoms, a biopsy was performed. The pathologic diagnosis was necrotizing cerebritis. After the commencement of steroid therapy, imaging findings and patient symptoms improved. This report will review the utility of advanced imaging for differentiating inflammatory from neoplastic appearing lesions on conventional imaging.     

MR imaging of flow with locally high spatial resolution.

Locally focused magnetic resonance imaging (LF MRI) allows imaging with variable spatial resolution within the field of view (FOV). Because LF MRI uses a priori information to provide locally high resolution in regions with rapid spatial variations in intensity (e.g., blood/tissue interface), it allows accurate reproduction of intense sharp edges in the specimen without blurring and truncation artifacts. This study employs LF MRI for 3D imaging of stationary and pulsatile flow. In the implemented version of LF MRI analytically defined basis functions are used to determine image intensity in regions depicted with low or high resolution. It is demonstrated that LF MRI of flow allows a significant (i.e. 3-4 times) reduction in scan time as compared to conventional FT MRI. It is also shown that LF images of pulsatile flow have a decreased appearance of ghosting artifacts as compared to the images reconstructed by using the conventional method.     

Phase-difference and spectroscopic imaging for monitoring of human brain temperature during cooling

Decrease of the human brain temperature was induced by intranasal cooling. The main purpose of this study was to compare the two magnetic resonance methods for monitoring brain temperature changes during cooling: phase-difference and magnetic resonance spectroscopic imaging (MRSI) with high spatial resolution. Ten healthy volunteers were measured. Selective brain cooling was performed through nasal cavities using saline-cooled balloon catheters. MRSI was based on a radiofrequency spoiled gradient echo sequence. The spectral information was encoded by incrementing the echo time of the subsequent eight image records. Reconstructed voxel size was 1×1×5 mm³. Relative brain temperature was computed from the positions of water spectral lines. Phase maps were obtained from the first image record of the MRSI sequence. Mild hypothermia was achieved in 15–20 min. Mean brain temperature reduction varied in the interval <−3.0; − 0.6>°C and <−2.7; − 0.7>°C as measured by the MRSI and phase-difference methods, respectively. Very good correlation was found in all locations between the temperatures measured by both techniques except in the frontal lobe. Measurements in the transversal slices were more robust to the movement artifacts than those in the sagittal planes. Good agreement was found between the MRSI and phase-difference techniques.     

Undersampled MRI reconstruction with patch-based directional wavelets

Compressed sensing has shown great potential in reducing data acquisition time in magnetic resonance imaging (MRI). In traditional compressed sensing MRI methods, an image is reconstructed by enforcing its sparse representation with respect to a preconstructed basis or dictionary. In this paper, patch-based directional wavelets are proposed to reconstruct images from undersampled k-space data. A parameter of patch-based directional wavelets, indicating the geometric direction of each patch, is trained from the reconstructed image using conventional compressed sensing MRI methods and incorporated into the sparsifying transform to provide the sparse representation for the image to be reconstructed. A reconstruction formulation is proposed and solved via an efficient alternating direction algorithm. Simulation results on phantom and in vivo data indicate that the proposed method outperforms conventional compressed sensing MRI methods in preserving the edges and suppressing the noise. Besides, the proposed method is not sensitive to the initial image when training directions.