Cancellation of streak artifacts in radial abdominal imaging using interference null space projection.

PURPOSE: In radial abdominal imaging, it has been commonly observed that signal from the arms cause streaks due to system imperfections. We previously introduced a streak removal technique (B-STAR), which is inherently spatially variant and limited to work in image space. In this work, we propose a spatially invariant streak cancellation technique (CACTUS), which can be applied in either image space or k-space and is compatible with iterative reconstructions. THEORY AND METHODS: Streak sources are typically spatially localized and can be represented using a low-dimensional subspace. CACTUS identifies the streak subspace by leveraging the spatial redundancy of receiver coils and projects the data onto the streak null space to eliminate the streaks. When applied in k-space, CACTUS can be combined with iterative reconstructions. CACTUS was tested in phantoms and in vivo abdominal imaging using a radial turbo spin-echo pulse sequence. RESULTS: In phantoms, CACTUS improved T2 estimation in comparison to previous de-streaking methods. In vivo experiments showed that CACTUS reduced streaks and yielded T2 estimation, in regions affected by streaks, closer to a streak-free reference. Evaluation using a clinical abdominal dataset (n = 20) showed that CACTUS is comparable to B-STAR and yields significantly better signal preservation and streak cancellation than coil removal and suppression methods. CONCLUSION: CACTUS provides superior signal preservation and streak reduction performance compared to coil removal and suppression methods. As a clear advantage over B-STAR, CACTUS can be integrated with iterative reconstruction methods. In abdominal T2 mapping, CACTUS improves the accuracy of parameter estimation in areas affected by streaks.

Abdominal T2-Weighted Imaging and T2 Mapping Using a Variable Flip Angle Radial Turbo Spin-Echo Technique.

BACKGROUND: T2 mapping is of great interest in abdominal imaging but current methods are limited by low resolution, slice coverage, motion sensitivity, or lengthy acquisitions. PURPOSE: Develop a radial turbo spin-echo technique with refocusing variable flip angles (RADTSE-VFA) for high spatiotemporal T2 mapping and efficient slice coverage within a breath-hold and compare to the constant flip angle counterpart (RADTSE-CFA). STUDY TYPE: Prospective technical efficacy. SUBJECTS: Testing performed on agarose phantoms and 12 patients. Focal liver lesion classification tested on malignant (N = 24) and benign (N = 11) lesions. FIELD STRENGTH/SEQUENCE: 1.5 T/RADTSE-VFA, RADTSE-CFA. ASSESSMENT: A constrained objective function was used to optimize the refocusing flip angles. Phantom and/or in vivo data were used to assess relative contrast, T2 estimation, specific absorption rate (SAR), and focal liver lesion classification. STATISTICAL TESTS: t-Tests or Mann-Whitney Rank Sum tests were used. RESULTS: Phantom data did not show significant differences in mean relative contrast (P = 0.10) and T2 accuracy (P = 0.99) between RADTSE-VFA and RADTSE-CFA. Adding noise caused T2 overestimation predominantly for RADTSE-CFA and low T2 values. In vivo results did not show significant differences in mean spleen-to-liver (P = 0.62) and kidney-to-liver (P = 0.49) relative contrast between RADTSE-VFA and RADTSE-CFA. Mean T2 values were not significantly different between the two techniques for spleen (T2VFA  = 109.2 ± 12.3 msec; T2CFA  = 110.7 ± 11.1 msec; P = 0.78) and kidney-medulla (T2VFA  = 113.0 ± 8.7 msec; T2CFA  = 114.0 ± 8.6 msec; P = 0.79). Liver T2 was significantly higher for RADTSE-CFA (T2VFA  = 52.6 ± 6.6 msec; T2CFA  = 60.4 ± 8.0 msec) consistent with T2 overestimation in the phantom study. Focal liver lesion classification had comparable T2 distributions for RADTSE-VFA and RADTSE-CFA for malignancies (P = 1.0) and benign lesions (P = 0.39). RADTSE-VFA had significantly lower SAR than RADTSE-CFA increasing slice coverage by 1.5. DATA CONCLUSION: RADTSE-VFA provided noise-robust T2 estimation compared to the constant flip angle counterpart while generating T2-weighted images with comparable contrast. The VFA scheme minimized SAR improving slice efficiency for breath-hold imaging. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 1.

Joint total variation-based reconstruction of multiparametric magnetic resonance images for mapping tissue types.

Multispectral analysis of coregistered multiparametric magnetic resonance (MR) images provides a powerful method for tissue phenotyping and segmentation. Acquisition of a sufficiently varied set of multicontrast MR images and parameter maps to objectively define multiple normal and pathologic tissue types can require long scan times. Accelerated MRI on clinical scanners with multichannel receivers exploits techniques such as parallel imaging, while accelerated preclinical MRI scanning must rely on alternate approaches. In this work, tumor-bearing mice were imaged at 7 T to acquire k-space data corresponding to a series of images with varying T1-, T2- and T2*-weighting. A joint reconstruction framework is proposed to reconstruct a series of T1-weighted images and corresponding T1 maps simultaneously from undersampled Cartesian k-space data. The ambiguity introduced by undersampling was resolved by using model-based constraints and structural information from a reference fully sampled image as the joint total variation prior. This process was repeated to reconstruct T2-weighted and T2*-weighted images and corresponding maps of T2 and T2* from undersampled Cartesian k-space data. Validation of the reconstructed images and parameter maps was carried out by computing tissue-type maps, as well as maps of the proton density fat fraction (PDFF), proton density water fraction (PDwF), fat relaxation rate ( R2f*) and water relaxation rate ( R2w*) from the reconstructed data, and comparing them with ground truth (GT) equivalents. Tissue-type maps computed using 18% k-space data were visually similar to GT tissue-type maps, with dice coefficients ranging from 0.43 to 0.73 for tumor, fluid adipose and muscle tissue types. The mean T1 and T2 values within each tissue type computed using only 18% k-space data were within 8%-10% of the GT values from fully sampled data. The PDFF and PDwF maps computed using 27% k-space data were within 3%-15% of GT values and showed good agreement with the expected values for the four tissue types.

Visibility of Quantization Errors in Reversible JPEG2000.

Image compression systems that exploit the properties of the human visual system have been studied extensively over the past few decades. For the JPEG2000 image compression standard, all previous methods that aim to optimize perceptual quality have considered the irreversible pipeline of the standard. In this work, we propose an approach for the reversible pipeline of the JPEG2000 standard. We introduce a new methodology to measure visibility of quantization errors when reversible color and wavelet transforms are employed. Incorporation of the visibility thresholds using this methodology into a JPEG2000 encoder enables creation of scalable codestreams that can provide both near-threshold and numerically lossless representations, which is desirable in applications where restoration of original image samples is required. Most importantly, this is the first work that quantifies the bitrate penalty incurred by the reversible transforms in near-threshold image compression compared to the irreversible transforms.

Radial streak artifact reduction using phased array beamforming.

PURPOSE: A new method for streak artifact reduction in radial MRI based on phased array filtering. THEORY: Radial imaging in applications that require large fields-of-view can be susceptible to streaking artifacts due to gradient nonlinearities. Coil removal methods prune the coils contributing the most to streaking artifacts at the expense of signal loss. Phased array beamforming is a form of spatial filtering used to suppress unwanted signals. The proposed method uses interference covariance generated from the streaking artifact samples which are manually extracted with phased array beamforming to suppress streaking in the images. METHODS: The performance of the proposed method was evaluated on abdomen radial fast spin echo images acquired on a 1.5T Siemens scanner and compared with previously proposed methods. RESULTS: Our results demonstrate that the proposed method can effectively suppress streaking artifacts without any noticeable loss in signal levels. Coil removal methods can suppress streaks as well but they may incur significant signal loss due to coil pruning. Quantitative metrics also demonstrate the superiority of the proposed method over earlier methods. CONCLUSION: The use of interference covariance with phased array beamforming can help reduce streaking artifacts.

An efficient 3D stack-of-stars turbo spin echo pulse sequence for simultaneous T2-weighted imaging and T2 mapping.

PURPOSE: To design a pulse sequence for efficient 3D T2-weighted imaging and T2 mapping. METHODS: A stack-of-stars turbo spin echo pulse sequence with variable refocusing flip angles and a flexible pseudorandom view ordering is proposed for simultaneous T2-weighted imaging and T2 mapping. An analytical framework is introduced for the selection of refocusing flip angles to maximize relative tissue contrast while minimizing T2 estimation errors and maintaining low specific absorption rate. Images at different echo times are generated using a subspace constrained iterative reconstruction algorithm. T2 maps are obtained by modeling the signal evolution using the extended phase graph model. The technique is evaluated using phantoms and demonstrated in vivo for brain, knee, and carotid imaging. RESULTS: Numerical simulations demonstrate an improved point spread function with the proposed pseudorandom view ordering compared to golden angle view ordering. Phantom experiments show that T2 values estimated from the stack-of-stars turbo spin echo pulse sequence with variable refocusing flip angles have good concordance with spin echo reference values. In vivo results show the proposed pulse sequence can generate qualitatively comparable T2-weighted images as conventional Cartesian 3D SPACE in addition to simultaneously generating 3D T2 maps. CONCLUSION: The proposed stack-of-stars turbo spin echo pulse sequence with pseudorandom view ordering and variable refocusing flip angles allows high resolution isotropic T2 mapping in clinically acceptable scan times. The optimization framework for the selection of refocusing flip angles improves T2 estimation accuracy while generating T2-weighted contrast comparable to conventional Cartesian imaging.

Rapid high-resolution T1 mapping using a highly accelerated radial steady-state free-precession technique.

BACKGROUND: T1 mapping is often used in some clinical protocols. Existing techniques are limited in slice coverage, and/or spatial-temporal resolution, or require long acquisitions. Here we present a multi-slice inversion-recovery (IR) radial steady-state free precession (radSSFP) pulse sequence combined with a principal component (PC) based reconstruction that overcomes these limitations. PURPOSE: To develop a fast technique for multi-slice high-resolution T1 mapping. STUDY TYPE: Technical efficacy study done prospectively. PHANTOM/SUBJECTS: IR-radSSFP was tested in phantoms, five healthy volunteers, and four patients with abdominal lesions. FIELD STRENGTH/SEQUENCE: IR-radSSFP was implemented at 3T. ASSESSMENT: Computer simulations were performed to optimize the flip angle for T1 estimation; testing was done in phantoms using as reference an IR spin-echo pulse sequence. T1 mapping with IR-radSSFP was also assessed in vivo (brain and abdomen) and T1 values were compared with literature. T1 maps were also compared with a radial IR-FLASH technique. STATISTICAL TESTS: A two-tailed t-test was used to compare T1 values in phantoms. A repeatability study was carried out in vivo using Bland-Altman analysis. RESULTS: Simulations and phantom experiments showed that a flip angle of 20˚ was optimal for T1 mapping. When comparing single to multi-slice experiments in phantoms there were no significant differences between the means T1 values (P = 0.0475). In vivo results show that T1 maps with spatial resolution as high as 0.69 mm × 0.69 mm × 2.00 mm (brain) and 0.83 mm × 0.83 mm × 3.00 mm (abdomen) can be generated for 84 brain slices in 3 min and 10 abdominal slices in a breath-hold; T1 values were comparable to those reported in literature. The coefficients of variation from the repeatability study were 1.7% for brain and 2.5-2.7% in the abdomen. DATA CONCLUSION: A multi-slice IR-radSSFP technique combined with a PC-based reconstruction was demonstrated for higher resolution T1 mapping. This technique is fast, motion-insensitive and yields repeatable T1 values comparable to those in literature. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:239-252.

Accelerated MR parameter mapping with a union of local subspaces constraint.

PURPOSE: A new reconstruction method for multi-contrast imaging and parameter mapping based on a union of local subspaces constraint is presented. THEORY: Subspace constrained reconstructions use a predetermined subspace to explicitly constrain the relaxation signals. The choice of subspace size ( K ) impacts the approximation error vs noise-amplification tradeoff associated with these methods. A different approach is used in the model consistency constraint (MOCCO) framework to leverage the subspace model to enforce a softer penalty. Our proposed method, MOCCO-LS, augments the MOCCO model with a union of local subspaces (LS) approach. The union of local subspaces model is coupled with spatial support constraints and incorporated into the MOCCO framework to regularize the contrast signals in the scene. METHODS: The performance of the MOCCO-LS method was evaluated in vivo on T1 and T2 mapping of the human brain and with Monte-Carlo simulations and compared against MOCCO and the explicit subspace constrained models. RESULTS: The results demonstrate a clear improvement in the multi-contrast images and parameter maps. We sweep across the model order space ( K ) to compare the different reconstructions and demonstrate that the reconstructions have different preferential operating points. Experiments on T2 mapping show that the proposed method yields substantial improvements in performance even when operating at very high acceleration rates. CONCLUSIONS: The use of a union of local subspace constraints coupled with a sparsity promoting penalty leads to improved reconstruction quality of multi-contrast images and parameter maps.

A multi-band double-inversion radial fast spin-echo technique for T2 cardiovascular magnetic resonance mapping of the heart.

BACKGROUND: Double inversion recovery (DIR) fast spin-echo (FSE) cardiovascular magnetic resonance (CMR) sequences are used clinically for black-blood T2-weighted imaging. However, these sequences suffer from slice inefficiency due to the non-selective inversion pulses. We propose a multi-band (MB) encoded DIR radial FSE (MB-DIR-RADFSE) technique to simultaneously excite two slices. This sequence has improved signal-to-noise ratio per unit time compared to a single slice excitation. It is also motion robust and enables the reconstruction of high-resolution black-blood T2-weighted images and T2 maps for the excited slices. METHODS: Hadamard encoded MB pulses were used in MB-DIR-RADFSE to simultaneously excite two slices. A principal component based iterative reconstruction was used to jointly reconstruct black-blood T2-weighted images and T2 maps. Phantom and in vivo experiments were performed to evaluate T2 mapping performance and results were compared to a T2-prepared balanced steady state free precession (bSSFP) method. The inter-segment variability of the T2 maps were assessed using data acquired on healthy subjects. A reproducibility study was performed to evaluate reproducibility of the proposed technique. RESULTS: Phantom experiments show that the T2 values estimated from MB-DIR-RADFSE are comparable to the spin-echo based reference, while T2-prepared bSSFP over-estimated T2 values. The relative contrast of the black-blood images from the multi-band scheme was comparable to those from a single slice acquisition. The myocardial segment analysis on 8 healthy subjects indicated a significant difference (p-value < 0.01) in the T2 estimates from the apical slice when compared to the mid-ventricular slice. The mean T2 estimate from 12 subjects obtained using T2-prepared bSSFP was significantly higher (p-value = 0.012) compared to MB-DIR-RADFSE, consistent with the phantom results. The Bland-Altman analysis showed excellent reproducibility between the MB-DIR-RADFSE measurements, with a mean T2 difference of 0.12 ms and coefficient of reproducibility of 2.07 in 15 clinical subjects. The utility of this technique is demonstrated in two subjects where the T2 maps show elevated values in regions of pathology. CONCLUSIONS: The use of multi-band pulses for excitation improves the slice efficiency of the double inversion fast spin-echo pulse sequence. The use of a radial trajectory and a joint reconstruction framework allows reconstruction of TE images and T2 maps for the excited slices.

Multiresolution imaging using golden angle stack-of-stars and compressed sensing for dynamic MR urography.

PURPOSE: To develop a novel multiresolution MRI methodology for accurate estimation of glomerular filtration rate (GFR) in vivo. MATERIALS AND METHODS: A three-dimensional golden-angle radial stack-of-stars (SoS) trajectory was used for data acquisition on a 3 Tesla MRI scanner. Multiresolution reconstruction and analysis was performed using arterial input function reconstructed at 1-s. temporal resolution and renal dynamic data reconstructed using compressed sensing (CS) with 4-s temporal resolution. The method was first validated using simulations and the clinical utility of the technique was evaluated by comparing the GFR estimates from the proposed method to the estimated GFR (eGFR) obtained from serum creatinine for 10 subjects. RESULTS: The 4-s temporal resolution CS images minimized streaking artifacts and noise while the 1-s temporal resolution AIF minimized errors in GFR estimates. A paired t-test showed that there was no statistically significant difference between MRI based total GFR values and serum creatinine based eGFR estimates (P = 0.92). CONCLUSION: We have demonstrated the feasibility of multiresolution MRI using a golden angle radial stack-of-stars scheme to accurately estimate GFR as well as produce diagnostic quality dynamic images in vivo. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2017;46:303-311.

Volumetric MRI of the lungs during forced expiration.

PURPOSE: Lung function is typically characterized by spirometer measurements, which do not offer spatially specific information. Imaging during exhalation provides spatial information but is challenging due to large movement over a short time. The purpose of this work is to provide a solution to lung imaging during forced expiration using accelerated magnetic resonance imaging. The method uses radial golden angle stack-of-stars gradient echo acquisition and compressed sensing reconstruction. METHODS: A technique for dynamic three-dimensional imaging of the lungs from highly undersampled data is developed and tested on six subjects. This method takes advantage of image sparsity, both spatially and temporally, including the use of reference frames called bookends. Sparsity, with respect to total variation, and residual from the bookends, enables reconstruction from an extremely limited amount of data. RESULTS: Dynamic three-dimensional images can be captured at sub-150 ms temporal resolution, using only three (or less) acquired radial lines per slice per timepoint. The images have a spatial resolution of 4.6×4.6×10 mm. Lung volume calculations based on image segmentation are compared to those from simultaneously acquired spirometer measurements. CONCLUSION: Dynamic lung imaging during forced expiration is made possible by compressed sensing accelerated dynamic three-dimensional radial magnetic resonance imaging. Magn Reson Med 75:2295-2302, 2016. © 2015 Wiley Periodicals, Inc.

Correcting partial volume effects in biexponential T2 estimation of small lesions.

PURPOSE: T2 mapping provides a quantitative approach for focal liver lesion characterization. For small lesions, a biexponential model should be used to account for partial volume effects (PVE). However, conventional biexponential fitting suffers from large uncertainty of the fitted parameters when noise is present. The purpose of this work is to develop a more robust method to correct for PVE affecting small lesions. METHODS: We developed a region of interest-based joint biexponential fitting (JBF) algorithm to estimate the T2 of lesions affected by PVE. JBF takes advantage of the lesion fraction variation among voxels within a region of interest. JBF is compared to conventional approaches using Cramér-Rao lower bound analysis, numerical simulations, phantom, and in vivo data. RESULTS: JBF provides more accurate and precise T2 estimates in the presence of PVE. Furthermore, JBF is less sensitive to region of interest drawing. Phantom and in vivo results show that JBF can be combined with a reconstruction method for highly undersampled data, enabling the characterization of small abdominal lesions from data acquired in a single breath hold. CONCLUSION: The JBF algorithm provides more accurate and stable T2 estimates for small structures than conventional techniques when PVE is present. It should be particularly useful for the characterization of small abdominal lesions.

High-resolution, large dynamic range field map estimation.

PURPOSE: We present a theory and a corresponding method to compute high-resolution field maps over a large dynamic range. THEORY AND METHODS: We derive a closed-form expression for the error in the field map value when computed from two echoes. We formulate an optimization problem to choose three echo times which result in a pair of maximally distinct error distributions. We use standard field mapping sequences at the prescribed echo times. We then design a corresponding estimation algorithm which takes advantage of the optimized echo times to disambiguate the field offset value. RESULTS: We validate our method using high-resolution images of a phantom at 7T. The resulting field maps demonstrate robust mapping over both a large dynamic range, and in low SNR regions. We also present high-resolution offset maps in vivo using both, GRE and multiecho gradient echo sequences. Even though the proposed echo time spacings are larger than the well known phase aliasing cutoff, the resulting field maps exhibit a large dynamic range without the use of phase unwrapping or spatial regularization techniques. CONCLUSION: We demonstrate a novel three-echo field map estimation method which overcomes the traditional noise-dynamic range trade-off.

Resveratrol: inhibitory effects on metastatic cell behaviors and voltage-gated Na⁺ channel activity in rat prostate cancer in vitro.

Resveratrol, a natural plant phenolic found at high concentration in red grapes, has been suggested to have a range of health benefits. Here, we tested its effects on metastatic cell behaviors. The strongly metastatic rat prostate MAT-LyLu cells were used as a model. At 20 µM, resveratrol had no effect on cellular proliferation or viability. However, it suppressed significantly 1) lateral motility by up to 25%; 2) transverse motility by 31%; and invasion by 37%. It also increased the cells' adhesion to substrate by 55%. Electrophysiologically, resveratrol inhibited voltage-gated Na(+) channel (VGSC) activity that has been shown previously to promote metastatic cell behaviors. This effect was dose-dependent with an IC50 of ∼50 µM. Voltage dependencies of current activation and peak were not affected but steady-state inactivation was shifted to more hyperpolarized potentials and recovery from inactivation was slowed. Coapplication of resveratrol with the highly specific VGSC blocker tetrodotoxin did not result in any additive effect on inhibition of both 1) VGSC activity and 2) metastatic cell behaviors. These results suggest 1) that a significant mode of action of resveratrol is VGSC blockage and 2) that resveratrol has promise as a natural antimetastatic agent.

Reducing annotation burden in MR: A novel MR-contrast guided contrastive learning approach for image segmentation.

BACKGROUND: Contrastive learning, a successful form of representational learning, has shown promising results in pretraining deep learning (DL) models for downstream tasks. When working with limited annotation data, as in medical image segmentation tasks, learning domain-specific local representations can further improve the performance of DL models. PURPOSE: In this work, we extend the contrastive learning framework to utilize domain-specific contrast information from unlabeled Magnetic Resonance (MR) images to improve the performance of downstream MR image segmentation tasks in the presence of limited labeled data. METHODS: The contrast in MR images is controlled by underlying tissue properties (e.g., T1 or T2) and image acquisition parameters. We hypothesize that learning to discriminate local representations based on underlying tissue properties should improve subsequent segmentation tasks on MR images. We propose a novel constrained contrastive learning (CCL) strategy that uses tissue-specific information via a constraint map to define positive and negative local neighborhoods for contrastive learning, embedding this information in the representational space during pretraining. For a given MR contrast image, the proposed strategy uses local signal characteristics (constraint map) across a set of related multi-contrast MR images as a surrogate for underlying tissue information. We demonstrate the utility of the approach for downstream: (1) multi-organ segmentation tasks in T2-weighted images where a DL model learns T2 information with constraint maps from a set of 2D multi-echo T2-weighted images (n = 101) and (2) tumor segmentation tasks in multi-parametric images from the public brain tumor segmentation (BraTS) (n = 80) dataset where DL models learn T1 and T2 information from multi-parametric BraTS images. Performance is evaluated on downstream multi-label segmentation tasks with limited data in (1) T2-weighted images of the abdomen from an in-house Radial-T2 (Train/Test = 30/20), (2) public Cartesian-T2 (Train/Test = 6/12) dataset, and (3) multi-parametric MR images from the public brain tumor segmentation dataset (BraTS) (Train/Test = 40/50). The performance of the proposed CCL strategy is compared to state-of-the-art self-supervised contrastive learning techniques. In each task, a model is also trained using all available labeled data for supervised baseline performance. RESULTS: The proposed CCL strategy consistently yielded improved Dice scores, Precision, and Recall metrics, and reduced HD95 values across all segmentation tasks. We also observed performance comparable to the baseline with reduced annotation effort. The t-SNE visualization of features for T2-weighted images demonstrates its ability to embed T2 information in the representational space. On the BraTS dataset, we also observed that using an appropriate multi-contrast space to learn T1+T2, T1, or T2 information during pretraining further improved the performance of tumor segmentation tasks. CONCLUSIONS: Learning to embed tissue-specific information that controls MR image contrast with the proposed constrained contrastive learning improved the performance of DL models on subsequent segmentation tasks compared to conventional self-supervised contrastive learning techniques. The use of such domain-specific local representations could help understand, improve performance, and mitigate the scarcity of labeled data in MR image segmentation tasks.

T2 relaxometry with indirect echo compensation from highly undersampled data.

PURPOSE: To develop an algorithm for fast and accurate T2 estimation from highly undersampled multi-echo spin-echo data. METHODS: The algorithm combines a model-based reconstruction with a signal decay based on the slice-resolved extended phase graph (SEPG) model with the goal of reconstructing T2 maps from highly undersampled radial multi-echo spin-echo data with indirect echo compensation. To avoid problems associated with the nonlinearity of the SEPG model, principal component decomposition is used to linearize the signal model. The proposed CUrve Reconstruction via principal component-based Linearization with Indirect Echo compensation (CURLIE) algorithm is used to estimate T2 curves from highly undersampled data. T2 maps are obtained by fitting the curves to the SEPG model. RESULTS: Results on phantoms showed T2 biases (1.9% to 18.4%) when indirect echoes are not taken into account. The T2 biases were reduced (< 3.2%) when the CURLIE reconstruction was performed along with SEPG fitting even for high degrees of undersampling (4% sampled). Experiments in vivo for brain, liver, and heart followed the same trend as the phantoms. CONCLUSION: The CURLIE reconstruction combined with SEPG fitting enables accurate T2 estimation from highly undersampled multi-echo spin-echo radial data thus, yielding a fast T2 mapping method without errors caused by indirect echoes.

T2 mapping from highly undersampled data by reconstruction of principal component coefficient maps using compressed sensing.

Recently, there has been an increased interest in quantitative MR parameters to improve diagnosis and treatment. Parameter mapping requires multiple images acquired with different timings usually resulting in long acquisition times. While acquisition time can be reduced by acquiring undersampled data, obtaining accurate estimates of parameters from undersampled data is a challenging problem, in particular for structures with high spatial frequency content. In this work, principal component analysis is combined with a model-based algorithm to reconstruct maps of selected principal component coefficients from highly undersampled radial MRI data. This novel approach linearizes the cost function of the optimization problem yielding a more accurate and reliable estimation of MR parameter maps. The proposed algorithm--reconstruction of principal component coefficient maps using compressed sensing--is demonstrated in phantoms and in vivo and compared with two other algorithms previously developed for undersampled data.

Convolutional Neural Network Based Frameworks for Fast Automatic Segmentation of Thalamic Nuclei from Native and Synthesized Contrast Structural MRI.

Thalamic nuclei have been implicated in several neurological diseases. Thalamic nuclei parcellation from structural MRI is challenging due to poor intra-thalamic nuclear contrast while methods based on diffusion and functional MRI are affected by limited spatial resolution and image distortion. Existing multi-atlas based techniques are often computationally intensive and time-consuming. In this work, we propose a 3D convolutional neural network (CNN) based framework for thalamic nuclei parcellation using T1-weighted Magnetization Prepared Rapid Gradient Echo (MPRAGE) images. Transformation of images to an efficient representation has been proposed to improve the performance of subsequent classification tasks especially when working with limited labeled data. We investigate this by transforming the MPRAGE images to White-Matter-nulled MPRAGE (WMn-MPRAGE) contrast, previously shown to exhibit good intra-thalamic nuclear contrast, prior to the segmentation step. We trained two 3D segmentation frameworks using MPRAGE images (n = 35 subjects): (a) a native contrast segmentation (NCS) on MPRAGE images and (b) a synthesized contrast segmentation (SCS) where synthesized WMn-MPRAGE representation generated by a contrast synthesis CNN were used. Thalamic nuclei labels were generated using THOMAS, a multi-atlas segmentation technique proposed for WMn-MPRAGE images. The segmentation accuracy and clinical utility were evaluated on a healthy cohort (n = 12) and a cohort (n = 45) comprising of healthy subjects and patients with alcohol use disorder (AUD), respectively. Both the segmentation CNNs yielded comparable performances on most thalamic nuclei with Dice scores greater than 0.84 for larger nuclei and at least 0.7 for smaller nuclei. However, for some nuclei, the SCS CNN yielded significant improvements in Dice scores (medial geniculate nucleus, P = 0.003, centromedian nucleus, P = 0.01) and percent volume difference (ventral anterior, P = 0.001, ventral posterior lateral, P = 0.01) over NCS. In the AUD cohort, the SCS CNN demonstrated a significant atrophy in ventral lateral posterior nucleus in AUD patients compared to healthy age-matched controls (P = 0.01), agreeing with previous studies on thalamic atrophy in alcoholism, whereas the NCS CNN showed spurious atrophy of the ventral posterior lateral nucleus. CNN-based segmentation of thalamic nuclei provides a fast and automated technique for thalamic nuclei prediction in MPRAGE images. The transformation of images to an efficient representation, such as WMn-MPRAGE, can provide further improvements in segmentation performance.

Improving subspace constrained radial fast spin echo MRI using block matching driven non-local low rank regularization.

Subspace-constrained reconstruction methods restrict the relaxation signals (of size M) in the scene to a pre-determined subspace (of size K≪M) and allow multi-contrast imaging and parameter mapping from accelerated acquisitions. However, these constraints yield poor image quality at some imaging contrasts, which can impact the parameter mapping performance. Additional regularization such as the use of joint-sparse (JS) or locally-low-rank (LLR) constraints can help improve the recovery of these images but are not sufficient when operating at high acceleration rates. We propose a method, non-local rank 3D (NLR3D), that is built on block matching and transform domain low rank constraints to allow high quality recovery of subspace-coefficient images (SCI) and subsequent multi-contrast imaging and parameter mapping. The performance of NLR3D was evaluated using Monte-Carlo (MC) simulations and compared against the JS and LLR methods. In vivo T 2 mapping results are presented on brain and knee datasets. MC results demonstrate improved bias, variance, and MSE behavior in both the multi-contrast images and parameter maps when compared to the JS and LLR methods. In vivo brain and knee results at moderate and high acceleration rates demonstrate improved recovery of high SNR early TE images as well as parameter maps. No significant difference was found in the T2 values measured in ROIs between the NLR3D reconstructions and the reference images (Wilcoxon signed rank test). The proposed method, NLR3D, enables recovery of high-quality SCI and, consequently, the associated multi-contrast images and parameter maps.

Rapid high-resolution volumetric T1 mapping using a highly accelerated stack-of-stars Look Locker technique.

PURPOSE: To develop a fast volumetric T1 mapping technique. MATERIALS AND METHODS: A stack-of-stars (SOS) Look Locker technique based on the acquisition of undersampled radial data (>30× relative to Nyquist) and an efficient multi-slab excitation scheme is presented. A principal-component based reconstruction is used to reconstruct T1 maps. Computer simulations were performed to determine the best choice of partitions per slab and degree of undersampling. The technique was validated in phantoms against reference T1 values measured with a 2D Cartesian inversion-recovery spin-echo technique. The SOS Look Locker technique was tested in brain (n = 4) and prostate (n = 5). Brain T1 mapping was carried out with and without kz acceleration and results between the two approaches were compared. Prostate T1 mapping was compared to standard techniques. A reproducibility study was conducted in brain and prostate. Statistical analyses were performed using linear regression and Bland Altman analysis. RESULTS: Phantom T1 values showed excellent correlations between SOS Look Locker and the inversion-recovery spin-echo reference (r2 = 0.9965; p < 0.0001) and between SOS Look Locker with slab-selective and non-slab selective inversion pulses (r2 = 0.9999; p < 0.0001). In vivo results showed that full brain T1 mapping (1 mm3) with kz acceleration is achieved in 4 min 21 s. Full prostate T1 mapping (0.9 × 0.9 × 4 mm3) is achieved in 2 min 43 s. T1 values for brain and prostate were in agreement with literature values. A reproducibility study showed coefficients of variation in the range of 0.18-0.2% (brain) and 0.15-0.18% (prostate). CONCLUSION: A rapid volumetric T1 mapping technique was developed. The technique enables high-resolution T1 mapping with adequate anatomical coverage in a clinically acceptable time.