Search PubMed⌕ Search

Biomedical subjects

G T Gullberg

Publications and source records attributed to G T Gullberg.

At least 37 records · Page 2Linked to original sources
PubMed1998

Projection space image reconstruction using strip functions to calculate pixels more "natural" for modeling the geometric response of the SPECT collimator.

The spatially varying geometric response of the collimator-detector system in single photon emission computed tomography (SPECT) causes loss in resolution, shape distortions, reconstructed density nonuniformity, and quantitative inaccuracies. A projection space image reconstruction algorithm is used to correct these reconstruction artifacts. The projectors F use strip functions to calculate pixels more "natural" for modeling the two-dimensional (2-D) geometric response of the SPECT collimator transaxially to the axis of rotation. These projectors are defined by summing the intersection of an array of multiple strips rotated at equal angles to approximate the ideal system geometric response of the collimator. Two projection models were evaluated for modeling the system geometric response function. For one projector each strip is of equal weight, for the other projector a Gaussian weighting is used. Parallel beam and fan beam projections of a physical three-dimensional (3-D) Hoffman brain phantom and a Jaszczak cold rod phantom were used to evaluate the geometric response correction. Reconstructions were obtained by using the singular value decomposition (SVD) method and the iterative conjugate gradient algorithm to solve for q in the imaging equation FGq = p, where p is the projection measurement. The projector F included the new models for the geometric response, whereas, the backprojector G did not always model the geometric response in order to increase the computational speed. The final reconstruction was obtained by sampling the backprojection Gq at a discrete array of points. Reconstructions produced by the two proposed projectors showed improved resolution when compared against a unit-strip "natural" pixel model, the conventional image pixelized model with ray tracing to calculate the geometric response, and the filtered backprojection algorithm. When the reconstruction is displayed on fine grid points, the continuity and resolution of the image is preserved without the ring artifacts seen in the unit-strip "natural" pixel model. With present computing power, the geometric response correction using the proposed projection space reconstruction approach is not yet feasible for routine clinical use.

Algorithms↗
PubMed1998

Iterative reconstruction of fluorine-18 SPECT using geometric point response correction.

UNLABELLED: This article demonstrates resolution recovery in 18F SPECT image reconstruction by using an iterative algorithm that corrects for the system geometric response. METHODS: Patient and phantom studies were performed using a Picker PRISM 3000 three-detector SPECT system (Picker International, Inc., Cleveland, OH) to image 18F with 511 keV collimators. A measured point response function of the imaging system was used in an iterative reconstruction algorithm in which the projector and backprojector modeled the system point response function by using an efficient layer-by-layer blurring technique. The blurring function was a five-element kernel in the shape of a cross. The iterative reconstruction algorithm was an ordered-subset maximum-likelihood expectation maximization algorithm. RESULTS: The iterative reconstruction algorithm with geometric response correction showed an improvement in resolution over the filtered backprojection reconstruction and the iterative reconstruction without correction. CONCLUSION: The proposed iterative reconstruction algorithm with geometric response correction is efficient and effective with significant resolution recovery.

Aged↗
PubMed1998

Effect of truncated projections on defect detection in attenuation-compensated fanbeam cardiac SPECT.

UNLABELLED: For some camera systems used in cardiac SPECT, the limited field of view of fanbeam-collimated detectors produces truncation in the projection data. This truncation may generate artifacts and distortions in the transmission CT images and in the attenuation-corrected myocardial SPECT images (that use the transmission CT images as attenuation maps), thus affecting clinical diagnosis. Concern over this problem stimulated us to evaluate the effect of truncation with human observer performance studies. METHODS: A three-dimensional mathematical cardiac-torso phantom that realistically models the attenuation and 201Tl uptake distributions in different organs of the upper torso was used for the investigation. Five degrees of truncation (from 0% to 40% truncation) in the projection data were simulated by the use of five different detector sizes collimated by fanbeam collimators. The attenuation maps were obtained by reconstructing the truncated transmission data generated from these fanbeam geometries. The emission data obtained with the same fanbeam geometries were reconstructed by using the maximum likelihood-expectation maximization algorithm and were corrected for attenuation using the reconstructed attenuation maps. Two observer performance studies were performed. In Study 1, the images were reconstructed without using the body contour as support, whereas in Study 2, the exact body contour was used in reconstructing both the attenuation maps and the attenuation-corrected SPECT images. RESULTS: The results of the receiver operating characteristics analysis indicate that there is very little difference between detection for the various degrees of truncation, and this difference only becomes noticeable for severe truncation of greater than 40%. CONCLUSION: Using the acquisition and processing methods, we found the use of the body contour as support in reconstructing truncated transmission CT and SPECT data reduces the loss of defect detectability in attenuation-compensated myocardial SPECT images due to truncated data.

Algorithms↗
PubMed1997

Reduction of geometric and intensity distortions in echo-planar imaging using a multireference scan.

Echo-planar imaging (EPI) is very sensitive to patient-induced field inhomogeneity caused by susceptibility changes between different anatomical regions. This results in geometric and intensity distortions in the image, especially near tissue/air and tissue/bone interfaces. A new approach is presented to reduce geometric and intensity distortions in EPI. A phase-encoded multireference scan is used to estimate the amplitude and phase errors in the measured signals due to the field inhomogeneity. The EPI data is corrected using both the amplitude and phase of the measured errors. This technique has been evaluated using EPI pulse sequences implemented with conventional gradients and implemented with imaging systems that have special resonating gradients and fast analog to digital converters. The results in both phantom and human studies show that in the absence of object motion the new correction technique can effectively reduce the geometric and intensity distortions.

Cerebellum↗
PubMed1997

An investigation into the effect of input function shape and image acquisition interval on estimates of washin for dynamic cardiac SPECT.

Dynamic cardiac SPECT and PET can be used to measure myocardial perfusion by estimating the kinetic rate constant describing the washing of radioactive-labelled tracers from the blood to the extravascular myocardial tissue. Because of differences in photon statistics and data acquisition techniques, protocols which produce optimal estimates of the washin for dynamic cardiac PET may give suboptimal estimates if applied in dynamic cardiac SPECT. Two important factors in the estimation of washin are the shape of the tracer input function and the image acquisition interval. This study uses computer simulations to investigate the effect of varying the tracer infusion length and image acquisition interval on the bias and variance of estimates of washin obtained with dynamic cardiac SPECT and 99mTc-labelled teboroxime. Bias in parameter estimates can be introduced by aliasing, integration of the time-varying radioactivity by the detector, and detector motion. This bias can be reduced by decreasing the acquisition interval and using a longer-duration input function. However, this results in poor photon statistics, which generate large variance, and can also introduce bias in the estimates of the washin. Our studies indicate that better estimates of the washin are obtained by using an acquisition interval that is of sufficient duration to obtain adequate photon statistics even if this is at the expense of temporal resolution. The increase in bias caused by using a 10 or 20 s acquisition interval instead of a 5 s acquisition interval is minimal when compared with the reduction in variance. Variance in estimates is also reduced by using a sharp input function, resulting in higher peak counts during washin. It is also shown that the variance of estimates of the washin increases generally when faster kinetics are observed. This variance can, however, be reduced by using longer acquisition intervals.

Fourier Analysis↗
PubMed1997

Toward accurate attenuation correction in SPECT without transmission measurements.

The current trend in attenuation correction for single photon emission computed tomography (SPECT) is to measure and reconstruct the attenuation coefficient map using a transmission scan, performed either sequentially or simultaneously with the emission scan. This approach requires dedicated hardware and increases the cost (and in some cases the scanning time) required to produce a clinical SPECT image. Furthermore, if short focal-length fan-beam collimators are used for transmission imaging, the projection data may be truncated, leading to errors in the attenuation coefficient map. Our goal is to obtain information about the attenuation distribution from only the measured emission data by exploiting the fact that only certain attenuation distributions are consistent with a given emission dataset. Ultimately this consistency information will either be used directly to compensate for attenuation or combined with the incomplete information from fan-beam transmission measurements to produce a more accurate attenuation coefficient map. In this manuscript the consistency conditions (which relate the measured SPECT data to the sinogram of the attenuation distribution) are used to find the uniform elliptical attenuation object which is most consistent with the measured emission data. This object is then used for attenuation correction during the reconstruction of the emission data. The method is tested using both simulated and experimentally acquired data from uniformly and nonuniformly attenuating objects. The results show that, for uniform elliptical attenuators, the consistency conditions of the SPECT data can be used to produce an accurate estimate of the attenuation map without performing any transmission measurements. The results also show that, in certain circumstances, the consistency conditions can prove useful for attenuation compensation with nonuniform attenuators.

Algorithms↗
PubMed1997

Implementation of a model-based nonuniform scatter correction scheme for SPECT.

In this paper four scatter-compensation schemes are considered. The four schemes are all based on a previously developed two-dimensional (2-D) scatter model. Reconstruction is achieved using the iterative expectation-maximization maximum-likelihood (EM-ML) algorithm. The schemes consist of: 1) including the model in both the forward and back projector; 2) just including the model in the forward projector; 3) and 4) implementing the model in a subtraction and addition scheme, respectively. Monte Carlo simulated projection data are used to test the accuracy, convergence properties, and noise properties of the four scatter-compensation schemes. Data are simulated for both uniformly and nonuniformly attenuating objects. The results show that all four correction schemes yield images which are similar in terms of accuracy to that obtained from reconstructing scatter-free data. The subtraction scheme is shown to converge faster than the other compensation schemes, both in terms of iterations and actual time required for reconstruction. The scheme in which the model is only used in the forward-projector and the scatter-addition scheme both performs slightly better, in terms of signal-to-noise ratio (SNR), than the subtraction scheme. However, the forward projector scheme requires significantly more CPU time for reconstruction. The correction scheme in which the scatter model was included in both the forward and backprojectors is shown to produce accurate images with SNR's higher than even a perfect scatter rejection scheme. While the scatter correction scheme with the model in both the forward projector and backprojector has superior noise properties to the other algorithms, the results suggest that the faster subtraction/addition schemes will probably prove most useful for routine clinical scatter compensation.

Image Processing, Computer-Assisted↗
PubMed1995

Reduction of phase error ghosting artifacts in thin slice fast spin-echo imaging.

Fast spin-echo (FSE) imaging techniques are very sensitive to the relative phase between the 90 degrees (excitation) RF pulse and the 180 degrees (refocusing) RF pulses. In this paper, it is demonstrated that a phase shift can be created between the excitation and refocusing pulses in such a manner that the received signal is divided into two components of distinctly different phase shifts. The nature of these two components is reviewed. It is demonstrated that ghosting artifacts will occur when images are reconstructed from this received signal. The ghosting is shown to be object dependent. A correction technique is presented which calculates the phase errors among different echoes based on measurements from a single echo train acquired without phase encoding gradients. The results in both phantom and human studies show that this method is capable of reducing the ghosting artifact in thin slice FSE images.

Algorithms↗
PubMed1995

Magnetic resonance imaging verification of a multi-compartment perfusion model for a chromatography gel phantom.

A model of MRI signal intensity was developed for chromatography column phantoms containing Sephadex gels, which were used to simulate tissue perfusion and the exchange of protons between extravascular and intravascular tissue compartments. Computer simulations of two-compartment and three-compartment models were compared in experiments that used two chromatography columns. One column contained a Sephadex G-25 gel that had a smaller pore size and a slower rate of exchange between extrabead and intrabead compartments than did a second column, which contained a G-50 gel with a larger pore size. In both columns the beads were of approximately the same size to simulate the same ratio of intrabead volume to extrabead volume. Slice-selective spin-echo experiments were performed with the columns oriented parallel to each other in a 1.5 T imager and parallel to the magnetic field, with water flowing through each column perpendicular to the slice at mean flow velocities ranging between 0 and 0.66 mm/s. The results of the experiments agreed with computer simulations, which showed that the MRI signal intensity in the perfused columns is a function of the rate of exchange between extrabead and intrabead compartments. The exchange process modifies the transit time of protons passing through an excited region. The simulations also showed that both two-compartment and three-compartment models could be used to fit the experimental data. Experiments with chromatography gel columns are a preliminary investigation into the potential of using MRI to measure the exchange between extravascular and intravascular compartments of endogenous protons during blood perfusion of biological tissue.

Body Fluid Compartments↗
PubMed1995

A transmission-map-based scatter correction technique for SPECT in inhomogeneous media.

In this paper a method of modeling the distribution of scattered events in emission projection data is developed and applied. This method is based on the use of a transmission map to define the inhomogeneous scattering object. The key point is the use of the set of line integrals calculated as part of the attenuation correction technique, as the basis of a model of the distribution of scattered events. The probability of a photon being scattered through a given angle and being detected in the emission energy window is approximated using a Gaussian function. The parameters of this Gaussian are determined using Monte Carlo generated parallel-beam scatter line spread functions from a nonuniformly attenuating phantom. The model is incorporated into a two-dimensional projector-backprojector and used with the Expectation-Maximization-Maximum-Likelihood algorithm for the reconstruction of fan-beam phantom data. The correction is shown to perform well for a phantom that varies slowly in the axial direction. For the more clinically realistic situation of a torso phantom, the method produces improvements in terms of blood pool to myocardium contrast, but does not restore the contrast to the level exhibited in a reconstruction from "scatter free" data.

Humans↗
PubMed1995

An investigation of the effect of finite system resolution and photon noise on the bias and precision of dynamic cardiac SPECT parameters.

The exchange of a radionuclide (and similarly of oxygen) between blood and the myocardium is a dynamic process. Dynamic SPECT offers the possibility of directly quantifying the kinetic parameters which describe this process. In this paper we investigate and quantify the effect of typical SPECT system resolution and photon counting statistics on the bias and precision of dynamic cardiac SPECT parameters. System resolution and photon noise are only two of the image degrading processes which occur in SPECT. Therefore, the bias and precision quoted should be viewed as a lower limit on those which could be expected with an experimental study. Dynamic SPECT projection data are simulated using a realistic human torso phantom. Data are simulated assuming both perfect system resolution and a system resolution typical of a clinical SPECT system. A triple detector SPECT system which acquires a full set of projection data in 10 s using continuous detector motion is modelled. Kinetic parameters are estimated using a number of myocardial regions of interest. The results show that the rate constant characterizing the washing of activity into the myocardium is more sensitive to region of interest position than is the washout rate constant. The bias and precision of the dynamic parameters are estimated from multiple realizations of projection data exhibiting various noise levels. The main effect of increased photon noise in the projection data is to decrease the precision of the estimated parameters. However, there is also some evidence that for high noise levels the bias of the parameters may also be affected.

Biophysical Phenomena↗
PubMed1994

Application of convergent-beam collimation and simultaneous transmission emission tomography to cardiac single-photon emission computed tomography.

Single-photon emission computed tomography (SPECT) is the most commonly performed imaging technique for perfusion studies of the heart and brain. However, these organs are much smaller than the crystal surface of gamma cameras. SPECT sensitivity and resolution can be improved by using fan- and cone-beam collimators to magnify the image of these organs over a larger portion of the crystal face. Special orbits and reconstruction algorithms must be used with convergent-beam acquisitions to prevent image distortion. Differential attenuation of source activity in the chest is one of the most vexing problems in cardiac SPECT, especially with Thallium-201. Multi-headed cameras equipped with convergent-beam collimators allow a transmission image to be obtained at the same time as emission images. Applying a transmission map of the chest attenuation values to the emission images produces a truer picture of source distribution in the heart. This article reviews the technical problems associated with convergent-beam geometry and simultaneous transmission emission tomography SPECT imaging of the heart.

Algorithms↗
PubMed1994

Implementation of Tuy's cone-beam inversion formula.

Tuy's cone-beam inversion formula was modified to develop a cone-beam reconstruction algorithm. The algorithm was implemented for a cone-beam vertex orbit consisting of a circle and two orthogonal lines. This orbit geometry satisfies the cone-beam data sufficiency condition and is easy to implement on commercial single photon emission computed tomography (PECT) systems. The algorithm, which consists of two derivative steps, one rebinning step, and one three-dimensional backprojection step, was verified by computer simulations and by reconstructing physical phantom data collected on a clinical SPECT system. The proposed algorithm gives equivalent results and is as efficient as other analytical cone-beam reconstruction algorithms.

Algorithms↗
PubMed1994

A cone beam SPECT reconstruction algorithm with a displaced center of rotation.

A filtered backprojection (FBP) algorithm is derived based on Feldkamp's FBP algorithm for a cone beam geometry that has a displaced center of rotation. In cone beam single photon emission computed tomography (CB-SPECT) the center of rotation displacement can degrade the reconstructed images. The center of rotation displacement of interest is mechanical shift, which is the displacement of the midplane of the cone beam collimator off the rotation center. Mechanical shift is characterized by two orthogonal components: the shift of the midplane of the cone beam collimator along the direction of the axis of rotation, and the distance between the midline of the cone beam collimator and the axis of rotation. This new algorithm corrects mechanical shift directly by incorporating mechanical shift into the algorithm. This new algorithm is evaluated using both Monte Carlo simulated data and experimentally acquired data. The results demonstrate that this algorithm is able to correct for blurring and the "doughnut" type artifacts caused by system mechanical shift and improve the image resolution.

Algorithms↗
PubMed1994

A simulation of emission and transmission noise propagation in cardiac SPECT imaging with nonuniform attenuation correction.

Transmission computed tomography provides information needed for nonuniform attenuation correction of cardiac single photon emission computed tomography (SPECT). Nonuniform attenuation correction is accomplished using an iterative ML-EM algorithm and a projection-backprojection operation that incorporates attenuation factors measured from the reconstructed transmission map. The precision and accuracy of the attenuation corrected emission reconstruction is a function of emission and transmission statistics. This paper presents an error propagation analysis that uses a mathematical cardiac chest phantom to simulate various combinations of total emission counts C and transmission flux I0 under ideal imaging conditions (without geometric response distortion and without scatter). The spatial average, spatial variance, and accuracy measures for a 4 x 4 pixel region in the heart are tabulated after 30 iterations of the ML-EM algorithm. The confidence intervals for these measures were determined from 1000 realizations of reconstructions from projections randomly generated with the same transmission and emission statistics. It can be shown empirically from the simulation results that the spatial %rms uncertainty for the simulated cardiac region has a simple expression: %rms2 = K1/C+K2/I0(2)+B2 where K1 and K2 are least-square estimates based on the simulation results, and B is the measured spatial %rms uncertainty for the simulation at infinite statistics. For a transmission incident flux of 1500 events per projection bin of 0.712 cm and typical clinical emission events totaling 1 x 10(5), the spatial %rms uncertainty is approximately 14%. At clinical transmission and emission statistics, the statistical noise in the simulated attenuation-corrected reconstructions are dominated by the emission statistics.

Algorithms↗
PubMed1994

Kinetic modeling of teboroxime using dynamic SPECT imaging of a canine model.

UNLABELLED: The tomographic utility of 99mTc-labeled teboroxime has been limited because of its fast washout from the heart, which requires rapid data acquisitions that have not been feasible until the recent development of multidetector SPECT systems. METHODS: Using a three-detector SPECT system to acquire dynamic tomographic data every 10.2 sec, we investigated the potential of modeling the kinetics of teboroxime to develop a sensitive and quantitative measure of cardiac perfusion. Seven studies were performed on four dogs; in three of the studies the LAD artery was occluded. The three-dimensional activity distributions were reconstructed and were corrected for attenuation using a transmission scan. Time-activity curves from the blood and tissue were fit to a two-compartment model with two-way exchange. RESULTS: Performing attenuation correction during the reconstruction process affected the washin parameter k21 significantly (p < 0.0001). The washin parameter k21 also decreased significantly (p < 0.002) when the LAD was occluded. CONCLUSIONS: The results indicate that the washin of teboroxime in myocardial tissue (k21) measured using dynamic SPECT imaging and kinetic modeling is an indicator of myocardial blood flow.

Algorithms↗
PubMed1993

The use of computer-assisted diagnosis in cardiac perfusion nuclear medicine studies: a review (Part 2).

Computer-assisted diagnosis (CAID) is commonly used to evaluate cardiac nuclear medicine studies such as thallium perfusion scans. Part I of this series (Journal of Digital Imaging, Vol 5, No 4, pp 209-222) reviewed the basic theory underlying CAID in nuclear medicine and its use in planar thallium imaging. This review discusses the application of CAID to SPECT perfusion studies.

Artifacts↗