Biomedical subjects
Publications and source records attributed to Donald L Snyder.
Two iterative methods are developed for forming a maximum-likelihood estimate of the attenuation density in a patient or object for transmission tomography when projection data are incomplete. The methods converge monotonically to the same limit points. Results of testing the methods with both simulated and real data are given.
We address the problem of image formation in transmission tomography when metal objects of known composition and shape, but unknown pose, are present in the scan subject. Using an alternating minimization (AM) algorithm, derived from a model in which the detected data are viewed as Poisson-distributed photon counts, we seek to eliminate the streaking artifacts commonly seen in filtered back projection images containing high-contrast objects. We show that this algorithm, which minimizes the I-divergence (or equivalently, maximizes the log-likelihood) between the measured data and model-based estimates of the means of the data, converges much faster when knowledge of the high-density materials (such as brachytherapy applicators or prosthetic implants) is exploited. The algorithm incorporates a steepest descent-based method to find the position and orientation (collectively called the pose) of the known objects. This pose is then used to constrain the image pixels to their known attenuation values, or, for example, to form a mask on the "missing" projection data in the shadow of the objects. Results from two-dimensional simulations are shown in this paper. The extension of the model and methods used to three dimensions is outlined.
The accurate determination of x-ray signal properties is important to several computed tomography (CT) research and development areas, notably for statistical reconstruction algorithms and dose-reduction simulation. The most commonly used model of CT signal formation, assuming monoenergetic x-ray sources with quantum counting detectors obeying simple Poisson statistics, does not reflect the actual physics of CT acquisition. This paper describes a more accurate model, taking into account the energy-integrating detection process, nonuniform flux profiles, and data-conditioning processes. Methods are developed to experimentally measure and theoretically calculate statistical distributions, as well as techniques to analyze CT signal properties. Results indicate the limitations of current models and suggest improvements for the description of CT signal properties.
X-ray computed tomography (CT) images of patients bearing metal intracavitary applicators or other metal foreign objects exhibit severe artifacts including streaks and aliasing. We have systematically evaluated via computer simulations the impact of scattered radiation, the polyenergetic spectrum, and measurement noise on the performance of three reconstruction algorithms: conventional filtered backprojection (FBP), deterministic iterative deblurring, and a new iterative algorithm, alternating minimization (AM), based on a CT detector model that includes noise, scatter, and polyenergetic spectra. Contrary to the dominant view of the literature, FBP streaking artifacts are due mostly to mismatches between FBP's simplified model of CT detector response and the physical process of signal acquisition. Artifacts on AM images are significantly mitigated as this algorithm substantially reduces detector-model mismatches. However, metal artifacts are reduced to acceptable levels only when prior knowledge of the metal object in the patient, including its pose, shape, and attenuation map, are used to constrain AM's iterations. AM image reconstruction, in combination with object-constrained CT to estimate the pose of metal objects in the patient, is a promising approach for effectively mitigating metal artifacts and making quantitative estimation of tissue attenuation coefficients a clinical possibility.