Our aim is to investigate the impact of respiratory motion on tumor quantification and delineation in static PET/CT imaging using a population of patient respiratory traces. A total of 1295 respiratory traces acquired during whole body PET/CT imaging were classified into three types according to the qualitative shape of their signal histograms. Each trace was scaled to three diaphragm motion amplitudes (6 mm, 11 mm and 16 mm) to drive a whole body PET/CT computer simulation that was validated with a physical phantom experiment. Three lung lesions and one liver lesion were simulated with diameters of 1 cm and 2 cm. PET data were reconstructed using the OS-EM algorithm with attenuation correction using CT images at the end-expiration phase and respiratory-averaged CT. The errors of the lesion maximum standardized uptake values (SUV max ) and lesion volumes between motion-free and motion-blurred PET/CT images were measured and analyzed. For respiration with 11 mm diaphragm motion and larger quiescent period fraction, respiratory motion can cause a mean lesion SUV max underestimation of 28% and a mean lesion volume overestimation of 130% in PET/CT images with 1 cm lesions. The errors of lesion SUV max and volume are larger for patient traces with larger motion amplitudes. Smaller lesions are more sensitive to respiratory motion than larger lesions for the same motion amplitude. Patient respiratory traces with relatively larger quiescent period fraction yield results less subject to respiratory motion than traces with long-term amplitude variability. Mismatched attenuation correction due to respiratory motion can cause SUV max overestimation for lesions in the lower lung region close to the liver dome. Using respiratoryaveraged CT for attenuation correction yields smaller mismatch errors than those using endexpiration CT. Respiratory motion can have a significant impact on static oncological PET/CT imaging where SUV and/or volume measurements are important. The impact is highly dependent upon motion amplitude, lesion location and size, attenuation map and respiratory pattern. To overcome the motion effect, motion compensation techniques may be necessary in clinical practice to improve the tumor quantification for determining the response to therapy or for radiation treatment planning.

Accurate system modeling in tomographic image reconstruction has been shown to reduce the spatial variance of resolution and improve quantitative accuracy. System modeling can be improved through analytic calculations, Monte Carlo simulations, and physical measurements. The purpose of this work is to improve clinical fully-3-D reconstruction without substantially increasing computation time. We present a practical method for measuring the detector blurring component of a whole-body positron emission tomography (PET) system to form an approximate system model for use with fully-3-D reconstruction. We employ Monte Carlo simulations to show that a non-collimated point source is acceptable for modeling the radial blurring present in a PET tomograph and we justify the use of a Na22 point source for collecting these measurements. We measure the system response on a whole-body scanner, simplify it to a 2-D function, and incorporate a parameterized version of this response into a modified fully-3-D OSEM algorithm. Empirical testing of the signal versus noise benefits reveal roughly a 15% improvement in spatial resolution and 10% improvement in contrast at matched image noise levels. Convergence analysis demonstrates improved resolution and contrast versus noise properties can be achieved with the proposed method with similar computation time as the conventional approach. Comparison of the measured spatially variant and invariant reconstruction revealed similar performance with conventional image metrics. Edge artifacts, which are a common artifact of resolution-modeled reconstruction methods, were less apparent in the spatially variant method than in the invariant method. With the proposed and other resolution-modeled reconstruction methods, edge artifacts need to be studied in more detail to determine the optimal tradeoff of resolution/contrast enhancement and edge fidelity.

The addition of accurate system modeling in PET image reconstruction results in images with distinct noise texture and characteristics. In particular, the incorporation of point spread functions (PSF) into the system model has been shown to visually reduce image noise, but the noise properties have not been thoroughly studied. This work offers a systematic evaluation of noise and signal properties in different combinations of reconstruction methods and parameters. We evaluate two fully 3D PET reconstruction algorithms: (1) OSEM with exact scanner line of response modeled (OSEM+LOR), (2) OSEM with line of response and a measured point spread function incorporated (OSEM+LOR +PSF), in combination with the effects of four post-reconstruction filtering parameters and 1-10 iterations, representing a range of clinically acceptable settings. We used a modified NEMA image quality (IQ) phantom, which was filled with 68 Ge and consisted of six hot spheres of different sizes with a target/background ratio of 4:1. The phantom was scanned 50 times in 3D mode on a clinical system to provide independent noise realizations. Data were reconstructed with OSEM+LOR and OSEM+LOR+PSF using different reconstruction parameters, and our implementations of the algorithms match the vendor's product algorithms. With access to multiple realizations, background noise characteristics were quantified with four metrics. Image roughness and the standard deviation image measured the pixel-to-pixel variation; background variability and ensemble noise quantified the region-to-region variation. Image roughness is the image noise perceived when viewing an individual image. At matched iterations, the addition of PSF leads to images with less noise defined as image roughness (reduced by 35% for unfiltered data) and as the standard deviation image, while it has no effect on background variability or ensemble noise. In terms of signal to noise performance, PSF-based reconstruction has a 7% improvement in contrast recovery at matched ensemble noise levels and 20% improvement of quantitation SNR in unfiltered data. In addition, the relations between different metrics are studied. A linear correlation is observed between background variability and ensemble noise for all different combinations of reconstruction methods and parameters, suggesting that background variability is a reasonable surrogate for ensemble noise when multiple realizations of scans are not available.

Abstract-Appropriate application of spatially variant system models can correct for degraded resolution response and mispositioning errors. This paper explores the detector blurring component of the system model for a whole body positron emission tomography (PET) system and extends this factor into a more general system response function to account for other system effects including the influence of Fourier rebinning (FORE). We model the system response function as a three-dimensional

Objective To compare myocardial blood flow (MBF) and myocardial flow reserve (MFR) estimates from 82Rb PET data using ten software packages (SPs): Carimas, Corridor4DM, FlowQuant, HOQUTO, ImagenQ, MunichHeart, PMOD, QPET, syngo MBF, and UW-QPP. Background It is unknown how MBF and MFR values from existing SPs agree for 82Rb PET. Methods Rest and stress 82Rb PET scans of 48 patients with suspected or known coronary artery disease (CAD) were analyzed in 10 centers. Each center used one of the 10 SPs to analyze global and regional MBF using the different kinetic models implemented. Values were considered to agree if they simultaneously had an intraclass correlation coefficient (ICC) > 0.75 and a difference < 20% of the median across all programs. Results The most common model evaluated was the one-tissue compartment model (1TCM) by Lortie et al. (2007). MBF values from seven of the eight software packages implementing this model agreed best (Carimas, Corridor4DM, FlowQuant, PMOD, QPET, syngoMBF, and UW-QPP). Values from two other models (El Fakhri et al. in Corridor4DM and Alessio et al. in UW-QPP) also agreed well, with occasional differences. The MBF results from other models (Sitek et al. 1TCM in Corridor4DM, Katoh et al. 1TCM in HOQUTO, Herrero et al. 2TCM in PMOD, Yoshida et al. retention in ImagenQ, and Lautamäki et al. retention in MunichHeart) were less in agreement with Lortie 1TCM values. Conclusions SPs using the same kinetic model, as described in Lortie et al. (2007), provided consistent results in measuring global and regional MBF values, suggesting they may be used interchangeably to process data acquired with a common imaging protocol.