Performance comparison of two resolution modeling PET reconstruction algorithms in terms of physical figures of merit used in quantitative imaging

Impact of list-mode reconstruction and image-space point spread function correction on PET image contrast and quantitative value using SiPM-based PET/CT system

We evaluate the effects of list-mode reconstruction and the image-space point spread function (iPSF) on the contrast and quantitative values of positron emission tomography (PET) images using a SiPM-PET/CT system. The evaluation is conducted on an NEMA body phantom and clinical images using a Cartesion Prime SiPM-PET/CT system. The signal-to-background ratio (SBR) of the phantom is set to 2, 4, 6, and 8, and all the PET image data are obtained and reconstructed using 3D-OSEM, time-of-flight, iPSF (-/ +), and a 4-mm Gaussian filter with several iterations. The evaluation criteria include % background variability (NB,10 mm), % contrast (QH,10 mm), iPSF change in QH,10 mm (ΔQH,10 mm) for edge artifact evaluation, profile curves, visual evaluation of edge artifacts, clinical imaging for the standardized uptake value (SUV) of lung nodules, and SNRliver. NB,10 mm demonstrates no significant difference in all SBRs with and without iPSF, whereas QH,10 mm is higher based on the SBR with and without iPSF. ΔQH,10 mm indicates increased iterations and a larger rate of change (> 5%) for small spheres of < 17 mm. The profile curves portrayed almost real concentrations, except for the 10-mm sphere of SBR2 without iPSF; however, with iPSF, an overshoot was observed in the 13-mm sphere of all SBRs. The degree of overshoot increased with increasing iteration and SBR. Edge artifacts were detected at values ≥ 17-22 mm in SBRs other than SBR2 with iPSF. Irrespective of the nodal size, SUV and SNRliver improved considerably after iPSF adjustment. Therefore, the effects of list-mode reconstruction and iPSF on PET image contrast were limited, and the overcorrection of the quantitative values was validated using iPSF.

The impact of time-of-flight, resolution recovery, and noise modelling in reconstruction algorithms in non-solid-state detectors PET/CT scanners: - multi-centric comparison of activity recovery in a 68Ge phantom

The reconstruction algorithms implemented on PET/CT scanners offer gain in activity recovery of small lesions at an extent that is not full known yet. METHODS: A cylindrical phantom with warm background and hot spheres filled with a ⁶⁸Ge epoxy was acquired with four non-state-solid-detectors PET/CT scanners: mCT, Ingenuity TF, Discovery 710, and IQ. Images were reconstructed switching on and off time-of-flight (TOF), point spread function (PSF) modelling, and Bayesian penalised likelihood (BPL). Images were reconstructed with the default parameters recommended by the manufacturers. The recovery coefficient (RC_max), defined as the ratio of the measured maximum activity concentration in each sphere and the actual one, and the coefficient of variation (CoV_BAC) defined as the ratio of the standard deviation and the average of background activity concentration were measured. RESULTS: While with IR alone, complete recovery of the activity concentration is achieved down to the 22 mm diameter's sphere, with TOF, TOF + PSF and BPL it is achieved down to the 17 mm diameter one. At smaller dimensions, the difference among the various studied reconstruction algorithms is substantial for the 13- and 17-mm diameters' spheres for all scanners and for all reconstructions with a considerable gain in RCmax when PSF and BPL are used. At 10 mm diameter's sphere the difference among the algorithms is significantly reduced, except for BPL which still guarantees a gain in RCmax.

Performance evaluation of a new time of flight PET/CT scanner: Results of a multicenter study

PURPOSE

The aim of this multicenter study was to evaluate the performance of the upgraded version of the Ingenuity TF PET/CT scanner, according to the NEMA NU-2 2012 standards.

METHODS

Spatial resolution, sensitivity, count rate response, scatter fraction, image quality and accuracy were evaluated on three Ingenuity TF scanners installed in Italian hospitals. Furthermore, energy and timing resolution were measured. A detailed image quality phantom analysis was performed to evaluate the effect of different clinical reconstruction parameters, including the application of PSF correction.

RESULTS

Results show an average spatial resolution of 4.7 mm and an average absolute system sensitivity of 7.9 cps/kBq. The average maximum NECR was 119.83 kcps at 20.67 kBq/ml, while the maximum true event rate was 322.62 kcps at the concentration of 24.51 kBq/ml. The average maximum bias below NECR peak was 12.58%. All the results of NEMA tests were in agreement with the values declared by the manufacturer. The estimated average energy and timing resolution were 10.83% and 536.2 ps, respectively. Image quality phantom analysis obtained with different reconstruction settings showed that PSF correction was the parameter that affected mainly on contrast recovery coefficient, while the iteration number and amplitude of Gaussian filter had no significant effect. Of relevance, the application of PSF correction never led to recovery coefficient values higher than 100% and to Gibbs or edge artifacts.

CONCLUSIONS

The new Ingenuity TF model shows physical performance similar to other scanners of the latest generation for all standard NEMA NU2-2012 measurements.

Task Group 174 Report: Utilization of [18 F]Fluorodeoxyglucose Positron Emission Tomography ([18 F]FDG-PET) in Radiation Therapy

The use of positron emission tomography (PET) in radiation therapy (RT) is rapidly increasing in the areas of staging, segmentation, treatment planning, and response assessment. The most common radiotracer is ¹⁸ F-fluorodeoxyglucose ([¹⁸ F]FDG), a glucose analog with demonstrated efficacy in cancer diagnosis and staging. However, diagnosis and RT planning are different endeavors with unique requirements, and very little literature is available for guiding physicists and clinicians in the utilization of [¹⁸ F]FDG-PET in RT. The two goals of this report are to educate and provide recommendations. The report provides background and education on current PET imaging systems, PET tracers, intensity quantification, and current utilization in RT (staging, segmentation, image registration, treatment planning, and therapy response assessment). Recommendations are provided on acceptance testing, annual and monthly quality assurance, scanning protocols to ensure consistency between interpatient scans and intrapatient longitudinal scans, reporting of patient and scan parameters in literature, requirements for incorporation of [¹⁸ F]FDG-PET in treatment planning systems, and image registration. The recommendations provided here are minimum requirements and are not meant to cover all aspects of the use of [¹⁸ F]FDG-PET for RT.

Differences in edge artifacts between 68Ga- and 18F-PET images reconstructed using point spread function correction

OBJECTIVE

Edge artifacts have been reported on in relation to F-PET using point spread function correction algorithms. The positron range of Ga is longer than F, and this difference is thought to result in different edge artifacts. The purpose of this study is to clarify the difference in edge artifacts in PET images using point spread function correction in Ga- and F-PET.

METHODS

We used a National Electrical Manufacturers Association International Electrotechnical Commission body phantom. The phantom was filled severally with Ga and F solution. The PET data were obtained over a 90 minutes period using a True Point Biograph 16 scanner. The images were then reconstructed with the ordered subset expectation maximization with point spread function correction. The phantom image analyses were performed by a visual assessment of the PET images and profiles, and an absolute recovery coefficient, which was the ratio of the maximum radioactivity of any given hot sphere to its true radioactivity.

RESULTS

The ring-like edge artifacts of Ga-PET were less prominent than those in F-PET. The relative radioactivity profiles of Ga-PET showed low overshoots of the maximum radioactivity although high overshoots did appear in F-PET. The absolute recovery coefficients of Ga-PET were smaller than those of F-PET.

CONCLUSION

The edge artifacts of Ga-PET were less prominent than those of F-PET, and their overshoots were smaller. The difference in the positron range between Ga and F may possibly result in the difference in edge artifacts of images reconstructed using the point spread function correction algorithm.

Multicenter study of quantitative PET system harmonization using NIST-traceable 68Ge/68Ga cross-calibration kit

PURPOSE

The present study aimed to define the errors in SUV and demonstrate the feasibility of SUV harmonization among contemporary PET/CT scanners using a novel National Institute of Standards and Technology (NIST)-traceable ⁶⁸Ge/⁶⁸Ga source as the reference standard.

METHODS

We used ⁶⁸Ge/⁶⁸Ga dose calibrator and PET sources made with same batch of ⁶⁸Ge/⁶⁸Ga embedded in epoxy that is traceable to the NIST standard. Bias in the amount of radioactivity and the radioactive concentrations measured by the dose calibrators and PET/CT scanners, respectively, was determined at five Japanese sites. We adjusted optimal dial setting of the dose calibrators and PET reconstruction parameters to close the actual amount of radioactivity and the radioactive concentration, respectively, of the NIST-traceable ⁶⁸Ge/⁶⁸Ga sources to harmonize SUV. Errors in SUV before and after harmonization were then calculated at each site.

RESULTS

The average bias in the amount of radioactivity and the radioactive concentrations measured by dose calibrator and PET scanner was -4.94% and -12.22%, respectively, before, and -0.14% and -4.81%, respectively, after harmonization. Corresponding averaged errors in SUV measured under clinical conditions were underestimated by 7.66%, but improved by -4.70% under optimal conditions.

CONCLUSION

Our proposed method using an NIST-traceable ⁶⁸Ge/⁶⁸Ga source identified bias in values obtained using dose calibrators and PET scanners, and reduced SUV variability to within 5% across different models of PET scanners at five sites. Our protocol using a standard source has considerable potential for harmonizing the SUV when contemporary PET scanners are involved in multicenter studies.

Bayesian penalized-likelihood reconstruction algorithm suppresses edge artifacts in PET reconstruction based on point-spread-function

PURPOSE

The Bayesian penalized-likelihood reconstruction algorithm (BPL), Q.Clear, uses relative difference penalty as a regularization function to control image noise and the degree of edge-preservation in PET images. The present study aimed to determine the effects of suppression on edge artifacts due to point-spread-function (PSF) correction using a Q.Clear.

METHODS

Spheres of a cylindrical phantom contained a background of 5.3 kBq/mL of [¹⁸F]FDG and sphere-to-background ratios (SBR) of 16, 8, 4 and 2. The background also contained water and spheres containing 21.2 kBq/mL of [¹⁸F]FDG as non-background. All data were acquired using a Discovery PET/CT 710 and were reconstructed using three-dimensional ordered-subset expectation maximization with time-of-flight (TOF) and PSF correction (3D-OSEM), and Q.Clear with TOF (BPL). We investigated β-values of 200-800 using BPL. The PET images were analyzed using visual assessment and profile curves, edge variability and contrast recovery coefficients were measured.

RESULTS

The 38- and 27-mm spheres were surrounded by higher radioactivity concentration when reconstructed with 3D-OSEM as opposed to BPL, which suppressed edge artifacts. Images of 10-mm spheres had sharper overshoot at high SBR and non-background when reconstructed with BPL. Although contrast recovery coefficients of 10-mm spheres in BPL decreased as a function of increasing β, higher penalty parameter decreased the overshoot.

CONCLUSIONS

BPL is a feasible method for the suppression of edge artifacts of PSF correction, although this depends on SBR and sphere size. Overshoot associated with BPL caused overestimation in small spheres at high SBR. Higher penalty parameter in BPL can suppress overshoot more effectively.

Optimization, evaluation and calibration of a cross-strip DOI detector

This study depicts the evaluation of a SiPM detector with depth of interaction (DOI) capability via a dual-sided readout that is suitable for high-resolution positron emission tomography and magnetic resonance (PET/MR) imaging. Two different 12  ×  12 pixelated LSO scintillator arrays with a crystal pitch of 1.60 mm are examined. One array is 20 mm-long with a crystal separation by the specular reflector Vikuiti enhanced specular reflector (ESR), and the other one is 18 mm-long and separated by the diffuse reflector Lumirror E60 (E60). An improvement in energy resolution from 22.6% to 15.5% for the scintillator array with the E60 reflector is achieved by taking a nonlinear light collection correction into account. The results are FWHM energy resolutions of 14.0% and 15.5%, average FWHM DOI resolutions of 2.96 mm and 1.83 mm, and FWHM coincidence resolving times of 1.09 ns and 1.48 ns for the scintillator array with the ESR and that with the E60 reflector, respectively. The measured DOI signal ratios need to be assigned to an interaction depth inside the scintillator crystal. A linear and a nonlinear method, using the intrinsic scintillator radiation from lutetium, are implemented for an easy to apply calibration and are compared to the conventional method, which exploits a setup with an externally collimated radiation beam. The deviation between the DOI functions of the linear or nonlinear method and the conventional method is determined. The resulting average of differences in DOI positions is 0.67 mm and 0.45 mm for the nonlinear calibration method for the scintillator array with the ESR and with the E60 reflector, respectively; Whereas the linear calibration method results in 0.51 mm and 0.32 mm for the scintillator array with the ESR and the E60 reflector, respectively; and is, due to its simplicity, also applicable in assembled detector systems.

Diagnostic performance of 18F-FDG PET/CT using point spread function reconstruction on initial staging of rectal cancer: a comparison study with conventional PET/CT and pelvic MRI

Background

Accurate staging is crucial for treatment selection and prognosis prediction in patients with rectal cancer. Point spread function (PSF) reconstruction can improve spatial resolution and signal-to-noise ratio of PET imaging. The aim of this study was to evaluate the effectiveness of ¹⁸F-FDG PET/CT with PSF reconstruction for initial staging in rectal cancer compared with conventional PET/CT and pelvic MRI.

Methods

A total of 59 patients with rectal cancer underwent preoperative ¹⁸F-FDG PET/CT and pelvic MRI. The maximum standardized uptake value (SUVmax) and lesion to background (L/B) ratio of possible metastatic lymph nodes, and metabolic tumor volumes (MTVs) of primary tumors were calculated. For N and T (T1-2 vs T3-4) staging, sensitivities, specificities, positive predictive values, negative predictive values, and accuracies were compared between conventional PET/CT [reconstructed with ordered subset expectation maximization (OSEM)], PSF-PET/CT (reconstructed with OSEM+PSF), and pelvic MRI. Histopathologic analysis was the reference standard.

Results

For N staging, PSF-PET/CT provided higher sensitivity (78.6%) than conventional PET/CT (64.3%), and pelvic MRI (57.1%), and all techniques showed high specificity (PSF-PET: 95.4%, conventional PET: 96.7%, pelvic MRI: 93.5%). SUVmax and L/B ratio were significantly higher in PSF-PET/CT than conventional-PET/CT (p < 0.001). The accuracy for T staging in PSF-PET/CT (69.4%) was not significantly different to conventional PET/CT (73.5%) and pelvic MRI (73.5%). MTVs of PSF and conventional PET showed a significant difference among T stages (p < 0.001), with higher values in advanced stages. In M staging, both PSF and conventional PET/CT diagnosed all distant metastases correctly.

Conclusions

PSF-PET/CT produced images with higher lesion-to-background contrast than conventional PET/CT, which allowed improved detection of lymph node metastasis without compromising specificity, and showed comparable diagnostic value to MRI in local staging. PSF-PET/CT is likely to have a great value for initial staging in rectal cancer.

Optimization of PET/CT image quality using the GE 'Sharp IR' point-spread function reconstruction algorithm

OBJECTIVE

The objective of this study was to quantify any improvement with the GE 'Sharp IR' point-spread function (PSF) reconstruction algorithm in addition to ordered subsets expectation maximum (OSEM) and time-of-flight (TOF) reconstruction algorithms and establish the optimum parameters to be used in clinical studies.

MATERIALS AND METHODS

We conducted a range of experiments using the National Electrical Manufacturers Association image quality phantom filled with a 4 : 1 signal-to-background ratio. We scanned the phantom using the GE Discovery 690 PET/CT scanner. We varied iteration number and Gaussian filtration. Results were compared for OSEM, OSEM+TOF and OSEM+TOF+PSF reconstructions. A sample of 15 whole-body fluorine-18-fluorodeoxyglucose were reconstructed with OSEM+TOF and OSEM+TOF+PSF using a selection of optimum reconstruction parameters determined in phantom studies. Clinicians qualitatively ranked their preferred images to choose optimum parameters.

RESULTS

The addition of PSF improved signal-to-noise ratios (SNRs), contrast, hot contrast recovery coefficients and noise over OSEM and OSEM+TOF reconstruction algorithms. SNRs were the highest at two iterations and with 0 or 2 mm filters with OSEM+TOF+PSF reconstruction in all phantom studies. Clinicians generally favoured OSEM+TOF+PSF reconstruction with three iterations and a 2 mm filter.

CONCLUSION

PSF reconstruction significantly improved image quality for both clinical and phantom studies. We recommended the optimum reconstruction parameters using three iterations, 24 subsets and a 2 mm filter, which improved SNRs by up to 28.8% for small lesions (P<0.05).

Image improvement method for positron emission mammography

PURPOSE

To evaluate in clinical use a rapidly converging, efficient iterative deconvolution algorithm (RSEMD) for improving the quantitative accuracy of previously reconstructed breast images by a commercial positron emission mammography (PEM) scanner.

MATERIALS AND METHODS

The RSEMD method was tested on imaging data from clinical Naviscan Flex Solo II PEM scanner. This method was applied to anthropomorphic like breast phantom data and patient breast images previously reconstructed with Naviscan software to determine improvements in image resolution, signal to noise ratio (SNR) and contrast to noise ratio (CNR).

RESULTS

In all of the patients' breast studies the improved images proved to have higher resolution, contrast and lower noise as compared with images reconstructed by conventional methods. In general, the values of CNR reached a plateau at an average of 6 iterations with an average improvement factor of about 2 for post-reconstructed Flex Solo II PEM images. Improvements in image resolution after the application of RSEMD have also been demonstrated.

CONCLUSIONS

A rapidly converging, iterative deconvolution algorithm with a resolution subsets-based approach (RSEMD) that operates on patient DICOM images has been used for quantitative improvement in breast imaging. The RSEMD method can be applied to PEM images to enhance the resolution and contrast in cancer diagnosis to monitor the tumor progression at the earliest stages.

Point-spread function reconstructed PET images of sub-centimeter lesions are not quantitative

Background

PET image reconstruction methods include modeling of resolution degrading phenomena, often referred to as point-spread function (PSF) reconstruction. The aim of this study was to develop a clinically relevant phantom and characterize the reproducibility and accuracy of high-resolution PSF reconstructed images of small lesions, which is a prerequisite for using PET in the prediction and evaluation of responses to treatment.

Sets of small homogeneous ¹⁸F-spheres (range 3–12 mm diameter, relevant for small lesions and lymph nodes) were suspended and covered by a ¹¹C-silicone, which provided a scattering medium and a varying sphere-to-background ratio. Repeated measurements were made on PET/CT scanners from two vendors using a wide range of reconstruction parameters. Recovery coefficients (RCs) were measured for clinically used volume-of-interest definitions.

Results

For non-PSF images, RCs were reproducible and fell monotonically as the sphere diameter decreased, which is the expected behavior. PSF images converged slower and had artifacts: RCs did not fall monotonically as sphere diameters decreased but had a maximum RC for sphere sizes around 8 mm, RCs could be greater than 1, and RCs were less reproducible. To some degree, post-reconstruction filters could suppress PSF artifacts.

Conclusions

High-resolution PSF images of small lesions showed artifacts that could lead to serious misinterpretations when used for monitoring treatment response. Thus, it could be safer to use non-PSF reconstruction for quantitative purposes unless PSF reconstruction parameters are optimized for the specific task.

Electronic supplementary material

The online version of this article (doi:10.1186/s40658-016-0169-9) contains supplementary material, which is available to authorized users.

Evaluation of spatial dependence of point spread function-based PET reconstruction using a traceable point-like 22Na source

Background

The point spread function (PSF) of positron emission tomography (PET) depends on the position across the field of view (FOV). Reconstruction based on PSF improves spatial resolution and quantitative accuracy. The present study aimed to quantify the effects of PSF correction as a function of the position of a traceable point-like ²²Na source over the FOV on two PET scanners with a different detector design.

Methods

We used Discovery 600 and Discovery 710 (GE Healthcare) PET scanners and traceable point-like ²²Na sources (<1 MBq) with a spherical absorber design that assures uniform angular distribution of the emitted annihilation photons. The source was moved in three directions at intervals of 1 cm from the center towards the peripheral FOV using a three-dimensional (3D)-positioning robot, and data were acquired over a period of 2 min per point. The PET data were reconstructed by filtered back projection (FBP), the ordered subset expectation maximization (OSEM), OSEM + PSF, and OSEM + PSF + time-of-flight (TOF). Full width at half maximum (FWHM) was determined according to the NEMA method, and total counts in regions of interest (ROI) for each reconstruction were quantified.

Results

The radial FWHM of FBP and OSEM increased towards the peripheral FOV, whereas PSF-based reconstruction recovered the FWHM at all points in the FOV of both scanners. The radial FWHM for PSF was 30–50 % lower than that of OSEM at the center of the FOV. The accuracy of PSF correction was independent of detector design. Quantitative values were stable across the FOV in all reconstruction methods. The effect of TOF on spatial resolution and quantitation accuracy was less noticeable.

Conclusions

The traceable ²²Na point-like source allowed the evaluation of spatial resolution and quantitative accuracy across the FOV using different reconstruction methods and scanners. PSF-based reconstruction reduces dependence of the spatial resolution on the position. The quantitative accuracy over the entire FOV of the PET system is good, regardless of the reconstruction methods, although it depends slightly on the position.

[Accuracy of Resolution Recovery in PSF-based Fully-3D PET Image Reconstruction: Simulation and Phantom Study in Multicenter Trial]

PURPOSE

Recently, the quality of positron emission tomography (PET) images has rapidly improved using resolution recovery algorithm with point spread function (PSF). The aim of this study was to investigate the accuracy of the resolution recovery algorithm using three different PET systems.

METHODS

Three PET scanner models, the GE Discovery 600 M (D600M), SIEMENS Biograph mCT (mCT), and SHIMADZU SET-3000GCT/X (3000GCT) were used in this study. The radial dependences of spatial resolution (full width at half maximum: FWHM) were obtained by point source measurements (0.9 mmφ). All PET data were acquired in three-dimensional (3D) mode and reconstructed using the filtered back projection (FBP) , 3D-ordered subsets expectation maximization (3D-OSEM or dynamic row-action maximum likelihood algorithm) , and 3D-OSEM+PSF (PSF) algorithms. Two indicators, aspect ratio (ASR) and resolution recovery ratio (RRR), were calculated from measured FWHMs and compared among the three PET scanners.

RESULTS

In D600 and 3000GCT, distortions of the radial direction were slightly increased at circumference of field of view (FOV). On the other hand, random distortions were occurred in both radial and tangential direction in mCT. ASRs calculated from 3D-OSEM images at circumference of FOV were 2.06, 1.22, and 2.04 on D600M, mCT, and 3000GCT, respectively. ASR improved with PSF in all PET scanners. On the other hand, RRR with PSF were calculated 57.6%, 61.4%, and 31.6%, respectively.

CONCLUSION

Our results suggest that the spatial resolutions of PET images could be improved with PSF algorithm in all PET systems; however, effect of PSF was different depending on PET systems. Furthermore, PSF algorithm could not completely improve spatial resolutions in circumference of FOV.