Tomographic reconstruction methods for decomposing directional components

Previous Article
Inverse source problems in electrodynamics
IPI Home
This Issue
Next Article

December 2018, 12(6): 1429-1442. doi: 10.3934/ipi.2018060

Tomographic reconstruction methods for decomposing directional components

Rasmus Dalgas Kongskov ^1, and Yiqiu Dong ^1,2,,

1.	Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
2.	College of Mathematics and Statistics, Shenzhen University, Shenzhen, Guangdong, China

* Corresponding author: Yiqiu Dong

Received February 2018 Revised May 2018 Published October 2018

Fund Project: The work was supported by Advanced Grant 291405 from the European Research Council and Grant 11701388 from the National Natural Science Foundation of China

X-ray computed tomography technique has been used in many different practical applications. Often after reconstruction we need segment or decompose objects into different components. In this paper, we propose two new reconstruction methods that can decompose objects at the same time. By incorporating direction information, the proposed methods can decompose objects into various directional components. Furthermore, we propose an algorithm to obtain the direction information of the object directly from its CT measurements. We demonstrate the proposed methods on simulated and real samples to show their practical applicability. The numerical results show the differences between the two methods and effectiveness as dealing with fibre-crack decomposition problems.

Keywords: X-ray computed tomography, reconstruction, decomposition, directional total variation regularization, variational methods.

Mathematics Subject Classification: Primary: 49N45; Secondary: 65F22.

Citation: Rasmus Dalgas Kongskov, Yiqiu Dong. Tomographic reconstruction methods for decomposing directional components. Inverse Problems & Imaging, 2018, 12 (6) : 1429-1442. doi: 10.3934/ipi.2018060

References:

[1]	G. Aubert and P. Kornprobst, Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations, Springer Science + Business Media, 2006.
[2]	J. F. Aujol, G. Aubert, L. Blanc-Féraud and A. Chambolle, Image decomposition into a bounded variation component and an oscillating component, J. Math. Imaging Vis., 22 (2005), 71-88. doi: 10.1007/s10851-005-4783-8.
[3]	J. F. Aujol, G. Gilboa, T. Chan and S. Osher, Structure-texture image decomposition-modeling, algorithms, and parameter selection, Int. J. Comput. Vis., 67 (2006), 111-136.
[4]	I. Bayram and M. E. Kamasak, A directional total variation, Eur. Signal Process. Conf., 19 (2012), 265-269.
[5]	A. Bovik, Handbook of Image and Video Processing, Academic Press, 2000.
[6]	T. M. Buzug, Computed Tomography : From Photon Statistics to Modern Cone-Beam CT, Springer, 2008.
[7]	A. Chambolle and P.-L. Lions, Image recovery via total variation minimization and related problems, Numer. Math., 76 (1997), 167-188. doi: 10.1007/s002110050258.
[8]	A. Chambolle and T. Pock, A first-order primal-dual algorithm for convex problems with applications to imaging, J. Math. Imaging Vis., 40 (2011), 120-145. doi: 10.1007/s10851-010-0251-1.
[9]	A. H. Delaney and Y. Bresler, Globally convergent edge-preserving regularized reconstruction: An application to limited-angle tomography, IEEE Trans. Image Process., 7 (1998), 204-221.
[10]	G. R. Easley, D. Labate and F. Colonna, Shearlet-based total variation diffusion for denoising, IEEE Trans. Image Process., 18 (2009), 260-268. doi: 10.1109/TIP.2008.2008070.
[11]	S. Esedoglu and S. J. Osher, Decomposition of images by the anisotropic Rudin-Osher-Fatemi model, Commun. Pure Appl. Math., 57 (2004), 1609-1626. doi: 10.1002/cpa.20045.
[12]	J. Gilles, Noisy image decomposition: A new structure, texture and noise model based on local adaptivity, J. Math. Imaging Vis., 28 (2007), 285-295. doi: 10.1007/s10851-007-0020-y.
[13]	M. Holler and K. Kunisch, On infimal convolution of total variation type functionals and applications, SIAM J. Imaging Sci., 7 (2014), 2258-2300. doi: 10.1137/130948793.
[14]	K. M. Jespersen, J. Zangenberg, T. Lowe, P. J. Withers and L. P. Mikkelsen, Fatigue damage assessment of uni-directional non-crimp fabric reinforced polyester composite using X-ray computed tomography, Compos. Sci. Technol., 136 (2016), 94-103.
[15]	R. D. Kongskov and Y. Dong, Directional total generalized variation regularization for impulse noise removal, Proceedings of Scale Sp. Var. Methods Comput. Vis. 2017, LNCS 10302 (2017), 221–231, URL http://link.springer.com/10.1007/978-3-642-24785-9.
[16]	R. D. Kongskov, Y. Dong and K. Knudsen, Directional total generalized variation regularization, arXiv: 1701.02675
[17]	V. P. Krishnan and E. T. Quinto, Handbook of Mathematical Methods in Imaging, Springer Science + Business Media, 2015.
[18]	J. Li, C. Miao, Z. Shen, G. Wang and H. Yu, Robust frame based X-ray CT reconstruction, J. Comput. Math., 34 (2016), 683-704. doi: 10.4208/jcm.1608-m2016-0499.
[19]	Y. Lou, T. Zeng, S. Osher and J. Xin, A weighted difference of anisotropic and isotropic total variation model for image processing, SIAM J. Imaging Sci., 8 (2015), 1798-1823. doi: 10.1137/14098435X.
[20]	Y. Meyer, Oscillating Patterns in Image Processing and Nonlinear Evolution Equations, vol. 22, American Mathematical Society, 2001. doi: 10.1090/ulect/022.
[21]	W. J. Palenstijn, K. J. Batenburg and J. Sijbers, Performance improvements for iterative electron tomography reconstruction using graphics processing units (GPUs), J. Struct. Biol., 176 (2011), 250-253.
[22]	E. T. Quinto, Singularities of the X-ray transform and limited data tomography in $ \mathbb{R}^2$ and $ \mathbb{R}^3$, SIAM J. Math. Anal., 24 (1993), 1215-1225. doi: 10.1137/0524069.
[23]	J. Radon, Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten, Akad. Wiss., 69 (1917), 262-277.
[24]	J. E. Rouse, Characterisation of Impact Damage in Carbon Fibre Reinforced Plastics by 3D X-Ray Tomography, PhD thesis, University of Manchester, 2012.
[25]	L. I. Rudin, S. Osher and E. Fatemi, Nonlinear total variation based noise removal algorithms, Phys. D Nonlinear Phenom., 60 (1992), 259-268. doi: 10.1016/0167-2789(92)90242-F.
[26]	S. R. Sandoghchi, G. T. Jasion, N. V. Wheeler, S. Jain, Z. Lian, J. P. Wooler, R. P. Boardman, N. K. Baddela, Y. Chen, J. R. Hayes, E. N. Fokoua, T. Bradley, D. R. Gray, S. M. Mousavi, M. N. Petrovich, F. Poletti and D. J. Richardson, X-ray tomography for structural analysis of microstructured and multimaterial optical fibers and preforms, Opt. Express, 22 (2014), 26181.
[27]	O. Scherzer, Handbook of Mathematical Methods in Imaging, Springer Science + Business Media, 2010.
[28]	E. Y. Sidky and X. Pan, Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization, Phys. Med. Biol., 53 (2008), 4777-4807.
[29]	J. L. Starck, M. Elad and D. L. Donoho, Image decomposition via the combination of sparse representation and a variational approach, IEEE Trans. Image Process., 14 (2005), 1570-1582. doi: 10.1109/TIP.2005.852206.
[30]	W. van Aarle, W. J. Palenstijn, J. Cant, E. Janssens, F. Bleichrodt, A. Dabravolski, J. De Beenhouwer, K. J. Batenburg and J. Sijbers, Fast and flexible X-ray tomography using the ASTRA toolbox, Opt. Express, 24 (2016), 25129-25147.
[31]	W. van Aarle, W. J. Palenstijn, J. De Beenhouwer, T. Altantzis, S. Bals, K. J. Batenburg and J. Sijbers, The ASTRA Toolbox: A platform for advanced algorithm development in electron tomography, Ultramicroscopy, 157 (2015), 35-47.
[32]	Z. Wang, A. C. Bovik, H. R. Sheikh and E. P. Simoncelli, Image quality assessment: from error visibility to structural similarity, IEEE Trans. Image Process., 13 (2004), 600-612.
[33]	J. Weickert, Anisotropic Diffusion in Image Processing, Teubner-Verlag, Stuttgart, Germany, 1998.
[34]	J. Weickert, Coherence-enhancing diffusion filtering, Int. J. Comput. Vis., 31 (1999), 111-127.

show all references

References:

[1]	G. Aubert and P. Kornprobst, Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations, Springer Science + Business Media, 2006.
[2]	J. F. Aujol, G. Aubert, L. Blanc-Féraud and A. Chambolle, Image decomposition into a bounded variation component and an oscillating component, J. Math. Imaging Vis., 22 (2005), 71-88. doi: 10.1007/s10851-005-4783-8.
[3]	J. F. Aujol, G. Gilboa, T. Chan and S. Osher, Structure-texture image decomposition-modeling, algorithms, and parameter selection, Int. J. Comput. Vis., 67 (2006), 111-136.
[4]	I. Bayram and M. E. Kamasak, A directional total variation, Eur. Signal Process. Conf., 19 (2012), 265-269.
[5]	A. Bovik, Handbook of Image and Video Processing, Academic Press, 2000.
[6]	T. M. Buzug, Computed Tomography : From Photon Statistics to Modern Cone-Beam CT, Springer, 2008.
[7]	A. Chambolle and P.-L. Lions, Image recovery via total variation minimization and related problems, Numer. Math., 76 (1997), 167-188. doi: 10.1007/s002110050258.
[8]	A. Chambolle and T. Pock, A first-order primal-dual algorithm for convex problems with applications to imaging, J. Math. Imaging Vis., 40 (2011), 120-145. doi: 10.1007/s10851-010-0251-1.
[9]	A. H. Delaney and Y. Bresler, Globally convergent edge-preserving regularized reconstruction: An application to limited-angle tomography, IEEE Trans. Image Process., 7 (1998), 204-221.
[10]	G. R. Easley, D. Labate and F. Colonna, Shearlet-based total variation diffusion for denoising, IEEE Trans. Image Process., 18 (2009), 260-268. doi: 10.1109/TIP.2008.2008070.
[11]	S. Esedoglu and S. J. Osher, Decomposition of images by the anisotropic Rudin-Osher-Fatemi model, Commun. Pure Appl. Math., 57 (2004), 1609-1626. doi: 10.1002/cpa.20045.
[12]	J. Gilles, Noisy image decomposition: A new structure, texture and noise model based on local adaptivity, J. Math. Imaging Vis., 28 (2007), 285-295. doi: 10.1007/s10851-007-0020-y.
[13]	M. Holler and K. Kunisch, On infimal convolution of total variation type functionals and applications, SIAM J. Imaging Sci., 7 (2014), 2258-2300. doi: 10.1137/130948793.
[14]	K. M. Jespersen, J. Zangenberg, T. Lowe, P. J. Withers and L. P. Mikkelsen, Fatigue damage assessment of uni-directional non-crimp fabric reinforced polyester composite using X-ray computed tomography, Compos. Sci. Technol., 136 (2016), 94-103.
[15]	R. D. Kongskov and Y. Dong, Directional total generalized variation regularization for impulse noise removal, Proceedings of Scale Sp. Var. Methods Comput. Vis. 2017, LNCS 10302 (2017), 221–231, URL http://link.springer.com/10.1007/978-3-642-24785-9.
[16]	R. D. Kongskov, Y. Dong and K. Knudsen, Directional total generalized variation regularization, arXiv: 1701.02675
[17]	V. P. Krishnan and E. T. Quinto, Handbook of Mathematical Methods in Imaging, Springer Science + Business Media, 2015.
[18]	J. Li, C. Miao, Z. Shen, G. Wang and H. Yu, Robust frame based X-ray CT reconstruction, J. Comput. Math., 34 (2016), 683-704. doi: 10.4208/jcm.1608-m2016-0499.
[19]	Y. Lou, T. Zeng, S. Osher and J. Xin, A weighted difference of anisotropic and isotropic total variation model for image processing, SIAM J. Imaging Sci., 8 (2015), 1798-1823. doi: 10.1137/14098435X.
[20]	Y. Meyer, Oscillating Patterns in Image Processing and Nonlinear Evolution Equations, vol. 22, American Mathematical Society, 2001. doi: 10.1090/ulect/022.
[21]	W. J. Palenstijn, K. J. Batenburg and J. Sijbers, Performance improvements for iterative electron tomography reconstruction using graphics processing units (GPUs), J. Struct. Biol., 176 (2011), 250-253.
[22]	E. T. Quinto, Singularities of the X-ray transform and limited data tomography in $ \mathbb{R}^2$ and $ \mathbb{R}^3$, SIAM J. Math. Anal., 24 (1993), 1215-1225. doi: 10.1137/0524069.
[23]	J. Radon, Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten, Akad. Wiss., 69 (1917), 262-277.
[24]	J. E. Rouse, Characterisation of Impact Damage in Carbon Fibre Reinforced Plastics by 3D X-Ray Tomography, PhD thesis, University of Manchester, 2012.
[25]	L. I. Rudin, S. Osher and E. Fatemi, Nonlinear total variation based noise removal algorithms, Phys. D Nonlinear Phenom., 60 (1992), 259-268. doi: 10.1016/0167-2789(92)90242-F.
[26]	S. R. Sandoghchi, G. T. Jasion, N. V. Wheeler, S. Jain, Z. Lian, J. P. Wooler, R. P. Boardman, N. K. Baddela, Y. Chen, J. R. Hayes, E. N. Fokoua, T. Bradley, D. R. Gray, S. M. Mousavi, M. N. Petrovich, F. Poletti and D. J. Richardson, X-ray tomography for structural analysis of microstructured and multimaterial optical fibers and preforms, Opt. Express, 22 (2014), 26181.
[27]	O. Scherzer, Handbook of Mathematical Methods in Imaging, Springer Science + Business Media, 2010.
[28]	E. Y. Sidky and X. Pan, Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization, Phys. Med. Biol., 53 (2008), 4777-4807.
[29]	J. L. Starck, M. Elad and D. L. Donoho, Image decomposition via the combination of sparse representation and a variational approach, IEEE Trans. Image Process., 14 (2005), 1570-1582. doi: 10.1109/TIP.2005.852206.
[30]	W. van Aarle, W. J. Palenstijn, J. Cant, E. Janssens, F. Bleichrodt, A. Dabravolski, J. De Beenhouwer, K. J. Batenburg and J. Sijbers, Fast and flexible X-ray tomography using the ASTRA toolbox, Opt. Express, 24 (2016), 25129-25147.
[31]	W. van Aarle, W. J. Palenstijn, J. De Beenhouwer, T. Altantzis, S. Bals, K. J. Batenburg and J. Sijbers, The ASTRA Toolbox: A platform for advanced algorithm development in electron tomography, Ultramicroscopy, 157 (2015), 35-47.
[32]	Z. Wang, A. C. Bovik, H. R. Sheikh and E. P. Simoncelli, Image quality assessment: from error visibility to structural similarity, IEEE Trans. Image Process., 13 (2004), 600-612.
[33]	J. Weickert, Anisotropic Diffusion in Image Processing, Teubner-Verlag, Stuttgart, Germany, 1998.
[34]	J. Weickert, Coherence-enhancing diffusion filtering, Int. J. Comput. Vis., 31 (1999), 111-127.

Figure 1. Comparison on simulated CT reconstruction problem. Regularization parameters in the $\ell_2$-TV and $\ell_2$-DTV models are tuned to maximize the PSNR values. The parameters in DTV are chosen as $a = 0.15$ and $\theta = 20^\circ$

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 2. Left: the objects with estimated direction from the noise-free sinogram. Right: the noise-free sinogram overlaid with the plot of the sum of the magnitudes

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 3. Left: fibre-crack phantom with fibres along the direction $20^\circ$ and cracks in a circular pattern. Right: simulated noise-free sinogram

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 4. Influence of the parameter $\alpha $ in the decomposition model (4) on the results. The SSIM values of $u+w$ are 0.9478, 0.9526, and 0.9513, respectively

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 5. Influence of the parameter $\beta $ in the decomposition model (4) on the results. Here, we set $\alpha = 0.7$ and $\lambda = 0.0038$. The SSIM values of $u+w$ are 0.9526 and 0.9545, respectively

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 6. Comparison the sinogram splitting method by applying FBP and the variational method, which are introduced in Section 3

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 7. Influence of $K$ in the sinogram splitting method on the results by applying the variational method

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 8. Carbon fibre sample from [24]

Figure Options

Download full-size image

Download as PowerPoint slide

Figure 9. Comparison of the sinogram splitting method with the image decomposition method on a real fibre sample

Figure Options

Download full-size image

Download as PowerPoint slide

Table 1. Direction estimation results for the phantom and the real object shown in Figure 2. Note that the correct main direction for the phantom is 20$^\circ$

$\rho $ (%)	0	1	3	5	10	20	30	40
Phantom	20.1	20.1	20.1	20.1	20.1	20.1	20.1	31.7
Real	81.5	81.7	81.5	80.9	81.7	79.5	-1.1	-34.9

Table Options

Download as excel

$\rho $ (%)	0	1	3	5	10	20	30	40
Phantom	20.1	20.1	20.1	20.1	20.1	20.1	20.1	31.7
Real	81.5	81.7	81.5	80.9	81.7	79.5	-1.1	-34.9

[1]	Zhenhua Zhao, Yining Zhu, Jiansheng Yang, Ming Jiang. Mumford-Shah-TV functional with application in X-ray interior tomography. Inverse Problems & Imaging, 2018, 12 (2) : 331-348. doi: 10.3934/ipi.2018015
[2]	Dan Jane, Gabriel P. Paternain. On the injectivity of the X-ray transform for Anosov thermostats. Discrete & Continuous Dynamical Systems - A, 2009, 24 (2) : 471-487. doi: 10.3934/dcds.2009.24.471
[3]	Raymond H. Chan, Haixia Liang, Suhua Wei, Mila Nikolova, Xue-Cheng Tai. High-order total variation regularization approach for axially symmetric object tomography from a single radiograph. Inverse Problems & Imaging, 2015, 9 (1) : 55-77. doi: 10.3934/ipi.2015.9.55
[4]	François Rouvière. X-ray transform on Damek-Ricci spaces. Inverse Problems & Imaging, 2010, 4 (4) : 713-720. doi: 10.3934/ipi.2010.4.713
[5]	Silvia Allavena, Michele Piana, Federico Benvenuto, Anna Maria Massone. An interpolation/extrapolation approach to X-ray imaging of solar flares. Inverse Problems & Imaging, 2012, 6 (2) : 147-162. doi: 10.3934/ipi.2012.6.147
[6]	Nuutti Hyvönen, Martti Kalke, Matti Lassas, Henri Setälä, Samuli Siltanen. Three-dimensional dental X-ray imaging by combination of panoramic and projection data. Inverse Problems & Imaging, 2010, 4 (2) : 257-271. doi: 10.3934/ipi.2010.4.257
[7]	Arun K. Kulshreshth, Andreas Alpers, Gabor T. Herman, Erik Knudsen, Lajos Rodek, Henning F. Poulsen. A greedy method for reconstructing polycrystals from three-dimensional X-ray diffraction data. Inverse Problems & Imaging, 2009, 3 (1) : 69-85. doi: 10.3934/ipi.2009.3.69
[8]	Jakob S. Jørgensen, Emil Y. Sidky, Per Christian Hansen, Xiaochuan Pan. Empirical average-case relation between undersampling and sparsity in X-ray CT. Inverse Problems & Imaging, 2015, 9 (2) : 431-446. doi: 10.3934/ipi.2015.9.431
[9]	Adriana González, Laurent Jacques, Christophe De Vleeschouwer, Philippe Antoine. Compressive optical deflectometric tomography: A constrained total-variation minimization approach. Inverse Problems & Imaging, 2014, 8 (2) : 421-457. doi: 10.3934/ipi.2014.8.421
[10]	Yuyuan Ouyang, Yunmei Chen, Ying Wu. Total variation and wavelet regularization of orientation distribution functions in diffusion MRI. Inverse Problems & Imaging, 2013, 7 (2) : 565-583. doi: 10.3934/ipi.2013.7.565
[11]	Yunhai Xiao, Junfeng Yang, Xiaoming Yuan. Alternating algorithms for total variation image reconstruction from random projections. Inverse Problems & Imaging, 2012, 6 (3) : 547-563. doi: 10.3934/ipi.2012.6.547
[12]	You-Wei Wen, Raymond Honfu Chan. Using generalized cross validation to select regularization parameter for total variation regularization problems. Inverse Problems & Imaging, 2018, 12 (5) : 1103-1120. doi: 10.3934/ipi.2018046
[13]	Florian Krügel. Some properties of minimizers of a variational problem involving the total variation functional. Communications on Pure & Applied Analysis, 2015, 14 (1) : 341-360. doi: 10.3934/cpaa.2015.14.341
[14]	Tim Kreutzmann, Andreas Rieder. Geometric reconstruction in bioluminescence tomography. Inverse Problems & Imaging, 2014, 8 (1) : 173-197. doi: 10.3934/ipi.2014.8.173
[15]	Xiaoqun Zhang, Tony F. Chan. Wavelet inpainting by nonlocal total variation. Inverse Problems & Imaging, 2010, 4 (1) : 191-210. doi: 10.3934/ipi.2010.4.191
[16]	Shousheng Luo, Tie Zhou. Superiorization of EM algorithm and its application in Single-Photon Emission Computed Tomography(SPECT). Inverse Problems & Imaging, 2014, 8 (1) : 223-246. doi: 10.3934/ipi.2014.8.223
[17]	Rinaldo M. Colombo, Francesca Monti. Solutions with large total variation to nonconservative hyperbolic systems. Communications on Pure & Applied Analysis, 2010, 9 (1) : 47-60. doi: 10.3934/cpaa.2010.9.47
[18]	Yunho Kim, Paul M. Thompson, Luminita A. Vese. HARDI data denoising using vectorial total variation and logarithmic barrier. Inverse Problems & Imaging, 2010, 4 (2) : 273-310. doi: 10.3934/ipi.2010.4.273
[19]	Sören Bartels, Marijo Milicevic. Iterative finite element solution of a constrained total variation regularized model problem. Discrete & Continuous Dynamical Systems - S, 2017, 10 (6) : 1207-1232. doi: 10.3934/dcdss.2017066
[20]	Juan C. Moreno, V. B. Surya Prasath, João C. Neves. Color image processing by vectorial total variation with gradient channels coupling. Inverse Problems & Imaging, 2016, 10 (2) : 461-497. doi: 10.3934/ipi.2016008

2017 Impact Factor: 1.465

Tools

Metrics

Tools

Article outline

Figures and Tables

2017 Impact Factor: 1.465

Tools

Figures and Tables

2017 Impact Factor: 1.465

American Institute of Mathematical Sciences