Simultaneous reconstruction of activity and attenuation in time-of-flight PET

Optimized MLAA for quantitative non-TOF PET/MR of the brain

For quantitative tracer distribution in positron emission tomography, attenuation correction is essential. In a hybrid PET/CT system the CT images serve as a basis for generation of the attenuation map, but in PET/MR, the MR images do not have a similarly simple relationship with the attenuation map. Hence attenuation correction in PET/MR systems is more challenging. Typically either of two MR sequences are used: the Dixon or the ultra-short time echo (UTE) techniques. However these sequences have some well-known limitations. In this study, a reconstruction technique based on a modified and optimized non-TOF MLAA is proposed for PET/MR brain imaging. The idea is to tune the parameters of the MLTR applying some information from an attenuation image computed from the UTE sequences and a T1w MR image. In this MLTR algorithm, an [Formula: see text] parameter is introduced and optimized in order to drive the algorithm to a final attenuation map most consistent with the emission data. Because the non-TOF MLAA is used, a technique to reduce the cross-talk effect is proposed. In this study, the proposed algorithm is compared to the common reconstruction methods such as OSEM using a CT attenuation map, considered as the reference, and OSEM using the Dixon and UTE attenuation maps. To show the robustness and the reproducibility of the proposed algorithm, a set of 204 [¹⁸F]FDG patients, 35 [¹¹C]PiB patients and 1 [¹⁸F]FET patient are used. The results show that by choosing an optimized value of [Formula: see text] in MLTR, the proposed algorithm improves the results compared to the standard MR-based attenuation correction methods (i.e. OSEM using the Dixon or the UTE attenuation maps), and the cross-talk and the scale problem are limited.

Recent developments in time-of-flight PET

While the first time-of-flight (TOF)-positron emission tomography (PET) systems were already built in the early 1980s, limited clinical studies were acquired on these scanners. PET was still a research tool, and the available TOF-PET systems were experimental. Due to a combination of low stopping power and limited spatial resolution (caused by limited light output of the scintillators), these systems could not compete with bismuth germanate (BGO)-based PET scanners. Developments on TOF system were limited for about a decade but started again around 2000. The combination of fast photomultipliers, scintillators with high density, modern electronics, and faster computing power for image reconstruction have made it possible to introduce this principle in clinical TOF-PET systems. This paper reviews recent developments in system design, image reconstruction, corrections, and the potential in new applications for TOF-PET. After explaining the basic principles of time-of-flight, the difficulties in detector technology and electronics to obtain a good and stable timing resolution are shortly explained. The available clinical systems and prototypes under development are described in detail. The development of this type of PET scanner also requires modified image reconstruction with accurate modeling and correction methods. The additional dimension introduced by the time difference motivates a shift from sinogram- to listmode-based reconstruction. This reconstruction is however rather slow and therefore rebinning techniques specific for TOF data have been proposed. The main motivation for TOF-PET remains the large potential for image quality improvement and more accurate quantification for a given number of counts. The gain is related to the ratio of object size and spatial extent of the TOF kernel and is therefore particularly relevant for heavy patients, where image quality degrades significantly due to increased attenuation (low counts) and high scatter fractions. The original calculations for the gain were based on analytical methods. Recent publications for iterative reconstruction have shown that it is difficult to quantify TOF gain into one factor. The gain depends on the measured distribution, the location within the object, and the count rate. In a clinical situation, the gain can be used to either increase the standardized uptake value (SUV) or reduce the image acquisition time or administered dose. The localized nature of the TOF kernel makes it possible to utilize local tomography reconstruction or to separate emission from transmission data. The introduction of TOF also improves the joint estimation of transmission and emission images from emission data only. TOF is also interesting for new applications of PET-like isotopes with low branching ratio for positron fraction. The local nature also reduces the need for fine angular sampling, which makes TOF interesting for limited angle situations like breast PET and online dose imaging in proton or hadron therapy. The aim of this review is to introduce the reader in an educational way into the topic of TOF-PET and to give an overview of the benefits and new opportunities in using this additional information.

Combined PET/MRI: from Status Quo to Status Go. Summary Report of the Fifth International Workshop on PET/MR Imaging; February 15-19, 2016; Tübingen, Germany

This article provides a collaborative perspective of the discussions and conclusions from the fifth international workshop of combined positron emission tomorgraphy (PET)/magnetic resonance imaging (MRI) that was held in Tübingen, Germany, from February 15 to 19, 2016. Specifically, we summarise the second part of the workshop made up of invited presentations from active researchers in the field of PET/MRI and associated fields augmented by round table discussions and dialogue boards with specific topics. This year, this included practical advice as to possible approaches to moving PET/MRI into clinical routine, the use of PET/MRI in brain receptor imaging, in assessing cardiovascular diseases, cancer, infection, and inflammatory diseases. To address perceived challenges still remaining to innovatively integrate PET and MRI system technologies, a dedicated round table session brought together key representatives from industry and academia who were engaged with either the conceptualisation or early adoption of hybrid PET/MRI systems. Discussions during the workshop highlighted that emerging unique applications of PET/MRI such as the ability to provide multi-parametric quantitative and visual information which will enable not only overall disease detection but also disease characterisation would eventually be regarded as compelling arguments for the adoption of PET/MR. However, as indicated by previous workshops, evidence in favour of this observation is only growing slowly, mainly due to the ongoing inability to pool data cohorts from independent trials as well as different systems and sites. The participants emphasised that moving from status quo to status go entails the need to adopt standardised imaging procedures and the readiness to act together prospectively across multiple PET/MRI sites and vendors.

Simultaneous reconstruction of attenuation and activity in cardiac PET can remove CT misalignment artifacts

BACKGROUND

Misalignment between positron emission tomography (PET) and computed tomography (CT) data is known to generate artifactual defects in cardiac PET images due to imprecise attenuation correction (AC). In this work, the use of a maximum likelihood attenuation and activity (MLAA) algorithm is proposed to avoid such artifacts in time-of-flight (TOF) PET.

METHODS

MLAA was implemented and tested using a thorax/heart phantom and retrospectively on fourteen (13)N-ammonia PET/CT perfusion studies. Global and local misalignments between PET and CT data were generated by shifting matched CT images or using CT data representative of the end-inspiration phase. PET images were reconstructed with MLAA and a 3D-ordered-subsets-expectation-maximization (OSEM)-TOF algorithm. Images obtained with 3D-OSEM-TOF and matched CT were used as references. These images were compared (qualitatively and semi-quantitatively) with those reconstructed with 3D-OSEM-TOF and MLAA for which a misaligned CT was used, respectively, for AC and initialization.

RESULTS

Phantom experiment proved the capability of MLAA to converge toward the correct emission and attenuation distributions using, as input, only PET emission data, but convergence was very slow. Initializing MLAA with phantom CT images markedly improved convergence speed. In patient studies, when shifted or end-inspiration CT images were used for AC, 3D-OSEM-TOF reconstructions showed artifacts of increasing severity, size, and frequency with increasing mismatch. Such artifacts were absent in the corresponding MLAA images.

CONCLUSION

The proposed implementation of the MLAA algorithm is a feasible and robust technique to avoid AC mismatch artifacts in cardiac PET studies provided that a CT of the source is available, even if poorly aligned.

Generation of a Four-Class Attenuation Map for MRI-Based Attenuation Correction of PET Data in the Head Area Using a Novel Combination of STE/Dixon-MRI and FCM Clustering

PURPOSE

The aim of this study is to generate a four-class magnetic resonance imaging (MRI)-based attenuation map (μ-map) for attenuation correction of positron emission tomography (PET) data of the head area using a novel combination of short echo time (STE)/Dixon-MRI and a dedicated image segmentation method.

PROCEDURES

MR images of the head area were acquired using STE and two-point Dixon sequences. μ-maps were derived from MRI images based on a fuzzy C-means (FCM) clustering method along with morphologic operations. Quantitative assessment was performed to evaluate generated MRI-based μ-maps compared to X-ray computed tomography (CT)-based μ-maps.

RESULTS

The voxel-by-voxel comparison of MR-based and CT-based segmentation results yielded an average of more than 95 % for accuracy and specificity in the cortical bone, soft tissue, and air region. MRI-based μ-maps show a high correlation with those derived from CT scans (R (2) > 0.95).

CONCLUSIONS

Results indicate that STE/Dixon-MRI data in combination with FCM-based segmentation yields precise MR-based μ-maps for PET attenuation correction in hybrid PET/MRI systems.

Atlas-guided generation of pseudo-CT images for MRI-only and hybrid PET-MRI-guided radiotherapy treatment planning

Magnetic resonance imaging (MRI)-guided attenuation correction (AC) of positron emission tomography (PET) data and/or radiation therapy (RT) treatment planning is challenged by the lack of a direct link between MRI voxel intensities and electron density. Therefore, even if this is not a trivial task, a pseudo-computed tomography (CT) image must be predicted from MRI alone. In this work, we propose a two-step (segmentation and fusion) atlas-based algorithm focusing on bone tissue identification to create a pseudo-CT image from conventional MRI sequences and evaluate its performance against the conventional MRI segmentation technique and a recently proposed multi-atlas approach. The clinical studies consisted of pelvic CT, PET and MRI scans of 12 patients with loco-regionally advanced rectal disease. In the first step, bone segmentation of the target image is optimized through local weighted atlas voting. The obtained bone map is then used to assess the quality of deformed atlases to perform voxel-wise weighted atlas fusion. To evaluate the performance of the method, a leave-one-out cross-validation (LOOCV) scheme was devised to find optimal parameters for the model. Geometric evaluation of the produced pseudo-CT images and quantitative analysis of the accuracy of PET AC were performed. Moreover, a dosimetric evaluation of volumetric modulated arc therapy photon treatment plans calculated using the different pseudo-CT images was carried out and compared to those produced using CT images serving as references. The pseudo-CT images produced using the proposed method exhibit bone identification accuracy of 0.89 based on the Dice similarity metric compared to 0.75 achieved by the other atlas-based method. The superior bone extraction resulted in a mean standard uptake value bias of  -1.5  ±  5.0% (mean  ±  SD) in bony structures compared to  -19.9  ±  11.8% and  -8.1  ±  8.2% achieved by MRI segmentation-based (water-only) and atlas-guided AC. Dosimetric evaluation using dose volume histograms and the average difference between minimum/maximum absorbed doses revealed a mean error of less than 1% for the both target volumes and organs at risk. Two-dimensional (2D) gamma analysis of the isocenter dose distributions at 1%/1 mm criterion revealed pass rates of 91.40  ±  7.56%, 96.00  ±  4.11% and 97.67  ±  3.6% for MRI segmentation, atlas-guided and the proposed methods, respectively. The proposed method generates accurate pseudo-CT images from conventional Dixon MRI sequences with improved bone extraction accuracy. The approach is promising for potential use in PET AC and MRI-only or hybrid PET/MRI-guided RT treatment planning.

One registration multi-atlas-based pseudo-CT generation for attenuation correction in PET/MRI

PURPOSE

The outcome of a detailed assessment of various strategies for atlas-based whole-body bone segmentation from magnetic resonance imaging (MRI) was exploited to select the optimal parameters and setting, with the aim of proposing a novel one-registration multi-atlas (ORMA) pseudo-CT generation approach.

METHODS

The proposed approach consists of only one online registration between the target and reference images, regardless of the number of atlas images (N), while for the remaining atlas images, the pre-computed transformation matrices to the reference image are used to align them to the target image. The performance characteristics of the proposed method were evaluated and compared with conventional atlas-based attenuation map generation strategies (direct registration of the entire atlas images followed by voxel-wise weighting (VWW) and arithmetic averaging atlas fusion). To this end, four different positron emission tomography (PET) attenuation maps were generated via arithmetic averaging and VWW scheme using both direct registration and ORMA approaches as well as the 3-class attenuation map obtained from the Philips Ingenuity TF PET/MRI scanner commonly used in the clinical setting. The evaluation was performed based on the accuracy of extracted whole-body bones by the different attenuation maps and by quantitative analysis of resulting PET images compared to CT-based attenuation-corrected PET images serving as reference.

RESULTS

The comparison of validation metrics regarding the accuracy of extracted bone using the different techniques demonstrated the superiority of the VWW atlas fusion algorithm achieving a Dice similarity measure of 0.82 ± 0.04 compared to arithmetic averaging atlas fusion (0.60 ± 0.02), which uses conventional direct registration. Application of the ORMA approach modestly compromised the accuracy, yielding a Dice similarity measure of 0.76 ± 0.05 for ORMA-VWW and 0.55 ± 0.03 for ORMA-averaging. The results of quantitative PET analysis followed the same trend with less significant differences in terms of SUV bias, whereas massive improvements were observed compared to PET images corrected for attenuation using the 3-class attenuation map. The maximum absolute bias achieved by VWW and VWW-ORMA methods was 06.4 ± 5.5 in the lung and 07.9 ± 4.8 in the bone, respectively.

CONCLUSIONS

The proposed algorithm is capable of generating decent attenuation maps. The quantitative analysis revealed a good correlation between PET images corrected for attenuation using the proposed pseudo-CT generation approach and the corresponding CT images. The computational time is reduced by a factor of 1/N at the expense of a modest decrease in quantitative accuracy, thus allowing us to achieve a reasonable compromise between computing time and quantitative performance.

Comparison between MRI-based attenuation correction methods for brain PET in dementia patients

INTRODUCTION

The combination of Positron Emission Tomography (PET) with magnetic resonance imaging (MRI) in hybrid PET/MRI scanners offers a number of advantages in investigating brain structure and function. A critical step of PET data reconstruction is attenuation correction (AC). Accounting for bone in attenuation maps (μ-map) was shown to be important in brain PET studies. While there are a number of MRI-based AC methods, no systematic comparison between them has been performed so far. The aim of this work was to study the different performance obtained by some of the recent methods presented in the literature. To perform such a comparison, we focused on [¹⁸F]-Fluorodeoxyglucose-PET/MRI neurodegenerative dementing disorders, which are known to exhibit reduced levels of glucose metabolism in certain brain regions.

METHODS

Four novel methods were used to calculate μ-maps from MRI data of 15 patients with Alzheimer's dementia (AD). The methods cover two atlas-based methods, a segmentation method, and a hybrid template/segmentation method. Additionally, the Dixon-based and a UTE-based method, offered by a vendor, were included in the comparison. Performance was assessed at three levels: tissue identification accuracy in the μ-map, quantitative accuracy of reconstructed PET data in specific brain regions, and precision in diagnostic images at identifying hypometabolic areas.

RESULTS

Quantitative regional errors of -20--10 % were obtained using the vendor's AC methods, whereas the novel methods produced errors in a margin of ±5 %. The obtained precision at identifying areas with abnormally low levels of glucose uptake, potentially regions affected by AD, were 62.9 and 79.5 % for the two vendor AC methods, the former ignoring bone and the latter including bone information. The precision increased to 87.5-93.3 % in average for the four new methods, exhibiting similar performances.

CONCLUSION

We confirm that the AC methods based on the Dixon and UTE sequences provided by the vendor are inferior to alternative techniques. As a novel finding, there was no substantial difference between the recently proposed atlas-based, template-based and segmentation-based methods.

Current image acquisition options in PET/MR

Whole-body PET/MR hybrid imaging combines excellent soft tissue contrast and various functional imaging parameters provided by MR with high sensitivity and quantification of radiotracer uptake provided by PET. Although clinical evaluation now is under way, PET/MR demands for new technologies and innovative solutions, currently subject to interdisciplinary research. Attenuation correction (AC) of human soft tissues and of hardware components has to be MR based to maintain quantification of PET imaging as CT attenuation information is missing. MR-based AC is inherently associated with the following challenges: patient tissues are segmented into only few tissue classes, providing discrete attenuation coefficients; bone is substituted as soft tissue in MR-based AC; the limited field of view in MRI leads to truncations in body imaging and, consequently, in MR-based AC; and correct segmentation of lung tissue may be hampered by breathing artifacts. Use of time of flight during PET image acquisition and reconstruction, however, may improve the accuracy of AC. This article provides a status of current image acquisition options in PET/MR hybrid imaging.

Combined PET/MR: The Real Work Has Just Started. Summary Report of the Third International Workshop on PET/MR Imaging; February 17-21, 2014, Tübingen, Germany

This paper summarises the proceedings and discussions at the third annual workshop held in Tübingen, Germany, dedicated to the advancement of the technical, scientific and clinical applications of combined PET/MRI systems in humans. Two days of basic scientific and technical instructions with "hands-on" tutorials were followed by 3 days of invited presentations from active researchers in this and associated fields augmented by round-table discussions and dialogue boards with specific themes. These included the use of PET/MRI in paediatric oncology and in adult neurology, oncology and cardiology, the development of multi-parametric analyses, and efforts to standardise PET/MRI examinations to allow pooling of data for evaluating the technology. A poll taken on the final day demonstrated that over 50 % of those present felt that while PET/MRI technology underwent an inevitable slump after its much-anticipated initial launch, it was now entering a period of slow, progressive development, with new key applications emerging. In particular, researchers are focusing on exploiting the complementary nature of the physiological (PET) and biochemical (MRI/MRS) data within the morphological framework (MRI) that these devices can provide. Much of the discussion was summed up on the final day when one speaker commented on the state of PET/MRI: "the real work has just started".

The Use of Anatomical Information for Molecular Image Reconstruction Algorithms: Attenuation/Scatter Correction, Motion Compensation, and Noise Reduction

PET and SPECT are important tools for providing valuable molecular information about patients to clinicians. Advances in nuclear medicine hardware technologies and statistical image reconstruction algorithms enabled significantly improved image quality. Sequentially or simultaneously acquired anatomical images such as CT and MRI from hybrid scanners are also important ingredients for improving the image quality of PET or SPECT further. High-quality anatomical information has been used and investigated for attenuation and scatter corrections, motion compensation, and noise reduction via post-reconstruction filtering and regularization in inverse problems. In this article, we will review works using anatomical information for molecular image reconstruction algorithms for better image quality by describing mathematical models, discussing sources of anatomical information for different cases, and showing some examples.

Quantitative analysis of MRI-guided attenuation correction techniques in time-of-flight brain PET/MRI

PURPOSE

In quantitative PET/MR imaging, attenuation correction (AC) of PET data is markedly challenged by the need of deriving accurate attenuation maps from MR images. A number of strategies have been developed for MRI-guided attenuation correction with different degrees of success. In this work, we compare the quantitative performance of three generic AC methods, including standard 3-class MR segmentation-based, advanced atlas-registration-based and emission-based approaches in the context of brain time-of-flight (TOF) PET/MRI.

MATERIALS AND METHODS

Fourteen patients referred for diagnostic MRI and (18)F-FDG PET/CT brain scans were included in this comparative study. For each study, PET images were reconstructed using four different attenuation maps derived from CT-based AC (CTAC) serving as reference, standard 3-class MR-segmentation, atlas-registration and emission-based AC methods. To generate 3-class attenuation maps, T1-weighted MRI images were segmented into background air, fat and soft-tissue classes followed by assignment of constant linear attenuation coefficients of 0, 0.0864 and 0.0975 cm(-1) to each class, respectively. A robust atlas-registration based AC method was developed for pseudo-CT generation using local weighted fusion of atlases based on their morphological similarity to target MR images. Our recently proposed MRI-guided maximum likelihood reconstruction of activity and attenuation (MLAA) algorithm was employed to estimate the attenuation map from TOF emission data. The performance of the different AC algorithms in terms of prediction of bones and quantification of PET tracer uptake was objectively evaluated with respect to reference CTAC maps and CTAC-PET images.

RESULTS

Qualitative evaluation showed that the MLAA-AC method could sparsely estimate bones and accurately differentiate them from air cavities. It was found that the atlas-AC method can accurately predict bones with variable errors in defining air cavities. Quantitative assessment of bone extraction accuracy based on Dice similarity coefficient (DSC) showed that MLAA-AC and atlas-AC resulted in DSC mean values of 0.79 and 0.92, respectively, in all patients. The MLAA-AC and atlas-AC methods predicted mean linear attenuation coefficients of 0.107 and 0.134 cm(-1), respectively, for the skull compared to reference CTAC mean value of 0.138cm(-1). The evaluation of the relative change in tracer uptake within 32 distinct regions of the brain with respect to CTAC PET images showed that the 3-class MRAC, MLAA-AC and atlas-AC methods resulted in quantification errors of -16.2 ± 3.6%, -13.3 ± 3.3% and 1.0 ± 3.4%, respectively. Linear regression and Bland-Altman concordance plots showed that both 3-class MRAC and MLAA-AC methods result in a significant systematic bias in PET tracer uptake, while the atlas-AC method results in a negligible bias.

CONCLUSION

The standard 3-class MRAC method significantly underestimated cerebral PET tracer uptake. While current state-of-the-art MLAA-AC methods look promising, they were unable to noticeably reduce quantification errors in the context of brain imaging. Conversely, the proposed atlas-AC method provided the most accurate attenuation maps, and thus the lowest quantification bias.

Attenuation correction in emission tomography using the emission data--A review

The problem of attenuation correction (AC) for quantitative positron emission tomography (PET) had been considered solved to a large extent after the commercial availability of devices combining PET with computed tomography (CT) in 2001; single photon emission computed tomography (SPECT) has seen a similar development. However, stimulated in particular by technical advances toward clinical systems combining PET and magnetic resonance imaging (MRI), research interest in alternative approaches for PET AC has grown substantially in the last years. In this comprehensive literature review, the authors first present theoretical results with relevance to simultaneous reconstruction of attenuation and activity. The authors then look back at the early history of this research area especially in PET; since this history is closely interwoven with that of similar approaches in SPECT, these will also be covered. We then review algorithmic advances in PET, including analytic and iterative algorithms. The analytic approaches are either based on the Helgason-Ludwig data consistency conditions of the Radon transform, or generalizations of John's partial differential equation; with respect to iterative methods, we discuss maximum likelihood reconstruction of attenuation and activity (MLAA), the maximum likelihood attenuation correction factors (MLACF) algorithm, and their offspring. The description of methods is followed by a structured account of applications for simultaneous reconstruction techniques: this discussion covers organ-specific applications, applications specific to PET/MRI, applications using supplemental transmission information, and motion-aware applications. After briefly summarizing SPECT applications, we consider recent developments using emission data other than unscattered photons. In summary, developments using time-of-flight (TOF) PET emission data for AC have shown promising advances and open a wide range of applications. These techniques may both remedy deficiencies of purely MRI-based AC approaches in PET/MRI and improve standalone PET imaging.

MR Imaging-Guided Attenuation Correction of PET Data in PET/MR Imaging

Attenuation correction (AC) is one of the most important challenges in the recently introduced combined PET/magnetic resonance (MR) scanners. PET/MR AC (MR-AC) approaches aim to develop methods that allow accurate estimation of the linear attenuation coefficients of the tissues and other components located in the PET field of view. MR-AC methods can be divided into 3 categories: segmentation, atlas, and PET based. This review provides a comprehensive list of the state-of-the-art MR-AC approaches and their pros and cons. The main sources of artifacts are presented. Finally, this review discusses the current status of MR-AC approaches for clinical applications.

Maximum-likelihood joint image reconstruction and motion estimation with misaligned attenuation in TOF-PET/CT

This work is an extension of our recent work on joint activity reconstruction/motion estimation (JRM) from positron emission tomography (PET) data. We performed JRM by maximization of the penalized log-likelihood in which the probabilistic model assumes that the same motion field affects both the activity distribution and the attenuation map. Our previous results showed that JRM can successfully reconstruct the activity distribution when the attenuation map is misaligned with the PET data, but converges slowly due to the significant cross-talk in the likelihood. In this paper, we utilize time-of-flight PET for JRM and demonstrate that the convergence speed is significantly improved compared to JRM with conventional PET data.

Advances in time-of-flight PET

This paper provides a review and an update on time-of-flight PET imaging with a focus on PET instrumentation, ranging from hardware design to software algorithms. We first present a short introduction to PET, followed by a description of TOF PET imaging and its history from the early days. Next, we introduce the current state-of-art in TOF PET technology and briefly summarize the benefits of TOF PET imaging. This is followed by a discussion of the various technological advancements in hardware (scintillators, photo-sensors, electronics) and software (image reconstruction) that have led to the current widespread use of TOF PET technology, and future developments that have the potential for further improvements in the TOF imaging performance. We conclude with a discussion of some new research areas that have opened up in PET imaging as a result of having good system timing resolution, ranging from new algorithms for attenuation correction, through efficient system calibration techniques, to potential for new PET system designs.

Maximum-Likelihood Joint Image Reconstruction/Motion Estimation in Attenuation-Corrected Respiratory Gated PET/CT Using a Single Attenuation Map

This work provides an insight into positron emission tomography (PET) joint image reconstruction/motion estimation (JRM) by maximization of the likelihood, where the probabilistic model accounts for warped attenuation. Our analysis shows that maximum-likelihood (ML) JRM returns the same reconstructed gates for any attenuation map (μ-map) that is a deformation of a given μ-map, regardless of its alignment with the PET gates. We derived a joint optimization algorithm accordingly, and applied it to simulated and patient gated PET data. We first evaluated the proposed algorithm on simulations of respiratory gated PET/CT data based on the XCAT phantom. Our results show that independently of which μ-map is used as input to JRM: (i) the warped μ-maps correspond to the gated μ-maps, (ii) JRM outperforms the traditional post-registration reconstruction and consolidation (PRRC) for hot lesion quantification and (iii) reconstructed gated PET images are similar to those obtained with gated μ-maps. This suggests that a breath-held μ-map can be used. We then applied JRM on patient data with a μ-map derived from a breath-held high resolution CT (HRCT), and compared the results with PRRC, where each reconstructed PET image was obtained with a corresponding cine-CT gated μ-map. Results show that JRM with breath-held HRCT achieves similar reconstruction to that using PRRC with cine-CT. This suggests a practical low-dose solution for implementation of motion-corrected respiratory gated PET/CT.

Sensitivity estimation in time-of-flight list-mode positron emission tomography

PURPOSE

An accurate quantification of the images in positron emission tomography (PET) requires knowing the actual sensitivity at each voxel, which represents the probability that a positron emitted in that voxel is finally detected as a coincidence of two gamma rays in a pair of detectors in the PET scanner. This sensitivity depends on the characteristics of the acquisition, as it is affected by the attenuation of the annihilation gamma rays in the body, and possible variations of the sensitivity of the scanner detectors. In this work, the authors propose a new approach to handle time-of-flight (TOF) list-mode PET data, which allows performing either or both, a self-attenuation correction, and self-normalization correction based on emission data only.

METHODS

The authors derive the theory using a fully Bayesian statistical model of complete data. The authors perform an initial evaluation of algorithms derived from that theory and proposed in this work using numerical 2D list-mode simulations with different TOF resolutions and total number of detected coincidences. Effects of randoms and scatter are not simulated.

RESULTS

The authors found that proposed algorithms successfully correct for unknown attenuation and scanner normalization for simulated 2D list-mode TOF-PET data.

CONCLUSIONS

A new method is presented that can be used for corrections for attenuation and normalization (sensitivity) using TOF list-mode data.

Joint Estimation of Activity and Attenuation in Whole-Body TOF PET/MRI Using Constrained Gaussian Mixture Models

It has recently been shown that the attenuation map can be estimated from time-of-flight (TOF) PET emission data using joint maximum likelihood reconstruction of attenuation and activity (MLAA). In this work, we propose a novel MRI-guided MLAA algorithm for emission-based attenuation correction in whole-body PET/MR imaging. The algorithm imposes MR spatial and CT statistical constraints on the MLAA estimation of attenuation maps using a constrained Gaussian mixture model (GMM) and a Markov random field smoothness prior. Dixon water and fat MR images were segmented into outside air, lung, fat and soft-tissue classes and an MR low-intensity (unknown) class corresponding to air cavities, cortical bone and susceptibility artifacts. The attenuation coefficients over the unknown class were estimated using a mixture of four Gaussians, and those over the known tissue classes using unimodal Gaussians, parameterized over a patient population. To eliminate misclassification of spongy bones with surrounding tissues, and thus include them in the unknown class, we heuristically suppressed fat in water images and also used a co-registered bone probability map. The proposed MLAA-GMM algorithm was compared with the MLAA algorithms proposed by Rezaei and Salomon using simulation and clinical studies with two different tracer distributions. The results showed that our proposed algorithm outperforms its counterparts in suppressing the cross-talk and scaling problems of activity and attenuation and thus produces PET images of improved quantitative accuracy. It can be concluded that the proposed algorithm effectively exploits the MR information and can pave the way toward accurate emission-based attenuation correction in TOF PET/MRI.

Transmission-less attenuation estimation from time-of-flight PET histo-images using consistency equations

In positron emission tomography (PET) imaging, attenuation correction with accurate attenuation estimation is crucial for quantitative patient studies. Recent research showed that the attenuation sinogram can be determined up to a scaling constant utilizing the time-of-flight information. The TOF-PET data can be naturally and efficiently stored in a histo-image without information loss, and the radioactive tracer distribution can be efficiently reconstructed using the DIRECT approaches. In this paper, we explore transmission-less attenuation estimation from TOF-PET histo-images. We first present the TOF-PET histo-image formation and the consistency equations in the histo-image parameterization, then we derive a least-squares solution for estimating the directional derivatives of the attenuation factors from the measured emission histo-images. Finally, we present a fast solver to estimate the attenuation factors from their directional derivatives using the discrete sine transform and fast Fourier transform while considering the boundary conditions. We find that the attenuation histo-images can be uniquely determined from the TOF-PET histo-images by considering boundary conditions. Since the estimate of the attenuation directional derivatives can be inaccurate for LORs tangent to the patient boundary, external sources, e.g. a ring or annulus source, might be needed to give an accurate estimate of the attenuation gradient for such LORs. The attenuation estimation from TOF-PET emission histo-images is demonstrated using simulated 2D TOF-PET data.

Emission-based estimation of lung attenuation coefficients for attenuation correction in time-of-flight PET/MR

In standard segmentation-based MRI-guided attenuation correction (MRAC) of PET data on hybrid PET/MRI systems, the inter/intra-patient variability of linear attenuation coefficients (LACs) is ignored owing to the assignment of a constant LAC to each tissue class. This can lead to PET quantification errors, especially in the lung regions. In this work, we aim to derive continuous and patient-specific lung LACs from time-of-flight (TOF) PET emission data using the maximum likelihood reconstruction of activity and attenuation (MLAA) algorithm. The MLAA algorithm was constrained for estimation of lung LACs only in the standard 4-class MR attenuation map using Gaussian lung tissue preference and Markov random field smoothness priors. MRAC maps were derived from segmentation of CT images of 19 TOF-PET/CT clinical studies into background air, lung, soft tissue and fat tissue classes, followed by assignment of predefined LACs of 0, 0.0224, 0.0864 and 0.0975 cm(-1), respectively. The lung LACs of the resulting attenuation maps were then estimated from emission data using the proposed MLAA algorithm. PET quantification accuracy of MRAC and MLAA methods was evaluated against the reference CT-based AC method in the lungs, lesions located in/near the lungs and neighbouring tissues. The results show that the proposed MLAA algorithm is capable of retrieving lung density gradients and compensate fairly for respiratory-phase mismatch between PET and corresponding attenuation maps. It was found that the mean of the estimated lung LACs generally follow the trend of the reference CT-based attenuation correction (CTAC) method. Quantitative analysis revealed that the MRAC method resulted in average relative errors of -5.2 ± 7.1% and -6.1 ± 6.7% in the lungs and lesions, respectively. These were reduced by the MLAA algorithm to -0.8 ± 6.3% and -3.3 ± 4.7%, respectively. In conclusion, we demonstrated the potential and capability of emission-based methods in deriving patient-specific lung LACs to improve the accuracy of attenuation correction in TOF PET/MR imaging, thus paving the way for their adaptation in the clinic.

Clinical Assessment of Emission- and Segmentation-Based MR-Guided Attenuation Correction in Whole-Body Time-of-Flight PET/MR Imaging

The joint maximum-likelihood reconstruction of activity and attenuation (MLAA) for emission-based attenuation correction has regained attention since the advent of time-of-flight PET/MR imaging. Recently, we improved the performance of the MLAA algorithm using an MR imaging-constrained gaussian mixture model (GMM). In this study, we compare the performance of our proposed algorithm with standard 4-class MR-based attenuation correction (MRAC) implemented on commercial systems.

METHODS

Five head and neck (18)F-FDG patients were scanned on PET/MR imaging and PET/CT scanners. Dixon fat and water MR images were registered to CT images. MRAC maps were derived by segmenting the MR images into 4 tissue classes and assigning predefined attenuation coefficients. For MLAA-GMM, MR images were segmented into known tissue classes, including fat, soft tissue, lung, background air, and an unknown MR low-intensity class encompassing cortical bones, air cavities, and metal artifacts. A coregistered bone probability map was also included in the unknown tissue class. Finally, the GMM prior was constrained over known tissue classes of attenuation maps using unimodal gaussians parameterized over a patient population.

RESULTS

The results showed that the MLAA-GMM algorithm outperformed the MRAC method by differentiating bones from air gaps and providing more accurate patient-specific attenuation coefficients of soft tissue and lungs. It was found that the MRAC and MLAA-GMM methods resulted in average standardized uptake value errors of -5.4% and -3.5% in the lungs, -7.4% and -5.0% in soft tissues/lesions, and -18.4% and -10.2% in bones, respectively.

CONCLUSION

The proposed MLAA algorithm is promising for accurate derivation of attenuation maps on time-of-flight PET/MR systems.

MR-based attenuation correction using ultrashort-echo-time pulse sequences in dementia patients

Attenuation correction (AC) is a critical requirement for quantitative PET reconstruction. Accounting for bone information in the attenuation map (μ map) is of paramount importance for accurate brain PET quantification. However, to measure the signal from bone structures represents a challenging task in MR. Recent (18)F-FDG PET/MR studies showed quantitative bias for the assessment of radiotracer concentration when bone was ignored. This work is focused on (18)F-FDG PET/MR neurodegenerative dementing disorders. These are known to lead to specific patterns of (18)F-FDG hypometabolism, mainly in superficial brain structures, which might suffer from attenuation artifacts and thus have immediate diagnostic consequences. A fully automatic method to estimate the μ map, including bone tissue using only MR information, is presented.

METHODS

The algorithm was based on a dual-echo ultrashort-echo-time MR imaging sequence to calculate the R2 map, from which the μ map was derived. The R2-based μ map was postprocessed to calculate an estimated distribution of the bone tissue. μ maps calculated from datasets of 9 patients were compared with their CT-based μ maps (μ mapCT) by determining the confusion matrix. Additionally, a region-of-interest comparison between reconstructed PET data, corrected using different μ maps, was performed. PET data were reconstructed using a Dixon-based μ map (μ mapDX) and a dual-echo ultrashort-echo-time-based μ map (μ mapUTE), which are both calculated by the scanner, and the R2-based μ map presented in this work was compared with reconstructed PET data using the μ mapCT as a reference.

RESULTS

Errors were approximately 20% higher using the μ mapDX and μ mapUTE for AC, compared with reconstructed PET data using the reference μ mapCT. However, PET AC using the R2-based μ map resulted, for all the patients and all the analyzed regions of interest, in a significant improvement, reducing the error to -5.8% to 2.5%.

CONCLUSION

The proposed method successfully showed significantly reduced errors in quantification, compared with the μ mapDX and μ mapUTE, and therefore delivered more accurate PET image quantification for an improved diagnostic workup in dementia patients.

PET-MRI: a review of challenges and solutions in the development of integrated multimodality imaging

The integration of positron emission tomography (PET) and magnetic resonance imaging (MRI) has been an ongoing research topic for the last 20 years. This paper gives an overview of the different developments and the technical problems associated with combining PET and MRI in one system. After explaining the different detector concepts for integrating PET-MRI and minimising interference the limitations and advantages of different solutions for the detector and system are described for preclinical and clinical imaging systems. The different integrated PET-MRI systems are described in detail. Besides detector concepts and system integration the challenges and proposed solutions for attenuation correction and the potential for motion correction and resolution recovery are also discussed in this topical review.

The dynamics of physics in PET

From a technical perspective, there are fundamentally two forces driving the evolution of instrumentation in positron emission tomography (PET) and nuclear medicine generally: clinical needs and technical innovation. This essay considers some of the dynamics of these forces as they act on physics-related developments in PET and suggests that progress will be greatest if these differing motivations are kept in balance as the field evolves.

Attenuation correction synthesis for hybrid PET-MR scanners: application to brain studies

Attenuation correction is an essential requirement for quantification of positron emission tomography (PET) data. In PET/CT acquisition systems, attenuation maps are derived from computed tomography (CT) images. However, in hybrid PET/MR scanners, magnetic resonance imaging (MRI) images do not directly provide a patient-specific attenuation map. The aim of the proposed work is to improve attenuation correction for PET/MR scanners by generating synthetic CTs and attenuation maps. The synthetic images are generated through a multi-atlas information propagation scheme, locally matching the MRI-derived patient's morphology to a database of MRI/CT pairs, using a local image similarity measure. Results show significant improvements in CT synthesis and PET reconstruction accuracy when compared to a segmentation method using an ultrashort-echo-time MRI sequence and to a simplified atlas-based method.

Quantitative simultaneous positron emission tomography and magnetic resonance imaging

Simultaneous positron emission tomography and magnetic resonance imaging (PET-MR) is an innovative and promising imaging modality that is generating substantial interest in the medical imaging community, while offering many challenges and opportunities. In this study, we investigated whether MR surface coils need to be accounted for in PET attenuation correction. Furthermore, we integrated motion correction, attenuation correction, and point spread function modeling into a single PET reconstruction framework. We applied our reconstruction framework to in vivo animal and patient PET-MR studies. We have demonstrated that our approach greatly improved PET image quality.

An SPM8-based approach for attenuation correction combining segmentation and nonrigid template formation: application to simultaneous PET/MR brain imaging

We present an approach for head MR-based attenuation correction (AC) based on the Statistical Parametric Mapping 8 (SPM8) software, which combines segmentation- and atlas-based features to provide a robust technique to generate attenuation maps (μ maps) from MR data in integrated PET/MR scanners.

METHODS

Coregistered anatomic MR and CT images of 15 glioblastoma subjects were used to generate the templates. The MR images from these subjects were first segmented into 6 tissue classes (gray matter, white matter, cerebrospinal fluid, bone, soft tissue, and air), which were then nonrigidly coregistered using a diffeomorphic approach. A similar procedure was used to coregister the anatomic MR data for a new subject to the template. Finally, the CT-like images obtained by applying the inverse transformations were converted to linear attenuation coefficients to be used for AC of PET data. The method was validated on 16 new subjects with brain tumors (n = 12) or mild cognitive impairment (n = 4) who underwent CT and PET/MR scans. The μ maps and corresponding reconstructed PET images were compared with those obtained using the gold standard CT-based approach and the Dixon-based method available on the Biograph mMR scanner. Relative change (RC) images were generated in each case, and voxel- and region-of-interest-based analyses were performed.

RESULTS

The leave-one-out cross-validation analysis of the data from the 15 atlas-generation subjects showed small errors in brain linear attenuation coefficients (RC, 1.38% ± 4.52%) compared with the gold standard. Similar results (RC, 1.86% ± 4.06%) were obtained from the analysis of the atlas-validation datasets. The voxel- and region-of-interest-based analysis of the corresponding reconstructed PET images revealed quantification errors of 3.87% ± 5.0% and 2.74% ± 2.28%, respectively. The Dixon-based method performed substantially worse (the mean RC values were 13.0% ± 10.25% and 9.38% ± 4.97%, respectively). Areas closer to the skull showed the largest improvement.

CONCLUSION

We have presented an SPM8-based approach for deriving the head μ map from MR data to be used for PET AC in integrated PET/MR scanners. Its implementation is straightforward and requires only the morphologic data acquired with a single MR sequence. The method is accurate and robust, combining the strengths of both segmentation- and atlas-based approaches while minimizing their drawbacks.

ML-reconstruction for TOF-PET with simultaneous estimation of the attenuation factors

In positron emission tomography (PET), attenuation correction is typically done based on information obtained from transmission tomography. Recent studies show that time-of-flight (TOF) PET emission data allow joint estimation of activity and attenuation images. Mathematical analysis revealed that the joint estimation problem is determined up to a scale factor. In this work, we propose a maximum likelihood reconstruction algorithm that jointly estimates the activity image together with the sinogram of the attenuation factors. The algorithm is evaluated with 2-D and 3-D simulations as well as clinical TOF-PET measurements of a patient scan and compared to reference reconstructions. The robustness of the algorithm to possible imperfect scanner calibration is demonstrated with reconstructions of the patient scan ignoring the varying detector sensitivities.

The contribution of physics to Nuclear Medicine: physicians' perspective on future directions

Background

Advances in Nuclear Medicine physics enabled the specialty of Nuclear Medicine and directed research in other aspects of radiotracer imaging, ultimately leading to Nuclear Medicine’s emergence as an important component of current medical practice.

Discussion

Nuclear Medicine’s unique ability to characterize in vivo biology without perturbing it will assure its ongoing role in a practice of medicine increasingly driven by molecular biology. However, in the future, it is likely that advances in molecular biology and radiopharmaceutical chemistry will increasingly direct future developments in Nuclear Medicine physics, rather than relying on physics as the primary driver of advances in Nuclear Medicine.

Summary

Working hand-in-hand with clinicians, chemists, and biologists, Nuclear Medicine physicists can greatly enhance the specialty by creating more sensitive and robust imaging devices, by enabling more facile and sophisticated image analysis to yield quantitative measures of regional in vivo biology, and by combining the strengths of radiotracer imaging with other imaging modalities in hybrid devices, with the overall goal to enhance Nuclear Medicine’s ability to characterize regional in vivo biology.

Transmission-less attenuation correction in time-of-flight PET: analysis of a discrete iterative algorithm

The maximum likelihood attenuation correction factors (MLACF) algorithm has been developed to calculate the maximum-likelihood estimate of the activity image and the attenuation sinogram in time-of-flight (TOF) positron emission tomography, using only emission data without prior information on the attenuation. We consider the case of a Poisson model of the data, in the absence of scatter or random background. In this case the maximization with respect to the attenuation factors can be achieved in a closed form and the MLACF algorithm works by updating the activity. Despite promising numerical results, the convergence of this algorithm has not been analysed. In this paper we derive the algorithm and demonstrate that the MLACF algorithm monotonically increases the likelihood, is asymptotically regular, and that the limit points of the iteration are stationary points of the likelihood. Because the problem is not convex, however, the limit points might be saddle points or local maxima. To obtain some empirical insight into the latter question, we present data obtained by applying MLACF to 2D simulated TOF data, using a large number of iterations and different initializations.

List-mode reconstruction for the Biograph mCT with physics modeling and event-by-event motion correction

Whole-body PET/CT scanners are important clinical and research tools to study tracer distribution throughout the body. In whole-body studies, respiratory motion results in image artifacts. We have previously demonstrated for brain imaging that, when provided with accurate motion data, event-by-event correction has better accuracy than frame-based methods. Therefore, the goal of this work was to develop a list-mode reconstruction with novel physics modeling for the Siemens Biograph mCT with event-by-event motion correction, based on the MOLAR platform (Motion-compensation OSEM List-mode Algorithm for Resolution-Recovery Reconstruction). Application of MOLAR for the mCT required two algorithmic developments. First, in routine studies, the mCT collects list-mode data in 32 bit packets, where averaging of lines-of-response (LORs) by axial span and angular mashing reduced the number of LORs so that 32 bits are sufficient to address all sinogram bins. This degrades spatial resolution. In this work, we proposed a probabilistic LOR (pLOR) position technique that addresses axial and transaxial LOR grouping in 32 bit data. Second, two simplified approaches for 3D time-of-flight (TOF) scatter estimation were developed to accelerate the computationally intensive calculation without compromising accuracy. The proposed list-mode reconstruction algorithm was compared to the manufacturer's point spread function + TOF (PSF+TOF) algorithm. Phantom, animal, and human studies demonstrated that MOLAR with pLOR gives slightly faster contrast recovery than the PSF+TOF algorithm that uses the average 32 bit LOR sinogram positioning. Moving phantom and a whole-body human study suggested that event-by-event motion correction reduces image blurring caused by respiratory motion. We conclude that list-mode reconstruction with pLOR positioning provides a platform to generate high quality images for the mCT, and to recover fine structures in whole-body PET scans through event-by-event motion correction.