Multi-GPU Acceleration of Branchless Distance Driven Projection and Backprojection for Clinical Helical CT.

Model-based image reconstruction (MBIR) techniques have the potential to generate high quality images from noisy measurements and a small number of projections which can reduce the x-ray dose in patients. These MBIR techniques rely on projection and backprojection to refine an image estimate. One of the widely used projectors for these modern MBIR based technique is called branchless distance driven (DD) projection and backprojection. While this method produces superior quality images, the computational cost of iterative updates keeps it from being ubiquitous in clinical applications. In this paper, we provide several new parallelization ideas for concurrent execution of the DD projectors in multi-GPU systems using CUDA programming tools. We have introduced some novel schemes for dividing the projection data and image voxels over multiple GPUs to avoid runtime overhead and inter-device synchronization issues. We have also reduced the complexity of overlap calculation of the algorithm by eliminating the common projection plane and directly projecting the detector boundaries onto image voxel boundaries. To reduce the time required for calculating the overlap between the detector edges and image voxel boundaries, we have proposed a pre-accumulation technique to accumulate image intensities in perpendicular 2D image slabs (from a 3D image) before projection and after backprojection to ensure our DD kernels run faster in parallel GPU threads. For the implementation of our iterative MBIR technique we use a parallel multi-GPU version of the alternating minimization (AM) algorithm with penalized likelihood update. The time performance using our proposed reconstruction method with Siemens Sensation 16 patient scan data shows an average of 24 times speedup using a single TITAN X GPU and 74 times speedup using 3 TITAN X GPUs in parallel for combined projection and backprojection.

Quantitative Computed Tomography Classification of Lung Nodules: Initial Comparison of 2- and 3-Dimensional Analysis.

OBJECTIVE: The aim of this study was to compare the performance of 2- (2D) and 3-dimensional (3D) quantitative computed tomography (CT) methods for classifying lung nodules as lung cancer, metastases, or benign. METHODS: Using semiautomated software and computerized analysis, we analyzed more than 50 quantitative CT features of 96 solid nodules in 94 patients, in 2D from a single slice and in 3D from the entire nodule volume. Multivariable logistic regression was used to classify nodule types. Model performance was assessed by the area under the receiver operating characteristic curve (AUC) using leave-one-out cross-validation. RESULTS: The AUC for distinguishing 53 primary lung cancers from 18 benign nodules and 25 metastases ranged from 0.79 to 0.83 and was not significantly different for 2D and 3D analyses (P = 0.29-0.78). Models distinguishing metastases from benign nodules were statistically significant only by 3D analysis (AUC = 0.84). CONCLUSIONS: Three-dimensional CT methods did not improve discrimination of lung cancer, but may help distinguish benign nodules from metastases.

Image gently campaign back to basics initiative: ten steps to help manage radiation dose in pediatric digital radiography.

OBJECTIVE: The purpose of this review is to summarize 10 steps a practice can take to manage radiation exposure in pediatric digital radiography. CONCLUSION: The Image Gently campaign raises awareness of opportunities for lowering radiation dose while maintaining diagnostic quality of images of children. The newest initiative in the campaign, Back to Basics, addresses methods for standardizing the approach to pediatric digital radiography, highlighting challenges related to the technology in imaging of patients of widely varying body sizes.

Motion-compensated scheme for sequential scanned statistical iterative dual-energy CT reconstruction.

Objective.Dual-energy computed tomography (DECT) has been widely used to reconstruct numerous types of images due its ability to better discriminate tissue properties. Sequential scanning is a popular dual-energy data acquisition method as it requires no specialized hardware. However, patient motion between two sequential scans may lead to severe motion artifacts in DECT statistical iterative reconstructions (SIR) images. The objective is to reduce the motion artifacts in such reconstructions.Approach.We propose a motion-compensation scheme that incorporates a deformation vector field into any DECT SIR. The deformation vector field is estimated via the multi-modality symmetric deformable registration method. The precalculated registration mapping and its inverse or adjoint are then embedded into each iteration of the iterative DECT algorithm.Main results.Results from a simulated and clinical case show that the proposed framework is capable of reducing motion artifacts in DECT SIRs. Percentage mean square errors in regions of interest in the simulated and clinical cases were reduced from 4.6% to 0.5% and 6.8% to 0.8%, respectively. A perturbation analysis was then performed to determine errors in approximating the continuous deformation by using the deformation field and interpolation. Our findings show that errors in our method are mostly propagated through the target image and amplified by the inverse matrix of the combination of the Fisher information and Hessian of the penalty term.Significance.We have proposed a novel motion-compensation scheme to incorporate a 3D registration method into the joint statistical iterative DECT algorithm in order to reduce motion artifacts caused by inter-scan motion, and successfully demonstrate that interscan motion corrections can be integrated into the DECT SIR process, enabling accurate imaging of radiological quantities on conventional SECT scanners, without significant loss of either computational efficiency or accuracy.

New exposure indicators for digital radiography simplified for radiologists and technologists.

OBJECTIVE: The purpose of this article is to educate radiologists and technologists about the clinically relevant portion of the new digital radiography standards. CONCLUSION: Both the International Electrotechnical Commission (IEC standard 62494-1) and the American Association of Physicists in Medicine (AAPM Task Group 116) have developed similar standards for monitoring exposure in digital radiography to eliminate proprietary and confusing terminology. Radiologists and technologists will need to learn three new terms--exposure index, target exposure index, and deviation index--to understand the new standards.

Towards subpercentage uncertainty proton stopping-power mapping via dual-energy CT: Direct experimental validation and uncertainty analysis of a statistical iterative image reconstruction method.

PURPOSE: To assess the potential of a joint dual-energy computerized tomography (CT) reconstruction process (statistical image reconstruction method built on a basis vector model (JSIR-BVM)) implemented on a 16-slice commercial CT scanner to measure high spatial resolution stopping-power ratio (SPR) maps with uncertainties of less than 1%. METHODS: JSIR-BVM was used to reconstruct images of effective electron density and mean excitation energy from dual-energy CT (DECT) sinograms for 10 high-purity samples of known density and atomic composition inserted into head and body phantoms. The measured DECT data consisted of 90 and 140 kVp axial sinograms serially acquired on a Philips Brilliance Big Bore CT scanner without beam-hardening corrections. The corresponding SPRs were subsequently measured directly via ion chamber measurements on a MEVION S250 superconducting synchrocyclotron and evaluated theoretically from the known sample compositions and densities. Deviations of JSIR-BVM SPR values from their theoretically calculated and directly measured ground-truth values were evaluated for our JSIR-BVM method and our implementation of the Hünemohr-Saito (H-S) DECT image-domain decomposition technique for SPR imaging. A thorough uncertainty analysis was then performed for five different scenarios (comparison of JSIR-BVM stopping-power ratio/stopping power (SPR/SP) to International Commission on Radiation Measurements and Units benchmarks; comparison of JSIR-BVM SPR to measured benchmarks; and uncertainties in JSIR-BVM SPR/SP maps for patients of unknown composition) per the Joint Committee for Guides in Metrology and the Guide to Expression of Uncertainty in Measurement, including the impact of uncertainties in measured photon spectra, sample composition and density, photon cross section and I-value models, and random measurement uncertainty. Estimated SPR uncertainty for three main tissue groups in patients of unknown composition and the weighted proportion of each tissue type for three proton treatment sites were then used to derive a composite range uncertainty for our method. RESULTS: Mean JSIR-BVM SPR estimates deviated by less than 1% from their theoretical and directly measured ground-truth values for most inserts and phantom geometries except for high-density Delrin and Teflon samples with SPR error relative to proton measurements of 1.1% and -1.0% (head phantom) and 1.1% and -1.1% (body phantom). The overall root-mean-square (RMS) deviations over all samples were 0.39% and 0.52% (head phantom) and 0.43% and 0.57% (body phantom) relative to theoretical and directly measured ground-truth SPRs, respectively. The corresponding RMS (maximum) errors for the image-domain decomposition method were 2.68% and 2.73% (4.68% and 4.99%) for the head phantom and 0.71% and 0.87% (1.37% and 1.66%) for the body phantom. Compared to H-S SPR maps, JSIR-BVM yielded 30% sharper and twofold sharper images for soft tissues and bone-like surrogates, respectively, while reducing noise by factors of 6 and 3, respectively. The uncertainty (coverage factor k = 1) of the DECT-to-benchmark values comparison ranged from 0.5% to 1.5% and is dominated by scanning-beam photon-spectra uncertainties. An analysis of the SPR uncertainty for patients of unknown composition showed a JSIR-BVM uncertainty of 0.65%, 1.21%, and 0.77% for soft-, lung-, and bony-tissue groups which led to a composite range uncertainty of 0.6-0.9%. CONCLUSIONS: Observed JSIR-BVM SPR estimation errors were all less than 50% of the estimated k = 1 total uncertainty of our benchmarking experiment, demonstrating that JSIR-BVM high spatial resolution, low-noise SPR mapping is feasible and is robust to variations in the geometry of the scanned object. In contrast, the much larger H-S SPR estimation errors are dominated by imaging noise and residual beam-hardening artifacts. While the uncertainties characteristic of our current JSIR-BVM implementation can be as large as 1.5%, achieving < 1% total uncertainty is feasible by improving the accuracy of scanner-specific scatter-profile and photon-spectrum estimates. With its robustness to beam-hardening artifact, image noise, and variations in phantom size and geometry, JSIR-BVM has the potential to achieve high spatial-resolution SPR mapping with subpercentage accuracy and estimated uncertainty in the clinical setting.

Noise-resolution tradeoffs in x-ray CT imaging: a comparison of penalized alternating minimization and filtered backprojection algorithms.

PURPOSE: In comparison with conventional filtered backprojection (FBP) algorithms for x-ray computed tomography (CT) image reconstruction, statistical algorithms directly incorporate the random nature of the data and do not assume CT data are linear, noiseless functions of the attenuation line integral. Thus, it has been hypothesized that statistical image reconstruction may support a more favorable tradeoff than FBP between image noise and spatial resolution in dose-limited applications. The purpose of this study is to evaluate the noise-resolution tradeoff for the alternating minimization (AM) algorithm regularized using a nonquadratic penalty function. METHODS: Idealized monoenergetic CT projection data with Poisson noise were simulated for two phantoms with inserts of varying contrast (7%-238%) and distance from the field-of-view (FOV) center (2-6.5 cm). Images were reconstructed for the simulated projection data by the FBP algorithm and two penalty function parameter values of the penalized AM algorithm. Each algorithm was run with a range of smoothing strengths to allow quantification of the noise-resolution tradeoff curve. Image noise is quantified as the standard deviation in the water background around each contrast insert. Modulation transfer functions (MTFs) were calculated from six-parameter model fits to oversampled edge-spread functions defined by the circular contrast-insert edges as a metric of local resolution. The integral of the MTF up to 0.5 1p/mm was adopted as a single-parameter measure of local spatial resolution. RESULTS: The penalized AM algorithm noise-resolution tradeoff curve was always more favorable than that of the FBP algorithm. While resolution and noise are found to vary as a function of distance from the FOV center differently for the two algorithms, the ratio of noises when matching the resolution metric is relatively uniform over the image. The ratio of AM-to-FBP image variances, a predictor of dose-reduction potential, was strongly dependent on the shape of the AM's nonquadratic penalty function and was also strongly influenced by the contrast of the insert for which resolution is quantified. Dose-reduction potential, reported here as the fraction (%) of FBP dose necessary for AM to reconstruct an image with comparable noise and resolution, for one penalty parameter value of the AM algorithm was found to vary from 70% to 50% for low-contrast and high-contrast structures, respectively, and from 70% to 10% for the second AM penalty parameter value. However, the second penalty, AM-700, was found to suffer from poor low-contrast resolution when matching the high-contrast resolution metric with FBP. CONCLUSIONS: The results of this simulation study imply that penalized AM has the potential to reconstruct images with similar noise and resolution using a fraction (10%-70%) of the FBP dose. However, this dose-reduction potential depends strongly on the AM penalty parameter and the contrast magnitude of the structures of interest. In addition, the authors' results imply that the advantage of AM can be maximized by optimizing the nonquadratic penalty function to the specific imaging task of interest. Future work will extend the methods used here to quantify noise and resolution in images reconstructed from real CT data.

Validation of CT dose-reduction simulation.

The objective of this research was to develop and validate a custom computed tomography dose-reduction simulation technique for producing images that have an appearance consistent with the same scan performed at a lower mAs (with fixed kVp, rotation time, and collimation). Synthetic noise is added to projection (sinogram) data, incorporating a stochastic noise model that includes energy-integrating detectors, tube-current modulation, bowtie beam filtering, and electronic system noise. Experimental methods were developed to determine the parameters required for each component of the noise model. As a validation, the outputs of the simulations were compared to measurements with cadavers in the image domain and with phantoms in both the sinogram and image domain, using an unbiased root-mean-square relative error metric to quantify agreement in noise processes. Four-alternative forced-choice (4AFC) observer studies were conducted to confirm the realistic appearance of simulated noise, and the effects of various system model components on visual noise were studied. The "just noticeable difference (JND)" in noise levels was analyzed to determine the sensitivity of observers to changes in noise level. Individual detector measurements were shown to be normally distributed (p > 0.54), justifying the use of a Gaussian random noise generator for simulations. Phantom tests showed the ability to match original and simulated noise variance in the sinogram domain to within 5.6% +/- 1.6% (standard deviation), which was then propagated into the image domain with errors less than 4.1% +/- 1.6%. Cadaver measurements indicated that image noise was matched to within 2.6% +/- 2.0%. More importantly, the 4AFC observer studies indicated that the simulated images were realistic, i.e., no detectable difference between simulated and original images (p = 0.86) was observed. JND studies indicated that observers' sensitivity to change in noise levels corresponded to a 25% difference in dose, which is far larger than the noise accuracy achieved by simulation. In summary, the dose-reduction simulation tool demonstrated excellent accuracy in providing realistic images. The methodology promises to be a useful tool for researchers and radiologists to explore dose reduction protocols in an effort to produce diagnostic images with radiation dose "as low as reasonably achievable".

Experimental implementation of a joint statistical image reconstruction method for proton stopping power mapping from dual-energy CT data.

PURPOSE: To experimentally commission a dual-energy CT (DECT) joint statistical image reconstruction (JSIR) method, which is built on a linear basis vector model (BVM) of material characterization, for proton stopping power ratio (SPR) estimation. METHODS: The JSIR-BVM method builds on the relationship between the energy-dependent photon attenuation coefficients and the proton stopping power via a pair of BVM component weights. The two BVM component images are simultaneously reconstructed from the acquired DECT sinograms and then used to predict the electron density and mean excitation energy (I-value), which are required by the Bethe equation for SPR computation. A post-reconstruction image-based DECT method, which utilizes the two separate CT images reconstructed via the scanner's software, was implemented for comparison. The DECT measurement data were acquired on a Philips Brilliance scanner at 90 and 140 kVp for two phantoms of different sizes. Each phantom contains 12 different soft and bony tissue surrogates with known compositions. The SPR estimation results were compared to the reference values computed from the known compositions. The difference of the computed water equivalent path lengths (WEPL) across the phantoms between the two methods was also compared. RESULTS: The overall root-mean-square (RMS) of SPR estimation error of the JSIR-BVM method are 0.33% and 0.37% for the head- and body-sized phantoms, respectively, and all SPR estimates of the test samples are within 0.7% of the reference ground truth. The image-based method achieves overall RMS errors of 2.35% and 2.50% for the head- and body-sized phantoms, respectively. The JSIR-BVM method also reduces the pixel-wise random variation by 4-fold to 6-fold within homogeneous regions compared to the image-based method. The average differences between the JSIR-BVM method and the image-based method are 0.54% and 1.02% for the head- and body-sized phantoms, respectively. CONCLUSION: By taking advantage of an accurate polychromatic CT data model and a model-based DECT statistical reconstruction algorithm, the JSIR-BVM method accounts for both systematic bias and random noise in the acquired DECT measurement data. Therefore, the JSIR-BVM method achieves good accuracy and precision on proton SPR estimation for various tissue surrogates and object sizes. In contrast, the experimentally achievable accuracy of the image-based method may be limited by the uncertainties in the image formation process. The result suggests that the JSIR-BVM method has the potential for more accurate SPR prediction compared to post-reconstruction image-based methods in clinical settings.

Use of computed tomography scans for cochlear implants.

While 3-dimensional (3D) imaging by computed tomography has long been desirable for research and treatment of cochlear-implant patients, technical challenges have limited its wide application. Recent developments in scanner hardware and image processing techniques now allow image quality improvements that make clinical applications feasible. Validation experiments were performed to characterize a new methodology and its imaging performance.

Neonatal chest computed radiography: image processing and optimal image display.

OBJECTIVE: The purpose of this study was to determine soft-copy image display preferences of brightness, latitude, and detail contrast for neonatal chest computed radiography to establish a baseline for future work on low-dose imaging. CONCLUSION: Observers preferred brighter images with higher detail contrast and narrow to middle latitude for soft-copy display compared with the typical screen-film hard-copy appearance. Future research on low-dose neonatal chest imaging will be facilitated by an understanding of optimal soft-copy image display.

Comparison of standard- and low-radiation-dose CT for quantification of emphysema.

OBJECTIVE: This study was performed to compare standard- and low-radiation-dose techniques in the CT quantification of emphysema. MATERIALS AND METHODS: The study population consisted of 36 men and 20 women who were current or former heavy smokers and underwent standard-dose (effective tube current, 100-250 mAs) chest CT at our institution within 6 months of having undergone low-dose (effective tube current, 30-60 mAs) chest CT. All CT scans were reconstructed at 5-mm slice thickness with a smooth filter. CT-measured lung volume, mean and median lung attenuation, and percentage of lung volume with attenuation lower than multiple thresholds (emphysema index values) were compared by Pearson correlation, two-tailed and paired Student's t tests, and regression analysis. RESULTS: There were no significant differences in mean attenuation (-848 vs -846 H, p > 0.35) for the low dose and the standard dose or in median lung attenuation (-879 vs -878 H, p > 0.66). Low- and standard-dose emphysema indexes were correlated at all attenuation thresholds (r = 0.86-0.97). Mean emphysema indexes were higher on the low-dose scans, but the mean difference at all thresholds was less than 3%. The differences were significant (p < 0.05) only at the lower index thresholds, correlated with differences in lung volume (r < or = 0.86), and increased with greater differences in dose. CONCLUSION: Low-dose technique has minimal effect on CT quantification of emphysema.

Statistical reconstruction for x-ray computed tomography using energy-integrating detectors.

Statistical image reconstruction (SR) algorithms have the potential to significantly reduce x-ray CT image artefacts because they use a more accurate model than conventional filtered backprojection and can incorporate effects such as noise, incomplete data and nonlinear detector response. Most SR algorithms assume that the CT detectors are photon-counting devices and generate Poisson-distributed signals. However, actual CT detectors integrate energy from the x-ray beam and exhibit compound Poisson-distributed signal statistics. This study presents the first assessment of the impact on image quality of the resultant mismatch between the detector and signal statistics models assumed by the sinogram data model and the reconstruction algorithm. A 2D CT projection simulator was created to generate synthetic polyenergetic transmission data assuming (i) photon-counting with simple Poisson-distributed signals and (ii) energy-weighted detection with compound Poisson-distributed signals. An alternating minimization (AM) algorithm was used to reconstruct images from the data models (i) and (ii) for a typical abdominal scan protocol with incident particle fluence levels ranging from 10(5) to 1.6 x 10(6) photons/detector. The images reconstructed from data models (i) and (ii) were compared by visual inspection and image-quality figures of merit. The reconstructed image quality degraded significantly when the means were mismatched from the assumed model. However, if the signal means are appropriately modified, images from data models (i) and (ii) did not differ significantly even when SNR is very low. While data-mean mismatches characteristic of the difference between particle-fluence and energy-fluence transmission can cause significant streaking and cupping artefacts, the mismatch between the actual and assumed CT detector signal statistics did not significantly degrade image quality once systematic data means mismatches were corrected.

In vivo estimates of the position of advanced bionics electrode arrays in the human cochlea.

OBJECTIVES: A new technique for determining the position of each electrode in the cochlea is described and applied to spiral computed tomography data from 15 patients implanted with Advanced Bionics HiFocus I, Ij, or Helix arrays. METHODS: ANALYZE imaging software was used to register 3-dimensional image volumes from patients' preoperative and postoperative scans and from a single body donor whose unimplanted ears were scanned clinically, with micro computed tomography and with orthogonal-plane fluorescence optical sectioning (OPFOS) microscopy. By use of this registration, we compared the atlas of OPFOS images of soft tissue within the body donor's cochlea with the bone and fluid/ tissue boundary available in patient scan data to choose the midmodiolar axis position and judge the electrode position in the scala tympani or scala vestibuli, including the distance to the medial and lateral scalar walls. The angular rotation 0 degrees start point is a line joining the midmodiolar axis and the middle of the cochlear canal entry from the vestibule. RESULTS: The group mean array insertion depth was 477 degrees (range, 286 degrees to 655 degrees). The word scores were negatively correlated (r = -0.59; p = .028) with the number of electrodes in the scala vestibuli. CONCLUSIONS: Although the individual variability in all measures was large, repeated patterns of suboptimal electrode placement were observed across subjects, underscoring the applicability of this technique.

On two-parameter models of photon cross sections: application to dual-energy CT imaging.

The goal of this study is to evaluate the theoretically achievable accuracy in estimating photon cross sections at low energies (20-1000 keV) from idealized dual-energy x-ray computed tomography (CT) images. Cross-section estimation from dual-energy measurements requires a model that can accurately represent photon cross sections of any biological material as a function of energy by specifying only two characteristic parameters of the underlying material, e.g., effective atomic number and density. This paper evaluates the accuracy of two commonly used two-parameter cross-section models for postprocessing idealized measurements derived from dual-energy CT images. The parametric fit model (PFM) accounts for electron-binding effects and photoelectric absorption by power functions in atomic number and energy and scattering by the Klein-Nishina cross section. The basis-vector model (BVM) assumes that attenuation coefficients of any biological substance can be approximated by a linear combination of mass attenuation coefficients of two dissimilar basis substances. Both PFM and BVM were fit to a modern cross-section library for a range of elements and mixtures representative of naturally occurring biological materials (Z = 2-20). The PFM model, in conjunction with the effective atomic number approximation, yields estimated the total linear cross-section estimates with mean absolute and maximum error ranges of 0.6%-2.2% and 1%-6%, respectively. The corresponding error ranges for BVM estimates were 0.02%-0.15% and 0.1%-0.5%. However, for photoelectric absorption frequency, the PFM absolute mean and maximum errors were 10.8%-22.4% and 29%-50%, compared with corresponding BVM errors of 0.4%-11.3% and 0.5%-17.0%, respectively. Both models were found to exhibit similar sensitivities to image-intensity measurement uncertainties. Of the two models, BVM is the most promising approach for realizing dual-energy CT cross-section measurement.

Properties of preprocessed sinogram data in x-ray computed tomography.

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.

Pose estimation of known objects during transmission tomographic image reconstruction.

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.

Image reconstruction for transmission tomography when projection data are incomplete.

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.

Line Integral Alternating Minimization Algorithm for Dual-Energy X-Ray CT Image Reconstruction.

We propose a new algorithm, called line integral alternating minimization (LIAM), for dual-energy X-ray CT image reconstruction. Instead of obtaining component images by minimizing the discrepancy between the data and the mean estimates, LIAM allows for a tunable discrepancy between the basis material projections and the basis sinograms. A parameter is introduced that controls the size of this discrepancy, and with this parameter the new algorithm can continuously go from a two-step approach to the joint estimation approach. LIAM alternates between iteratively updating the line integrals of the component images and reconstruction of the component images using an image iterative deblurring algorithm. An edge-preserving penalty function can be incorporated in the iterative deblurring step to decrease the roughness in component images. Images from both simulated and experimentally acquired sinograms from a clinical scanner were reconstructed by LIAM while varying the regularization parameters to identify good choices. The results from the dual-energy alternating minimization algorithm applied to the same data were used for comparison. Using a small fraction of the computation time of dual-energy alternating minimization, LIAM achieves better accuracy of the component images in the presence of Poisson noise for simulated data reconstruction and achieves the same level of accuracy for real data reconstruction.

Reduction of motion blurring artifacts using respiratory gated CT in sinogram space: a quantitative evaluation.

Techniques have been developed for reducing motion blurring artifacts by using respiratory gated computed tomography (CT) in sinogram space and quantitatively evaluating the artifact reduction. A synthetic sinogram was built from multiple scans intercepting a respiratory gating window. A gated CT image was then reconstructed using the filtered back-projection algorithm. Wedge phantoms, developed for quantifying the motion artifact reduction, were scanned while being moved using a computer-controlled linear stage. The resulting artifacts appeared between the high and low density regions as an apparent feature with a Hounsfield value that was the average of the two regions. A CT profile through these regions was fit using two error functions, each modeling the partial-volume averaging characteristics for the unmoving phantom. The motion artifact was quantified by determining the apparent distance between the two functions. The blurring artifact had a linear relationship with both the speed and the tangent of the wedge angles. When gating was employed, the blurring artifact was reduced systematically at the air-phantom interface. The gated image of phantoms moving at 20 mm/s showed similar blurring artifacts as the nongated image of phantoms moving at 10 mm/s. Nine patients were also scanned using the synchronized respiratory motion technique. Image artifacts were evaluated in the diaphragm, where high contrast interfaces intercepted the imaging plane. For patients, this respiratory gating technique reduced the blurring artifacts by 9%-41% at the lung-diaphragm interface.