|
|
Llosa, G., Barrillon, P., Barrio, J., Bisogni, M. G., Cabello, J., Del Guerra, A., et al. (2013). High performance detector head for PET and PET/MR with continuous crystals and SiPMs. Nucl. Instrum. Methods Phys. Res. A, 702, 3–5.
Abstract: A high resolution PET detector head for small animal PET applications has been developed. The detector is composed of a 12 mm x 12 mm continuous LYSO crystal coupled to a 64-channel monolithic SiPM matrix from FBK-irst. Crystal thicknesses of 5 mm and 10 mm have been tested, both yielding an intrinsic spatial resolution around 0.7 mm FWHM with a position determination algorithm that can also provide depth-of-interaction information. The detectors have been tested in a rotating system that makes it possible to acquire tomographic data and reconstruct images of Na-22 sources. An image reconstruction method specifically adapted for continuous crystals has been employed. The Full Width at Half Maximum measured from a point source reconstructed with ML-EM was 0.7 mm with the 5 mm crystal and 0.8 mm with the 10 mm crystal.
|
|
|
|
Gillam, J. E., Solevi, P., Oliver, J. F., & Rafecas, M. (2013). Simulated one-pass list-mode: an approach to on-the-fly system matrix calculation. Phys. Med. Biol., 58(7), 2377–2394.
Abstract: In the development of prototype systems for positron emission tomography a valid and robust image reconstruction algorithm is required. However, prototypes often employ novel detector and system geometries which may change rapidly under optimization. In addition, developing systems generally produce highly granular, or possibly continuous detection domains which require some level of on-the-fly calculation for retention of measurement precision. In this investigation a new method of on-the-fly system matrix calculation is proposed that provides advantages in application to such list-mode systems in terms of flexibility in system modeling. The new method is easily adaptable to complicated system geometries and available computational resources. Detection uncertainty models are used as random number generators to produce ensembles of possible photon trajectories at image reconstruction time for each datum in the measurement list. However, the result of this approach is that the system matrix elements change at each iteration in a non-repetitive manner. The resulting algorithm is considered the simulation of a one-pass list (SOPL) which is generated and the list traversed during image reconstruction. SOPL alters the system matrix in use at each iteration and so behavior within the maximum likelihood-expectation maximization algorithm was investigated. A two-pixel system and a small two dimensional imaging model are used to illustrate the process and quantify aspects of the algorithm. The two-dimensional imaging system showed that, while incurring a penalty in image resolution, in comparison to a non-random equal-computation counterpart, SOPL provides much enhanced noise properties. In addition, enhancement in system matrix quality is straightforward (by increasing the number of samples in the ensemble) so that the resolution penalty can be recovered when desired while retaining improvement in noise properties. Finally the approach is tested and validated against a standard (highly accurate) system matrix using experimental data from a prototype system-the AX-PET.
|
|
|
|
Solevi, P. et al, Oliver, J. F., Gillam, J. E., & Rafecas, M. (2013). A Monte-Carlo based model of the AX-PET demonstrator and its experimental validation. Phys. Med. Biol., 58(16), 5495–5510.
Abstract: AX-PET is a novel PET detector based on axially oriented crystals and orthogonal wavelength shifter (WLS) strips, both individually read out by silicon photo-multipliers. Its design decouples sensitivity and spatial resolution, by reducing the parallax error due to the layered arrangement of the crystals. Additionally the granularity of AX-PET enhances the capability to track photons within the detector yielding a large fraction of inter-crystal scatter events. These events, if properly processed, can be included in the reconstruction stage further increasing the sensitivity. Its unique features require dedicated Monte-Carlo simulations, enabling the development of the device, interpreting data and allowing the development of reconstruction codes. At the same time the non-conventional design of AX-PET poses several challenges to the simulation and modeling tasks, mostly related to the light transport and distribution within the crystals and WLS strips, as well as the electronics readout. In this work we present a hybrid simulation tool based on an analytical model and a Monte-Carlo based description of the AX-PET demonstrator. It was extensively validated against experimental data, providing excellent agreement.
|
|
|
|
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.
|
|
|
|
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.
|
|
