Abstract
MR-based attenuation correction is instrumental for integrated PET/MR imaging. It is generally achieved by segmenting MR images into a set of tissue classes with known attenuation properties (e.g., air, lung, bone, fat, soft tissue). Bone identification with MR imaging is, however, quite challenging, because of the low proton density and fast decay time of bone tissue. The clinical evaluation of a novel, recently published method for zero-echo-time (ZTE)–based MR bone depiction and segmentation in the head is presented here. Methods: A new paradigm for MR imaging bone segmentation, based on proton density–weighted ZTE imaging, was disclosed earlier in 2014. In this study, we reviewed the bone maps obtained with this method on 15 clinical datasets acquired with a PET/CT/MR trimodality setup. The CT scans acquired for PET attenuation-correction purposes were used as reference for the evaluation. Quantitative measurements based on the Jaccard distance between ZTE and CT bone masks and qualitative scoring of anatomic accuracy by an experienced radiologist and nuclear medicine physician were performed. Results: The average Jaccard distance between ZTE and CT bone masks evaluated over the entire head was 52% ± 6% (range, 38%–63%). When only the cranium was considered, the distance was 39% ± 4% (range, 32%–49%). These results surpass previously reported attempts with dual-echo ultrashort echo time, for which the Jaccard distance was in the 47%–79% range (parietal and nasal regions, respectively). Anatomically, the calvaria is consistently well segmented, with frequent but isolated voxel misclassifications. Air cavity walls and bone/fluid interfaces with high anatomic detail, such as the inner ear, remain a challenge. Conclusion: This is the first, to our knowledge, clinical evaluation of skull bone identification based on a ZTE sequence. The results suggest that proton density–weighted ZTE imaging is an efficient means of obtaining high-resolution maps of bone tissue with sufficient anatomic accuracy for, for example, PET attenuation correction.
Hybrid PET/MR imaging is an emerging modality with potential applications in oncology, cardiology, and neurology. Since its commercial introduction in 2011, this new technology has been steadily expanding its installed base. Currently, 4 models—2 fully integrated and 2 sequential—are commercially available for clinical applications.
One of the main challenges of hybrid PET/MR imaging is attenuation correction. Indeed, a critical step for quantitative PET is determining and compensating for the signal attenuation introduced by the patient. In PET/CT scanners, this attenuation correction can be achieved adapting the x-ray attenuation information provided by the CT in PET/CT scanners. In the case of PET/MR scanners, the attenuation is estimated by segmenting MR data to identify different tissue classes (e.g., fat, soft tissue, lung) (1,2). It follows that the accuracy of reconstructed PET data will be determined by how many tissue classes can be correctly identified.
Of particular relevance is the correct identification of bone tissue, due to its distinctively high attenuation value. However, in contrast to CT, MR is intrinsically much less suited for the depiction of cortical bone structures: the low proton density (PD) (∼20% of water) and heterogeneous structure lead to weak and short-lived signals (T2, ∼390 μs at 3T) (3). Atlas methods have been proposed to overcome this obstacle (4–6), but they are not well suited for patients with off-norm anatomy, a common occurrence in PET indications such as oncology, for which postsurgical examinations are frequent.
Conventional sequences with echo times in the millisecond range are therefore too slow for direct bone signal detection. In contrast, ultrashort-echo-time (UTE) sequences with center-out k-space acquisition enable sufficiently fast data acquisition to capture the bone signals (7,8). Several studies have been published discussing the technical feasibility of UTE imaging for PET/MR attenuation correction. Generally, long T2 suppression methods (e.g., echo subtraction, saturation prepulses, or multiple sequences) are necessary to separate the bone signal from soft tissue (9). Arguably, the most widely accepted approach of bone segmentation for MR-based attenuation correction is based on the postprocessing of dual-echo UTE images (10–17).
A recent development has been the publication of a new bone identification technique, based on 3-dimensional radial zero-echo-time (ZTE) imaging (18). This sequence (similar to pointwise encoding time reduction with radial acquisition (19) and, to a certain extent, SWeep Imaging with Fourier Transformation (20)) has been proven to provide high-resolution isotropic images, suitable for bone segmentation—1.4 mm3 in this study. Furthermore, bone segmentation is achieved without the need of preparation pulses or multiple echoes, making it a time-efficient acquisition. Rather than obtaining contrast from T2 relaxation time differences, this method explores the use of PD differences.
The goal of the present study was to evaluate ZTE acquisitions for skull bone identification in PET/MR. For this purpose, a PET/CT/MR trimodality setup has been used to acquire clinical patient datasets, providing both the MR sequences under study and the gold-standard CT reference.
MATERIALS AND METHODS
Data Acquisition
Data were acquired using a trimodality setup consisting of a Discovery 750w 3T MR system (GE Healthcare) located in a room adjacent to a Discovery 690 time-of-flight PET/CT system (GE Healthcare). Patients were transported between the 2 systems using a dedicated transfer device, enabling a consistent patient placement between the PET/CT and MR imaging systems (GE Healthcare) (21).
Fifteen patients referred for a clinical oncology PET/CT examination were enrolled in this study. The average patient age was 59 ± 17 y (range, 32–88 y), the average weight was 75 ± 9 kg (range, 66–100 kg), and the average body mass index was 25 ± 3 kg/m2 (range, 20–33 kg/m2). Sixty percent of the patients were men and 40% women. Eleven of the patients displayed artificial alterations of the skull (craniotomies and craniectomies at various locations). The present study did not involve any extra radiation dose delivered to the patients, since the used CT was part of the clinical routine PET/CT examination. This study was approved by the institutional ethics committee, and written informed consent was obtained from all patients before the examination.
The voluntary MR examination was performed during the resting time after the radiotracer injection, so the total time at the site for the patient (PET/CT/MR, compared with standard PET/CT) was not altered, allowing approximately 30 min of MR scan time. The MR protocol included a ZTE acquisition, as well as further sequences used for the anatomic referencing of whole-body PET findings and not relevant for the current study. Thirteen patients were scanned with a 32-channel head coil (MR Instruments Inc.) and the remaining 2 with a head and neck coil (geometry embracing method head and neck unit; GE Healthcare).
The present study was based on the new bone imaging paradigm recently proposed by Wiesinger et al., based on proton density (PD)–weighted ZTE (22,23). The acquisition lasted 172 s, covering a 26-cm transaxial and axial field of view with a resolution of 1.4 × 1.4 × 1.4 mm. Four excitations were acquired, with a flip angle of 1° and 62.5-kHz bandwidth. The k-space was parsed using a 3-dimensional radial trajectory. The missing region at the center of k-space was filled with a second radial acquisition, with a reduced number of spokes and bandwidth, similar to the so-called water- and fat-suppressed proton projection MRI method (24). Before reconstruction, the 2 acquisitions were merged with a linear transition in the overlap region, to prevent the formation of artifacts. The raw data were reconstructed offline using standard 3-dimensional gridding (25,26).
The PET/CT acquisition followed the standard protocol for a clinical oncology study. The helical CT scan acquired for PET attenuation correction (120 kV; 15–80 mA with automatic dose modulation; rotation time, 0.5 s; helical thickness, 3.75 mm; pitch, 39.37 mm/rotation; matrix size, 512 × 512, 1.4 × 1.4 × 3.3 mm3) was used in this study as ground truth for the evaluation of MR-based bone imaging. The subsequent PET acquisition is not relevant for the present study.
Postprocessing
A detailed description of the segmentation procedure to generate bone and tissue masks from the reconstructed ZTE volumes can be found in Wiesinger et al. (22). The procedure essentially consists of a first intensity equalization step, followed by a set of histogram-based thresholding operations: first, a logarithmic rescaling of the normalized image intensities is applied to enhance the bone tissue. Then, a mask of the main body object is obtained to eliminate the background and coil elements. Within this body region, the bone and air regions are defined by thresholding operations based on a gaussian fitting of the main histogram peaks. Simple thresholding was also used to identify internal air cavities. The precise threshold values were defined based on the mean and full width at half maximum values of the noise and signal peaks. The formulas to generate the thresholds based on these were defined empirically (±2.25 × full width at half maximum) but remained constant for all patients.
Data Analysis
The segmentation results were quantitatively evaluated using the Jaccard distance between CT and ZTE bone masks. This overlap distance is equal to 1 (or 100%) for maximal dissimilarity and to 0 (or 0%) for complete equality:The CT datasets were automatically registered to the segmented ZTE volumes using the rigid registration tool of the Integrated Registration plugin in an Advantage Workstation (GE Healthcare).
Notice that the results of the Jaccard distance are dependent on the threshold used to generate the CT bone mask. We determined this threshold by plotting, for each patient, the Jaccard distance for a wide range of possible threshold values. The optimal threshold for each patient was recorded, and the average threshold over all patients was used in the subsequent analysis.
This study of the quality of ZTE bone masks, compared with CT, was performed both considering the entire head image and considering only the cranium (to be precise, the axial span of the neurocranium). Studying both areas separately is useful for 2 reasons: first, it removes from the measure effects such as dental artifacts and jaw and vertebral (nonrigid) misalignment. Second, it provides a better idea of the different performance that can be expected when using ZTE for the attenuation correction of brain PET studies, as opposed to head and neck imaging.
In terms of clinical evaluation, each CT and ZTE image pair was inspected by a dual board-certified radiologist and nuclear medicine physician using the Advantage Workstation Volume Viewer software. Axial, sagittal, and coronal 2-dimensional views, as well as volume renderings, were used to determine the anatomic accuracy of ZTE bone maps in comparison to the CT gold standard. Selected areas of interest were assessed: the calvaria was subdivided into frontal, parietal, temporal, and occipital areas. The inner ear (bone/fluid interfaces), the mastoid (bone/air/fluid interfaces), and the orbit (bone/muscle/fat interfaces) were assessed. Calcifications within the brain found consistently in adults were addressed as well (choroid plexus of lateral ventricles, pineal gland). Midline brain/fluid interfaces (third ventricle, fourth ventricle) were also part of the assessment. Another defined area of interest was artificial bone alterations, such as drilling holes or osteoplastic craniectomy sites.
RESULTS
Figures 1A and 1C illustrate a typical single-echo ZTE dataset (in this case acquired with the geometry embracing method head and neck coil). Notice the flat contrast, uncommon in diagnostic imaging, intended to facilitate the segmentation. The use of very low flip angles, an intrinsic property of ZTE imaging, has in this case the benefit of minimizing saturation effects in tissues with long relaxation times. In consequence, the low-intensity range of the image is occupied by little more than air and bone tissue, which is highlighted by the logarithmic intensity rescaling of the images (Figs. 1B and 1D).
Axial (A and B) and sagittal views (C and D) of reconstructed ZTE dataset. Logarithmic intensity rescaling (B and D) was applied to enhance bone tissue and facilitate its separation from background and internal air cavities.
An example of the rescaled intensity histogram of a patient ZTE dataset can be found in Figure 2A. After intensity equalization, all datasets present clear peaks for the air and soft-tissue classes, easily identified by the automated gaussian fitting. Figures 2B and 2C show a typical segmentation result, obtained using the average and SD of these peaks to define the expected intensity range of bone tissue.
(A) Typical intensity histogram of ZTE dataset after logarithmic rescaling. Gaussian fitting result for soft-tissue and background peaks is indicated by dashed lines. Expected intensity range of bone tissue used in segmentation is indicated by dotted lines. (B) Axial view of corresponding ZTE volume. (C) Soft tissue (blue) and bone mask (green) obtained with segmentation.
Figure 3 shows a side-by-side comparison of a segmented ZTE volume with the corresponding (registered) CT. Supplemental Figures 1 and 2 (supplemental materials are available at http://jnm.snmjournals.org) provide further examples, including artificial alterations of the skull.
Volume rendering (top) and sagittal views (bottom) of CT dataset (left) and corresponding segmented ZTE dataset (right). Dental artifacts on both images, minor misclassification of cartilage and auditory canal air, and oversegmentation around sinusal cavities exist.
Before the quantitative evaluation of the segmentation results, a suitable threshold had to be determined for the generation of bone maps from the reference CT datasets. Figure 4 shows the plot of Jaccard distances as a function of the selected CT threshold, for all patients. The optimal CT threshold was found to be 471 ± 88 Hounsfield units (HUs) (range, 356–636 HUs). As can be appreciated in the figure, it was a relatively wide and flat minimum in all cases. The difference in the Jaccard distance, when using this average threshold (471 HUs) rather than the individual optimum for each patient, was found to be no larger than +1%.
Jaccard distance between CT and ZTE bone masks, as function of CT threshold value. Each line in plot corresponds to 1 of the patients.
Table 1 shows the Jaccard distance measured between the CT and ZTE bone masks, for all patients. The measurements were performed twice, once considering the entire imaged area and again considering only the cranium. The average distance for the entire area was 52% ± 6% (range, 38%–63%) and for the cranium only 39% ± 4% (range, 32%–49%). A visual representation of the agreement between CT and ZTE bone masks can be found in Figure 5.
Jaccard Overlap Distances Measured Between ZTE and CT Bone Masks
Logarithmic rescaling of ZTE dataset (A and B) and corresponding map of agreement between ZTE and CT bone masks (C and D). Green pixels indicate true-positive bone identification, red pixels indicate false-negative (missed bone), and yellow pixels indicate false-positive. Residual misregistration can be perceived in, for example, occipital region. Notice surgical bone alteration.
As with UTE segmentation approaches, partial-volume effects on tissue/air interfaces led to false-positive bone voxels. These are accounted for by the segmentation procedure and eliminated using a Euclidean distance transform (similar to an erosion operation). However, some false-positive voxels remain in isolated areas such as the auditory canal and nasal fossae. The fraction of false-positive voxels was measured to be 80% ± 6%. Most of these corresponded to high-density CT regions, close to but below the segmentation threshold. The percentage of bone voxels being classified as air and vice versa was in all cases below 1%.
The results of the comparative clinical scoring of the ZTE bone maps are summarized in Table 2. Notice that these observations are specific to segmented bone maps, and some of the structures found missing there are still visible in the unsegmented ZTE data.
Clinical Scoring of ZTE Bone Masks in Comparison with CT Gold Standard
The calvaria yielded mainly minor issues, irrespective of the site. However, a few major issues were found in all calvaria sites, mainly frontally, consisting of apparent osteolysis on ZTE without correlative on CT. The apparent osteolysis was found to be caused by bone being misclassified as soft tissue. In most cases, this misclassification was limited to clusters of less than 5 voxels, often in regions with thin cortical bone layers separated by an air cavity or thick layer of trabecular bone. The misclassification seldom affected both layers of cortical bone (creating a perforation), and the overall morphology and continuity of the cranium were preserved in all cases. Bone union ridges and vessel foramina were visible in several cases.
Bone/fluid interfaces with high anatomic detail, as found in the inner ear, were consistently found to be a major issue for ZTE. Furthermore, localized misclassifications of bone as air were noticed in some regions in which bone tissue was too dense to yield sufficient MR signal. The mastoid air cells, although being complex anatomic structures as well, were found to pose fewer issues.
Almost no major issues were found in the orbit. Only in 1 subject, pseudoossifications were found in the conal space on ZTE. Localized false-positive bone voxels were occasionally found on the cornea.
The physiologic intracranial calcifications in the pineal gland (respectively, choroid plexus of the lateral ventricles) were depicted on ZTE in almost all subjects (respectively, half of the subjects). Voxels incorrectly labeled as hard material were displayed on ZTE images in the third ventricle more often than in the fourth ventricle. No cases were found in which brain structures were misclassified as air.
Contrary to previously reported skull segmentation studies based on UTE sequences (27), the presence of tendons on the bone masks obtained with ZTE was minor. Nasal cartilage was partially included in the bone masks.
In artificial osseous alterations (11/15 subjects), minor issues were found in 6 subjects and major issues in 3. Those major issues consisted of apparent osteolysis (bone to tissue misclassification, 2 subjects) and air (bone to air misclassification, 1 subject), for which reimplanted bone had a normal appearance on CT. Metal artifacts due to dental implants led to unwanted structures in the bone mask, corresponding to the edges of the MR signal void.
DISCUSSION
To our knowledge, this is the first clinical evaluation of skull bone identification based on a ZTE sequence, comparing the obtained bone masks with corresponding CT data. We have shown that single-echo PD-weighted acquisitions are well suited for automated bone identification. Furthermore, the quantitative and qualitative evaluation of the obtained bone masks compare positively with alternative methods found in the literature.
Our results are in agreement with those reported in the methodologic paper by Wiesinger et al. (22). They confirm that ZTE MR imaging is well suited for the depiction of tissues with fast decaying resonance signals and, in particular, for the efficient acquisition of PD-weighted datasets. These can be robustly segmented into air, soft-tissue, and bone compartments.
Relying on single-echo PD imaging is more time-efficient than the multiecho/multisequence T2* estimation approaches commonly found in the prior work (10–17). The saved time can be reinvested in the acquisition of higher resolution datasets, further improving the resulting bone maps. The ability to obtain images of higher resolution is restricted by the intrinsic signal-to-noise limitations of low flip-angle ZTE acquisition, leading to a tradeoff between the achievable resolution and the need for additional averaging.
The inverse logarithmic scaling of ZTE datasets provides a suitable means of direct inspection of the bone tissue information captured by the sequence. The image background and partial-volume effects in the air–tissue interfaces must be removed before rendering tools can be used.
The bone masks yielded by ZTE postprocessing were evaluated relative to CT, the current gold standard for bone imaging. Quantitative evaluation using the Jaccard distance shows better performance than previously reported for dual-echo UTE (27), which ranged between 47% and 79%, as opposed to the 38%–63% obtained in this study.
The wide minima found during the optimization of the threshold for CT bone segmentation suggest that the bone structures detectable with ZTE are relatively well defined in terms of CT intensity (i.e., there are no undetected structures with similar attenuation properties, nor an excess of gradual intensity transitions). Hence, the ZTE bone mask is well suited for the purpose of attenuation modeling. Still, future work will aim at improving the processing of ZTE datasets to yield nonbinary estimations of tissue density.
Overall, the anatomic detail of ZTE bone masks was insufficient for direct clinical inspection but acceptable for attenuation-correction purposes. Despite air cavity differentiation working remarkably well because of the PD contrast of the sequence, the problem of false-positive bone identification in certain interfaces affected by partial-volume effects remains.
The general morphology of the cranium was correctly captured in the ZTE masks. Apparent osteolyses were found, especially in the frontal calvaria. These were due to local (typically <1 cm) misclassification of thin trabecular bone layers as soft tissue. Often, only 1 layer was affected, suggesting that the misclassification would have a limited impact on attenuation estimation. Block effects related to the implementation of the intensity inhomogeneity correction were also identified as an occasional source of local bone misclassification. This problem will need to be addressed in future studies.
Structural detail in the inner ear was almost completely lost in ZTE bone masks of all patients, compared with CT (severely limiting their use for clinical diagnosis). Although structural detail was also diminished in the mastoid, with loss of some thin air cell septa, this did not affect the overall interpretation of the scans, because no relevant information was lost. More problematic, from the point of view of attenuation estimation, would be the misclassification of dense bone tissue as air. This is probably an inherent limitation of the acquisition, due to the low density of signal-emitting tissue. Regardless of the highly localized nature of this issue, further study will be required to assess its impact in attenuation correction.
Fat/muscle/bone interfaces such as in the retrobulbar region yielded satisfactory results, except in 1 subject in whom apparent conal ossification on ZTE was possibly inferred by eyeball movements.
We could show that even small calcifications, such as in the pineal gland, are reliably interpreted as hard structures by ZTE. The incorrect hard structures displayed by ZTE in the third ventricle were likely induced by a combination of magnetization saturation and the thin, slitlike configuration of this cerebrospinal fluid space, with 2 brain/fluid interfaces being close, as opposed to the fourth ventricle, for which pseudocalcifications were found less commonly. Classification issues in the ventricles have also been reported for UTE-based methods (28). Anatomy-based postprocessing may need to be included in the segmentation to mitigate this issue.
Similarly, the throat and vertebrae area would benefit from dedicated processing, which could account both for the lower density bone tissue and for the false-positives in the oropharyngeal cavity interfaces.
Apparent lytic or even air-containing lesions in reimplanted bone on ZTE, being observed in 3 of 11 subjects with postoperative changes, do have grave clinical influence. Extreme caution should therefore be taken when interpreting ZTE datasets around such areas. Signal voids around implanted metallic objects are a serious limitation of the proposed method, and correction approaches are currently being investigated.
Further work will be directed at extending the use of ZTE bone segmentation to whole-body applications. Other potential applications of the ZTE sequence for PET/MR imaging would be lung visualization, to improve the estimation of parenchyma attenuation (29) and the detection of local coils.
CONCLUSION
This is the first, to our knowledge, clinical evaluation of skull bone identification based on a ZTE sequence. The results suggest that single-echo, PD-weighted ZTE imaging is an efficient means of obtaining anatomically accurate maps of bone tissue for the purpose of PET attenuation correction. The segmentation results show better bone depiction and separation from air cavities and collagenous tissue than previously reported methods. Further work will be aimed at obtaining nonbinary estimations of bone density.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. No potential conflict of interest relevant to this article was reported.
Footnotes
- © 2015 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication October 14, 2014.
- Accepted for publication December 15, 2014.