|
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.
|
|
|
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.
|
|
|
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.
|
|
|
Blume, M., Navab, N., & Rafecas, M. (2012). Joint image and motion reconstruction for PET using a B-spline motion model. Phys. Med. Biol., 57(24), 22pp.
Abstract: We present a novel joint image and motion reconstruction method for PET. The method is based on gated data and reconstructs an image together with amotion function. The motion function can be used to transform the reconstructed image to any of the input gates. All available events (from all gates) are used in the reconstruction. The presented method uses a B-spline motion model, together with a novel motion regularization procedure that does not need a regularization parameter (which is usually extremely difficult to adjust). Several image and motion grid levels are used in order to reduce the reconstruction time. In a simulation study, the presented method is compared to a recently proposed joint reconstruction method. While the presented method provides comparable reconstruction quality, it is much easier to use since no regularization parameter has to be chosen. Furthermore, since the B-spline discretization of the motion function depends on fewer parameters than a displacement field, the presented method is considerably faster and consumes less memory than its counterpart. The method is also applied to clinical data, for which a novel purely data-driven gating approach is presented.
|
|
|
Bolle, E., Casella, C., Chesi, E., De Leo, R., Dissertori, G., Fanti, V., et al. (2012). AX-PET: A novel PET concept with G-APD readout. Nucl. Instrum. Methods Phys. Res. A, 695, 129–134.
Abstract: The AX-PET collaboration has developed a novel concept for high resolution PET imaging to overcome some of the performance limitations of classical PET cameras, in particular the compromise between spatial resolution and sensitivity introduced by the parallax error. The detector consists of an arrangement of long LYSO scintillating crystals axially oriented around the field of view together with arrays of wave length shifter strips orthogonal to the crystals. This matrix allows a precise 3D measurement of the photon interaction point. This is valid both for photoelectric absorption at 511 key and for Compton scattering down to deposited energies of about 100 keV. Crystals and WLS strips are individually read out using Geiger-mode Avalanche Photo Diodes (G-APDs). The sensitivity of such a detector can be adjusted by changing the number of layers and the resolution is defined by the crystal and strip dimensions. Two AX-PET modules were built and fully characterized in dedicated test set-ups at CERN, with point-like Na-22 sources. Their performance in terms of energy (Renew approximate to 11.8% (FWMH) at 511 key) and spatial resolution was assessed (sigma(axial) approximate to 0.65 mm), both individually and for the two modules in coincidence. Test campaigns at ETH Zurich and at the company AAA allowed the tomographic reconstructions of more complex phantoms validating the 3D reconstruction algorithms. The concept of the AX-PET modules will be presented together with some characterization results. We describe a count rate model which allows to optimize the planing of the tomographic scans.
|
|
|
Brzezinski, K., Oliver, J. F., Gillam, J., & Rafecas, M. (2014). Study of a high-resolution PET system using a Silicon detector probe. Phys. Med. Biol., 59(20), 6117–6140.
Abstract: A high-resolution silicon detector probe, in coincidence with a conventional PET scanner, is expected to provide images of higher quality than those achievable using the scanner alone. Spatial resolution should improve due to the finer pixelization of the probe detector, while increased sensitivity in the probe vicinity is expected to decrease noise. A PET-probe prototype is being developed utilizing this principle. The system includes a probe consisting of ten layers of silicon detectors, each a 80 x 52 array of 1 x 1 x 1 mm(3) pixels, to be operated in coincidence with a modern clinical PET scanner. Detailed simulation studies of this system have been performed to assess the effect of the additional probe information on the quality of the reconstructed images. A grid of point sources was simulated to study the contribution of the probe to the system resolution at different locations over the field of view (FOV). A resolution phantom was used to demonstrate the effect on image resolution for two probe positions. A homogeneous source distribution with hot and cold regions was used to demonstrate that the localized improvement in resolution does not come at the expense of the overall quality of the image. Since the improvement is constrained to an area close to the probe, breast imaging is proposed as a potential application for the novel geometry. In this sense, a simplified breast phantom, adjacent to heart and torso compartments, was simulated and the effect of the probe on lesion detectability, through measurements of the local contrast recovery coefficient-to-noise ratio (CNR), was observed. The list-mode ML-EM algorithm was used for image reconstruction in all cases. As expected, the point spread function of the PET-probe system was found to be non-isotropic and vary with position, offering improvement in specific regions. Increase in resolution, of factors of up to 2, was observed in the region close to the probe. Images of the resolution phantom showed visible improvement in resolution when including the probe in the simulations. The image quality study demonstrated that contrast and spill-over ratio in other areas of the FOV were not sacrificed for this enhancement. The CNR study performed on the breast phantom indicates increased lesion detectability provided by the probe.
|
|
|
Brzezinski, K., Oliver, J. F., Gillam, J., Rafecas, M., Studen, A., Grkovski, M., et al. (2016). Experimental evaluation of the resolution improvement provided by a silicon PET probe. J. Instrum., 11, P09016–13pp.
Abstract: A high-resolution PET system, which incorporates a silicon detector probe into a conventional PET scanner, has been proposed to obtain increased image quality in a limited region of interest. Detailed simulation studies have previously shown that the additional probe information improves the spatial resolution of the reconstructed image and increases lesion detectability, with no cost to other image quality measures. The current study expands on the previous work by using a laboratory prototype of the silicon PET-probe system to examine the resolution improvement in an experimental setting. Two different versions of the probe prototype were assessed, both consisting of a back-to-back pair of 1-mm thick silicon pad detectors, one arranged in 32 x 16 arrays of 1.4mm x 1.4mm pixels and the other in 40 x 26 arrays of 1.0mm x 1.0mm pixels. Each detector was read out by a set of VATAGP7 ASICs and a custom-designed data acquisition board which allowed trigger and data interfacing with the PET scanner, itself consisting of BGO block detectors segmented into 8 x 6 arrays of 6mm x 12mm x 30mm crystals. Limited-angle probe data was acquired from a group of Na-22 point-like sources in order to observe the maximum resolution achievable using the probe system. Data from a Derenzo-like resolution phantom was acquired, then scaled to obtain similar statistical quality as that of previous simulation studies. In this case, images were reconstructed using measurements of the PET ring alone and with the inclusion of the probe data. Images of the Na-22 source demonstrated a resolution of 1.5mm FWHM in the probe data, the PET ring resolution being approximately 6 mm. Profiles taken through the image of the Derenzo-like phantom showed a clear increase in spatial resolution. Improvements in peak-to-valley ratios of 50% and 38%, in the 4.8mm and 4.0mm phantom features respectively, were observed, while previously unresolvable 3.2mm features were brought to light by the addition of the probe. These results support the possibility of improving the image resolution of a clinical PET scanner using the silicon PET-probe.
|
|
|
Cabello, J., & Rafecas, M. (2012). Comparison of basis functions for 3D PET reconstruction using a Monte Carlo system matrix. Phys. Med. Biol., 57(7), 1759–1777.
Abstract: In emission tomography, iterative statistical methods are accepted as the reconstruction algorithms that achieve the best image quality. The accuracy of these methods relies partly on the quality of the system response matrix (SRM) that characterizes the scanner. The more physical phenomena included in the SRM, the higher the SRM quality, and therefore higher image quality is obtained from the reconstruction process. High-resolution small animal scanners contain as many as 10(3)-10(4) small crystal pairs, while the field of view (FOV) is divided into hundreds of thousands of small voxels. These two characteristics have a significant impact on the number of elements to be calculated in the SRM. Monte Carlo (MC) methods have gained popularity as a way of calculating the SRM, due to the increased accuracy achievable, at the cost of introducing some statistical noise and long simulation times. In the work presented here the SRM is calculated using MC methods exploiting the cylindrical symmetries of the scanner, significantly reducing the simulation time necessary to calculate a high statistical quality SRM and the storage space necessary. The use of cylindrical symmetries makes polar voxels a convenient basis function. Alternatively, spherically symmetric basis functions result in improved noise properties compared to cubic and polar basis functions. The quality of reconstructed images using polar voxels, spherically symmetric basis functions on a polar grid, cubic voxels and post-reconstruction filtered polar and cubic voxels is compared from a noise and spatial resolution perspective. This study demonstrates that polar voxels perform as well as cubic voxels, reducing the simulation time necessary to calculate the SRM and the disk space necessary to store it. Results showed that spherically symmetric functions outperform polar and cubic basis functions in terms of noise properties, at the cost of slightly degraded spatial resolution, larger SRM file size and longer reconstruction times. However, we demonstrate that post-reconstruction smoothing, usually applied in emission imaging to reduce the level of noise, can produce a spatial resolution degradation of similar to 50%, while spherically symmetric basis functions produce a degradation of only similar to 6%, compared to polar and cubic voxels, at the same noise level. Therefore, the image quality trade-off obtained with blobs is higher than that obtained with cubic or polar voxels.
|
|
|
Cabello, J., Torres-Espallardo, I., Gillam, J. E., & Rafecas, M. (2013). PET Reconstruction From Truncated Projections Using Total-Variation Regularization for Hadron Therapy Monitoring. IEEE Trans. Nucl. Sci., 60(5), 3364–3372.
Abstract: Hadron therapy exploits the properties of ion beams to treat tumors by maximizing the dose released to the target and sparing healthy tissue. With hadron beams, the dose distribution shows a relatively low entrance dose which rises sharply at the end of the range, providing the characteristic Bragg peak that drops quickly thereafter. It is of critical importance in order not to damage surrounding healthy tissues and/or avoid targeting underdosage to know where the delivered dose profile ends-the location of the Bragg peak. During hadron therapy, short-lived beta(+)-emitters are produced along the beam path, their distribution being correlated with the delivered dose. Following positron annihilation, two photons are emitted, which can be detected using a positron emission tomography (PET) scanner. The low yield of emitters, their short half-life, and the wash out from the target region make the use of PET, even only a few minutes after hadron irradiation, a challenging application. In-beam PET represents a potential candidate to estimate the distribution of beta(+)-emitters during or immediately after irradiation, at the cost of truncation effects and degraded image quality due to the partial rings required of the PET scanner. Time-of-flight (ToF) information can potentially be used to compensate for truncation effects and to enhance image contrast. However, the highly demanding timing performance required in ToF-PET makes this option costly. Alternatively, the use of maximum-a-posteriori-expectation-maximization (MAP-EM), including total variation (TV) in the cost function, produces images with low noise, while preserving spatial resolution. In this paper, we compare data reconstructed with maximum-likelihood-expectation-maximization (ML-EM) and MAP-EM using TV as prior, and the impact of including ToF information, from data acquired with a complete and a partial-ring PET scanner, of simulated hadron beams interacting with a polymethyl methacrylate (PMMA) target. The results show that MAP-EM, in the absence of ToF information, produces lower noise images and more similar data compared to the simulated beta(+) distributions than ML-EM with ToF information in the order of 200-600 ps. The investigation is extended to the combination of MAP-EM and ToF information to study the limit of performance using both approaches.
|
|
|
Gillam, J. E., Solevi, P., Oliver, J. F., Casella, C., Heller, M., Joram, C., et al. (2014). Sensitivity recovery for the AX-PET prototype using inter-crystal scattering events. Phys. Med. Biol., 59(15), 4065–4083.
Abstract: The development of novel detection devices and systems such as the AX-positron emission tomography (PET) demonstrator often introduce or increase the measurement of atypical coincidence events such as inter-crystal scattering (ICS). In more standard systems, ICS events often go undetected and the small measured fraction may be ignored. As the measured quantity of such events in the data increases, so too does the importance of considering them during image reconstruction. Generally, treatment of ICS events will attempt to determine which of the possible candidate lines of response (LoRs) correctly determine the annihilation photon trajectory. However, methods of assessment often have low success rates or are computationally demanding. In this investigation alternative approaches are considered. Experimental data was taken using the AX-PET prototype and a NEMA phantom. Three methods of ICS treatment were assessed-each of which considered all possible candidate LoRs during image reconstruction. Maximum likelihood expectation maximization was used in conjunction with both standard (line-like) and novel (V-like in this investigation) detection responses modeled within the system matrix. The investigation assumed that no information other than interaction locations was available to distinguish between candidates, yet the methods assessed all provided means by which such information could be included. In all cases it was shown that the signal to noise ratio is increased using ICS events. However, only one method, which used full modeling of the ICS response in the system matrix-the V-like model-provided enhancement in all figures of merit assessed in this investigation. Finally, the optimal method of ICS incorporation was demonstrated using data from two small animals measured using the AX-PET demonstrator.
|
|