Search

Purpose: In recent years, the authors and others have been exploring the use of penalized-likelihood sinogram-domain smoothing and restoration approaches for emission and transmission tomography. The motivation for this strategy was initially pragmatic: to provide a more computationally feasible alternative to fully iterative penalized-likelihood image reconstruction involving expensive backprojections and reprojections, while still obtaining some of the benefits of the statistical modeling employed in penalized-likelihood approaches. In this work, the authors seek to compare the two approaches in greater detail., Methods: The sinogram-domain strategy entails estimating the "ideal" line integrals needed for reconstruction of an activity or attenuation distribution from the set of noisy, potentially degraded tomographic measurements by maximizing a penalized-likelihood objective function. The objective function models the data statistics as well as any degradation that can be represented in the sinogram domain. The estimated line integrals can then be input to analytic reconstruction algorithms such as filtered backprojection (FBP). The authors compare this to fully iterative approaches maximizing similar objective functions., Results: The authors present mathematical analyses based on so-called equivalent optimization problems that establish that the approaches can be made precisely equivalent under certain restrictive conditions. More significantly, by use of resolution-variance tradeoff studies, the authors show that they can yield very similar performance under more relaxed, realistic conditions., Conclusions: The sinogram- and image-domain approaches are equivalent under certain restrictive conditions and can perform very similarly under more relaxed conditions. The match is particularly good for fully sampled, high-resolution CT geometries. One limitation of the sinogram-domain approach relative to the image-domain approach is the difficulty of imposing additional constraints, such as image non-negativity.

Most X-ray tubes comprise a rotating anode that is bombarded with electrons to produce X-rays. A substantial amount of heat is generated, and to increase the area of the anode exposed to the electrons, without increasing the apparent size of the focal spot, the focal track of the anode is generally beveled with a very shallow angle (typically 5 degrees-7 degrees in a computed tomography (CT) tube). Due to the line focus principle, this allows a fairly large area of the focal track to be exposed to electrons while retaining a fairly small effective projected focal spot. One side effect of anode angulation is that the focal spot appears different from different positions in the detector array; the effective focal spot size at a constant distance from the tube will be larger for a peripheral detector channel than for a central one. These differences in the effective size of the focal spot across the field-of-view lead to worse resolution in the periphery than in the center of reconstructed images. In this work we describe a method for achieving more uniform resolution in fanbeam CT images by correcting for these focal spot angulation effects. We do so by modeling the effects as a series of local blurrings in the space of transmitted CT intensities and determining the effective coefficients of the corresponding discrete convolutions. The effect of these blurrings can then be compensated for in the sinogram domain through the use of a penalized-likelihood sinogram restoration model we have recently developed.

In this paper, we derive a monotonic penalized-likelihood algorithm for image reconstruction in X-ray fluorescence computed tomography (XFCT) when the attenuation maps at the energies of the fluorescence X-rays are unknown. In XFCT, a sample is irradiated with pencil beams of monochromatic synchrotron radiation that stimulate the emission of fluorescence X-rays from atoms of elements whose K- or L-edges lie below the energy of the stimulating beam. Scanning and rotating the object through the beam allows for acquisition of a tomographic dataset that can be used to reconstruct images of the distribution of the elements in question. XFCT is a stimulated emission tomography modality, and it is thus necessary to correct for attenuation of the incident and fluorescence photons. The attenuation map is, however, generally known only at the stimulating beam energy and not at the energies of the various fluorescence X-rays of interest. We have developed a penalized-likelihood image reconstruction strategy for this problem. The approach alternates between updating the distribution of a given element and updating the attenuation map for that element's fluorescence X-rays. The approach is guaranteed to increase the penalized likelihood at each iteration. Because the joint objective function is not necessarily concave, the approach may drive the solution to a local maximum. To encourage the algorithm to seek out a reasonable local maximum, we include in the objective function a prior that encourages a relationship, based on physical considerations, between the fluorescence attenuation map and the distribution of the element being reconstructed.

We present a statistically principled sinogram smoothing approach for X-ray computed tomography (CT) with the intent of reducing noise-induced streak artifacts. These artifacts arise in CT when some subset of the transmission measurements capture relatively few photons because of high attenuation along the measurement lines. Attempts to reduce these artifacts have focused on the use of adaptive filters that strive to tailor the degree of smoothing to the local noise levels in the measurements. While these approaches involve loose consideration of the measurement statistics to determine smoothing levels, they do not explicitly model the statistical distributions of the measurement data. In this paper, we present an explicitly statistical approach to sinogram smoothing in the presence of photon-starved measurements. It is an extension of a nonparametric sinogram smoothing approach using penalized Poisson-likelihood functions that we have previously developed for emission tomography. Because the approach explicitly models the data statistics, it is naturally adaptive--it will smooth more variable measurements more heavily than it does less variable measurements. We find that it significantly reduces streak artifacts and noise levels without comprising image resolution.

X-ray fluorescence computed tomography (XFCT) is an emerging imaging modality that allows for the reconstruction of the distribution of nonradioactive elements within a sample from measurements of fluorescence x-rays produced by irradiation of the sample with monochromatic synchrotron radiation. XFCT is not a transmission tomography modality, but rather a stimulated emission tomography modality and thus correction for attenuation of the incident and fluorescence photons is essential if accurate images are to be obtained. In this work, we develop and characterize an approximate analytic approach to image reconstruction with attenuation correction in XFCT that is applicable when the incident beam attenuation is uniform and when a factor involving fluorescence attenuation and solid angle effects satisfies a certain approximation. When these conditions hold, we demonstrate that the XFCT imaging equation reduces to the exponential Radon transform, which can be readily inverted. The necessary approximation worsens as the total fluorescence attenuation in the sample grows, but the approach is found to be relatively robust as the approximation breaks down. In a long-axis, small solid angle geometry the proposed approach performs comparably to a previously proposed, more computationally expensive approximate method across a range of attenuation levels. In a short-axis, large solid angle geometry, the proposed approach is found to outperform this previous method.

In this paper, the sampling and aliasing consequences of employing a quarter-detector-offset (QDO) in helical computed tomography (CT) are analyzed. QDO is often used in conventional CT to reduce in-plane aliasing by eliminating data redundancies to improve radial sampling. In helical CT, these same redundancies are exploited to improve longitudinal sampling and so it might seem ill-advised to employ QDO. The relative merit of the two geometries for helical CT is studied by conducting a multidimensional sampling analysis of projection-space sampling as well as a Fourier crosstalk analysis of crosstalk among the object's Fourier basis components. Both a standard fanbeam helical CT geometry and a hypothetical parallel-beam CT geometry, which helps illuminate the more complicated fanbeam results, are analyzed. Using the sampling analysis, it was found that the use of QDO leads to very different spectral tiling than arise when not using QDO. However, due to the shape of the essential support of the projection data spectra that arises in practice, both configurations lead to very similar or identical amounts of spectral overlap. This perspective also predicts the spatially variant longitudinal aliasing that has been observed in helical CT. The crosstalk results were consistent with those of the multidimensional sampling analysis. Thus, from the standpoint of aliasing and crosstalk, no compelling difference is found between the two geometries.

We have investigated the effect of computed tomography (CT) image reconstruction algorithm on the performance of our automated lung nodule detection method. Commercial CT scanners offer a choice of several algorithms for the reconstruction of projection data into transaxial images. Different algorithms produce images with substantially different properties that are apparent not only quantitatively, but also through visual assessment. During some clinical thoracic CT examinations, patient scans are reconstructed with multiple reconstruction algorithms. Thirty-eight such cases were collected to form two databases: one with patient projection data reconstructed with the "standard" reconstruction algorithm and the other with the same patient projection data reconstructed with the "lung" reconstruction algorithm. The automated nodule detection method was applied to both databases. This method is based on gray-level-thresholding techniques to segment the lung regions from each CT section to create a segmented lung volume. Further gray-level-thresholding techniques are applied within the segmented lung volume to identify a set of lung nodule candidates. Rule-based and linear discriminant classifiers are used to differentiate between lung nodule candidates that correspond to actual nodules and those that correspond to non-nodules. The automated method that was applied to both databases was exactly the same, except that the classifiers were calibrated separately for each database. For comparison, the classifier then was trained on one database and tested independently on the other database. When applied to the databases in this manner, the automated method demonstrated overall a similar level of performance, indicating an encouraging degree of robustness.