Practical joint reconstruction of activity and attenuation with autonomous scaling for time-of-flight PET

Supplemental Transmission Aided Attenuation Correction for Quantitative Cardiac PET

Quantitative PET attenuation correction (AC) for cardiac PET/CT and PET/MR is a challenging problem. We propose and evaluate an AC approach that uses coincidences from a relatively weak and physically fixed sparse external source, in combination with that from the patient, to reconstruct μ -maps based on physics principles alone. The low 30 cm3 volume of the source makes it easy to fill and place, and the method does not use prior image data or attenuation map assumptions. Our supplemental transmission aided maximum likelihood reconstruction of attenuation and activity (sTX-MLAA) algorithm contains an attenuation map update that maximizes the likelihood of terms representing coincidences originating from tracer in the patient and a weighted expression of counts segmented from the external source alone. Both external source and patient scatter and randoms are fully corrected. We evaluated performance of sTX-MLAA compared to reference standard CT-based AC with FDG PET/CT phantom studies; including modeling a patient with myocardial inflammation. Through an ROI analysis we measured ≤ 5 % bias in activity concentrations for PET images generated with sTX-MLAA and a TX source strength ≥ 12.7 MBq, relative to CT-AC. PET background variability (from noise and sparse sampling) was substantially reduced with sTX-MLAA compared to using counts segmented from the transmission source alone for AC. Results suggest that sTX-MLAA will enable quantitative PET during cardiac PET/CT and PET/MR of human patients.

The SNR of time-of-flight positron emission tomography data for joint reconstruction of the activity and attenuation images

Objective.Measurement of the time-of-flight (TOF) difference of each coincident pair of photons increases the effective sensitivity of positron emission tomography (PET). Many authors have analyzed the benefit of TOF for quantification and hot spot detection in the reconstructed activity images. However, TOF not only improves the effective sensitivity, it also enables the joint reconstruction of the tracer concentration and attenuation images. This can be used to correct for errors in CT- or MR-derived attenuation maps, or to apply attenuation correction without the help of a second modality. This paper presents an analysis of the effect of TOF on the variance of the jointly reconstructed attenuation and (attenuation corrected) tracer concentration images.Approach.The analysis is performed for PET systems that have a distribution of possibly non-Gaussian TOF-kernels, and includes the conventional Gaussian TOF-kernel as a special case. Non-Gaussian TOF-kernels are often observed in novel detector designs, which make use of two (or more) different mechanisms to convert the incoming 511 keV photon to optical photons. The analytical result is validated with a simple 2D simulation.Main results.We show that if two different TOF-kernels are equivalent for image reconstruction with known attenuation, then they are also equivalent for joint reconstruction of the activity and the attenuation images. The variance increase in the activity, caused by also jointly reconstructing the attenuation image, vanishes when the TOF-resolution approaches perfection.Significance.These results are of interest for PET detector development and for the development of stand-alone PET systems.

Total-body PET: a new paradigm for molecular imaging

Total body (TB) positron emission tomography (PET) instruments have dramatically changed the paradigm of PET clinical and research studies due to their very high sensitivity and capability to image dynamic radiopharmaceutical distributions in the major organs of the body simultaneously. In this manuscript, we review the design of these systems and discuss general challenges and trade-offs to maximize the performance gains of current TB-PET systems. We then describe new concepts and technology that may impact future TB-PET systems. The manuscript summarizes what has been learned from the initial sites with TB-PET and explores potential research and clinical applications of TB-PET. The current generation of TB-PET systems range in axial field-of-view (AFOV) from 1 to 2 m and serve to illustrate the benefits and opportunities of a longer AFOV for various applications in PET. In only a few years of use these new TB-PET systems have shown that they will play an important role in expanding the field of molecular imaging and benefiting clinical practice.

Abbreviated scan protocols to capture 18F-FDG kinetics for long axial FOV PET scanners

PURPOSE

Kinetic parameters from dynamic 18F-fluorodeoxyglucose (FDG) imaging offer complementary insights to the study of disease compared to static clinical imaging. However, dynamic imaging protocols are cumbersome due to the long acquisition time. Long axial field-of-view (LAFOV) PET scanners (> 70 cm) have two advantages for dynamic imaging over clinical PET scanners with a standard axial field-of-view (SAFOV; 16-30 cm). The large axial coverage enables multi-organ dynamic imaging in a single bed position, and the high sensitivity may enable clinically routine abbreviated dynamic imaging protocols.

METHODS

In this work, we studied two abbreviated protocols using data from a 65-min dynamic 18F-FDG scan: (A) dynamic imaging immediately post-injection (p.i.) for variable durations, and (B) dynamic imaging immediately p.i. for variable durations plus a 1-h p.i. (5-min-long) datapoint. Nine cancer patients were imaged on the Biograph Vision Quadra (Siemens Healthineers). Time-activity curves over the lesions (N = 39) were fitted using the Patlak graphical analysis and a 2-tissue-compartment (2C, k4 = 0) model for variable scan durations (5-60 min). Kinetic parameters from the complete dataset served as the reference. Lesions from all cancers were grouped into low, medium, and high flux groups, and bias and precision of Ki (Patlak) and Ki, K1, k2, and k3 (2C) were calculated for each group.

RESULTS

Using only early dynamic data with the 2C (or Patlak) model, accurate quantification of Ki required at least 50 (or 55) min of dynamic data for low flux lesions, at least 30 (or 40) min for medium flux lesions, and at least 15 (or 20) min for high flux lesions to achieve both 10% bias and precision. The addition of the final (5-min) datapoint allowed for accurate quantification of Ki with a bias and precision of 10% using only 10-15 min of early dynamic data for either model.

CONCLUSION

Dynamic imaging for 10-15 min immediately p.i. followed by a 5-min scan at 1-h p.i can accurately and precisely quantify 18F-FDG on a long axial FOV scanner, potentially allowing for more widespread use of dynamic 18F-FDG imaging.

Initial evaluation of a new maximum-likelihood attenuation correction factor-based attenuation correction for time-of-flight brain PET

AIM

The aim of this study was to evaluate an image reconstruction algorithm, including a new maximum-likelihood attenuation correction factor (ML-ACF) for time of flight (TOF) brain positron emission tomography (PET).

METHODS

The implemented algorithm combines an ML-ACF method that simultaneously estimates both the emission image and attenuation sinogram from TOF emission data, and a scaling method based on anatomical features. To evaluate the algorithm's quantitative accuracy, three-dimensional brain phantom images were acquired and soft-tissue attenuation coefficients and emission values were analyzed.

RESULTS

The heterogeneous distributions of attenuation coefficients in soft tissue, skull, and nasal cavity were sufficiently visualized. The attenuation coefficient of soft tissue remained within 5% of theoretical value. Attenuation-corrected emission showed no lateral differences, and significant differences among soft tissue were within the error range.

CONCLUSION

The ML-ACF-based attenuation correction implemented for TOF brain PET worked well and obtained practical levels of accuracy.

Can a ToF-PET photon attenuation reconstruction test stopping-power estimations in proton therapy? A phantom study

Objective. The aim of the phantom study was to validate and to improve the computed tomography (CT) images used for the dose computation in proton therapy. It was tested, if the joint reconstruction of activity and attenuation images of time-of-flight PET (ToF-PET) scans could improve the estimation of the proton stopping-power.Approach. The attenuation images, i.e. CT images with 511 keV gamma-rays (γCTs), were jointly reconstructed with activity maps from ToF-PET scans. Theβ⁺activity was produced with FDG and in a separate experiment with proton-induced radioactivation. The phantoms contained slabs of tissue substitutes. The use of theγCTs for the prediction of the beam stopping in proton therapy was based on a linear relationship between theγ-ray attenuation, the electron density, and the stopping-power of fast protons.Main results. The FDG based experiment showed sufficient linearity to detect a bias of bony tissue in the heuristic look-up table, which maps between x-ray CT images and proton stopping-power.γCTs can be used for dose computation, if the electron density of one type of tissue is provided as a scaling factor. A possible limitation is imposed by the spatial resolution, which is inferior by a factor of 2.5 compared to the one of the x-ray CT.γCTs can also be derived from off-line, ToF-PET scans subsequent to the application of a proton field with a hypofractionated dose level.Significance. γCTs are a viable tool to support the estimation of proton stopping with radiotracer-based ToF-PET data from diagnosis or staging. This could be of higher potential relevance in MRI-guided proton therapy.γCTs could form an alternative approach to make use of in-beam or off-line PET scans of proton-inducedβ⁺activity with possible clinical limitations due to the low number of coincidence counts.

2D feasibility study of joint reconstruction of attenuation and activity in limited angle TOF-PET

Several research groups are studying organ-dedicated limited angle positron emission tomography (PET) systems to optimize performance-cost ratio, sensitivity, access to the patient and/or flexibility. Often open systems are considered, typically consisting of two detector panels of various sizes. Such systems provide incomplete sampling due to limited angular coverage and/or truncation, which leads to artefacts in the reconstructed activity images. In addition, these organ-dedicated PET systems are usually stand-alone systems, and as a result, no attenuation information can be obtained from anatomical images acquired in the same imaging session. It has been shown that the use of time-of-flight information reduces incomplete data artefacts and enables the joint estimation of the activity and the attenuation factors. In this work, we explore with simple 2D simulations the performance and stability of a joint reconstruction algorithm, for imaging with a limited angle PET system. The reconstruction is based on the so-called MLACF (Maximum Likelihood Attenuation Correction Factors) algorithm and uses linear attenuation coefficients in a known-tissue-class region to obtain absolute quantification. Different panel sizes and different time-of-flight (TOF) resolutions are considered. The noise propagation is compared to that of MLEM reconstruction with exact attenuation correction (AC) for the same PET system. The results show that with good TOF resolution, images of good visual quality can be obtained. If also a good scatter correction can be implemented, quantitative PET imaging will be possible. Further research, in particular on scatter correction, is required.

Modified kernel MLAA using autoencoder for PET-enabled dual-energy CT

Combined use of PET and dual-energy CT provides complementary information for multi-parametric imaging. PET-enabled dual-energy CT combines a low-energy X-ray CT image with a high-energy γ-ray CT (GCT) image reconstructed from time-of-flight PET emission data to enable dual-energy CT material decomposition on a PET/CT scanner. The maximum-likelihood attenuation and activity (MLAA) algorithm has been used for GCT reconstruction but suffers from noise. Kernel MLAA exploits an X-ray CT image prior through the kernel framework to guide GCT reconstruction and has demonstrated substantial improvements in noise suppression. However, similar to other kernel methods for image reconstruction, the existing kernel MLAA uses image intensity-based features to construct the kernel representation, which is not always robust and may lead to suboptimal reconstruction with artefacts. In this paper, we propose a modified kernel method by using an autoencoder convolutional neural network (CNN) to extract an intrinsic feature set from the X-ray CT image prior. A computer simulation study was conducted to compare the autoencoder CNN-derived feature representation with raw image patches for evaluation of kernel MLAA for GCT image reconstruction and dual-energy multi-material decomposition. The results show that the autoencoder kernel MLAA method can achieve a significant image quality improvement for GCT and material decomposition as compared to the existing kernel MLAA algorithm. A weakness of the proposed method is its potential over-smoothness in a bone region, indicating the importance of further optimization in future work. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Pitfalls on PET/CT Due to Artifacts and Instrumentation

PET/CT has become a preferred imaging modality over PET-only scanners in clinical practice. However, along with the significant improvement in diagnostic accuracy and patient throughput, pitfalls on PET/CT are reported as well. This review provides a general overview on the potential influence of the limitations with respect to PET/CT instrumentation and artifacts associated with the modality integration on the image appearance and quantitative accuracy of PET. Approaches proposed in literature to address the limitations or minimize the artifacts are discussed as well as their current challenges for clinical applications. Although the CT component can play an important role in assisting clinical diagnosis, we concentrate on the imaging scenarios where CT is used to provide auxiliary information for attenuation compensation and scatter correction in PET.

Data-driven respiratory phase-matched PET attenuation correction without CT

We propose a deep learning-based data-driven respiratory phase-matched gated-PET attenuation correction (AC) method that does not need a gated-CT. The proposed method is a multi-step process that consists of data-driven respiratory gating, gated attenuation map estimation using maximum-likelihood reconstruction of attenuation and activity (MLAA) algorithm, and enhancement of the gated attenuation maps using convolutional neural network (CNN). The gated MLAA attenuation maps enhanced by the CNN allowed for the phase-matched AC of gated-PET images. We conducted a non-rigid registration of the gated-PET images to generate motion-free PET images. We trained the CNN by conducting a 3D patch-based learning with 80 oncologic whole-body¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) PET/CT scan data and applied it to seven regional PET/CT scans that cover the lower lung and upper liver. We investigated the impact of the proposed respiratory phase-matched AC of PET without utilizing CT on tumor size and standard uptake value (SUV) assessment, and PET image quality (%STD). The attenuation corrected gated and motion-free PET images generated using the proposed method yielded sharper organ boundaries and better noise characteristics than conventional gated and ungated PET images. A banana artifact observed in a phase-mismatched CT-based AC was not observed in the proposed approach. By employing the proposed method, the size of tumor was reduced by 12.3% and SUV_90%was increased by 13.3% in tumors with larger movements than 5 mm. %STD of liver uptake was reduced by 11.1%. The deep learning-based data-driven respiratory phase-matched AC method improved the PET image quality and reduced the motion artifacts.

Update on the PennPET Explorer: A Whole-body Imager with Scalable Axial Field-of-View

Following successful performance testing and human imaging of a prototype PennPET Explorer, the scanner has been expanded to a current axial field of view of 1.12 m. Initial studies on this instrument have demonstrated encouraging results for total-body positron emission tomography imaging. Planned studies will test the capabilities of the PennPET Explorer further and inform the design of further human imaging protocols.

Practical joint reconstruction of activity and attenuation with autonomous scaling for time-of-flight PET

Recent research has showed that attenuation images can be determined from emission data, jointly with activity images, up to a scaling constant when utilizing the time-of-flight (TOF) information. We aim to develop practical CT-less joint reconstruction for clinical TOF PET scanners to obtain quantitatively accurate activity and attenuation images. In this work, we present a joint reconstruction of activity and attenuation based on MLAA (maximum likelihood reconstruction of attenuation and activity) with autonomous scaling determination and joint TOF scatter estimation from TOF PET data. Our idea for scaling is to use a selected volume of interest (VOI) in a reconstructed attenuation image with known attenuation, e.g. a liver in patient imaging. First, we construct a unit attenuation medium which has a similar, though not necessarily the same, support to the imaged emission object. All detectable LORs intersecting the unit medium have an attenuation factor of e -1≈ 0.3679, i.e. the line integral of linear attenuation coefficients is one. The scaling factor can then be determined from the difference between the reconstructed attenuation image and the known attenuation within the selected VOI normalized by the unit attenuation medium. A four-step iterative joint reconstruction algorithm is developed. In each iteration, (1) first the activity is updated using TOF OSEM from TOF list-mode data; (2) then the attenuation image is updated using XMLTR-a extended MLTR from non-TOF LOR sinograms; (3) a scaling factor is determined based on the selected VOI and both activity and attenuation images are updated using the estimated scaling; and (4) scatter is estimated using TOF single scatter simulation with the jointly reconstructed activity and attenuation images. The performance of joint reconstruction is studied using simulated data from a generic whole-body clinical TOF PET scanner and a long axial FOV research PET scanner as well as 3D experimental data from the PennPET Explorer scanner. We show that the proposed joint reconstruction with proper autonomous scaling provides low bias results comparable to the reference reconstruction with known attenuation.