ZT-AC:  Zero TE based PET/MR attenuation correction
Florian Wiesinger1, Sandeep Kaushik2, Dattesh Shanbhag2, Venkata Chebrolu2, Vivek Vaidya2, Sangtae Ahn3, Lishui Cheng3, Andrew Leynes4, Jaewon Yang4, Thomas Hope4, and Peder Larson4

1GE Global Research, Munich, Germany, 2GE Global Research, Bangalore, India, 3GE Global Research, Niskayuna, NY, United States, 4University of California San Francisco, San Francisco, CA, United States

Synopsis

We describe a novel PET/MR attenuation correction (AC) method based on zero TE MR imaging that is fast, robust and accounts for bone. The method, termed ZT-AC, was tested for N=10 PET/MR patients in the head and compared relative to gold-standard CT-AC and atlas-based AC methods.

Purpose

PET/MR is a novel hybrid imaging modality which promises to influence diagnostic imaging in the near and long-term future. Despite dedicated and focused research efforts, PET/MR attenuation correction (AC) as required for quantitative PET is still considered challenging; especially when compared to PET/CT. A particular challenge is the correct characterization of bone; i.e. the tissue with the highest PET attenuation but which appears invisible in standard MRI. Recently, zero TE was demonstrated for fast and reliable MR bone depiction [1]. Here we extend that work and demonstrate zero TE for PET/MR attenuation correction.

Methods and Materials

N=10 oncology patients were scanned using a 3-tesla, time-of-flight GE Signa PET/MR scanner (GE Healthcare, Waukesha, WI). Zero TE images were acquired using proton density (PD) weighted parameter settings (FA=1deg, BW=±62.5kHz, FOV=26.4cm, res=2.4mm, 48400 spokes acquired in 41sec) and a GEM head array coil for signal reception. All zero TE images were bias correction using N4ITK [2] and normalized relative to soft tissue. Additionally, outside patient structures (e.g. plastic components from the coil housing and the table) were removed using connected component analysis. For each patient, available CT datasets were aligned to the PET/MR by rigidly registering the whole-body CT to zero TE data followed by diffeomorphic non-rigid registration [3,4]. The zero TE images were converted into pseudo CTs (ZT) using either: 1) image segmentation based on local, pixel-wise thresholding [ZTSEG(AIR)=-1000HU for zero TE < 0.2 and ZTSEG(SOFT-TISSUE)=40HU for zero TE > 0.8] and linear scaling [ZTSEG(BONE)=-2350*(ZTE-1)+40) for 0.2 < zero TE < 0.8] (cf. Fig. 1), or 2) machine learning based on a random forests regression model with each voxel in the pseudo CT represented by a feature vector derived from the corresponding voxel and its surrounding local neighborhood in the zero TE image. The corresponding pseudo CTs are referred to as ZTSEG and ZTML, respectively. PET images were reconstructed using a time-of-flight, ordered subset expectation maximization (TOF-OSEM) algorithm (2 iterations, 24 subsets). For PET/MR attenuation correction, existing atlas-based AC (AT-AC) [4] and gold standard CT-AC were compared relative to ZTSEG-AC and ZTML-AC.

Results

The left subplot of Fig. 1 compares zero TE in linear (left), and negative linear scale (middle) vs. CT (right). The inverted zero TE appears similar to the CT by depicting soft-tissue dark and bone grayish against a bright background corresponding to air. The scatter plot, depicted on the right, justifies the linear scaling for ZTSEG(BONE) as described above. Figure 2 illustrates a close correspondence between the true CT (left) and the two zero TE derived pseudo CTs; i.e. ZTSEG (middle) and ZTML (right). By taking into account neighborhood information and local non-linearity, the Machine Learning based approach better captures tissue / air interfaces otherwise affected by partial volume effects. Figure 3 illustrates the same comparison after rescaling the CTs towards 511keV PET attenuation maps [1/mm] and pasting hardware templates to account for attenuation of auxiliary MR equipment (i.e. coil and table). The energy rescaling and the down-sampling to PET resolution minimize / smooth differences. Figure 4 shows the same comparison for the reconstructed PET images. The error maps are calculated relative to gold-standard PETCT-AC (left) i.e. (PETX-PETCT-AC)/max(PETCT-AC) and demonstrate significantly improved PET quantitation using ZT-derived pseudo CTs. Aforementioned advantages of the Machine Learning based approach also translate into more accurate PET images using ZTML-AC; especially around the trachea, the sinuses and implants.

Discussion and Conclusions

The flat PDw soft-tissue contrast and the efficient capture of short-lived bone signals renders zero TE well suited for pseudo CT conversion as required for PET/MR-AC. This becomes especially apparent when inverting the color scale (cf. Fig. 1). In comparison to a pixel-based pseudo CT conversion (using thresholding and linear scaling), the Machine Learning based framework is more flexible in the sense that it allows a local non-parametric mapping by including local neighborhood and high-dimensional information. This becomes apparent around air cavities often affected by partial volume effects. Figure 5 exemplarily demonstrates improved robustness of ZTML for implants (less pronounced and more localized metal artifact) and its applicability for other anatomies such as the pelvis. In contrast to previous Machine Learning attempts requiring multi-contrast (PD, T1, T2, …) input data, the presented method achieves the same with only a single channel input information (i.e. PDw Zero TE) acquired in a scan of only 41sec. From a PET quantitation perspective, the performed N=10 patient study indicates significant improved PET quantitation (relative to gold standard CT-AC) for ZT-AC in the head when compared to existing alternatives.

Acknowledgements

No acknowledgement found.

References

[1] F. Wiesinger et al. "Zero TE MR bone imaging in the head." Magnetic Resonance in Medicine (2015). [2] H. Patel et al. "Automatic Determination of Anatomical Correspondences for Multimodal Field of View Correction." Pattern Recognition. Springer International Publishing, 2014. 432-442. [3] BB. Avants et al. "The Insight ToolKit image registration framework." Frontiers in neuroinformatics 8 (2014). [4] SD Wollenweber et al. "Evaluation of an atlas-based PET head attenuation correction using PET/CT & MR patient data." Nuclear Science, IEEE Transactions on 60.5 (2013): 3383-3390.

Figures

Figure 1 (patient #01): Zero TE images in linear (left) and inverted (middle) grayscale compared to the registered CT (right). The scatterplot demonstrates a linear relationship for converting zero TE bone signals into pseudo CT according to: ZTSEG(BONE)=-2350*(ZT-1)+40. The exact shape and slope depends on parameter settings used.

Figure 2 (patient #10): True CT (left) compared to zero TE derived pseudo CTs (ZTSEG(middle) and ZTML(right)).

Figure 3 (patient #04): PET photon attenuation maps: CT-AC (left), AT-AC, ZTSEG-AC, ZTML-AC (right).

Figure 4 (patient #04): CT-AC PET reconstruction (left) and relative error maps (±10% windowing) corresponding to AT-AC, ZTSEG-AC and ZTML-AC.

Figure 5: Zero TE (left), ZTML (middle), CT (right) for dental implants (top) and the pevis (bottom).



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3559