Performance of three-photon PET imaging: Monte Carlo simulations

We have recently introduced the idea of making use of three-photon positron annihilations in positron emission tomography. In this paper the basic characteristics of the three-gamma imaging in PET are studied by means of Monte Carlo simulations and analytical computations. Two typical configurations of human and small animal scanners are considered. Three-photon imaging requires high energy resolution detectors. Parameters currently attainable by CdZnTe semiconductor detectors, the technology of choice for the future development of radiation imaging, are assumed. Spatial resolution is calculated as a function of detector energy resolution and size, position in the field of view, scanner size, and the energies of the three gamma annihilation photons. Possible ways to improve the spatial resolution obtained for nominal parameters: 1.5 cm and 3.2 mm FWHM for human and small animal scanners, respectively, are indicated. Counting rates of true and random three-photon events for typical human and small animal scanning configurations are assessed. A simple formula for minimum size of lesions detectable in the three-gamma based images is derived. Depending on the contrast and total number of registered counts, lesions of a few mm size for human and sub mm for small animal scanners can be detected

Feasibility study of the positronium imaging with the J-PET tomograph

A detection system of the conventional PET tomograph is set-up to record data from e+ e- annihilation into two photons with energy of 511 keV, and it gives information on the density distribution of a radiopharmaceutical in the body of the object. In this paper we explore the possibility of performing the three gamma photons imaging based on ortho-positronium annihilation, as well as the possibility of positronium mean lifetime imaging with the J-PET tomograph constructed from plastic scintillators. For this purposes simulations of the ortho-positronium formation and its annihilation into three photons were performed taking into account distributions of photons' momenta as predicted by the theory of quantum electrodynamics and the response of the J-PET tomograph. In order to test the proposed ortho-positronium lifetime image reconstruction method, we concentrate on the decay of the ortho-positronium into three photons and applications of radiopharmaceuticals labeled with isotopes emitting a prompt gamma quantum. The proposed method of imaging is based on the determination of hit-times and hit-positions of registered photons which enables the reconstruction of the time and position of the annihilation point as well as the lifetime of the ortho-positronium on an event-by-event basis. We have simulated the production of the positronium in a cylindrical phantom composed of a set of different materials in which the ortho-positronium lifetime varied from 2.0 ns to 3.0 ns, as expected for ortho-positronium created in the human body. The presented reconstruction method for total-body J-PET like detector allows to achieve a mean lifetime resolution of about 40 ps. Recent Positron Annihilation Lifetime Spectroscopy measurements of cancerous and healthy uterine tissues show that this sensitivity may allow to study the morphological changes in cell structures.Comment: accepted in PMB (http://iopscience.iop.org/article/10.1088/1361-6560/aafe20

Incorporating accurate statistical modeling in PET: reconstruction for whole-body imaging

Tese de doutoramento em Biofísica, apresentada à Universidade de Lisboa através da Faculdade de Ciências, 2007The thesis is devoted to image reconstruction in 3D whole-body PET imaging. OSEM ( Ordered Subsets Expectation maximization ) is a statistical algorithm that assumes Poisson data. However, corrections for physical effects (attenuation, scattered and random coincidences) and detector efficiency remove the Poisson characteristics of these data. The Fourier Rebinning (FORE), that combines 3D imaging with fast 2D reconstructions, requires corrected data. Thus, if it will be used or whenever data are corrected prior to OSEM, the need to restore the Poisson-like characteristics is present. Restoring Poisson-like data, i.e., making the variance equal to the mean, was achieved through the use of weighted OSEM algorithms. One of them is the NECOSEM, relying on the NEC weighting transformation. The distinctive feature of this algorithm is the NEC multiplicative factor, defined as the ratio between the mean and the variance. With real clinical data this is critical, since there is only one value collected for each bin the data value itself. For simulated data, if we keep track of the values for these two statistical moments, the exact values for the NEC weights can be calculated. We have compared the performance of five different weighted algorithms (FORE+AWOSEM, FORE+NECOSEM, ANWOSEM3D, SPOSEM3D and NECOSEM3D) on the basis of tumor detectablity. The comparison was done for simulated and clinical data. In the former case an analytical simulator was used. This is the ideal situation, since all the weighting factors can be exactly determined. For comparing the performance of the algorithms, we used the Non-Prewhitening Matched Filter (NPWMF) numerical observer. With some knowledge obtained from the simulation study we proceeded to the reconstruction of clinical data. In that case, it was necessary to devise a strategy for estimating the NEC weighting factors. The comparison between reconstructed images was done by a physician largely familiar with whole-body PET imaging

Imaging, Detection, and Identification Algorithms for Position-Sensitive Gamma-Ray Detectors.

Three-dimensional-position-sensitive semiconductors record both the locations and energies of gamma-ray interactions with high resolution, enabling spectroscopy and imaging of gamma-ray-emitting materials. Imaging enables the detection of point sources of gamma rays in an otherwise extended-source background, even when the background spectrum is unknown and may share the point source's spectrum. The generalized likelihood ratio test (GLRT) and source-intensity test (SIT) are applied to this situation to detect one-or-more unshielded point sources from a library of isotopes in a spectrally unknown or known background when the background intensity varies spatially by a factor of two or less. In addition to estimating the number of sources present, their activities, isotopes, and directions from the detector are estimated. Experimental and some simulated results are presented for a single detector and an 18-detector array of 2 cm by 2 cm by 1.5 cm CdZnTe crystals and compared with the performance of spectral-only detection when the background and source are assumed to be spectrally different. Furthermore, the expected detection performance of the 18-detector array system is investigated statistically using experimental data in the case where the background is distinct spectrally from the point source and the possible source location and isotopic identity are known. Including imaging gave at least 7% higher SNR compared to ignoring the image dimension. Also, imaging methods based on the maximum-likelihood expectation-maximization method are introduced to determine the spatial distribution of isotopes and to find the activity distributions within targets moving with known motion through a radioactive background. Software has also been developed to support the analysis of the data from 3D-position-sensitive spectroscopic systems, for a range of detector designs and applications. The software design and unique features that allow fast multidimensional data analysis are presented, along with parallel computing performance.Ph.D.Nuclear Engineering & Radiological SciencesUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/89797/1/cgwahl_1.pd