|
|
Oliver, J. F., & Rafecas, M. (2010). Improving the singles rate method for modeling accidental coincidences in high-resolution PET. Phys. Med. Biol., 55(22), 6951–6971.
Abstract: Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.
|
|
|
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
|
|
|
Blume, M., Martinez-Moller, A., Keil, A., Navab, N., & Rafecas, M. (2010). Joint Reconstruction of Image and Motion in Gated Positron Emission Tomography. IEEE Trans. Med. Imaging, 29(11), 1892–1906.
Abstract: We present a novel intrinsic method for joint reconstruction of both image and motion in positron emission tomography (PET). Intrinsic motion compensation methods exclusively work on the measured data, without any external motion measurements. Most of these methods separate image from motion estimation: They use deformable image registration/optical flow techniques in order to estimate the motion from individually reconstructed gates. Then, the image is estimated based on this motion information. With these methods, a main problem lies in the motion estimation step, which is based on the noisy gated frames. The more noise is present, the more inaccurate the image registration becomes. As we show both visually and quantitatively, joint reconstruction using a simple deformation field motion model can compete with state-of-the-art image registration methods which use robust multilevel B-spline motion models.
|
|
|
|
Llosa, G., Barrio, J., Cabello, J., Crespo, A., Lacasta, C., Rafecas, M., et al. (2012). Detector characterization and first coincidence tests of a Compton telescope based on LaBr3 crystals and SiPMs. Nucl. Instrum. Methods Phys. Res. A, 695, 105–108.
Abstract: A Compton telescope for dose monitoring in hadron therapy consisting of several layers of continuous LaBr3 crystals coupled to silicon photomultiplier (SiPM) arrays is under development within the ENVISION project. In order to test the possibility of employing such detectors for the telescope, a detector head consisting of a continuous 16 mm x 18 mm x 5 mm LaBr3 crystal coupled to a SiPM array has been assembled and characterized, employing the SPIROC1 ASIC as readout electronics. The best energy resolution obtained at 511 key is 6.5% FWHM and the timing resolution is 3.1 ns FWHM. A position determination method for continuous crystals is being tested, with promising results. In addition, the detector has been operated in time coincidence with a second detector layer, to determine the coincidence capabilities of the system. The first tests are satisfactory, and encourage the development of larger detectors that will compose the telescope prototype.
|
|
|
|
Beltrame, P. et al, Oliver, J. F., Rafecas, M., & Solevi, P. (2011). The AX-PET demonstrator-Design, construction and characterization. Nucl. Instrum. Methods Phys. Res. A, 654(1), 546–559.
Abstract: Axial PET is a novel geometrical concept for Positron Emission Tomography (PET), based on layers of long scintillating crystals axially aligned with the bore axis. The axial coordinate is obtained from arrays of wavelength shifting (WLS) plastic strips placed orthogonally to the crystals. This article describes the design, construction and performance evaluation of a demonstrator set-up which consists of two identical detector modules, used in coincidence. Each module comprises 48 LYSO crystals of 100 mm length and 156 WLS strips. Crystals and strips are readout by Geiger-mode Avalanche Photo Diodes (G-APDs). The signals from the two modules are processed by fully analog front-end electronics and recorded in coincidence by a VME-based data acquisition system. Measurements with point-like (22)Na sources, with the modules used both individually and in coincidence mode, allowed for a complete performance evaluation up to the focal plane reconstruction of point sources. The results obtained are in good agreement with expectations and proved the set-up to be ready for the next evaluation phase with PET phantoms filled with radiotracers.
|
|
