Search

PURPOSE: Accurate photon attenuation assessment from MR data remains an unmet challenge in the thorax due to tissue heterogeneity and the difficulty of MR lung imaging. As thoracic tissues encompass the whole physiologic range of photon absorption, large errors can occur when using, for example, a uniform, water‐equivalent or a soft‐tissue‐only approximation. The purpose of this study was to introduce a method for voxel‐wise thoracic synthetic CT (sCT) generation from MR data attenuation correction (AC) for PET/MR or for MR‐only radiation treatment planning (RTP). METHODS: Acquisition: A radial stack‐of‐stars combining ultra‐short‐echo time (UTE) and modified Dixon (mDixon) sequence was optimized for thoracic imaging. The UTE‐mDixon pulse sequence collects MR signals at three TE times denoted as UTE, Echo1, and Echo2. Three‐point mDixon processing was used to reconstruct water and fat images. Bias field correction was applied in order to avoid artifacts caused by inhomogeneity of the MR magnetic field. Analysis: Water fraction and R2* maps were estimated using the UTE‐mDixon data to produce a total of seven MR features, that is UTE, Echo1, Echo2, Dixon water, Dixon fat, Water fraction, and R2*. A feature selection process was performed to determine the optimal feature combination for the proposed automatic, 6‐tissue classification for sCT generation. Fuzzy c‐means was used for the automatic classification which was followed by voxel‐wise attenuation coefficient assignment as a weighted sum of those of the component tissues. Performance evaluation: MR data collected using the proposed pulse sequence were compared to those using a traditional two‐point Dixon approach. Image quality measures, including image resolution and uniformity, were evaluated using an MR ACR phantom. Data collected from 25 normal volunteers were used to evaluate the accuracy of the proposed method compared to the template‐based approach. Notably, the template approach is applicable here, that is normal volunteers, but may not be robust enough for patients with pathologies. RESULTS: The free breathing UTE‐mDixon pulse sequence yielded images with quality comparable to those using the traditional breath holding mDixon sequence. Furthermore, by capturing the signal before T2* decay, the UTE‐mDixon image provided lung and bone information which the mDixon image did not. The combination of Dixon water, Dixon fat, and the Water fraction was the most robust for tissue clustering and supported the classification of six tissues, that is, air, lung, fat, soft tissue, low‐density bone, and dense bone, used to generate the sCT. The thoracic sCT had a mean absolute difference from the template‐based (reference) CT of less than 50 HU and which was better agreement with the reference CT than the results produced using the traditional Dixon‐based data. CONCLUSION: MR thoracic acquisition and analyses have been established to automatically provide six distinguishable tissue types to generate sCT for MR‐based AC of PET/MR and for MR‐only RTP.

Purpose: MR-based pseudo-CT has an important role in MR-based radiation therapy planning and PET attenuation correction. The purpose of this study is to establish a clinically feasible approach, including image acquisition, correction, and CT formation, for pseudo-CT generation of the brain using a single-acquisition, undersampled ultrashort echo time (UTE)-mDixon pulse sequence. Methods: Nine patients were recruited for this study. For each patient, a 190-s, undersampled, single acquisition UTE-mDixon sequence of the brain was acquired (TE = 0.1, 1.5, and 2.8 ms). A novel method of retrospective trajectory correction of the free induction decay (FID) signal was performed based on point-spread functions of three external MR markers. Two-point Dixon images were reconstructed using the first and second echo data (TE = 1.5 and 2.8 ms). R2∗ images (1/T2∗) were then estimated and were used to provide bone information. Three image features, i.e., Dixon-fat, Dixon-water, and R2∗, were used for unsupervised clustering. Five tissue clusters, i.e., air, brain, fat, fluid, and bone, were estimated using the fuzzy c-means (FCM) algorithm. A two-step, automatic tissue-assignment approach was proposed and designed according to the prior information of the given feature space. Pseudo-CTs were generated by a voxelwise linear combination of the membership functions of the FCM. A low-dose CT was acquired for each patient and was used as the gold standard for comparison. Results: The contrast and sharpness of the FID images were improved after trajectory correction was applied. The mean of the estimated trajectory delay was 0.774 μs (max: 1.350 μs; min: 0.180 μs). The FCM-estimated centroids of different tissue types showed a distinguishable pattern for different tissues, and significant differences were found between the centroid locations of different tissue types. Pseudo-CT can provide additional skull detail and has low bias and absolute error of estimated CT numbers of voxels (−22 ± 29 HU and 130 ± 16 HU) when compared to low-dose CT. Conclusions: The MR features generated by the proposed acquisition, correction, and processing methods may provide representative clustering information and could thus be used for clinical pseudo-CT generation.

Purpose Inaccuracy of conventional four-dimensional (4D) flow MR imaging in the presence of random unsteady and turbulent blood flow distal to a narrowing has been an important challenge. Previous investigations have revealed that shorter echo times (TE) decrease the errors, leading to more accurate flow assessments. Methods In this study, as part of a 4D flow acquisition, an Ultra-Short TE (UTE) method was adopted. UTE works based on a center-out radial k-space trajectory that inherently has a short TE. By employing free induction decay sampling starting from read-out gradient ramp-up, and by combining the refocusing lobe of the slice select gradient with the bipolar flow encoding gradient, TEs of ≈1 msec may be achieved. Results Both steady and pulsatile flow regimes, and in each case a range of Reynolds numbers, were studied in an in-vitro model. Flow assessment at low and medium flow rates demonstrated a good agreement between 4D UTE and conventional 4D flow techniques. However, 4D UTE flow significantly outperformed conventional 4D flow, at high flow rates for both steady and pulsatile flow regimes. Feasibility of the method in one patient with Aortic Stenosis was also demonstrated. Conclusion For both steady and pulsatile high flow rates, the measured flow distal to the stenotic narrowing using conventional 4D flow revealed more than 20% error compared to the ground-truth flow. This error was reduced to less than 5% using the 4D UTE flow technique. Magn Reson Med 73:939–950, 2015. © 2014 Wiley Periodicals, Inc.

Purpose: Accurate and robust photon attenuation derived from MR is essential for PET/MR and MR-based radiation treatment planning applications. Although the fuzzy C-means (FCM) algorithm has been applied for pseudo-CT generation, the input feature combination and the number of clusters have not been optimized. This study aims to optimize both for clinically practical pseudo-CT generation. Methods: Nine volunteers were recruited. A 190-second, single-acquisition UTE-mDixon with 25% (angular) sampling and 3D radial readout was performed to acquire three primitive MR features at TEs of 0.1, 1.5, and 2.8 ms: the free-induction-decay (FID), the first and the second echo images. Three derived images, Dixon-fat and Dixon-water generated by two-point Dixon water/fat separation, and R2* (1/T2*) map, were also created. To identify informative inputs for generating a pseudo-CT image volume, all 63 combinations, choosing one to six of the feature images, were used as inputs to FCM for pseudo-CT generation. Further, the number of clusters was varied from four to seven to find the optimal approach. Mean prediction deviation (MPD), mean absolute prediction deviation (MAPD), and correlation coefficient (R) of different combinations were compared for feature selection. Results: Among the 63 feature combinations, the four that resulted in the best MAPD and R were further compared along with the set containing all six features. The results suggested that R2* and Dixon-water are the most informative features. Further, including FID also improved the performance of pseudo-CT generation. Consequently, the set containing FID, Dixon-water, and R2* resulted in the most accurate, robust pseudo-CT when the number of cluster equals to five (5C). The clusters were interpreted as air, fat, bone, brain, and fluid. The six-cluster Result additionally included bone marrow. Conclusion: The results suggested that FID, Dixon-water, R2* are the most important features. The findings can be used to facilitate pseudo-CT generation for unsupervised clustering. Please note that the project was completed with partial funding from the Ohio Department of Development grant TECH 11-063 and a sponsored research agreement with Philips Healthcare that is managed by Case Western Reserve University. As noted in the affiliations, some of the authors are Philips employees

Purpose: Emerging technologies such as dedicated PET/MRI and MR-therapy systems require robust and clinically practical methods for determining photon attenuation. Herein, we propose using novel MR acquisition methods and processing for the generation of pseudo-CTs. Methods: A single acquisition, 190-second UTE-mDixon sequence with 25% (angular) sampling density and 3D radial readout was performed on nine volunteers. Three water-filled tubes were placed in the FOV for trajectory-delay correction. The MR data were reconstructed to generate three primitive images acquired at TEs of 0.1, 1.5 and 2.8 ms. In addition, three derived MR images were generated, i.e. two-point Dixon water/fat separation and R2* (1/T2*) map. Furthermore, two spatial features, i.e. local binary pattern (S-1) and relative spatial coordinates (S-2), were incorporated. A direct-mapping operator was generated using Artificial Neural Networks (ANNs) for transforming the MR features to a pseudo-CT. CT images served as the training data and, using a leave-one-out method, for performance evaluation using mean prediction deviation (MPD), mean absolute prediction deviation (MAPD), and correlation coefficient (R). Results: The errors between measured CT and pseudo-CT declined dramatically when the spatial features, i.e. S-1 and S-2, were included. The MPD, MAPD, and R were, respectively, 5±57 HU, 141±41 HU, and 0.815±0.066 for results generated by the ANN trained without the spatial features and were 32±26 HU, 115±18 HU, and 0.869±0.035 with the spatial features. The estimation errors of the pseudo-CT were smaller when both the S-1 and S-2 were used together than when either the S-1 or the S-2 was used. Pseudo-CT generation (256×256×256 voxels) by ANN took < 0.5 s using a computer having an Intel i7 3.4GHz CPU and 16 GB RAM. Conclusion: The proposed direct-mapping ANN approach is a technically accurate, clinically practical method for pseudo-CT generation and can potentially help improve the accuracy of MR-AC and MR-RTP applications. Please note that the project was completed with partial funding from the Ohio Department of Development grant TECH 11-063 and a sponsored research agreement with Philips Healthcare that is managed by Case Western Reserve University. As noted in the affiliations, some of the authors are Philips employees

Purpose: The ultrashort echo-time (UTE) sequence is a promising MR pulse sequence for imaging cortical bone which is otherwise difficult to image using conventional MR sequences and also poses strong attenuation for photons in radiation therapy and PET imaging. The authors report here a systematic characterization of cortical bone signal decay and a scanning time optimization strategy for the UTE sequence through k-space undersampling, which can result in up to a 75% reduction in acquisition time. Using the undersampled UTE imaging sequence, the authors also attempted to quantitatively investigate the MR properties of cortical bone in healthy volunteers, thus demonstrating the feasibility of using such a technique for generating bone-enhanced images which can be used for radiation therapy planning and attenuation correction with PET/MR. Methods: An angularly undersampled, radially encoded UTE sequence was used for scanning the brains of healthy volunteers. Quantitative MR characterization of tissue properties, including water fraction and R2∗ = 1/T2∗, was performed by analyzing the UTE images acquired at multiple echo times. The impact of different sampling rates was evaluated through systematic comparison of the MR image quality, bone-enhanced image quality, image noise, water fraction, and R2∗ of cortical bone. Results: A reduced angular sampling rate of the UTE trajectory achieves acquisition durations in proportion to the sampling rate and in as short as 25% of the time required for full sampling using a standard Cartesian acquisition, while preserving unique MR contrast within the skull at the cost of a minimal increase in noise level. The R2∗ of human skull was measured as 0.2–0.3 ms−1 depending on the specific region, which is more than ten times greater than the R2∗ of soft tissue. The water fraction in human skull was measured to be 60%–80%, which is significantly less than the >90% water fraction in brain. High-quality, bone-enhanced images can be generated using a reduced sampled UTE sequence with no visible compromise in image quality and they preserved bone-to-air contrast with as low as a 25% sampling rate. Conclusions: This UTE strategy with angular undersampling preserves the image quality and contrast of cortical bone, while reducing the total scanning time by as much as 75%. The quantitative results of R2∗ and the water fraction of skull based on Dixon analysis of UTE images acquired at multiple echo times provide guidance for the clinical adoption and further parameter optimization of the UTE sequence when used for radiation therapy and MR-based PET attenuation correction.

Purpose: To develop semi‐automatic methods for robust segmentation of head and neck anatomy through structure‐specific image features selection from multiple MRI sequences. Methods: We developed a semi‐automatic approach for probabilistic segmentation of head and neck anatomy that combines several candidate segmentations using random forests and structure‐specific image features selection from multiple MRI sequences. Starting from user‐drawn line segments on a single axial image, in each tree, structure segmentations are generated from composite feature images computed at the tree split nodes. Composite features are computed through random combinations of the MRI sequences and image features, including texture and Laplacian of Gaussian gradients. We used seven different MRI pulse sequences obtained from a Phillips scanner, including Dixon water only T1, Dixon fat only T1, Dixon in‐phase, Dixon opposed‐phase, T1‐weighted, T2‐weighted and diffusion weighted whole body imaging with background body signal suppression (DWIBS). The highest quality anatomical structure segmentation in a tree is extracted from all available split node segmentations with segmentation quality measured using within‐segment statistics including kurtosis, edge quality, and inter‐segment dissimilarity. The joint segmentation with probabilistic score is obtained by combining all the generated tree segmentations. Results: We used our approach to segment anatomical structures in head and neck images using the seven MRI sequences listed above. The algorithm returned the segmentations, their probabilistic scores, and the highest quality feature images for each structure. We evaluated our approach by computing DICE overlap scores with manually segmented ground truths and obtained an overall segmentation accuracy of 80% for the various radiotherapeutically‐relevant anatomic structures. Conclusion: We developed a semi‐automatic, random forest‐based probabilistic segmentation of head and neck anatomy by combining multiple MRI pulse sequences. Our approach achieves accurate segmentation for various radiotherapeutically‐relevant anatomic structures. Additionally, it extracts the highest quality composite feature images for obtaining optimal segmentation for the same structures. Phillips

Purpose: To generate electron density (ED) maps and digitally reconstructed radiographs (DRRs) of the pelvis based only on magnetic resonance imaging (MRI). Methods: A conventional 3D fast‐field dual‐echo sequence was used to acquire image data of 8 healthy subjects with a Philips 3.0T Ingenia TX system in approximately 1:50 min scan time per volunteer. Dixon reconstruction of the nearly out‐of‐phase (echo 1) and in‐phase images (echo 2) allowed for water and fat classification. A bone‐enhanced image was generated by automatically thresholding the noise level of the in‐phase image with subsequent background removal. ED maps were then produced by assigning known bulk electron densities to the classified bone and tissue fractions. A bone probability atlas derived from CT data was registered to the ED map in order to filter out misclassified voxels. Finally, DRRs were reconstructed from bone‐enhanced images as well as from ED maps. Results: The proposed MRI sequence with subsequent Dixon reconstruction and probabilistic filtering makes it possible to classify cortical bone, soft tissue and adipose tissue in the pelvis and yields ED maps and corresponding DRRs. Bowel content was misclassified as cortical bone or air and compromised the segmentation in some slices as well as in the DRRs. Automatic probabilistic atlas filtering can significantly reduce artifacts induced by bowel content without affecting pelvic bone structures markedly. In total, the artifact/bone fraction dropped from 1.7 before filtering to 0.2 after filtering. The average reduction of artifact volume is 87%, and the average bone preservation is 99%. Remaining artifacts are spatially close to the true bone in areas of positive bone probability. Conclusion: This study demonstrated the feasibility of generating realistic ED maps of the pelvis by using MRI only. The method has the potential to become an essential component of emerging applications such as MR‐only‐based radiation therapy planning. All authors have the following relevant financial interest or relationship to disclose with regard to the subject matter of this presentation: Company name: Philips Research; Type of relationship: Employee