Analytical study of the performance of a multilayer positron computed tomography scanner

Efficient 3-D TOF PET reconstruction using view-grouped histo-images: DIRECT-direct image reconstruction for TOF

For modern time-of-flight (TOF) positron emission tomography (PET) systems, in which the number of possible lines of response and TOF bins is much larger than the number of acquired events, the most appropriate reconstruction approaches are considered to be list-mode methods. However, their shortcomings are relatively high computational costs for reconstruction and for sensitivity matrix calculation. Efficient treatment of TOF data within the proposed DIRECT approach is obtained by 1) angular (azimuthal and co-polar) grouping of TOF events to a set of views as given by the angular sampling requirements for the TOF resolution, and 2) deposition (weighted-histogramming) of these grouped events, and correction data, into a set of "histo-images," one histo-image per view. The histo-images have the same geometry (voxel grid, size and orientation) as the reconstructed image. The concept is similar to the approach involving binning of the TOF data into angularly subsampled histo-projections-projections expanded in the TOF directions. However, unlike binning into histo-projections, the deposition of TOF events directly into the image voxels eliminates the need for tracing and/or interpolation operations during the reconstruction. Together with the performance of reconstruction operations directly in image space, this leads to a very efficient implementation of TOF reconstruction algorithms. Furthermore, the resolution properties are not compromised either, since events are placed into the image elements of the desired size from the beginning. Concepts and efficiency of the proposed data partitioning scheme are demonstrated in this work by using the DIRECT approach in conjunction with the row-action maximum-likelihood (RAMLA) algorithm.

A 2D/3D hybrid PET scanner with rotating partial slice-septa and its quantitative procedures

This paper presents a PET scanner capable of acquiring projection data in three-dimensional (3D) and two-dimensional (2D) modes simultaneously. The scanner has rotating partial slice-septa, and coincidence events are stored as 2D data or as 3D data depending on whether the lines of response are collimated by the septa or not. 68Ge/Ga rod sources can be set on the rotating septa, and a transmission scan for attenuation correction is performed in the 2D mode. The scanner allows simultaneous 3D-emission/2D-transmission scanning or post-injection transmission scanning with little cross-talk. A blank scan for detector normalization is also performed with the rotating rod sources in the 2D mode, from which we can derive the normalizing factors in both modes. The 3D/2D difference method is available for scatter correction, even in a dynamic study where the source distribution is changing. A 'summation method' is proposed as a new image reconstruction algorithm, in which the high- and low-frequency components of images are reconstructed from the 3D and 2D data respectively. In this method, most of the scatter contribution in the 3D data is removed by high-pass filtering, not by subtracting estimated scatter distribution, and hence the method is expected to be robust for scatter from outside the axial field of view. Computer simulations revealed that the rotating partial septa offer a single-scatter to true ratio similar to that of the conventional full septa if the depth of the partial septa is properly lengthened, with a small increase in multiple scattering.

Assessment of response function in two PET scanners with and without interplane septa

The authors have assessed the response function both experimentally and theoretically for two commercial tomographs: CTI 931/08-12 and CTI 953B with and without interplane septa. Monte Carlo simulations were undertaken using the GEANT package from CERN. Spatial resolution (tomographic and axial) was calculated for line sources at various positions in the field of view. Sensitivity and scatter fraction (SF) were calculated for various source geometries as a function of energy discrimination. A very realistic response function in positron emission tomography (PET) is obtained by Monte Carlo methods, using global parameters to account for unsimulated phenomena such as scintillation light transport inside a detector block and its sharing among the various phototubes. Minor discrepancies remain for sensitivity and SF at high energy thresholds and may probably be explained by introducing the observed dispersion in the energy response for the various crystals within a detector block.

Scatter fraction: measurement and correction

The concept of scatter in Positron Emission Tomography is reviewed regarding origin and influence on data. Different ways to measure and correct for scatter are discussed.

Sensitivity and cross calibration

In Positron Emission Tomography it is possible to obtain quantitative information on the distribution of positron emitting radiopharmaceuticals. The different components influencing the sensitivity and the cross calibration of PET systems are discussed especially with respect to the quantitative nature of the PET technique.

On estimating the loss of quantification in PET due to finite detector resolution

In positron emission computed tomography, images are reconstructed from projection data of the slice. Such projection data are inherently subject to bias owing to the detection of coincident annihilation photons with detectors of finite resolution. The influence of these unavoidable effects on reconstructed data is reviewed and a straight forward method is outlined for calculating the loss of quantification of small objects owing to finite detector resolution. A sample calculation using idealized point spread functions is presented.

A fast reconstruction algorthm for stationary positron emission tomography based on a modified em algorithm

An efficient iterative reconstruction method for positron emission tomography (PET) is presented. The algorithm is basically an enhanced EM (expectation maximization) algorithm with improved frequency response. High-frequency components of the ratio of measured to calculated projections are extracted and are taken into account for the iterative correction of image density in such a way that the correction is performed with a uniform efficiency over the image plane and with a flat frequency response. As a result, the convergence speed is not so sensitive to the image pattern or matrix size as the standard EM algorithm, and nonuniformity of the spatial resolution is significantly improved. Nonnegativity of the reconstructed image is preserved. Simulation studies have been made assuming two PET systems: a scanning PET with ideal sampling and a stationary PET with sparse sampling. In the latter, a "bank array" of detectors is employed to improve the sampling in the object plane. The new algorithm provides satisfactory images by two or three iterations starting from a flat image in either case. The behavior of convergence is monitored by evaluating the root mean square of C(b)-1 where C(b) is the correction factor for pixel b in the EM algorithm. The value decreases rapidly and monotonically with iteration number. Although the theory is not accurate enough to assure the stability of convergence, the algorithm is promising to achieve significant saving in computation compared to the standard EM algorithm.

Imaging Performance of a Dynamic Positron Emission Tomograph: Positome IIIp

This paper documents modifications to an older PET system to improve its dead time, scatter fraction, and spatial resolution in high count rate, short duration studies. A new dual-tapered collimator reduces scatter by 33 percent while providing excellent resolution uniformity in all slices. A data encoding scheme produces uniformly sampled parallel projections from the coincidence data in real time while the detector array executes an orbital motion. The image uniformity, scatter compensation and high count-rate performance have been validated up to 40 kBq/cc in a 20 cm flood source. The errors in image quantification due to counting statistics, live time, and random counts are estimated from repeated measurements on a contrast phantom at high count-rates. The effects of two methods of scatter compensation on image contrast are shown in contrast phantoms and a typical glucose utilization study. Blood flow measurements using 0-15 labeled water bolus method, made under different physiological conditions, reflected the changes expected. The true count efficiency of 75 kcps/(uCi/cc) permits these studies to be done with only 500-900 MBq (13-25 mCi) injected activity.

Positron emission tomography imaging--technical considerations

Positron imaging instrumentation has improved rapidly in the last few years. Scanners currently under development are beginning to approach fundamental limits set by positron range and noncolinearity effects. This report reviews the latest developments in positron emission tomography (PET) instrumentation, emphasizing the development of coding schemes that reduce the complexity and cost of high-resolution scanners. The relative benefits of using time-of-flight (TOF) information is discussed as well.

Axial resolution and the value of interpolating scan in multislice positron computed tomography

We have calculated the aperture function of a positron computed tomograph (PCT) with computer simulation, and evaluated the axial resolution of a multislice PCT, Positologica III, both theoretically and experimentally. The axial point spread function (PSF) was approximately a triangle at or near the center of the field, and the sensitivity for the slice decreased significantly as the source moved away off the image plane. Accordingly, there were low sensitivity areas between an in-plane and the adjacent cross-plane. This invisible region was clinically significant if the object was thin enough in the z-axis. In order to fill up the gaps between adjacent slices, it is valuable to move the patient half the slice interval in the z-axis and perform an " interpolating scan."

Software correction of scatter coincidence in positron CT

This paper describes a software correction of scatter coincidence in positron CT which is based on an estimation of scatter projections from true projections by an integral transform. Kernels for the integral transform are projected distributions of scatter coincidences for a line source at different positions in a water phantom and are calculated by Klein-Nishina's formula. True projections of any composite object can be determined from measured projections by iterative applications of the integral transform. The correction method was tested in computer simulations and phantom experiments with Positologica. The results showed that effects of scatter coincidence are not negligible in the quantitation of images, but the correction reduces them significantly.

Performance of positron imaging systems as a function of energy threshold and shielding depth

The effect of energy discrimination and shielding on positron imaging data quality was investigated using a detector pair to simulate a ring positron emission tomograph. Formulas are presented relating the sensitivity, random fraction, and scatter fraction for a detector pair to the same parameters for a ring system. Data were fitted to detector pair expressions for the variation of the above parameters with shielding depth in order to obtain information on the effect of energy threshold level. These fitted curves were used to determine the sensitivity, random fraction, and scatter fraction, as well as an overall data quality factor as a function of energy threshold level and shielding depth. Data were obtained for both NaI(T1) and BGO detector types using activity levels in the range of 1.5 muCi/cm3. Results show that for NaI(T1) detectors, the lowest possible energy threshold level is optimal, with the corresponding optimal shielding depth determined by the level of activity to be imaged. For BGO detectors, a tradeoff exists between energy thresholds of 100-400 keV and shielding depths of 15-30 cm with smaller shielding depths requiring higher energy thresholds.

Performance study of single-slice positron emission tomography scanners by monte carlo techniques

A Monte Carlo simulation of the gamma ray transport within a single-slice positron emission tomograph has been generated to study the effects of system parameters on performance. Included in the simulation are the radioactive source distribution, collimators, and detectors with intercrystal septa. Data are first presented to show the coincidence and singles sensitivities as a function of ring radius. Then, for a fixed radius of 26 cm, the variation of sensitivities are shown as a function of the following variables: slice thickness, patient port size, intercrystal septum dimensions, lower energy discriminator level, and coincidence fan angle. Simulation-generated sensitivity data are compared with experimental values for several tomographs andgood agreement is obtained. Discrepancies between two definitions used in experimentally determining scatter fractions are discussed. The Monte Carlo simulation shows that small radii rings have an effective count rate (quality factor) that is more than 90 percent of that for larger rings at low and moderate activity levels (=/< 0.25 muCi cm(-3)), contrary to what is predicted from analytical calculations. It is concluded that small radius rings are better suited for low dose-rate static studies, while larger radius rings are preferred for high dose-rate dynamic studies.