Bryan W. Reutter

Neural correlates of woman face processing by 2-month-old infants.

The age of 2 months marks a turn in the development of face processing in humans with the emergence of recognition based on internal feature configuration. We studied the neural bases of this early cognitive expertise, critical for adaptive behavior in the social world, by mapping with positron emission tomography the brain activity of 2-month-old alert infants while looking at unknown woman faces. We observed the activation of a distributed network of cortical areas that largely overlapped the adult face-processing network, including the so-called fusiform face area. We also evidenced the activation of left superior temporal and inferior frontal gyri, regions associated, in adults, with language processing. These findings demonstrates that cognitive development proceeds early in functionally active interconnected cortical areas despite the fact they have not all yet reached full metabolic maturation.

List-mode maximum-likelihood reconstruction applied to positron emission mammography (PEM) with irregular sampling

Presents a preliminary study of list-mode likelihood reconstruction of images for a rectangular positron emission tomograph (PET) specifically designed to image the human breast. The prospective device consists of small arrays of scintillation crystals for which depth of interaction is estimated. Except in very rare instances, the number of annihilation events detected is expected to be far less than the number of distinguishable events. If one were to histogram the acquired data, most histogram bins would remain vacant. Therefore, it seems natural to investigate the efficacy of processing events one at a time rather than processing the data in histogram format. From a reconstruction perspective, the new tomograph presents a challenge in that the rectangular geometry leads to irregular radial and angular sampling, and the field of view extends completely to the detector faces. Simulations are presented that indicate that the proposed tomograph can detect 8-mm-diameter spherical tumors with a tumor-to-background tracer density ratio of 3:1 using realistic image acquisition parameters. Spherical tumors of 4-mm diameter are near the limit of detectability with the image acquisition parameters used. Expressions are presented to estimate the loss of image contrast due to Compton scattering.

Real-time system for respiratory-cardiac gating in positron tomography

A Macintosh-based signal processing system has been developed to support simultaneous respiratory and cardiac gating on the ECAT EXACT HR PET scanner. Using the Lab-View real-time software environment, the system reads analog inputs from a pneumatic respiratory bellows and an ECG monitor to compute an appropriate histogram memory location for the PET data. Respiratory state is determined by the bellows signal amplitude; cardiac state is based on the time since the last R-wave. These two states are used in a 2D lookup table to determine a combined respiratory-cardiac state. A 4-bit address encoding the selected histogram is directed from the system to the ECAT scanner, which dynamically switches the destination of tomograph events as respiratory-cardiac state changes. To test the switching efficiency of the combined Macintosh/ECAT system, a rotating emission phantom was built. Acquisitions with 25 msec states while the phantom was rotating at 240 rpm demonstrate the system could effectively stop motion at this rate, with approximately 5 msec switching time between states.

Non-rigid summing of gated PET via optical flow

A method for summing together datasets from gated cardiac PET acquisitions is described. Optical flow techniques are used to accurately model non-rigid motion present during the cardiac cycle so that a one-to-one mapping is found between each voxel of two gated volumes. Using this mapping, image summing can take place, producing a composite dataset with improved statistics and reduced motion-induced blur. Results using a data from a gated cardiac study on a dog are presented.

Kinetic parameter estimation from SPECT cone-beam projection measurements.

Kinetic parameters are commonly estimated from dynamically acquired nuclear medicine data by first reconstructing a dynamic sequence of images and subsequently fitting the parameters to time-activity curves generated from regions of interest overlaid upon the image sequence. Biased estimates can result from images reconstructed using inconsistent projections of a time-varying distribution of radiopharmaceutical acquired by a rotating SPECT system. If the SPECT data are acquired using cone-beam collimators wherein the gantry rotates so that the focal point of the collimators always remains in a plane, additional biases can arise from images reconstructed using insufficient, as well as truncated, projection samples. To overcome these problems we have investigated the estimation of kinetic parameters directly from SPECT cone-beam projection data by modelling the data acquisition process. To accomplish this it was necessary to parametrize the spatial and temporal distribution of the radiopharmaceutical within the SPECT field of view. In a simulated chest image volume, kinetic parameters were estimated for simple one-compartment models for four myocardial regions of interest. Myocardial uptake and washout parameters estimated by conventional analysis of noiseless simulated cone-beam data had biases ranging between 3-26% and 0-28%, respectively. Parameters estimated directly from the noiseless projection data were unbiased as expected, since the model used for fitting was faithful to the simulation. Statistical uncertainties of parameter estimates for 10,000,000 events ranged between 0.2-9% for the uptake parameters and between 0.3-6% for the washout parameters.

Dynamic single photon emission computed tomography—basic principles and cardiac applications

The very nature of nuclear medicine, the visual representation of injected radiopharmaceuticals, implies imaging of dynamic processes such as the uptake and wash-out of radiotracers from body organs. For years, nuclear medicine has been touted as the modality of choice for evaluating function in health and disease. This evaluation is greatly enhanced using single photon emission computed tomography (SPECT), which permits three-dimensional (3D) visualization of tracer distributions in the body. However, to fully realize the potential of the technique requires the imaging of in vivo dynamic processes of flow and metabolism. Tissue motion and deformation must also be addressed. Absolute quantification of these dynamic processes in the body has the potential to improve diagnosis. This paper presents a review of advancements toward the realization of the potential of dynamic SPECT imaging and a brief history of the development of the instrumentation. A major portion of the paper is devoted to the review of special data processing methods that have been developed for extracting kinetics from dynamic cardiac SPECT data acquired using rotating detector heads that move as radiopharmaceuticals exchange between biological compartments. Recent developments in multi-resolution spatiotemporal methods enable one to estimate kinetic parameters of compartment models of dynamic processes using data acquired from a single camera head with slow gantry rotation. The estimation of kinetic parameters directly from projection measurements improves bias and variance over the conventional method of first reconstructing 3D dynamic images, generating time-activity curves from selected regions of interest and then estimating the kinetic parameters from the generated time-activity curves. Although the potential applications of SPECT for imaging dynamic processes have not been fully realized in the clinic, it is hoped that this review illuminates the potential of SPECT for dynamic imaging, especially in light of new developments that enable measurement of dynamic processes directly from projection measurements.

Scintillator Non-Proportionality: Present Understanding and Future Challenges

Scintillator non-proportionality (the fact that the conversion factor between the energy deposited in a scintillator and the number of visible photons produced is not constant) has been studied both experimentally and theoretically for ~50 years. Early research centered on the dependence of the conversion factor on the species of the ionizing radiation (gamma, alpha, beta, proton, etc.), and researchers during the 1960s discovered a strong correlation between the scintillation efficiency and the ionization density. In more recent years, non-proportionality has been proposed as the reason why the energy resolution of most scintillators is worse than that predicted by counting statistics. While much progress has been made, there are still major gaps in our understanding of both the fundamental causes of non-proportionality and their quantitative link to scintillator energy resolution. This paper summarizes the present state of knowledge on the nature of the light-yield non-proportionality and its effect on energy resolution.

Automated 3-D segmentation of respiratory-gated PET transmission images

As a preliminary step toward performing respiration compensated attenuation correction of respiratory-gated cardiac PET data, we acquired and automatically segmented respiratory-gated transmission data for a dog breathing on its own under gas anesthesia. Transmission data were acquired for 20 min on a CTI/Siemens ECAT EXACT HR (47-slice) scanner. Two respiratory gates were obtained using data from a pneumatic bellows placed around the dogs chest. For each respiratory gate, torso and lung surfaces were segmented automatically using a differential 3-D image edge detection algorithm. Three-dimensional visualizations showed that during inspiration the heart translated about 4 mm transversely and the diaphragm translated about 9 mm inferiorly. The observed respiratory motion of the canine heart and diaphragm suggests that respiration compensated attenuation correction may be necessary for accurate quantitation of high-resolution respiratory-gated human cardiac PET data. Our automated image segmentation results suggest that respiration compensated segmented attenuation correction may be possible using respiratory-gated transmission data obtained with as little as 3 min of acquisition time per gate.

Effects of temporal modelling on the statistical uncertainty of spatiotemporal distributions estimated directly from dynamic SPECT projections

Artefacts can result when reconstructing a dynamic image sequence from inconsistent single photon emission computed tomography (SPECT) projection data acquired by a slowly rotating gantry. The artefacts can lead to biases in kinetic parameters estimated from time-activity curves generated by overlaying volumes of interest on the images. Insufficient sampling and truncation of projections by cone-beam collimators can cause additional artefacts. To overcome these sources of bias in conventional image based dynamic data analysis, we have been investigating the estimation of time-activity curves and kinetic model parameters directly from dynamic SPECT projection data by modelling the spatial and temporal distribution of the radiopharmaceutical throughout the projected field of view. In the present work, we perform Monte Carlo simulations to study the effects of the temporal modelling on the statistical variability of the reconstructed spatiotemporal distributions. The simulations utilize fast methods for fully four-dimensional (4D) direct estimation of spatiotemporal distributions and their statistical uncertainties, using a spatial segmentation and temporal B-splines. The simulation results suggest that there is benefit in modelling higher orders of temporal spline continuity. In addition, the accuracy of the time modelling can be increased substantially without unduly increasing the statistical uncertainty, by using relatively fine initial time sampling to capture rapidly changing activity distributions.

A methodology for specifying PET VOIs using multimodality techniques

Volume-of-interest (VOI) extraction for radionuclide and anatomical measurements requires correct identification and delineation of the anatomical feature being studied. The authors have developed a toolset for specifying three dimensional (3-D) VOIs on a multislice positron emission tomography (PET) dataset. The software is particularly suited for specifying cerebral cortex VOIs which represent a particular gyrus or deep brain structure. A registered 3-D magnetic resonance image (MRI) dataset is used to provide high-resolution anatomical information, both as oblique two-dimensional (2-D) sections and as volume renderings of a segmented cortical surface. VOIs are specified indirectly in two dimensions by drawing a stack of 2-D regions on the MRI data. The regions are tiled together to form closed triangular mesh surface models, which are subsequently transformed into the observation space of the PET scanner. Quantification by this method allows calculation of radionuclide activity in the VOIs, as well as their statistical uncertainties and correlations. The methodology for this type of analysis and validation results are presented.