M. J. Rodriguez-Alvarez

Performance tests of two portable mini gamma cameras for medical applications.

We have developed two prototypes of portable gamma cameras for medical applications based on a previous prototype designed and tested by our group. These cameras use a CsI(Na) continuous scintillation crystal coupled to the new flat-panel-type multianode position-sensitive photomultiplier tube, H8500 from Hamamatsu Photonics. One of the prototypes, mainly intended for intrasurgical use, has a field of view of 44×44mm2, and weighs 1.2kg. Its intrinsic resolution is better than 1.5mm and its energy resolution is about 13% at 140keV. The second prototype, mainly intended for osteological, renal, mammary, and endocrine (thyroid, parathyroid, and suprarenal) scintigraphies, weighs a total of 2kg. Its average spatial resolution is 2mm; it has a field of view of 95×95mm2, with an energy resolution of about 15% at 140keV. The main advantages of these gamma camera prototypes with respect to those previously reported in the literature are high portability and low weight, with no significant loss of sensitivity and spatial resolution. All the electronic components are packed inside the mini gamma cameras, and no external electronic devices are required. The cameras are only connected through the universal serial bus port to a portable PC. In this paper, we present the design of the cameras and describe the procedures that have led us to choose their configuration together with the most important performance features of the cameras. For one of the prototypes, clinical tests on melanoma patients are presented and images are compared with those obtained with a conventional camera.

Small animal PET scanner based on monolithic LYSO crystals: Performance evaluation

PURPOSE The authors have developed a small animal Positron emission tomography (PET) scanner based on monolithic LYSO crystals coupled to multi-anode photomultiplier tubes (MA-PMTs). In this study, the authors report on the design, calibration procedure, and performance evaluation of a PET system that the authors have developed using this innovative nonpixelated detector design. METHODS The scanner is made up of eight compact modules forming an octagon with an axial field of view (FOV) of 40 mm and a transaxial FOV of 80 mm diameter. In order to fully determine its performance, a recently issued National Electrical Manufacturers Association (NEMA) NU-4 protocol, specifically developed for small animal PET scanners, has been followed. By measuring the width of light distribution collected in the MA-PMT the authors are able to determine depth of interaction (DOI), thus making the proper identification of lines of response (LORs) with large incidence angles possible. PET performances are compared with those obtained with currently commercially available small animal PET scanners. RESULTS At axial center when the point-like source is located at 5 mm from the radial center, the spatial resolution measured was 1.65, 1.80, and 1.86 mm full width at half maximum (FWHM) for radial, tangential, and axial image profiles, respectively. A system scatter fraction of 7.5% (mouse-like phantom) and 13% (rat-like phantom) was obtained, while the maximum noise equivalent count rate (NECR) was 16.9 kcps at 12.7 MBq (0.37 MBq/ml) for mouse-like phantom and 12.8 kcps at 12.4 MBq (0.042 MBq/ml) for rat-like phantom The peak absolute sensitivity in the center of the FOV is 2% for a 30% peak energy window. Several animal images are also presented. CONCLUSIONS The overall performance of our small animal PET is comparable to that obtained with much more complex crystal pixelated PET systems. Moreover, the new proposed PET produces high-quality images suitable for studies with small animals.

ALBIRA: A small animal PET/SPECT/CT imaging system

PURPOSE The authors have developed a trimodal PET∕SPECT∕CT scanner for small animal imaging. The gamma ray subsystems are based on monolithic crystals coupled to multianode photomultiplier tubes (MA-PMTs), while computed tomography (CT) comprises a commercially available microfocus x-ray tube and a CsI scintillator 2D pixelated flat panel x-ray detector. In this study the authors will report on the design and performance evaluation of the multimodal system. METHODS X-ray transmission measurements are performed based on cone-beam geometry. Individual projections were acquired by rotating the x-ray tube and the 2D flat panel detector, thus making possible a transaxial field of view (FOV) of roughly 80 mm in diameter and an axial FOV of 65 mm for the CT system. The single photon emission computed tomography (SPECT) component has a dual head detector geometry mounted on a rotating gantry. The distance between the SPECT module detectors can be varied in order to optimize specific user requirements, including variable FOV. The positron emission tomography (PET) system is made up of eight compact modules forming an octagon with an axial FOV of 40 mm and a transaxial FOV of 80 mm in diameter. The main CT image quality parameters (spatial resolution and uniformity) have been determined. In the case of the SPECT, the tomographic spatial resolution and system sensitivity have been evaluated with a (99m)Tc solution using single-pinhole and multi-pinhole collimators. PET and SPECT images were reconstructed using three-dimensional (3D) maximum likelihood and ordered subset expectation maximization (MLEM and OSEM) algorithms developed by the authors, whereas the CT images were obtained using a 3D based FBP algorithm. RESULTS CT spatial resolution was 85 μm while a uniformity of 2.7% was obtained for a water filled phantom at 45 kV. The SPECT spatial resolution was better than 0.8 mm measured with a Derenzo-like phantom for a FOV of 20 mm using a 1-mm pinhole aperture collimator. The full width at half-maximum PET radial spatial resolution at the center of the field of view was 1.55 mm. The SPECT system sensitivity for a FOV of 20 mm and 15% energy window was 700 cps∕MBq (7.8 × 10(-2)%) using a multi-pinhole equipped with five apertures 1 mm in diameter, whereas the PET absolute sensitivity was 2% for a 350-650 keV energy window and a 5 ns timing window. Several animal images are also presented. CONCLUSIONS The new small animal PET∕SPECT∕CT proposed here exhibits high performance, producing high-quality images suitable for studies with small animals. Monolithic design for PET and SPECT scintillator crystals reduces cost and complexity without significant performance degradation.

Calibration and Performance Tests of Detectors for Laser-Accelerated Protons

We present the calibration and performance tests carried out with two detectors for intense proton pulses accelerated by lasers. Most of the procedures were realized with proton beams of 0.46-5.60 MeV from a tandem accelerator. One approach made use of radiochromic films, for which we calibrated the relation between optical density and energy deposition over more than three orders of magnitude. The validity of these results and of our analysis algorithms has been confirmed by controlled irradiation of film stacks and reconstruction of the total beam charge for strongly non-uniform beam profiles. For the spectral analysis of protons from repeated laser shots, we have designed an online monitor based on a plastic scintillator. The resulting signal from a photomultiplier directly measured on a fast oscilloscope is especially useful for time-of-flight applications. Variable optical filters allow for suppression of saturation and an extension of the dynamic range. With pulsed proton beams we have tested the detector response to a wide range of beam intensities from single particles to 3 ×105 protons per 100 ns time interval.

Expectation maximization (EM) algorithms using polar symmetries for computed tomography (CT) image reconstruction

We suggest a symmetric-polar pixellation scheme which makes possible a reduction of the computational cost for expectation maximization (EM) iterative algorithms. The proposed symmetric-polar pixellation allows us to deal with 3D images as a whole problem without dividing the 3D problem into 2D slices approach. Performance evaluation of each approach in terms of stability and image quality is presented. Exhaustive comparisons between all approaches were conducted in a 2D based image reconstruction model. From these 2D approaches, that showing the best performances were finally implemented and evaluated in a 3D based image reconstruction model. Comparison to 3D images reconstructed with FBP is also presented. Although the algorithm is presented in the context of computed tomography (CT) image reconstruction, it can be applied to any other tomographic technique as well, due to the fact that the only requirement is a scanning geometry involving measurements of an object under different projection angles. Real data have been acquired with a small animal (CT) scanner to verify the proposed mathematical description of the CT system.

Characterization of protons accelerated from a 3 TW table-top laser system

This project has been funded by Centro para el Desarrollo Tecnologico Industrial (CDTI, Spain) within the INNPRONTA program, grant no. IPT-20111027, by EUROSTARS project E9113, and by the Spanish Ministry for Economy and Competitiveness within the Retos-Colaboracion 2015 initiative, ref. RTC-2015-3278-1.

Detailed requirements for a laser-based proton/ion accelerator for radioisotope production

Background and objectives: Laser-plasma acceleration of protons and ions is often considered a promising technique for compact applications of highly intense beams of multi-MeV particles. A remarkable example is the on-site production of short-lived radioisotopes for medical and preclinical interventions. We study quantitatively the activity of four important PET isotopes which may be obtained by irradiation of suitable target nuclei with laser-accelerated protons and deuterons. These simulations allow for confining the range of useful parameters of a laser-based production system. Methods: We choose a total of ten p- and d-induced reaction channels for the production of F-18, C-11, O-15, and N-13 from suitable target nuclei. We calculate the activity yield as a function of projectile energy starting from the corresponding, known cross sections. In order to simulate typical laser-plasma particle spectra we generate exponentially decaying distributions spread over a wide range, up to a maximum energy between 6 and 16 MeV. From the yield curves and the spectra we obtain the single-shot activation and the total activity after a realistic production time at 100 Hz pulse rate, taking into account saturation effects due to decay during irradiation. Results: We present numerical results for ten reaction channels and six realistic projectile spectra. With single laser shots, the highest activities are generated for O-15 (up to 11 (20) kBq for 16 MeV maximum p (d) energy). After prolonged irradiation at 100 Hz pulse rate, useful quantities of C-11 and O-15 may be obtained from spectra with 10 MeV maximum energies. The production of N-13 and F-18, to the contrary, requires higher energies and/or shot rates. Conclusions: 10 MeV particle energy and 100 Hz pulse rate are realistic benchmarks for a laser-based PET isotope production system. Experimental work to achieve these demanding objectives is in progress.

Sparse Givens resolution of large system of linear equations: Applications to image reconstruction

In medicine, computed tomographic images are reconstructed from a large number of measurements of X-ray transmission through the patient (projection data). The mathematical model used to describe a computed tomography device is a large system of linear equations of the form AX=B. In this paper we propose the QR decomposition as a direct method to solve the linear system. QR decomposition can be a large computational procedure. However, once it has been calculated for a specific system, matrices Q and R are stored and used for any acquired projection on that system. Implementation of the QR decomposition in order to take more advantage of the sparsity of the system matrix is discussed.

Determination of the Interaction Position of Gamma Photons in Monolithic Scintillators Using Neural Network Fitting

In Positron Emission Tomography (PET) detectors based on monolithic scintillators, the photon interaction position needs to be estimated from the light distribution (LD) on the photodetector pixels. Due to the finite size of the scintillator volume, the symmetry of the LD is truncated everywhere except for the crystal center. This effect produces a poor estimation of the interaction positions towards the edges, an especially critical situation when linear algorithms, such as Center of Gravity (CoG), are used. When all the crystal faces are painted black, except the one in contact with the photodetector, the LD can be assumed to behave as the inverse square law, providing a simple theoretical model. Using this LD model, the interaction coordinates can be determined by means of fitting each event to a theoretical distribution. In that sense, the use of neural networks (NNs) has been shown to be an effective alternative to more traditional fitting techniques as nonlinear least squares (LS). The multilayer perceptron is one type of NN which can model non-linear functions well and can be trained to accurately generalize when presented with new data. In this work we have shown the capability of NNs to approximate the LD and provide the interaction coordinates of γ-photons with two different photodetector setups. One experimental setup was based on analog Silicon Photomultipliers (SiPMs) and a charge division diode network, whereas the second setup was based on digital SiPMs (dSiPMs). In both experiments NNs minimized border effects. Average spatial resolutions of 1.9 ±0.2 mm and 1.7 ±0.2 mm for the entire crystal surface were obtained for the analog and dSiPMs approaches, respectively.

Analysis of the Statistical Moments of the Scintillation Light Distribution With dSiPMs

In γ-ray detectors, monolithic scintillation crystals offer the possibility of preserving the scintillation light distribution especially when painted black. The statistical moments of this distribution provide accurate information on the three photon impact coordinates, including their depth of interaction (DOI). Digital SiPMs (dSiPMs) return digital information based on pixels about the collected light distribution, since the signal is a digital sum of the trigger bins. In this work we present, for the first time, an accurate analysis of the statistical moments of the light distribution using monolithic painted black crystals and state-of-the-art dSiPMs. Two 32.6 ×32.6 mm2 monolithic LYSO crystals covering the entire photodetectors area have been used in coincidence with 10 mm in thickness. The photosensor tiles were kept at a stable temperature of T = 20 °C. Energy resolution of about 18% was reached in relation to the zeroth moment. The first moment, related to the impact position, determined a spatial resolution of about 3 mm near the crystal center, but quadratically degrading towards the crystal borders. The DOI resolution, measured by means of the second moment, was found to be nearing 4 mm in the crystal center region. The third order moment, the so-called skewness, is related to the degree of truncation and once calibrated minimizes the compression effects. A corrected spatial resolution of about 3 mm was then measured for the entire crystal surface. DOI resolution improved at the crystals center, reaching 3.5 mm, but a degradation towards the borders remained due to truncation of the scintillation light distribution.