Abstract
This paper introduces a new decomposition of the 3D X-ray transform based on the shearlet representation, a multiscale directional representation which is optimally efficient in handling 3D data containing edge singularities. Using this decomposition, we derive a highly effective reconstruction algorithm yielding a near-optimal rate of convergence in estimating piecewise smooth objects from 3D X-ray tomographic data which are corrupted by white Gaussian noise. This algorithm is achieved by applying a thresholding scheme on the 3D shearlet transform coefficients of the noisy data which, for a given noise level ε, can be tuned so that the estimator attains the essentially optimal mean square error rate O(log(ε − 1)ε 2/3), as ε→0. This is the first published result to achieve this type of error estimate, outperforming methods based on Wavelet-Vaguelettes decomposition and on SVD, which can only achieve MSE rates of O(ε 1/2) and O(ε 1/3), respectively.
Similar content being viewed by others
References
Candès, E.J., Donoho, D.L.: Recovering edges in ill-posed inverse problems: optimality of curvelet frames. Ann. Stat. 30(3), 784–842 (2002)
Candès, E.J., Demanet, L., Donoho, D.L., Ying, L.: Fast discrete curvelet transforms. SIAM Multiscale Model. Simul. 5, 861–899 (2006)
Colonna, F., Easley, G.R., Guo, K., Labate, D.: Radon transform inversion using the shearlet representation. Appl. Comput. Harmon. Anal. 29(2), 232–250 (2010)
Cormack, A.M.: Representation of a function by its line integrals, with some radiological applications. J. Appl. Phys. 34, 2722–2727 (1963)
Cormack, A.M.: Representation of a function by its line integrals, with some radiological applications II. J. Appl. Phys. 35, 195–207 (1964)
Dahlke, S., Steidl, G., Teschke, G.: Shearlet coorbit spaces: compactly supported analyzing shearlets, traces and embeddings. J. Fourier Anal. Appl. 17(6), 1232–1255 (2011)
Dahlke, S., Häuser, S., Steidl, G., Teschke, G.: Shearlet coorbit spaces: traces and embeddings in higher dimensions—extended version. Monatshefte Math. (2012, to appear)
Donoho, D.L.: Nonlinear solution of linear inverse problem by Wavelet–Vaguelette Decomposition. Appl. Comput. Harmon. Anal. 2, 101–126 (1995)
Donoho, D.L., Johnstone, I.M.: Ideal spatial adaptation via wavelet shrinkage. Biometrika 81, 425–455 (1994)
Easley, G.R., Labate, D., Lim, W.: Sparse directional image representations using the discrete shearlet transform. Appl. Comput. Harmon. Anal. 25(1), 25–46 (2008)
Easley, G.R., Colonna, F., Labate, D.: Improved radon based imaging using the shearlet transform. In: Proc. SPIE, Independent Component Analyses, Wavelets, Unsupervised Smart Sensors, and Neural Networks II, vol. 7343. Orlando (2009)
Easley, G.R., Labate, D., Colonna, F.: Shearlet-based total variation diffusion for denoising. IEEE Trans. Image Process. 18(2), 260–268 (2009)
Guo, K., Labate, D.: Characterization and analysis of edges using the continuous shearlet transform. SIAM J. Imag. Sci. 2, 959–986 (2009)
Guo, K., Labate, D.: Optimally sparse 3D approximations using shearlet representations. Electron. Res. Announc. Math. Sci 17, 126–138 (2010)
Guo, K., Labate, D.: The Construction of Smooth Parseval Frames of Shearlets. Preprint (2011)
Guo, K., Labate, D.: Optimally sparse representations of 3D data with C 2 surface singularities using parseval frames of shearlets. SIAM J Math. Anal. 44, 851–886 (2012)
Guo, K., Lim, W., Labate, D., Weiss, G., Wilson, E.: Wavelets with composite dilations. Electron. Res. Announc. Am. Math. Soc. 10, 78–87 (2004)
Guo, K., Lim, W., Labate, D., Weiss, G., Wilson, E.: Wavelets with composite dilations and their MRA properties. Appl. Comput. Harmon. Anal. 20, 231–249 (2006)
Han, B., Kutyniok, G., Shen, Z.: Adaptive multiresolution analysis structures and shearlet systems. SIAM J. Numer. Anal. 49, 1921–1946 (2011)
Kolaczyk, E.D.: A wavelet shrinkage approach to tomographic image reconstruction. J. Am. Stat. Assoc. 91, 1079–1090 (1996)
Lee, N.Y., Lucier, B.J.: Wavelet methods for inverting the radon transform with noisy data. IEEE Trans. Image Process. 10(1), 79–94 (2001)
Natterer, F., Wübbeling, F.: Mathematical methods in image reconstruction. In: SIAM Monographs on Mathematical Modeling and Computation. SIAM, Philadelphia (2001)
Neelamani, R., Nowak, R., Baraniuk, R.: Model-based inverse halftoning with Wavelet–Vaguelette deconvolution. In: Proc. IEEE Int. Conf. Image Processing ICIP 2000, (Vancouver, Canada), pp. 973–976 (2000)
Negi, P.S., Labate, D.: 3D discrete shearlet transform and video processing. IEEE Trans. Image Process. 21(6), 2944–2954 (2012)
Patel, V.M., Easley, G., Healy, D.M.: Shearlet-based deconvolution. IEEE Trans. Image Process. 18(12), 2673–2685 (2009)
Radon, J.: Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten. In: Berichte Sächsische Akademie der Wissenschaften, Leipzig, Math.–Phys. Kl., vol. 69, pp. 262–267 (1917). Reprinted in: Helgason, S.: The Radon Transform, pp. 177–192. Birkhäuser, Boston (1980)
Yi, S., Labate, D., Easley, G.R., Krim, H.: A shearlet approach to edge analysis and detection. IEEE Trans. Image Process. 18(5), 929–941 (2009)
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by: Yang Wang.
Rights and permissions
About this article
Cite this article
Guo, K., Labate, D. Optimal recovery of 3D X-ray tomographic data via shearlet decomposition. Adv Comput Math 39, 227–255 (2013). https://doi.org/10.1007/s10444-012-9276-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10444-012-9276-x