<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.0 20120330//EN" "JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.0" xml:lang="en" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">pmp</journal-id>
<journal-title-group>
<journal-title>Progress in Medical Physics</journal-title>
<abbrev-journal-title>Progress in Medical Physics</abbrev-journal-title>
</journal-title-group>
<issn pub-type="ppub">2508-4445</issn>
<issn pub-type="epub">2508-4453</issn>
<publisher>
<publisher-name>Korean Society of Medical Physics</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.14316/pmp.2013.24.4.303</article-id>
<article-id pub-id-type="publisher-id">pmp-24-303</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Theoretical Investigation of Metal Artifact Reduction Based on Sinogram Normalization in Computed Tomography</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name name-style="western" xml:lang="en">
<surname>Jeon</surname><given-names>Hosang</given-names></name>
<xref ref-type="aff" rid="aff1-pmp-24-303"><sup>&#x2217;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name name-style="western" xml:lang="en">
<surname>Youn</surname><given-names>Hanbean</given-names></name>
<xref ref-type="aff" rid="aff2-pmp-24-303"><sup>&#x2020;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name name-style="western" xml:lang="en">
<surname>Nam</surname><given-names>Jiho</given-names></name>
<xref ref-type="aff" rid="aff1-pmp-24-303"><sup>&#x2217;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name name-style="western" xml:lang="en">
<surname>Kim</surname><given-names>Ho Kyung</given-names></name>
<xref ref-type="aff" rid="aff2-pmp-24-303"><sup>&#x2020;</sup></xref>
<xref ref-type="aff" rid="aff3-pmp-24-303"><sup>&#x2021;</sup></xref>
</contrib>
<aff id="aff1-pmp-24-303" xml:lang="en"><label>&#x2217;</label>Department of Radiation Oncology, Pusan National University Yangsan Hospital, Yangsan, <country>Korea</country></aff>
<aff id="aff2-pmp-24-303" xml:lang="en"><label>&#x2020;</label>School of Mechanical Engineering, Pusan National University, Busan, <country>Korea</country></aff>
<aff id="aff3-pmp-24-303" xml:lang="en"><label>&#x2021;</label>Advanced Medical Engineering Research Center, Pusan National University, Busan, <country>Korea</country></aff>
</contrib-group>
<pub-date pub-type="ppub"><month>12</month><year>2013</year></pub-date>
<pub-date pub-type="epub"><day>17</day><month>12</month><year>2013</year></pub-date>
<volume>24</volume><issue>4</issue>
<fpage>303</fpage><lpage>314</lpage>
<permissions>
<copyright-statement>Copyright &#x00A9; 2013 Korean Society of Medical Physics</copyright-statement>
<copyright-year>2013</copyright-year>
<license><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/3.0/">http://creativecommons.org/licenses/by-nc/3.0/</ext-link>) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p></license>
</permissions>
<abstract xml:lang="en">
<title>Abstract</title>
<p>Image quality of computed tomography (CT) is very vulnerable to metal artifacts. Recently, the thickness and background normalization techniques have been introduced. Since they provide flat sinograms, it is easy to determine metal traces and a simple linear interpolation would be enough to describe the missing data in sinograms. In this study, we have developed a theory describing two normalization methods and compared two methods with respect to various sizes and numbers of metal inserts by using simple numerical simulations. The developed theory showed that the background normalization provide flatter sinograms than the thickness normalization, which was validated with the simulation results. Numerical simulation results with respect to various sizes and numbers of metal inserts showed that the background normalization was better than the thickness normalization for metal artifact corrections. Although the residual artifacts still existed, we have showed that the background normalization without the segmentation procedure was better than the thickness normalization for metal artifact corrections. Since the background normalization without the segmentation procedure is simple and it does not require any users&#x0027; intervention, it can be readily installed in conventional CT systems.</p>
</abstract>
<kwd-group xml:lang="en">
<kwd>Computed tomography</kwd>
<kwd>Metal artifact</kwd>
<kwd>Sinogram</kwd>
<kwd>Linear interpolation</kwd>
</kwd-group>
</article-meta>
</front>
<back>
<ref-list xml:lang="en">
<title>REFERENCES</title>
<ref id="b1"><label>1.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parker</surname> <given-names>RP</given-names></name> <name><surname>Hobday</surname> <given-names>PA</given-names></name> <name><surname>Cassell</surname> <given-names>KJ</given-names></name></person-group> <article-title>The direct use of CT numbers in radiotherapy dosage calculations for inhomogeneous media</article-title> <source>Phys Med Biol</source> <volume>24</volume><issue>4</issue> <fpage>802</fpage><lpage>809</lpage> <year>1979</year></element-citation></ref>
<ref id="b2"><label>2.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guerrero</surname> <given-names>ME</given-names></name> <name><surname>Jacobs</surname> <given-names>R</given-names></name> <name><surname>Loubele</surname> <given-names>M</given-names></name> <etal/></person-group> <article-title>State-ofthe- art on cone beam CT imaging for preoperative planning of implant placement</article-title> <source>Clin Oral Invest</source> <volume>10</volume><issue>1</issue> <fpage>1</fpage><lpage>7</lpage> <year>2006</year></element-citation></ref>
<ref id="b3"><label>3.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cho</surname> <given-names>MK</given-names></name> <name><surname>Kim</surname> <given-names>HK</given-names></name> <name><surname>Youn</surname> <given-names>H</given-names></name> <name><surname>Kim</surname> <given-names>SS</given-names></name></person-group> <article-title>A feasibility study of digital tomosynthesis for volumetric dental imaging</article-title> <source>J Instrum</source> <volume>7</volume><issue>3</issue> <fpage>P03007</fpage> <year>2012</year></element-citation></ref>
<ref id="b4"><label>4.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Glover</surname> <given-names>GH</given-names></name> <name><surname>Pelc</surname> <given-names>NJ</given-names></name></person-group> <article-title>An algorithm for the reduction of metal clip artifacts in CT reconstructions</article-title> <source>Med Phys</source> <volume>8</volume><issue>6</issue> <fpage>799</fpage><lpage>807</lpage> <year>1981</year></element-citation></ref>
<ref id="b5"><label>5.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kalender</surname> <given-names>WA</given-names></name> <name><surname>Hebel</surname> <given-names>R</given-names></name> <name><surname>Ebersberger</surname> <given-names>J</given-names></name></person-group> <article-title>Reduction of CT artifacts caused by metallic implants</article-title> <source>Radiology</source> <volume>164</volume><issue>2</issue> <fpage>576</fpage><lpage>577</lpage> <year>1987</year></element-citation></ref>
<ref id="b6"><label>6.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>G</given-names></name> <name><surname>Snyder</surname> <given-names>DL</given-names></name> <name><surname>O&#x0027;Sullivan</surname> <given-names>JA</given-names></name> <etal/></person-group> <article-title>Iterative deblurring for CT metal artifact reduction</article-title> <source>IEEE Trans Med Imaging</source> <volume>15</volume><issue>5</issue> <fpage>657</fpage><lpage>664</lpage> <year>1996</year></element-citation></ref>
<ref id="b7"><label>7.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>De Man</surname> <given-names>B</given-names></name> <name><surname>Nuyts</surname> <given-names>J</given-names></name> <name><surname>Dupont</surname> <given-names>P</given-names></name> <etal/></person-group> <article-title>An iterative maximum- likelihood polychromatic algorithm for CT</article-title> <source>IEEE Trans Med Imaging</source> <volume>20</volume><issue>10</issue> <fpage>999</fpage><lpage>1008</lpage> <year>2001</year></element-citation></ref>
<ref id="b8"><label>8.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kachelrie&#x03B2;</surname> <given-names>M</given-names></name> <name><surname>Watzke</surname> <given-names>O</given-names></name> <name><surname>Kalender</surname> <given-names>WA</given-names></name></person-group> <article-title>Generalized multi-dimensional adaptive filtering for conventional and spiral single-slice, multislice, and cone-beam CT</article-title> <source>Med Phys</source> <volume>28</volume><issue>4</issue> <fpage>475</fpage><lpage>490</lpage> <year>2001</year></element-citation></ref>
<ref id="b9"><label>9.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nuyts</surname> <given-names>J</given-names></name> <name><surname>De Man</surname> <given-names>B</given-names></name> <name><surname>Fessler</surname> <given-names>JA</given-names></name> <etal/></person-group> <article-title>Modelling the physics in the iterative reconstruction for transmission computed tomography</article-title> <source>Phys Med Biol</source> <volume>58</volume><issue>12</issue> <fpage>R63</fpage><lpage>R96</lpage> <year>2013</year></element-citation></ref>
<ref id="b10"><label>10.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>JC</given-names></name> <name><surname>Song</surname> <given-names>B</given-names></name> <name><surname>Kim</surname> <given-names>JS</given-names></name> <etal/></person-group> <article-title>Fast compressed sensing- based CBCT reconstruction using Barzilai-Borwein formulation for application to online IGRT</article-title> <source>Med Phys</source> <volume>39</volume> <fpage>1207</fpage><lpage>1217</lpage> <year>2012</year></element-citation></ref>
<ref id="b11"><label>11.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Choi</surname> <given-names>J</given-names></name> <name><surname>Kim</surname> <given-names>KS</given-names></name> <name><surname>Kim</surname> <given-names>MW</given-names></name> <etal/></person-group> <article-title>Sparsity driven metal part reconstruction for artifact removal in dental CT</article-title> <source>J X-ray Sci Tech</source> <volume>19</volume><issue>4</issue> <fpage>457</fpage><lpage>475</lpage> <year>2011</year></element-citation></ref>
<ref id="b12"><label>12.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>M&#x00FC;ller</surname> <given-names>J</given-names></name> <name><surname>Buzug</surname> <given-names>TM</given-names></name></person-group> <article-title>Spurious structures created by interpolation- based CT metal artifact reduction</article-title> <source>Proc SPIE</source> <volume>7258</volume> <fpage>72581Y(8pp)</fpage> <year>2009</year></element-citation></ref>
<ref id="b13"><label>13.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>M&#x00FC;ller</surname> <given-names>J</given-names></name> <name><surname>Buzug</surname> <given-names>TM</given-names></name></person-group> <source>Intersection Line Length Normalization in CT Projection Data</source> <publisher-name>Bildverarbeitung fur die Medizin</publisher-name> <publisher-loc>Berlin</publisher-loc> <year>2008</year> <fpage>77</fpage><lpage>81</lpage></element-citation></ref>
<ref id="b14"><label>14.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meyer</surname> <given-names>E</given-names></name> <name><surname>Raupach</surname> <given-names>R</given-names></name> <name><surname>Lell</surname> <given-names>M</given-names></name> <name><surname>Schmidt</surname> <given-names>B</given-names></name> <name><surname>Kachelrie&#x03B2;</surname> <given-names>M</given-names></name></person-group> <article-title>Normalized metal artifact reduction (NMAR) in computed tomography</article-title> <source>Med Phys</source> <volume>37</volume><issue>10</issue> <fpage>5482</fpage><lpage>5493</lpage> <year>2010</year></element-citation></ref>
<ref id="b15"><label>15.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Jain</surname> <given-names>AK</given-names></name></person-group> <source>Fundamentals of Digital Image Processing</source> <publisher-name>Prentice-Hall Inc</publisher-name> <publisher-loc>Englewood Cliffs, NJ</publisher-loc> <year>1989</year><comment>), pp.</comment> <fpage>438</fpage><lpage>439</lpage></element-citation></ref>
<ref id="b16"><label>16.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Otsu</surname> <given-names>N</given-names></name></person-group> <article-title>A threshold selection method from gray-level histograms</article-title> <source>IEEE Trans Syst Man Cybern</source> <volume>9</volume><issue>1</issue> <fpage>62</fpage><lpage>66</lpage> <year>1979</year></element-citation></ref>
<ref id="b17"><label>17.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sezgin</surname> <given-names>M</given-names></name> <name><surname>Sankur</surname> <given-names>B</given-names></name></person-group> <article-title>Survey over image thresholding techniques and quantitative performance evaluation</article-title> <source>J Electron Imaging</source> <volume>13</volume><issue>1</issue> <fpage>146</fpage><lpage>165</lpage> <year>2004</year></element-citation></ref>
</ref-list>
<sec sec-type="display-objects">
<title>Figures and Tables</title>
<fig id="F1-pmp-24-303" position="float">
<label>Fig. 1.</label>
<caption xml:lang="en"><p>Sketch describing the projection of an object <italic>f</italic>(<italic>x, y</italic>) containing a circular metal insert by a single ray at a projection angle &#x03B8;. The lengths intersected by the ray with the object and the metal insert are <italic>L</italic>(<italic>t</italic>,&#x03B8;) and <italic>l</italic>(<italic>t</italic>,&#x03B8;), respectively.</p></caption>
<graphic xlink:href="pmp-24-303f1.tif"/>
</fig>
<fig id="F2-pmp-24-303" position="float">
<label>Fig. 2.</label>
<caption xml:lang="en"><p>Illustration of the metal artifact reduction algorithm based on the thickness and background normalizations. 1) Generation of a sinogram from the CT scan. 2) Image reconstruction using the sinogram 1). Generation of sinograms 3) and 3&#x0027;) by forward-projecting the uncorrected image after replacing its internal components by air and redefining the internal components, except the metal insert, by the segmentation procedure, respectively. 4) Extraction of metal traces from the original sinogram. Generation of normalized sinograms 5) and 5&#x0027;) dividing 1) by 3) and 3&#x0027;), respectively. Interpolated sinograms 6) and 6&#x0027;) by applying the linear interpolation to 5) and 5&#x0027;), respectively, with the known metal traces. Generation of corrected sinograms 7) and 7&#x0027;) by multiplying 6) and 6&#x0027;) by 3) and 3&#x0027;), respectively (the denormalization procedure). Generation of corrected images 8) and 8&#x0027;).</p></caption>
<graphic xlink:href="pmp-24-303f2.tif"/>
</fig>
<fig id="F3-pmp-24-303" position="float">
<label>Fig. 3.</label>
<caption xml:lang="en"><p>(a) Simple circular numerical phantom consisting of water (&#x03BC;<italic><sub>water</sub></italic>), bone (&#x03BC;<italic><sub>bone</sub></italic>) and gold (&#x03BC;<italic><sub>Au</sub></italic>). Radii of each material <italic>r<sub>j</sub></italic>, where <italic>j</italic>=water, bone, and Au, are 200, 100 and 20 mm, respectively. The center of gold material is located at <italic>t</italic><sub>0</sub>=120 mm. (b) Comparison of normalized sinograms obtained at &#x03B8;=0<sup>o</sup>. Symbols describe the simulation results, whereas solid lines describe theoretical calculations.</p></caption>
<graphic xlink:href="pmp-24-303f3.tif"/>
</fig>
<fig id="F4-pmp-24-303" position="float">
<label>Fig. 4.</label>
<caption xml:lang="en"><p>(a) Interpolated sinograms with various functions which include simple linear function, 2nd-, 3rd-, and 4th-order Lagrange functions, and cubic spline. For comparison, the thickness- and background-normalized sinograms with linear interpolation are also plotted. (b) Sinograms denormalized from the thickness and background normalizations. For comparison, the original sinogram with linear interpolation is also plotted. All the profiles are extracted from sinograms obtained at &#x03B8;=263<sup>o</sup> at which the effect of interpolation can be well addressed.</p></caption>
<graphic xlink:href="pmp-24-303f4.tif"/>
</fig>
<fig id="F5-pmp-24-303" position="float">
<label>Fig. 5.</label>
<caption xml:lang="en"><p>Restored images by using various interpolation methods for (a) an original image: (b) linear interpolation, (c) 2nd-order Lagrange interpolation, (d) 3rd-order Lagrange interpolation, (e) 4th-order Lagrange interpolation, (f) cubic spline interpolation, (g) linear interpolation in the thickness-normalized sinogram, and (h) linear interpolation in the background-normalized sinogram.</p></caption>
<graphic xlink:href="pmp-24-303f5.tif"/>
</fig>
<fig id="F6-pmp-24-303" position="float">
<label>Fig. 6.</label>
<caption xml:lang="en"><p>Column (a) shows the reconstructed images with metal objects having various sizes. Column (b) and (c) shows the restored images using the thickness and background normalization methods, respectively. Column (d) shows the restored images using the background normalization method without the segmentation procedure. The insets are the profiles extracted from sinograms obtained at &#x03B8;=30<sup>o</sup>.</p></caption>
<graphic xlink:href="pmp-24-303f6.tif"/>
</fig>
<fig id="F7-pmp-24-303" position="float">
<label>Fig. 7.</label>
<caption xml:lang="en"><p>Column (a) shows the reconstructed images with various numbers of metal objects. Column (b) and (c) shows the restored images using the thickness and background normalization methods, respectively. Column (d) shows the restored images using the background normalization method without the segmentation procedure. The insets are the profiles extracted from sinograms obtained at &#x03B8;=30<sup>o</sup>.</p></caption>
<graphic xlink:href="pmp-24-303f7.tif"/>
</fig>
<table-wrap id="T1-pmp-24-303" position="float">
<label>Table 1.</label>
<caption xml:lang="en"><p>Linear attenuation coefficients of tissue and metal, and ratios between metal and tissue at diagnostic energy.</p></caption>
<table frame="hsides" rules="all">
<thead>
<tr>
<th valign="middle" align="center" rowspan="2">Energy (keV)</th>
<th valign="middle" align="center" colspan="3">&#x03BC;<italic><sub>b</sub></italic> (cm<sup>&#xFF0D;1</sup>)</th>
<th valign="middle" align="center" colspan="2">&#x03BC;<italic><sub>m</sub></italic> (cm<sup>&#xFF0D;1</sup>)</th>
<th valign="middle" align="center" colspan="2">&#x03BC;<italic><sub>m</sub></italic>/&#x03BC;<italic><sub>b</sub></italic></th>
</tr>
<tr>
<th valign="middle" align="left">Bone &#x03BC;<italic><sub>bone</sub></italic></th>
<th valign="middle" align="center">Tissue &#x03BC;<italic><sub>tissue</sub></italic></th>
<th valign="middle" align="center">Water &#x03BC;<italic><sub>water</sub></italic></th>
<th valign="middle" align="center">Tungsten &#x03BC;<italic><sub>W</sub></italic></th>
<th valign="middle" align="center">Gold &#x03BC;<italic><sub>Au</sub></italic></th>
<th valign="middle" align="center">&#x03BC;<italic><sub>W</sub></italic>/&#x03BC;<italic><sub>water</sub></italic></th>
<th valign="middle" align="center">&#x03BC;<italic><sub>Au</sub></italic>/&#x03BC;<italic><sub>water</sub></italic></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">10</td>
<td valign="middle" align="center">54.74</td>
<td valign="middle" align="center">5.70</td>
<td valign="middle" align="center">5.33</td>
<td valign="middle" align="center">1,870.17</td>
<td valign="middle" align="center">2,281.69</td>
<td valign="middle" align="center">351</td>
<td valign="middle" align="center">428</td>
</tr>
<tr>
<td valign="middle" align="left">20</td>
<td valign="middle" align="center">7.68</td>
<td valign="middle" align="center">0.87</td>
<td valign="middle" align="center">0.81</td>
<td valign="middle" align="center">1,268.40</td>
<td valign="middle" align="center">1,522.61</td>
<td valign="middle" align="center">1,567</td>
<td valign="middle" align="center">1,881</td>
</tr>
<tr>
<td valign="middle" align="left">30</td>
<td valign="middle" align="center">2.56</td>
<td valign="middle" align="center">0.40</td>
<td valign="middle" align="center">0.38</td>
<td valign="middle" align="center">438.69</td>
<td valign="middle" align="center">531.69</td>
<td valign="middle" align="center">1,168</td>
<td valign="middle" align="center">1,416</td>
</tr>
<tr>
<td valign="middle" align="left">40</td>
<td valign="middle" align="center">1.28</td>
<td valign="middle" align="center">0.28</td>
<td valign="middle" align="center">0.27</td>
<td valign="middle" align="center">205.93</td>
<td valign="middle" align="center">250.77</td>
<td valign="middle" align="center">768</td>
<td valign="middle" align="center">935</td>
</tr>
<tr>
<td valign="middle" align="left">50</td>
<td valign="middle" align="center">0.81</td>
<td valign="middle" align="center">0.24</td>
<td valign="middle" align="center">0.23</td>
<td valign="middle" align="center">114.84</td>
<td valign="middle" align="center">140.22</td>
<td valign="middle" align="center">506</td>
<td valign="middle" align="center">618</td>
</tr>
<tr>
<td valign="middle" align="left">60</td>
<td valign="middle" align="center">0.60</td>
<td valign="middle" align="center">0.22</td>
<td valign="middle" align="center">0.21</td>
<td valign="middle" align="center">71.66</td>
<td valign="middle" align="center">87.50</td>
<td valign="middle" align="center">348</td>
<td valign="middle" align="center">425</td>
</tr>
<tr>
<td valign="middle" align="left">80</td>
<td valign="middle" align="center">0.43</td>
<td valign="middle" align="center">0.19</td>
<td valign="middle" align="center">0.18</td>
<td valign="middle" align="center">150.71</td>
<td valign="middle" align="center">42.21</td>
<td valign="middle" align="center">820</td>
<td valign="middle" align="center">230</td>
</tr>
<tr>
<td valign="middle" align="left">100</td>
<td valign="middle" align="center">0.36</td>
<td valign="middle" align="center">0.18</td>
<td valign="middle" align="center">0.17</td>
<td valign="middle" align="center">85.63</td>
<td valign="middle" align="center">99.65</td>
<td valign="middle" align="center">502</td>
<td valign="middle" align="center">584</td>
</tr>
<tr>
<td valign="middle" align="left">150</td>
<td valign="middle" align="center">0.28</td>
<td valign="middle" align="center">0.16</td>
<td valign="middle" align="center">0.15</td>
<td valign="middle" align="center">30.51</td>
<td valign="middle" align="center">35.92</td>
<td valign="middle" align="center">203</td>
<td valign="middle" align="center">239</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</back>
</article>