Search

Dual-panel PET system configuration can lead to spatially variable point-spread functions (PSF) of considerable deformations due to depth-of-interaction effects and limited angular coverage. If not modelled properly, these effects result in decreased and inconsistent recovery of lesion activity across the field-of-view (FOV), as well as mispositioning of lesions in the reconstructed image caused by strong PSF asymmetries. We implemented and evaluated models of such PSF deformations with spatially-variant image-based resolution modeling (IRM) within reconstruction (varRM) using the Direct Image REConstruction for Time-of-flight (DIRECT) method and within post-reconstruction deconvolution methods. In addition, DIRECT reconstruction was performed with a spatially-invariant IRM (invRM) and without resolution modeling (noRM) for comparison. The methods were evaluated using simulated data for a realistic breast model with a set of 5 mm lesions located throughout the FOV of a dual-panel Breast-PET scanner. We simulated high-count data to focus on the ability of each method to correctly recover the PSF deformations, and a clinically realistic count level to assess the impact of low count data on the quantitative performance of the evaluated techniques. Performance of the methods evaluated herein was assessed by comparing lesion activity recovery (%BIAS), consistency (%SD) across the FOV, overall error (%RMSE), and recovery of each lesion location. As expected, all techniques using IRM provide considerable improvement over the noRM reconstruction. For the high-count cases, the overall quantitative performance of all IRM techniques, whether within reconstruction or within post-reconstruction, is similar if the lesion location misplacements are ignored. However, invRM provides less consistent performance on activity across lesions and is not able to recover accurate lesion locations. For a clinically realistic count level, varRM reconstruction consistently outperforms all compared approaches, while the post-reconstruction IRM approaches exhibit higher %SD and %RMSE values due to being more affected by the data noise than the within-reconstruction IRM approaches. [ABSTRACT FROM AUTHOR]

This work explores the feasibility and impact of including both the motion correction and the image registration transformation parameters from positron emission tomography (PET) image space to magnetic resonance (MR), or stereotaxic, image space within the systemmatrix of PET image reconstruction. This approach is motivated by the fields of neuroscience and psychiatry, where PET is used to investigate differences in activation patterns between different groups of participants, requiring all images to be registered to a common spatial atlas. Currently, image registration is performed after image reconstruction which introduces interpolation effects into the final image. Furthermore, motion correction (also requiring registration) introduces a further level of interpolation, and the overall result of these operations can lead to resolution degradation and possibly artifacts. It is important to note that performing such operations on a post-reconstruction basis means, strictly speaking, that the final images are not ones which maximize the desired objective function (e.g. maximum likelihood (ML), or maximum a posteriori reconstruction (MAP)). To correctly seek parameter estimates in the desired spatial atlas which are in accordance with the chosen reconstruction objective function, it is necessary to include the transformation parameters for both motion correction and registration within the system modeling stage of image reconstruction. Such an approach not only respects the statistically chosen objective function (e.g. ML or MAP), but furthermore should serve to reduce the interpolation effects. To evaluate the proposed method, this work investigates registration (including motion correction) using 2D and 3D simulations based on the high resolution research tomograph (HRRT) PET scanner geometry, with and without resolutionmodeling, using the ML expectation maximization (MLEM) reconstruction algorithm. The quality of reconstructionwas assessed using biasvariance and root mean squared error analyses, comparing the proposed method to conventional post-reconstruction registration methods. An overall reduction in bias (for a cold region: from 41% down to 31% (2D) and 97% down to 65% (3D), and for a hot region: from 11% down to 8% (2D) and from 16% down to 14% (3D)) and in root mean squared error analyses (for a cold region: from 43% to 37% (2D) and from 97% to 65% (3D), and for a hot region: from 11% to 9% (2D) and from 16% down to 14% (3D)) in reconstructed regional mean activities (full regions of interest; all with statistical significance: p < 5 ✕ 10-10) is found when including the motion correction and registration in the system matrix of the MLEM reconstruction, with resolution modeling. However, this improvement in performance comes with an extra computational cost of about 40 min. In this context, this work constitutes an important step toward the goal of estimating parameters of interest directly from the raw Poisson-distributed PET data, and hence toward the complete elimination of post-processing steps. [ABSTRACT FROM AUTHOR]