Suppression of metal artifacts in CT using a reconstruction procedure that combines MAP and projection completion

Deep-Learning-Based Metal Artefact Reduction With Unsupervised Domain Adaptation Regularization for Practical CT Images

CT metal artefact reduction (MAR) methods based on supervised deep learning are often troubled by domain gap between simulated training dataset and real-application dataset, i.e., methods trained on simulation cannot generalize well to practical data. Unsupervised MAR methods can be trained directly on practical data, but they learn MAR with indirect metrics and often perform unsatisfactorily. To tackle the domain gap problem, we propose a novel MAR method called UDAMAR based on unsupervised domain adaptation (UDA). Specifically, we introduce a UDA regularization loss into a typical image-domain supervised MAR method, which mitigates the domain discrepancy between simulated and practical artefacts by feature-space alignment. Our adversarial-based UDA focuses on a low-level feature space where the domain difference of metal artefacts mainly lies. UDAMAR can simultaneously learn MAR from simulated data with known labels and extract critical information from unlabeled practical data. Experiments on both clinical dental and torso datasets show the superiority of UDAMAR by outperforming its supervised backbone and two state-of-the-art unsupervised methods. We carefully analyze UDAMAR by both experiments on simulated metal artefacts and various ablation studies. On simulation, its close performance to the supervised methods and advantages over the unsupervised methods justify its efficacy. Ablation studies on the influence from the weight of UDA regularization loss, UDA feature layers, and the amount of practical data used for training further demonstrate the robustness of UDAMAR. UDAMAR provides a simple and clean design and is easy to implement. These advantages make it a very feasible solution for practical CT MAR.

InDuDoNet+: A deep unfolding dual domain network for metal artifact reduction in CT images

During the computed tomography (CT) imaging process, metallic implants within patients often cause harmful artifacts, which adversely degrade the visual quality of reconstructed CT images and negatively affect the subsequent clinical diagnosis. For the metal artifact reduction (MAR) task, current deep learning based methods have achieved promising performance. However, most of them share two main common limitations: (1) the CT physical imaging geometry constraint is not comprehensively incorporated into deep network structures; (2) the entire framework has weak interpretability for the specific MAR task; hence, the role of each network module is difficult to be evaluated. To alleviate these issues, in the paper, we construct a novel deep unfolding dual domain network, termed InDuDoNet+, into which CT imaging process is finely embedded. Concretely, we derive a joint spatial and Radon domain reconstruction model and propose an optimization algorithm with only simple operators for solving it. By unfolding the iterative steps involved in the proposed algorithm into the corresponding network modules, we easily build the InDuDoNet+ with clear interpretability. Furthermore, we analyze the CT values among different tissues, and merge the prior observations into a prior network for our InDuDoNet+, which significantly improve its generalization performance. Comprehensive experiments on synthesized data and clinical data substantiate the superiority of the proposed methods as well as the superior generalization performance beyond the current state-of-the-art (SOTA) MAR methods. Code is available at https://github.com/hongwang01/InDuDoNet_plus.

Physics-informed sinogram completion for metal artifact reduction in CT imaging

Objective.Metal artifacts in the computed tomography (CT) imaging are unavoidably adverse to the clinical diagnosis and treatment outcomes. Most metal artifact reduction (MAR) methods easily result in the over-smoothing problem and loss of structure details near the metal implants, especially for these metal implants with irregular elongated shapes. To address this problem, we present the physics-informed sinogram completion (PISC) method for MAR in CT imaging, to reduce metal artifacts and recover more structural textures.Approach.Specifically, the original uncorrected sinogram is firstly completed by the normalized linear interpolation algorithm to reduce metal artifacts. Simultaneously, the uncorrected sinogram is also corrected based on the beam-hardening correction physical model, to recover the latent structure information in metal trajectory region by leveraging the attenuation characteristics of different materials. Both corrected sinograms are fused with the pixel-wise adaptive weights, which are manually designed according to the shape and material information of metal implants. To furtherly reduce artifacts and improve the CT image quality, a post-processing frequency split algorithm is adopted to yield the final corrected CT image after reconstructing the fused sinogram.Main results.We qualitatively and quantitatively evaluated the presented PISC method on two simulated datasets and three real datasets. All results demonstrate that the presented PISC method can effectively correct the metal implants with various shapes and materials, in terms of artifact suppression and structure preservation.Significance.We proposed a sinogram-domain MAR method to compensate for the over-smoothing problem existing in most MAR methods by taking advantage of the physical prior knowledge, which has the potential to improve the performance of the deep learning based MAR approaches.

Metal Artifact Reduction in Dental CBCT Images Using Direct Sinogram Correction Combined with Metal Path-Length Weighting

Metal artifacts in dental computed tomography (CT) images, caused by highly X-ray absorbing objects, such as dental implants or crowns, often more severely compromise image readability than in medical CT images. Since lower tube voltages are used for dental CTs in spite of the more frequent presence of metallic objects in the patient, metal artifacts appear more severely in dental CT images, and the artifacts often persist even after metal artifact correction. The direct sinogram correction (DSC) method, which directly corrects the sinogram using the mapping function derived by minimizing the sinogram inconsistency, works well in the case of mild metal artifacts, but it often fails to correct severe metal artifacts. We propose a modified DSC method to reduce severe metal artifacts, and we have tested it on human dental images. We first segment the metallic objects in the CT image, and then we forward-project the segmented metal mask to identify the metal traces in the projection data with computing the metal path length for the rays penetrating the metal mask. In the sinogram correction with the DSC mapping function, we apply the weighting proportional to the metal path length. We have applied the proposed method to the phantom and patient images taken at the X-ray tube voltage of 90 kVp. We observed that the proposed method outperforms the original DSC method when metal artifacts were severe. However, we need further extensive studies to verify the proposed method for various CT scan conditions with many more patient images.

Metal artifact correction in photon-counting detector computed tomography: metal trace replacement using high-energy data

BACKGROUND

Metal artifacts have been an outstanding issue in computed tomography (CT) since its first uses in the clinic and continue to interfere. Metal artifact reduction (MAR) methods continue to be proposed and photon-counting detectors (PCDs) have recently been the subject of research toward this purpose. PCDs offer the ability to distinguish the energy of incident x-rays and sort them in a set number of energy bins. High-energy data captured using PCDs have been shown to reduce metal artifacts in reconstructions due to reduced beam hardening.

PURPOSE

High-energy reconstructions using PCD-CT have their drawbacks, such as reduced image contrast and increased noise. Here, we demonstrate a MAR algorithm, trace replacement MAR (TRMAR), in which the data corrupted by metal artifacts in full energy spectrum projections are corrected using the high-energy data captured during the same scan. The resulting reconstructions offer similar MAR to that seen in high-energy reconstructions, but with improved image quality.

METHODS

Experimental data were collected using a bench-top PCD-CT system with a cadmium zinc telluride PCD. Simulations were performed to determine the optimal high-energy threshold and to test TRMAR in simulations using the XCAT phantom and a biological sample. For experiments a 100-mm diameter cylindrical phantom containing vials of water, two screws, various densities of Ca(ClO₄ )₂ , and a spatial resolution phantom was imaged with and without the screws. The screws were segmented in the initial reconstruction and forward projected to identify them in the sinogram space in order to perform TRMAR. The resulting reconstructions were compared to the control and to reconstructions corrected using normalized metal artifact reduction (NMAR). Additionally, a beef short rib was imaged with and without metal to provide a more realistic phantom.

RESULTS

XCAT simulations showed a reduction in the streak artifact from -978 HU in uncorrected images to -10 HU with TRMAR. The magnitude of the metal artifact in uncorrected images of the 100-mm phantom was -442 HU, compared to the desired -81 HU with no metal. TRMAR reduced the magnitude of the artifact to -142 HU, with NMAR reducing the magnitude to -96 HU. Relative image noise was reduced from 176% in the high-energy image to 56% using TRMAR. Density quantification was better with NMAR, with the Ca(ClO₄ )₂ vial affected most by metal artifacts showing 0.8% error compared to 2.1% with TRMAR. Small features were preserved to a greater extent with TRMAR, with the limiting spatial frequency at 20% of the MTF fully maintained at 1.31 lp/mm, while with NMAR it was reduced to 1.22 lp/mm. Images of the beef short rib showed better delineation of the shape of the metal using TRMAR.

CONCLUSIONS

NMAR offers slightly better performance compared to TRMAR in streak reduction and image quality metrics. However, TRMAR is less susceptible to metal segmentation errors and can closely approximate the reduction in the streak metal artifact seen in NMAR at 1/3 the computation time. With the recent introduction of PCD-CT into the clinic, TRMAR offers notable potential for fast, effective MAR.

Metal artifact reduction in kV CT images throughout two-step sequential deep convolutional neural networks by combining multi-modal imaging (MARTIAN)

This work attempted to construct a new metal artifact reduction (MAR) framework in kilo-voltage (kV) computed tomography (CT) images by combining (1) deep learning and (2) multi-modal imaging, defined as MARTIAN (Metal Artifact Reduction throughout Two-step sequentIAl deep convolutional neural Networks). Most CNNs under supervised learning require artifact-free images to artifact-contaminated images for artifact correction. Mega-voltage (MV) CT is insensitive to metal artifacts, unlike kV CT due to different physical characteristics, which can facilitate the generation of artifact-free synthetic kV CT images throughout the first network (Network 1). The pairs of true kV CT and artifact-free kV CT images after post-processing constructed a subsequent network (Network 2) to conduct the actual MAR process. The proposed framework was implemented by GAN from 90 scans for head-and-neck and brain radiotherapy and validated with 10 independent cases against commercial MAR software. The artifact-free kV CT images following Network 1 and post-processing led to structural similarity (SSIM) of 0.997, and mean-absolute-error (MAE) of 10.2 HU, relative to true kV CT. Network 2 in charge of actual MAR successfully suppressed metal artifacts, relative to commercial MAR, while retaining the detailed imaging information, yielding the SSIM of 0.995 against 0.997 from the commercial MAR.

One-step iterative reconstruction approach based on eigentissue decomposition for spectral photon-counting computed tomography

Purpose: We propose a one-step tissue characterization method for spectral photon-counting computed tomography (SPCCT) using eigentissue decomposition (ETD), tailored for highly accurate human tissue characterization in radiotherapy. Methods: The approach combines a Poisson likelihood, a spatial prior, and a quantitative prior constraining eigentissue fractions based on expected values for tabulated tissues. There are two regularization parameters: α for the quantitative prior, and β for the spatial prior. The approach is validated in a realistic simulation environment for SPCCT. The impact of α and β is evaluated on a virtual phantom. The framework is tested on a virtual patient and compared with two sinogram-based two-step methods [using respectively filtered backprojection (FBP) and an iterative method for the second step] and a post-reconstruction approach with the same quantitative prior. All methods use ETD. Results: Optimal performance with respect to bias or RMSE is achieved with different combinations of α and β on the cylindrical phantom. Evaluated in tissues of the virtual patient, the one-step framework outperforms two-step and post-reconstruction approaches to quantify proton-stopping power (SPR). The mean absolute bias on the SPR is 0.6% (two-step FBP), 0.6% (two-step iterative), 0.6% (post-reconstruction), and 0.2% (one-step optimized for low bias). Following the same order, the RMSE on the SPR is 13.3%, 2.5%, 3.2%, and 1.5%. Conclusions: Accurate and precise characterization with ETD can be achieved with noisy SPCCT data without the need to rely on post-reconstruction methods. The one-step framework is more accurate and precise than two-step methods for human tissue characterization.

[An adaptive CT metal artifact reduction algorithm that combines projection interpolation and physical correction]

OBJECTIVE

To propose an adaptive weighted CT metal artifact reduce algorithm that combines projection interpolation and physical correction.

METHODS

A normalized metal projection interpolation algorithm was used to obtain the initial corrected projection data. A metal physical correction model was then introduced to obtain the physically corrected projection data. To verify the effectiveness of the method, we conducted experiments using simulation data and clinical data. For the simulation data, the quantitative indicators PSNR and SSIM were used for evaluation, while for the clinical data, the resultant images were evaluated by imaging experts to compare the artifact-reducing performance of different methods.

RESULTS

For the simulation data, the proposed method improved the PSNR value by at least 0.2 dB and resulted in the highest SSIM value among the methods for comparison. The experiment with the clinical data showed that the imaging experts gave the highest scores of 3.616±0.338 (in a 5-point scale) to the images processed using the proposed method, which had significant better artifact-reducing performance than the other methods (P < 0.001).

CONCLUSION

The metal artifact reduction algorithm proposed herein can effectively reduce metal artifacts while preserving the tissue structure information and reducing the generation of new artifacts.

Multi-modal feature-fusion for CT metal artifact reduction using edge-enhanced generative adversarial networks

Computed Tomography (CT) imaging is one of the most widely-used and cost-effective technology for organ screening and diseases diagnosis. Because of existence of metallic implants in some patients, the CT images acquired from these patients are often corrupted by undesirable metal artifacts, which causes severe problem of metal artifact. Although there have been proposed many methods to reduce metal artifact, reduction is still challenging and inadequate, and results are suffering from symptom variance, second artifact and poor subjective evaluation. To address these problems, we propose a novel metal artifact reduction method based on generative adversarial networks to simultaneously reduce metal artifacts and enhance texture structure of corrected CT images. Specifically, we firstly incorporate interactive information (text) and imaging CT (image) into a comprehensive feature to yield multi-modal feature-fusion representation, which overcomes the representative ability limitation of single-modal data. The incorporation of interaction information constrains the feature generation to ensure symptom consistency between corrected and target CT. Then, we design an edge-enhance sub-network to avoid second artifact and suppress noise. Besides, we invite three professional physicians to evaluate corrected CT image subjectively. In this paper, We achieved average increment of 11.3% PSNR and 12.1% SSIM on DeepLesion dataset. The subjective evaluations by physicians show that ours outperforms over 6.3%, 7.1%, 5.50% and 6.9% in term of sharpness, resolution, invariance and acceptability, respectively. Our proposed method can achieve high-quality metal artifact reduction results.

DICDNet: Deep Interpretable Convolutional Dictionary Network for Metal Artifact Reduction in CT Images

Computed tomography (CT) images are often impaired by unfavorable artifacts caused by metallic implants within patients, which would adversely affect the subsequent clinical diagnosis and treatment. Although the existing deep-learning-based approaches have achieved promising success on metal artifact reduction (MAR) for CT images, most of them treated the task as a general image restoration problem and utilized off-the-shelf network modules for image quality enhancement. Hence, such frameworks always suffer from lack of sufficient model interpretability for the specific task. Besides, the existing MAR techniques largely neglect the intrinsic prior knowledge underlying metal-corrupted CT images which is beneficial for the MAR performance improvement. In this paper, we specifically propose a deep interpretable convolutional dictionary network (DICDNet) for the MAR task. Particularly, we first explore that the metal artifacts always present non-local streaking and star-shape patterns in CT images. Based on such observations, a convolutional dictionary model is deployed to encode the metal artifacts. To solve the model, we propose a novel optimization algorithm based on the proximal gradient technique. With only simple operators, the iterative steps of the proposed algorithm can be easily unfolded into corresponding network modules with specific physical meanings. Comprehensive experiments on synthesized and clinical datasets substantiate the effectiveness of the proposed DICDNet as well as its superior interpretability, compared to current state-of-the-art MAR methods. Code is available at https://github.com/hongwang01/DICDNet.

Inpainting-filtering for metal artifact reduction (IMIF-MAR) in computed tomography

The reduction of metal artifacts remains a challenge in computed tomography because they decrease image quality, and consequently might affect the medical diagnosis. The objective of this study is to present a novel method to correct metal artifacts based solely on the CT-slices. The proposed method consists of four steps. First, metal implants in the original CT-slice are segmented using an entropy based method, producing a metal image. Second, a prior image is acquired using three transformations: Gaussian filter, Parisotto and Schoenlieb inpainting method with the Mumford-Shah image model and L0 Gradient Minimization method (L0GM). Next, based on the projections from the original CT-slice, prior image and metal image, the sinogram is corrected in the traces affected by metal in the process called normalization and denormalization. Finally, the reconstructed image is obtained by FBP and a Nonlocal Means (NLM) filtering. The efficacy of the algorithm is evaluated by comparing five image quality metrics of the images and by inspecting regions of interest (ROI). Phantom data as well as clinical datasets are included. The proposed method is compared with three established metal artifact reduction (MAR) methods. The results from a phantom and clinical dataset show the visible reduction of artifacts. The conclusion is that IMIF-MAR method can reduce streak metal artifacts effectively and avoid new artifacts around metal implants, while preserving the anatomical structures. Considering both clinical and phantom studies, the proposed MAR algorithm improves the quality of clinical images affected by metal artifacts, and could be integrated in clinical setting.

Comprehensive assessments of the open mouth dynamic maneuver and metal artifact reduction algorithm on computed tomography images of the oral cavity and oropharynx

Objectives

To determine the optimal utility of the open mouth maneuver and Metal Artifact Reduction for the Orthopedic Implants (O-MAR) technique for CT of the oral cavity and oropharynx.

Methods

Between July 2017 and May 2019, 59 subjects who underwent both conventional and open mouth head and neck CT scans were included in this retrospective study. All images were reconstructed using the O-MAR algorithm. With conventional CT with/without the O-MAR (CTc_O/CTc) and open mouth CT with/without O-MAR (CTo_O/CTo), one reader measured the noise level in multiple anatomic regions of the oral cavity and oropharynx. Visual scores for the streak artifact and overall subjective image quality were assessed by two independent readers.

Results

For the mobile tongue, retromolar trigone, and palatine tonsil, the mean noise was significantly lower, and the mean visual scores were significantly higher, with CTo than with CTc or CTc_O (all, P < 0.001). The mean visual scores were higher with CTo_O than with CTo for the mobile tongue and palatine tonsil (all, P < 0.001). Contrarily, for the mouth floor and tongue base, the mean noise was significantly higher with CTo_O than with CTc or CTc_O, and the mean visual scores were significantly higher with CTc than with CTo or CTo_O (all, P < 0.001).

Conclusions

The open mouth maneuver and O-MAR technique can have different influences on the CT image quality according to the anatomical subsites of the oral cavity and oropharynx.

Metal artifact reduction method based on a constrained beam-hardening estimator for polychromatic x-ray CT

Beam hardening in x-ray computed tomography (CT) is inevitable because of the polychromatic x-ray spectrum and energy-dependent attenuation coefficients of materials, leading to the underestimation of artifacts arising from projection data, especially on metal regions. State-of-the-art research on beam-hardening artifacts is based on a numerical method that recursively performs CT reconstruction, which leads to a heavy computational burden. To address this computational issue, we propose a constrained beam-hardening estimator that provides an efficient numerical solution via a linear combination of two images reconstructed only once during the entire process. The proposed estimator reflects the geometry of metal objects and physical characteristics of beam hardening during the transmission of polychromatic x-rays through a material. Most of the associated parameters are numerically obtained from an initial uncorrected CT image and forward projection transformation without additional optimization procedures. Only the unknown parameter related to beam-hardening artifacts is fine-tuned by linear optimization, which is performed only in the reconstruction image domain. The proposed approach was systematically assessed using numerical simulations and phantom data for qualitative and quantitative comparisons. Compared with existing sinogram inpainting-based and model-based approaches, the proposed scheme in conjunction with the constrained beam-hardening estimator not only provided improved image quality in areas surrounding the metal but also achieved fast beam-hardening correction owing to the analytical reconstruction structure. This work may have significant implications in improving dose calculation accuracy or target volume delineation for treatment planning in radiotherapy.

Deep learning-based metal artefact reduction in PET/CT imaging

Objectives

The susceptibility of CT imaging to metallic objects gives rise to strong streak artefacts and skewed information about the attenuation medium around the metallic implants. This metal-induced artefact in CT images leads to inaccurate attenuation correction in PET/CT imaging. This study investigates the potential of deep learning–based metal artefact reduction (MAR) in quantitative PET/CT imaging.

Methods

Deep learning–based metal artefact reduction approaches were implemented in the image (DLI-MAR) and projection (DLP-MAR) domains. The proposed algorithms were quantitatively compared to the normalized MAR (NMAR) method using simulated and clinical studies. Eighty metal-free CT images were employed for simulation of metal artefact as well as training and evaluation of the aforementioned MAR approaches. Thirty ¹⁸F-FDG PET/CT images affected by the presence of metallic implants were retrospectively employed for clinical assessment of the MAR techniques.

Results

The evaluation of MAR techniques on the simulation dataset demonstrated the superior performance of the DLI-MAR approach (structural similarity (SSIM) = 0.95 ± 0.2 compared to 0.94 ± 0.2 and 0.93 ± 0.3 obtained using DLP-MAR and NMAR, respectively) in minimizing metal artefacts in CT images. The presence of metallic artefacts in CT images or PET attenuation correction maps led to quantitative bias, image artefacts and under- and overestimation of scatter correction of PET images. The DLI-MAR technique led to a quantitative PET bias of 1.3 ± 3% compared to 10.5 ± 6% without MAR and 3.2 ± 0.5% achieved by NMAR.

Conclusion

The DLI-MAR technique was able to reduce the adverse effects of metal artefacts on PET images through the generation of accurate attenuation maps from corrupted CT images.

Key Points

• The presence of metallic objects, such as dental implants, gives rise to severe photon starvation, beam hardening and scattering, thus leading to adverse artefacts in reconstructed CT images.

• The aim of this work is to develop and evaluate a deep learning–based MAR to improve CT-based attenuation and scatter correction in PET/CT imaging.

• Deep learning–based MAR in the image (DLI-MAR) domain outperformed its counterpart implemented in the projection (DLP-MAR) domain. The DLI-MAR approach minimized the adverse impact of metal artefacts on whole-body PET images through generating accurate attenuation maps from corrupted CT images.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-021-07709-z.

Deep Sinogram Completion With Image Prior for Metal Artifact Reduction in CT Images

Computed tomography (CT) has been widely used for medical diagnosis, assessment, and therapy planning and guidance. In reality, CT images may be affected adversely in the presence of metallic objects, which could lead to severe metal artifacts and influence clinical diagnosis or dose calculation in radiation therapy. In this article, we propose a generalizable framework for metal artifact reduction (MAR) by simultaneously leveraging the advantages of image domain and sinogram domain-based MAR techniques. We formulate our framework as a sinogram completion problem and train a neural network (SinoNet) to restore the metal-affected projections. To improve the continuity of the completed projections at the boundary of metal trace and thus alleviate new artifacts in the reconstructed CT images, we train another neural network (PriorNet) to generate a good prior image to guide sinogram learning, and further design a novel residual sinogram learning strategy to effectively utilize the prior image information for better sinogram completion. The two networks are jointly trained in an end-to-end fashion with a differentiable forward projection (FP) operation so that the prior image generation and deep sinogram completion procedures can benefit from each other. Finally, the artifact-reduced CT images are reconstructed using the filtered backward projection (FBP) from the completed sinogram. Extensive experiments on simulated and real artifacts data demonstrate that our method produces superior artifact-reduced results while preserving the anatomical structures and outperforms other MAR methods.

Potentials and caveats of AI in hybrid imaging

State-of-the-art patient management frequently mandates the investigation of both anatomy and physiology of the patients. Hybrid imaging modalities such as the PET/MRI, PET/CT and SPECT/CT have the ability to provide both structural and functional information of the investigated tissues in a single examination. With the introduction of such advanced hardware fusion, new problems arise such as the exceedingly large amount of multi-modality data that requires novel approaches of how to extract a maximum of clinical information from large sets of multi-dimensional imaging data. Artificial intelligence (AI) has emerged as one of the leading technologies that has shown promise in facilitating highly integrative analysis of multi-parametric data. Specifically, the usefulness of AI algorithms in the medical imaging field has been heavily investigated in the realms of (1) image acquisition and reconstruction, (2) post-processing and (3) data mining and modelling. Here, we aim to provide an overview of the challenges encountered in hybrid imaging and discuss how AI algorithms can facilitate potential solutions. In addition, we highlight the pitfalls and challenges in using advanced AI algorithms in the context of hybrid imaging and provide suggestions for building robust AI solutions that enable reproducible and transparent research.

Superiorized method for metal artifact reduction

PURPOSE

Metal artifact reduction (MAR) is a challenging problem in computed tomography (CT) imaging. A popular class of MAR methods replace sinogram measurements that are corrupted by metal with artificial data, typically generated from some combination of interpolation along with other heuristics. While these "projection completion" approaches are successful in eliminating severe artifacts, secondary artifacts may be introduced by the artificial data. In this paper, we propose an approach which uses projection completion to generate a prior image, which is then incorporated into an iterative reconstruction algorithm based on the superiorization framework. The rationale is that the image produced by the iterative algorithm can inherit the desirable properties of the prior image, while also reducing secondary artifacts.

METHODS

The prior image is reconstructed using normalized metal artifact reduction (NMAR), a popular projection completion approach. The iterative algorithm is a modified version of the simultaneous algebraic reconstruction technique (SART), which reduces artifacts by incorporating a polyenergetic forward model, least-squares weighting, and superiorization. The penalty function used for superiorization is a weighted average between a total variation (TV) term and a term promoting similarity with the prior image, similar to penalty functions used in prior image constrained compressive sensing (PICCS). Because the prior is largely free of severe metal artifacts, these artifacts are discouraged from arising during iterative reconstruction; additionally, because the iterative approach uses the original projection data, it is able to recover information that is lost during the NMAR process.

RESULTS

We perform numerical experiments modeling a simple geometric object, as well as several more realistic scenarios such as metal pins, bilateral hip implants, and dental fillings placed within an anatomical phantom. The proposed iterative algorithm is largely successful at eliminating severe metal artifacts as well as secondary artifacts introduced by the NMAR process, especially lost edges of bone structures in the neighborhood of the metal regions. In one case modeling severe photon starvation, the NMAR algorithm is found to provide better results.

CONCLUSION

The proposed algorithm is effective in applying the superiorization methodology to the problem of MAR, providing better results than both NMAR and a purely total variation-based superiorization approach in nearly all cases.

CT metal artifact reduction algorithms: Toward a framework for objective performance assessment

Purpose

Although several metal artifact reduction (MAR) algorithms for computed tomography (CT) scanning are commercially available, no quantitative, rigorous, and reproducible method exists for assessing their performance. The lack of assessment methods poses a challenge to regulators, consumers, and industry. We explored a phantom‐based framework for assessing an important aspect of MAR performance: how applying MAR in the presence of metal affects model observer performance at a low‐contrast detectability (LCD) task This work is, to our knowledge, the first model observer–based framework for the evaluation of MAR algorithms in the published literature.

Methods

We designed a numerical head phantom with metal implants. In order to incorporate an element of randomness, the phantom included a rotatable inset with an inhomogeneous background. We generated simulated projection data for the phantom. We applied two variants of a simple MAR algorithm, sinogram inpainting, to the projection data, that we reconstructed using filtered backprojection. To assess how MAR affected observer performance, we examined the detectability of a signal at the center of a region of interest (ROI) by a channelized Hotelling observer (CHO). As a figure of merit, we used the area under the ROC curve (AUC).

Results

We used simulation to test our framework on two variants of the MAR technique of sinogram inpainting. We found that our method was able to resolve the difference in two different MAR algorithms’ effect on LCD task performance, as well as the difference in task performances when MAR was applied, vs not.

Conclusion

We laid out a phantom‐based framework for objective assessment of how MAR impacts low‐contrast detectability, that we tested on two MAR algorithms. Our results demonstrate the importance of testing MAR performance over a range of object and imaging parameters, since applying MAR does not always improve the quality of an image for a given diagnostic task. Our framework is an initial step toward developing a more comprehensive objective assessment method for MAR, which would require developing additional phantoms and methods specific to various clinical applications of MAR, and increasing study efficiency.

Impact of different iterative metal artifact reduction (iMAR) algorithms on PET/CT attenuation correction after port implementation

PURPOSE

To evaluate the effect of various interactive metal artifact reduction (iMAR) algorithms on attenuation correction in the vicinity of port chambers in PET/CT.

MATERIAL AND METHODS

In this prospective study, 30 oncological patients (12 female, 18 male, mean age 59.6 ± 10.5y) with implanted port chambers undergoing 18F-FDG PET/CT were included. CT images were reconstructed with standard weighted filtered back projection (WFBP) and three different iMAR algorithms (hip, dental filling (DF) and pacemaker (PM)). PET attenuation correction was performed with all four CT datasets. SUVmean, SUVmax and HU measurements were performed in fat and muscle tissue in the vicinity of the port chamber at the location of the strongest bright and dark band artifacts. Differences between HU and SUV values across all CT- and PET-images were investigated using a paired t-test. Bonferroni correction was used to prevent alpha-error accumulation (p < 0.008).

RESULTS

In comparison to WFBP (fat: 94.2 ± 53.9 HU, muscle: 197.6 ± 49.2 HU) all three iMAR algorithms led to a decrease of HU in bright band artifacts. iMAR-DF led to a decrease of 159.2 % (fat: -51.9 ± 58.5 HU, muscle: 94.5 ± 55.3 HU), iMAR-hip of 138.3 % (fat: -30.3 ± 58.5, muscle: 70.4 ± 28.8) and iMAR-PM of 122.3 % (fat: -21.2 ± 47.2 HU, muscle: 72.5 ± 25.1 HU; for all p < 0.008). There was no significant effect of iMAR on SUV measurements in comparison to WFBP.

CONCLUSION

iMAR leads to a significant change of HU values in artifacts caused by port catheter chambers in comparison to WFBP. However, no significant differences in attenuation correction and consecutive changes in SUV measurements can be observed.

Reduction of beam hardening artifacts on real C-arm CT data using polychromatic statistical image reconstruction

PURPOSE

This work aims at the compensation of beam hardening artifacts by the means of an extended three-dimensional polychromatic statistical reconstruction to be applied for flat panel cone-beam CT.

METHODS

We implemented this reconstruction technique as being introduced by Elbakri et al. (2002) [1] for a multi-GPU system, assuming the underlying object consists of several well-defined materials. Furthermore, we assume one voxel can only contain an overlap of at most two materials, depending on its density value. Given the X-ray spectrum, the procedure enables to reconstruct the energy-dependent attenuation values of the volume.

RESULTS

We evaluated the method by using flat-panel cone-beam CT measurements of structures containing small metal objects and clinical head scan data. In comparison with the water-corrected filtered backprojection, as well as a maximum likelihood reconstruction with a consistency-based beam hardening correction, our method features clearly reduced beam hardening artifacts and a more accurate shape of metal objects.

CONCLUSIONS

Our multi-GPU implementation of the polychromatic reconstruction, which does not require any image pre-segmentation, clearly outperforms the standard reconstructions of objects, with respect to beam hardening even in the presence of metal objects inside the volume. However, remaining artifacts, caused mainly by the limited dynamic range of the detector, may have to be addressed in future work.

A dual-stream deep convolutional network for reducing metal streak artifacts in CT images

Machine learning and deep learning are rapidly finding applications in the medical imaging field. In this paper, we address the long-standing problem of metal artifacts in computed tomography (CT) images by training a dual-stream deep convolutional neural network for streak removal. While many metal artifact reduction methods exist, even state-of-the-art algorithms fall short in some clinical applications. Specifically, proton therapy planning requires high image quality with accurate tumor volumes to ensure treatment success. We explore a dual-stream deep network structure with residual learning to correct metal streak artifacts after a first-pass by a state-of-the-art interpolation-based algorithm, NMAR. We provide the network with a mask of the streaks in order to focus attention on those areas. Our experiments compare a mean squared error loss function with a perceptual loss function to emphasize preservation of image features and texture. Both visual and quantitative metrics are used to assess the resulting image quality for metal implant cases. Success may be due to the duality of information processing, with one network stream performing local structure correction, while the other stream provides an attention mechanism to destreak effectively. This study shows that image-domain deep learning can be highly effective for metal artifact reduction (MAR), and highlights the benefits and drawbacks of different loss functions for solving a major CT reconstruction challenge.

MR-based CT metal artifact reduction for head-and-neck photon, electron, and proton radiotherapy

PURPOSE

We investigated the impact on computed tomography (CT) image quality and photon, electron, and proton head-and-neck (H&N) radiotherapy (RT) dose calculations of three CT metal artifact reduction (MAR) approaches: A CT-based algorithm (oMAR Philips Healthcare), manual water override, and our recently presented, Magnetic Resonance (MR)-based kerMAR algorithm. We considered the following three hypotheses: I: Manual water override improves MAR over the CT- and MR-based alternatives; II: The automatic algorithms (oMAR and kerMAR) improve MAR over the uncorrected CT; III: kerMAR improves MAR over oMAR.

METHODS

We included a veal shank phantom with/without six metal inserts and nine H&N RT patients with dental implants. We quantified the MAR capabilities by the reduction of outliers in the CT value distribution in regions of interest, and the change in particle range and photon depth at maximum dose.

RESULTS

Water override provided apparent image improvements in the soft tissue region but insignificantly or negatively influenced the dose calculations. We however found significant improvements in image quality and particle range impact, compared to the uncorrected CT, when using oMAR and kerMAR. kerMAR in turn provided superior improvements in terms of high intensity streak suppression compared to oMAR, again with associated impacts on the particle range estimates.

CONCLUSION

We found no benefits of the water override compared to the rest, and tentatively reject hypothesis I. We however found improvements in the automatic algorithms, and thus support for hypothesis II, and found the MR-based kerMAR to improve upon oMAR, supporting hypothesis III.

Known-component metal artifact reduction (KC-MAR) for cone-beam CT

Intraoperative cone-beam CT (CBCT) is increasingly used for surgical navigation and validation of device placement. In spinal deformity correction, CBCT provides visualization of pedicle screws and fixation rods in relation to adjacent anatomy. This work reports and evaluates a method that uses prior information regarding such surgical instrumentation for improved metal artifact reduction (MAR). The known-component MAR (KC-MAR) approach achieves precise localization of instrumentation in projection images using rigid or deformable 3D-2D registration of component models, thereby overcoming residual errors associated with segmentation-based methods. Projection data containing metal components are processed via 2D inpainting of the detector signal, followed by 3D filtered back-projection (FBP). Phantom studies were performed to identify nominal algorithm parameters and quantitatively investigate performance over a range of component material composition and size. A cadaver study emulating screw and rod placement in spinal deformity correction was conducted to evaluate performance under realistic clinical imaging conditions. KC-MAR demonstrated reduction in artifacts (standard deviation in voxel values) across a range of component types and dose levels, reducing the artifact to 5-10 HU. Accurate component delineation was demonstrated for rigid (screw) and deformable (rod) models with sub-mm registration errors, and a single-pixel dilation of the projected components was found to compensate for partial-volume effects. Artifacts associated with spine screws and rods were reduced by 40%-80% in cadaver studies, and the resulting images demonstrated markedly improved visualization of instrumentation (e.g. screw threads) within cortical margins. The KC-MAR algorithm combines knowledge of surgical instrumentation with 3D image reconstruction in a manner that overcomes potential pitfalls of segmentation. The approach is compatible with FBP-thereby maintaining simplicity in a manner that is consistent with surgical workflow-or more sophisticated model-based reconstruction methods that could further improve image quality and/or help reduce radiation dose.

The effects of the orthopedic metal artifact reduction (O-MAR) algorithm on contouring and dosimetry of head and neck radiotherapy patients

Metallic objects, such as dental fillings, cause artifacts in computed tomography (CT) scans. We quantify the contouring and dosimetric effects of Orthopedic Metal Artifact Reduction (O-MAR), in head and neck radiotherapy. The ease of organ contouring was assessed by having a radiation oncologist identify the CT data set with or without O-MAR for each of 28 patients that was easier to contour. The effect on contouring was quantified further by having the physician recontour parotid glands, previously drawn by him on the O-MAR scans, on uncorrected scans, and calculating the Dice coefficent (a measure of overlap) for the contours. Radiotherapy plans originally generated on scans reconstructed with O-MAR were recalculated on scans without metal artifact correction. The study was done using the Analytical Anisotropic Algorithm (AAA) dose calculation algorithm. The 15 patients with a planning target volume (PTV) extending to the same slice as the artifacts were used for this part of the study. The normal tissue doses were not significantly affected. The PTV mean dose and V95 were not affected, but the cold spots became less severe in the O-MAR corrected plans, with the minimum point dose on average being 4.1% higher. In 79% of the cases, the radiation oncologist identified the O-MAR scan as easier to contour; in 11% he chose the uncorrected scan and in 11% the scans were judged to have equal quality. A total of nine parotid glands (on both scans-18 contours in total) in 5 patients were recontoured. The average Dice coefficient for parotids drawn with and without O-MAR was found to be 0.775 +/- 0.045. The O-MAR algorithm does not produce a significant dosimetric effect in head and neck plans when using the AAA dose calculation algorithm. It can therefore be used for improved contouring accuracy without updating the critical structure tolerance doses and target coverage expectations.

Prior-Guided Metal Artifact Reduction for Iterative X-Ray Computed Tomography

High-attenuation materials pose significant challenges to computed tomographic imaging. Formed of high mass-density and high atomic number elements, they cause more severe beam hardening and scattering artifacts than do water-like materials. Pre-corrected line-integral density measurements are no longer linearly proportional to the path lengths, leading to reconstructed image suffering from streaking artifacts extending from metal, often along highest-density directions. In this paper, a novel prior-based iterative approach is proposed to reduce metal artifacts. It combines the superiority of statistical methods with the benefits of sinogram completion methods to estimate and correct metal-induced biases. Preliminary results show minimized residual artifacts and significantly improved image quality.

Model Image-Based Metal Artifact Reduction for Computed Tomography

Metal implants often produce severe artifacts in the reconstructed computed tomography (CT) images, causing information and image detail loss and making the CT images diagnostically unusable. In order to eliminate the metal artifacts and enhance the diagnostic value of the reconstructed CT images, a post-processing metal artifact reduction algorithm, based on a tissue-class model segmented by thresholding and k-means clustering with spatial information, is proposed. The image inpainting technique is incorporated into the algorithm to improve the segmentation accuracy for CT images severely corrupted by metal artifacts. A study of a water phantom and of two sets of clinical CT images was performed to test the algorithm performance. The proposed method effectively eliminates typical metal artifacts, restores the average CT numbers of different tissues to the proper levels, and preserves the edge and contrast information, thus allowing the accurate reconstruction of the tissue attenuation map. The quality of the artifact-corrected CT images allows them to be subsequently used in other clinical applications, such as three-dimensional rendering, dose estimation for radiotherapy, attenuation correction for PET and SPECT, etc. The algorithm does not rely on the use of the raw sinogram and so is not limited by the proprietary format restrictions.

An additional tilted-scan-based CT metal-artifact-reduction method for radiation therapy planning

Purpose

As computed tomography (CT) imaging is the most commonly used modality for treatment planning in radiation therapy, metal artifacts in the planning CT images may complicate the target delineation and reduce the dose calculation accuracy. Although current CT scanners do provide certain correction steps, it is a common understanding that there is not a universal solution yet to the metal artifact reduction (MAR) in general. Particularly noting the importance of MAR for radiation treatment planning, we propose a novel MAR method in this work that recruits an additional tilted CT scan and synthesizes nearly metal‐artifact‐free CT images.

Methods

The proposed method is based on the facts that the most pronounced metal artifacts in CT images show up along the x‐ray beam direction traversing multiple metallic objects and that a tilted CT scan can provide complementary information free of such metal artifacts in the earlier scan. Although the tilted CT scan would contain its own metal artifacts in the images, the artifacts may manifest in a different fashion leaving a chance to concatenate the two CT images with the metal artifacts much suppressed. We developed an image processing technique that uses the structural similarity (SSIM) for suppressing the metal artifacts. On top of the additional scan, we proposed to use an existing MAR method for each scan if necessary to further suppress the metal artifacts.

Results

The proposed method was validated by a simulation study using the pelvic region of an XCAT numerical phantom and also by an experimental study using the head part of the Rando phantom. The proposed method was found to effectively reduce the metal artifacts. Quantitative analyses revealed that the proposed method reduced the mean absolute percentages of the error by up to 86% and 89% in the simulation and experimental studies, respectively.

Conclusions

It was confirmed that the proposed method, using complementary information acquired from an additional tilted CT scan, can provide nearly metal‐artifact‐free images for the treatment planning.

Reduction of Angularly-Varying-Data Truncation in C-Arm CBCT Imaging

C-arm cone-beam computed tomography (CBCT) has been used increasingly as an imaging tool for yielding 3D anatomical information about the subjects in surgical and interventional procedures. In the clinical applications, the limited field-of-view (FOV) of C-arm CBCT can lead to significant data truncation, resulting in image artifacts that can obscure low contrast tumor embedded within soft-tissue background, thus limiting the utility of C-arm CBCT. The truncation issue can become serious as most of the surgical and interventional procedures would involve devices and tubes that are placed outside the FOV of C-arm CBCT and thus can engender angularly-varying-data truncation. Existing methods may not be adequately applicable to dealing with the angularly-varying truncation. In this work, we seek to reduce truncation artifacts by tailoring optimization-based reconstruction directly from truncated data, without performing pre-reconstruction data compensation, collected from physical phantoms and human subjects. The reconstruction problem is formulated as a constrained optimization program in which a data-derivative-ℓ2-norm fidelity is included for effectively suppressing image artifacts caused by the angularly-varying-data truncation, and the generic Chambolle-Pock algorithm is tailored to solve the optimization program. The results of the study suggest that an appropriately designed optimization-based reconstruction can be exploited for yielding images with reduced artifacts caused by angularly-varying-data truncation.

Convolutional Neural Network Based Metal Artifact Reduction in X-Ray Computed Tomography

In the presence of metal implants, metal artifacts are introduced to x-ray computed tomography CT images. Although a large number of metal artifact reduction (MAR) methods have been proposed in the past decades, MAR is still one of the major problems in clinical x-ray CT. In this paper, we develop a convolutional neural network (CNN)-based open MAR framework, which fuses the information from the original and corrected images to suppress artifacts. The proposed approach consists of two phases. In the CNN training phase, we build a database consisting of metal-free, metal-inserted and pre-corrected CT images, and image patches are extracted and used for CNN training. In the MAR phase, the uncorrected and pre-corrected images are used as the input of the trained CNN to generate a CNN image with reduced artifacts. To further reduce the remaining artifacts, water equivalent tissues in a CNN image are set to a uniform value to yield a CNN prior, whose forward projections are used to replace the metal-affected projections, followed by the FBP reconstruction. The effectiveness of the proposed method is validated on both simulated and real data. Experimental results demonstrate the superior MAR capability of the proposed method to its competitors in terms of artifact suppression and preservation of anatomical structures in the vicinity of metal implants.

Metal Artifact Reduction in X-ray Computed Tomography Using Computer-Aided Design Data of Implants as Prior Information

OBJECTIVES

The performance of metal artifact reduction (MAR) methods in x-ray computed tomography (CT) suffers from incorrect identification of metallic implants in the artifact-affected volumetric images. The aim of this study was to investigate potential improvements of state-of-the-art MAR methods by using prior information on geometry and material of the implant.

MATERIALS AND METHODS

The influence of a novel prior knowledge-based segmentation (PS) compared with threshold-based segmentation (TS) on 2 MAR methods (linear interpolation [LI] and normalized-MAR [NORMAR]) was investigated. The segmentation is the initial step of both MAR methods. Prior knowledge-based segmentation uses 3-dimensional registered computer-aided design (CAD) data as prior knowledge to estimate the correct position and orientation of the metallic objects. Threshold-based segmentation uses an adaptive threshold to identify metal. Subsequently, for LI and NORMAR, the selected voxels are projected into the raw data domain to mark metal areas. Attenuation values in these areas are replaced by different interpolation schemes followed by a second reconstruction. Finally, the previously selected metal voxels are replaced by the metal voxels determined by PS or TS in the initial reconstruction. First, we investigated in an elaborate phantom study if the knowledge of the exact implant shape extracted from the CAD data provided by the manufacturer of the implant can improve the MAR result. Second, the leg of a human cadaver was scanned using a clinical CT system before and after the implantation of an artificial knee joint. The results were compared regarding segmentation accuracy, CT number accuracy, and the restoration of distorted structures.

RESULTS

The use of PS improved the efficacy of LI and NORMAR compared with TS. Artifacts caused by insufficient segmentation were reduced, and additional information was made available within the projection data. The estimation of the implant shape was more exact and not dependent on a threshold value. Consequently, the visibility of structures was improved when comparing the new approach to the standard method. This was further confirmed by improved CT value accuracy and reduced image noise.

CONCLUSIONS

The PS approach based on prior implant information provides image quality which is superior to TS-based MAR, especially when the shape of the metallic implant is complex. The new approach can be useful for improving MAR methods and dose calculations within radiation therapy based on the MAR corrected CT images.

Dosimetry Effects Caused by Unilateral and Bilateral Hip Prostheses: A Monte Carlo Case Study in Megavoltage Photon Radiotherapy for Computed Tomography Data without Metal Artifacts

Background:

Hip prostheses (HPs) are routinely used in hip augmentation to replace painful or dysfunctional hip joints. However, high-density and high-atomic-number (Z) inserts may cause dose perturbations in the target volume and interface regions.

Aim:

To evaluate the dosimetric influence of various HPs during megavoltage conformal radiotherapy (RT) of the prostate using Monte Carlo (MC) simulations.

Materials and Methods:

BEAMnrc and DOSXYZnrc MC user-codes were respectively used to simulate the linac head and to calculate 3D absorbed dose distributions in a computed tomography (CT)-based phantom. A novel technique was used to synthetically introduce HPs into the raw patient CT dataset. The prosthesis materials evaluated were stainless steel (SS316L), titanium (Ti6Al4V), and ultra-high-molecular-weight polyethylene (UHMWPE). Four, five, and six conformal photon fields of 6–20 MV were used.

Results:

The absorbed dose within and beyond metallic prostheses dropped significantly due to beam attenuation. For bilateral HPs, the target dose reduction ranged up to 23% and 17% for SS316L and Ti6Al4V, respectively. For unilateral HP, the respective dose reductions were 19% and 12%. Dose enhancement was always <1% for UHMWPE. The 6-field plan produced the best target coverage. Up to 38% dose increase was found at the bone–SS316L proximal interface.

Conclusions:

The novel technique used enabled the complete exclusion of metal artifacts in the CT dataset. High-energy plans with more oblique beams can help minimize dose attenuation through HPs. Shadowing and interface effects are density dependent and greatest for SS316L, while UHMWPE poses negligible dose perturbation.

A new approach for reducing beam hardening artifacts in polychromatic X-ray computed tomography using more accurate prior image

PURPOSE

Metal artifacts severely degrade CT image quality in clinical diagnosis, which are difficult to removed, especially for the beam hardening artifacts. The metal artifact reduction (MAR) based on prior images are the most frequently-used methods. However, there exists a lot misclassification in most prior images caused by absence of prior information such as the spectrum distribution of X-ray beam source, especially many or big metal included. The purpose of this work is to find a more accurate prior image to improve image quality.

METHODS

The proposed method comprise of following four steps. First, the metal image is segmented by thresholding an initial image, where the metal traces are identified in the initial projection data using the forward projection of the metal image. Second, the accurate absorbent model of certain metal image is calculated according to the spectrum distribution of certain X-ray beam source and energy-dependent attenuation coefficients of metal. Then, a new metal image is reconstructed by the general analytical reconstruction algorithm such as filtered back projection (FPB). The prior image is obtained by segmenting the difference image between the initial image and the new metal image into air, tissue and bone. Finally, the initial projection data are normalized by dividing the projection data of prior image pixel to pixel, the corrected image is obtained by interpolation, denormalization and reconstruction.

RESULTS

Some clinical images with dental fillings and knee prostheses are used to evaluate the proposed algorithm and normalized metal artifact reduction (NMAR) and linear interpolation (LI) method. The results demonstrate the artifacts can be reduced efficiently by the proposed method.

CONCLUSIONS

The proposed method could obtain an exact prior image using the prior information about X-ray beam source and energy-dependent attenuation coefficients of metal. As a result, the better performance of reducing beam hardening artifacts can be improved, even though there were many or big implants. Moreover, the process of the proposed method is rather simple and little extra calculation burden is necessary. It has superiorities over other algorithms when include big or many implants.

Reduction of Metallic Artifacts of the Post-treatment Intracranial Aneurysms: Effects of Single Energy Metal Artifact Reduction Algorithm

PURPOSE

This study evaluated the quality of computed tomography (CT) and CT angiography images generated using the single-energy metal artifact reduction (SEMAR) algorithm during perfusion examination in patients who had undergone reconstruction with neurosurgical clipping or endovascular coiling for treatment of aneurysms.

METHODS

A total of 55 patients with implanted intracranial clips or coils (24 men and 31 women; mean age 60.15 ± 15.86 years) underwent perfusion studies evaluated by CT and CT angiography with a 320-row CT scanner. Images were reconstructed with either the SEMAR algorithm combined with iterative reconstruction (SEMAR group), or by iterative reconstruction only (non-SEMAR group control). The SEMAR and control images were compared for artifacts (index and maximum diameter), image quality, cerebral perfusion parameters, noise (images with the worst artifacts), and contrast-to-noise ratio. The metallic artifacts were visually evaluated by two radiologists using a four-point scale in a double-blinded manner.

RESULTS

The noise, artifact diameter, and artifact index of the SEMAR images were significantly lower than that of the control images, and the subjective image quality score and contrast-to-noise ratio were significantly higher (P < 0.01, all). The cerebral perfusion parameters of the SEMAR and control images were comparable (i. e. blood flow, blood volume, and mean transit time).

CONCLUSION

For imaging intracranial metallic implants, the SEMAR algorithm produced images with significantly fewer artifacts than the iterative reconstruction alone, with no statistical changes in perfusion parameters. Thus, SEMAR reconstruction can be instrumental in improving CT image quality and may ultimately improve the detection of postoperative complications and patient prognosis.

A metal artifact reduction algorithm in CT using multiple prior images by recursive active contour segmentation

We propose a novel metal artifact reduction (MAR) algorithm for CT images that completes a corrupted sinogram along the metal trace region. When metal implants are located inside a field of view, they create a barrier to the transmitted X-ray beam due to the high attenuation of metals, which significantly degrades the image quality. To fill in the metal trace region efficiently, the proposed algorithm uses multiple prior images with residual error compensation in sinogram space. Multiple prior images are generated by applying a recursive active contour (RAC) segmentation algorithm to the pre-corrected image acquired by MAR with linear interpolation, where the number of prior image is controlled by RAC depending on the object complexity. A sinogram basis is then acquired by forward projection of the prior images. The metal trace region of the original sinogram is replaced by the linearly combined sinogram of the prior images. Then, the additional correction in the metal trace region is performed to compensate the residual errors occurred by non-ideal data acquisition condition. The performance of the proposed MAR algorithm is compared with MAR with linear interpolation and the normalized MAR algorithm using simulated and experimental data. The results show that the proposed algorithm outperforms other MAR algorithms, especially when the object is complex with multiple bone objects.

Dual-Energy CT for the Musculoskeletal System

The principal advantages of dual-energy computed tomography (CT) over conventional CT in the musculoskeletal setting relate to the additional information provided regarding tissue composition, artifact reduction, and image optimization. This article discusses the manifestations of these in clinical practice-urate and bone marrow edema detection, metal artifact reduction, and tendon analysis, with potential in arthrography, bone densitometry, and metastases surveillance. The basic principles of dual-energy CT physics and scanner design will also be discussed. ^© RSNA, 2016.

Metal Artifact Reduction of CT Scans to Improve PET/CT

In recent years, different metal artifact reduction methods have been developed for CT. These methods have only recently been introduced for PET/CT even though they could be beneficial for interpretation, segmentation, and quantification of the PET/CT images. In this study, phantom and patient scans were analyzed visually and quantitatively to measure the effect on PET images of iterative metal artifact reduction (iMAR) of CT data. Methods: The phantom consisted of 2 types of hip prostheses in a solution of ¹⁸F-FDG and water. ¹⁸F-FDG PET/CT scans of 14 patients with metal implants (either dental implants, hip prostheses, shoulder prostheses, or pedicle screws) and ⁶⁸Ga-labeled prostate-specific membrane antigen (⁶⁸Ga-PSMA) PET/CT scans of 7 patients with hip prostheses were scored by 2 experienced nuclear medicine physicians to analyze clinical relevance. For all patients, a lesion was located in the field of view of the metal implant. Phantom and patients were scanned in a PET/CT scanner. The standard low-dose CT scans were processed with the iMAR algorithm. The PET data were reconstructed using attenuation correction provided by both standard CT and iMAR-processed CT. Results: For the phantom scans, cold artifacts were visible on the PET image. There was a 30% deficit in ¹⁸F-FDG concentration, which was restored by iMAR processing, indicating that metal artifacts on CT images induce quantification errors in PET data. The iMAR algorithm was useful for most patients. When iMAR was used, the confidence in interpretation increased or stayed the same, with an average improvement of 28% ± 20% (scored on a scale of 0%-100% confidence). The SUV increase or decrease depended on the type of metal artifact. The mean difference in absolute values of SUV_mean of the lesions was 3.5% ± 3.3%. Conclusion: The iMAR algorithm increases the confidence of the interpretation of the PET/CT scan and influences the SUV. The added value of iMAR depends on the indication for the PET/CT scan, location and size/type of the prosthesis, and location and extent of the disease.

Polyenergetic known-component CT reconstruction with unknown material compositions and unknown x-ray spectra

Metal artifacts can cause substantial image quality issues in computed tomography. This is particularly true in interventional imaging where surgical tools or metal implants are in the field-of-view. Moreover, the region-of-interest is often near such devices which is exactly where image quality degradations are largest. Previous work on known-component reconstruction (KCR) has shown the incorporation of a physical model (e.g. shape, material composition, etc) of the metal component into the reconstruction algorithm can significantly reduce artifacts even near the edge of a metal component. However, for such approaches to be effective, they must have an accurate model of the component that include energy-dependent properties of both the metal device and the CT scanner, placing a burden on system characterization and component material knowledge. In this work, we propose a modified KCR approach that adopts a mixed forward model with a polyenergetic model for the component and a monoenergetic model for the background anatomy. This new approach called Poly-KCR jointly estimates a spectral transfer function associated with known components in addition to the background attenuation values. Thus, this approach eliminates both the need to know component material composition a prior as well as the requirement for an energy-dependent characterization of the CT scanner. We demonstrate the efficacy of this novel approach and illustrate its improved performance over traditional and model-based iterative reconstruction methods in both simulation studies and in physical data including an implanted cadaver sample.

Image-based Metal Artifact Reduction in X-ray Computed Tomography utilizing Local Anatomical Similarity

X-ray computed tomography (CT) is widely used in radiation therapy treatment planning in recent years. However, metal implants such as dental fillings and hip prostheses can cause severe bright and dark streaking artifacts in reconstructed CT images. These artifacts decrease image contrast and degrade HU accuracy, leading to inaccuracies in target delineation and dose calculation. In this work, a metal artifact reduction method is proposed based on the intrinsic anatomical similarity between neighboring CT slices. Neighboring CT slices from the same patient exhibit similar anatomical features. Exploiting this anatomical similarity, a gamma map is calculated as a weighted summation of relative HU error and distance error for each pixel in an artifact-corrupted CT image relative to a neighboring, artifact-free image. The minimum value in the gamma map for each pixel is used to identify an appropriate pixel from the artifact-free CT slice to replace the corresponding artifact-corrupted pixel. With the proposed method, the mean CT HU error was reduced from 360 HU and 460 HU to 24 HU and 34 HU on head and pelvis CT images, respectively. Dose calculation accuracy also improved, as the dose difference was reduced from greater than 20% to less than 4%. Using 3%/3mm criteria, the gamma analysis failure rate was reduced from 23.25% to 0.02%. An image-based metal artifact reduction method is proposed that replaces corrupted image pixels with pixels from neighboring CT slices free of metal artifacts. This method is shown to be capable of suppressing streaking artifacts, thereby improving HU and dose calculation accuracy.

Iterative metal artifact reduction for x-ray computed tomography using unmatched projector/backprojector pairs

PURPOSE

Metal artifact reduction (MAR) is a major problem and a challenging issue in x-ray computed tomography (CT) examinations. Iterative reconstruction from sinograms unaffected by metals shows promising potential in detail recovery. This reconstruction has been the subject of much research in recent years. However, conventional iterative reconstruction methods easily introduce new artifacts around metal implants because of incomplete data reconstruction and inconsistencies in practical data acquisition. Hence, this work aims at developing a method to suppress newly introduced artifacts and improve the image quality around metal implants for the iterative MAR scheme.

METHODS

The proposed method consists of two steps based on the general iterative MAR framework. An uncorrected image is initially reconstructed, and the corresponding metal trace is obtained. The iterative reconstruction method is then used to reconstruct images from the unaffected sinogram. In the reconstruction step of this work, an iterative strategy utilizing unmatched projector/backprojector pairs is used. A ramp filter is introduced into the back-projection procedure to restrain the inconsistency components in low frequencies and generate more reliable images of the regions around metals. Furthermore, a constrained total variation (TV) minimization model is also incorporated to enhance efficiency. The proposed strategy is implemented based on an iterative FBP and an alternating direction minimization (ADM) scheme, respectively. The developed algorithms are referred to as "iFBP-TV" and "TV-FADM," respectively. Two projection-completion-based MAR methods and three iterative MAR methods are performed simultaneously for comparison.

RESULTS

The proposed method performs reasonably on both simulation and real CT-scanned datasets. This approach could reduce streak metal artifacts effectively and avoid the mentioned effects in the vicinity of the metals. The improvements are evaluated by inspecting regions of interest and by comparing the root-mean-square errors, normalized mean absolute distance, and universal quality index metrics of the images. Both iFBP-TV and TV-FADM methods outperform other counterparts in all cases. Unlike the conventional iterative methods, the proposed strategy utilizing unmatched projector/backprojector pairs shows excellent performance in detail preservation and prevention of the introduction of new artifacts.

CONCLUSIONS

Qualitative and quantitative evaluations of experimental results indicate that the developed method outperforms classical MAR algorithms in suppressing streak artifacts and preserving the edge structural information of the object. In particular, structures lying close to metals can be gradually recovered because of the reduction of artifacts caused by inconsistency effects.

PET/MRI in the Presence of Metal Implants: Completion of the Attenuation Map from PET Emission Data

We present a novel technique for accurate whole-body attenuation correction in the presence of metallic endoprosthesis, on integrated non-time-of-flight (non-TOF) PET/MRI scanners. The proposed implant PET-based attenuation map completion (IPAC) method performs a joint reconstruction of radioactivity and attenuation from the emission data to determine the position, shape, and linear attenuation coefficient (LAC) of metallic implants. Methods: The initial estimate of the attenuation map was obtained using the MR Dixon method currently available on the Siemens Biograph mMR scanner. The attenuation coefficients in the area of the MR image subjected to metal susceptibility artifacts are then reconstructed from the PET emission data using the IPAC algorithm. The method was tested on 11 subjects presenting 13 different metallic implants, who underwent CT and PET/MR scans. Relative mean LACs and Dice similarity coefficients were calculated to determine the accuracy of the reconstructed attenuation values and the shape of the metal implant, respectively. The reconstructed PET images were compared with those obtained using the reference CT-based approach and the Dixon-based method. Absolute relative change (aRC) images were generated in each case, and voxel-based analyses were performed. Results: The error in implant LAC estimation, using the proposed IPAC algorithm, was 15.7% ± 7.8%, which was significantly smaller than the Dixon- (100%) and CT- (39%) derived values. A mean Dice similarity coefficient of 73% ± 9% was obtained when comparing the IPAC- with the CT-derived implant shape. The voxel-based analysis of the reconstructed PET images revealed quantification errors (aRC) of 13.2% ± 22.1% for the IPAC- with respect to CT-corrected images. The Dixon-based method performed substantially worse, with a mean aRC of 23.1% ± 38.4%. Conclusion: We have presented a non-TOF emission-based approach for estimating the attenuation map in the presence of metallic implants, to be used for whole-body attenuation correction in integrated PET/MR scanners. The Graphics Processing Unit implementation of the algorithm will be included in the open-source reconstruction toolbox Occiput.io.

Modeling and Pre-Treatment of Photon-Starved CT Data for Iterative Reconstruction

An increasing number of X-ray CT procedures are being conducted with drastically reduced dosage, due at least in part to advances in statistical reconstruction methods that can deal more effectively with noise than can traditional techniques. As data become photon-limited, more detailed models are necessary to deal with count rates that drop to the levels of system electronic noise. We present two options for sinogram pre-treatment that can improve the performance of photon-starved measurements, with the intent of following with model-based image reconstruction. Both the local linear minimum mean-squared error (LLMMSE) filter and pointwise Bayesian restoration (PBR) show promise in extracting useful, quantitative information from very low-count data by reducing local bias while maintaining the lower noise variance of statistical methods. Results from clinical data demonstrate the potential of both techniques.

Reduce beam hardening artifacts of polychromatic X-ray computed tomography by an iterative approximation approach

BACKGROUND

The beam hardening artifact is one of most important modalities of metal artifact for polychromatic X-ray computed tomography (CT), which can impair the image quality seriously.

OBJECTIVE

An iterative approach is proposed to reduce beam hardening artifact caused by metallic components in polychromatic X-ray CT.

METHODS

According to Lambert-Beer law, the (detected) projections can be expressed as monotonic nonlinear functions of element geometry projections, which are the theoretical projections produced only by the pixel intensities (image grayscale) of certain element (component). With help of a prior knowledge on spectrum distribution of X-ray beam source and energy-dependent attenuation coefficients, the functions have explicit expressions. Newton-Raphson algorithm is employed to solve the functions. The solutions are named as the synthetical geometry projections, which are the nearly linear weighted sum of element geometry projections with respect to mean of each attenuation coefficient. In this process, the attenuation coefficients are modified to make Newton-Raphson iterative functions satisfy the convergence conditions of fixed pointed iteration(FPI) so that the solutions will approach the true synthetical geometry projections stably. The underlying images are obtained using the projections by general reconstruction algorithms such as the filtered back projection (FBP). The image gray values are adjusted according to the attenuation coefficient means to obtain proper CT numbers.

RESULTS

Several examples demonstrate the proposed approach is efficient in reducing beam hardening artifacts and has satisfactory performance in the term of some general criteria. In a simulation example, the normalized root mean square difference (NRMSD) can be reduced 17.52% compared to a newest algorithm.

CONCLUSIONS

Since the element geometry projections are free from the effect of beam hardening, the nearly linear weighted sum of them, the synthetical geometry projections, are almost free from the effect of beam hardening. By working out the synthetical geometry projections, the proposed approach becomes quite efficient in reducing beam hardening artifacts.

Optimization-based image reconstruction with artifact reduction in C-arm CBCT

We investigate an optimization-based reconstruction, with an emphasis on image-artifact reduction, from data collected in C-arm cone-beam computed tomography (CBCT) employed in image-guided interventional procedures. In the study, an image to be reconstructed is formulated as a solution to a convex optimization program in which a weighted data divergence is minimized subject to a constraint on the image total variation (TV); a data-derivative fidelity is introduced in the program specifically for effectively suppressing dominant, low-frequency data artifact caused by, e.g. data truncation; and the Chambolle-Pock (CP) algorithm is tailored to reconstruct an image through solving the program. Like any other reconstructions, the optimization-based reconstruction considered depends upon numerous parameters. We elucidate the parameters, illustrate their determination, and demonstrate their impact on the reconstruction. The optimization-based reconstruction, when applied to data collected from swine and patient subjects, yields images with visibly reduced artifacts in contrast to the reference reconstruction, and it also appears to exhibit a high degree of robustness against distinctively different anatomies of imaged subjects and scanning conditions of clinical significance. Knowledge and insights gained in the study may be exploited for aiding in the design of practical reconstructions of truly clinical-application utility.

High-kVp Assisted Metal Artifact Reduction for X-ray Computed Tomography

In X-ray computed tomography (CT), the presence of metallic parts in patients causes serious artifacts and degrades image quality. Many algorithms were published for metal artifact reduction (MAR) over the past decades with various degrees of success but without a perfect solution. Some MAR algorithms are based on the assumption that metal artifacts are due only to strong beam hardening and may fail in the case of serious photon starvation. Iterative methods handle photon starvation by discarding or underweighting corrupted data, but the results are not always stable and they come with high computational cost. In this paper, we propose a high-kVp-assisted CT scan mode combining a standard CT scan with a few projection views at a high-kVp value to obtain critical projection information near the metal parts. This method only requires minor hardware modifications on a modern CT scanner. Two MAR algorithms are proposed: dual-energy normalized MAR (DNMAR) and high-energy embedded MAR (HEMAR), aiming at situations without and with photon starvation respectively. Simulation results obtained with the CT simulator CatSim demonstrate that the proposed DNMAR and HEMAR methods can eliminate metal artifacts effectively.

Performances of low-dose dual-energy CT in reducing artifacts from implanted metallic orthopedic devices

OBJECTIVES

The objective was to evaluate the performances of dose-reduced dual-energy computed tomography (DECT) in decreasing metallic artifacts from orthopedic devices compared with dose-neutral DECT, dose-neutral single-energy computed tomography (SECT), and dose-reduced SECT.

MATERIALS AND METHODS

Thirty implants in 20 consecutive cadavers underwent both SECT and DECT at three fixed CT dose indexes (CTDI): 20.0, 10.0, and 5.0 mGy. Extrapolated monoenergetic DECT images at 64, 69, 88, 105, 120, and 130 keV, and individually adjusted monoenergy for optimized image quality (OPTkeV) were generated. In each group, the image quality of the seven monoenergetic images and of the SECT image was assessed qualitatively and quantitatively by visually rating and by measuring the maximum streak artifact respectively.

RESULTS

The comparison between SECT and OPTkeV evaluated overall within all groups showed a significant difference (p <0.001), with OPTkeV images providing better images. Comparing OPTkeV with the other DECT images, a significant difference was shown (p <0.001), with OPTkeV and 130-keV images providing the qualitatively best results. The OPTkeV images of 5.0-mGy acquisitions provided percentages of images with scores 1 and 2 of 36 % and 30 % respectively, compared with 0 % and 33.3 % of the corresponding SECT images of 10- and 20-mGy acquisitions. Moreover, DECT reconstructions at the OPTkeV of the low-dose group showed higher CT numbers than the SECT images of dose groups 1 and 2.

CONCLUSIONS

This study demonstrates that low-dose DECT permits a reduction of artifacts due to metallic implants to be obtained in a similar manner to neutral-dose DECT and better than reduced or neutral-dose SECT.

Quantifying metal-induced susceptibility artifacts of the instrumented spine at 1.5T using fast-spin echo and 3D-multispectral MRI

PURPOSE

To evaluate magnetic resonance imaging (MRI) artifacts near metallic spinal instrumentation using both conventional metal artifact reduction sequences (MARS) and 3D multispectral imaging sequences (3D-MSI).

MATERIALS AND METHODS

Both MARS and 3D-MSI images were acquired in 10 subjects with titanium spinal hardware on a 1.5T GE 450W scanner. Clinical computed tomography (CT) images were used to measure the volume of the implant using seed-based region growing. Using 30-40 landmarks, the MARS and 3D-MSI images were coregistered to the CT images. Three independent users manually segmented the artifact volume from both MR sequences. For five L-spine subjects, one user independently segmented the nerve root in both MARS and 3D-MSI images.

RESULTS

For all 10 subjects, the measured artifact volume for the 3D-MSI images closely matched that of the CT implant volume (absolute error: 4.3 ± 2.0 cm³ ). The MARS artifact volume was ∼8-fold higher than that of the 3D-MSI images (30.7 ± 20.2, P = 0.002). The average nerve root volume for the MARS images was 24 ± 7.3% lower than the 3D-MSI images (P = 0.06).

CONCLUSION

Compared to 3D-MSI images, the higher-resolution MARS images may help study features farther away from the implant surface. However, the MARS images retained substantial artifacts in the slice-dimension that result in a larger artifact volume. These artifacts have the potential to obscure physiologically relevant features, and can be mitigated with 3D-MSI sequences. Hence, MR study protocols may benefit with the inclusion both MARS and 3D-MSI sequences to accurately study pathology near the spine.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:51-58.