Experimental evaluation of the resolution improvement provided by a silicon PET probe.

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 × 16 arrays of 1.4 mm × 1.4 mm pixels and the other in 40 × 26 arrays of 1.0 mm × 1.0 mm 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 × 6 arrays of 6 mm × 12 mm × 30 mm 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.5 mm 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.8 mm and 4.0 mm phantom features respectively, were observed, while previously unresolvable 3.2 mm 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.

Silicon detectors for combined MR-PET and MR-SPECT imaging.

Silicon based devices can extend PET-MR and SPECT-MR imaging to applications, where their advantages in performance outweigh benefits of high statistical counts.Silicon is in many ways an excellent detector material with numerous advantages, among others: excellent energy and spatial resolution, mature processing technology, large signal to noise ratio, relatively low price, availability, versatility and malleability. The signal in silicon is also immune to effects of magnetic field at the level normally used in MR devices. Tests in fields up to 7 T were performed in a study to determine effects of magnetic field on positron range in a silicon PET device. The curvature of positron tracks in direction perpendicular to the field's orientation shortens the distance between emission and annihilation point of the positron. The effect can be fully appreciated for a rotation of the sample for a fixed field direction, compressing range in all dimensions. A popular Ga-68 source was used showing a factor of 2 improvement in image noise compared to zero field operation. There was also a little increase in noise as the reconstructed resolution varied between 2.5 and 1.5 mm.A speculative applications can be recognized in both emission modalities, SPECT and PET.Compton camera is a subspecies of SPECT, where a silicon based scatter as a MR compatible part could inserted into the MR bore and the secondary detector could operate in less constrained environment away from the magnet. Introducing a Compton camera also relaxes requirements of the radiotracers used, extending the range of conceivable photon energies beyond 140.5 keV of the Tc-99m.In PET, one could exploit the compressed sub-millimeter range of positrons in the magnetic field. To exploit the advantage, detectors with spatial resolution commensurate to the effect must be used with silicon being an excellent candidate. Measurements performed outside of the MR achieving spatial resolution below 1 mm are reported.

Silicon as an Unconventional Detector in Positron Emission Tomography.

Positron emission tomography (PET) is a widely used technique in medical imaging and in studying small animal models of human disease. In the conventional approach, the 511 keV annihilation photons emitted from a patient or small animal are detected by a ring of scintillators such as LYSO read out by arrays of photodetectors. Although this has been a successful in achieving ~5mm FWHM spatial resolution in human studies and ~1mm resolution in dedicated small animal instruments, there is interest in significantly improving these figures. Silicon, although its stopping power is modest for 511 keV photons, offers a number of potential advantages over more conventional approaches. Foremost is its high spatial resolution in 3D: our past studies show that there is little diffculty in localizing 511 keV photon interactions to ~0.3mm. Since spatial resolution and reconstructed image noise trade off in a highly non-linear manner that depends on the PET instrument response, if high spatial resolution is the goal, silicon may outperform standard PET detectors even though it has lower sensitivity to 511 keV photons. To evaluate silicon in a variety of PET "magnifying glass" configurations, an instrument has been constructed that consists of an outer partial-ring of PET scintillation detectors into which various arrangements of silicon detectors can be inserted to emulate dual-ring or imaging probe geometries. Recent results have demonstrated 0.7 mm FWHM resolution using pad detectors having 16×32 arrays of 1.4mm square pads and setups have shown promising results in both small animal and PET imaging probe configurations. Although many challenges remain, silicon has potential to become the PET detector of choice when spatial resolution is the primary consideration.

Timing performance of the silicon PET insert probe.

Simulation indicates that PET image could be improved by upgrading a conventional ring with a probe placed close to the imaged object. In this paper, timing issues related to a PET probe using high-resistivity silicon as a detector material are addressed. The final probe will consist of several (four to eight) 1-mm thick layers of silicon detectors, segmented into 1 x 1 mm(2) pads, each pad equivalent to an independent p + nn+ diode. A proper matching of events in silicon with events of the external ring can be achieved with a good timing resolution. To estimate the timing performance, measurements were performed on a simplified model probe, consisting of a single 1-mm thick detector with 256 square pads (1.4 mm side), coupled with two VATAGP7s, application-specific integrated circuits. The detector material and electronics are the same that will be used for the final probe. The model was exposed to 511 keV annihilation photons from an (22)Na source, and a scintillator (LYSO)-PMT assembly was used as a timing reference. Results were compared with the simulation, consisting of four parts: (i) GEANT4 implemented realistic tracking of electrons excited by annihilation photon interactions in silicon, (ii) calculation of propagation of secondary ionisation (electron-hole pairs) in the sensor, (iii) estimation of the shape of the current pulse induced on surface electrodes and (iv) simulation of the first electronics stage. A very good agreement between the simulation and the measurements were found. Both indicate reliable performance of the final probe at timing windows down to 20 ns.

Minimising activity and dose with enhanced image quality by radiopharmaceutical administrations.

Owing to the introduction of new diagnostic procedures, such as computed tomography (CT), positron emission tomography (PET) and single photon emission computed tomography (SPECT), the individual dose caused by medical exposures has grown rapidly in the last years. This is especially a subject to radiation protection for nuclear medical diagnosis, since in this case radiopharmaceuticals are administered to the patient, meaning not only a radiation exposure to the diseased tissue but also to the healthy tissues of large parts of the body. 'Minimizing Activity and Dose with Enhanced Image quality by Radiopharmaceutical Administrations' (MADEIRA) is a project cofunded by the European Commission within the Seventh Euratom Framework Programme that aims to improve three-dimensional (3D) nuclear medical imaging technologies significantly. MADEIRA is aiming to improve the efficacy and safety of 3D PET and SPECT functional imaging by optimising the spatial resolution and the signal-to-noise ratio, improving the knowledge of the temporal variation of the radiopharmaceuticals' uptake in and clearance from tumourous and healthy tissues, and evaluation of the corresponding patient dose. Using an optimised imaging procedure that improves the information gained per unit administered dose, MADEIRA aims especially to reduce the dose to healthy tissues of the patient. In this paper, an overall summary of the current achievements will be presented.

Probe-guided localization of cancer deposits using [18F]fluorodeoxyglucose.

In recent years, several probes have been developed to allow for the intraoperative detection of tumour tissue using [18F]fluorodeoxyglucose (FDG). Detector designs include high-energy gamma and beta probes, as well as combination devices with background rejection capabilities. Some laboratory prototypes and commercialized systems have demonstrated reasonable sensitivities for 511 keV photons and /or b particles emitted from 18F for in vivo use. This review focuses on the ability of these devices to detect tumour deposits in the low-contrast environment of the operating room . Important technical and biological factors that influence tumour-to-background contrast are discussed and potential future applications and developments are highlighted. In addition, we evaluate the limited data on absorbed doses resulting from [18F] FDG administration immediately prior to surgery that indicate acceptable levels of radiation exposure to operating room personnel.

Timing in thick silicon detectors - an update.

Thick silicon detectors are becoming widely used with reliable detectors available. They are often considered as a part of coincidence setup, where the timing resolution is of a crucial importance. Since over-biasing of thick detectors is sometimes unpractical, the timing resolution can be compromised in thick detectors. For this article the electric and Ramo fields in a 1.4 by 1.4 mm(2) pad size and 1 mm thick pad detector were calculated. GEANT4 was used to determine the tracks of the interaction electron produced in photon interactions, and the drift of ionized carriers in the detector was simulated. The signals were processed using a virtual preamplifier, a CR-RC shaper with a shaping time of 200 ns and a leading edge discriminator. The distributions of delay of the trigger after the event were compared to the measurements and a good agreement was found, allowing for additional noise in experimental setup. We proposed and evaluated an alternative readout strategy reading signals from 9 nearest pads which greatly reduces the effects of inhomogeneous Ramo field on the timing resolution.

Design and Feasibility Study of a Single Photon Emission Microscope System for Small Animal I-125 Imaging.

This paper presents a design study of a single photon emission microscope (SPEM) system for small animal imaging using I-125 labelled radiotracers. This system is based on the use of a very-high resolution gamma camera coupled to a converging non-multiplexing multiple pinhole collimator. This enables one to "zoom" into a small local region inside the object to extract imaging information with a very high spatial resolution and a reasonable sensitivity for gamma rays emitted from this local region. The SPEM system also includes a pinhole (or multiple pinhole) gamma camera that has a full angular coverage of the entire object. The designed imaging spatial resolution for the SPEM system is between 250 µm to 500 µm FWHM.

A modified uniform Cramer-Rao bound for multiple pinhole aperture design.

This paper presents a modified Uniform Cramer-Rao bound (UCRB) for studying estimator spatial resolution and variance tradeoffs. We proposed to use a resolution constraint that is imposed on mean gradient vectors of achieved estimators and derived the minimum achievable variance for any estimator satisfies this resolution constraint. This approach partially overcomes the limitations of the former UCRB approach based on a bias-gradient norm constraint. We applied this method in a feasibility study of using multiple pinhole apertures for small animal SPECT imaging applications. The SPECT system studied was based on an existing gamma camera. The achievable spatial resolution and variance tradeoffs for systems with different design parameters, such as number of pinholes and pinhole size, were studied.

Feasibility Study of Compton Scattering Enhanced Multiple Pinhole Imager for Nuclear Medicine.

This paper presents a feasibility study of a Compton scattering enhanced (CSE) multiple pinhole imaging system for gamma rays with energy of 140keV or higher. This system consists of a multiple-pinhole collimator, a position sensitive scintillation detector as used in standard Gamma camera, and a Silicon pad detector array, inserted between the collimator and the scintillation detector. The problem of multiplexing, normally associated with multiple pinhole system, is reduced by using the extra information from the detected Compton scattering events. In order to compensate for the sensitivity loss, due to the low probability of detecting Compton scattered events, the proposed detector is designed to collect both Compton scattering and Non-Compton events. It has been shown that with properly selected pinhole spacing, the proposed detector design leads to an improved image quality.

Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography.

We present a dual-energy (DE) transmission computed tomography (CT) reconstruction method. It is statistically motivated and features nonnegativity constraints in the density domain. A penalized weighted least squares (PWLS) objective function has been chosen to handle the non-Poisson noise added by amorphous silicon (aSi:H) detectors. A Gauss-Seidel algorithm has been used to minimize the objective function. The behavior of the method in terms of bias/standard deviation tradeoff has been compared to that of a DE method that is based on filtered back projection (FBP). The advantages of the DE PWLS method are largest for high noise and/or low flux cases. Qualitative results suggest this as well. Also, the reconstructed images of an object with opaque regions are presented. Possible applications of the method are: attenuation correction for positron emission tomography (PET) images, various quantitative computed tomography (QCT) methods such as bone mineral densitometry (BMD), and the removal of metal streak artifacts.

Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction.

This paper presents a new class of algorithms for penalized-likelihood reconstruction of attenuation maps from low-count transmission scans. We derive the algorithms by applying to the transmission log-likelihood a version of the convexity technique developed by De Pierro for emission tomography. The new class includes the single-coordinate ascent (SCA) algorithm and Lange's convex algorithm for transmission tomography as special cases. The new grouped-coordinate ascent (GCA) algorithms in the class overcome several limitations associated with previous algorithms. 1) Fewer exponentiations are required than in the transmission maximum likelihood-expectation maximization (ML-EM) algorithm or in the SCA algorithm. 2) The algorithms intrinsically accommodate nonnegativity constraints, unlike many gradient-based methods. 3) The algorithms are easily parallelizable, unlike the SCA algorithm and perhaps line-search algorithms. We show that the GCA algorithms converge faster than the SCA algorithm, even on conventional workstations. An example from a low-count positron emission tomography (PET) transmission scan illustrates the method.

Model-based estimation for dynamic cardiac studies using ECT.

The authors develop a strategy for joint estimation of physiological parameters and myocardial boundaries using ECT (emission computed tomography). They construct an observation model to relate parameters of interest to the projection data and to account for limited ECT system resolution and measurement noise. The authors then use a maximum likelihood (ML) estimator to jointly estimate all the parameters directly from the projection data without reconstruction of intermediate images. They also simulate myocardial perfusion studies based on a simplified heart model to evaluate the performance of the model-based joint ML estimator and compare this performance to the Cramer-Rao lower bound. Finally, the authors discuss model assumptions and potential uses of the joint estimation strategy.

Model-based estimation with boundary side information or boundary regularization [cardiac emission CT].

The authors have previously developed a model-based strategy for joint estimation of myocardial perfusion and boundaries using ECT (emission computed tomography). They have also reported difficulties with boundary estimation in low contrast and low count rate situations. Here they propose using boundary side information (obtainable from high resolution MRI and CT images) or boundary regularization to improve both perfusion and boundary estimation in these situations. To fuse boundary side information into the emission measurements, the authors formulate a joint log-likelihood function to include auxiliary boundary measurements as well as ECT projection measurements. In addition, they introduce registration parameters to align auxiliary boundary measurements with ECT measurements and jointly estimate these parameters with other parameters of interest from the composite measurements. In simulated PET O-15 water myocardial perfusion studies using a simplified model, the authors show that the joint estimation improves perfusion estimation performance and gives boundary alignment accuracy of <0.5 mm even at 0.2 million counts. They implement boundary regularization through formulating a penalized log-likelihood function. They also demonstrate in simulations that simultaneous regularization of the epicardial boundary and myocardial thickness gives comparable perfusion estimation accuracy with the use of boundary side information.

Thin-film, flat-panel, composite imagers for projection and tomographic imaging.

The recent development of large-area, flat-panel a-Si:H imaging arrays is generally expected to lead to real-time diagnostic and megavoltage X-ray projection imagers with film-cassette-like profiles. While such flat-panel imagers offer numerous advantages over existing fluoroscopic and radiographic imaging devices, the unique properties of the arrays also offer the prospect of detector configurations not previously possible with other real-time technologies. The thin, highly uniform profile of the arrays allows the creation of composite imaging devices in which a flat-panel detector overlies a second imaging detector. A dual-energy (diagnostic and megavoltage) composite imager consisting of a pair of stacked, flat-panel imagers would provide unique information helping to resolve the patient localization and verification problem in megavoltage radiotherapy. In PET or SPECT, attenuation corrections could be obtained by placing a flat-panel array for transmission measurements directly in front of the main emission detector. In this article, the concept of such real-time flat-panel composite imagers is proposed. Specific embodiments of this concept applied toward the resolution of outstanding problems in radiotherapy, PET and SPECT are outlined and calculations and data supporting the feasibility of the concept are presented.

Preconditioning methods for improved convergence rates in iterative reconstructions.

Because of the characteristics of the tomographic inversion problem, iterative reconstruction techniques often suffer from poor convergence rates-especially at high spatial frequencies. By using preconditioning methods, the convergence properties of most iterative methods can be greatly enhanced without changing their ultimate solution. To increase reconstruction speed, spatially invariant preconditioning filters that can be designed using the tomographic system response and implemented using 2-D frequency-domain filtering techniques have been applied. In a sample application, reconstructions from noiseless, simulated projection data, were performed using preconditioned and conventional steepest-descent algorithms. The preconditioned methods demonstrated residuals that were up to a factor of 30 lower than the assisted algorithms at the same iteration. Applications of these methods to regularized reconstructions from projection data containing Poisson noise showed similar, although not as dramatic, behavior.

Acceleration and filtering in the generalized Landweber iteration using a variable shaping matrix.

The generalized Landweber iteration with a variable shaping matrix is used to solve the large linear system of equations arising in the image reconstruction problem of emission tomography. The method is based on the property that once a spatial frequency image component is almost recovered within in in the generalized Landweber iteration, this component will still stay within in during subsequent iterations with a different shaping matrix, as long as this shaping matrix satisfies the convergence criterion for the component. Two different shaping matrices are used: the first recovers low-frequency image components; and the second may be used either to accelerate the reconstruction of high-frequency image components, or to attenuate these components to filter the image. The variable shaping matrix gives results similar to truncated inverse filtering, but requires much less computation and memory, since it does not rely on the singular value decomposition.

Clinical SPRINT imaging. Preliminary results compared to conventional SPECT brain scanning using Tc-99m HMPAO.

We have performed initial clinical studies using the high resolution single photon ring tomograph (SPRINT) and Tc-99m HMPAO. To determine what additional anatomic detail can be depicted using this high resolution, dedicated head, multidetector SPECT device compared to conventional SPECT, six patients with stroke and one normal volunteer were imaged after the injection of 20 mCi Tc-99m HMPAO on a conventional rotating Anger gamma camera (GE-400AC), followed immediately by imaging on SPRINT. Imaging acquisition on the GE-400AC gamma camera was performed using 360 degrees rotation, 64 stops, at 30 sec/stop, yielding an average of 985,714 counts for a 10.0 mm thick slice. GE-400AC images were of good quality, having in-plane full width half maximum (FWHM) resolution of approximately 15 mm. SPRINT acquisition of image data was performed using both the high resolution and high sensitivity apertures, with data collection over 15 or 20 minutes of imaging time accumulating approximately 500,000 counts and 1,000,000 counts, respectively, from patients in a 10.0 mm thick slice, achieving an in-plane FWHM resolution of approximately 8 mm and 10 mm for the two apertures, respectively. Both image resolution and contrast for visualization of gray, white, and cerebral spinal fluid filled brain structures were improved using SPRINT compared with the GE-400AC. We conclude that SPRINT is well suited for brain imaging with Tc-99m HMPAO and is of particular value for applications demanding high resolution.

Testing of local gamma-ray scatter fractions determined by spectral fitting.

The spectral-fitting method of correction for gamma-ray Compton scattering within objects separates the unscattered and scattered components of locally measured energy spectra. Here, we employ a third-order polynomial for the scattering and an approximately constant fitting window. A scatter fraction, defined as total scattered over total unscattered counts within a 20% window, is calculated for each point in our Anger camera images. These scatter fractions are tested against those from Monte-Carlo simulation for 99mTc and against results from semiconductor detector measurements for 131I. A radioactive sphere at several locations within a non-radioactive cylinder and the inverse are imaged for the testing. For one case, reproducibility of the spectral-fitting scatter fraction as a function of the number of unscattered counts within the 20% acceptance window was also determined. With 99mTc, for all cases, the agreement between spectral fitting and the standard estimation method is within 16%. With 131I, for the 'hot' sphere at two locations, the agreement is within 21%. For the 'hot' sphere at the third location (off the cylinder axis towards the camera), the dependence of scatter fraction on transverse distance is good although the absolute values are too large. Scatter fraction reproducibility is within 10% for 1000 or more counts. Therefore, further testing of spectral fitting and initial application to realistic clinical images seem to be in order.