Simultaneous reconstruction of activity and attenuation for PET/MR

Feasibility of PET-enabled dual-energy CT imaging: First physical phantom and initial patient study results

PURPOSE

Dual-energy (DE) CT enables material decomposition by using two different x-ray energies and may be combined with PET for improved multimodality imaging. However, this increases radiation dose and may require a hardware upgrade due to the added second x-ray CT scan. The recently proposed PET-enabled DECT method allows dual-energy imaging using a conventional PET/CT scanner without the need to change scanner hardware or increase radiation exposure. Here we demonstrate the first-time physical phantom and patient data evaluation of this method.

METHODS

The PET-enabled DECT method reconstructs a gamma-ray CT (gCT) image at 511 keV from the time-of-flight PET data with the maximum-likelihood attenuation and activity (MLAA) approach and then combines this image with the low-energy x-ray CT images to form a dual-energy image pair for material decomposition. To improve the image quality of gCT, a kernel MLAA method was developed using the x-ray CT as a priori information. Here we developed a general open-source implementation for gCT reconstruction and used this implementation for the first real data validation using both physical phantom study and human-subject study. Results from PET-enabled DECT were compared using x-ray DECT as the reference. Further, we applied the PET-enabled DECT method in another patient study to evaluate bone lesions.

RESULTS

Compared to the standard MLAA, results from the kernel MLAA showed significantly improved image quality. PET-enabled DECT with the kernel MLAA was able to generate fractional images that were comparable to the x-ray DECT, with high correlation coefficients for both the phantom study and human subject study (R > 0.99). The application study also indicates that PET-enabled DECT has potential to characterize bone lesions.

CONCLUSION

Results from this study have demonstrated the feasibility of this PET-enabled method for CT imaging and material decomposition. PET-enabled DECT shows promise to provide comparable results to x-ray DECT.

openSSS: an open-source implementation of scatter estimation for 3D TOF-PET

BACKGROUND

Scatter correction is essential for quantitative and accurate time-of-flight (TOF) PET imaging. It is implemented by an accurate scatter estimation algorithm, to calculate the statistical distribution of scattered photons among the measured coincidences. However, to our knowledge, scatter estimation algorithms that account for TOF and that are compatible with custom geometries are not available in open-source reconstruction libraries, such as CASToR and STIR. To this end, we have developed an open-source implementation of the TOF-aware single-scatter-simulation (SSS) algorithm: openSSS.

RESULTS

openSSS is validated on NEMA phantoms and patient data, for three PET geometries, compared to Monte-Carlo simulations and two proprietary vendor-specific reconstruction platforms. The reconstructed images have similar contrast recovery and background variability, deviating by up to 3.7%-point on contrast recovery and 1.8 on background variability and looking visually similar.

CONCLUSION

We have developed and validated an open-source scatter estimation library to complement reconstruction frameworks. By enabling vendor-independent clinical-grade reconstructions on custom scanner geometries, openSSS represents a crucial step in transparent research on quantitative PET and novel PET scanner designs.

Initial results of the Hyperion IIDPET insert for simultaneous PET-MRI applied to atherosclerotic plaque imaging in New-Zealand white rabbits

Objective.In preclinical research,in vivoimaging of mice and rats is more common than any other animal species, since their physiopathology is very well-known and many genetically altered disease models exist. Animal studies based on small rodents are usually performed using dedicated preclinical imaging systems with high spatial resolution. For studies that require animal models such as mini-pigs or New-Zealand White (NZW) rabbits, imaging systems with larger bore sizes are required. In case of hybrid imaging using positron emission tomography (PET) and magnetic resonance imaging (MRI), clinical systems have to be used, as these animal models do not typically fit in preclinical simultaneous PET-MRI scanners.Approach.In this paper, we present initial imaging results obtained with the Hyperion IIDPET insert which can accommodate NZW rabbits when combined with a large volume MRI RF coil. First, we developed a rabbit-sized image quality phantom of comparable size to a NZW rabbit in order to evaluate the PET imaging performance of the insert under high count rates. For this phantom, radioactive spheres with inner diameters between 3.95 and7.86mm were visible in a warm background with a tracer activity ratio of 4.1 to 1 and with a total18F activity in the phantom of58MBq at measurement start. Second, we performed simultaneous PET-MR imaging of atherosclerotic plaques in a rabbitin vivousing a single injection containing18F-FDG for detection of inflammatory activity, and Gd-ESMA for visualization of the aortic vessel wall and plaques with MRI.Main results.The fused PET-MR images reveal18F-FDG uptake within an active plaques with plaque thicknesses in the sub-millimeter range. Histology showed colocalization of18F-FDG uptake with macrophages in the aortic vessel wall lesions.Significance.Our initial results demonstrate that this PET insert is a promising system for simultaneous high-resolution PET-MR atherosclerotic plaque imaging studies in NZW rabbits.

Simultaneous Activity and Attenuation Estimation in TOF-PET With TV-Constrained Nonconvex Optimization

An alternating direction method of multipliers (ADMM) framework is developed for nonsmooth biconvex optimization for inverse problems in imaging. In particular, the simultaneous estimation of activity and attenuation (SAA) problem in time-of-flight positron emission tomography (TOF-PET) has such a structure when maximum likelihood estimation (MLE) is employed. The ADMM framework is applied to MLE for SAA in TOF-PET, resulting in the ADMM-SAA algorithm. This algorithm is extended by imposing total variation (TV) constraints on both the activity and attenuation map, resulting in the ADMM-TVSAA algorithm. The performance of this algorithm is illustrated using the penalized maximum likelihood activity and attenuation estimation (P-MLAA) algorithm as a reference.

Deep learning techniques in PET/CT imaging: A comprehensive review from sinogram to image space

Positron emission tomography/computed tomography (PET/CT) is increasingly used in oncology, neurology, cardiology, and emerging medical fields. The success stems from the cohesive information that hybrid PET/CT imaging offers, surpassing the capabilities of individual modalities when used in isolation for different malignancies. However, manual image interpretation requires extensive disease-specific knowledge, and it is a time-consuming aspect of physicians' daily routines. Deep learning algorithms, akin to a practitioner during training, extract knowledge from images to facilitate the diagnosis process by detecting symptoms and enhancing images. This acquired knowledge aids in supporting the diagnosis process through symptom detection and image enhancement. The available review papers on PET/CT imaging have a drawback as they either included additional modalities or examined various types of AI applications. However, there has been a lack of comprehensive investigation specifically focused on the highly specific use of AI, and deep learning, on PET/CT images. This review aims to fill that gap by investigating the characteristics of approaches used in papers that employed deep learning for PET/CT imaging. Within the review, we identified 99 studies published between 2017 and 2022 that applied deep learning to PET/CT images. We also identified the best pre-processing algorithms and the most effective deep learning models reported for PET/CT while highlighting the current limitations. Our review underscores the potential of deep learning (DL) in PET/CT imaging, with successful applications in lesion detection, tumor segmentation, and disease classification in both sinogram and image spaces. Common and specific pre-processing techniques are also discussed. DL algorithms excel at extracting meaningful features, and enhancing accuracy and efficiency in diagnosis. However, limitations arise from the scarcity of annotated datasets and challenges in explainability and uncertainty. Recent DL models, such as attention-based models, generative models, multi-modal models, graph convolutional networks, and transformers, are promising for improving PET/CT studies. Additionally, radiomics has garnered attention for tumor classification and predicting patient outcomes. Ongoing research is crucial to explore new applications and improve the accuracy of DL models in this rapidly evolving field.

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

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

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.

Using domain knowledge for robust and generalizable deep learning-based CT-free PET attenuation and scatter correction

Despite the potential of deep learning (DL)-based methods in substituting CT-based PET attenuation and scatter correction for CT-free PET imaging, a critical bottleneck is their limited capability in handling large heterogeneity of tracers and scanners of PET imaging. This study employs a simple way to integrate domain knowledge in DL for CT-free PET imaging. In contrast to conventional direct DL methods, we simplify the complex problem by a domain decomposition so that the learning of anatomy-dependent attenuation correction can be achieved robustly in a low-frequency domain while the original anatomy-independent high-frequency texture can be preserved during the processing. Even with the training from one tracer on one scanner, the effectiveness and robustness of our proposed approach are confirmed in tests of various external imaging tracers on different scanners. The robust, generalizable, and transparent DL development may enhance the potential of clinical translation.

Deep learning-based methods have been proposed to substitute CT-based PET attenuation and scatter correction to achieve CT-free PET imaging. Here, the authors present a simple way to integrate domain knowledge in deep learning for CT-free PET imaging.

Lutetium background radiation in total-body PET-A simulation study on opportunities and challenges in PET attenuation correction

The current generation of total-body positron emission tomography (PET) scanners offer significant sensitivity increase with an extended axial imaging extent. With the large volume of lutetium-based scintillation crystals that are used as detector elements in these scanners, there is an increased flux of background radiation originating from ¹⁷⁶Lu decay in the crystals and higher sensitivity for detecting it. Combined with the ability of scanning the entire body in a single bed position, this allows more effective utilization of the lutetium background as a transmission source for estimating 511 keV attenuation coefficients. In this study, utilization of the lutetium background radiation for attenuation correction in total-body PET was studied using Monte Carlo simulations of a 3D whole-body XCAT phantom in the uEXPLORER PET scanner, with particular focus on ultralow-dose PET scans that are now made possible with these scanners. Effects of an increased acceptance angle, reduced scan durations, and Compton scattering on PET quantification were studied. Furthermore, quantification accuracy of lutetium-based attenuation correction was compared for a 20-min scan of the whole body on the uEXPLORER, a one-meter-long, and a conventional 24-cm-long scanner. Quantification and lesion contrast were minimally affected in both long axial field-of-view scanners and in a whole-body 20-min scan, the mean bias in all analyzed organs of interest were within a ±10% range compared to ground-truth activity maps. Quantification was affected in certain organs, when scan duration was reduced to 5 min or a reduced acceptance angle of 17° was used. Analysis of the Compton scattered events suggests that implementing a scatter correction method for the transmission data will be required, and increasing the energy threshold from 250 keV to 290 keV can reduce the computational costs and data rates, with negligible effects on PET quantification. Finally, the current results can serve as groundwork for transferring lutetium-based attenuation correction into research and clinical practice.

Comparison of deep learning-based emission-only attenuation correction methods for positron emission tomography

PURPOSE

This study aims to compare two approaches using only emission PET data and a convolution neural network (CNN) to correct the attenuation (μ) of the annihilation photons in PET.

METHODS

One of the approaches uses a CNN to generate μ-maps from the non-attenuation-corrected (NAC) PET images (μ-CNN_NAC). In the other method, CNN is used to improve the accuracy of μ-maps generated using maximum likelihood estimation of activity and attenuation (MLAA) reconstruction (μ-CNN_MLAA). We investigated the improvement in the CNN performance by combining the two methods (μ-CNN_MLAA+NAC) and the suitability of μ-CNN_NAC for providing the scatter distribution required for MLAA reconstruction. Image data from ¹⁸F-FDG (n = 100) or ⁶⁸ Ga-DOTATOC (n = 50) PET/CT scans were used for neural network training and testing.

RESULTS

The error of the attenuation correction factors estimated using μ-CT and μ-CNN_NAC was over 7%, but that of scatter estimates was only 2.5%, indicating the validity of the scatter estimation from μ-CNN_NAC. However, CNN_NAC provided less accurate bone structures in the μ-maps, while the best results in recovering the fine bone structures were obtained by applying CNN_MLAA+NAC. Additionally, the μ-values in the lungs were overestimated by CNN_NAC. Activity images (λ) corrected for attenuation using μ-CNN_MLAA and μ-CNN_MLAA+NAC were superior to those corrected using μ-CNN_NAC, in terms of their similarity to λ-CT. However, the improvement in the similarity with λ-CT by combining the CNN_NAC and CNN_MLAA approaches was insignificant (percent error for lung cancer lesions, λ-CNN_NAC = 5.45% ± 7.88%; λ-CNN_MLAA = 1.21% ± 5.74%; λ-CNN_MLAA+NAC = 1.91% ± 4.78%; percent error for bone cancer lesions, λ-CNN_NAC = 1.37% ± 5.16%; λ-CNN_MLAA = 0.23% ± 3.81%; λ-CNN_MLAA+NAC = 0.05% ± 3.49%).

CONCLUSION

The use of CNN_NAC was feasible for scatter estimation to address the chicken-egg dilemma in MLAA reconstruction, but CNN_MLAA outperformed CNN_NAC.

A CT-less approach to quantitative PET imaging using the LSO intrinsic radiation for long-axial FOV PET scanners

Purpose

Long‐axial field‐of‐view (FOV) positron emission tomography (PET) scanners have gained a lot of interest in the recent years. Such scanners provide increased sensitivity and enable unique imaging opportunities that were not previously feasible. Benefiting from the high sensitivity of a long‐axial FOV PET scanner, we studied a computed tomography (CT)–less reconstruction algorithm for the Siemens Biograph Vision Quadra with an axial FOV of 106 cm.

Methods

In this work, the background radiation from radioisotope lutetium‐176 in the scintillators was used to create an initial estimate of the attenuation maps. Then, joint activity and attenuation reconstruction algorithms were used to create an improved attenuation map of the object. The final attenuation maps were then used to reconstruct quantitative PET images, which were compared against CT‐based PET images. The proposed method was evaluated on data from three patients who underwent a flurodeoxyglucouse PET scan.

Results

Segmentation of the PET images of the three studied patients showed an average quantitative error of 6.5%–8.3% across all studied organs when using attenuation maps from maximum likelihood estimation of attenuation and activity and 5.3%–6.6% when using attenuation maps from maximum likelihood estimation of activity and attenuation correction coefficients.

Conclusions

Benefiting from the background radiation of lutetium‐based scintillators, a quantitative CT‐less PET imaging technique was evaluated in this work.

Accurate Transmission-Less Attenuation Correction Method for Amyloid-β Brain PET Using Deep Neural Network

The lack of physically measured attenuation maps (μ-maps) for attenuation and scatter correction is an important technical challenge in brain-dedicated stand-alone positron emission tomography (PET) scanners. The accuracy of the calculated attenuation correction is limited by the nonuniformity of tissue composition due to pathologic conditions and the complex structure of facial bones. The aim of this study is to develop an accurate transmission-less attenuation correction method for amyloid-β (Aβ) brain PET studies. We investigated the validity of a deep convolutional neural network trained to produce a CT-derived μ-map (μ-CT) from simultaneously reconstructed activity and attenuation maps using the MLAA (maximum likelihood reconstruction of activity and attenuation) algorithm for Aβ brain PET. The performance of three different structures of U-net models (2D, 2.5D, and 3D) were compared. The U-net models generated less noisy and more uniform μ-maps than MLAA μ-maps. Among the three different U-net models, the patch-based 3D U-net model reduced noise and cross-talk artifacts more effectively. The Dice similarity coefficients between the μ-map generated using 3D U-net and μ-CT in bone and air segments were 0.83 and 0.67. All three U-net models showed better voxel-wise correlation of the μ-maps compared to MLAA. The patch-based 3D U-net model was the best. While the uptake value of MLAA yielded a high percentage error of 20% or more, the uptake value of 3D U-nets yielded the lowest percentage error within 5%. The proposed deep learning approach that requires no transmission data, anatomic image, or atlas/template for PET attenuation correction remarkably enhanced the quantitative accuracy of the simultaneously estimated MLAA μ-maps from Aβ brain PET.

Imitation learning for improved 3D PET/MR attenuation correction

The assessment of the quality of synthesised/pseudo Computed Tomography (pCT) images is commonly measured by an intensity-wise similarity between the ground truth CT and the pCT. However, when using the pCT as an attenuation map (μ-map) for PET reconstruction in Positron Emission Tomography Magnetic Resonance Imaging (PET/MRI) minimising the error between pCT and CT neglects the main objective of predicting a pCT that when used as μ-map reconstructs a pseudo PET (pPET) which is as similar as possible to the gold standard CT-derived PET reconstruction. This observation motivated us to propose a novel multi-hypothesis deep learning framework explicitly aimed at PET reconstruction application. A convolutional neural network (CNN) synthesises pCTs by minimising a combination of the pixel-wise error between pCT and CT and a novel metric-loss that itself is defined by a CNN and aims to minimise consequent PET residuals. Training is performed on a database of twenty 3D MR/CT/PET brain image pairs. Quantitative results on a fully independent dataset of twenty-three 3D MR/CT/PET image pairs show that the network is able to synthesise more accurate pCTs. The Mean Absolute Error on the pCT (110.98 HU ± 19.22 HU) compared to a baseline CNN (172.12 HU ± 19.61 HU) and a multi-atlas propagation approach (153.40 HU ± 18.68 HU), and subsequently lead to a significant improvement in the PET reconstruction error (4.74% ± 1.52% compared to baseline 13.72% ± 2.48% and multi-atlas propagation 6.68% ± 2.06%).

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

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

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.

Attenuation correction for human PET/MRI studies

Attenuation correction has been one of the main methodological challenges in the integrated positron emission tomography and magnetic resonance imaging (PET/MRI) field. As standard transmission or computed tomography approaches are not available in integrated PET/MRI scanners, MR-based attenuation correction approaches had to be developed. Aspects that have to be considered for implementing accurate methods include the need to account for attenuation in bone tissue, normal and pathological lung and the MR hardware present in the PET field-of-view, to reduce the impact of subject motion, to minimize truncation and susceptibility artifacts, and to address issues related to the data acquisition and processing both on the PET and MRI sides. The standard MR-based attenuation correction techniques implemented by the PET/MRI equipment manufacturers and their impact on clinical and research PET data interpretation and quantification are first discussed. Next, the more advanced methods, including the latest generation deep learning-based approaches that have been proposed for further minimizing the attenuation correction related bias are described. Finally, a future perspective focused on the needed developments in the field is given.

Intrascanner Reproducibility of an SPM-based Head MR-based Attenuation Correction Method

Recently, an exhaustive examination of 11 state of the art MR-based attenuation correction (AC) concluded that there are currently a few methods showing similar results compared to the gold-standard, CT-based AC. While the study presented a thorough portfolio of metrics to quantify accuracy (bias) and quality, it lacked one of the most important metrics to quantify robustness that is critical for its clinical applicability: intrascanner reproducibility (repeatability). In this work, we provide for the first time a study of the repeatability of one of the outperforming brain MR-based AC methods: the SPM-based pseudo-CT approach. 22 subjects undergoing 3 ¹⁸F-FDG PET/MRI visits within 2 months were retrospectively analyzed in this study. Pseudo-CT mu-maps were obtained from the coregistered MR images for all 3 visits and the PET data from visit 1 was reconstructed using all three mu-maps. Relative changes (RC), Intraclass correlation coefficient (ICC), Reproducibility coefficient (RDC_95%) and Bland-Altman Limits of Agreement (LoA) were used to measure repeatability. Voxel-based and ROI-based results showed that absolute RC for the reconstructed PET images are within ~2%. The brain cortex and the cerebellum were the regions with the largest variability (~3%). The differences across visits were not statistically significant (p=0.90). In conclusion this study shows for the first time the repeatability of the SPM-based pseudo-CT approach for brain MR-AC. These results, in addition to the ease of implementation and the quality and robustness previously demonstrated, confer this SPM-based method an ideal candidate for routine brain PET/MRI research and clinical studies.

Practical issues and limitations of brain attenuation correction on a simultaneous PET-MR scanner

BACKGROUND

Despite the advent of clinical PET-MR imaging for routine use in 2011 and the development of several methods to address the problem of attenuation correction, some challenges remain. We have identified and investigated several issues that might affect the reliability and accuracy of current attenuation correction methods when these are implemented for clinical and research studies of the brain. These are (1) the accuracy of converting CT Hounsfield units, obtained from an independently acquired CT scan, to 511 keV linear attenuation coefficients; (2) the effect of padding used in the MR head coil; (3) the presence of close-packed hair; (4) the effect of headphones. For each of these, we have examined the effect on reconstructed PET images and evaluated practical mitigating measures.

RESULTS

Our major findings were (1) for both Siemens and GE PET-MR systems, CT data from either a Siemens or a GE PET-CT scanner may be used, provided the conversion to 511 keV μ-map is performed by the PET-MR vendor's own method, as implemented on their PET-CT scanner; (2) the effect of the head coil pads is minimal; (3) the effect of dense hair in the field of view is marked (> 10% error in reconstructed PET images); and (4) using headphones and not including them in the attenuation map causes significant errors in reconstructed PET images, but the risk of scanning without them may be acceptable following sound level measurements.

CONCLUSIONS

It is important that the limitations of attenuation correction in PET-MR are considered when designing research and clinical PET-MR protocols in order to enable accurate quantification of brain PET scans. Whilst the effect of pads is not significant, dense hair, the use of headphones and the use of an independently acquired CT-scan can all lead to non-negligible effects on PET quantification. Although seemingly trivial, these effects add complications to setting up protocols for clinical and research PET-MR studies that do not occur with PET-CT. In the absence of more sophisticated PET-MR brain attenuation correction, the effect of all of the issues above can be minimised if the pragmatic approaches presented in this work are followed.

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.

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.

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.

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.

A head coil system with an integrated orbiting transmission point source mechanism for attenuation correction in PET/MRI

The combination of positron emission tomography (PET) and magnetic resonance imaging (MRI) provides a benefit for diagnostic imaging. Still, attenuation correction (AC) is a challenge in PET/MRI compared to stand-alone PET and PET-computed tomography (PET/CT). In the absence of photonic transmission sources, AC in PET/MRI is usually based on retrospective segmentation of MR images or complex additional MR-sequences. However, most methods available today are still challenged by either the incorporation of cortical bone or substantial anatomical variations of subjects. This leads to a bias in quantification of tracer concentration in PET. Therefore, we have developed a fully integrated transmission source system for PET/MRI of the head to enable direct measurement of attenuation coefficients using external positron emitters, which is the reference standard in AC. Based on a setup called the 'liquid drive' presented by Jones et al (1995) two decades ago, we built a head coil system consisting of an MR-compatible hydraulic system driving a point source on a helical path around a 24-channel MR-receiver coil to perform a transmission scan. Sinogram windowing of the moving source allows for post-injection measurements. The prototype was tested in the Siemens Biograph mMR using a homogeneous water phantom and a phantom with air cavities and a Teflon (PTFE) cylinder. The second phantom was measured both with and without emission activity. For both measurements air, water and Teflon were clearly distinguishable and homogeneous regions of the phantom were successfully reproduced in the AC map. For water the linear attenuation coefficient was measured as (0.096  ±  0.005) cm^-1 in accordance with the true physical value. This combined MR head coil and transmission source system is, to our knowledge, the first working example to use an orbiting point source in PET/MRI and may be helpful in providing fully-quantitative PET data in neuro-PET/MRI.

Accurate hybrid template-based and MR-based attenuation correction using UTE images for simultaneous PET/MR brain imaging applications

Background

Attenuation correction is one of the most crucial correction factors for accurate PET data quantitation in hybrid PET/MR scanners, and computing accurate attenuation coefficient maps from MR brain acquisitions is challenging. Here, we develop a method for accurate bone and air segmentation using MR ultrashort echo time (UTE) images.

Methods

MR UTE images from simultaneous MR and PET imaging of five healthy volunteers was used to generate a whole head, bone and air template image for inclusion into an improved MR derived attenuation correction map, and applied to PET image data for quantitative analysis. Bone, air and soft tissue were segmented based on Gaussian Mixture Models with probabilistic tissue maps as a priori information. We present results for two approaches for bone attenuation coefficient assignments: one using a constant attenuation correction value; and another using an estimated continuous attenuation value based on a calibration fit. Quantitative comparisons were performed to evaluate the accuracy of the reconstructed PET images, with respect to a reference image reconstructed with manually segmented attenuation maps.

Results

The DICE coefficient analysis for the air and bone regions in the images demonstrated improvements compared to the UTE approach, and other state-of-the-art techniques. The most accurate whole brain and regional brain analyses were obtained using constant bone attenuation coefficient values.

Conclusions

A novel attenuation correction method for PET data reconstruction is proposed. Analyses show improvements in the quantitative accuracy of the reconstructed PET images compared to other state-of-the-art AC methods for simultaneous PET/MR scanners. Further evaluation is needed with radiopharmaceuticals other than FDG, and in larger cohorts of participants.

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.

Hybrid cardiac imaging using PET/MRI: a joint position statement by the European Society of Cardiovascular Radiology (ESCR) and the European Association of Nuclear Medicine (EANM)

Positron emission tomography (PET) and magnetic resonance imaging (MRI) have both been used for decades in cardiovascular imaging. Since 2010, hybrid PET/MRI using sequential and integrated scanner platforms has been available, with hybrid cardiac PET/MR imaging protocols increasingly incorporated into clinical workflows. Given the range of complementary information provided by each method, the use of hybrid PET/MRI may be justified and beneficial in particular clinical settings for the evaluation of different disease entities. In the present joint position statement, we critically review the role and value of integrated PET/MRI in cardiovascular imaging, provide a technical overview of cardiac PET/MRI and practical advice related to the cardiac PET/MRI workflow, identify cardiovascular applications that can potentially benefit from hybrid PET/MRI, and describe the needs for future development and research. In order to encourage its wide dissemination, this article is freely accessible on the European Radiology and European Journal of Hybrid Imaging web sites.

KEY POINTS

• Studies and case-reports indicate that PET/MRI is a feasible and robust technology. • Promising fields of application include a variety of cardiac conditions. • Larger studies are required to demonstrate its incremental and cost-effective value. • The translation of novel radiopharmaceuticals and MR-sequences will provide exciting new opportunities.

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.

Feasibility of Deep Learning-Based PET/MR Attenuation Correction in the Pelvis Using Only Diagnostic MR Images

This study evaluated the feasibility of using only diagnostically relevant magnetic resonance (MR) images together with deep learning for positron emission tomography (PET)/MR attenuation correction (deepMRAC) in the pelvis. Such an approach could eliminate dedicated MRAC sequences that have limited diagnostic utility but can substantially lengthen acquisition times for multibed position scans. We used axial T2 and T1 LAVA Flex magnetic resonance imaging images that were acquired for diagnostic purposes as inputs to a 3D deep convolutional neural network. The network was trained to produce a discretized (air, water, fat, and bone) substitute computed tomography (CT) (CTsub). Discretized (CTref-discrete) and continuously valued (CTref) reference CT images were created to serve as ground truth for network training and attenuation correction, respectively. Training was performed with data from 12 subjects. CTsub, CTref, and the system MRAC were used for PET/MR attenuation correction, and quantitative PET values of the resulting images were compared in 6 test subjects. Overall, the network produced CTsub with Dice coefficients of 0.79 ± 0.03 for cortical bone, 0.98 ± 0.01 for soft tissue (fat: 0.94 ± 0.0; water: 0.88 ± 0.02), and 0.49 ± 0.17 for bowel gas when compared with CTref-discrete. The root mean square error of the whole PET image was 4.9% by using deepMRAC and 11.6% by using the system MRAC. In evaluating 16 soft tissue lesions, the distribution of errors for maximum standardized uptake value was significantly narrower using deepMRAC (-1.0% ± 1.3%) than using system MRAC method (0.0% ± 6.4%) according to the Brown-Forsy the test (P < .05). These results indicate that improved PET/MR attenuation correction can be achieved in the pelvis using only diagnostically relevant MR images.

Hybrid PET/MRI Methodology

The hybrid PET/MR scanner represents the first implementation of the effective integration of two modalities allowing truly synchronous/simultaneous acquisition of their imaging signals. This integration, resulting from the innovation and development of specific hardware components has paved the way for new approaches in the study of neurodegenerative diseases. This chapter will describe the hardware development that has led to the availability of different clinical solutions for PET/MR imaging as well as the still-open technological challenges and opportunities related to the processing and exploitation of the simultaneous acquisition in neurological studies.

From simultaneous to synergistic MR-PET brain imaging: A review of hybrid MR-PET imaging methodologies

Simultaneous Magnetic Resonance Imaging (MRI) and Positron Emission Tomography (PET) scanning is a recent major development in biomedical imaging. The full integration of the PET detector ring and electronics within the MR system has been a technologically challenging design to develop but provides capacity for simultaneous imaging and the potential for new diagnostic and research capability. This article reviews state-of-the-art MR-PET hardware and software, and discusses future developments focusing on neuroimaging methodologies for MR-PET scanning. We particularly focus on the methodologies that lead to an improved synergy between MRI and PET, including optimal data acquisition, PET attenuation and motion correction, and joint image reconstruction and processing methods based on the underlying complementary and mutual information. We further review the current and potential future applications of simultaneous MR-PET in both systems neuroscience and clinical neuroimaging research. We demonstrate a simultaneous data acquisition protocol to highlight new applications of MR-PET neuroimaging research studies.

Zero TE-based pseudo-CT image conversion in the head and its application in PET/MR attenuation correction and MR-guided radiation therapy planning

PURPOSE

To describe a method for converting Zero TE (ZTE) MR images into X-ray attenuation information in the form of pseudo-CT images and demonstrate its performance for (1) attenuation correction (AC) in PET/MR and (2) dose planning in MR-guided radiation therapy planning (RTP).

METHODS

Proton density-weighted ZTE images were acquired as input for MR-based pseudo-CT conversion, providing (1) efficient capture of short-lived bone signals, (2) flat soft-tissue contrast, and (3) fast and robust 3D MR imaging. After bias correction and normalization, the images were segmented into bone, soft-tissue, and air by means of thresholding and morphological refinements. Fixed Hounsfield replacement values were assigned for air (-1000 HU) and soft-tissue (+42 HU), whereas continuous linear mapping was used for bone.

RESULTS

The obtained ZTE-derived pseudo-CT images accurately resembled the true CT images (i.e., Dice coefficient for bone overlap of 0.73 ± 0.08 and mean absolute error of 123 ± 25 HU evaluated over the whole head, including errors from residual registration mismatches in the neck and mouth regions). The linear bone mapping accounted for bone density variations. Averaged across five patients, ZTE-based AC demonstrated a PET error of -0.04 ± 1.68% relative to CT-based AC. Similarly, for RTP assessed in eight patients, the absolute dose difference over the target volume was found to be 0.23 ± 0.42%.

CONCLUSION

The described method enables MR to pseudo-CT image conversion for the head in an accurate, robust, and fast manner without relying on anatomical prior knowledge. Potential applications include PET/MR-AC, and MR-guided RTP.

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.

Comparison of volumetric and functional parameters in simultaneous cardiac PET/MR: feasibility of volumetric assessment with residual activity from prior PET/CT

OBJECTIVE

To compare cardiac left ventricular (LV) parameters in simultaneously acquired hybrid fluorine-18-fluorodeoxyglucose ([18F] FDG) positron emission tomography/magnetic resonance imaging (PET/MRI) in patients with residual tracer activity of upstream PET/CT.

METHODS

Twenty-nine patients (23 men, age 58±17 years) underwent cardiac PET/MRI either directly after a non-cardiac PET/CT with homogenous cardiac [18F] FDG uptake (n=20) or for viability assessment (n=9). Gated cardiac [18F] FDG PET and cine MR sequences were acquired simultaneously and evaluated blinded to the cross-imaging results. Image quality (IQ), end-diastolic (LVEDV), end-systolic volume (LVESV), ejection fraction (LVEF) and myocardial mass (LVMM) were measured. Pearson correlation and intraclass correlation coefficient (ICC), regression and a Bland-Altman analysis were assessed.

RESULTS

Except LVMM, volumetric and functional LV parameters demonstrated high correlations (LVESV: r=0.97, LVEDV: r=0.95, LVEF: r=0.91, LVMM: r=0.87, each p<0.05), but wide limits of agreement (LOA) for LVEDV (-25.3-82.5ml); LVESV (-33.1-72.7ml); LVEF (-18.9-14.8%) and LVMM (-78.2-43.2g). Intra- and interobserver reliability were very high (ICC≥0.95) for all parameters, except for MR-LVEF (ICC=0.87). PET-IQ (0-3) was high (mean: 2.2±0.9) with significant influence on LVMM calculations only.

CONCLUSION

In simultaneously acquired cardiac PET/MRI data, LVEDV, LVESV and LVEF show good agreement. However, the agreement seems to be limited if cardiac PET/MRI follows PET/CT and only the residual activity is used.

KEY POINTS

• [ 18 F] FDG PET-MRI is feasible with residual [ 18 F] FDG activity in patients with homogenous cardiac uptake. • Cardiac volumes and function assessed by PET/MRI show good agreement. • LVEDV and LVESV are underestimated; PET overestimates LVMM and LVEF. • Cardiac PET and MRI data correlate better when acquired simultaneously than sequentially. • PET and MRI should not assess LV parameters interchangeably.

Dual nano-sized contrast agents in PET/MRI: a systematic review

Nowadays molecular imaging plays a vital role in achieving a successful targeted and personalized treatment. Hence, the approach of combining two or more medical imaging modalities was developed. The objective of this review is to systematically compare recent dual contrast agents in Positron Emission Tomography (PET)/Magnetic Resonance Imaging (MRI) and in some cases Single photon emission computed tomography (SPECT)/MRI in terms of some their characteristics, such as tumor uptake, and reticuloendothelial system uptake (especially liver) and their relaxivity rates for early detection of primary cancer tumor. To the best of our knowledge, this is the first systematic and integrated overview of this field. Two reviewers individually directed the systematic review search using PubMed, MEDLINE and Google Scholar. Two other reviewers directed quality assessment, using the criteria checklist from the CAMARADES (Collaborative Approach to Meta-Analysis and Review of Animal Data from Experimental Studies) tool, and differences were resolved by consensus. After reviewing all 49 studies, we concluded that a size range of 20-200 nm can be used for molecular imaging, although it is better to try to achieve as small a size as it is possible. Also, small nanoparticles with a hydrophilic coating and positive charge are suitable as a T₂ contrast agent. According to our selected data, the most successful dual probes in terms of high targeting were with an average size of 40 nm, PEGylated using peptides as a biomarker and radiolabeled with copper 64 and gallium 68. Copyright © 2017 John Wiley & Sons, Ltd.

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.

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.