Joint Reconstruction of Activity and Attenuation in Time-of-Flight PET: A Quantitative Analysis

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.

Attenuation correction and truncation completion for breast PET/MR imaging using deep learning

Objective. Simultaneous PET/MR scanners combine the high sensitivity of MR imaging with the functional imaging of PET. However, attenuation correction of breast PET/MR imaging is technically challenging. The purpose of this study is to establish a robust attenuation correction algorithm for breast PET/MR images that relies on deep learning (DL) to recreate the missing portions of the patient's anatomy (truncation completion), as well as to provide bone information for attenuation correction from only the PET data.Approach. Data acquired from 23 female subjects with invasive breast cancer scanned with18F-fluorodeoxyglucose PET/CT and PET/MR localized to the breast region were used for this study. Three DL models, U-Net with mean absolute error loss (DLMAE) model, U-Net with mean squared error loss (DLMSE) model, and U-Net with perceptual loss (DLPerceptual) model, were trained to predict synthetic CT images (sCT) for PET attenuation correction (AC) given non-attenuation corrected (NAC) PETPET/MRimages as inputs. The DL and Dixon-based sCT reconstructed PET images were compared against those reconstructed from CT images by calculating the percent error of the standardized uptake value (SUV) and conducting Wilcoxon signed rank statistical tests.Main results. sCT images from the DLMAEmodel, the DLMSEmodel, and the DLPerceptualmodel were similar in mean absolute error (MAE), peak-signal-to-noise ratio, and normalized cross-correlation. No significant difference in SUV was found between the PET images reconstructed using the DLMSEand DLPerceptualsCTs compared to the reference CT for AC in all tissue regions. All DL methods performed better than the Dixon-based method according to SUV analysis.Significance. A 3D U-Net with MSE or perceptual loss model can be implemented into a reconstruction workflow, and the derived sCT images allow successful truncation completion and attenuation correction for breast PET/MR images.

A review of PET attenuation correction methods for PET-MR

Despite being thirteen years since the installation of the first PET-MR system, the scanners constitute a very small proportion of the total hybrid PET systems installed. This is in stark contrast to the rapid expansion of the PET-CT scanner, which quickly established its importance in patient diagnosis within a similar timeframe. One of the main hurdles is the development of an accurate, reproducible and easy-to-use method for attenuation correction. Quantitative discrepancies in PET images between the manufacturer-provided MR methods and the more established CT- or transmission-based attenuation correction methods have led the scientific community in a continuous effort to develop a robust and accurate alternative. These can be divided into four broad categories: (i) MR-based, (ii) emission-based, (iii) atlas-based and the (iv) machine learning-based attenuation correction, which is rapidly gaining momentum. The first is based on segmenting the MR images in various tissues and allocating a predefined attenuation coefficient for each tissue. Emission-based attenuation correction methods aim in utilising the PET emission data by simultaneously reconstructing the radioactivity distribution and the attenuation image. Atlas-based attenuation correction methods aim to predict a CT or transmission image given an MR image of a new patient, by using databases containing CT or transmission images from the general population. Finally, in machine learning methods, a model that could predict the required image given the acquired MR or non-attenuation-corrected PET image is developed by exploiting the underlying features of the images. Deep learning methods are the dominant approach in this category. Compared to the more traditional machine learning, which uses structured data for building a model, deep learning makes direct use of the acquired images to identify underlying features. This up-to-date review goes through the literature of attenuation correction approaches in PET-MR after categorising them. The various approaches in each category are described and discussed. After exploring each category separately, a general overview is given of the current status and potential future approaches along with a comparison of the four outlined categories.

Synthetic PET via Domain Translation of 3-D MRI

Historically, patient datasets have been used to develop and validate various reconstruction algorithms for PET/MRI and PET/CT. To enable such algorithm development, without the need for acquiring hundreds of patient exams, in this article we demonstrate a deep learning technique to generate synthetic but realistic whole-body PET sinograms from abundantly available whole-body MRI. Specifically, we use a dataset of 56 18F-FDG-PET/MRI exams to train a 3-D residual UNet to predict physiologic PET uptake from whole-body T1-weighted MRI. In training, we implemented a balanced loss function to generate realistic uptake across a large dynamic range and computed losses along tomographic lines of response to mimic the PET acquisition. The predicted PET images are forward projected to produce synthetic PET (sPET) time-of-flight (ToF) sinograms that can be used with vendor-provided PET reconstruction algorithms, including using CT-based attenuation correction (CTAC) and MR-based attenuation correction (MRAC). The resulting synthetic data recapitulates physiologic 18F-FDG uptake, e.g., high uptake localized to the brain and bladder, as well as uptake in liver, kidneys, heart, and muscle. To simulate abnormalities with high uptake, we also insert synthetic lesions. We demonstrate that this sPET data can be used interchangeably with real PET data for the PET quantification task of comparing CTAC and MRAC methods, achieving ≤ 7.6% error in mean-SUV compared to using real data. These results together show that the proposed sPET data pipeline can be reasonably used for development, evaluation, and validation of PET/MRI reconstruction methods.

A quantitative clinical evaluation of simultaneous reconstruction of attenuation and activity in time-of-flight PET

BACKGROUND

The maximum likelihood activity and attenuation (MLAA) reconstruction algorithm has been proposed to jointly estimate tracer activity and attenuation at the same time, and proven to be a promising solution to the CT attenuation correction (CT-AC) artifacts in PET images. This study aimed to perform a quantitative evaluation and clinical validation of the MLAA method.

METHODS

A uniform cylinder phantom filled with ¹⁸F-FDG solution was scanned to optimize the reconstruction parameters for the implemented MLAA algorithm. 67 patients who underwent whole-body ¹⁸F-FDG PET/CT scan were retrospectively recruited. PET images were reconstructed using MLAA and clinical standard OSEM algorithm with CT-AC (CT-OSEM). The mean and maximum standardized uptake values (SUVmean and SUVmax) in regions of interest (ROIs) of organs, high uptake lesions and areas affected by metal implants and respiration motion artifacts were quantitatively analyzed.

RESULTS

In quantitative analysis, SUVs in patient's organ ROIs between two methods showed R² ranging from 0.91 to 0.98 and k ranging from 0.90 to 1.06, and the average SUVmax and SUVmean differences between two methods were within 10% range, except for the lung ROI, which was 10.5% and 16.73% respectively. The average SUVmax and SUVmean differences of a total of 117 high uptake lesions were 7.25% and 7.10% respectively. 20 patients were identified to have apparent respiration motion artifacts in the liver in CT-OSEM images, and the SUVs differences between two methods measured at dome of the liver were significantly larger than measured at middle part of the liver. 10 regions with obvious metal artifacts were identified in CT-OSEM images and the average SUVmean and SUVmax differences in metal implants affected regions were reported to be 52.90% and 56.20% respectively.

CONCLUSIONS

PET images reconstructed using MLAA are clinically acceptable in terms of image quality as well as quantification and it is a useful tool in clinical practice, especially when CT-AC may cause respiration motion and metal artifacts. Moreover, this study also provides technical reference and data support for the future iteration and development of PET reconstruction technology of SUV accurate quantification.

Deep learning-based attenuation map generation with simultaneously reconstructed PET activity and attenuation and low-dose application

Objective. In PET/CT imaging, CT is used for positron emission tomography (PET) attenuation correction (AC). CT artifacts or misalignment between PET and CT can cause AC artifacts and quantification errors in PET. Simultaneous reconstruction (MLAA) of PET activity (λ-MLAA) and attenuation (μ-MLAA) maps was proposed to solve those issues using the time-of-flight PET raw data only. However,λ-MLAA still suffers from quantification error as compared to reconstruction using the gold-standard CT-based attenuation map (μ-CT). Recently, a deep learning (DL)-based framework was proposed to improve MLAA by predictingμ-DL fromλ-MLAA andμ-MLAA using an image domain loss function (IM-loss). However, IM-loss does not directly measure the AC errors according to the PET attenuation physics. Our preliminary studies showed that an additional physics-based loss function can lead to more accurate PET AC. The main objective of this study is to optimize the attenuation map generation framework for clinical full-dose¹⁸F-FDG studies. We also investigate the effectiveness of the optimized network on predicting attenuation maps for synthetic low-dose oncological PET studies.Approach. We optimized the proposed DL framework by applying different preprocessing steps and hyperparameter optimization, including patch size, weights of the loss terms and number of angles in the projection-domain loss term. The optimization was performed based on 100 skull-to-toe¹⁸F-FDG PET/CT scans with minimal misalignment. The optimized framework was further evaluated on 85 clinical full-dose neck-to-thigh¹⁸F-FDG cancer datasets as well as synthetic low-dose studies with only 10% of the full-dose raw data.Main results. Clinical evaluation of tumor quantification as well as physics-based figure-of-merit metric evaluation validated the promising performance of our proposed method. For both full-dose and low-dose studies, the proposed framework achieved <1% error in tumor standardized uptake value measures.Significance. It is of great clinical interest to achieve CT-less PET reconstruction, especially for low-dose PET studies.

Deep learning-based attenuation correction for whole-body PET - a multi-tracer study with 18F-FDG, 68 Ga-DOTATATE, and 18F-Fluciclovine

A novel deep learning (DL)-based attenuation correction (AC) framework was applied to clinical whole-body oncology studies using 18F-FDG, 68 Ga-DOTATATE, and 18F-Fluciclovine. The framework used activity (λ-MLAA) and attenuation (µ-MLAA) maps estimated by the maximum likelihood reconstruction of activity and attenuation (MLAA) algorithm as inputs to a modified U-net neural network with a novel imaging physics-based loss function to learn a CT-derived attenuation map (µ-CT).

METHODS

Clinical whole-body PET/CT datasets of 18F-FDG (N = 113), 68 Ga-DOTATATE (N = 76), and 18F-Fluciclovine (N = 90) were used to train and test tracer-specific neural networks. For each tracer, forty subjects were used to train the neural network to predict attenuation maps (µ-DL). µ-DL and µ-MLAA were compared to the gold-standard µ-CT. PET images reconstructed using the OSEM algorithm with µ-DL (OSEMDL) and µ-MLAA (OSEMMLAA) were compared to the CT-based reconstruction (OSEMCT). Tumor regions of interest were segmented by two radiologists and tumor SUV and volume measures were reported, as well as evaluation using conventional image analysis metrics.

RESULTS

µ-DL yielded high resolution and fine detail recovery of the attenuation map, which was superior in quality as compared to µ-MLAA in all metrics for all tracers. Using OSEMCT as the gold-standard, OSEMDL provided more accurate tumor quantification than OSEMMLAA for all three tracers, e.g., error in SUVmax for OSEMMLAA vs. OSEMDL: - 3.6 ± 4.4% vs. - 1.7 ± 4.5% for 18F-FDG (N = 152), - 4.3 ± 5.1% vs. 0.4 ± 2.8% for 68 Ga-DOTATATE (N = 70), and - 7.3 ± 2.9% vs. - 2.8 ± 2.3% for 18F-Fluciclovine (N = 44). OSEMDL also yielded more accurate tumor volume measures than OSEMMLAA, i.e., - 8.4 ± 14.5% (OSEMMLAA) vs. - 3.0 ± 15.0% for 18F-FDG, - 14.1 ± 19.7% vs. 1.8 ± 11.6% for 68 Ga-DOTATATE, and - 15.9 ± 9.1% vs. - 6.4 ± 6.4% for 18F-Fluciclovine.

CONCLUSIONS

The proposed framework provides accurate and robust attenuation correction for whole-body 18F-FDG, 68 Ga-DOTATATE and 18F-Fluciclovine in tumor SUV measures as well as tumor volume estimation. The proposed method provides clinically equivalent quality as compared to CT in attenuation correction for the three tracers.

Data-driven gated (DDG) CT: An automated respiratory gating method to enable DDG PET/CT

BACKGROUND

The accuracy of PET quantification and localization can be compromised if a misregistered CT is used for attenuation correction (AC) in PET/CT. As data-driven gating (DDG) continues to grow in clinical use, these issues are becoming more relevant with respect to solutions for gated CT.

PURPOSE

In this work, a new automated data-driven gated (DDG) CT method was developed to provide average CT and DDG CT for AC of PET and DDG PET, respectively.

METHODS

An automatic DDG CT was developed to provide the end-expiratory (EE) and end-inspiratory (EI) phases of images from low-dose cine CT images, with all phases being averaged to generate an average CT. The respiratory phases of EE and EI were determined according to lung region Hounsfield unit (HU) values and body outline contours. The average CT was used for AC of baseline PET and DDG CT at EE phase was used for AC of DDG PET at the quiescent or EE phase. The EI and EE phases obtained with DDG CT were used for assessing the magnitude of respiratory motion. The proposed DDG CT was compared to two commercial CT gating methods: 1) 4D CT (external device based) and 2) D4D CT (DDG based) in 38 patient data sets with respect to respiratory phase image selection, lung HU, lung volume, and image artifacts. In a separate set of twenty consecutive PET/CT studies containing a mix of ¹⁸ F-FDG, ⁶⁸ Ga-Dotatate, and ⁶⁴ Cu-Dotatate scans, the proposed DDG CT was compared with D4D CT for impacts on registration and quantification in DDG PET/CT.

RESULTS

In the EE phase, the images selected by DDG CT and 4D CT were identical 62.5±21.6% of the time, while DDG CT and D4D CT were 6.5±9.7%, and 4D CT and D4D CT were 8.6±12.2%. These differences in EE phase image selection were significant (p<0.0001). In the EI phase, the images selected by DDG CT and 4D CT were identical 68.2±18.9% of the time, DDG CT and D4D CT were 63.9±18.8%, and 4D CT and D4D CT were 61.2±19.8%. These differences were not significant. The mean lung HU and volumes were not statistically different (p > 0.1) among the three methods. In some studies, DDG CT was better than D4D or 4D CT in appropriate selection of the EE and EI phases, and D4D CT was found to reverse the EE and EI phases or not select the correct images by visual inspection. A statistically significant improvement of DDG CT over D4D CT for AC of DDG PET was also demonstrated with PET quantification analysis. When irregular breath cycles were present in the cine CT, DDG CT could be used to replace average CT for improved AC of baseline PET.

CONCLUSION

A new automatic DDG CT was developed to tackle the issues of misregistration and tumor motion in PET/CT imaging. DDG CT was significantly more consistent than D4D CT in selecting the EE phase images as the clinical standard of 4D CT. When compared to both commercial gated CT methods of 4D CT and D4D CT, DDG CT appeared to be more robust in the lower lung and upper diaphragm regions where misregistration and tumor motion often occur. DDG CT offered improved AC for DDG PET relative to D4D CT. In cases with irregular respiratory motion, DDG CT improved AC over average CT for baseline PET. The new DDG CT provides the benefits of 4D CT without the need for external device gating. This article is protected by copyright. All rights reserved.

MRI-guided attenuation correction in torso PET/MRI: Assessment of segmentation-, atlas-, and deep learning-based approaches in the presence of outliers

Purpose

We compare the performance of three commonly used MRI‐guided attenuation correction approaches in torso PET/MRI, namely segmentation‐, atlas‐, and deep learning‐based algorithms.

Methods

Twenty‐five co‐registered torso ¹⁸F‐FDG PET/CT and PET/MR images were enrolled. PET attenuation maps were generated from in‐phase Dixon MRI using a three‐tissue class segmentation‐based approach (soft‐tissue, lung, and background air), voxel‐wise weighting atlas‐based approach, and a residual convolutional neural network. The bias in standardized uptake value (SUV) was calculated for each approach considering CT‐based attenuation corrected PET images as reference. In addition to the overall performance assessment of these approaches, the primary focus of this work was on recognizing the origins of potential outliers, notably body truncation, metal‐artifacts, abnormal anatomy, and small malignant lesions in the lungs.

Results

The deep learning approach outperformed both atlas‐ and segmentation‐based methods resulting in less than 4% SUV bias across 25 patients compared to the segmentation‐based method with up to 20% SUV bias in bony structures and the atlas‐based method with 9% bias in the lung. The deep learning‐based method exhibited superior performance. Yet, in case of sever truncation and metallic‐artifacts in the input MRI, this approach was outperformed by the atlas‐based method, exhibiting suboptimal performance in the affected regions. Conversely, for abnormal anatomies, such as a patient presenting with one lung or small malignant lesion in the lung, the deep learning algorithm exhibited promising performance compared to other methods.

Conclusion

The deep learning‐based method provides promising outcome for synthetic CT generation from MRI. However, metal‐artifact and body truncation should be specifically addressed.

Time of Flight in Perspective: Instrumental and Computational Aspects of Time Resolution in Positron Emission Tomography

The first time-of-flight positron emission tomography (TOF-PET) scanners were developed as early as in the 1980s. However, the poor light output and low detection efficiency of TOF-capable detectors available at the time limited any gain in image quality achieved with these TOF-PET scanners over the traditional non-TOF PET scanners. The discovery of LSO and other Lu-based scintillators revived interest in TOF-PET and led to the development of a second generation of scanners with high sensitivity and spatial resolution in the mid-2000s. The introduction of the silicon photomultiplier (SiPM) has recently yielded a third generation of TOF-PET systems with unprecedented imaging performance. Parallel to these instrumentation developments, much progress has been made in the development of image reconstruction algorithms that better utilize the additional information provided by TOF. Overall, the benefits range from a reduction in image variance (SNR increase), through allowing joint estimation of activity and attenuation, to better reconstructing data from limited angle systems. In this work, we review these developments, focusing on three broad areas: 1) timing theory and factors affecting the time resolution of a TOF-PET system; 2) utilization of TOF information for improved image reconstruction; and 3) quantification of the benefits of TOF compared to non-TOF PET. Finally, we offer a brief outlook on the TOF-PET developments anticipated in the short and longer term. Throughout this work, we aim to maintain a clinically driven perspective, treating TOF as one of multiple (and sometimes competitive) factors that can aid in the optimization of PET imaging performance.

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.

Deep learning-based T1-enhanced selection of linear attenuation coefficients (DL-TESLA) for PET/MR attenuation correction in dementia neuroimaging

PURPOSE

The accuracy of existing PET/MR attenuation correction (AC) has been limited by a lack of correlation between MR signal and tissue electron density. Based on our finding that longitudinal relaxation rate, or R1 , is associated with CT Hounsfield unit in bone and soft tissues in the brain, we propose a deep learning T1 -enhanced selection of linear attenuation coefficients (DL-TESLA) method to incorporate quantitative R1 for PET/MR AC and evaluate its accuracy and longitudinal test-retest repeatability in brain PET/MR imaging.

METHODS

DL-TESLA uses a 3D residual UNet (ResUNet) for pseudo-CT (pCT) estimation. With a total of 174 participants, we compared PET AC accuracy of DL-TESLA to 3 other methods adopting similar 3D ResUNet structures but using UTE R 2 ∗ , or Dixon, or T1 -MPRAGE as input. With images from 23 additional participants repeatedly scanned, the test-retest differences and within-subject coefficient of variation of standardized uptake value ratios (SUVR) were compared between PET images reconstructed using either DL-TESLA or CT for AC.

RESULTS

DL-TESLA had (1) significantly lower mean absolute error in pCT, (2) the highest Dice coefficients in both bone and air, (3) significantly lower PET relative absolute error in whole brain and various brain regions, (4) the highest percentage of voxels with a PET relative error within both ±3% and ±5%, (5) similar to CT test-retest differences in SUVRs from the cerebrum and mean cortical (MC) region, and (6) similar to CT within-subject coefficient of variation in cerebrum and MC.

CONCLUSION

DL-TESLA demonstrates excellent PET/MR AC accuracy and test-retest repeatability.

Quantitative PET in the 2020s: a roadmap

Positron emission tomography (PET) plays an increasingly important role in research and clinical applications, catalysed by remarkable technical advances and a growing appreciation of the need for reliable, sensitive biomarkers of human function in health and disease. Over the last 30 years, a large amount of the physics and engineering effort in PET has been motivated by the dominant clinical application during that period, oncology. This has led to important developments such as PET/CT, whole-body PET, 3D PET, accelerated statistical image reconstruction, and time-of-flight PET. Despite impressive improvements in image quality as a result of these advances, the emphasis on static, semi-quantitative 'hot spot' imaging for oncologic applications has meant that the capability of PET to quantify biologically relevant parameters based on tracer kinetics has not been fully exploited. More recent advances, such as PET/MR and total-body PET, have opened up the ability to address a vast range of new research questions, from which a future expansion of applications and radiotracers appears highly likely. Many of these new applications and tracers will, at least initially, require quantitative analyses that more fully exploit the exquisite sensitivity of PET and the tracer principle on which it is based. It is also expected that they will require more sophisticated quantitative analysis methods than those that are currently available. At the same time, artificial intelligence is revolutionizing data analysis and impacting the relationship between the statistical quality of the acquired data and the information we can extract from the data. In this roadmap, leaders of the key sub-disciplines of the field identify the challenges and opportunities to be addressed over the next ten years that will enable PET to realise its full quantitative potential, initially in research laboratories and, ultimately, in clinical practice.

A Role of PET/MR Imaging in Dementia?

Since many years, magnetic resonance imaging (MRI) and positron emission tomography (PET) have a prominent role in neurodegenerative disorders and dementia, not only in a research setting but also in a clinical setting. For several decades, information from both modalities is combined ranging from individual visual assessments to fully integrating all images. Several tools are available to coregister images from MRI and PET and to covisualize these images. When studying neurodegenerative disorders with PET it is important to perform a partial volume correction and this can be done using the structural information obtained by MRI. With the advent of PET/MR, the question arises in how far this hybrid imaging modality is an added value compared to combining PET and MRI data from two separate modalities. One issue in PET/MR is still not yet completely settled, that is, the attenuation correction. This is of less importance for visual assessments but it can become an issue when combining data from PET/CT and PET/MR scanners in multicenter studies or when using cut-off values to classify patients. Simultaneous imaging has clearly some advantages: for the patient it is beneficial to have only one scan session instead of two but also in cases in which PET data are related to functional of physiological data acquired with MRI (such as functional MRI or arterial spin labeling). However, the most important benefit is currently the more integrated use of PET and MRI. This is also possible with separate measurements but requires more streamlining of the whole process. In that case coregistration of images is mandatory. It needs to be determined in which cases simultaneous PET/MRI leads to new insights or improved diagnosis compared to multimodal imaging using dedicated scanners.

PET-enabled dual-energy CT: image reconstruction and a proof-of-concept computer simulation study

Standard dual-energy computed tomography (CT) uses two different x-ray energies to obtain energy-dependent tissue attenuation information to allow quantitative material decomposition. The combined use of dual-energy CT and positron emission tomography (PET) may provide a more comprehensive characterization of disease states in cancer and other diseases. However, the integration of dual-energy CT with PET is not trivial, either requiring costly hardware upgrades or increasing radiation exposure. This paper proposes a different dual-energy CT imaging method that is enabled by PET. Instead of using a second x-ray CT scan with a different energy, this method exploits time-of-flight PET image reconstruction via the maximum likelihood attenuation and activity (MLAA) algorithm to obtain a 511 keV gamma-ray attenuation image from PET emission data. The high-energy gamma-ray attenuation image is then combined with the low-energy x-ray CT of PET/CT to provide a pair of dual-energy CT images. A major challenge with the standard MLAA reconstruction is the high noise present in the reconstructed 511 keV attenuation map, which would not compromise the PET activity reconstruction too much but may significantly affect the performance of the gamma-ray attenuation image for material decomposition. To overcome the problem, we further propose a kernel MLAA algorithm to exploit the prior information from the available x-ray CT image. We conducted a computer simulation to test the concept and algorithm for the task of material decomposition. The simulation results demonstrate that this PET-enabled dual-energy CT method is promising for quantitative material decomposition. The proposed method can be readily implemented on time-of-flight PET/CT scanners to enable simultaneous PET and dual-energy CT imaging.

Simultaneous emission and attenuation reconstruction in time-of-flight PET using a reference object

Background

Simultaneous reconstruction of emission and attenuation images in time-of-flight (TOF) positron emission tomography (PET) does not provide a unique solution. In this study, we propose to solve this limitation by including additional information given by a reference object with known attenuation placed outside the patient. Different configurations of the reference object were studied including geometry, material composition, and activity, and an optimal configuration was defined. In addition, this configuration was tested for different timing resolutions and noise levels.

Results

The proposed strategy was tested in 2D simulations obtained by forward projection of available PET/CT data and noise was included using Monte Carlo techniques. Obtained results suggest that the optimal configuration corresponds to a water cylinder inserted in the patient table and filled with activity. In that case, mean differences between reconstructed and true images were below 10%. However, better results can be obtained by increasing the activity of the reference object.

Conclusion

This study shows promising results that might allow to obtain an accurate attenuation map from pure TOF-PET data without prior knowledge obtained from CT, MRI, or transmission scans.

PET/MRI attenuation estimation in the lung: A review of past, present, and potential techniques

Positron emission tomography/magnetic resonance imaging (PET/MRI) potentially offers several advantages over positron emission tomography/computed tomography (PET/CT), for example, no CT radiation dose and soft tissue images from MR acquired at the same time as the PET. However, obtaining accurate linear attenuation correction (LAC) factors for the lung remains difficult in PET/MRI. LACs depend on electron density and in the lung, these vary significantly both within an individual and from person to person. Current commercial practice is to use a single‐valued population‐based lung LAC, and better estimation is needed to improve quantification. Given the under‐appreciation of lung attenuation estimation as an issue, the inaccuracy of PET quantification due to the use of single‐valued lung LACs, the unique challenges of lung estimation, and the emerging status of PET/MRI scanners in lung disease, a review is timely. This paper highlights past and present methods, categorizing them into segmentation, atlas/mapping, and emission‐based schemes. Potential strategies for future developments are also presented.

Algorithms for joint activity-attenuation estimation from positron emission tomography scatter

Background

Attenuation correction in positron emission tomography remains challenging in the absence of measured transmission data. Scattered emission data may contribute missing information, but quantitative scatter-to-attenuation (S2A) reconstruction needs to input the reconstructed activity image. Here, we study S2A reconstruction as a building block for joint estimation of activity and attenuation.

Methods

We study two S2A reconstruction algorithms, maximum-likelihood expectation maximization (MLEM) with one-step-late attenuation (MLEM-OSL) and a maximum-likelihood gradient ascent (MLGA). We study theoretical properties of these algorithms with a focus on convergence and convergence speed and compare convergence speeds and the impact of object size in simulations using different spatial scale factors. Then, we propose joint estimation of activity and attenuation from scattered and nonscattered (true) emission data, combining MLEM-OSL or MLGA with scatter-MLEM as well as trues-MLEM and the maximum-likelihood transmission (MLTR) algorithm.

Results

Shortcomings of MLEM-OSL inhibit convergence to the true solution with high attenuation; these shortcomings are related to the linearization of a nonlinear measurement equation and can be linked to a new numerical criterion allowing geometrical interpretations in terms of low and high attenuation. Comparisons using simulated data confirm that while MLGA converges largely independent of the attenuation scale, MLEM-OSL converges if low-attenuation data dominate, but not with high attenuation. Convergence of MLEM-OSL can be improved by isolating data satisfying the aforementioned low-attenuation criterion. In joint estimation of activity and attenuation, scattered data helps avoid local minima that nonscattered data alone cannot. Combining MLEM-OSL with trues-MLEM may be sufficient for low-attenuation objects, while MLGA, scatter-MLEM, and MLTR may additionally be needed with higher attenuation.

Conclusions

The performance of S2A algorithms depends on spatial scales. MLGA provides lower computational complexity and convergence in more diverse setups than MLEM-OSL. Finally, scattered data may provide additional information to joint estimation of activity and attenuation through S2A reconstruction.

A Quantitative Evaluation of Joint Activity and Attenuation Reconstruction in TOF PET/MR Brain Imaging

Time-of-flight (TOF) PET data provide an effective means for attenuation correction (AC) when no (or incomplete or inaccurate) attenuation information is available. Since MR scanners provide little information on photon attenuation of different tissue types, AC in hybrid PET/MR scanners has always been challenging. In this contribution, we aim at validating the activity reconstructions of the maximum-likelihood ordered-subsets activity and attenuation (OSAA) reconstruction algorithm on a patient brain data set. We present a quantitative comparison of joint reconstructions with the current clinical gold standard-ordered-subsets expectation maximization-using CT-based AC in PET/CT, as well as the current state of the art in PET/MR, that is, zero time echo (ZTE)-based AC. Methods: The TOF PET emission data were initially used in a preprocessing stage to estimate crystal maps of efficiencies, timing offsets, and timing resolutions. Applying these additional corrections during reconstructions, OSAA, ZTE-based, and the vendor-provided atlas-based AC techniques were analyzed and compared with CT-based AC. In our initial study, we used the CT-based estimate of the expected scatter and later used the ZTE-based and OSAA attenuation estimates to compute the expected scatter contribution of the data during reconstructions. In all reconstructions, a maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results were analyzed in the 86 segmented regions of interest of the Hammers atlas. Results: Our quantitative analysis showed that, in practice, a tracer activity difference of +0.5% (±2.1%) and +0.1% (±2.3%) could be expected for the state-of-the-art ZTE-based and OSAA AC methods, respectively, in PET/MR compared with the clinical gold standard in PET/CT. Conclusion: Joint activity and attenuation estimation methods can provide an effective solution to the challenging AC problem for brain studies in hybrid TOF PET/MR scanners. With an accurate TOF-based (timing offsets and timing resolutions) calibration, and similar to the results of the state-of-the-art method in PET/MR, regional errors of joint TOF PET reconstructions are within a few percentage points.

Impact of Tissue Classification in MRI-Guided Attenuation Correction on Whole-Body Patlak PET/MRI

PURPOSE

The aim of this work is to investigate the impact of tissue classification in magnetic resonance imaging (MRI)-guided positron emission tomography (PET) attenuation correction (AC) for whole-body (WB) Patlak net uptake rate constant (K_i) imaging in PET/MRI studies.

PROCEDURES

WB dynamic PET/CT data were acquired for 14 patients. The CT images were utilized to generate attenuation maps (μ-map_CTAC) of continuous attenuation coefficient values (A_coeff). The μ-map_CTAC were then segmented into four tissue classes (μ-map_4-classes), namely background (air), lung, fat, and soft tissue, where a predefined A_coeff was assigned to each class. To assess the impact of bone for AC, the bones in the μ-map_CTAC were then assigned a predefined soft tissue A_coeff (0.1 cm^-1) to produce an AC μ-map without bones (μ-map_no-bones). Thereafter, both WB static SUV and dynamic PET images were reconstructed using μ-map_CTAC, μ-map_4-classes, and μ-map_no-bones (PET_CTAC, PET_4-classes, and PET_no-bones), respectively. WB indirect and direct parametric K_i images were generated using Patlak graphical analysis. Malignant lesions were delineated on PET images with an automatic segmentation method that uses an active contour model (MASAC). Then, the quantitative metrics of the metabolically active tumor volume (MATV), target-to-background (TBR), contrast-to-noise ratio (CNR), peak region-of-interest (ROI_peak), maximum region-of-interest (ROI_max), mean region-of-interest (ROI_mean), and metabolic volume product (MVP) were analyzed. The Wilcoxon test was conducted to assess the difference between PET_4-classes and PET_no-bones against PET_CTAC for all images. The same test was also adopted to compare the differences between SUV, indirect K_i, and direct K_i images for each evaluated AC method.

RESULTS

No significant differences in MATV, TBR, and CNR were observed between PET_4-classes and PET_CTAC for either SUV or K_i images. PET_4-classes significantly overestimated ROI_peak, ROI_max, ROI_mean, as well as MVP scores compared with PET_CTAC in both SUV and K_i images. SUV images exhibited the highest median relative errors for PET_4-classes with respect to PET_CTAC (RE_4-classes): 6.91 %, 6.55 %, 5.90 %, and 6.56 % for ROI_peak, ROI_max, ROI_mean, and MVP, respectively. On the contrary, K_i images showed slightly reduced RE_4-classes (indirect 5.52 %, 5.95 %, 4.43 %, and 5.70 %, direct 6.61 %, 6.33 %, 5.53 %, and 4.96 %) for ROI_peak, ROI_max, ROI_mean, and MVP, respectively. A higher TBR was observed on indirect and direct K_i images relative to SUV, while direct K_i images demonstrated the highest CNR.

CONCLUSIONS

Four-tissue class AC may impact SUV and K_i parameter estimation but only to a limited extent, thereby suggesting that WB Patlak K_i imaging for dynamic WB PET/MRI studies is feasible. Patlak K_i imaging can enhance TBR, thereby facilitating lesion segmentation and quantification. However, patient-specific A_coeff for each tissue class should be used when possible to address the high inter-patient variability of A_coeff distributions.

Magnetic Resonance Fingerprinting: Implications and Opportunities for PET/MR

Magnetic Resonance Imaging (MRI) can be used to assess anatomical structure, and its sensitivity to a variety of tissue properties enables superb contrast between tissues as well as the ability to characterize these tissues. However, despite vast potential for quantitative and functional evaluation, MRI is typically used qualitatively, in which the underlying tissue properties are not measured, and thus the brightness of each pixel is not quantitatively meaningful. Positron Emission Tomography (PET) is an inherently quantitative imaging modality that interrogates functional activity within a tissue, probed by a molecule of interest coupled with an appropriate tracer. These modalities can complement one another to provide clinical information regarding both structure and function, but there are still technical and practical hurdles in the way of the integrated use of both modalities. Recent advances in MRI have moved the field in an increasingly quantitative direction, which is complementary to PET, and could also potentially help solve some of the challenges in PET/MR. Magnetic Resonance Fingerprinting (MRF) is a recently described MRI-based technique which can efficiently and simultaneously quantitatively map several tissue properties in a single exam. Here, the basic principles behind the quantitative approach of MRF are laid out, and the potential implications for combined PET/MR are discussed.