Impact of polychromatic x-ray sources on helical, cone-beam computed tomography and dual-energy methods

Image Domain Multi-Material Decomposition Noise Suppression Through Basis Transformation and Selective Filtering

Spectral CT can provide material characterization ability to offer more precise material information for diagnosis purposes. However, the material decomposition process generally leads to amplification of noise which significantly limits the utility of the material basis images. To mitigate such problem, an image domain noise suppression method was proposed in this work. The method performs basis transformation of the material basis images based on a singular value decomposition. The noise variances of the original spectral CT images were incorporated in the matrix to be decomposed to ensure that the transformed basis images are statistically uncorrelated. Due to the difference in noise amplitudes in the transformed basis images, a selective filtering method was proposed with the low-noise transformed basis image as guidance. The method was evaluated using both numerical simulation and real clinical dual-energy CT data. Results demonstrated that compared with existing methods, the proposed method performs better in preserving the spatial resolution and the soft tissue contrast while suppressing the image noise. The proposed method is also computationally efficient and can realize real-time noise suppression for clinical spectral CT images.

Dual-energy head cone-beam CT using a dual-layer flat-panel detector: Hybrid material decomposition and a feasibility study

BACKGROUND

Flat panel detector (FPD) based cone-beam computed tomography (CT) has made tremendous progress in the last two decades, with many new and advanced medical and industrial applications keeping emerging from diagnostic imaging and image guidance for radiotherapy and interventional surgery. The current cone-beam CT (CBCT), however, is still suboptimal for head CT scan which requires a high standard of image quality. While the dual-layer FPD technology is under extensive development and is promising to further advance CBCT from qualitative anatomic imaging to quantitative dual-energy CT, its potential of enabling head CBCT applications has not yet been fully investigated.

PURPOSE

The relatively moderate energy separation from the dual-layer FPD and the overall low signal level especially at the bottom-layer detector, could raise significant challenges in performing high-quality dual-energy material decomposition (MD). In this work, we propose a hybrid, physics and model guided, MD algorithm that attempts to fully use the detected x-ray signals and prior-knowledge behind head CBCT using dual-layer FPD.

METHODS

Firstly, a regular projection-domain MD is performed as initial results of our approach and for comparison as conventional method. Secondly, based on the combined projection, a dual-layer multi-material spectral correction (dMMSC) is applied to generate beam hardening free images. Thirdly, the dMMSC corrected projections are adopted as a physics-model based guidance to generate a hybrid MD. A set of physics experiments including fan-beam scan and cone-beam scan using a head phantom and a Gammex Multi-Energy CT phantom are conducted to validate our proposed approach.

RESULTS

The combined reconstruction could reduce noise by about 10% with no visible resolution degradation. The fan-beam studies on the Gammex phantom demonstrated an improved MD performance, with the averaged iodine quantification error for the 5-15 mg/ml iodine inserts reduced from about 5.6% to 3.0% by the hybrid method. On fan-beam scan of the head phantom, our proposed hybrid MD could significantly reduce the streak artifacts, with CT number nonuniformity (NU) in the selected regions of interest (ROIs) reduced from 23 Hounsfield Units (HU) to 4.2 HU, and the corresponding noise suppressed from 31 to 6.5 HU. For cone-beam scan, after scatter correction (SC) and cone-beam artifact reduction (CBAR), our approach can also significantly improve image quality, with CT number NU in the selected ROI reduced from 24.2 to 6.6 HU and the noise level suppressed from 22.1 to 8.2 HU.

CONCLUSIONS

Our proposed physics and model guided hybrid MD for dual-layer FPD based head CBCT can significantly improve the robustness of MD and suppress the low-signal artifact. This preliminary feasibility study also demonstrated that the dual-layer FPD is promising to enable head CBCT spectral imaging.

Detection performance of pixelated lutetium-yttrium oxyorthosilicate (LYSO) scintillators for high-resolution photon-counting CT imaging

High-resolution photon-counting detector (PCD) computed tomography (CT) imaging is increasingly used for several applications. Recent technological advances in CT instrumentation have introduced various types of radiation detectors. Therefore, this work aims at evaluating the lutetium-yttrium oxyorthosilicate (LYSO) scintillator for use in PCD CT from a detector point of view. To do so, a mini-CT prototype was designed and constructed based on the pixelated LYSO blocks. The detector comprises four 10 × 10 linearly arranged LYSO blocks coupled with four position-sensitive photomultiplier tubes. The prototype utilizes a point gamma-ray source along with a cone-beam collimator. An in-home MATLAB-based data processing software package was also developed for storing the list-mode data, event positioning, and energy windowing. A set of experiments were conducted to assess the performance of the constructed energy-resolved LYSO:Ce detector for mini-CT imaging. The results show good crystal identification for all blocks with a maximum peak-to-valley ratio of 3.48. In addition, the findings confirm that the developed detector is position-sensitive. The 20% energy window provides an optimal performance by simultaneously providing good crystal identification and a scatter removal factor of 0.71. A 96% uniformity was also observed when the detector was irradiated with a uniform flood. The spatial resolution of the mini-CT prototype in the x- and y-directions was calculated to be 0.9 and 0.93 mm, respectively, corrected for a magnification factor of 2.5. It is concluded that the pixelated LYSO crystal is a promising alternative to the current detectors and would be the scintillator of choice for high-resolution PCD CT imaging tasks.

Multi-Material Decomposition for Single Energy CT Using Material Sparsity Constraint

Multi-material decomposition (MMD) decomposes CT images into basis material images, and is a promising technique in clinical diagnostic CT to identify material compositions within the human body. MMD could be implemented on measurements obtained from spectral CT protocol, although spectral CT data acquisition is not readily available in most clinical environments. MMD methods using single energy CT (SECT), broadly applied in radiological departments of most hospitals, have been proposed in the literature while challenged by the inferior decomposition accuracy and the limited number of material bases due to the constrained material information in the SECT measurement. In this paper, we propose an image-domain SECT MMD method using material sparsity as an assistance under the condition that each voxel of the CT image contains at most two different elemental materials. L₀ norm represents the material sparsity constraint (MSC) and is integrated into the decomposition objective function with a least-square data fidelity term, total variation term, and a sum-to-one constraint of material volume fractions. An accelerated primal-dual (APD) algorithm with line-search scheme is applied to solve the problem. The pixelwise direct inversion method with the two-material assumption (TMA) is applied to estimate the initials. We validate the proposed method on phantom and patient data. Compared with the TMA method, the proposed MSC method increases the volume fraction accuracy (VFA) from 92.0% to 98.5% in the phantom study. In the patient study, the calcification area can be clearly visualized in the virtual non-contrast image generated by the proposed method, and has a similar shape to that in the ground-truth contrast-free CT image. The high decomposition image quality from the proposed method substantially facilitates the SECT-based MMD clinical applications.

Image domain multi-material decomposition using single energy CT

Multi-material decomposition (MMD) technique decomposes the CT images into basis material images and has been promising in clinical practice for material composition quantification within the human body. MMD could be implemented using the image data acquired from spectral CT or its special case, dual-energy CT (DECT) while the spectral CT data acquisition usually requires a hardware modification. In this paper, we propose an image domain MMD method using single energy CT (SECT). The proposed objective function applies a least square data fidelity term to enforce the minimization between the linear combination of decomposed material image and the measured SECT image, and an edge-preserving (EP) regularization term to meet the piecewise constant property of the material image. We apply the optimization transfer principle to form a pixel-wise separable quadratic surrogate (PWSQS) function in each iteration to decrease the objective function. The pixelwise direct inversion method assisted by the two-material assumption (TMA) is applied to obtain a good initial value. The proposed method is evaluated using a digital phantom, a Catphan phantom and the clinical data. A low-pass filtration method is implemented for a comparison purpose. In the phantom study, the proposed TMA method achieves high volume fraction accuracy (VFA) of 79.64% and the proposed EP method further increases the VFA by 15.56% and decreases the decomposition standard deviation (STD) by 81.51% compared with the TMA method. At the comparable noise level, the proposed EP method increases spatial resolution by an overall factor of 1.01 when the modulation transfer function magnitude is decreased to 50% compared with the low-pass filtration method. In clinical data study, the virtual non-contrast image generated by the proposed method achieves the root-mean-squared-relative error of 2.93% compared with the contrast-free ground-truth image.

Reconstruction method for DECT with one half-scan plus a second limited-angle scan using prior knowledge of complementary support set (Pri-CSS)

Dual-energy computed tomography (DECT) has capability to improve material differentiation, but most scanning schemes require two sets of full-scan measurements at different x-ray spectra, limiting its application to imaging system with incomplete scan. In this study, using one half-scan and a second limited-angle scan, we propose a DECT reconstruction method by exploiting the consistent information of gradient images at high- and low-energy spectra, which relaxes the requirement of data acquisition of DECT. Based on the theory of sampling condition analysis, the complementary support set of gradient images plays an important role in image reconstruction because it constitutes the sufficient and necessary condition for accurate CT reconstruction. For DECT, the gradient images of high- and low-energy CT images ideally share the same complementary support set for the same object. Inspired by this idea, we extract the prior knowledge of complementary support set (Pri-CSS) from the gradient image of the first half-scan CT image to promote the second limited-angle CT reconstruction. Pri-CSS will be incorporated into total variation regularization model in the form of constrains. Alternative direction method is applied to iteratively solve the modified optimization model, thereby deriving the proposed algorithm to recover low-energy CT image from limited-angle measurements. The qualitative and quantitative experiments on digital and real data are performed to validate the proposed method. The results show that the proposed method outperforms its counterparts and achieve high reconstruction quality for the designed scanning configuration.

Statistical image-based material decomposition for triple-energy computed tomography using total variation regularization

BACKGROUND

Triple-energy computed tomography (TECT) can obtain x-ray attenuation measurements at three energy spectra, thereby allowing identification of different material compositions with same or very similar attenuation coefficients. This ability is known as material decomposition, which can decompose TECT images into different basis material image. However, the basis material image would be severely degraded when material decomposition is directly performed on the noisy TECT measurements using a matrix inversion method.

OBJECTIVE

To achieve high quality basis material image, we present a statistical image-based material decomposition method for TECT, which uses the penalized weighted least-squares (PWLS) criteria with total variation (TV) regularization (PWLS-TV).

METHODS

The weighted least-squares term involves the noise statistical properties of the material decomposition process, and the TV regularization penalizes differences between local neighboring pixels in a decomposed image, thereby contributing to improving the quality of the basis material image. Subsequently, an alternating optimization method is used to minimize the objective function.

RESULTS

The performance of PWLS-TV is quantitatively evaluated using digital and mouse thorax phantoms. The experimental results show that PWLS-TV material decomposition method can greatly improve the quality of decomposed basis material image compared to the quality of images obtained using the competing methods in terms of suppressing noise and preserving edge and fine structure details.

CONCLUSIONS

The PWLS-TV method can simultaneously perform noise reduction and material decomposition in one iterative step, and it results in a considerable improvement of basis material image quality.

Single-Scan Dual-Energy CT Using Primary Modulation

Compared with conventional computed tomography (CT), dual-energy CT (DECT) provides better material differentiation but requires projection data acquired with two different effective x-ray spectra, limiting DECT applications to specialized scanners. We propose a hardware-based method, known as PM-DECT, which utilizes primary beam modulation to enable single-scan DECT on a conventional CT scanner. PM-DECT inserts an attenuation sheet with a spatially varying pattern-primary beam modulator-between the x-ray source and imaged object. During a CT scan, the modulator selectively hardens the x-ray beam, thereby increasing the average photon energy at specific detector pixel locations. Thus, PM-DECT simultaneously acquires high and low energy data at each projection angle. From the sparse projection data, high and low energy CT images are jointly reconstructed and simultaneously decomposed into basis materials via an iterative CT reconstruction algorithm with gradient weighting and an improved version of similarity based regularization. Studies on Catphan 600 and anthropomorphic head phantoms demonstrate that PM-DECT retains a high level of spatial resolution compared with conventional CT scans. Electron density values calculated from decomposed images indicate a limited error of 1.12% for PM-DECT. Comparison against a two-scan DECT technique shows that PM-DECT's image reconstruction from sparse data sets contributes only 0.66% error. By granting the opportunity for high-quality single-scan DECT on conventional CT scanners via limited hardware modification, PM-DECT has the potential to liberate DECT from specialized scanners, extending clinical availability, and implementation.

Statistical image-domain multimaterial decomposition for dual-energy CT

PURPOSE

Dual-energy CT (DECT) enhances tissue characterization because of its basis material decomposition capability. In addition to conventional two-material decomposition from DECT measurements, multimaterial decomposition (MMD) is required in many clinical applications. To solve the ill-posed problem of reconstructing multi-material images from dual-energy measurements, additional constraints are incorporated into the formulation, including volume and mass conservation and the assumptions that there are at most three materials in each pixel and various material types among pixels. The recently proposed flexible image-domain MMD method decomposes pixels sequentially into multiple basis materials using a direct inversion scheme which leads to magnified noise in the material images. In this paper, we propose a statistical image-domain MMD method for DECT to suppress the noise.

METHODS

The proposed method applies penalized weighted least-square (PWLS) reconstruction with a negative log-likelihood term and edge-preserving regularization for each material. The statistical weight is determined by a data-based method accounting for the noise variance of high- and low-energy CT images. We apply the optimization transfer principles to design a serial of pixel-wise separable quadratic surrogates (PWSQS) functions which monotonically decrease the cost function. The separability in each pixel enables the simultaneous update of all pixels.

RESULTS

The proposed method is evaluated on a digital phantom, Catphan©600 phantom and three patients (pelvis, head, and thigh). We also implement the direct inversion and low-pass filtration methods for a comparison purpose. Compared with the direct inversion method, the proposed method reduces noise standard deviation (STD) in soft tissue by 95.35% in the digital phantom study, by 88.01% in the Catphan©600 phantom study, by 92.45% in the pelvis patient study, by 60.21% in the head patient study, and by 81.22% in the thigh patient study, respectively. The overall volume fraction accuracy is improved by around 6.85%. Compared with the low-pass filtration method, the root-mean-square percentage error (RMSE(%)) of electron densities in the Catphan©600 phantom is decreased by 20.89%. As modulation transfer function (MTF) magnitude decreased to 50%, the proposed method increases the spatial resolution by an overall factor of 1.64 on the digital phantom, and 2.16 on the Catphan©600 phantom. The overall volume fraction accuracy is increased by 6.15%.

CONCLUSIONS

We proposed a statistical image-domain MMD method using DECT measurements. The method successfully suppresses the magnified noise while faithfully retaining the quantification accuracy and anatomical structure in the decomposed material images. The proposed method is practical and promising for advanced clinical applications using DECT imaging.

Noise suppression for dual-energy CT via penalized weighted least-square optimization with similarity-based regularization

PURPOSE

Dual-energy CT (DECT) expands applications of CT imaging in its capability to decompose CT images into material images. However, decomposition via direct matrix inversion leads to large noise amplification and limits quantitative use of DECT. Their group has previously developed a noise suppression algorithm via penalized weighted least-square optimization with edge-preservation regularization (PWLS-EPR). In this paper, the authors improve method performance using the same framework of penalized weighted least-square optimization but with similarity-based regularization (PWLS-SBR), which substantially enhances the quality of decomposed images by retaining a more uniform noise power spectrum (NPS).

METHODS

The design of PWLS-SBR is based on the fact that averaging pixels of similar materials gives a low-noise image. For each pixel, the authors calculate the similarity to other pixels in its neighborhood by comparing CT values. Using an empirical Gaussian model, the authors assign high/low similarity value to one neighboring pixel if its CT value is close/far to the CT value of the pixel of interest. These similarity values are organized in matrix form, such that multiplication of the similarity matrix to the image vector reduces image noise. The similarity matrices are calculated on both high- and low-energy CT images and averaged. In PWLS-SBR, the authors include a regularization term to minimize the L-2 norm of the difference between the images without and with noise suppression via similarity matrix multiplication. By using all pixel information of the initial CT images rather than just those lying on or near edges, PWLS-SBR is superior to the previously developed PWLS-EPR, as supported by comparison studies on phantoms and a head-and-neck patient.

RESULTS

On the line-pair slice of the Catphan(©)600 phantom, PWLS-SBR outperforms PWLS-EPR and retains spatial resolution of 8 lp/cm, comparable to the original CT images, even at 90% reduction in noise standard deviation (STD). Similar performance on spatial resolution is observed on an anthropomorphic head phantom. In addition, results of PWLS-SBR show substantially improved image quality due to preservation of image NPS. On the Catphan(©)600 phantom, NPS using PWLS-SBR has a correlation of 93% with that via direct matrix inversion, while the correlation drops to -52% for PWLS-EPR. Electron density measurement studies indicate high accuracy of PWLS-SBR. On seven different materials, the measured electron densities calculated from the decomposed material images using PWLS-SBR have a root-mean-square error (RMSE) of 1.20%, while the results of PWLS-EPR have a RMSE of 2.21%. In the study on a head-and-neck patient, PWLS-SBR is shown to reduce noise STD by a factor of 3 on material images with image qualities comparable to CT images, whereas fine structures are lost in the PWLS-EPR result. Additionally, PWLS-SBR better preserves low contrast on the tissue image.

CONCLUSIONS

The authors propose improvements to the regularization term of an optimization framework which performs iterative image-domain decomposition for DECT with noise suppression. The regularization term avoids calculation of image gradient and is based on pixel similarity. The proposed method not only achieves a high decomposition accuracy, but also improves over the previous algorithm on NPS as well as spatial resolution.

Dual energy CT with one full scan and a second sparse-view scan using structure preserving iterative reconstruction (SPIR)

Conventional dual-energy CT (DECT) reconstruction requires two full-size projection datasets with two different energy spectra. In this study, we propose an iterative algorithm to enable a new data acquisition scheme which requires one full scan and a second sparse-view scan for potential reduction in imaging dose and engineering cost of DECT. A bilateral filter is calculated as a similarity matrix from the first full-scan CT image to quantify the similarity between any two pixels, which is assumed unchanged on a second CT image since DECT scans are performed on the same object. The second CT image from reduced projections is reconstructed by an iterative algorithm which updates the image by minimizing the total variation of the difference between the image and its filtered image by the similarity matrix under data fidelity constraint. As the redundant structural information of the two CT images is contained in the similarity matrix for CT reconstruction, we refer to the algorithm as structure preserving iterative reconstruction (SPIR). The proposed method is evaluated on both digital and physical phantoms, and is compared with the filtered-backprojection (FBP) method, the conventional total-variation-regularization-based algorithm (TVR) and prior-image-constrained-compressed-sensing (PICCS). SPIR with a second 10-view scan reduces the image noise STD by a factor of one order of magnitude with same spatial resolution as full-view FBP image. SPIR substantially improves over TVR on the reconstruction accuracy of a 10-view scan by decreasing the reconstruction error from 6.18% to 1.33%, and outperforms TVR at 50 and 20-view scans on spatial resolution with a higher frequency at the modulation transfer function value of 10% by an average factor of 4. Compared with the 20-view scan PICCS result, the SPIR image has 7 times lower noise STD with similar spatial resolution. The electron density map obtained from the SPIR-based DECT images with a second 10-view scan has an average error of less than 1%.

A general framework of noise suppression in material decomposition for dual-energy CT

PURPOSE

As a general problem of dual-energy CT (DECT), noise amplification in material decomposition severely reduces the signal-to-noise ratio on the decomposed images compared to that on the original CT images. In this work, the authors propose a general framework of noise suppression in material decomposition for DECT. The method is based on an iterative algorithm recently developed in their group for image-domain decomposition of DECT, with an extension to include nonlinear decomposition models. The generalized framework of iterative DECT decomposition enables beam-hardening correction with simultaneous noise suppression, which improves the clinical benefits of DECT.

METHODS

The authors propose to suppress noise on the decomposed images of DECT using convex optimization, which is formulated in the form of least-squares estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, the authors include the inverse of the estimated variance-covariance matrix of the decomposed images as the penalty weight in the least-squares term. Analytical formulas are derived to compute the variance-covariance matrix for decomposed images with general-form numerical or analytical decomposition. As a demonstration, the authors implement the proposed algorithm on phantom data using an empirical polynomial function of decomposition measured on a calibration scan. The polynomial coefficients are determined from the projection data acquired on a wedge phantom, and the signal decomposition is performed in the projection domain.

RESULTS

On the Catphan(®)600 phantom, the proposed noise suppression method reduces the average noise standard deviation of basis material images by one to two orders of magnitude, with a superior performance on spatial resolution as shown in comparisons of line-pair images and modulation transfer function measurements. On the synthesized monoenergetic CT images, the noise standard deviation is reduced by a factor of 2-3. By using nonlinear decomposition on projections, the authors' method effectively suppresses the streaking artifacts of beam hardening and obtains more uniform images than their previous approach based on a linear model. Similar performance of noise suppression is observed in the results of an anthropomorphic head phantom and a pediatric chest phantom generated by the proposed method. With beam-hardening correction enabled by their approach, the image spatial nonuniformity on the head phantom is reduced from around 10% on the original CT images to 4.9% on the synthesized monoenergetic CT image. On the pediatric chest phantom, their method suppresses image noise standard deviation by a factor of around 7.5, and compared with linear decomposition, it reduces the estimation error of electron densities from 33.3% to 8.6%.

CONCLUSIONS

The authors propose a general framework of noise suppression in material decomposition for DECT. Phantom studies have shown the proposed method improves the image uniformity and the accuracy of electron density measurements by effective beam-hardening correction and reduces noise level without noticeable resolution loss.

Noise Suppression for Dual-Energy CT Through Entropy Minimization

In dual energy CT (DECT), noise amplification during signal decomposition significantly limits the utility of basis material images. Since clinically relevant objects typically contain a limited number of different materials, we propose an Image-domain Decomposition method through Entropy Minimization (IDEM) for noise suppression in DECT. Pixels of decomposed images are first linearly transformed into 2D clusters of data points, which are highly asymmetric due to strong signal correlation. An optimal axis is identified in the 2D space via numerical search such that the projection of data clusters onto the axis has minimum entropy. Noise suppression is performed on each image pixel by estimating the center-of-mass value of each data cluster along the direction perpendicular to the projection axis. The IDEM method is distinct from other noise suppression techniques in that it does not suppress pixel noise by reducing spatial variation between neighboring pixels. As supported by studies on Catphan©600 and anthropomorphic head phantoms, this feature endows our algorithm with a unique capability of reducing noise standard deviation on DECT decomposed images by approximately one order of magnitude while preserving spatial resolution and image noise power spectra (NPS). Compared with a filtering method and recently developed iterative method at the same level of noise suppression, the IDEM algorithm obtains high-resolution images with less artifacts. It also maintains accuracy of electron density measurements with less than 2% bias error. The IDEM method effectively suppresses noise of DECT for quantitative use, with appealing features on preservation of image spatial resolution and NPS.

Combined iterative reconstruction and image-domain decomposition for dual energy CT using total-variation regularization

PURPOSE

Dual-energy CT (DECT) is being increasingly used for its capability of material decomposition and energy-selective imaging. A generic problem of DECT, however, is that the decomposition process is unstable in the sense that the relative magnitude of decomposed signals is reduced due to signal cancellation while the image noise is accumulating from the two CT images of independent scans. Direct image decomposition, therefore, leads to severe degradation of signal-to-noise ratio on the resultant images. Existing noise suppression techniques are typically implemented in DECT with the procedures of reconstruction and decomposition performed independently, which do not explore the statistical properties of decomposed images during the reconstruction for noise reduction. In this work, the authors propose an iterative approach that combines the reconstruction and the signal decomposition procedures to minimize the DECT image noise without noticeable loss of resolution.

METHODS

The proposed algorithm is formulated as an optimization problem, which balances the data fidelity and total variation of decomposed images in one framework, and the decomposition step is carried out iteratively together with reconstruction. The noise in the CT images from the proposed algorithm becomes well correlated even though the noise of the raw projections is independent on the two CT scans. Due to this feature, the proposed algorithm avoids noise accumulation during the decomposition process. The authors evaluate the method performance on noise suppression and spatial resolution using phantom studies and compare the algorithm with conventional denoising approaches as well as combined iterative reconstruction methods with different forms of regularization.

RESULTS

On the Catphan©600 phantom, the proposed method outperforms the existing denoising methods on preserving spatial resolution at the same level of noise suppression, i.e., a reduction of noise standard deviation by one order of magnitude. This improvement is mainly attributed to the high noise correlation in the CT images reconstructed by the proposed algorithm. Iterative reconstruction using different regularization, including quadratic orq-generalized Gaussian Markov random field regularization, achieves similar noise suppression from high noise correlation. However, the proposed TV regularization obtains a better edge preserving performance. Studies of electron density measurement also show that our method reduces the average estimation error from 9.5% to 7.1%. On the anthropomorphic head phantom, the proposed method suppresses the noise standard deviation of the decomposed images by a factor of ∼14 without blurring the fine structures in the sinus area.

CONCLUSIONS

The authors propose a practical method for DECT imaging reconstruction, which combines the image reconstruction and material decomposition into one optimization framework. Compared to the existing approaches, our method achieves a superior performance on DECT imaging with respect to decomposition accuracy, noise reduction, and spatial resolution.

Iterative image-domain decomposition for dual-energy CT

PURPOSE

Dual energy CT (DECT) imaging plays an important role in advanced imaging applications due to its capability of material decomposition. Direct decomposition via matrix inversion suffers from significant degradation of image signal-to-noise ratios, which reduces clinical values of DECT. Existing denoising algorithms achieve suboptimal performance since they suppress image noise either before or after the decomposition and do not fully explore the noise statistical properties of the decomposition process. In this work, the authors propose an iterative image-domain decomposition method for noise suppression in DECT, using the full variance-covariance matrix of the decomposed images.

METHODS

The proposed algorithm is formulated in the form of least-square estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, the authors include the inverse of the estimated variance-covariance matrix of the decomposed images as the penalty weight in the least-square term. The regularization term enforces the image smoothness by calculating the square sum of neighboring pixel value differences. To retain the boundary sharpness of the decomposed images, the authors detect the edges in the CT images before decomposition. These edge pixels have small weights in the calculation of the regularization term. Distinct from the existing denoising algorithms applied on the images before or after decomposition, the method has an iterative process for noise suppression, with decomposition performed in each iteration. The authors implement the proposed algorithm using a standard conjugate gradient algorithm. The method performance is evaluated using an evaluation phantom (Catphan©600) and an anthropomorphic head phantom. The results are compared with those generated using direct matrix inversion with no noise suppression, a denoising method applied on the decomposed images, and an existing algorithm with similar formulation as the proposed method but with an edge-preserving regularization term.

RESULTS

On the Catphan phantom, the method maintains the same spatial resolution on the decomposed images as that of the CT images before decomposition (8 pairs/cm) while significantly reducing their noise standard deviation. Compared to that obtained by the direct matrix inversion, the noise standard deviation in the images decomposed by the proposed algorithm is reduced by over 98%. Without considering the noise correlation properties in the formulation, the denoising scheme degrades the spatial resolution to 6 pairs/cm for the same level of noise suppression. Compared to the edge-preserving algorithm, the method achieves better low-contrast detectability. A quantitative study is performed on the contrast-rod slice of Catphan phantom. The proposed method achieves lower electron density measurement error as compared to that by the direct matrix inversion, and significantly reduces the error variation by over 97%. On the head phantom, the method reduces the noise standard deviation of decomposed images by over 97% without blurring the sinus structures.

CONCLUSIONS

The authors propose an iterative image-domain decomposition method for DECT. The method combines noise suppression and material decomposition into an iterative process and achieves both goals simultaneously. By exploring the full variance-covariance properties of the decomposed images and utilizing the edge predetection, the proposed algorithm shows superior performance on noise suppression with high image spatial resolution and low-contrast detectability.

Quantitative material decomposition using spectral computed tomography with an energy-resolved photon-counting detector

Dual-energy computed tomography (CT) techniques have been used to decompose materials and characterize tissues according to their physical and chemical compositions. However, these techniques are hampered by the limitations of conventional x-ray detectors operated in charge integrating mode. Energy-resolved photon-counting detectors provide spectral information from polychromatic x-rays using multiple energy thresholds. These detectors allow simultaneous acquisition of data in different energy ranges without spectral overlap, resulting in more efficient material decomposition and quantification for dual-energy CT. In this study, a pre-reconstruction dual-energy CT technique based on volume conservation was proposed for three-material decomposition. The technique was combined with iterative reconstruction algorithms by using a ray-driven projector in order to improve the quality of decomposition images and reduce radiation dose. A spectral CT system equipped with a CZT-based photon-counting detector was used to implement the proposed dual-energy CT technique. We obtained dual-energy images of calibration and three-material phantoms consisting of low atomic number materials from the optimal energy bins determined by Monte Carlo simulations. The material decomposition process was accomplished by both the proposed and post-reconstruction dual-energy CT techniques. Linear regression and normalized root-mean-square error (NRMSE) analyses were performed to evaluate the quantitative accuracy of decomposition images. The calibration accuracy of the proposed dual-energy CT technique was higher than that of the post-reconstruction dual-energy CT technique, with fitted slopes of 0.97-1.01 and NRMSEs of 0.20-4.50% for all basis materials. In the three-material phantom study, the proposed dual-energy CT technique decreased the NRMSEs of measured volume fractions by factors of 0.17-0.28 compared to the post-reconstruction dual-energy CT technique. It was concluded that the proposed dual-energy CT technique can potentially be used to decompose mixtures into basis materials and characterize tissues according to their composition.

Kilovoltage dependence of the attenuation of a potassium iodide/water solution on CT: presentation of a computer model implementing polychromatic character of the X-ray photon beam

PURPOSE

To describe a program that simulates a computed tomographic scan with the polychromatic aspect of the output of the x-ray tube implemented. This program can be used in the study of the attenuation of different solutions/solutes. These results can subsequently guide the radiologist to obtain a satisfying contrast enhancement at lower tube voltages and eventually lower contrast volumes.

MATERIALS AND METHODS

A Matlab program was written to simulate a computed tomographic scan. The spectrum of the x-ray tube at different kilovoltages was generated with another program (XOP) and used as input. Beam-hardening correction and zero padding were added. The results were validated with attenuation measurements of a corresponding potassium iodide solution in water.

RESULTS

There was a good agreement between the calculated and measured attenuations; the calculated results matched with the measured values and fell within a 5% deviation.

CONCLUSION

It is possible to simulate correctly the attenuation of a potassium iodide in water solution in silico. This can be helpful to determine kilovoltages, administered contrast medium volumes and concentrations to reduce the irradiation of the patients, and obtain equally good contrast enhancement on a basis other than empirical.

Correlation between CT attenuation value and iodine concentration in vitro: discrepancy between gemstone spectral imaging on single-source dual-energy CT and traditional polychromatic X-ray imaging

INTRODUCTION

To assess the relation between CT attenuation value and iodine concentration in vitro, using gemstone spectral imaging (GSI) with single-source dual-energy CT and traditional polychromatic X-ray imaging (TPXI), respectively.

METHODS

A polypropylene phantom with eight test tubes in which iodine concentrations of solution were 0.4, 0.7, 2, 5, 10, 20, 30 and 50 mg/mL underwent GSI and traditional polychromatic X-ray scans (80, 100, 120 and 140 kV(p)), using single-source dual-energy spectral CT (Discovery CT750HD; GE Healthcare Technologies, Milwaukee, WI, USA) at the same tube speed of 0.8 s/rotation. All spectral imaging data were analysed with GSI viewer to obtain monochromatic images (50-140 keV, interval of 10 keV). Computed tomography attenuation values of iodine solution were measured with the same size of regions of interest and at the exact same level for both monochromatic and polychromatic images. The relation between CT attenuation value and iodine concentration was examined.

RESULTS

A linear correlation was found between CT attenuation value and iodine concentration for both monochromatic and polychromatic images. Moreover, the fitting coefficients for CT attenuation values and iodine concentrations were closer to one with GSI (r(2) = 0.99824-0.99996) than that with TPXI (r(2) = 0.99640-0.99736).

CONCLUSIONS

Owing to the better correlation coefficients between CT attenuation value and iodine concentration, GSI may be a preferred method for quantitative measurement compared with TPXI.

Imaging of nanoparticles with dual-energy computed tomography

A simulation study was performed to determine the feasibility and performance of imaging nanoparticles as contrast agents in dual-energy computed tomography. An analytical simulation model was used to model the relevant signal-to-noise ratio (SNR) in dual-energy imaging for the specific case of a three-material patient phantom consisting of water, calcium hydroxyapatite and contrast agent. Elemental gold and iodine were both considered as contrast agents. Simulations were performed for a range of monoenergetic (20-150 keV) and polyenergetic (20-150 kVp) beam spectra. A reference configuration was defined with beam energies of 80 and 140 kVp to match current clinical practice. The effect of adding a silver filter to the high-energy beam was also studied. A figure of merit (FOM), which normalized the dual-energy SNR to the square root of the patient integral dose, was calculated for all cases. The units of the FOM were keV(-1/2). A simple Rose model of detectability was used to estimate the minimum concentration of either elements needed to be detected (SNR > 5). For monoenergetic beams, the peak FOM of gold was 6.4 × 10(-6) keV(-1/2), while the peak FOM of iodine was 3.1 × 10(-6) keV(-1/2), a factor of approximately 2 greater for gold. For polyenergetic spectra, at the reference energies of 80 and 140 kVp, the FOM for gold and iodine was 1.65 × 10(-6) and 5.0 × 10(-7) keV(-1/2), respectively, a factor of approximately 3.3 greater. Also at these energies, the minimum detectable concentration of gold was estimated to be 58.5 mg mL(-1), while iodine was estimated to be 117.5 mg mL(-1). The results suggest that the imaging of a gold nanoparticle contrast agent is well suited to current conditions used in clinical imaging. The addition of a silver filter of 800 µm further increased the image quality of the gold signal by approximately 50% for the same absorbed dose to the patient.

Image quality optimization and evaluation of linearly mixed images in dual-source, dual-energy CT

In dual-source dual-energy CT, the images reconstructed from the low- and high-energy scans (typically at 80 and 140 kV, respectively) can be mixed together to provide a single set of nonmaterial-specific images for the purpose of routine diagnostic interpretation. Different from the material-specific information that may be obtained from the dual-energy scan data, the mixed images are created with the purpose of providing the interpreting physician a single set of images that have an appearance similar to that in single-energy images acquired at the same total radiation dose. In this work, the authors used a phantom study to evaluate the image quality of linearly mixed images in comparison to single-energy CT images, assuming the same total radiation dose and taking into account the effect of patient size and the dose partitioning between the low-and high-energy scans. The authors first developed a method to optimize the quality of the linearly mixed images such that the single-energy image quality was compared to the best-case image quality of the dual-energy mixed images. Compared to 80 kV single-energy images for the same radiation dose, the iodine CNR in dual-energy mixed images was worse for smaller phantom sizes. However, similar noise and similar or improved iodine CNR relative to 120 kV images could be achieved for dual-energy mixed images using the same total radiation dose over a wide range of patient sizes (up to 45 cm lateral thorax dimension). Thus, for adult CT practices, which primarily use 120 kV scanning, the use of dual-energy CT for the purpose of material-specific imaging can also produce a set of non-material-specific images for routine diagnostic interpretation that are of similar or improved quality relative to single-energy 120 kV scans.