Dataset Information

Sequentially reweighted TV minimization for CT metal artifact reduction.

ABSTRACT: Metal artifact reduction has long been an important topic in x-ray CT image reconstruction. In this work, the authors propose an iterative method that sequentially minimizes a reweighted total variation (TV) of the image and produces substantially artifact-reduced reconstructions.A sequentially reweighted TV minimization algorithm is proposed to fully exploit the sparseness of image gradients (IG). The authors first formulate a constrained optimization model that minimizes a weighted TV of the image, subject to the constraint that the estimated projection data are within a specified tolerance of the available projection measurements, with image non-negativity enforced. The authors then solve a sequence of weighted TV minimization problems where weights used for the next iteration are computed from the current solution. Using the complete projection data, the algorithm first reconstructs an image from which a binary metal image can be extracted. Forward projection of the binary image identifies metal traces in the projection space. The metal-free background image is then reconstructed from the metal-trace-excluded projection data by employing a different set of weights. Each minimization problem is solved using a gradient method that alternates projection-onto-convex-sets and steepest descent. A series of simulation and experimental studies are performed to evaluate the proposed approach.Our study shows that the sequentially reweighted scheme, by altering a single parameter in the weighting function, flexibly controls the sparsity of the IG and reconstructs artifacts-free images in a two-stage process. It successfully produces images with significantly reduced streak artifacts, suppressed noise and well-preserved contrast and edge properties.The sequentially reweighed TV minimization provides a systematic approach for suppressing CT metal artifacts. The technique can also be generalized to other "missing data" problems in CT image reconstruction.

SUBMITTER: Zhang X

PROVIDER: S-EPMC3710254 | biostudies-other | 2013 Jul

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:PurposeAs 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.MethodsThe 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.ResultsThe 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.ConclusionsIt 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.

Project description:BackgroundMetal artifacts due to spinal implants may affect the accuracy of dose calculation for radiotherapy. However, the dosimetric impact of metal artifact reduction (MAR) for spinal implants in stereotactic body radiotherapy (SBRT) plans has not been well studied. The objective of this study was to evaluate the dosimetric impact of MAR in spinal SBRT planning with three clinically common dose calculation algorithms.MethodsGammex phantom and 10 patients' computed tomography (CT) images were studied to investigate the effects of titanium implants. A commercial orthopedic MAR algorithm was employed to reduce artifacts. Dose calculations for SBRT were conducted on both artifact-corrected and uncorrected images using three commercial algorithms [analytical anisotropic algorithm (AAA), Acuros XB (AXB), and Monte Carlo (MC)]. Dose discrepancies between artifact-corrected and uncorrected cases were appraised using a dose-volume histogram (DVH) and 3-dimensional (3D) gamma analysis with different distance to agreement (DTA) and dose difference criteria. The gamma agreement index (GAI) was denoted as G(∆D, DTA). Statistical analysis of t-test was utilized to evaluate the dose differences of different algorithms.ResultsThe phantom study demonstrated that titanium metal artifacts can be effectively reduced. The patient cases study showed that dose differences between the artifact-corrected and uncorrected datasets were small evaluated by gamma index and DVH. Gamma analysis found that even the strict criterion local G(1,1) had average values ≥93.9% for the three algorithms. For all DVH metrics, average differences did not exceed 0.7% in planning target volume (PTV) and 2.1% in planning risk volume of spinal cord (PRV-SC). Statistical analysis showed that the observed dose differences of MC method were significantly larger than those of AAA (P<0.01 for D98% of PTV and P<0.001 for D0.1cc of spinal cord) and AXB methods (P<0.001 for D98% and P<0.0001 for D0.1cc).ConclusionsDosimetric impact of artifacts caused by titanium implants is not significant in spinal SBRT planning, which indicates that dose calculation algorithms might not be very sensitive to CT number variation caused by titanium inserts.

Project description:Severe artifacts in kilovoltage-CT simulation images caused by large metallic implants can significantly degrade the conspicuity and apparent CT Hounsfield number of targets and anatomic structures, jeopardize the confidence of anatomical segmentation, and introduce inaccuracies into the radiation therapy treatment planning process. This study evaluated the performance of the first commercial orthopedic metal artifact reduction function (O-MAR) for radiation therapy, and investigated its clinical applications in treatment planning.Both phantom and clinical data were used for the evaluation. The CIRS electron density phantom with known physical (and electron) density plugs and removable titanium implants was scanned on a Philips Brilliance Big Bore 16-slice CT simulator. The CT Hounsfield numbers of density plugs on both uncorrected and O-MAR corrected images were compared. Treatment planning accuracy was evaluated by comparing simulated dose distributions computed using the true density images, uncorrected images, and O-MAR corrected images. Ten CT image sets of patients with large hip implants were processed with the O-MAR function and evaluated by two radiation oncologists using a five-point score for overall image quality, anatomical conspicuity, and CT Hounsfield number accuracy. By utilizing the same structure contours delineated from the O-MAR corrected images, clinical IMRT treatment plans for five patients were computed on the uncorrected and O-MAR corrected images, respectively, and compared.Results of the phantom study indicated that CT Hounsfield number accuracy and noise were improved on the O-MAR corrected images, especially for images with bilateral metal implants. The γ pass rates of the simulated dose distributions computed on the uncorrected and O-MAR corrected images referenced to those of the true densities were higher than 99.9% (even when using 1% and 3 mm distance-to-agreement criterion), suggesting that dose distributions were clinically identical. In all patient cases, radiation oncologists rated O-MAR corrected images as higher quality. Formerly obscured critical structures were able to be visualized. The overall image quality and the conspicuity in critical organs were significantly improved compared with the uncorrected images: overall quality score (1.35 vs 3.25, P = 0.0022); bladder (2.15 vs 3.7, P = 0.0023); prostate and seminal vesicles∕vagina (1.3 vs 3.275, P = 0.0020); rectum (2.8 vs 3.9, P = 0.0021). The noise levels of the selected ROIs were reduced from 93.7 to 38.2 HU. On most cases (8∕10), the average CT Hounsfield numbers of the prostate∕vagina on the O-MAR corrected images were closer to the referenced value (41.2 HU, an average measured from patients without metal implants) than those on the uncorrected images. High γ pass rates of the five IMRT dose distribution pairs indicated that the dose distributions were not significantly affected by the CT image improvements.Overall, this study indicated that the O-MAR function can remarkably reduce metal artifacts and improve both CT Hounsfield number accuracy and target and critical structure visualization. Although there was no significant impact of the O-MAR algorithm on the calculated dose distributions, we suggest that O-MAR corrected images are more suitable for the entire treatment planning process by offering better anatomical structure visualization, improving radiation oncologists' confidence in target delineation, and by avoiding subjective density overrides of artifact regions on uncorrected images.

Dataset Information

Sequentially reweighted TV minimization for CT metal artifact reduction.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets