Metal Artifact Reduction in Cone-Beam X-Ray CT via Ray Profile Correction

08/06/2018 ∙ by Sungsoo Ha, et al. ∙ Stony Brook University 0

In computed tomography (CT), metal implants increase the inconsistencies between the measured data and the linear attenuation assumption made by analytic CT reconstruction algorithms. The inconsistencies give rise to dark and bright bands and streaks in the reconstructed image, collectively called metal artifacts. These artifacts make it difficult for radiologists to render correct diagnostic decisions. We describe a data-driven metal artifact reduction (MAR) algorithm for image-guided spine surgery that applies to scenarios in which a prior CT scan of the patient is available. We tested the proposed method with two clinical datasets that were both obtained during spine surgery. Using the proposed method, we were not only able to remove the dark and bright streaks caused by the implanted screws but we also recovered the anatomical structures hidden by these artifacts. This results in an improved capability of surgeons to confirm the correctness of the implanted pedicle screw placements.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

X-ray computed tomography (CT) is a leading cross-sectional imaging technique lauded for its high image resolution and rapid speed of acquisition. CT reconstruction algorithms based on the Radon transform reconstruct the internal structures of human bodies by utilizing the X-ray photon interactions with matter which are following the Beer-Lambert law Kak and Slaney (1988). However, this procedure takes the assumption that the X-ray source is monochromatic. When it comes to (more realistic) polychromatic X-ray sources, low-energy photons are attenuated more easily than high-energy ones and, as a result, the X-ray beam is hardened as it passes through the material, shifting the energy spectrum toward higher energies. This becomes a particular problem when CT scanning is conducted with patients bearing metallic implants, which have dramatically increased attenuation properties for lower energies. The selective photon absorption not only increases the amount of dose absorbed by the object, it also amplifies the X-ray beam’s hardness. Moreover, the implanted metals can severely change photon direction (and energy). A failure to consider all of these kinds of non-linear behaviors in the X-ray photon interactions (also include Poisson noise, photon starvation, motion, partial volume effect and etc) results in various artifacts, for example, dark streaks along the lines of greatest attenuation Brooks and Di Chiro (1976); Hsieh (2009). The high pass filter used in Filtered-back projection (FBP) Feldkamp et al. (1984) then further exaggerates the differences between adjacent detector elements where one element has received a hardened beam and the other has not. This unintended contrast produces bright streaks in other directions. As a consequence, due to these adverse mechanisms, metal artifacts obscure information about anatomical structures, making it difficult for radiologists to correctly interpret the affected CT images.

There have been extensive efforts in developing metal artifact reduction (MAR) algorithms to compensate the approximation errors caused by implanted metals or high density objects. These efforts can be largely divided into two types of approaches: iterative reconstruction and sinogram correction. The iterative reconstruction methods adapt existing X-ray CT systems by incorporating one or more of the following types of a-priori knowledge: (1) low-level information of the images to be reconstructed Xue et al. (2009); Duan et al. (2008); Zhang et al. (2011); Mehranian et al. (2013), (2) the X-ray spectrum of the source De Man et al. (2001), (3) the attenuation functions of the base materials Elbakri and Fessler (2003); Srivastava and Fessler (2005); Abella and Fessler (2009), and (4) the composition of the metal components Stayman et al. (2012). More recently proposed iterative algorithms attempt to reduce beam hardening effects without the need of any prior knowledge by decomposing the image to be reconstructed into low and high density components Kyriakou et al. (2010); Jin et al. (2015).

On the other hand, the sinogram correction methods aim to directly correct the metal shadow regions in the projection data in which the corresponding rays have interacted with metal objects. One early approach replaces the corrupted data with their neighbors using linear Kalender et al. (1987)

or high-order interpolation schemes

Bazalova et al. (2007); Abdoli et al. (2010); Zhao et al. (2000). However, the interpolation-based MARs often suffer from loss of detail around the metal objects, and they also have high propensity to introduce new streak artifacts Müller and Buzug (2009). To address the lack of structural information, Prell et al. Prell et al. (2009) and Meyer et al. Meyer et al. (2010)

attempted to build prior CT images by roughly segmenting the uncorrected or pre-corrected CT image into soft-tissue, air, and bone equivalent materials. This has been a promising idea and further efforts have emerged that seek to produce a better prior image with the help of advanced computer vision techniques

Li et al. (2014); Karimi et al. (2012). Recently, for example, Zhang et al. Zhang and Yu (2017)

utilized a convolutional neural network to generate a more sophisticated prior and used it to correct the sinogram contaminated by metal artifacts.

In this paper we present a new MAR method that also uses the general approach of correcting a contaminated sinogram by substituting corrupted data with cleaner data available in prior images. The synthesis process we propose is not unlike the one often used in image-guided surgery (IGS) Grimson et al. (1999). These methods perform a real-time correlation of the operative field with a preoperative imaging dataset to show the precise location of a selected surgical instrument in the surrounding anatomic structures. To realize this, before the surgery, the patient undergoes a series of CT scans that reveal the soft tissues and bony structures. In our scenario, these preoperative CT images serve as prior images to help remove the metal artifacts that appear in intra- or post-surgery CT scans due to the implanted metal objects. Since such prior images have been acquired from the same patient, they will likely contain very similar internal structures, especially around metal implants. Furthermore, as these regions are often at least partially surrounded by bone it is unlikely that they are markedly deformed during the surgery. Thus, to find surrogate values to replace unreliable data in the sinogram, we first search ray paths in the prior images that have very similar density profiles along the ray passing through the metal objects. Then, the best matched prior ray profiles are used to correct the ray paths profiles that are corrupted by metal artifacts. Finally, the unreliable data are replaced with the re-projections of corrected ray profiles. We explored this general idea in Ha and Mueller (2016) but this preliminary work was limited to 2D fan-beam CT geometry. In this paper, we generalize our method to 3D cone-beam CT geometry and also present a significantly refined and mature framework.

In the following, Section 2 describes the proposed method and its technical details. Then, in Section 3, we show metal artifact reduction results. Section 4 concludes the paper with a discussion on future research directions for the proposed method.

2 Methods

In the following we use spine surgery as an example where a prior patient scan is available and an immediate post-surgery scan is required to confirm the correctness of the placements of the implanted metals screws.

Our method requires two CT data sets (or one that has a sufficient portion free of metal). One data set is an artifact-free prior CT scan taken before the (spine) surgery. The other is obtained during the surgery, containing metal artifacts due to the implanted pedicle screws. Although the two CT scans are taken from the same patient, they might be obtained in different conditions (e.g. patient’s pose, X-ray dosage amount, field-of-view and etc.). Therefore, it is necessary to register one volume to the other before applying the proposed ray profile correction method. This registration step is described in Section 2.1. Ray profile correction is also required to know which parts of a ray path belong to metal objects, and which ones do not. For this, we extract the implanted pedicle screws from the uncorrected CT volume (the volume suffering from metal artifacts). This metal localization step is described in Section 2.2. In the following we shall call a set of sample points along a ray the ray profile. A line integral is then computed as the weighted sum of all sample points of a ray profile. Also, we will define the regions traversed by rays passing through metal objects as metal shadow. The projection values under the metal shadow are unreliable because of beam-hardening, photon starvation, and so on, and they will result in metal artifacts. Our goal is to compute surrogate values in the metal shadow regions by correcting the corresponding ray profiles using the aligned prior CT volume and geometric information of the implanted metal (here the pedicle screws). This new correction scheme is explained in Section 2.3. Finally, the corrected metal shadow is combined smoothly with the original CT projection data as described in Section 2.4. The overall process is illustrated in Figure 1.

Figure 1: Overview of the proposed MAR method

2.1 Rigid volume registration

For the ray profile correction, we need to find matched prior profiles from the prior CT volume generated in the pre-operative CT scan. This prior CT volume is usually significantly misaligned with the CT volume obtained during or after the surgery. The patient may be in a different pose or the CT scan may cover a different range of the spine region (or have a different field-of-view). All that these two volumes might have in common is the surgical region itself. One naive approach for finding matched ray profiles would be to exhaustively search the prior CT volume with the ray profiles extracted from the uncorrected CT volume. This approach would be computationally very demanding as the searching space is almost infinite.

Instead, we align the two volumes and then find the set of matched prior profiles. The challenge in aligning two CT volumes taken at different times is that there can be large discrepancies in the soft internal structures (e.g. tissues). To resolve this problem, we first extract bone structures which are quite robust to deformation and also less contaminated by metal artifacts, in contrast to soft structures (see Figure 2 for a visualization). For the bone structure extraction, assuming there are only low and high density materials, we employ the balanced histogram thresholding (BHT) method Anjos and Shahbazkia (2008). Figure 3 shows two CT volumes obtained before and after surgery along with bone structures extracted using the BHT method. The prior volume is then rigidly registered to the uncorrected one by solving the following minimization problem:

(1)

Here, is the bone-only CT volume and its super-script and indicate the uncorrected and the prior CT volumes, respectively, which consist of voxels in total. The

is a rigid volume transformation operator with parameter vector,

, which includes translations and rotations in a 3D Cartesian coordinate system. Using Eq.(1) we find the optimal parameter vector, , which has a minimum weighted sum of the voxel-wise absolute difference between the two volumes. The weight term, , is there to further penalize the data mismatch term if two voxels came from different anatomical structures (materials). It is formalized as follows:

(2)

Here, we give more penalty, if the two voxels are not in same material (i.e. two voxels are from bone structures or not). This minimizes the contribution of mis-categorized voxels in the uncorrected volume, such as voxels in the bright band or in implanted metals that are regarded as bone after applying the BHT segmentation.

Figure 2: Similarity in bone structures between CT scans (a) before and (b) after a spine surgery. The visualizations are manually adjusted and colored to show the best view (Blue for bone and red for implanted pedicle screws and jaw that have higher density than bone).

We use a GPU-accelerated Hybrid-PSO (particle swarm optimization) algorithm to solve the minimization problem in Eq.(

1) which avoids a convergence to a local minima Kennedy (2011); Wachowiak et al. (2004)

. More specifically, in each generation, we randomly choose half of the particles and randomly adjust either a translation or rotation parameter with uniform probability. In every third generation, we pick half of the worst particles. The first half of these are replaced with completely new random values. Among the remainder, three-fifth of the particles are randomized as we do in each generation and the crossover is applied to the others. These types of randomization strategies have proven effective in finding the global solution in different optimization tasks

Wachowiak et al. (2004); Sharp et al. (2015). Figure 3 shows an example of the rigid registration result.

Figure 3: Rigid CT volume registration. (a) uncorrected CT volume, (b) prior CT volume, (c) registered prior CT volume to (a), (d) extracted bone structures of (a), (e) extracted bone structures of (b), (f) pseudo-colored overlap view of (a) and (c). Note that the presented CT slices, (a) and (b), initially show a large discrepancy although they have been sampled from the same z-index before applying the registration.

2.2 Localization of implanted metal objects

Before we can correct the ray profiles (see Section 2.3), we need to know whether a sample point in a given ray profile originates from metal (or not). For this purpose, we segment the screws from the uncorrected CT volume as follow. In the first step, we use the balanced histogram thresholding (BHT) method Anjos and Shahbazkia (2008) to get a coarse segmentation. After that we apply the DBSCAN (Density-based spatial clustering of applications with noise) algorithm Ester et al. (1996)

to the segmented structures. DBSCAN is a popular clustering algorithm which classifies points that are not well connected to a cluster as outliers. We found that DBSCAN did very well to remove any remaining noise and obtain an accurate segmentation of the metal objects (in our case, the screws). The clean segmentation also allows us to to precisely determine how many screws were implanted. Optionally, we might also include prior geometric knowledge to accelerate the process and to further improve the clustering accuracy. Figure

4 shows a screw extracted from an uncorrected CT volume. Figure 4c is what we will refer to as the metal-only CT volume.

Figure 4: Metal localization. (a) uncorrected CT volume, (b) extracted high density materials using the BHT, and (c) final results with DBSCAN

2.3 Ray profile correction

To compute the profile’s surrogate values (the values subject to replacement), we use the observation that metal artifacts usually appear around implanted metals and that the degree of corruption tends to decrease with distance from the metal region. Using this observation, the noisy ray profiles are corrected by dividing it into two regions, metal and non-metal regions. A metal region is the part of a profile that traverses a metal-only CT volume. For these regions, since we usually know the material of the implanted metals and their linear attenuation coefficients, the surrogate values are replaced with the linearly interpolated values of the two profiles extracted from the prior and the metal-only CT volume. We use linear interpolation to take into account the partial volume effect around the metal boundaries and so avoid any binarization artifacts. For the non-metal regions, the surrogate values are computed by linear interpolating between the noisy and prior profiles extracted from the uncorrected and the prior CT volume, respectively. The interpolation weight is given by the distance from the nearby metal boundaries along the ray path. It takes into account that when a point in a ray profile is close to a metal it is more likely deteriorated by metal artifacts and thus we put more emphasis on prior information. Vice versa, when a profile point is sufficiently far away from metal we can safely rely on its value in the uncorrected volume. As such, our method smoothly blends prior image information into the currently acquired imagery but only at locations where the current image information is likely unreliable due to metal artifacts.

Figure 5: Ray profile correction. [left] ray profiles extracted from noisy, prior, metal-only and corrected CT volumes and interpolation weights in [right-top] metal region and [right-bottom] non-metal region. In these figures the x-axis denotes the sampled locations of the ray profiles, while the y-axis shows the intensity (or weight) values at these positions.

Our ray profile correction scheme is described in Figure 5 and can be formulated as follows:

(3)

where is the linear interpolation operator such that . In this equation, represents the sampled value of a ray profile at position while its superscripts , , and , indicate the corrected profile and the profiles extracted from the metal-only, uncorrected, and aligned prior volumes, respectively. The superscript, , denotes the distance transform of and henceforth, is the distance from the position to the closest metal boundary along the ray path Maurer et al. (2003). The value is the linear attenuation coefficient of the implanted metals while is a scalar that controls the smoothness of the weight factor and is a constant value representing how similar the aligned prior volume is to the uncorrected one. We experimentally determined it as in this work.

2.4 Seamless in-painting

Recall that sample profiles are only computed for rays that traverse the segmented metal in the metal-only CT volume. We integrate these rays and store them in a corrected sinogram. Figure 6a shows a portion of an uncorrected sinogram while Figure 6b shows the same region with the corrected profiles only. The final task is to replace the metal shadow regions of the uncorrected sinogram with these corrected regions. However, a direct replacement of the data can lead to undesired discontinuities around the boundary of the metal shadow, resulting in the generation of new artifacts Müller and Buzug (2009). Therefore, it is important to seamlessly combine the new data with the existing ones at the boundary while internally keeping the relative contrast and the details of the data. To achieve such a seamless in-painting, we solve the following minimization problem Pérez et al. (2003):

(4)

where is projection data and its super-scripts, , and , represent the in-painted, corrected, and original projection data, respectively. denotes the metal shadow region and is 8-connected neighborhood of a pixel, . In this equation, the first term aims to preserve the gradients of the original (uncorrected) projection data in the metal shadow regions while lowering the intensities to non-metal values. Preserving the original gradients ensures that the detail and contrast in the projection data is maintained for the subsequent reconstruction. The second term in the equation affects the metal shadow region’s boundary only and ensures a smooth transition to the outside regions. Figure 6c shows a result of this seamless in-painting process using the same region than in panel a and b.

Figure 6: An example of seamless in-painting. (a) original, (b) corrected and (c) in-painted projection data

3 Results and Discussions

To test the proposed MAR framework, we used clinical CT projection data obtained during an image-guided surgery procedure on a cervical spine region using a Medtronic O-Arm surgical imaging CT scanner. The scanner has a source to axis distance of and a source to detector distance of . It is equipped with a flat X-ray detector with bins and an active area of

. During the scan, 360 projections were collected uniformly distributed over 360

. The 3D reconstruction used the filtered back-projection algorithm Kak and Slaney (1988) and produced a volume with a voxel size of .

Figures 7 and 8 show some results we obtained using our metal artifact reduction algorithm. The spine has two pedicle screws implanted. Figure 7 shows one of them in transverse, sagittal and coronal views. Note that the sagittal and coronal views are horizontal and vertical cut slices passing through the screws, respectively. Figure 8 shows the other implanted screws in the same manner. Overall, the proposed method effectively removes the metal artifacts (dark/bright bands and streaks) and reveals clear outlines of the implanted pedicle screws which are suitable for evaluating their placements after the surgery. For example, in Figure 7, top and bottom row, the yellow arrow indicates a pedicle screw where only the corrected image (column b) can reveal that is has been correctly inserted into the bone without extending into the tissue. Likewise, the yellow arrow in the middle row in Figure 7) shows a volume feature that was previously hidden by the beam hardening artifacts (column a) but is now readily visible.

Figure 7: Case study I. From top to bottom, transverse, coronal and sigittal views.

One side effect of the proposed method is the tendency of blurring the anatomical structures near metal objects. The difference images (column c) between the uncorrected and corrected images show (1) the removed artifacts, (2) a bright version of the metal pieces (as mentioned, our method lowers their projection values in the sinogram), and (3) some incorrectly removed details. The latter causes the blurring effects (annotated by the red arrow in Figure 7 and Figure 8). We think this is primarily because of the distance-based artifact region prediction model in Eq.(3

) where the model estimates the artifact regions based on the distance (

) from a point to the nearest metal boundary along a ray path regardless of whether the point is corrupted by metal artifacts or not. One way to mitigate the blurring effect could be to introduce an additional stage at the end of our MAR framework that would exploit the information hidden in the low- and high-pass filtered sinograms Jeong and Ra (2009) or images Meyer et al. (2012) to control the correction process.

Figure 8: Case study II. From top to bottom, transverse, coronal and sagittal views.

4 Conclusion

We have presented a new method for metal artifact reduction (MAR). It assumes that a prior CT scan taken before implanting the metal objects into the patient is available. Using this prior scan and a segmentation or a model of the metal implant, we utilize a novel ray profile correction scheme which computes an accurate estimate of the rays traversing the projection regions affected by the metal artifacts. Our experiments with clinical CT data indicate that the proposed method can clearly reveal the placements of implanted pedicle screws that were ambiguous before due the significant metal artifacts.

Future work will extend this work along three directions. Firstly, we will further investigate the behavior of the phenomenon of beam hardening for better metal artifact region prediction. Secondly, using a CAD model available for commercial pedicle screws, we believe a better localization of the implanted metal can be accomplished, leading to a better ray profile correction. Lastly, we also believe that utilizing existing algorithms that can reveal features hidden by metal artifacts in conjunction with a MAR method will be able to help resolve the occasional blurring effects around metals. Finally, we also plan to conduct more clinical case studies which will contain a wide selection of different implanted metals and different amounts of metal artifacts.

Acknowledgement

This research was partially supported by NSF grant IIS 1527200 and the Ministry of Science, ICT and Future Planning, Korea, under the “IT Consilience Creative Program (ITCCP)” supervised by NIPA. We also thank Medtronic, Inc. for the data and support.

References

References

References

  • Kak and Slaney (1988) Avinash C.. Kak and Malcolm Slaney, Principles of computerized tomographic imaging (IEEE press, 1988).
  • Brooks and Di Chiro (1976) Rodney A Brooks and Giovanni Di Chiro, “Beam hardening in x-ray reconstructive tomography,” Physics in medicine and biology 21, 390 (1976).
  • Hsieh (2009) Jiang Hsieh, “Computed tomography: principles, design, artifacts, and recent advances,”  (SPIE Bellingham, WA, 2009).
  • Feldkamp et al. (1984) LA Feldkamp, LC Davis,  and JW Kress, “Practical cone-beam algorithm,” JOSA A 1, 612–619 (1984).
  • Xue et al. (2009) Hui Xue, Li Zhang, Yongshun Xiao, Zhiqiang Chen,  and Yuxiang Xing, “Metal artifact reduction in dual energy ct by sinogram segmentation based on active contour model and tv inpainting,” in 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC) (IEEE, 2009) pp. 904–908.
  • Duan et al. (2008) Xinhui Duan, Li Zhang, Yongshun Xiao, Jianping Cheng, Zhiqiang Chen,  and Yuxiang Xing, “Metal artifact reduction in ct images by sinogram tv inpainting,” in 2008 IEEE Nuclear Science Symposium Conference Record (IEEE, 2008) pp. 4175–4177.
  • Zhang et al. (2011) Yi Zhang, Yi-Fei Pu, Jin-Rong Hu, Yan Liu,  and Ji-Liu Zhou, “A new ct metal artifacts reduction algorithm based on fractional-order sinogram inpainting,” Journal of X-ray science and technology 19, 373–384 (2011).
  • Mehranian et al. (2013) Abolfazl Mehranian, Mohammad Reza Ay, Arman Rahmim,  and Habib Zaidi, “3d prior image constrained projection completion for x-ray ct metal artifact reduction,” IEEE Transactions on Nuclear Science 60, 3318–3332 (2013).
  • De Man et al. (2001) Bruno De Man, Johan Nuyts, Patrick Dupont, Guy Marchal,  and Paul Suetens, “An iterative maximum-likelihood polychromatic algorithm for ct,” IEEE transactions on medical imaging 20, 999–1008 (2001).
  • Elbakri and Fessler (2003) Idris A Elbakri and Jeffrey A Fessler, “Segmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation,” Physics in medicine and biology 48, 2453 (2003).
  • Srivastava and Fessler (2005) Somesh Srivastava and Jeffrey A Fessler, “Simplified statistical image reconstruction algorithm for polyenergetic x-ray ct,” in IEEE Nuclear Science Symposium Conference Record, Vol. 3 (Citeseer, 2005) pp. 1551–1555.
  • Abella and Fessler (2009) Mónica Abella and Jeffrey A Fessler, “A new statistical image reconstruction algorithm for polyenergetic x-ray ct,” in 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro (IEEE, 2009) pp. 165–168.
  • Stayman et al. (2012) J Webster Stayman, Yoshito Otake, Jerry L Prince, A Jay Khanna,  and Jeffrey H Siewerdsen, “Model-based tomographic reconstruction of objects containing known components,” IEEE transactions on medical imaging 31, 1837–1848 (2012).
  • Kyriakou et al. (2010) Yiannis Kyriakou, Esther Meyer, Daniel Prell,  and Marc Kachelrieß, “Empirical beam hardening correction (ebhc) for ct,” Medical physics 37, 5179–5187 (2010).
  • Jin et al. (2015) Pengchong Jin, Charles A Bouman,  and Ken D Sauer, “A model-based image reconstruction algorithm with simultaneous beam hardening correction for x-ray ct,” IEEE Transactions on Computational Imaging 1, 200–216 (2015).
  • Kalender et al. (1987) Willi A Kalender, Robert Hebel,  and Johannes Ebersberger, “Reduction of ct artifacts caused by metallic implants.” Radiology 164, 576–577 (1987).
  • Bazalova et al. (2007) Magdalena Bazalova, Luc Beaulieu, Steven Palefsky,  and Frank Verhaegen, “Correction of ct artifacts and its influence on monte carlo dose calculations,” Medical physics 34, 2119–2132 (2007).
  • Abdoli et al. (2010)

    Mehrsima Abdoli, Mohammad Reza Ay, Alireza Ahmadian, Rudi AJO Dierckx,  and Habib Zaidi, “Reduction of dental filling metallic artifacts in ct-based attenuation correction of pet data using weighted virtual sinograms optimized by a genetic algorithm,” Medical physics 

    37, 6166–6177 (2010).
  • Zhao et al. (2000) Shiying Zhao, DD Robeltson, Ge Wang, Bruce Whiting,  and Kyongtae T Bae, “X-ray ct metal artifact reduction using wavelets: an application for imaging total hip prostheses,” IEEE transactions on medical imaging 19, 1238–1247 (2000).
  • Müller and Buzug (2009) Jan Müller and Thorsten M Buzug, “Spurious structures created by interpolation-based ct metal artifact reduction,” in SPIE Medical Imaging (International Society for Optics and Photonics, 2009) pp. 72581Y–72581Y.
  • Prell et al. (2009) Daniel Prell, Yiannis Kyriakou, Marcel Beister,  and Willi A Kalender, “A novel forward projection-based metal artifact reduction method for flat-detector computed tomography,” Physics in medicine and biology 54, 6575 (2009).
  • Meyer et al. (2010) Esther Meyer, Rainer Raupach, Michael Lell, Bernhard Schmidt,  and Marc Kachelrieß, “Normalized metal artifact reduction (nmar) in computed tomography,” Medical physics 37, 5482–5493 (2010).
  • Li et al. (2014) Ming Li, Zhaobang Liu, Yuefang Dong, Kai Su,  and Kangxin Luo, “A prior-interpolation based metal artifact reduction algorithm in computed tomography,” in 2014 7th International Conference on Biomedical Engineering and Informatics (IEEE, 2014) pp. 24–28.
  • Karimi et al. (2012) Seemeen Karimi, Pamela Cosman, Christoph Wald,  and Harry Martz, “Segmentation of artifacts and anatomy in ct metal artifact reduction,” Medical physics 39, 5857–5868 (2012).
  • Zhang and Yu (2017) Yanbo Zhang and Hengyong Yu, “Convolutional neural network based metal artifact reduction in x-ray computed tomography,” arXiv preprint arXiv:1709.01581  (2017).
  • Grimson et al. (1999) WEL Grimson, RJFA Kikinis, Ferenc A Jolesz,  and PM Black, “Image-guided surgery,” Scientific American 280, 54–61 (1999).
  • Ha and Mueller (2016) Sungsoo Ha and Klaus Mueller, “Metal artifact reduction in ct via ray profile correction,” in SPIE Medical Imaging (International Society for Optics and Photonics, 2016) pp. 978334–978334.
  • Anjos and Shahbazkia (2008) António Anjos and Hamid Reza Shahbazkia, “Bi-level image thresholding - a fast method,” BIOSIGNALS 2, 70–76 (2008).
  • Kennedy (2011) James Kennedy, “Particle swarm optimization,” in 

    Encyclopedia of machine learning

     (Springer, 2011) pp. 760–766.
  • Wachowiak et al. (2004)

    Mark P Wachowiak, Renata Smolíková, Yufeng Zheng, Jacek M Zurada,  and Adel Said Elmaghraby, “An approach to multimodal biomedical image registration utilizing particle swarm optimization,” IEEE Transactions on evolutionary computation 

    8, 289–301 (2004).
  • Sharp et al. (2015) Toby Sharp, Cem Keskin, Duncan Robertson, Jonathan Taylor, Jamie Shotton, David Kim, Christoph Rhemann, Ido Leichter, Alon Vinnikov, Yichen Wei, et al., “Accurate, robust, and flexible real-time hand tracking,” in Proceedings of the 33rd Annual ACM Conference on Human Factors in Computing Systems (ACM, 2015) pp. 3633–3642.
  • Ester et al. (1996) Martin Ester, Hans-Peter Kriegel, Jörg Sander, Xiaowei Xu, et al., “A density-based algorithm for discovering clusters in large spatial databases with noise.” in Kdd, Vol. 96 (1996) pp. 226–231.
  • Maurer et al. (2003) Calvin R Maurer, Rensheng Qi,  and Vijay Raghavan, “A linear time algorithm for computing exact euclidean distance transforms of binary images in arbitrary dimensions,” IEEE Transactions on Pattern Analysis and Machine Intelligence 25, 265–270 (2003).
  • Pérez et al. (2003) Patrick Pérez, Michel Gangnet,  and Andrew Blake, “Poisson image editing,” in ACM Transactions on graphics (TOG), Vol. 22 (ACM, 2003) pp. 313–318.
  • Jeong and Ra (2009) Kye Young Jeong and Jong Beom Ra, “Metal artifact reduction based on sinogram correction in ct,” in 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC) (IEEE, 2009) pp. 3480–3483.
  • Meyer et al. (2012) Esther Meyer, Rainer Raupach, Michael Lell, Bernhard Schmidt,  and Marc Kachelrieß, “Frequency split metal artifact reduction (fsmar) in computed tomography,” Medical physics 39, 1904–1916 (2012).