Comprehensive Survey on Metal Artifact Reduction Methods in Computed Tomography Images

1. Introduction

Computed tomography (CT) has become an invaluable imaging tool for the diagnosis or further evaluation of various medical conditions. With evolving technology, decreasing acquisition times, and improving resolution, the usage of CT has increased tremendously over recent years.^1.2 It is estimated that more than 85 million CT scans were performed in the United States in 2011.³ The increasing use of CT has been documented by a number of investigators.^{4, 5} Trauma is one of the medical domains to attract wide usage of CT modality. In India more than 1.2 million patients are reported with trauma cases every year. Among these trauma cases one common problem reported is streaking artifacts in CT images near the metal implant.

Metallic implants inside the legs, head, femur or teeth, aneurysm clips/coils(brain), dental fillings, pacer wires, surgical clips, bullet/shrapnel/lead shot, shoulder/hip/knee replacement, temporal neck screw, spinal rods, orthopedic plates, pedicle screws, depth electrodes (brain), cryoblation probes, iodinated contrast which cannot be removed from the body during CT study are reported to produces metal artifact.⁶ Metal artifact comes from the fact that X-ray in the diagnostic energy range is strongly attenuated after passing through metal pieces, and much fewer signals are detected by the detectors in these regions. Beam hardening, scatter, noise, exponential edge gradient effect and under-sampling are reported to be the main causes of this metal artifact.⁷ These artifacts deteriorate image quality and often superimpose structures of interest leading to prevention of recognition and management of complications of areas close to the implants.

Filtered Back Projection (FBP) algorithm is the most widely used reconstruction technique in X-ray CT, where transformations between the sinogram (projection) and image domains are performed using the forward and inverse radon transforms respectively.⁸ The radon transform of a function , denoted as , is defined as its line integral along a line inclined at an angle from the y-axis and at a distance “s” from the origin. It is also called as sinogram and is written mathematically as:

Given the sonogram , we would like to reconstruct the image described in (x, y) coordinates. Basically it works by subsequently “smearing” the acquired , across an empty plane. This is referred as simple back projection, and is given as:

The simple back projection yields a blurry laminogram (reconstructed image from sinogram) that must be deconvolved by a filter to yield the original image. This shall be done by filtering the projections in the fourier domain and back projecting the inverse fourier transform of the result. A formal derivation of the FBP method uses the fourier-slice theorem as follows:

Thus can be obtained by back projection of . The filtered projection profile is obtained by applying the ramp filter , defined in the frequency space as .