An automatic algorithm to exploit the symmetries of the system response matrix in positron emission tomography iterative reconstruction.

Positron emission tomography (PET) iterative 3D reconstruction is a very computational demanding task. One of the main issues of the iterative reconstruction concerns the management of the system response matrix (SRM). The SRM models the relationship between the projection and the voxel space and its memory footprint can easily exceed hundreds of GB. Moreover, in order to make the reconstruction fast enough not to hinder its practical application, the SRM must be stored in the random access memory of the workstation used for the reconstruction. This issue is normally solved by implementing efficient storage schemes and by reducing the number of redundant patterns in the SRM through symmetries. However, finding a sufficient number of symmetries is often non-trivial and is typically performed using dedicated solutions that cannot be exported to different detectors and geometries. In this paper, an automatic approach to reduce the memory footprint of a pre-computed SRM is described. The proposed approach was named symmetry search algorithm (SSA) and consists in an algorithm that searches for some of the redundant patterns present in the SRM, leading to its lossy compression. This approach was built to detect translations, reflections and coordinates swap in voxel space. Therefore, it is particularly well suited for those scanners where some of the rotational symmetries are broken, e.g. small animal scanner where the modules are arranged in a polygonal ring made of few elements, and dual head planar PET systems. In order to validate this approach, the SSA is applied to the SRM of a preclinical scanner (the IRIS PET/CT). The data acquired by the scanner were reconstructed with a dedicated maximum likelihood estimation maximization algorithm with both the uncompressed and the compressed SRMs. The results achieved show that the information lost due to the SSA compression is negligible. Compression factors up to 52 when using the SSA together with manually inserted symmetries and up to 204 when using the SSA alone, can be obtained for the IRIS SRM. These results come without significant differences in the values and in the main quality metrics of the reconstructed images, i.e. spatial resolution and noise. Although the compression factors depend on the system considered, the SSA is applicable to any SRM and therefore it can be considered a general tool to reduce the footprint of a pre-computed SRM.

Cherenkov luminescence measurements with digital silicon photomultipliers: a feasibility study.

BACKGROUND: A feasibility study was done to assess the capability of digital silicon photomultipliers to measure the Cherenkov luminescence emitted by a ß source. Cherenkov luminescence imaging (CLI) is possible with a charge coupled device (CCD) based technology, but a stand-alone technique for quantitative activity measurements based on Cherenkov luminescence has not yet been developed. Silicon photomultipliers (SiPMs) are photon counting devices with a fast impulse response and can potentially be used to quantify ß-emitting radiotracer distributions by CLI. METHODS: In this study, a Philips digital photon counting (PDPC) silicon photomultiplier detector was evaluated for measuring Cherenkov luminescence. The PDPC detector is a matrix of avalanche photodiodes, which were read one at a time in a dark count map (DCM) measurement mode (much like a CCD). This reduces the device active area but allows the information from a single avalanche photodiode to be preserved, which is not possible with analog SiPMs. An algorithm to reject the noisiest photodiodes and to correct the measured count rate for the dark current was developed. RESULTS: The results show that, in DCM mode and at (10-13) °C, the PDPC has a dynamic response to different levels of Cherenkov luminescence emitted by a ß source and transmitted through an opaque medium. This suggests the potential for this approach to provide quantitative activity measurements. Interestingly, the potential use of the PDPC in DCM mode for direct imaging of Cherenkov luminescence, as a opposed to a scalar measurement device, was also apparent. CONCLUSIONS: We showed that a PDPC tile in DCM mode is able to detect and image a ß source through its Cherenkov radiation emission. The detector's dynamic response to different levels of radiation suggests its potential quantitative capabilities, and the DCM mode allows imaging with a better spatial resolution than the conventional event-triggered mode. Finally, the same acquisition procedure and data processing could be employed also for other low light levels applications, such as bioluminescence.

Quantitative micro-CT based coronary artery profiling using interactive local thresholding and cylindrical coordinates.

BACKGROUND: Micro-CT is an established imaging technique for high-resolution non-destructive assessment of vascular samples, which is gaining growing interest for investigations of atherosclerotic arteries both in humans and in animal models. However, there is still a lack in the definition of micro-CT image metrics suitable for comprehensive evaluation and quantification of features of interest in the field of experimental atherosclerosis (ATS). OBJECTIVE: A novel approach to micro-CT image processing for profiling of coronary ATS is described, providing comprehensive visualization and quantification of contrast agent-free 3D high-resolution reconstruction of full-length artery walls. METHODS: Accelerated coronary ATS has been induced by high fat cholesterol-enriched diet in swine and left coronary artery (LCA) harvested en bloc for micro-CT scanning and histologic processing. A cylindrical coordinate system has been defined on the image space after curved multiplanar reformation of the coronary vessel for the comprehensive visualization of the main vessel features such as wall thickening and calcium content. A novel semi-automatic segmentation procedure based on 2D histograms has been implemented and the quantitative results validated by histology. RESULTS: The potentiality of attenuation-based micro-CT at low kV to reliably separate arterial wall layers from adjacent tissue as well as identify wall and plaque contours and major tissue components has been validated by histology. Morphometric indexes from histological data corresponding to several micro-CT slices have been derived (double observer evaluation at different coronary ATS stages) and highly significant correlations (R2 > 0.90) evidenced. Semi-automatic morphometry has been validated by double observer manual morphometry of micro-CT slices and highly significant correlations were found (R2 > 0.92). CONCLUSION: The micro-CT methodology described represents a handy and reliable tool for quantitative high resolution and contrast agent free full length coronary wall profiling, able to assist atherosclerotic vessels morphometry in a preclinical experimental model of coronary ATS and providing a link between in vivo imaging and histology.

Multi-ray-based system matrix generation for 3D PET reconstruction.

Iterative image reconstruction algorithms for positron emission tomography (PET) require a sophisticated system matrix (model) of the scanner. Our aim is to set up such a model offline for the YAP-(S)PET II small animal imaging tomograph in order to use it subsequently with standard ML-EM (maximum-likelihood expectation maximization) and OSEM (ordered subset expectation maximization) for fully three-dimensional image reconstruction. In general, the system model can be obtained analytically, via measurements or via Monte Carlo simulations. In this paper, we present the multi-ray method, which can be considered as a hybrid method to set up the system model offline. It incorporates accurate analytical (geometric) considerations as well as crystal depth and crystal scatter effects. At the same time, it has the potential to model seamlessly other physical aspects such as the positron range. The proposed method is based on multiple rays which are traced from/to the detector crystals through the image volume. Such a ray-tracing approach itself is not new; however, we derive a novel mathematical formulation of the approach and investigate the positioning of the integration (ray-end) points. First, we study single system matrix entries and show that the positioning and weighting of the ray-end points according to Gaussian integration give better results compared to equally spaced integration points (trapezoidal integration), especially if only a small number of integration points (rays) are used. Additionally, we show that, for a given variance of the single matrix entries, the number of rays (events) required to calculate the whole matrix is a factor of 20 larger when using a pure Monte-Carlo-based method. Finally, we analyse the quality of the model by reconstructing phantom data from the YAP-(S)PET II scanner.

Development of the YAP-PEM scanner for breast cancer imaging.

A prototype for positron emission mammography is under development at the Department of Physics of Pisa University. The device will be composed of two opposing detectors (parallel plane geometry). The active part of each detector head consists of a matrix of 900 YAP: Ce pixel scintillators, with a 2x2 mm(2) pitch and a 30 mm thickness. The read out is performed by an array of nine metal channel dynode PSPMTS (mod. R8520-00-C12) from Hamamatsu. In the previous version of the head, the PSPMTS were independently read out. For the clinical implementation of the prototype we have designed a simplified circuitry for the readout of the nine tubes based on a multiplexed resistive divider, reducing the number of channels from 36 to 4. A simulation study for an optimised amplifier has been carried out. The housing for each of the two yap-pem detectors has been fully engineered and is in the assembly stage.