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

Roadmap toward the 10 ps time-of-flight PET challenge

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

What scans we will read: imaging instrumentation trends in clinical oncology

Oncological diseases account for a significant portion of the burden on public healthcare systems with associated costs driven primarily by complex and long-lasting therapies. Through the visualization of patient-specific morphology and functional-molecular pathways, cancerous tissue can be detected and characterized non-invasively, so as to provide referring oncologists with essential information to support therapy management decisions. Following the onset of stand-alone anatomical and functional imaging, we witness a push towards integrating molecular image information through various methods, including anato-metabolic imaging (e.g., PET/CT), advanced MRI, optical or ultrasound imaging.This perspective paper highlights a number of key technological and methodological advances in imaging instrumentation related to anatomical, functional, molecular medicine and hybrid imaging, that is understood as the hardware-based combination of complementary anatomical and molecular imaging. These include novel detector technologies for ionizing radiation used in CT and nuclear medicine imaging, and novel system developments in MRI and optical as well as opto-acoustic imaging. We will also highlight new data processing methods for improved non-invasive tissue characterization. Following a general introduction to the role of imaging in oncology patient management we introduce imaging methods with well-defined clinical applications and potential for clinical translation. For each modality, we report first on the status quo and, then point to perceived technological and methodological advances in a subsequent status go section. Considering the breadth and dynamics of these developments, this perspective ends with a critical reflection on where the authors, with the majority of them being imaging experts with a background in physics and engineering, believe imaging methods will be in a few years from now.Overall, methodological and technological medical imaging advances are geared towards increased image contrast, the derivation of reproducible quantitative parameters, an increase in volume sensitivity and a reduction in overall examination time. To ensure full translation to the clinic, this progress in technologies and instrumentation is complemented by advances in relevant acquisition and image-processing protocols and improved data analysis. To this end, we should accept diagnostic images as "data", and - through the wider adoption of advanced analysis, including machine learning approaches and a "big data" concept - move to the next stage of non-invasive tumour phenotyping. The scans we will be reading in 10 years from now will likely be composed of highly diverse multi-dimensional data from multiple sources, which mandate the use of advanced and interactive visualization and analysis platforms powered by Artificial Intelligence (AI) for real-time data handling by cross-specialty clinical experts with a domain knowledge that will need to go beyond that of plain imaging.

Estimation of Crystal Timing Properties and Efficiencies for the Improvement of (Joint) Maximum-Likelihood Reconstructions in TOF-PET

With increasing improvements in the time of flight (TOF) resolution of positron emission tomography (PET) scanners, an accurate model of the TOF measurements is becoming increasingly important. This work considers two parameters of the TOF kernel; the relative positioning of the timing data-bins and the timing resolution along each line of response (LOR). Similar to an existing data-driven method, we assume that any shifts of data-bins along lines of response can be modelled as differences between crystal timing offsets. Inspired by this, timing resolutions of all LORs are modelled as the hypotenuse of timing resolutions of the crystal-pairs in coincidence. Furthermore, in order to mitigate the influence of potential inaccuracies of detector-pair sensitivities on crystal timing resolutions, relative LOR sensitivities are modelled as the product of efficiency factors for the two crystals in coincidence. We validate estimating maps of crystal timing offsets, timing resolutions and efficiencies from the emission data using noisy simulations of a brain phantom. Results are shown for phantom and patient data scanned on clinically available TOF-PET scanners. We find that the estimation of crystal timing resolutions is more sensitive to the data statistics than the estimation of crystal timing offsets. As a result, estimation of crystal timing properties could either be limited to high count emission data, or be obtained utilizing additional regularizations on the estimates. Using a more accurate model of the TOF acquisition, improvements are observed in standard activity reconstructions as well as joint reconstructions of activity and attenuation.

Evaluation of improved attenuation correction in whole-body PET/MR on patients with bone metastasis using various radiotracers

Purpose

This study evaluates the quantitative effect of improved MR-based attenuation correction (AC), including bone segmentation and the HUGE method for truncation correction in PET/MR whole-body hybrid imaging specifically of oncologic patients with bone metastasis and using various radiotracers.

Methods

Twenty-three patients that underwent altogether 28 whole-body PET/MR examinations with findings of bone metastasis were included in this study. Different radiotracers (¹⁸F-FDG, ⁶⁸Ga-PSMA, ⁶⁸Ga-DOTATOC, ¹²⁴I–MIBG) were injected according to appropriate clinical indications. Each of the 28 whole-body PET datasets was reconstructed three times using AC with (1) standard four-compartment μ-maps (background air, lung, muscle, and soft tissue), (2) five-compartment μ-maps (adding bone), and (3) six-compartment μ-maps (adding bone and HUGE truncation correction). The SUV_max of each detected bone lesion was measured in each reconstruction to evaluate the quantitative impact of improved MR-based AC. Relative difference images between four- and six-compartment μ-maps were calculated. MR-based HUGE truncation correction was compared with the PET-based MLAA truncation correction method in all patients.

Results

Overall, 69 bone lesions were detected and evaluated. The mean increase in relative difference over all 69 lesions in SUV_max was 5.4 ± 6.4% when comparing the improved six-compartment AC with the standard four-compartment AC. Maximal relative difference of 28.4% was measured in one lesion. Truncation correction with HUGE worked robust and resulted in realistic body contouring in all 28 exams and for all 4 different radiotracers. Truncation correction with MLAA revealed overestimations of arm tissue volume in all PET/MR exams with ¹⁸F-FDG radiotracer and failed in all other exams with radiotracers ⁶⁸Ga-PSMA, ⁶⁸Ga-DOTATOC, and ¹²⁴I- MIBG due to limitations in body contour detection.

Conclusion

Improved MR-based AC, including bone segmentation and HUGE truncation correction in whole-body PET/MR on patients with bone lesions and using various radiotracers, is important to ensure best possible diagnostic image quality and accurate PET quantification. The HUGE method for truncation correction based on MR worked robust and results in realistic body contouring, independent of the radiotracers used.

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.

Machine Learning in PET: from Photon Detection to Quantitative Image Reconstruction

Machine learning has found unique applications in nuclear medicine from photon detection to quantitative image reconstruction. While there have been impressive strides in detector development for time-of-flight positron emission tomography, most detectors still make use of simple signal processing methods to extract the time and position information from the detector signals. Now with the availability of fast waveform digitizers, machine learning techniques have been applied to estimate the position and arrival time of high-energy photons. In quantitative image reconstruction, machine learning has been used to estimate various corrections factors, including scattered events and attenuation images, as well as to reduce statistical noise in reconstructed images. Here machine learning either provides a faster alternative to an existing time-consuming computation, such as in the case of scatter estimation, or creates a data-driven approach to map an implicitly defined function, such as in the case of estimating the attenuation map for PET/MR scans. In this article, we will review the abovementioned applications of machine learning in nuclear medicine.

Metal artifact correction strategies in MRI-based attenuation correction in PET/MRI

In hybrid positron emission tomography (PET) and MRI systems, attenuation correction for PET image reconstruction is commonly based on processing of dedicated MR images. The image quality of the latter is strongly affected by metallic objects inside the body, such as e.g. dental implants, endoprostheses, or surgical clips which all lead to substantial artifacts that propagate into MRI-based attenuation images. In this work, we review publications about metal artifact correction strategies in MRI-based attenuation correction in PET/MRI. Moreover, we also give an overview about publications investigating the impact of MRI-based attenuation correction metal artifacts on the reconstructed PET image quality and quantification.

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.

Design, Implementation, and Evaluation of a Head and Neck MRI RF Array Integrated with a 511 keV Transmission Source for Attenuation Correction in PET/MR

The goal of this work is to further improve positron emission tomography (PET) attenuation correction and magnetic resonance (MR) sensitivity for head and neck applications of PET/MR. A dedicated 24-channel receive-only array, fully-integrated with a hydraulic system to move a transmission source helically around the patient and radiofrequency (RF) coil array, is designed, implemented, and evaluated. The device enables the calculation of attenuation coefficients from PET measurements at 511 keV including the RF coil and the particular patient. The RF coil design is PET-optimized by minimizing photon attenuation from coil components and housing. The functionality of the presented device is successfully demonstrated by calculating the attenuation map of a water bottle based on PET transmission measurements; results are in excellent agreement with reference values. It is shown that the device itself has marginal influence on the static magnetic field B₀ and the radiofrequency transmit field B₁ of the 3T PET/MR system. Furthermore, the developed RF array is shown to outperform a standard commercial 16-channel head and neck coil in terms of signal-to-noise ratio (SNR) and parallel imaging performance. In conclusion, the presented hardware enables accurate calculation of attenuation maps for PET/MR systems while improving the SNR of corresponding MR images in a single device without degrading the B₀ and B₁ homogeneity of the scanner.

Zero-Extra-Dose PET Delayed Imaging with Data-Driven Attenuation Correction Estimation

PURPOSE

Delayed positron emission tomography (PET) imaging may improve sensitivity and specificity in lesion detection. We proposed a PET data-driven method to estimate the attenuation map (AM) for the delayed scan without an additional x-ray computed tomography (CT).

PROCEDURES

An emission-attenuation-scatter joint estimation framework was developed. Several practical issues for clinical datasets were addressed. Particularly, the unknown scatter correction was incorporated in the joint estimation algorithm. The scaling problem was solved using prior information from the early CT scan. Fourteen patient datasets were added to evaluate the method. These patients went through two separate PET/CT scans. The delayed CT-based AM served as ground truth for the delayed scan. Standard uptake values (SUVmean and SUVmax) of lesion and normal tissue regions of interests (ROIs) in the early and delayed phase and the respective %DSUV (percentage change of SUVmean at two different time points) were analyzed, all with estimated and the true AM. Three radiologists participated in lesion detection tasks with images reconstructed with both AMs and rated scores for detectability.

RESULTS

The mean relative difference of SUVmean in lesion and normal liver tissue were 3.30 and 6.69 %. The average lesion-to-background contrast (detectability) with delayed PET images using CT AM was 60 % higher than that of the earlier PET image, and was 64 % higher when using the data-based AM. %DSUV for lesions and liver backgrounds with CT-based AM were - 0.058 ± 0.25 and - 0.33 ± 0.08 while with data-based AM were - 0.00 ± 0.26 and - 0.28 ± 0.08. Only slight significance difference was found between using CT-based AM and using the data-based AM reconstruction delay phase on %DSUV of lesion. The scores associated with the two AMs matched well consistently.

CONCLUSIONS

Our method may be used in delayed PET imaging, which allows no secondary CT radiation in delayed phase. The quantitative analysis for lesion detection purpose could be ensured.

Direct attenuation correction of brain PET images using only emission data via a deep convolutional encoder-decoder (Deep-DAC)

OBJECTIVE

To obtain attenuation-corrected PET images directly from non-attenuation-corrected images using a convolutional encoder-decoder network.

METHODS

Brain PET images from 129 patients were evaluated. The network was designed to map non-attenuation-corrected (NAC) images to pixel-wise continuously valued measured attenuation-corrected (MAC) PET images via an encoder-decoder architecture. Image quality was evaluated using various evaluation metrics. Image quantification was assessed for 19 radiomic features in 83 brain regions as delineated using the Hammersmith atlas (n30r83). Reliability of measurements was determined using pixel-wise relative errors (RE; %) for radiomic feature values in reference MAC PET images.

RESULTS

Peak signal-to-noise ratio (PSNR) and structural similarity index metric (SSIM) values were 39.2 ± 3.65 and 0.989 ± 0.006 for the external validation set, respectively. RE (%) of SUV_mean was - 0.10 ± 2.14 for all regions, and only 3 of 83 regions depicted significant differences. However, the mean RE (%) of this region was 0.02 (range, - 0.83 to 1.18). SUV_max had mean RE (%) of - 3.87 ± 2.84 for all brain regions, and 17 regions in the brain depicted significant differences with respect to MAC images with a mean RE of - 3.99 ± 2.11 (range, - 8.46 to 0.76). Homogeneity amongst Haralick-based radiomic features had the highest number (20) of regions with significant differences with a mean RE (%) of 7.22 ± 2.99.

CONCLUSIONS

Direct AC of PET images using deep convolutional encoder-decoder networks is a promising technique for brain PET images. The proposed deep learning method shows significant potential for emission-based AC in PET images with applications in PET/MRI and dedicated brain PET scanners.

KEY POINTS

• We demonstrate direct emission-based attenuation correction of PET images without using anatomical information. • We performed radiomics analysis of 83 brain regions to show robustness of direct attenuation correction of PET images. • Deep learning methods have significant promise for emission-based attenuation correction in PET images with potential applications in PET/MRI and dedicated brain PET scanners.

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 improved attenuation correction on 18F-FDG PET/MR hybrid imaging of the heart

PURPOSE

The aim of this study was to evaluate and quantify the effect of improved attenuation correction (AC) including bone segmentation and truncation correction on 18F-Fluordesoxyglucose cardiac positron emission tomography/magnetic resonance (PET/MR) imaging.

METHODS

PET data of 32 cardiac PET/MR datasets were reconstructed with three different AC-maps (1. Dixon-VIBE only, 2. HUGE truncation correction and bone segmentation, 3. MLAA). The Dixon-VIBE AC-maps served as reference of reconstructed PET data. 17-segment short-axis polar plots of the left ventricle were analyzed regarding the impact of each of the three AC methods on PET quantification in cardiac PET/MR imaging. Non-AC PET images were segmented to specify the amount of truncation in the Dixon-VIBE AC-map serving as a reference. All AC-maps were evaluated for artifacts.

RESULTS

Using HUGE + bone AC results in a homogeneous gain of ca. 6% and for MLAA 8% of PET signal distribution across the myocardium of the left ventricle over all patients compared to Dixon-VIBE AC only. Maximal relative differences up to 18% were observed in segment 17 (apex). The body volume truncation of -12.7 ± 7.1% compared to the segmented non-AC PET images using the Dixon-VIBE AC method was reduced to -1.9 ± 3.9% using HUGE and 7.8 ± 8.3% using MLAA. In each patient, a systematic overestimation in AC-map volume was observed when applying MLAA. Quantitative impact of artifacts showed regional differences up to 6% within single segments of the myocardium.

CONCLUSIONS

Improved AC including bone segmentation and truncation correction in cardiac PET/MR imaging is important to ensure best possible diagnostic quality and PET quantification. The results exhibited an overestimation of AC-map volume using MLAA, while HUGE resulted in a more realistic body contouring. Incorporation of bone segmentation into the Dixon-VIBE AC-map resulted in homogeneous gain in PET signal distribution across the myocardium. The majority of observed AC-map artifacts did not significantly affect the quantitative assessment of the myocardium.

Generation of PET Attenuation Map for Whole-Body Time-of-Flight 18F-FDG PET/MRI Using a Deep Neural Network Trained with Simultaneously Reconstructed Activity and Attenuation Maps

We propose a new deep learning-based approach to provide more accurate whole-body PET/MRI attenuation correction than is possible with the Dixon-based 4-segment method. We use activity and attenuation maps estimated using the maximum-likelihood reconstruction of activity and attenuation (MLAA) algorithm as inputs to a convolutional neural network (CNN) to learn a CT-derived attenuation map. Methods: The whole-body ¹⁸F-FDG PET/CT scan data of 100 cancer patients (38 men and 62 women; age, 57.3 ± 14.1 y) were retrospectively used for training and testing the CNN. A modified U-net was trained to predict a CT-derived μ-map (μ-CT) from the MLAA-generated activity distribution (λ-MLAA) and μ-map (μ-MLAA). We used 1.3 million patches derived from 60 patients' data for training the CNN, data of 20 others were used as a validation set to prevent overfitting, and the data of the other 20 were used as a test set for the CNN performance analysis. The attenuation maps generated using the proposed method (μ-CNN), μ-MLAA, and 4-segment method (μ-segment) were compared with the μ-CT, a ground truth. We also compared the voxelwise correlation between the activity images reconstructed using ordered-subset expectation maximization with the μ-maps, and the SUVs of primary and metastatic bone lesions obtained by drawing regions of interest on the activity images. Results: The CNN generates less noisy attenuation maps and achieves better bone identification than MLAA. The average Dice similarity coefficient for bone regions between μ-CNN and μ-CT was 0.77, which was significantly higher than that between μ-MLAA and μ-CT (0.36). Also, the CNN result showed the best pixel-by-pixel correlation with the CT-based results and remarkably reduced differences in activity maps in comparison to CT-based attenuation correction. Conclusion: The proposed deep neural network produced a more reliable attenuation map for 511-keV photons than the 4-segment method currently used in whole-body PET/MRI studies.

Feasibility study of a point-of-care positron emission tomography system with interactive imaging capability

PURPOSE

We investigated the feasibility of a novel positron emission tomography (PET) system that provides near real-time feedback to an operator who can interactively scan a patient to optimize image quality. The system should be compact and mobile to support point-of-care (POC) molecular imaging applications. In this study, we present the key technologies required and discuss the potential benefits of such new capability.

METHODS

The core of this novel PET technology includes trackable PET detectors and a fully three-dimensional, fast image reconstruction engine implemented on multiple graphics processing units (GPUs) to support dynamically changing geometry by calculating the system matrix on-the-fly using a tube-of-response approach. With near real-time image reconstruction capability, a POC-PET system may comprise a maneuverable front PET detector and a second detector panel which can be stationary or moved synchronously with the front detector such that both panels face the region-of-interest (ROI) with the detector trajectory contoured around a patient's body. We built a proof-of-concept prototype using two planar detectors each consisting of a photomultiplier tube (PMT) optically coupled to an array of 48 × 48 lutetium-yttrium oxyorthosilicate (LYSO) crystals (1.0 × 1.0 × 10.0 mm³ each). Only 38 × 38 crystals in each arrays can be clearly re-solved and used for coincidence detection. One detector was mounted to a robotic arm which can position it at arbitrary locations, and the other detector was mounted on a rotational stage. A cylindrical phantom (102 mm in diameter, 150 mm long) with nine spherical lesions (8:1 tumor-to-background activity concentration ratio) was imaged from 27 sampling angles. List-mode events were reconstructed to form images without or with time-of-flight (TOF) information. We conducted two Monte Carlo simulations using two POC-PET systems. The first one uses the same phantom and detector setup as our experiment, with the detector coincidence re-solving time (CRT) ranging from 100 to 700 ps full-width-at-half-maximum (FWHM). The second study simulates a body-size phantom (316 × 228 × 160 mm³ ) imaged by a larger POC-PET system that has 4 × 6 modules (32 × 32 LYSO crystals/module, four in axial and six in transaxial directions) in the front panel and 3 × 8 modules (16 × 16 LYSO crystals/module, three in axial and eight in transaxial directions) in the back panel. We also evaluated an interactive scanning strategy by progressively increasing the number of data sets used for image reconstruction. The updated images were analyzed based on the number of data sets and the detector CRT.

RESULTS

The proof-of-concept prototype re-solves most of the spherical lesions despite a limited number of coincidence events and incomplete sampling. TOF information reduces artifacts in the reconstructed images. Systems with better timing resolution exhibit improved image quality and reduced artifacts. We observed a reconstruction speed of 0.96 × 10⁶ events/s/iteration for 600 × 600 × 224 voxel rectilinear space using four GPUs. A POC-PET system with significantly higher sensitivity can interactively image a body-size object from four angles in less than 7 min.

CONCLUSIONS

We have developed GPU-based fast image reconstruction capability to support a PET system with arbitrary and dynamically changing geometry. Using TOF PET detectors, we demonstrated the feasibility of a PET system that can provide timely visual feedback to an operator who can scan a patient interactively to support POC imaging applications.

Evaluation of zero-echo-time attenuation correction for integrated PET/MR brain imaging-comparison to head atlas and 68Ge-transmission-based attenuation correction

BACKGROUND

MRI does not offer a direct method to obtain attenuation correction maps as its predecessors (stand-alone PET and PET/CT), and bone visualisation is particularly challenging. Recently, zero-echo-time (ZTE) was suggested for MR-based attenuation correction (AC). The aim of this work was to evaluate ZTE- and atlas-AC by comparison to 68Ge-transmission scan-based AC. Nine patients underwent brain PET/MR and stand-alone PET scanning using the dopamine transporter ligand 11C-PE2I. For each of them, two AC maps were obtained from the MR images: an atlas-based, obtained from T1-weighted LAVA-FLEX imaging with cortical bone inserted using a CT-based atlas, and an AC map generated from proton-density-weighted ZTE images. Stand-alone PET 68Ge-transmission AC map was used as gold standard. PET images were reconstructed using the three AC methods and standardised uptake value (SUV) values for the striatal, limbic and cortical regions, as well as the cerebellum (VOIs) were compared. SUV ratio (SUVR) values normalised for the cerebellum were also assessed. Bias, precision and agreement were calculated; statistical significance was evaluated using Wilcoxon matched-pairs signed-rank test.

RESULTS

Both ZTE- and atlas-AC showed a similar bias of 6-8% in SUV values across the regions. Correlation coefficients with 68Ge-AC were consistently high for ZTE-AC (r 0.99 for all regions), whereas they were lower for atlas-AC, varying from 0.99 in the striatum to 0.88 in the posterior cortical regions. SUVR showed an overall bias of 2.9 and 0.5% for atlas-AC and ZTE-AC, respectively. Correlations with 68Ge-AC were higher for ZTE-AC, varying from 0.99 in the striatum to 0.96 in the limbic regions, compared to atlas-AC (0.99 striatum to 0.77 posterior cortex).

CONCLUSIONS

Absolute SUV values showed less variability for ZTE-AC than for atlas-AC when compared to 68Ge-AC, but bias was similar for both methods. This bias is largely caused by higher linear attenuation coefficients in atlas- and ZTE-AC image compared to 68Ge-images. For SUVR, bias was lower when using ZTE-AC than for atlas-AC. ZTE-AC shows to be a more robust technique than atlas-AC in terms of both intra- and inter-patient variability.

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

Methods for joint activity reconstruction and attenuation reconstruction of time-of-flight (TOF) PET data provide an effective solution to attenuation correction when no (or incomplete or inaccurate) information on attenuation is available. One of the main barriers limiting use of these methods in clinical practice is their lack of validation in a relatively large patient database. In this contribution, we aim to validate reconstruction performed with maximum-likelihood activity reconstruction and attenuation registration (MLRR) in a whole-body patient dataset. Furthermore, a partial validation (because the scale problem of the algorithm is avoided for now) of reconstruction performed with maximum-likelihood activity and attenuation (MLAA) is also provided. We present a quantitative comparison between these 2 methods of joint reconstruction and the current clinical gold standard, maximum-likelihood expectation maximization (MLEM) with CT-based attenuation correction. Methods: The whole-body TOF PET emission data of each patient dataset were processed as a whole to reconstruct an activity volume covering all the acquired bed positions, helping reduce the problem of a scale per bed position in MLAA to a global scale for the entire activity volume. Three reconstruction algorithms were used: MLEM, MLRR, and MLAA. A maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results for various regions of interest were then analyzed. Results: The joint reconstructions of the whole-body patient dataset provided better quantification than the gold standard in cases of PET and CT misalignment caused by patient or organ motion. Our quantitative analysis showed a difference of -4.2% ± 2.3% between MLRR and MLEM and a difference of -7.5% ± 4.6% between MLAA and MLEM, averaged over all regions of interest. Conclusion: Joint reconstruction of activity and attenuation provides a useful means to estimate tracer distribution when CT-based-attenuation images are subject to misalignment or are not available. With an accurate estimate of the scatter contribution in the emission measurements, the joint reconstructions of TOF PET data are within clinically acceptable accuracy.

Simultaneous Emission-Transmission Tomography in an MRI Hardware Framework

Multi-modality imaging is essential for diagnosis and therapy in challenging cases. A Holy Grail of medical imaging is a hybrid imaging system combining computed tomography (CT), nuclear imaging, and magnetic resonance imaging (MRI) to deliver registered morphological, functional, and cellular/molecular information simultaneously and quantitatively for precision medicine. Recently, a unique imaging approach was demonstrated that combines nuclear imaging with polarized radiotracers and MRI-based spatial encoding. The detection scheme exploits the directional preference of γ-rays emitted from the polarized nuclei, and the result is a concentration image with resolution that can outperform standard nuclear imaging at a sensitivity significantly higher than that of MRI. However, the method does not calculate the attenuation image. Here we propose to obtain MRI-modulated γ-ray data for simultaneous image reconstruction of emission and transmission parameters, which could serve as a stepping stone toward simultaneous CT-SPECT-MRI. This method acquires synchronized datasets to provide insight into morphological features and molecular activities with accurate spatiotemporal registration. We present a complete overview of the system design and the formulation for tomographic reconstruction when the distribution of polarized radiotracers is either global or limited to a region of interest (ROI). Numerical results support the feasibility of our approach and suggest further research topics.

PET/MRI: a frontier in era of complementary hybrid imaging

With primitive approaches, the diagnosis and therapy were operated at the cellular, molecular, or even at the genetic level. As the diagnostic techniques are more concentrated towards molecular level, multi modal imaging becomes specifically essential. Multi-modal imaging has extensive applications in clinical as well as in pre-clinical studies. Positron Emission Tomography (PET) has flourished in the field of nuclear medicine, which has motivated it to fuse with Computed Tomography (CT) and Magnetic Resonance Imaging (MRI) for PET/CT and PET/MRI respectively. However, the challenges in PET/CT are due to the inability of simultaneous acquisition and reduced soft tissue contrast, which has led to the development of PET/MRI. Also, MRI offers the better soft tissue contrast over CT. Hence, fusion of PET and MRI results in combining structural information with functional image from PET. Yet, it has many technical challenges due to the interference between the modalities. Also, it must be resolved with various approaches for addressing the shortcomings of each system and improvise on the image quantification system. This review elaborates on the various challenges in the present PET/MRI system and the future directions of the hybrid modality. Also, the different data acquisition and analysis techniques of PET/MRI system are discussed with enhanced details on the software tools.

Reconstruction/segmentation of attenuation map in TOF-PET based on mixture models

Attenuation correction is known as a necessary step in positron emission tomography (PET) system to have accurate and quantitative activity images. Emission-based method is known as a promising approach for attenuation map estimation on TOF-PET scanners. The proposed method in this study imposes additional histogram-based information as a mixture model prior on the emission-based approach using maximum a posteriori (MAP) framework to improve its performance and make such a nearly segmented attenuation map. To eliminate misclassification of histogram modeling, a Median root prior is incorporated on the proposed approach to reduce the noise between neighbor voxels and encourage spatial smoothness in the reconstructed attenuation map. The joint-MAP optimization is carried out as an iterative approach wherein an alteration of the activity and attenuation updates is followed by a mixture decomposition of the attenuation map histogram. Also, the proposed method can segment attenuation map during the reconstruction. The evaluation of the proposed method on the numerical, simulation and real contexts indicate that the presented method has the potential to be used as a stand-alone method or even combined with other methods for attenuation correction on PET/MR systems.

Attenuation correction for brain PET imaging using deep neural network based on Dixon and ZTE MR images

Positron emission tomography (PET) is a functional imaging modality widely used in neuroscience studies. To obtain meaningful quantitative results from PET images, attenuation correction is necessary during image reconstruction. For PET/MR hybrid systems, PET attenuation is challenging as magnetic resonance (MR) images do not reflect attenuation coefficients directly. To address this issue, we present deep neural network methods to derive the continuous attenuation coefficients for brain PET imaging from MR images. With only Dixon MR images as the network input, the existing U-net structure was adopted and analysis using forty patient data sets shows it is superior to other Dixon-based methods. When both Dixon and zero echo time (ZTE) images are available, we have proposed a modified U-net structure, named GroupU-net, to efficiently make use of both Dixon and ZTE information through group convolution modules when the network goes deeper. Quantitative analysis based on fourteen real patient data sets demonstrates that both network approaches can perform better than the standard methods, and the proposed network structure can further reduce the PET quantification error compared to the U-net structure.

PET/MRI Hybrid Systems

Over the last decade, the combination of PET and MRI in one system has proven to be highly successful in basic preclinical research, as well as in clinical research. Nowadays, PET/MRI systems are well established in preclinical imaging and are progressing into clinical applications to provide further insights into specific diseases, therapeutic assessments, and biological pathways. Certain challenges in terms of hardware had to be resolved concurrently with the development of new techniques to be able to reach the full potential of both combined techniques. This review provides an overview of these challenges and describes the opportunities that simultaneous PET/MRI systems can exploit in comparison with stand-alone or other combined hybrid systems. New approaches were developed for simultaneous PET/MRI systems to correct for attenuation of 511 keV photons because MRI does not provide direct information on gamma photon attenuation properties. Furthermore, new algorithms to correct for motion were developed, because MRI can accurately detect motion with high temporal resolution. The additional information gained by the MRI can be employed to correct for partial volume effects as well. The development of new detector designs in combination with fast-decaying scintillator crystal materials enabled time-of-flight detection and incorporation in the reconstruction algorithms. Furthermore, this review lists the currently commercially available systems both for preclinical and clinical imaging and provides an overview of applications in both fields. In this regard, special emphasis has been placed on data analysis and the potential for both modalities to evolve with advanced image analysis tools, such as cluster analysis and machine learning.

Joint estimation of activity and attenuation for PET using pragmatic MR-based prior: application to clinical TOF PET/MR whole-body data for FDG and non-FDG tracers

Accurate and robust attenuation correction remains challenging in hybrid PET/MR particularly for torsos because it is difficult to segment bones, lungs and internal air in MR images. Additionally, MR suffers from susceptibility artifacts when a metallic implant is present. Recently, joint estimation (JE) of activity and attenuation based on PET data, also known as maximum likelihood reconstruction of activity and attenuation, has gained considerable interest because of (1) its promise to address the challenges in MR-based attenuation correction (MRAC), and (2) recent advances in time-of-flight (TOF) technology, which is known to be the key to the success of JE. In this paper, we implement a JE algorithm using an MR-based prior and evaluate the algorithm using whole-body PET/MR patient data, for both FDG and non-FDG tracers, acquired from GE SIGNA PET/MR scanners with TOF capability. The weight of the MR-based prior is spatially modulated, based on MR signal strength, to control the balance between MRAC and JE. Large prior weights are used in strong MR signal regions such as soft tissue and fat (i.e. MR tissue classification with a high degree of certainty) and small weights are used in low MR signal regions (i.e. MR tissue classification with a low degree of certainty). The MR-based prior is pragmatic in the sense that it is convex and does not require training or population statistics while exploiting synergies between MRAC and JE. We demonstrate the JE algorithm has the potential to improve the robustness and accuracy of MRAC by recovering the attenuation of metallic implants, internal air and some bones and by better delineating lung boundaries, not only for FDG but also for more specific non-FDG tracers such as 68Ga-DOTATOC and 18F-Fluoride.

Respiratory Motion Compensation for PET/CT with Motion Information Derived from Matched Attenuation-Corrected Gated PET Data

Respiratory motion degrades the detection and quantification capabilities of PET/CT imaging. Moreover, mismatch between a fast helical CT image and a time-averaged PET image due to respiratory motion results in additional attenuation correction artifacts and inaccurate localization. Current motion compensation approaches typically have 3 limitations: the mismatch among respiration-gated PET images and the CT attenuation correction (CTAC) map can introduce artifacts in the gated PET reconstructions that can subsequently affect the accuracy of the motion estimation; sinogram-based correction approaches do not correct for intragate motion due to intracycle and intercycle breathing variations; and the mismatch between the PET motion compensation reference gate and the CT image can cause an additional CT-mismatch artifact. In this study, we established a motion correction framework to address these limitations. Methods: In the proposed framework, the combined emission-transmission reconstruction algorithm was used for phase-matched gated PET reconstructions to facilitate the motion model building. An event-by-event nonrigid respiratory motion compensation method with correlations between internal organ motion and external respiratory signals was used to correct both intracycle and intercycle breathing variations. The PET reference gate was automatically determined by a newly proposed CT-matching algorithm. We applied the new framework to 13 human datasets with 3 different radiotracers and 323 lesions and compared its performance with CTAC and non-attenuation correction (NAC) approaches. Validation using 4-dimensional CT was performed for one lung cancer dataset. Results: For the 10 ¹⁸F-FDG studies, the proposed method outperformed (P < 0.006) both the CTAC and the NAC methods in terms of region-of-interest-based SUV_mean, SUV_max, and SUV ratio improvements over no motion correction (SUV_mean: 19.9% vs. 14.0% vs. 13.2%; SUV_max: 15.5% vs. 10.8% vs. 10.6%; SUV ratio: 24.1% vs. 17.6% vs. 16.2%, for the proposed, CTAC, and NAC methods, respectively). The proposed method increased SUV ratios over no motion correction for 94.4% of lesions, compared with 84.8% and 86.4% using the CTAC and NAC methods, respectively. For the 2 ¹⁸F-fluoropropyl-(+)-dihydrotetrabenazine studies, the proposed method reduced the CT-mismatch artifacts in the lower lung where the CTAC approach failed and maintained the quantification accuracy of bone marrow where the NAC approach failed. For the ¹⁸F-FMISO study, the proposed method outperformed both the CTAC and the NAC methods in terms of motion estimation accuracy at 2 lung lesion locations. Conclusion: The proposed PET/CT respiratory event-by-event motion-correction framework with motion information derived from matched attenuation-corrected PET data provides image quality superior to that of the CTAC and NAC methods for multiple tracers.

Combined PET/MRI: Global Warming-Summary Report of the 6th International Workshop on PET/MRI, March 27-29, 2017, Tübingen, Germany

The 6th annual meeting to address key issues in positron emission tomography (PET)/magnetic resonance imaging (MRI) was held again in Tübingen, Germany, from March 27 to 29, 2017. Over three days of invited plenary lectures, round table discussions and dialogue board deliberations, participants critically assessed the current state of PET/MRI, both clinically and as a research tool, and attempted to chart future directions. The meeting addressed the use of PET/MRI and workflows in oncology, neurosciences, infection, inflammation and chronic pain syndromes, as well as deeper discussions about how best to characterise the tumour microenvironment, optimise the complementary information available from PET and MRI, and how advanced data mining and bioinformatics, as well as information from liquid biomarkers (circulating tumour cells and nucleic acids) and pathology, can be integrated to give a more complete characterisation of disease phenotype. Some issues that have dominated previous meetings, such as the accuracy of MR-based attenuation correction (AC) of the PET scan, were finally put to rest as having been adequately addressed for the majority of clinical situations. Likewise, the ability to standardise PET systems for use in multicentre trials was confirmed, thus removing a perceived barrier to larger clinical imaging trials. The meeting openly questioned whether PET/MRI should, in all cases, be used as a whole-body imaging modality or whether in many circumstances it would best be employed to give an in-depth study of previously identified disease in a single organ or region. The meeting concluded that there is still much work to be done in the integration of data from different fields and in developing a common language for all stakeholders involved. In addition, the participants advocated joint training and education for individuals who engage in routine PET/MRI. It was agreed that PET/MRI can enhance our understanding of normal and disrupted biology, and we are in a position to describe the in vivo nature of disease processes, metabolism, evolution of cancer and the monitoring of response to pharmacological interventions and therapies. As such, PET/MRI is a key to advancing medicine and patient care.

Principles of Simultaneous PET/MR Imaging

Combined PET/MR imaging scanners capable of acquiring simultaneously the complementary information provided by the 2 imaging modalities are now available for human use. After addressing the hardware challenges for integrating the 2 imaging modalities, most of the efforts in the field have focused on developing MR-based attenuation correction methods for neurologic and whole-body applications, implementing approaches for improving one modality by using the data provided by the other and exploring research and clinical applications that could benefit from the synergistic use of the multimodal data.

MR-guided joint reconstruction of activity and attenuation in brain PET-MR

With the advent of time-of-flight (TOF) PET scanners, joint maximum-likelihood reconstruction of activity and attenuation (MLAA) maps has recently regained attention for the estimation of PET attenuation maps from emission data. However, the estimated attenuation and activity maps are scaled by unknown scaling factors. We recently demonstrated that in hybrid PET-MR, the scaling issue of this algorithm can be effectively addressed by imposing MR spatial constraints on the estimation of attenuation maps using a penalized MLAA (P-MLAA⁺) algorithm. With the advent of simultaneous PET-MR systems, MRI-guided PET image reconstruction has also gained attention for improving the quantitative accuracy of PET images, usually degraded by noise and partial volume effects. The aim of this study is therefore to increase the benefits of MRI information for improving the quantitative accuracy of PET images by exploiting MRI-based anatomical penalty functions to guide the reconstruction of both activity and attenuation maps during their joint estimation. We employed an anato-functional joint entropy penalty function for the reconstruction of activity and an anatomical quadratic penalty function for the reconstruction of attenuation. The resulting algorithm was referred to as P-MLAA⁺⁺ since it exploits both activity and attenuation penalty functions. The performance of the P-MLAA algorithms were compared with MLAA and the widely used activity reconstruction algorithms such as maximum likelihood expectation maximization (MLEM) and penalized MLEM (P-MLEM) both corrected for attenuation using a conventional MRI segmentation-based attenuation correction (MRAC) method. The studied methods were evaluated using simulations and clinical studies taking the PET image reconstructed using reference CT-based attenuation maps as a reference. The simulation results showed that the proposed method can notably improve the visual quality of the PET images by reducing noise while preserving structural boundaries and at the same time improving the quantitative accuracy of the PET images. Our clinical reconstruction results showed that the MLEM-MRAC, P-MLEM-MRAC, MLAA, P-MLAA⁺ and P-MLAA⁺⁺ algorithms result in, on average, quantification errors of -13.5 ± 3.1%, -13.4 ± 3.1%, -2.0 ± 6.5%, -3.0 ± 3.5% and -4.2 ± 3.6%, respectively, in different regions of the brain. In conclusion, whilst the P-MLAA⁺ algorithm showed the best overall quantification performance, the proposed P-MLAA⁺⁺ algorithm provided simultaneous partial volume and attenuation corrections with only a minor compromise of PET quantification.

Numerical algorithms for scatter-to-attenuation reconstruction in PET: empirical comparison of convergence, acceleration, and the effect of subsets

The use of scattered coincidences for attenuation correction of positron emission tomography (PET) data has recently been proposed. For practical applications, convergence speeds require further improvement, yet there exists a trade-off between convergence speed and the risk of non-convergence. In this respect, a maximum-likelihood gradient-ascent (MLGA) algorithm and a two-branch back-projection (2BP), which was previously proposed, were evaluated.

METHODS

MLGA was combined with the Armijo step size rule; and accelerated using conjugate gradients, Nesterov's momentum method, and data subsets of different sizes. In 2BP, we varied the subset size, an important determinant of convergence speed and computational burden. We used three sets of simulation data to evaluate the impact of a spatial scale factor.

RESULTS AND DISCUSSION

The Armijo step size allowed 10-fold increased step sizes compared to native MLGA. Conjugate gradients and Nesterov momentum lead to slightly faster, yet non-uniform convergence; improvements were mostly confined to later iterations, possibly due to the non-linearity of the problem. MLGA with data subsets achieved faster, uniform, and predictable convergence, with a speed-up factor equivalent to the number of subsets and no increase in computational burden. By contrast, 2BP computational burden increased linearly with the number of subsets due to repeated evaluation of the objective function, and convergence was limited to the case of many (and therefore small) subsets, which resulted in high computational burden.

CONCLUSION

Possibilities of improving 2BP appear limited. While general-purpose acceleration methods appear insufficient for MLGA, results suggest that data subsets are a promising way of improving MLGA performance.

Plane-dependent ML scatter scaling: 3D extension of the 2D simulated single scatter (SSS) estimate

Scatter correction is typically done using a simulation of the single scatter, which is then scaled to account for multiple scatters and other possible model mismatches. This scaling factor is determined by fitting the simulated scatter sinogram to the measured sinogram, using only counts measured along LORs that do not intersect the patient body, i.e. 'scatter-tails'. Extending previous work, we propose to scale the scatter with a plane dependent factor, which is determined as an additional unknown in the maximum likelihood (ML) reconstructions, using counts in the entire sinogram rather than only the 'scatter-tails'. The ML-scaled scatter estimates are validated using a Monte-Carlo simulation of a NEMA-like phantom, a phantom scan with typical contrast ratios of a ⁶⁸Ga-PSMA scan, and 23 whole-body ¹⁸F-FDG patient scans. On average, we observe a 12.2% change in the total amount of tracer activity of the MLEM reconstructions of our whole-body patient database when the proposed ML scatter scales are used. Furthermore, reconstructions using the ML-scaled scatter estimates are found to eliminate the typical 'halo' artifacts that are often observed in the vicinity of high focal uptake regions.

Joint estimation of activity image and attenuation sinogram using time-of-flight positron emission tomography data consistency condition filtering

Attenuation correction is essential for quantitative reliability of positron emission tomography (PET) imaging. In time-of-flight (TOF) PET, attenuation sinogram can be determined up to a global constant from noiseless emission data due to the TOF PET data consistency condition. This provides the theoretical basis for jointly estimating both activity image and attenuation sinogram/image directly from TOF PET emission data. Multiple joint estimation methods, such as maximum likelihood activity and attenuation (MLAA) and maximum likelihood attenuation correction factor (MLACF), have already been shown that can produce improved reconstruction results in TOF cases. However, due to the nonconcavity of the joint log-likelihood function and Poisson noise presented in PET data, the iterative method still requires proper initialization and well-designed regularization to prevent convergence to local maxima. To address this issue, we propose a joint estimation of activity image and attenuation sinogram using the TOF PET data consistency condition as an attenuation sinogram filter, and then evaluate the performance of the proposed method using computer simulations.

MLAA-based attenuation correction of flexible hardware components in hybrid PET/MR imaging

BACKGROUND

Accurate PET quantification demands attenuation correction (AC) for both patient and hardware attenuation of the 511 keV annihilation photons. In hybrid PET/MR imaging, AC for stationary hardware components such as patient table and MR head coil is straightforward, employing CT-derived attenuation templates. AC for flexible hardware components such as MR-safe headphones and MR radiofrequency (RF) surface coils is more challenging. Registration-based approaches, aligning CT-based attenuation templates with the current patient position, have been proposed but are not used in clinical routine. Ignoring headphone or RF coil attenuation has been shown to result in regional activity underestimation values of up to 18%. We propose to employ the maximum-likelihood reconstruction of attenuation and activity (MLAA) algorithm to estimate the attenuation of flexible hardware components. Starting with an initial attenuation map not including flexible hardware components, the attenuation update of MLAA is applied outside the body outline only, allowing to estimate hardware attenuation without modifying the patient attenuation map. Appropriate prior expectations on the attenuation coefficients are incorporated into MLAA. The proposed method is investigated for non-TOF PET phantom and ¹⁸F-FDG patient data acquired with a clinical PET/MR device, using headphones or RF surface coils as flexible hardware components.

RESULTS

Although MLAA cannot recover the exact physical shape of the hardware attenuation maps, the overall attenuation of the hardware components is accurately estimated. Therefore, the proposed algorithm significantly improves PET quantification. Using the phantom data, local activity underestimation when neglecting hardware attenuation was reduced from up to 25% to less than 3% under- or overestimation as compared to reference scans without hardware present or to CT-derived AC. For the patient data, we found an average activity underestimation of 7.9% evaluated in the full brain and of 6.1% for the abdominal region comparing the uncorrected case with MLAA.

CONCLUSIONS

MLAA is able to provide accurate estimations of the attenuation of flexible hardware components and can therefore be used to significantly improve PET quantification. The proposed approach can be readily incorporated into clinical workflow.

Attenuation correction in 4D-PET using a single-phase attenuation map and rigidity-adaptive deformable registration

PURPOSE

Four-dimensional positron emission tomography (4D-PET) imaging is a potential solution to the respiratory motion effect in the thoracic region. Computed tomography (CT)-based attenuation correction (AC) is an essential step toward quantitative imaging for PET. However, due to the temporal difference between 4D-PET and a single attenuation map from CT, typically available in routine clinical scanning, motion artifacts are observed in the attenuation-corrected PET images, leading to errors in tumor shape and uptake. We introduced a practical method to align single-phase CT with all other 4D-PET phases for AC.

METHODS

A penalized non-rigid Demons registration between individual 4D-PET frames without AC provides the motion vectors to be used for warping single-phase attenuation map. The non-rigid Demons registration was used to derive deformation vector fields (DVFs) between PET matched with the CT phase and other 4D-PET images. While attenuated PET images provide useful data for organ borders such as those of the lung and the liver, tumors cannot be distinguished from the background due to loss of contrast. To preserve the tumor shape in different phases, an ROI-covering tumor was excluded from nonrigid transformation. Instead the mean DVF of the central region of the tumor was assigned to all voxels in the ROI. This process mimics a rigid transformation of the tumor along with a nonrigid transformation of other organs. A 4D-XCAT phantom with spherical lung tumors, with diameters ranging from 10 to 40 mm, was used to evaluate the algorithm. The performance of the proposed hybrid method for attenuation map estimation was compared to (a) the Demons nonrigid registration only and (b) a single attenuation map based on quantitative parameters in individual PET frames.

RESULTS

Motion-related artifacts were significantly reduced in the attenuation-corrected 4D-PET images. When a single attenuation map was used for all individual PET frames, the normalized root-mean-square error (NRMSE) values in tumor region were 49.3% (STD: 8.3%), 50.5% (STD: 9.3%), 51.8% (STD: 10.8%) and 51.5% (STD: 12.1%) for 10-mm, 20-mm, 30-mm, and 40-mm tumors, respectively. These errors were reduced to 11.9% (STD: 2.9%), 13.6% (STD: 3.9%), 13.8% (STD: 4.8%), and 16.7% (STD: 9.3%) by our proposed method for deforming the attenuation map. The relative errors in total lesion glycolysis (TLG) values were -0.25% (STD: 2.87%) and 3.19% (STD: 2.35%) for 30-mm and 40-mm tumors, respectively, in proposed method. The corresponding values for Demons method were 25.22% (STD: 14.79%) and 18.42% (STD: 7.06%). Our proposed hybrid method outperforms the Demons method especially for larger tumors. For tumors smaller than 20 mm, nonrigid transformation could also provide quantitative results.

CONCLUSION

Although non-AC 4D-PET frames include insignificant anatomical information, they are still useful to estimate the DVFs to align the attenuation map for accurate AC. The proposed hybrid method can recover the AC-related artifacts and provide quantitative AC-PET images.

PET/MRI in the Presence of Metal Implants: Completion of the Attenuation Map from PET Emission Data

We present a novel technique for accurate whole-body attenuation correction in the presence of metallic endoprosthesis, on integrated non-time-of-flight (non-TOF) PET/MRI scanners. The proposed implant PET-based attenuation map completion (IPAC) method performs a joint reconstruction of radioactivity and attenuation from the emission data to determine the position, shape, and linear attenuation coefficient (LAC) of metallic implants. Methods: The initial estimate of the attenuation map was obtained using the MR Dixon method currently available on the Siemens Biograph mMR scanner. The attenuation coefficients in the area of the MR image subjected to metal susceptibility artifacts are then reconstructed from the PET emission data using the IPAC algorithm. The method was tested on 11 subjects presenting 13 different metallic implants, who underwent CT and PET/MR scans. Relative mean LACs and Dice similarity coefficients were calculated to determine the accuracy of the reconstructed attenuation values and the shape of the metal implant, respectively. The reconstructed PET images were compared with those obtained using the reference CT-based approach and the Dixon-based method. Absolute relative change (aRC) images were generated in each case, and voxel-based analyses were performed. Results: The error in implant LAC estimation, using the proposed IPAC algorithm, was 15.7% ± 7.8%, which was significantly smaller than the Dixon- (100%) and CT- (39%) derived values. A mean Dice similarity coefficient of 73% ± 9% was obtained when comparing the IPAC- with the CT-derived implant shape. The voxel-based analysis of the reconstructed PET images revealed quantification errors (aRC) of 13.2% ± 22.1% for the IPAC- with respect to CT-corrected images. The Dixon-based method performed substantially worse, with a mean aRC of 23.1% ± 38.4%. Conclusion: We have presented a non-TOF emission-based approach for estimating the attenuation map in the presence of metallic implants, to be used for whole-body attenuation correction in integrated PET/MR scanners. The Graphics Processing Unit implementation of the algorithm will be included in the open-source reconstruction toolbox Occiput.io.

Clinical evaluation of TOF versus non-TOF on PET artifacts in simultaneous PET/MR: a dual centre experience

PURPOSE

Our objective was to determine clinically the value of time-of-flight (TOF) information in reducing PET artifacts and improving PET image quality and accuracy in simultaneous TOF PET/MR scanning.

METHODS

A total 65 patients who underwent a comparative scan in a simultaneous TOF PET/MR scanner were included. TOF and non-TOF PET images were reconstructed, clinically examined, compared and scored. PET imaging artifacts were categorized as large or small implant-related artifacts, as dental implant-related artifacts, and as implant-unrelated artifacts. Differences in image quality, especially those related to (implant) artifacts, were assessed using a scale ranging from 0 (no artifact) to 4 (severe artifact).

RESULTS

A total of 87 image artifacts were found and evaluated. Four patients had large and eight patients small implant-related artifacts, 27 patients had dental implants/fillings, and 48 patients had implant-unrelated artifacts. The average score was 1.14 ± 0.82 for non-TOF PET images and 0.53 ± 0.66 for TOF images (p < 0.01) indicating that artifacts were less noticeable when TOF information was included.

CONCLUSION

Our study indicates that PET image artifacts are significantly mitigated with integration of TOF information in simultaneous PET/MR. The impact is predominantly seen in patients with significant artifacts due to metal implants.

The Effect of Susceptibility Artifacts Related to Metallic Implants on Adjacent-Lesion Assessment in Simultaneous TOF PET/MR

Metalic implants may affect attenuation correction (AC) in PET/MR imaging. The purpose of this study was to evaluate the effect of susceptibility artifacts related to metallic implants on adjacent metabolically active lesions in clinical simultaneous PET/MR scanning for both time-of-flight (TOF) and non-TOF reconstructed PET images. Methods: We included 27 patients without implants but with confirmed ¹⁸F-FDG-avid lesions adjacent to common implant locations. In all patients, a clinically indicated whole-body ¹⁸F-FDG PET/MR scan was acquired. Baseline non-TOF and TOF PET images were reconstructed. Reconstruction was repeated after the introduction of artificial signal voids in the AC map to simulate metallic implants in standard anatomic areas. All reconstructed images were qualitatively and quantitatively assessed and compared with the baseline images. Results: In total, 51 lesions were assessed. In 40 and 50 of these cases (non-TOF and TOF, respectively), the detectability of the lesions did not change; in 9 and 1 cases, the detectability changed; and in 2 non-TOF cases, the lesions were no longer visible after the introduction of metallic artifacts. The inclusion of TOF information significantly reduced artifacts due to simulated implants in the femoral head, sternum, and spine (P = 0.01, 0.01, and 0.03, respectively). It also improved image quality in these locations (P = 0.02, 0.01, and 0.01, respectively). The mean percentage error was -3.5% for TOF and -4.8% for non-TOF reconstructions, meaning that the inclusion of TOF information reduced the percentage error in SUV_max by 28.5% (P < 0.01). Conclusion: Qualitatively, there was a significant reduction of artifacts in the femoral head, sternum, and spine. There was also a significant qualitative improvement in image quality in these locations. Furthermore, our study indicated that simulated susceptibility artifacts related to metallic implants have a significant effect on small, moderately ¹⁸F-FDG-avid lesions near the implant site that possibly may go unnoticed without TOF information. On larger, highly ¹⁸F-FDG-avid lesions, the metallic implants had only a limited effect. The largest significant quantitative difference was found in artifacts of the sternum. There was only a weak inverse correlation between lesions affected by artifacts and distance from the implant.

An Expectation Maximization Method for Joint Estimation of Emission Activity Distribution and Photon Attenuation Map in PET

A maximum likelihood expectation maximization (MLEM) method is proposed for joint estimation of emission activity distribution and photon attenuation map from positron emission tomography (PET) emission data alone. The method is appealing since: (i) it guarantees monotonic likelihood increase to a local extremum, (ii) does not require arbitrary parameters, and (iii) guarantees the positivity of the estimated distributions. Moreover, we propose a discrete Poisson data acquisition model and numerical algorithm for: (i) efficient graphics processing unit (GPU) based formulation, and (ii) a closed form exact solution for the MLEM update equations, which is essential for accurate and robust estimation. Numerical experiments indicate that in the presence of noise, joint EMAA estimation converges to the true emission activity distribution with root mean square errors of 4% and 0.5% respectively in estimation of lung- and myocardial emission activity distributions for a computational XCAT thorax phantom.

A multi-centre evaluation of eleven clinically feasible brain PET/MRI attenuation correction techniques using a large cohort of patients

AIM

To accurately quantify the radioactivity concentration measured by PET, emission data need to be corrected for photon attenuation; however, the MRI signal cannot easily be converted into attenuation values, making attenuation correction (AC) in PET/MRI challenging. In order to further improve the current vendor-implemented MR-AC methods for absolute quantification, a number of prototype methods have been proposed in the literature. These can be categorized into three types: template/atlas-based, segmentation-based, and reconstruction-based. These proposed methods in general demonstrated improvements compared to vendor-implemented AC, and many studies report deviations in PET uptake after AC of only a few percent from a gold standard CT-AC. Using a unified quantitative evaluation with identical metrics, subject cohort, and common CT-based reference, the aims of this study were to evaluate a selection of novel methods proposed in the literature, and identify the ones suitable for clinical use.

METHODS

In total, 11 AC methods were evaluated: two vendor-implemented (MR-AC_DIXON and MR-AC_UTE), five based on template/atlas information (MR-AC_SEGBONE (Koesters et al., 2016), MR-AC_ONTARIO (Anazodo et al., 2014), MR-AC_BOSTON (Izquierdo-Garcia et al., 2014), MR-AC_UCL (Burgos et al., 2014), and MR-AC_MAXPROB (Merida et al., 2015)), one based on simultaneous reconstruction of attenuation and emission (MR-AC_MLAA (Benoit et al., 2015)), and three based on image-segmentation (MR-AC_MUNICH (Cabello et al., 2015), MR-AC_CAR-RiDR (Juttukonda et al., 2015), and MR-AC_RESOLUTE (Ladefoged et al., 2015)). We selected 359 subjects who were scanned using one of the following radiotracers: [¹⁸F]FDG (210), [¹¹C]PiB (51), and [¹⁸F]florbetapir (98). The comparison to AC with a gold standard CT was performed both globally and regionally, with a special focus on robustness and outlier analysis.

RESULTS

The average performance in PET tracer uptake was within ±5% of CT for all of the proposed methods, with the average±SD global percentage bias in PET FDG uptake for each method being: MR-AC_DIXON (-11.3±3.5)%, MR-AC_UTE (-5.7±2.0)%, MR-AC_ONTARIO (-4.3±3.6)%, MR-AC_MUNICH (3.7±2.1)%, MR-AC_MLAA (-1.9±2.6)%, MR-AC_SEGBONE (-1.7±3.6)%, MR-AC_UCL (0.8±1.2)%, MR-AC_CAR-RiDR (-0.4±1.9)%, MR-AC_MAXPROB (-0.4±1.6)%, MR-AC_BOSTON (-0.3±1.8)%, and MR-AC_RESOLUTE (0.3±1.7)%, ordered by average bias. The overall best performing methods (MR-AC_BOSTON, MR-AC_MAXPROB, MR-AC_RESOLUTE and MR-AC_UCL, ordered alphabetically) showed regional average errors within ±3% of PET with CT-AC in all regions of the brain with FDG, and the same four methods, as well as MR-AC_CAR-RiDR, showed that for 95% of the patients, 95% of brain voxels had an uptake that deviated by less than 15% from the reference. Comparable performance was obtained with PiB and florbetapir.

CONCLUSIONS

All of the proposed novel methods have an average global performance within likely acceptable limits (±5% of CT-based reference), and the main difference among the methods was found in the robustness, outlier analysis, and clinical feasibility. Overall, the best performing methods were MR-ACBOSTON, MR-ACMAXPROB, MR-ACRESOLUTE and MR-ACUCL, ordered alphabetically. These methods all minimized the number of outliers, standard deviation, and average global and local error. The methods MR-ACMUNICH and MR-ACCAR-RiDR were both within acceptable quantitative limits, so these methods should be considered if processing time is a factor. The method MR-ACSEGBONE also demonstrates promising results, and performs well within the likely acceptable quantitative limits. For clinical routine scans where processing time can be a key factor, this vendor-provided solution currently outperforms most methods. With the performance of the methods presented here, it may be concluded that the challenge of improving the accuracy of MR-AC in adult brains with normal anatomy has been solved to a quantitatively acceptable degree, which is smaller than the quantification reproducibility in PET imaging.

Vision 20/20: Magnetic resonance imaging-guided attenuation correction in PET/MRI: Challenges, solutions, and opportunities

Attenuation correction is an essential component of the long chain of data correction techniques required to achieve the full potential of quantitative positron emission tomography (PET) imaging. The development of combined PET/magnetic resonance imaging (MRI) systems mandated the widespread interest in developing novel strategies for deriving accurate attenuation maps with the aim to improve the quantitative accuracy of these emerging hybrid imaging systems. The attenuation map in PET/MRI should ideally be derived from anatomical MR images; however, MRI intensities reflect proton density and relaxation time properties of biological tissues rather than their electron density and photon attenuation properties. Therefore, in contrast to PET/computed tomography, there is a lack of standardized global mapping between the intensities of MRI signal and linear attenuation coefficients at 511 keV. Moreover, in standard MRI sequences, bones and lung tissues do not produce measurable signals owing to their low proton density and short transverse relaxation times. MR images are also inevitably subject to artifacts that degrade their quality, thus compromising their applicability for the task of attenuation correction in PET/MRI. MRI-guided attenuation correction strategies can be classified in three broad categories: (i) segmentation-based approaches, (ii) atlas-registration and machine learning methods, and (iii) emission/transmission-based approaches. This paper summarizes past and current state-of-the-art developments and latest advances in PET/MRI attenuation correction. The advantages and drawbacks of each approach for addressing the challenges of MR-based attenuation correction are comprehensively described. The opportunities brought by both MRI and PET imaging modalities for deriving accurate attenuation maps and improving PET quantification will be elaborated. Future prospects and potential clinical applications of these techniques and their integration in commercial systems will also be discussed.