Next Article in Journal
An Iterative Algorithm for Semisupervised Classification of Hotspots on Bone Scintigraphies of Patients with Prostate Cancer
Next Article in Special Issue
Micro-CT for Biological and Biomedical Studies: A Comparison of Imaging Techniques
Previous Article in Journal
Investigating Semantic Augmentation in Virtual Environments for Image Segmentation Using Convolutional Neural Networks
Previous Article in Special Issue
Imaging with Coherent X-rays: From the Early Synchrotron Tests to SYNAPSE
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Suppression of Cone-Beam Artefacts with Direct Iterative Reconstruction Computed Tomography Trajectories (DIRECTT)

1
Bundesanstalt für Materialforschung und -prüfung (BAM), Unter den Eichen 87, 12205 Berlin, Germany
2
Institute of Physics and Astronomy, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
*
Author to whom correspondence should be addressed.
Submission received: 5 July 2021 / Revised: 11 August 2021 / Accepted: 13 August 2021 / Published: 15 August 2021
(This article belongs to the Special Issue X-ray Digital Radiography and Computed Tomography)

Abstract

:
The reconstruction of cone-beam computed tomography data using filtered back-projection algorithms unavoidably results in severe artefacts. We describe how the Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) algorithm can be combined with a model of the artefacts for the reconstruction of such data. The implementation of DIRECTT results in reconstructed volumes of superior quality compared to the conventional algorithms.

1. Introduction

Computed tomography (CT) systems for non-destructive testing and material analysis generally use a cone beam on a sample that rotates in a circular orbit [1], with cylindrical samples being among the most common [1,2,3,4]. The exact reconstruction of data acquired during such a measurement is not possible because the geometry does not satisfy Tuy’s sufficiency condition [5]. This is demonstrated in Figure 1 with the reconstruction of a concrete rod by the commonly used algorithm developed by Feldkamp, Davis and Kress (FDK) [6]. For higher cone angles, there is a decrease of the grey values, which represent the attenuation coefficient μ, and of the image quality in the direction of the rotation axis (z-axis in Figure 1).
Several algorithms have been proposed to reduce such artefacts. Hsieh proposed a two-pass algorithm that estimates the cone-beam artefacts from the segmented high-density material and then subtracts them from the FDK reconstruction [7]. Han and Baek went further by devising a multi-pass approach that they tested for larger cone angles and different material densities [8]. Maaß et al. proposed an iterative algorithm that also subtracts the estimated artefacts from the FDK reconstruction without requiring segmentation [9].
Here, we will describe how we have adjusted the Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) algorithm [10,11,12] to estimate such artefacts and compensate for them.

2. Materials and Methods

2.1. Sample Images

A set of 3000 cone-beam projections of the concrete rod of Figure 1 was acquired over 360° on an in-house GE v|tome|x L 300 scanner. A 2024 × 2024 PerkinElmer detector with a pixel size of 0.2 mm was used. Source-object and source-detector distances of 81 mm and 1018 mm, respectively, resulted in a magnification of 12.5 for a voxel size of 0.016 mm. The voltage and current settings of the source were set to 140 kV and 80 μA, respectively. A 0.5 mm Cu prefilter was used. The acquisition time per projection was 6 s.
The geometry of the CT scan of the cylindrical sample is represented, not to scale, by Figure 2. The orange cone represents the field of view (FoV), while the blue dashed lines represent rays that traverse the front and rear edges of the sample. Near the lower edge of the FoV, an inverse conical area of the sample is defined by the solid and dashed orange lines. This is the part of the sample that lies within the FoV during only some of the projections and, therefore, is not fully reconstructible by FDK [13].
We consider the case that dimensions of the sample are not known precisely. It is possible to determine them from the projections. The total height h of the sample within the FoV is the sum of its parts h1 and h2 that extend respectively above and below the plane SOD. The plane is defined by the source (S), the centre of rotation (O) and the centre of the detector (D). The plane SOD is assumed to be perpendicular to the detector.
We can see from Figure 2 that,
SO r SD = h 1 DF
and,
SO + r SD = h 1 DR
where r is the radius of the sample, and DF and DR the distance between the central detector row and the detector rows where the front (F) and rear (R) edge of the concrete rod are respectively projected.
Dividing Equation (1) by Equation (2) and rearranging, we obtain:
r = DF / DR 1 DF / DR + 1 SO .
The length h1 can be calculated now from either Equation (1) or Equation (2), while the length h2 is:
h 2 = SO + r SD H 2
where H is the height of the detector. The part of the sample that, as mentioned above, does not always lie within the FoV has been accounted for through the inclusion of the radius r in Equation (4).

2.2. The Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) Algorithm

The DIRECTT algorithm was first proposed for the reconstruction of two-dimensional (2D) images by Lange et al. and, in a previous article [12], we introduced a new, more efficient, and fully 3D version. The algorithm operates on finding the best solution possible by mimicking the actual physical projection process, instead of directly solving the inverse problem. It only reconstructs certain voxels during each iteration, simulating the projection of the partial reconstruction, and repeating the workflow for the difference between measured and simulated projections until this difference is sufficiently close to zero [10,11,12]. Although the concrete rod of Figure 1 was one of the two datasets that were used to showcase the performance of DIRECTT, the algorithm was implemented only on the slice that corresponds to the cross section of the sample with the plane SOD [12]. That slice will be hereafter referred to as the central slice. Attempting to implement DIRECTT on the whole dataset does not lead to an improvement over the FDK reconstruction of Figure 1. On the contrary, it results in severe artefacts and missing data because the algorithm fails to predict the decreasing grey values along the z-axis. However, these artefacts can be reduced by modelling them based on the shape of the sample.

2.3. Software

For this work, the forward- and back-projection operations involved in DIRECTT were performed using the Python programming language and the open-source ASTRA (All Scale Tomographic Reconstruction Antwerp) toolbox [14]. Via ASTRA, computationally demanding operations are offloaded to a graphics processing unit using the CUDA (Compute Unified Device Architecture) language. The toolbox also includes several reconstruction algorithms, such as the FDK, the simultaneous iterative reconstruction technique (SIRT) [15], and a conjugate gradient (CG) method based on the Krylov subspace [16], that can run with little input from the user [14].

3. Results

A single measured projection of the concrete and the corresponding simulated projection of a virtual homogeneous cylinder of height h and radius r are shown in Figure 3a,b, respectively. The two horizontal lines in the former indicate the detector rows that correspond to points F and R of Figure 2. The two rows are identified automatically from the projections. Specifically, F is the lowermost detector row containing exclusively values that correspond to the background. Similarly, R is the lowermost row containing exclusively values lower than the mode of all absorption projections, which roughly corresponds to the absorption of rays that penetrate the sample perpendicularly.
By simulating the geometry of the scan, projecting and back-projecting the virtual cylinder, and normalizing for the central slice, a 3D model M of the artefacts that arise during the reconstruction is computed. Normalizing for the central slice ensures that hardly any artefacts are considered for the parts of the volume that always lie within the FoV and can be accurately reconstructed by FDK. Note that the model M needs to be computed just once. The slices of M that correspond to those in Figure 1 are shown in the top row of Figure 4. The volume resulting from the back-projection of the unfiltered projections of the concrete rod, the first step during each iteration of DIRECTT [12], is shown in the bottom row of Figure 4. Both volumes of Figure 4 comprise more slices along the z-axis than that of Figure 1. The extra slices have been included because it is essential for the implementation of any iterative reconstruction algorithm, such as DIRECTT, that the parts of the sample that do not lie within the FoV for every projection are reconstructed too.
There is an obvious qualitative relation between the two rows of Figure 4. The pixel-by-pixel division of the reconstructed volume by the model M results in a “corrected” volume, on which DIRECTT can be successfully implemented. The threshold values, based on which the voxels to be reconstructed during each iteration are selected, are calculated from the central slice as described in [12] but are applied simultaneously on the whole volume. The algorithm terminates when the projections of the reconstruction array match the measured ones. The final reconstructed volume is shown in Figure 5.

4. Discussion and Conclusions

The reconstruction by DIRECTT is an improvement over that by FDK in Figure 1 as there are no artefacts linked to the increasing cone angle. Lacking the ground truth for the reconstructed volume, the evaluation of the two algorithms using a full-reference metric is meaningless. Nevertheless, a quantitative evaluation of the algorithms can be undertaken by calculating the histogram entropy of the respective volumes. The histogram entropy (HE) is a global metric defined according to the relation:
HE = p ( μ ) log [ p ( μ ) ] d μ ,
where p(μ) is the distribution function of the grey values. The value of the HE increases if homogeneously distributed noise is present in the image, and decreases if sharp edges are present [17]. Therefore, a low value is an indication of a good balance between noise and blur. The FDK- and DIRECTT-reconstructed volumes (Figure 1 and Figure 5, respectively) have histogram entropies of 2.86 and 1.89, respectively.
An alternative way to evaluate the results of the algorithms is to simulate the projection of the reconstructed volumes and calculate the Pearson correlation coefficient [12,18] between these projections and the measured ones. The Pearson correlation coefficient (PCC), the value of which can range between 1 for total linear correlation and −1 for total linear anti-correlation, is calculated according to the relation:
PCC M , P = σ M , P σ M σ P ,
where σM,P is the covariance and σM and σP are the standard deviation of the measured and simulated projections, respectively. Such coefficients are calculated for each detector pixel and are plotted in Figure 6. While the PCC values decrease for large cone angles in the case of FDK, they remain near 1 regardless of the angle in the case of DIRECTT.
The comparison between the two algorithms can also be done locally on parts of the volume that are affected in a greater degree by the cone-beam artefacts. For instance, the profiles along the x-axis through the centres of the xy-slices in both Figure 1 and Figure 5 are plotted in Figure 7. It is evident that, in the case of FDK, the attenuation values away from the centre of the sample are underestimated.
Finally, the reconstruction of the volume has also been performed using SIRT and CG, which were programmed to perform a fixed number of 600 iterations. The results are shown in Figure 8. The values of the metrics for each reconstruction algorithm are listed in Table 1. While the simulated projections in the case of both SIRT and CG have a high correlation to the measured ones, the respective histogram entropy values are significantly higher than that of DIRECTT. Moreover, although CG has, for the better part, suppressed the artefacts that are linked to the increasing cone angle, it has, at the same time, resulted in ring-like artefacts exactly where the regions that lie fully within the FoV border the regions that are only partially within it.
To sum up, we have described how the DIRECTT algorithm is adjusted for the reconstruction of cone-beam CT data. The artefacts that normally arise during such an operation have been suppressed resulting in a reconstruction that is a clear improvement over that computed by FDK. The performance of DIRECTT has been quantified by both global and local metrics and compared to the performance of other iterative reconstruction algorithms.

Author Contributions

Conceptualization, S.M. and A.K.; methodology, S.M. and A.K.; software, S.M.; validation, S.M., A.K. and G.B.; formal analysis, S.M.; investigation, S.M.; resources, G.B.; data curation, S.M.; writing—original draft preparation, S.M.; writing—review and editing, S.M., A.K. and G.B.; visualization, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was part of the MUMMERING Innovative Training Network. The project has received funding from the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie Grant Agreement No. 765604.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from BAM and are available from the authors with the permission of BAM.

Acknowledgments

The authors kindly acknowledge Dietmar Meinel (BAM) for providing the experimental data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Chiffre, L.; Carmignato, S.; Kruth, J.-P.; Schmitt, R.; Weckenmann, A. Industrial applications of computed tomography. CIRP Ann. 2014, 63, 655–677. [Google Scholar] [CrossRef]
  2. Bartscher, M.; Hilpert, U.; Fiedler, D. Ermittlung der Messunsicherheit von Computertomographie-Messungen am Beispiel eines Zylinderkopfs (Determination of the Measurement Uncertainty of Computed Tomography Measurements Using a Cylinder Head as an Example). TM Tech. Mess. 2008, 75, 178–186. [Google Scholar] [CrossRef]
  3. Carmignato, S. Accuracy of industrial computed tomography measurements: Experimental results from an international comparison. CIRP Ann. 2012, 61, 491–494. [Google Scholar] [CrossRef]
  4. Brunke, O.; Lübbehüsen, J.; Hansen, F.; Butz, F.F. A new Concept for High-Speed atline and inlineCT for up to 100% Mass Production Process Control. In International Congress of Metrology; EDP Sciences: Les Ulis, France, 2013. [Google Scholar]
  5. Tuy, H.K. An Inversion Formula for Cone-Beam Reconstruction. SIAM J. Appl. Math. 1983, 43, 546–552. [Google Scholar] [CrossRef]
  6. Feldkamp, L.A.; Davis, L.C.; Kress, J.W. Practical cone-beam algorithm. J. Opt. Soc. Am. A 1984, 1, 612–619. [Google Scholar] [CrossRef] [Green Version]
  7. Hsieh, J. Two-pass algorithm for cone-beam reconstruction. In Medical Imaging 2000: Image Processing; International Society for Optics and Photonics: Bellingham, WA, USA, 2000. [Google Scholar]
  8. Han, C.; Baek, J. Multi-pass approach to reduce cone-beam artifacts in a circular orbit cone-beam CT system. Opt. Express 2019, 27, 10108–10126. [Google Scholar] [CrossRef] [PubMed]
  9. Maaß, C.; Dennerlein, F.; Noo, F.; Kachelrieß, M. Comparing short scan CT reconstruction algorithms regarding cone-beam artifact performance. In Proceedings of the IEEE Nuclear Science Symposuim & Medical Imaging Conference, Knoxville, TN, USA, 30 October–6 November 2010. [Google Scholar]
  10. Lange, A.; Hentschel, M.P.; Kupsch, A. Computertomographische Rekonstruktion mit DIRECTT. Mater. Test. 2008, 50, 272–277. [Google Scholar] [CrossRef]
  11. Kupsch, A.; Lange, A.; Hentschel, M.P.; Lück, S.; Schmidt, V.; Grothausmann, R.; Hilger, A.; Manke, I. Missing wedge computed tomography by iterative algorithm DIRECTT. J. Microsc. 2016, 261, 36–45. [Google Scholar] [CrossRef] [PubMed]
  12. Magkos, S.; Kupsch, A.; Bruno, G. Direct Iterative Reconstruction of Computed Tomography Trajectories Reconstruction from limited number of projections with DIRECTT. Rev. Sci. Instrum. 2020, 91, 103107. [Google Scholar] [CrossRef] [PubMed]
  13. Grass, M.; Köhler, T.; Proksa, R. 3D cone-beam CT reconstruction for circular trajectories. Phys. Med. Biol. 1999, 45, 329–347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Van Aarle, W.; Palenstijn, W.J.; Cant, J.; Janssens, E.; Bleichrodt, F.; Dabravolski, A.; Sijbers, J. Fast and flexible X-ray tomography using the ASTRA toolbox. Opt. Express 2016, 24, 25129–25147. [Google Scholar] [CrossRef] [PubMed]
  15. Gregor, J.; Benson, T. Computational Analysis and Improvement of SIRT. IEEE Trans. Med. Imaging 2008, 27, 918–924. [Google Scholar] [CrossRef] [PubMed]
  16. Hestenes, M.R.; Stiefel, E. Methods of conjugate gradients for solving linear systems. J. Res. Natl. Bur. Stand. 1952, 49, 409–436. [Google Scholar] [CrossRef]
  17. Allner, S.; Gustschin, A.; Fehringer, A.; Noël, P.B.; Pfeiffer, F. Metric-guided regularisation parameter selection for statistical iterative reconstruction in computed tomography. Sci. Rep. 2019, 9, 6016. [Google Scholar] [CrossRef] [PubMed]
  18. Rodgers, J.L.; Nicewander, W.A. Thirteen Ways to Look at the Correlation Coefficient. Am. Stat. 1988, 42, 59–66. [Google Scholar] [CrossRef]
Figure 1. Orthogonal slices through the volume of a concrete rod as reconstructed by Feldkamp, Davis and Kress (FDK). The orange lines indicate the relative position of the cross sections.
Figure 1. Orthogonal slices through the volume of a concrete rod as reconstructed by Feldkamp, Davis and Kress (FDK). The orange lines indicate the relative position of the cross sections.
Jimaging 07 00147 g001
Figure 2. Geometric representation of the computed tomography (CT) scan (symmetric with respect to the central plane SOD).
Figure 2. Geometric representation of the computed tomography (CT) scan (symmetric with respect to the central plane SOD).
Jimaging 07 00147 g002
Figure 3. (a) Projection of the concrete rod. The detector rows on which the front (F) and rear (R) edge of the concrete rod are respectively projected are marked by horizontal lines. (b) Simulated projection of a homogeneous cylinder of dimensions equal to the rod. The greyscale values in both figures correspond to the attenuation integral values μd.
Figure 3. (a) Projection of the concrete rod. The detector rows on which the front (F) and rear (R) edge of the concrete rod are respectively projected are marked by horizontal lines. (b) Simulated projection of a homogeneous cylinder of dimensions equal to the rod. The greyscale values in both figures correspond to the attenuation integral values μd.
Jimaging 07 00147 g003
Figure 4. Top row: Orthogonal slices through the model M of the artefacts that arise during the reconstruction; Bottom row: Orthogonal slices through the volume resulting from the back-projection of the unfiltered projections of the concrete rod. The extent of the volume of Figure 1 along the z-axis is indicated by the red dashed lines. Apart from that, the position of the slices shown in either row corresponds precisely to that of the slices in Figure 1.
Figure 4. Top row: Orthogonal slices through the model M of the artefacts that arise during the reconstruction; Bottom row: Orthogonal slices through the volume resulting from the back-projection of the unfiltered projections of the concrete rod. The extent of the volume of Figure 1 along the z-axis is indicated by the red dashed lines. Apart from that, the position of the slices shown in either row corresponds precisely to that of the slices in Figure 1.
Jimaging 07 00147 g004
Figure 5. Final reconstruction of the concrete rod as computed by DIRECTT. The volume has the same size as the volume shown in Figure 1 through omission of its upper and lower slices. The position of the slices shown corresponds precisely to that of the slices in Figure 1.
Figure 5. Final reconstruction of the concrete rod as computed by DIRECTT. The volume has the same size as the volume shown in Figure 1 through omission of its upper and lower slices. The position of the slices shown corresponds precisely to that of the slices in Figure 1.
Jimaging 07 00147 g005
Figure 6. Pearson correlation coefficient between measured projections and simulated projections of the volumes reconstructed by: (a) FDK, and (b) Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT). The Pearson correlation coefficient (PCC) has been calculated for each detector pixel.
Figure 6. Pearson correlation coefficient between measured projections and simulated projections of the volumes reconstructed by: (a) FDK, and (b) Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT). The Pearson correlation coefficient (PCC) has been calculated for each detector pixel.
Jimaging 07 00147 g006
Figure 7. Profile along the x-axis of the attenuation coefficient values of the first slice shown in both Figure 1 and Figure 5.
Figure 7. Profile along the x-axis of the attenuation coefficient values of the first slice shown in both Figure 1 and Figure 5.
Jimaging 07 00147 g007
Figure 8. Orthogonal slices through the volume of a concrete rod as reconstructed by simultaneous iterative reconstruction technique (SIRT, top row) and conjugate gradient (CG, bottom row). The volume has the same size as the volume shown in Figure 1 through omission of its upper and lower slices. The position of the slices shown corresponds precisely to that of the slices in Figure 1.
Figure 8. Orthogonal slices through the volume of a concrete rod as reconstructed by simultaneous iterative reconstruction technique (SIRT, top row) and conjugate gradient (CG, bottom row). The volume has the same size as the volume shown in Figure 1 through omission of its upper and lower slices. The position of the slices shown corresponds precisely to that of the slices in Figure 1.
Jimaging 07 00147 g008
Table 1. Comparison of the performance of the reconstruction algorithms.
Table 1. Comparison of the performance of the reconstruction algorithms.
Reconstruction AlgorithmNumber of IterationsAverage Time per Iteration (s) 1Histogram EntropyMean Value of PCC
FDK154.85 ± 0.262.860.75
DIRECTT56214.51 ± 0.181.890.92
SIRT6006.15 ± 0.052.610.97
CG6006.41 ± 0.012.330.97
1 On a computer equipped with an NVIDIA GeForce GTX 1080 Ti GPU.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Magkos, S.; Kupsch, A.; Bruno, G. Suppression of Cone-Beam Artefacts with Direct Iterative Reconstruction Computed Tomography Trajectories (DIRECTT). J. Imaging 2021, 7, 147. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7080147

AMA Style

Magkos S, Kupsch A, Bruno G. Suppression of Cone-Beam Artefacts with Direct Iterative Reconstruction Computed Tomography Trajectories (DIRECTT). Journal of Imaging. 2021; 7(8):147. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7080147

Chicago/Turabian Style

Magkos, Sotirios, Andreas Kupsch, and Giovanni Bruno. 2021. "Suppression of Cone-Beam Artefacts with Direct Iterative Reconstruction Computed Tomography Trajectories (DIRECTT)" Journal of Imaging 7, no. 8: 147. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7080147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop