Search

We present the performance of a digital phoswich positron emission tomography (PET) detector, composed by layers of pixilated scintillator arrays, read out by solid state light detectors and an application specific integrated circuit (ASIC). We investigated the use of integrated charge from the scintillation pulses along with time-over-threshold (ToT) to determine the layer of interaction (DOI) in the scintillator. Simulations were performed to assess the effectiveness of the ToT measurements for separating the scintillator events and identifying cross-layer-crystal-scatter (CLCS) events. These simulations indicate that ToT and charge integration from such a detector provide sufficient information to determine the layer of interaction. To demonstrate this in practice, we used a pair of prototype LYSO/BGO detectors. One detector consisted of a 19 × 19 array of 7 mm long LYSO crystals (1.36 mm pitch) coupled to a 16 × 16 array of 8 mm long BGO crystals (1.63 mm pitch). The other detector was similar except the LYSO crystal pitch was 1.63 mm. These detectors were coupled to an 8 × 8 multi-pixel photon counter mounted on a PETsys TOFPET2 ASIC. This high performance ASIC provided digital readout of the integrated charge and ToT from these detectors. We present a method to separate the events from the two scintillator layers using the ToT, and also investigate the performance of this detector. All the crystals within the proposed detector were clearly resolved, and the peak to valley ratio was 11.8 ± 4.0 and 10.1 ± 2.9 for the LYSO and BGO flood images. The measured energy resolution was 9.9% ± 1.3% and 28.5% ± 5.0% respectively for the LYSO and BGO crystals in the phoswich layers. The timing resolution between the LYSO–LYSO, LYSO–BGO and BGO–BGO coincidences was 468 ps, 1.33 ns and 2.14 ns respectively. Results show ToT can be used to identify the crystal layer where events occurred and also identify and reject the majority of CLCS events between layers.

A beta camera has been developed that allows planar imaging of the spatial and temporal distribution of beta particles using a 14 × 14 mm^2 position sensitive avalanche photodiode (PSAPD). This camera system, which we call Betabox, can be directly coupled to microfluidic chips designed for cell incubation or other biological applications. Betabox allows for imaging the cellular uptake of molecular imaging probes labeled with charged particle emitters such as ^(18)F inside these chips. In this work, we investigate the quantitative imaging capabilities of Betabox for ^(18)F beta particles, in terms of background rate, efficiency, spatial resolution, and count rate. Measurements of background and spatial resolution are considered both at room temperature (21 °C ± 1 °C) and at an elevated operating temperature (37 °C ± 1 °C), as is often required for biological assays. The background rate measured with a 4 keV energy cutoff is below 2 cph mm^(−2) at both 21 and 37 °C. The absolute efficiency of Betabox for the detection of ^(18)F positron sources in contact with a PSAPD with the surface passivated from ambient light and damage is 46% ± 1%. The lower detection limit is estimated using the Rose Criterion to be 0.2 cps mm^(−2) for 1 min acquisitions and a 62 × 62 µm^2 pixel size. The upper detection limit is approximately 21 000 cps. The spatial resolution at both 21 and 37 °C ranges from 0.4 mm FWHM at the center of the field of view (FOV), and degrades to 1 mm at a distance of 5 mm away from center yielding a useful FOV of approximately 10 × 10 mm^2. We also investigate the effects on spatial resolution and sensitivity that result from the use of a polymer based microfluidic chip. For these studies we place varying layers of low-density polyethylene (LDPE) between the detector and the source and find that the spatial resolution degrades by ~180 µm for every 100 µm of LDPE film. Sensitivity is reduced by half with the inclusion of ~200 µm of additional LDPE film. Lastly, we demonstrate the practical utilization of Betabox, with an imaging test of its linearity, when coupled to a polydimethylsiloxane microfluidic chip designed for cell based assays.

This paper introduces a mouse atlas registration system (MARS), composed of a stationary top-view x-ray projector and a side-view optical camera, coupled to a mouse atlas registration algorithm. This system uses the x-ray and optical images to guide a fully automatic co-registration of a mouse atlas with each subject, in order to provide anatomical reference for small animal molecular imaging systems such as positron emission tomography (PET). To facilitate the registration, a statistical atlas that accounts for inter-subject anatomical variations was constructed based on 83 organ-labeled mouse micro-computed tomography (CT) images. The statistical shape model and conditional Gaussian model techniques were used to register the atlas with the x-ray image and optical photo. The accuracy of the atlas registration was evaluated by comparing the registered atlas with the organ-labeled micro-CT images of the test subjects. The results showed excellent registration accuracy of the whole-body region, and good accuracy for the brain, liver, heart, lungs and kidneys. In its implementation, the MARS was integrated with a preclinical PET scanner to deliver combined PET/MARS imaging, and to facilitate atlas-assisted analysis of the preclinical PET images.

Dual-energy x-ray computed tomography (DECT) has the capability to decompose attenuation coefficients using two basis functions and has proved its potential in reducing beam-hardening artifacts from reconstructed images. The method typically involves two successive scans with different x-ray tube voltage settings. This work proposes an approach to dual-energy imaging through x-ray beam filtration that requires only one scan and a single tube voltage setting. It has been implemented in a preclinical microCT tomograph with minor modifications. Retrofitting of the microCT scanner involved the addition of an automated filter wheel and modifications to the acquisition and reconstruction software. Results show that beam-hardening artifacts are reduced to noise level. Acquisition of a mu-Compton image is well suited for attenuation-correction of PET images while dynamic energy selection (4D viewing) offers flexibility in image viewing by adjusting contrast and noise levels to suit the task at hand. All dual-energy and single energy reference scans were acquired at the same soft tissue dose level of 50 mGy.

It has been observed that microfluidic chips used for synthesizing (18)F-labeled compounds demonstrate visible light emission without nearby scintillators or fluorescent materials. The origin of the light was investigated and found to be consistent with the emission characteristics from Cerenkov radiation. Since (18)F decays through the emission of high-energy positrons, the energy threshold for beta particles, i.e. electrons or positrons, to generate Cerenkov radiation was calculated for water and polydimethylsiloxane (PDMS), the most commonly used polymer-based material for microfluidic chips. Beta particles emitted from (18)F have a continuous energy spectrum, with a maximum energy that exceeds this energy threshold for both water and PDMS. In addition, the spectral characteristics of the emitted light from (18)F in distilled water were also measured, yielding a broad distribution from 300 nm to 700 nm, with higher intensity at shorter wavelengths. A photograph of the (18)F solution showed a bluish-white light emitted from the solution, further suggesting Cerenkov radiation. In this study, the feasibility of using this Cerenkov light emission as a method for quantitative measurements of the radioactivity within the microfluidic chip in situ was evaluated. A detector previously developed for imaging microfluidic platforms was used. The detector consisted of a charge-coupled device (CCD) optically coupled to a lens. The system spatial resolution, minimum detectable activity and dynamic range were evaluated. In addition, the calibration of a Cerenkov signal versus activity concentration in the microfluidic chip was determined. This novel method of Cerenkov radiation measurements will provide researchers with a simple yet robust quantitative imaging tool for microfluidic applications utilizing beta particles.

Inevitable discrepancies between the mouse tissue optical properties assumed by an experimenter and the actual physiological values may affect the tomographic localization of bioluminescent sources. In a previous work, the simplifying assumption of optically homogeneous tissues led to inaccurate localization of deep sources. Improved results may be obtained if a mouse anatomical map is provided by a high-resolution imaging modality and optical properties are assigned to segmented tissues. In this work, the feasibility of this approach was explored by simulating the effect of different magnitude optical property errors on the image formation process of a combined optical–PET system. Some comparisons were made with corresponding simulations using higher spatial resolution data that are typically attainable by CCD cameras. In addition, simulation results provided insights on some of the experimental conditions that could lead to poor localization of bioluminescent sources. They also provided a rough guide on how accurately tissue optical properties need to be known in order to achieve correct localization of point sources with increasing tissue depth under low background noise conditions.

Many research institutions have a full suite of preclinical tomographic scanners to answer biomedical questions in vivo. Routine multi-modality imaging requires robust registration of images generated by various tomographs. We have implemented a hardware registration method for preclinical imaging that is similar to that used in the combined positron emission tomography (PET)/computed tomography (CT) scanners in the clinic. We designed an imaging chamber which can be rigidly and reproducibly mounted on separate microPET and microCT scanners. We have also designed a three-dimensional grid phantom with 1288 lines that is used to generate the spatial transformation matrix from software registration using a 15-parameter perspective model. The imaging chamber works in combination with the registration phantom synergistically to achieve the image registration goal. We verified that the average registration error between two imaging modalities is 0.335 mm using an in vivo mouse bone scan. This paper also estimates the impact of image misalignment on PET quantitation using attenuation corrections generated from misregistered images. Our technique is expected to produce PET quantitation errors of less than 5%. The methods presented are robust and appropriate for routine use in high throughput animal imaging facilities.

The feasibility and limits in performing tomographic bioluminescence imaging with a combined optical-PET (OPET) system were explored by simulating its image formation process. A micro-MRI based virtual mouse phantom was assigned appropriate tissue optical properties to each of its segmented internal organs at wavelengths spanning the emission spectrum of the firefly luciferase at 37 degrees C. The TOAST finite-element code was employed to simulate the diffuse transport of photons emitted from bioluminescence sources in the mouse. OPET measurements were simulated for single-point, two-point and distributed bioluminescence sources located in different organs such as the liver, the kidneys and the gut. An expectation maximization code was employed to recover the intensity and location of these simulated sources. It was found that spectrally resolved measurements were necessary in order to perform tomographic bioluminescence imaging. The true location of emission sources could be recovered if the mouse background optical properties were known a priori. The assumption of a homogeneous optical property background proved inadequate for describing photon transport in optically heterogeneous tissues and led to inaccurate source localization in the reconstructed images. The simulation results pointed out specific methodological challenges that need to be addressed before a practical implementation of OPET-based bioluminescence tomography is achieved.

MicroPET II is a second-generation animal PET scanner designed for high-resolution imaging of small laboratory rodents. The system consists of 90 scintillation detector modules arranged in three contiguous axial rings with a ring diameter of 16.0 cm and an axial length of 4.9 cm. Each detector module consists of a 14 x 14 array of lutetium oxyorthosilicate (LSO) crystals coupled to a multi-channel photomultiplier tube (MC-PMT) through a coherent optical fibre bundle. Each LSO crystal element measures 0.975 mm x 0.975 mm in cross section by 12.5 mm in length. A barium sulphate reflector material was used between LSO elements leading to a detector pitch of 1.15 mm in both axial and transverse directions. Fused optical fibre bundles were made from 90 microm diameter glass fibres with a numerical aperture of 0.56. Interstitial extramural absorber was added between the fibres to reduce optical cross talk. A charge-division readout circuit was implemented on printed circuit boards to decode the 196 crystals in each array from the outputs of the 64 anode signals of the MC-PMT. Electronics from Concorde Microsystems Inc. (Knoxville, TN) were used for signal amplification, digitization, event qualification, coincidence processing and data capture. Coincidence data were passed to a host PC that recorded events in list mode. Following acquisition, data were sorted into sinograms and reconstructed using Fourier rebinning and filtered hackprojection algorithms. Basic evaluation of the system has been completed. The absolute sensitivity of the microPET II scanner was 2.26% at the centre of the field of view (CFOV) for an energy window of 250-750 keV and a timing window of 10 ns. The intrinsic spatial resolution of the detectors in the system averaged 1.21 mm full width at half maximum (FWHM) when measured with a 22Na point source 0.5 mm in diameter. Reconstructed image resolution ranged from 0.83 mm FWHM at the CFOV to 1.47 mm FWHM in the radial direction, 1.17 mm FWHM in the tangential direction and 1.42 mm FWHM in the axial direction at 1 cm offset from the CFOV. These values represent highly significant improvements over our earlier microPET scanner (approximately fourfold sensitivity increase and 25-35% improvement in linear spatial resolution under equivalent operating conditions) and are expected to be further improved when the system is fully optimized.

We are currently developing a small animal positron emission tomography (PET) scanner with a design goal of 1 microlitre (1 mm3) image resolution. The detectors consist of a 12 x 12 array of 1 x 1 x 10 mm lutetium oxyorthosilicate (LSO) scintillator crystals coupled to a 64-channel photomultiplier tube (PMT) via 5 cm long optical fibre bundles. The optical fibre connection allows a high detector packing fraction despite the dead space surrounding the active region of the PMT. Optical fibre bundles made from different types of glass were tested for light transmission, and also their effects on crystal identification and energy resolution, and compared to direct coupling of the LSO arrays to the PMTs. We also investigated the effects of extramural absorber (EMA) in the fibre bundles. Based on these results, fibre bundles manufactured from F2 glass were selected. We built three pairs of prototype detectors (directly coupled LSO array, fibre bundle without EMA and fibre bundle with EMA) and measured flood histograms, energy resolution, intrinsic spatial resolution and timing resolution. The results demonstrated an intrinsic spatial resolution (FWHM) of 1.12 mm (directly coupled), 1.23 mm (fibre bundle without EMA coupling) and 1.27 mm (fibre bundle with EMA coupling) using an approximately 500 microm diameter Na-22 point source. Using a 330 microm outer diameter steel needle line source filled with F-18, spatial resolution for the detector with the EMA optical fibre bundle improved to 1.05 mm. The respective timing and energy FWHM values were 1.96 ns, 21% (directly coupled), 2.20 ns, 23% (fibre bundle without EMA) and 2.99 ns, 30% (fibre bundle with EMA). The peak-to-valley ratio in the flood histograms was better with EMA (5:1) compared to the optical fibre bundle without EMA (2.5:1), due to the decreased optical cross-talk. In comparison to the detectors used in our current generation microPET scanner, these detectors substantially improve on the spatial resolution, preserve the timing resolution and provide adequate energy resolution for a modern high-resolution animal PET tomograph.

A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 mm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and 18FDG animal studies appear to confirm this observation.

Attenuation correction is one of the important corrections required for quantitative positron emission tomography (PET). This work will compare the quantitative accuracy of attenuation correction using a simple global scale factor with traditional transmission-based methods acquired either with a small animal PET or a small animal x-ray computed tomography (CT) scanner. Two phantoms (one mouse-sized and one rat-sized) and two animal subjects (one mouse and one rat) were scanned in CTI Concorde Microsystem's MicroPET(R) Focus(TM) for emission and transmission data and in ImTek's MicroCAT(TM) II for transmission data. PET emission image values were calibrated against a scintillation well counter. Results indicate that the scale factor method of attenuation correction places the average measured activity concentration about the expected value, without correcting for the cupping artefact from attenuation. Noise analysis in the phantom studies with the PET-based method shows that noise in the transmission data increases the noise in the corrected emission data. The CT-based method was accurate and delivered low-noise images suitable for both PET data correction and PET tracer localization.

Monte Carlo simulation is an essential tool in emission tomography that can assist in the design of new medical imaging devices, the optimization of acquisition protocols and the development or assessment of image reconstruction algorithms and correction techniques. GATE, the Geant4 Application for Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE allows the description of time-dependent phenomena such as source or detector movement, and source decay kinetics. This feature makes it possible to simulate time curves under realistic acquisition conditions and to test dynamic reconstruction algorithms. This paper gives a detailed description of the design and development of GATE by the OpenGATE collaboration, whose continuing objective is to improve, document and validate GATE by simulating commercially available imaging systems for PET and SPECT. Large effort is also invested in the ability and the flexibility to model novel detection systems or systems still under design. A public release of GATE licensed under the GNU Lesser General Public License can be downloaded at http:/www-lphe.epfl.ch/GATE/. Two benchmarks developed for PET and SPECT to test the installation of GATE and to serve as a tutorial for the users are presented. Extensive validation of the GATE simulation platform has been started, comparing simulations and measurements on commercially available acquisition systems. References to those results are listed. The future prospects towards the gridification of GATE and its extension to other domains such as dosimetry are also discussed. ispartof: Physics in medicine and biology vol:49 issue:19 pages:4543-4561 ispartof: location:England status: published

We describe a method for normalization in 3D PET for use with maximum a posteriori (MAP) or other iterative model-based image reconstruction methods. This approach is an extension of previous factored normalization methods in which we include separate factors for detector sensitivity, geometric response, block effects and deadtime. Since our MAP reconstruction approach already models some of the geometric factors in the forward projection, the normalization factors must be modified to account only for effects not already included in the model. We describe a maximum likelihood approach to joint estimation of the count-rate independent normalization factors, which we apply to data from a uniform cylindrical source. We then compute block-wise and block-profile deadtime correction factors using singles and coincidence data, respectively, from a multiframe cylindrical source. We have applied this method for reconstruction of data from the Concorde microPET P4 scanner. Quantitative evaluation of this method using well-counter measurements of activity in a multicompartment phantom compares favourably with normalization based directly on cylindrical source measurements.

The microPET Primate 4-ring system (P4) is an animal PET tomograph with a 7.8 cm axial extent, a 19 cm diameter transaxial field of view (FOV) and a 22 cm animal port. The system is composed of 168 detector modules, each with an 8 x 8 array of 2.2 x 2.2 x 10 mm3 lutetium oxyorthosilicate crystals, arranged as 32 crystal rings 26 cm in diameter. The detector crystals are coupled to a Hamamatsu R5900-C8 PS-PMT via a 10 cm long optical fibre bundle. The detectors have a timing resolution of 3.2 ns, an average energy resolution of 26%, and an average intrinsic spatial resolution of 1.75 mm. The system operates in 3D mode without inter-plane septa, acquiring data in list mode. The reconstructed image spatial resolution ranges from 1.8 mm at the centre to 3 mm at 4 cm radial offset. The tomograph has a peak system sensitivity of 2.25% at the centre of the FOV with a 250-750 keV energy window. The noise equivalent count rate peaks at 100-290 kcps for representative object sizes. Images from two phantoms and three different types of laboratory animal demonstrate the advantage of the P4 system over the original prototype microPET. including its threefold improvement in sensitivity and a large axial FOV sufficient to image an entire mouse in a single bed position.