Processing of radial fast spin-echo data for obtaining T2 estimates from a single k-space data set

Concurrent water T2 and fat fraction mapping of the breast using the radial gradient and spin echo (RADGRASE) pulse sequence

This work describes and evaluates an efficient radial gradient- and spin-echo (RADGRASE) pulse sequence and reconstruction algorithm for concurrent measurement of proton-density weighted fat fraction (FF) and water component T2 (T2w) within breast tissues. The ability to estimate T2w in breast tissues, where fat can be highly abundant, is demonstrated using oil/gel phantoms across a wide range of FF values (0.1-0.7). Successful T2w mapping of breast tissues is also demonstrated in vivo by comparison with fat suppressed T2 values. The sensitivity of RADGRASE to detect changes in the breast was assessed by tracking T2w in 3 healthy volunteers through their menstrual cycle, demonstrating T2w values in the late luteal phase to be 18-29 ms higher than in the follicular phase. The technique is also applied to a cohort of 68 patients taking tamoxifen for breast cancer risk reduction, where significant positive correlation between the FF parameter Frac50 and T2w (p = 0.035) was observed in premenopausal subjects (n = 20). Our findings demonstrate the ability and efficacy of RADGRASE for simultaneously mapping FF and T2w within breast tissues, and the potential utility of the technique in studying breast tissue changes in clinical applications.

Whole-brain T2 mapping with radial sampling and retrospective motion correction at 3T

PURPOSE

Several barriers prevent the use of whole-brain T2 mapping in routine use despite increasing interest in this parameter. One of the main barriers is the long scan time resulting in patient discomfort and motion corrupted data. To address this challenge, a method for accurate whole-brain T2 mapping with a limited acquisition time and motion correction capabilities is investigated.

METHODS

A 3D radial multi-echo spin-echo sequence was implemented with optimized sampling trajectory enabling the estimation of intra-scan motion, subsequently used to correct the raw data. Motion corrected echo images are then reconstructed with linear subspace constrained reconstruction. Experiments were carried out on phantom and volunteers at 3T to evaluate the accuracy of the T2 estimation, the sensitivity to lesions and the efficiency of the correction on motion corrupted data.

RESULTS

Whole-brain T2 mapping acquired in less than 7 min enabled the depiction of lesions in the white matter with longer T2. Data retrospectively corrupted with typical motion traces of pediatric patients highly benefited from the motion correction by reducing the error in T2 estimates within the lesions. All datasets acquired on seven volunteers, with deliberate motion, also showed that motion corrupted T2 maps could be improved with the retrospective motion correction both at the voxel level and the structure level.

CONCLUSION

A whole-brain T2 mapping sequence with retrospective intra-scan motion correction and reasonable acquisition time is proposed. The method necessitates advanced iterative reconstruction strategies but no additional navigator, external device, or increased scan time is required.

Quantitative T2 Mapping of Acute Pancreatitis

BACKGROUND

Quantification of the T2 signal by means of T2 mapping in acute pancreatitis (AP) has the potential to quantify the parenchymal edema. Quantitative T2 mapping may overcome the limitations of previously reported scoring systems for reliable assessment of AP.

PURPOSE

To evaluate MR-derived pancreatic T2 mapping values in AP and correlate them with markers of disease severity.

STUDY TYPE

Prospective single-center study.

POPULATION

76 adults with AP (20-91 years, females/males: 39/37).

FIELD STRENGTH/SEQUENCE

Fat suppressed multiecho spin-echo prototype sequence to quantify T2 signal at 3T MRI.

ASSESSMENT

The severity of AP was assessed clinically, biologically, and by contrast-enhanced CT (CECT) performed 48-72 hours after symptom onset. MRI was then performed ≤24 hours after CT. Two readers blinded to any clinical information independently evaluated the T2 values by placing three regions of interest inside the pancreatic head, body, and tail on the T2 mapping MR sequence. Results were compared with corresponding CECT images as the standard and clinical severity parameters, using the length of hospital stay as our primary endpoint.

STATISTICAL TESTS

Continuous variables were compared using the Spearman's rank correlation coefficient, analysis of variance (ANOVA) or Student's t-test.

RESULTS

T2 values significantly correlated with the length of hospital stay (rs(74) = 0.29), CT severity index (CTSI) (rs(73) = 0.61; CTSI 0-3: 72 ± 14 msec, CTSI 4-10: 88 ± 15), intensive care unit (ICU) admission (t(2.77) = -3.41) and presence of organ failure (t(6.72) = -3.42), whereas the CTSI and Ranson score were not significantly related with ICU admission (CTSI: P = 0.24; Ranson score: P = 0.24) and organ failure (CTSI: P = 0.11; Ranson score P = 0.11).

CONCLUSION

T2 mapping correlates with AP severity parameters and is useful for assessing the severity of AP with higher sensitivity than the usual clinical and radiological scoring systems.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 2.

Feasibility and Reproducibility of T2 Mapping Compared with Diffusion-Weighted Imaging in Solid Renal Masses

Accurate prediction of renal mass subtypes, along with the WHO/ISUP grade and pathological T (pT) stage of clear cell renal cell carcinoma (ccRCC), is crucial for optimal decision making. Our study aimed to investigate the feasibility and reproducibility of motion-robust radial T2 mapping in differentiating lipid-poor angiomyolipoma (MFAML) from RCC and characterizing the WHO/ISUP grade and pT stage of ccRCC. Finally, 92 patients undergoing renal radial T2 mapping and ZOOMit DWI were recruited. The T2 values and apparent diffusion coefficient (ADC) were analyzed. Correlation coefficients were calculated between ADC and T2 values. Notably, ccRCC exhibited higher T2 and ADC values than MFAML (p < 0.05). T2 values were lower in the higher WHO/ISUP grade and pT stage of ccRCC (all p < 0.05). ADC showed no significant difference for pT stage (p = 0.056). T2 values revealed a higher area under the curve (AUC) in evaluating the WHO/ISUP grade compared to ADC (0.936 vs. 0.817, p = 0.027). T2 values moderately positively correlated with ADC (r = 0.675, p < 0.001). In conclusion, quantitative motion-robust radial T2 mapping is feasible for characterizing solid renal masses and could provide additional value for multiparametric imaging in predicting WHO/ISUP grade and pT stage of ccRCC.

Constrained alternating minimization for parameter mapping (CAMP)

OBJECTIVE

To improve image quality in highly accelerated parameter mapping by incorporating a linear constraint that relates consecutive images.

APPROACH

In multi-echo T1 or T2 mapping, scan time is often shortened by acquiring undersampled but complementary measures of k-space at each TE or TI. However, residual undersampling artifacts from the individual images can then degrade the quality of the final parameter maps. In this work, a new reconstruction method, dubbed Constrained Alternating Minimization for Parameter mapping (CAMP), is introduced. This method simultaneously extracts T2 or T1* maps in addition to an image for each TE or TI from accelerated datasets, leveraging the constraints of the decay to improve the reconstructed image quality. The model enforces exponential decay through a linear constraint, resulting in a biconvex objective function that lends itself to alternating minimization. The method was tested in four in vivo volunteer experiments and validated in phantom studies and healthy subjects, using T2 and T1 mapping, with accelerations of up to 12.

MAIN RESULTS

CAMP is demonstrated for accelerated radial and Cartesian acquisitions in T2 and T1 mapping. The method is even applied to generate an entire T2 weighted image series from a single TSE dataset, despite the blockwise k-space sampling at each echo time. Experimental undersampled phantom and in vivo results processed with CAMP exhibit reduced artifacts without introducing bias.

SIGNIFICANCE

For a wide array of applications, CAMP linearizes the model cost function without sacrificing model accuracy so that the well-conditioned and highly efficient reconstruction algorithm improves the image quality of accelerated parameter maps.

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.

The CALIPR framework for highly accelerated myelin water imaging with improved precision and sensitivity

Quantitative magnetic resonance imaging (MRI) techniques are powerful tools for the study of human tissue, but, in practice, their utility has been limited by lengthy acquisition times. Here, we introduce the Constrained, Adaptive, Low-dimensional, Intrinsically Precise Reconstruction (CALIPR) framework in the context of myelin water imaging (MWI); a quantitative MRI technique generally regarded as the most rigorous approach for noninvasive, in vivo measurement of myelin content. The CALIPR framework exploits data redundancy to recover high-quality images from a small fraction of an imaging dataset, which allowed MWI to be acquired with a previously unattainable sequence (fully sampled acquisition 2 hours:57 min:20 s) in 7 min:26 s (4.2% of the dataset, acceleration factor 23.9). CALIPR quantitative metrics had excellent precision (myelin water fraction mean coefficient of variation 3.2% for the brain and 3.0% for the spinal cord) and markedly increased sensitivity to demyelinating disease pathology compared to a current, widely used technique. The CALIPR framework facilitates drastically improved MWI and could be similarly transformative for other quantitative MRI applications.

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.

Comparison of T2 Quantification Strategies in the Abdominal-Pelvic Region for Clinical Use

OBJECTIVES

The aim of the study was to compare different magnetic resonance imaging (MRI) acquisition strategies appropriate for T2 quantification in the abdominal-pelvic area. The different techniques targeted in the study were chosen according to 2 main considerations: performing T2 measurement in an acceptable time for clinical use and preventing/correcting respiratory motion.

MATERIALS AND METHODS

Acquisitions were performed at 3 T. To select sequences for in vivo measurements, a phantom experiment was conducted, for which the T2 values obtained with the different techniques of interest were compared with the criterion standard (single-echo SE sequence, multiple acquisitions with varying echo time). Repeatability and temporal reproducibility studies for the different techniques were also conducted on the phantom. Finally, an in vivo study was conducted on 12 volunteers to compare the techniques that offer acceptable acquisition time for clinical use and either address or correct respiratory motion.

RESULTS

For the phantom study, the DESS and T2-preparation techniques presented the lowest precision (ρ2 = 0.9504 and ρ2 = 0.9849 respectively), and showed a poor repeatability/reproducibility compared with the other techniques. The strategy relying on SE-EPI showed the best precision and accuracy (ρ2 = 0.9994 and Cb = 0.9995). GRAPPATINI exhibited a very good precision (ρ2 = 0.9984). For the technique relying on radial TSE, the precision was not as good as GRAPPATINI (ρ2 = 0.9872). The in vivo study demonstrated good respiratory motion management for all of the selected techniques. It also showed that T2 estimate ranges were different from one method to another. For GRAPPATINI and radial TSE techniques, there were significant differences between all the different types of organs of interest.

CONCLUSIONS

To perform T2 measurement in the abdominal-pelvic region, one should favor a technique with acceptable acquisition time for clinical use, with proper respiratory motion management, with good repeatability, reproducibility, and precision. In this study, the techniques relying respectively on SE-EPI, radial TSE, and GRAPPATINI appeared as good candidates.

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.

Evaluation of a single-breath-hold radial turbo-spin-echo sequence for T2 mapping of the liver at 3T

OBJECTIVES

T2 mapping of the liver is a potential diagnostic tool, but conventional techniques are difficult to perform in clinical practice due to long scan time. We aimed to evaluate the accuracy of a prototype radial turbo-spin-echo (rTSE) sequence, optimized for multi-slice T2 mapping in the abdomen during one breath-hold at 3 T.

METHODS

A multi-sample (fat: 0-35%) agarose phantom doped with MnCl₂ and 80 subjects (73 patients undergoing abdomen MR examination and 7 healthy volunteers) were investigated. A radial turbo-spin-echo (rTSE) sequence with and without fat suppression, a Cartesian turbo-spin-echo (Cart-TSE) sequence, and a single-voxel multi-echo STEAM spectroscopy (HISTO) were performed in phantom, and fat-suppressed rTSE and HISTO sequences were performed in in vivo measurements. Two approaches were used to sample T2 values: manually selected circular ROIs and whole liver analysis with Gaussian mixture models (GMM).

RESULTS

The rTSE-T2s values exhibited a strong correlation with Cart-TSE-T2s (R² = 0.988) and with HISTO-T2s of water (R² = 0.972) in phantom with an offset between rTSE and Cart-TSE maps (mean difference = 3.17 ± 1.18 ms). The application of fat suppression decreased T2 values, and the effect was directly proportional to the amount of fat. Measurements in patients yielded a linear relationship between rTSE- and HISTO-T2s (R² = 0.546 and R² = 0.580 for ROI and GMM, respectively).

CONCLUSION

The fat-suppressed rTSE sequence allows for fast and accurate determination of T2 values of the liver, and appears to be suitable for further large cohort studies.

KEY POINTS

•Radial turbo-spin-echo T2 mapping performs comparably to Cartesian TSE-T2 mapping, but an offset in values is observed in phantom measurements. •Fat-suppressed radial turbo-spin-echo T2 mapping is consistent with T2 of water as assessed by MRS in phantom measurements. •Fat-suppressed radial turbo-spin-echo sequence allows fast T2 mapping of the liver in a single breath-hold and is correlated with MRS-based T2 of water.

T2 mapping in gadoxetic acid-enhanced MRI: utility for predicting decompensation and death in cirrhosis

OBJECTIVES

To determine whether T2 mapping in liver MRI can predict decompensation and death in cirrhotic patients.

METHODS

This retrospective study included 292 cirrhotic patients who underwent gadoxetic acid-enhanced MRI, including T1 and T2 mapping at 10-min hepatobiliary phase by using the Look-Locker and radial turbo spin-echo sequences, respectively. T1 and T2 values of the liver and spleen were measured. The association of MR parameters and serum markers with decompensation and death was investigated. Risk models combining T2_Liver, serum albumin level, and Model for End-Stage Liver Disease (MELD) score were created for predicting decompensation (T2_Liver, < 49.3 versus ≥ 49.3 ms) and death (< 57.4 versus ≥ 57.4 ms).

RESULTS

In patients with compensated cirrhosis at baseline and in the full patient cohort, 9.6% (19 of 197) and 5.1% (15 of 292) developed decompensation and died during the mean follow-up periods of 18.7 and 19.2 months, respectively. A prolonged T2_Liver (hazard ratio (HR), 2.59; 95% confidence interval (CI), 1.26, 5.31) was independently predictive of decompensation along with the serum albumin level (HR, 0.28; 95% CI, 0.12, 0.68) and MELD score (HR, 1.34; 95% CI, 1.08, 1.66). T2_Liver (HR, 2.61; 95% CI, 1.19, 5.72) and serum albumin level (HR, 0.46; 95% CI, 0.19, 1.14) were independent predictors of death. The mean times to decompensation (12.9 versus 29.2 months) and death (16.5 versus 29.6 months) were significantly different between the high- and low-risk groups (p < 0.001).

CONCLUSION

T2_Liver from T2 mapping can predict decompensation and death in patients with cirrhosis.

KEY POINTS

• Liver T2 values from the radial turbo spin-echo (TSE) T2 mapping sequence with tiered echo sharing and pseudo golden-angle (pGA) reordering were significantly higher in decompensated cirrhosis than compensated cirrhosis. • Liver T2 values from the radial TSE T2 mapping sequence with tiered echo sharing and pGA reordering can predict decompensation and death in patients with cirrhosis. • T2 mapping is recommended as part of liver MRI examinations for cirrhotic patients because it can provide a noninvasive prognostic marker for the development of decompensation and death.

Evaluation of feasibility and image quality of a new radial quantitative T2 weighted imaging sequence for liver MRI

OBJECTIVES

To evaluate the clinical feasibility of a new T2 weighted sequence to calculate T2 relaxation times (T2RT) of liver lesions using two-dimensional radial turbo spin echo (2DRTSE) and to evaluate this sequence by performing image quality and relaxation time comparison of multiple liver lesions.

MATERIALS AND METHODS

This prospective analysis of 2DRTSE sequences (using 22 echoes) was performed in 19 patients with 36 liver lesions. Two radiologists independently obtained T2RTs for liver lesions and scored image quality and image artifacts. Lesions were classified as cyst, hemangioma, solid, or necrotic. T2RT values were compared. Inter-reader agreement was evaluated.

RESULTS

The 2DRTSE images were considered good quality with few artifacts by both radiologists. Nineteen patients were included in the study, with a total of 36 liver lesions. Two of the liver lesions were classified as cysts, 7 as hemangiomas, 4 as necrotic lesions, and 23 as solid lesions. The concordance correlation coefficient was 0.996 for the calculated T2RT of each liver lesion between the two readers, indicating good agreement. There was statically significant difference of the calculated T2RT for each lesion type.

CONCLUSION

The 2DRTSE sequence can be performed and provides good T2W image quality and a quantitative T2RT map of the entire abdomen. The liver lesions can be distinguished based on the calculated T2RT using this technique. 2DRTSE could potentially supplant the current T2-weighted imaging sequence with the benefit of quantitative T2RTs.

T1 mapping, T2 mapping and MR elastography of the liver for detection and staging of liver fibrosis

PURPOSE

To compare liver stiffness measurements obtained from MR elastography with liver T1 relaxation times obtained from T1 mapping and T2 relaxation times obtained from T2 mapping for detection and staging of liver fibrosis.

MATERIALS AND METHODS

223 patients with known or suspected liver disease underwent MRI of the liver with T1 mapping (Look-Locker sequence) and 2D SE-EPI MR elastography (MRE) sequences. 139 of these patients also underwent T2 mapping with radial T2 TSE sequence. Two readers (R1 & R2) measured liver stiffness, T1 relaxation times and T2 relaxation times. T1 and T2 times were correlated with stiffness measurements. ROC analysis was used to compare the performance of both techniques in discriminating fibrosis stage in 23 patients who underwent liver biopsy.

RESULTS

For each reader there was significant moderate positive correlation between liver MRE and liver T1 mapping (r = 0.49 and 0.36). There was significant moderate positive correlation between liver T2 mapping and each of MRE and T1 mapping for one of the readers (r = 0.40 and 0.27). AUC for differentiating early (F0-F2) from advanced (F3-F4) fibrosis in biopsied patients was 0.975 (R1) and 0.925 (R2) for MRE, 0.671 (R1) and 0.642 (R2) for T1 mapping and 0.671 (R1) and 0.743 (R2) for T2 mapping. Inter-reader agreement was good for MRE (ICC = 0.84) substantial for T1 mapping (0.94) and T2 mapping (0.96).

CONCLUSIONS

Liver T1 and T2 mapping showed moderate positive correlation with MR elastography. Accuracy of MRE is however superior to T1 and T2 mapping in the subset of patients who underwent liver biopsy. Accuracy of combination of MRE and T1 mapping/T2 mapping was not superior to MRE alone.

Quantification of MRI T2-weighted High Signal Volume in Cystic Fibrosis: A Pilot Study

Background In patients with cystic fibrosis (CF), pulmonary structures with high MRI T2 signal intensity relate to inflammatory changes in the lung and bronchi. These areas of pathologic abnormalities can serve as imaging biomarkers. The feasibility of automated quantification is unknown. Purpose To quantify the MRI T2 high-signal-intensity lung volume and T2-weighted volume-intensity product (VIP) by using a black-blood T2-weighted radial fast spin-echo sequence in participants with CF. Materials and Methods Healthy individuals and study participants with CF were prospectively enrolled between January 2017 and November 2017. All participants underwent a lung MRI protocol including T2-weighted radial fast spin-echo sequence. Participants with CF also underwent pulmonary function tests the same day. Participants with CF exacerbation underwent repeat MRI after their treatment with antibiotics. Two observers supervised automated quantification of T2-weighted high-signal-intensity volume (HSV) and T2-weighted VIP independently, and the average score was chosen as consensus. Statistical analysis used the Mann-Whitney test for comparison of medians, correlations used the Spearman test, comparison of paired medians used the Wilcoxon signed rank test, and reproducibility was evaluated by using intraclass correlation coefficient. Results In 10 healthy study participants (median age, 21 years [age range, 18-27 years]; six men) and 12 participants with CF (median age, 18 years [age range, 9-40 years]; eight men), T2-weighted HSV was equal to 0% and 4.1% (range, 0.1%-17%), respectively, and T2-weighted VIP was equal to 0 msec and 303 msec (range, 39-1012 msec), respectively (P < .001). In participants with CF, T2-weighted HSV or T2-weighted VIP were associated with forced expiratory volume in 1 second percentage predicted (ρ = -0.88 and ρ = -0.94, respectively; P < .001). In six participants with CF exacerbation and follow-up after treatment, a decrease in both T2-weighted HSV and T2-weighted VIP was observed (P = .03). The intra- and interobserver reproducibility of MRI were good (intraclass correlation coefficients, >0.99 and >0.99, respectively). Conclusion In patients with cystic fibrosis (CF), automated quantification of lung MRI high-signal-intensity volume was reproducible and correlated with pulmonary function testing severity, and it improved after treatment for CF exacerbation. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Revel and Chassagnon in this issue.

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.

Accelerated and high-resolution cardiac T2 mapping through peripheral k-space sharing

PURPOSE

To develop high-spatial-resolution cardiac T₂ mapping that allows for a reduced acquisition time while maintaining its precision. We implemented and optimized a new golden-angle radial T₂ mapping technique named SKRATCH (Shared k-space Radial T₂ Characterization of the Heart) that shares k-space peripheries of T₂ -weighted images while preserving their contrasts.

METHODS

Six SKRATCH variants (gradient-recalled echo and balanced SSFP, free-breathing and breath-held, with and without a saturation preparation) were implemented, and their precision was compared with a navigator-gated reference technique in phantoms and 22 healthy volunteers at 3 T. The optimal breath-held SKRATCH technique was applied in a small cohort of patients with subacute myocardial infarction.

RESULTS

The faster free-breathing SKRATCH technique reduced the acquisition time by 52.4%, while maintaining the precision and spatial resolution of the reference technique. Similarly, the most precise and robust breath-held SKRATCH technique demonstrated homogenous T₂ values that did not significantly differ from the navigator-gated reference (T₂  = 39.9 ± 3.4 ms versus 39.5 ± 3.4 ms, P > .20, respectively). All infarct patients demonstrated a large T₂ elevation in the ischemic regions of the myocardium.

CONCLUSION

The optimized SKRATCH technique enabled the accelerated acquisition of high-spatial-resolution T₂ maps, was validated in healthy adult volunteers, and was successfully applied to a small initial group of patients.

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 T₁ and T₂ 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 T₂ 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.

Ultrafast compartmentalized relaxation time mapping with linear algebraic modeling

PURPOSE

To dramatically accelerate compartmental-average longitudinal (T₁ ) and transverse (T₂ ) relaxation measurements using the minimal-acquisition linear algebraic modeling (SLAM) method, and to validate it in phantoms and humans.

METHODS

Relaxation times were imaged at 3 Tesla in phantoms, in the abdomens of six volunteers, and in six brain tumor patients using standard inversion recovery and multi-spin-echo sequences. k-space was fully sampled to provide reference T₁ and T₂ measurements, and SLAM was performed using a limited set of phase encodes from central k-space. Anatomical compartments were segmented on scout images post-acquisition, and SLAM reconstruction was implemented using two algorithms. Compartment-average T₁ and T₂ measurements were determined retroactively from fully sampled data sets, and proactively from SLAM data sets at acceleration factors of up to 16. Values were compared with reference measurements. The compartment's localization properties were analyzed using the discrete spatial response function.

RESULTS

At 16-fold acceleration, compartment-average SLAM T₁ measurements agreed with the full k-space compartment-average results to within 0.0% ± 0.7%, 1.4% ± 3.4%, and 0.5% ± 2.9% for phantom, abdominal, and brain T₁ measurements, respectively. The corresponding T₂ measurements agreed within 0.2% ± 1.9%, 0.9% ± 7.9%, and 0.4% ± 5.8%, respectively.

CONCLUSION

SLAM can dramatically accelerate relaxation time measurements when compartmental or lesion-average values can suffice, or when standard relaxometry is precluded by scan-time limitations. Magn Reson Med 79:286-297, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Experimental O-space turbo spin echo imaging

PURPOSE

Turbo spin echo (TSE) imaging reduces imaging time by acquiring multiple echoes per repetition (TR), requiring fewer TRs. O-space can also require fewer TRs by using a combination of nonlinear magnetic gradient fields and surface coil arrays. Although to date, O-space has only been demonstrated for gradient echo imaging, it is valuable to combine these two techniques. However, collecting multiple O-space echoes per TR is difficult because of the different local k-space trajectories and variable T2-weighting.

THEORY AND METHODS

A practical scheme is demonstrated to combine the benefits of TSE and O-space for highly accelerated T2-weighted images. The scheme uses a modified acquisition order and filtered projection reconstruction to reduce artifacts caused by T2 decay, while retaining T2 contrast that corresponds to a specific echo time.

RESULTS

The experiments revealed that the proposed method can produce highly accelerated T2-weighted images. Moreover, the method can generate multiple images with different T2 contrasts from a single dataset.

CONCLUSIONS

The proposed O-space TSE imaging method requires fewer echoes than conventional TSE and fewer repetitions than conventional O-space imaging. It retains resilience to undersampling, clearly outperforming Cartesian SENSE at high levels of undersampling, and can generate undistorted images with a range of T2 contrast from a single acquired dataset.

Accelerated and motion-robust in vivo T2 mapping from radially undersampled data using bloch-simulation-based iterative reconstruction

PURPOSE

Development of a quantitative transverse relaxation time (T2)-mapping platform that operates at clinically feasible timescales by employing advanced image reconstruction of radially undersampled multi spin-echo (MSE) datasets.

METHODS

Data was acquired on phantom and in vivo at 3 Tesla using MSE protocols employing radial k-space sampling trajectories. In order to overcome the nontrivial spin evolution associated with MSE protocols, a numerical signal model was precalculated based on Bloch simulations of the actual pulse-sequence scheme used in the acquisition process. This signal model was subsequently incorporated into an iterative model-based image reconstruction process, producing T2 and proton-density maps.

RESULTS

T2 maps of phantom and in vivo brain were successfully constructed, closely matching values produced by a single spin-echo reference scan. High-resolution mapping was also performed for the spinal cord in vivo, differentiating the underlying gray/white matter morphology.

CONCLUSION

The presented MSE data-processing framework offers reliable mapping of T2 relaxation values in a ∼ 5-minute timescale, free of user- and scanner-dependent variations. The use of radial k-space sampling provides further advantages in the form of high immunity to irregular physiological motion, as well as enhanced spatial resolutions, owing to its inherent ability to perform alias-free limited field-of-view imaging.

T2 mapping of the heart with a double-inversion radial fast spin-echo method with indirect echo compensation

BACKGROUND

The abnormal signal intensity in cardiac T2-weighted images is associated with various pathologies including myocardial edema. However, the assessment of pathologies based on signal intensity is affected by the acquisition parameters and the sensitivities of the receiver coils. T2 mapping has been proposed to overcome limitations of T2-weighted imaging, but most methods are limited in spatial and/or temporal resolution. Here we present and evaluate a double inversion recovery radial fast spin-echo (DIR-RADFSE) technique that yields data with high spatiotemporal resolution for cardiac T2 mapping.

METHODS

DIR-RADFSE data were collected at 1.5 T on phantoms and subjects with echo train length (ETL) = 16, receiver bandwidth (BW) = ±32 kHz, TR = 1RR, matrix size = 256 × 256. Since only 16 views per echo time (TE) are collected, two algorithms designed to reconstruct highly undersampled radial data were used to generate images for 16 time points: the Echo-Sharing (ES) and the CUrve Reconstruction via pca-based Linearization with Indirect Echo compensation (CURLIE) algorithm. T2 maps were generated via least-squares fitting or the Slice-resolved Extended Phase Graph (SEPG) model fitting. The CURLIE-SEPG algorithm accounts for the effect of indirect echoes. The algorithms were compared based on reproducibility, using Bland-Altman analysis on data from 7 healthy volunteers, and T2 accuracy (against a single-echo spin-echo technique) using phantoms.

RESULTS

Both reconstruction algorithms generated in vivo images with high spatiotemporal resolution and showed good reproducibility. Mean T2 difference between repeated measures and the coefficient of repeatability were 0.58 ms and 2.97 for ES and 0.09 ms and 4.85 for CURLIE-SEPG. In vivo T2 estimates from ES were higher than those from CURLIE-SEPG. In phantoms, CURLIE-SEPG yielded more accurate T2s compared to reference values (error was 7.5-13.9% for ES and 0.6-2.1% for CURLIE-SEPG), consistent with the fact that CURLIE-SEPG compensates for the effects of indirect echoes. The potential of T2 mapping with CURLIE-SEPG is demonstrated in two subjects with known heart disease. Elevated T2 values were observed in areas of suspected pathology.

CONCLUSIONS

DIR-RADFSE yielded TE images with high spatiotemporal resolution. Two algorithms for generating T2 maps from highly undersampled data were evaluated in terms of accuracy and reproducibility. Results showed that CURLIE-SEPG yields T2 estimates that are reproducible and more accurate than ES.

Reducing contrast contamination in radial turbo-spin-echo acquisitions by combining a narrow-band KWIC filter with parallel imaging

PURPOSE

Cartesian turbo spin-echo (TSE) and radial TSE images are usually reconstructed by assembling data containing different contrast information into a single k-space. This approach results in mixed contrast contributions in the images, which may reduce their diagnostic value. The goal of this work is to improve the image contrast from radial TSE acquisitions by reducing the contribution of signals with undesired contrast information.

METHODS

Radial TSE acquisitions allow the reconstruction of multiple images with different T2 contrasts using the k-space weighted image contrast (KWIC) filter. In this work, the image contrast is improved by reducing the band-width of the KWIC filter. Data for the reconstruction of a single image are selected from within a small temporal range around the desired echo time. The resulting dataset is undersampled and, therefore, an iterative parallel imaging algorithm is applied to remove aliasing artifacts.

RESULTS

Radial TSE images of the human brain reconstructed with the proposed method show an improved contrast when compared with Cartesian TSE images or radial TSE images with conventional KWIC reconstructions.

CONCLUSION

The proposed method provides multi-contrast images from radial TSE data with contrasts similar to multi spin-echo images. Contaminations from unwanted contrast weightings are strongly reduced.

Quantifying iron-oxide nanoparticles at high concentration based on longitudinal relaxation using a three-dimensional SWIFT Look-Locker sequence

PURPOSE

Iron-oxide nanoparticles (IONPs) have proven utility as contrast agents in many MRI applications. Previous quantitative IONP mapping has been performed using mainly T2 * mapping methods. However, in applications requiring high IONP concentrations, such as magnetic nanoparticles based thermal therapies, conventional pulse sequences are unable to map T2 * because the signal decays too rapidly. In this article, sweep imaging with Fourier transformation (SWIFT) sequence is combined with the Look-Locker method to map T1 of IONPs in high concentrations.

METHODS

T1 values of agar containing IONPs in different concentrations were measured with the SWIFT Look-Locker method and with inversion recovery spectroscopy. Precisions of Look-Locker and variable flip angle (VFA) methods were compared in simulations.

RESULTS

The measured R1 (=1/T1 ) has a linear relationship with IONP concentration up to 53.6 mM of Fe. This concentration exceeds concentrations measured in previous work by almost an order of magnitude. Simulations show SWIFT Look-Locker method is also much less sensitive to B1 inhomogeneity than the VFA method.

CONCLUSION

SWIFT Look-Locker can accurately measure T1 of IONP concentrations ≤53.6 mM. By mapping T1 as a function of IONP concentration, IONP distribution maps might be used in the future to plan effective magnetic nanoparticle hyperthermia therapy.

A comparison of radial keyhole strategies for high spatial and temporal resolution 4D contrast-enhanced MRI in small animal tumor models

PURPOSE

Dynamic contrast-enhanced (DCE) MRI has been widely used as a quantitative imaging method for monitoring tumor response to therapy. The simultaneous challenges of increasing temporal and spatial resolution in a setting where the signal from the much smaller voxel is weaker have made this MR technique difficult to implement in small-animal imaging. Existing protocols employed in preclinical DCE-MRI acquire a limited number of slices resulting in potentially lost information in the third dimension. This study describes and compares a family of four-dimensional (3D spatial + time), projection acquisition, radial keyhole-sampling strategies that support high spatial and temporal resolution.

METHODS

The 4D method is based on a RF-spoiled, steady-state, gradient-recalled sequence with minimal echo time. An interleaved 3D radial trajectory with a quasi-uniform distribution of points in k-space was used for sampling temporally resolved datasets. These volumes were reconstructed with three different k-space filters encompassing a range of possible radial keyhole strategies. The effect of k-space filtering on spatial and temporal resolution was studied in a 5 mM CuSO(4) phantom consisting of a meshgrid with 350-μm spacing and in 12 tumors from three cell lines (HT-29, LoVo, MX-1) and a primary mouse sarcoma model (three tumors∕group). The time-to-peak signal intensity was used to assess the effect of the reconstruction filters on temporal resolution. As a measure of heterogeneity in the third dimension, the authors analyzed the spatial distribution of the rate of transport (K(trans)) of the contrast agent across the endothelium barrier for several different types of tumors.

RESULTS

Four-dimensional radial keyhole imaging does not degrade the system spatial resolution. Phantom studies indicate there is a maximum 40% decrease in signal-to-noise ratio as compared to a fully sampled dataset. T1 measurements obtained with the interleaved radial technique do not differ significantly from those made with a conventional Cartesian spin-echo sequence. A bin-by-bin comparison of the distribution of the time-to-peak parameter shows that 4D radial keyhole reconstruction does not cause significant temporal blurring when a temporal resolution of 9.9 s is used for the subsamples of the keyhole data. In vivo studies reveal substantial tumor heterogeneity in the third spatial dimension that may be missed with lower resolution imaging protocols.

CONCLUSIONS

Volumetric keyhole imaging with projection acquisition provides a means to increase spatiotemporal resolution and coverage over that provided by existing 2D Cartesian protocols. Furthermore, there is no difference in temporal resolution between the higher spatial resolution keyhole reconstruction and the undersampled projection data. The technique allows one to measure complex heterogeneity of kinetic parameters with isotropic, microscopic spatial resolution.

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.

A novel mouse model of Niemann-Pick type C disease carrying a D1005G-Npc1 mutation comparable to commonly observed human mutations

We have identified a point mutation in Npc1 that creates a novel mouse model (Npc1(nmf164)) of Niemann-Pick type C1 (NPC) disease: a single nucleotide change (A to G at cDNA bp 3163) that results in an aspartate to glycine change at position 1005 (D1005G). This change is in the cysteine-rich luminal loop of the NPC1 protein and is highly similar to commonly occurring human mutations. Genetic and molecular biological analyses, including sequencing the Npc1(spm) allele and identifying a truncating mutation, confirm that the mutation in Npc1(nmf164) mice is distinct from those in other existing mouse models of NPC disease (Npc1(nih), Npc1(spm)). Analyses of lifespan, body and spleen weight, gait and other motor activities, as well as acoustic startle responses all reveal a more slowly developing phenotype in Npc1(nmf164) mutant mice than in mice with the null mutations (Npc1(nih), Npc1(spm)). Although Npc1 mRNA levels appear relatively normal, Npc1(nmf164) brain and liver display dramatic reductions in Npc1 protein, as well as abnormal cholesterol metabolism and altered glycolipid expression. Furthermore, histological analyses of liver, spleen, hippocampus, cortex and cerebellum reveal abnormal cholesterol accumulation, glial activation and Purkinje cell loss at a slower rate than in the Npc1(nih) mouse model. Magnetic resonance imaging studies also reveal significantly less demyelination/dysmyelination than in the null alleles. Thus, although prior mouse models may correspond to the severe infantile onset forms of NPC disease, Npc1(nmf164) mice offer many advantages as a model for the late-onset, more slowly progressing forms of NPC disease that comprise the large majority of human cases.

In vivo assessment of neurodegeneration in Niemann-Pick type C mice by quantitative T2 mapping and diffusion tensor imaging

PURPOSE

To quantitatively and noninvasively assess neurological disease progression in a mouse model of Niemann-Pick type C (NPC) disease by measuring white matter status with magnetic resonance imaging (MRI) techniques of T2 mapping and diffusion tensor imaging (DTI).

MATERIALS AND METHODS

Quantitative T2 and DTI experiments were performed in vivo in NPC disease model and control mice at three timepoints to quantify differences and changes in white matter with measurements of T2 relaxation and DTI parameters. Histological staining for myelin content was also performed at two timepoints to compare with the MRI findings.

RESULTS

The results of the T2 and DTI measurements show significant differences in white matter areas of the brain in the NPC disease model compared to control mice at several timepoints, and were seen to change over time in both groups.

CONCLUSION

The findings of this study suggest that quantitative MRI measurements may be suitable in vivo biomarkers of disease status for future studies of NPC disease models. The changes in white matter measurements between timepoints in both control and NPC disease groups suggest that white matter structures continue to change and develop over time in the NPC model and can be tracked with MRI techniques.

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.

Model-based iterative reconstruction for radial fast spin-echo MRI

In radial fast spin-echo magnetic resonance imaging (MRI), a set of overlapping spokes with an inconsistent T2 weighting is acquired, which results in an averaged image contrast when employing conventional image reconstruction techniques. This work demonstrates that the problem may be overcome with the use of a dedicated reconstruction method that further allows for T2 quantification by extracting the embedded relaxation information. Thus, the proposed reconstruction method directly yields a spin-density and relaxivity map from only a single radial data set. The method is based on an inverse formulation of the problem and involves a modeling of the received MRI signal. Because the solution is found by numerical optimization, the approach exploits all data acquired. Further, it handles multicoil data and optionally allows for the incorporation of additional prior knowledge. Simulations and experimental results for a phantom and human brain in vivo demonstrate that the method yields spin-density and relaxivity maps that are neither affected by the typical artifacts from TE mixing, nor by streaking artifacts from the incomplete k-space coverage at individual echo times.

Rapid water and lipid imaging with T2 mapping using a radial IDEAL-GRASE technique

Three-point Dixon methods have been investigated as a means to generate water and fat images without the effects of field inhomogeneities. Recently, an iterative algorithm (IDEAL, iterative decomposition of water and fat with echo asymmetry and least squares estimation) was combined with a gradient and spin-echo acquisition strategy (IDEAL-GRASE) to provide a time-efficient method for lipid-water imaging with correction for the effects of field inhomogeneities. The method presented in this work combines IDEAL-GRASE with radial data acquisition. Radial data sampling offers robustness to motion over Cartesian trajectories as well as the possibility of generating high-resolution T(2) maps in addition to the water and fat images. The radial IDEAL-GRASE technique is demonstrated in phantoms and in vivo for various applications including abdominal, pelvic, and cardiac imaging.

Cardiac MRI is Complementary to Echocardiography in the Assessment of Cardiac Masses

Despite the fact that the incidence of cardiac tumors is low, the prompt evaluation and adequate intervention of these is highly important. Although most tumors of the heart are considered histologically benign, there are significant risks associated with these "benign" tumors. These are associated with significant morbidity and mortality due to obstruction of blood flow, alterations of conduction, propagation of arrhythmias, and thromboembolism, depending on their size, location, and nature. With the advent of noninvasive imaging modalities--traditionally echocardiography; but more recently using cross-sectional imaging with cardiac computed tomography and magnetic resonance imaging--cardiac tumors can be optimally assessed providing a greater opportunity for curative treatments by cardiothoracic surgery.

An optimal radial profile order based on the Golden Ratio for time-resolved MRI

In dynamic magnetic resonance imaging (MRI) studies, the motion kinetics or the contrast variability are often hard to predict, hampering an appropriate choice of the image update rate or the temporal resolution. A constant azimuthal profile spacing (111.246 degrees), based on the Golden Ratio, is investigated as optimal for image reconstruction from an arbitrary number of profiles in radial MRI. The profile order is evaluated and compared with a uniform profile distribution in terms of signal-to-noise ratio (SNR) and artifact level. The favorable characteristics of such a profile order are exemplified in two applications on healthy volunteers. First, an advanced sliding window reconstruction scheme is applied to dynamic cardiac imaging, with a reconstruction window that can be flexibly adjusted according to the extent of cardiac motion that is acceptable. Second, a contrast-enhancing k-space filter is presented that permits reconstructing an arbitrary number of images at arbitrary time points from one raw data set. The filter was utilized to depict the T1-relaxation in the brain after a single inversion prepulse. While a uniform profile distribution with a constant angle increment is optimal for a fixed and predetermined number of profiles, a profile distribution based on the Golden Ratio proved to be an appropriate solution for an arbitrary number of profiles.

Inversion recovery radial MRI with interleaved projection sets

The radial trajectory has found applications in cardiac imaging because of its resilience to undersampling and motion artifacts. Recent work has shown that interleaved and weighted radial imaging can produce images with multiple contrasts from a single data set. This feature was investigated for inversion recovery imaging of scar using a radial technique. The 2D radial imaging method was modified to acquire quadruply interleaved projection sets within each acquisition window of the cardiac cycle. These data were reconstructed using k-space weightings that used a smaller segment of the acquisition window for the central k-space data, the determinant of image contrast. This method generates four images with different T1 weightings. The novel approach was compared with noninterleaved radial imaging, interleaved radial without weightings, and Cartesian imaging in simulations, phantoms, and seven subjects with clinical myocardial infarction. The results show that during a typical acquisition window after an inversion pulse, magnetization changes rapidly. The interleaved acquisition provided better image quality than the noninterleaved radial acquisition. Interleaving with weighting provided better quality when the inversion time (TI) was shorter than optimal; otherwise, interleaving without weighting was superior. These methods enable a radial trajectory to be employed in conjunction with preparation pulses for viability imaging.