February  2013, 7(1): 267-282. doi: 10.3934/ipi.2013.7.267

Absorption and phase retrieval with Tikhonov and joint sparsity regularizations

1. 

CREATIS, CNRS UMR 5220, Inserm U630, INSA Lyon, Université Lyon 1, F-69621 Villeurbanne Cedex, France, France

2. 

European Synchrotron Radiation Facility, 6 rue Jules Horowitz, F-38043, Grenoble Cedex, France, France

Received  June 2012 Revised  October 2012 Published  February 2013

The X-ray phase contrast imaging technique relies on the measurement of the Fresnel diffraction intensity patterns associated to a phase shift induced by the object. The simultaneous recovery of the phase and of the absorption is an ill-posed nonlinear inverse problem. In this work, we investigate the resolution of this problem with nonlinear Tikhonov regularization and with a joint sparsity constraint regularization. The regularization functionals are minimized with a Gauss-Newton method and with a fixed point iterative method based on a surrogate functional. The algorithms are evalutated using simulated noisy data. The joint sparsity regularization gives better reconstructions for high noise levels.
Citation: Bruno Sixou, Valentina Davidoiu, Max Langer, Francoise Peyrin. Absorption and phase retrieval with Tikhonov and joint sparsity regularizations. Inverse Problems & Imaging, 2013, 7 (1) : 267-282. doi: 10.3934/ipi.2013.7.267
References:
[1]

S. Bayat, L. Apostol, E. Boller, T. Brochard and F. Peyrin, In vivo imaging of bone micro-architecture in mice with 3D synchrotron radiation micro-tomography,, Nucl. Instrum. Methods. Phys. Res., 548 (2005), 247.   Google Scholar

[2]

M. Born and E. Wolf, "Principles of Optics,", Cambridge University Press, (1997).   Google Scholar

[3]

J. H.Bramble, A. Cohen and W. Dahmen, "Multiscale Problems and Methods in Numerical Simulations,", Lectures given at the C.I.M.E Summer School held in Martina Franca, (2001), 9.   Google Scholar

[4]

E. J. Candès, J. Romberg and T. Tao, Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information,, IEEE Trans. Inf. Theory, 52 (2006), 489.  doi: 10.1109/TIT.2005.862083.  Google Scholar

[5]

C. Chappard, A. Basillais, L. Benhamou, A. Bonassie, N. Bonnet, B. Brunet-Imbault and F. Peyrin, Comparison of synchrotron radiation and conventional X-ray microcomputed tomography for assessing trabecular bone microarchitecture of human femoral heads,, Med. Phys., 33 (2006), 3568.   Google Scholar

[6]

P. Cloetens, R. Barrett, J. Baruchel, J. P. Guigay and M. Schlenker, Phase objects in synchrotron radiation hard X-ray imaging,, J. Phys. D, 29 (1996), 133.   Google Scholar

[7]

I. Daubechies, M. Fornasier and I. Loris, Accelerated projected gradient method for linear inverse problems with sparsity constraints,, J. Fourier Anal. Appl., 14 (2008), 764.  doi: 10.1007/s00041-008-9039-8.  Google Scholar

[8]

V. Davidoiu, B. Sixou, M. Langer and F. Peyrin, Nonlinear iterative phase retrieval based on Fréchet derivative,, Optic Express, 23 (2011), 22809.   Google Scholar

[9]

G. R. Davis and S. L. Wong, X-ray microtomography of bones and teeth,, Physiol. Meas., 17 (1996), 121.   Google Scholar

[10]

M. Defrise, I. Daubechies and C. De Mol, An iterative thresholding algorithm for linear inverse problems with sparsity constraint,, Commun. Pure. Appl. Math., 57 (2004), 1413.  doi: 10.1002/cpa.20042.  Google Scholar

[11]

V. Dicken, A new approach towards simultaneous activity and attenuation reconstruction in emission tomography,, Inverse Problems, 15 (1999), 931.  doi: 10.1088/0266-5611/15/4/307.  Google Scholar

[12]

D. L. Donoho, Compressed sensing,, IEEE Trans. Inf. Theory, 52 (2006), 1289.  doi: 10.1109/TIT.2006.871582.  Google Scholar

[13]

H. Engl, M. Hanke and A. Neubauer, "Regularization of Inverse Problems,", Mathematics and its Applications, 375 (1996).  doi: 10.1007/978-94-009-1740-8.  Google Scholar

[14]

H. Engl, K. Kunisch and A. Neubauer, Convergence rates for Tikhonov regularization of nonlinear ill-posed problems,, Inverse Problems, 5 (1989), 523.   Google Scholar

[15]

M. Fornasier and H. Rauhut, Recovery algorithms for vector-valued data with joint sparsity constraints,, SIAM J. Numer. Anal., 46 (2008), 577.  doi: 10.1137/0606668909.  Google Scholar

[16]

J. P. Guigay, M. Langer, R. Boistel and P. Cloetens, A mixed contrast transfer and transport of intensity approach for phase retrieval in the Fresnel region,, Opt. Lett., 32 (2007), 1617.   Google Scholar

[17]

T. E. Gureyev, Composite techniques for phase retrieval in the Fresnel region,, Opt. Commun., 220 (2003), 49.   Google Scholar

[18]

T. E. Gureyev and K. A. Nugent, Phase retrieval with the transport of intensity equation: Orthogonal series solution for non uniform illumination,, Opt. Commun., 13 (1996), 1670.   Google Scholar

[19]

B. Han and Z. Shen, Dual wavelet frames and Riesz bases in Sobolev spaces,, Constructive Approximation, 29 (2009), 369.  doi: 10.1007/s00365-008-9027-x.  Google Scholar

[20]

M. Langer, P. Cloetens and F. Peyrin, Regularization of phase retrieval with phase attenuation duality prior for 3-D holotomography,, IEEE Trans. Image Process, 19 (2010), 2425.  doi: 10.1109/TIP.2010.2048608.  Google Scholar

[21]

M. Langer, P. Cloetens, J. P. Guigay and F. Peyrin, Quantitative comparison of direct phase retrieval algorithms in in-line phase tomography,, Medical Physics, 35 (2008), 4556.   Google Scholar

[22]

A. Momose, T. Takeda, Y. Tai, A. Yoneyama and K. Hirano, Phase-contrast tomographic imaging using an X-ray interferometer,, J. Synchrotron. Rad., 5 (1998), 309.   Google Scholar

[23]

R. D. Nowak, S. J. Wright and M. A. T. Figueiredo, Sparse reconstruction by separable approximation,, IEEE Trans. Sig. Proc., 57 (2009), 2479.  doi: 10.1109/TSP.2009.2016892.  Google Scholar

[24]

K. A. Nugent, Coherent mehtods in the X-rays science,, Advances in Physics, 59 (2010), 1.   Google Scholar

[25]

S. Nuzzo, F. Peyrin, P. Cloetens, J. Baruchel and G. Boivin, Quantification of the degree of mineralization of bone in three dimensions using synchrotron radiation microtomography,, Med. Phys., 29 (2002), 2672.   Google Scholar

[26]

D. M. Paganin, "Coherent X-Ray Optics,", Oxford University Press, (2006).   Google Scholar

[27]

R. Ramlau, Morozov's discrepancy principle for Tikhonov-regularization of nonlinear operators,, J. Num. Funct. Anal. Opt., 23 (2002), 147.  doi: 10.1081/NFA-120003676.  Google Scholar

[28]

M. Salome, F. Peyrin, P. Cloetens, C. Odet, A. M. Laval-Jeantet, J. Baruchel and P. Spanne, A synchrotron radiation microtomography system for the analysis of trabecular bone samples,, Med. Phys., 26 (1999), 2194.   Google Scholar

[29]

O. Scherzer, M. Grasmair, H. Grossauer, M. Haltmeier and F. Lenzen, "Variational Methods in Imaging,", Applied Mathematical Sciences, 167 (2009).   Google Scholar

[30]

G. Teschke and C. Borries, Accelerated projected steepest descent method for nonlinear inverse problems with sparsity constraints,, Inverse Problems, 26 (2010).  doi: 10.1088/0266-5611/26/2/025007.  Google Scholar

[31]

G. Teschke and R. Ramlau, An iterative algorithm for nonlinear inverse problems with joint sparsity constraints in vector-valued regimes and an application to color image impainting,, Inverse Problems, 23 (2007), 1851.  doi: 10.1088/0266-5611/23/5/005.  Google Scholar

[32]

R. Ramlau and G. Teschke, A Tikhonov-based projection iteration for nonlinear ill-posed problems with sparsity constraints,, Numer. Math., 104 (2006), 177.  doi: 10.1007/s00211-006-0016-3.  Google Scholar

[33]

J. Tropp, Algorithm for simultaneous sparse approximation. Part II: Convex relaxation,, IEEE Transactions on Signal Processing, 86 (2006), 589.   Google Scholar

[34]

T. Weikamp, C. David, O. Bunk, J. Bruder, P. Cloetens and F. Pfeiffer, X-ray phase radiography and tomography of soft tissue using grating interferometry,, Eur. J. Radiol., 68 (2008).   Google Scholar

[35]

S. W. Wilkins, T. E. Gureyev, D. Gao, A. Pogany and A. W. Stevenson, Phase contrast imaging using polychromatic X-rays,, Nature, 384 (1996), 335.   Google Scholar

[36]

E. Zeidler, "Nonlinear Functional Analysis and its Applications. II/B. Nonlinear Monotone Operators,", Springer-Verlag, (1990).   Google Scholar

show all references

References:
[1]

S. Bayat, L. Apostol, E. Boller, T. Brochard and F. Peyrin, In vivo imaging of bone micro-architecture in mice with 3D synchrotron radiation micro-tomography,, Nucl. Instrum. Methods. Phys. Res., 548 (2005), 247.   Google Scholar

[2]

M. Born and E. Wolf, "Principles of Optics,", Cambridge University Press, (1997).   Google Scholar

[3]

J. H.Bramble, A. Cohen and W. Dahmen, "Multiscale Problems and Methods in Numerical Simulations,", Lectures given at the C.I.M.E Summer School held in Martina Franca, (2001), 9.   Google Scholar

[4]

E. J. Candès, J. Romberg and T. Tao, Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information,, IEEE Trans. Inf. Theory, 52 (2006), 489.  doi: 10.1109/TIT.2005.862083.  Google Scholar

[5]

C. Chappard, A. Basillais, L. Benhamou, A. Bonassie, N. Bonnet, B. Brunet-Imbault and F. Peyrin, Comparison of synchrotron radiation and conventional X-ray microcomputed tomography for assessing trabecular bone microarchitecture of human femoral heads,, Med. Phys., 33 (2006), 3568.   Google Scholar

[6]

P. Cloetens, R. Barrett, J. Baruchel, J. P. Guigay and M. Schlenker, Phase objects in synchrotron radiation hard X-ray imaging,, J. Phys. D, 29 (1996), 133.   Google Scholar

[7]

I. Daubechies, M. Fornasier and I. Loris, Accelerated projected gradient method for linear inverse problems with sparsity constraints,, J. Fourier Anal. Appl., 14 (2008), 764.  doi: 10.1007/s00041-008-9039-8.  Google Scholar

[8]

V. Davidoiu, B. Sixou, M. Langer and F. Peyrin, Nonlinear iterative phase retrieval based on Fréchet derivative,, Optic Express, 23 (2011), 22809.   Google Scholar

[9]

G. R. Davis and S. L. Wong, X-ray microtomography of bones and teeth,, Physiol. Meas., 17 (1996), 121.   Google Scholar

[10]

M. Defrise, I. Daubechies and C. De Mol, An iterative thresholding algorithm for linear inverse problems with sparsity constraint,, Commun. Pure. Appl. Math., 57 (2004), 1413.  doi: 10.1002/cpa.20042.  Google Scholar

[11]

V. Dicken, A new approach towards simultaneous activity and attenuation reconstruction in emission tomography,, Inverse Problems, 15 (1999), 931.  doi: 10.1088/0266-5611/15/4/307.  Google Scholar

[12]

D. L. Donoho, Compressed sensing,, IEEE Trans. Inf. Theory, 52 (2006), 1289.  doi: 10.1109/TIT.2006.871582.  Google Scholar

[13]

H. Engl, M. Hanke and A. Neubauer, "Regularization of Inverse Problems,", Mathematics and its Applications, 375 (1996).  doi: 10.1007/978-94-009-1740-8.  Google Scholar

[14]

H. Engl, K. Kunisch and A. Neubauer, Convergence rates for Tikhonov regularization of nonlinear ill-posed problems,, Inverse Problems, 5 (1989), 523.   Google Scholar

[15]

M. Fornasier and H. Rauhut, Recovery algorithms for vector-valued data with joint sparsity constraints,, SIAM J. Numer. Anal., 46 (2008), 577.  doi: 10.1137/0606668909.  Google Scholar

[16]

J. P. Guigay, M. Langer, R. Boistel and P. Cloetens, A mixed contrast transfer and transport of intensity approach for phase retrieval in the Fresnel region,, Opt. Lett., 32 (2007), 1617.   Google Scholar

[17]

T. E. Gureyev, Composite techniques for phase retrieval in the Fresnel region,, Opt. Commun., 220 (2003), 49.   Google Scholar

[18]

T. E. Gureyev and K. A. Nugent, Phase retrieval with the transport of intensity equation: Orthogonal series solution for non uniform illumination,, Opt. Commun., 13 (1996), 1670.   Google Scholar

[19]

B. Han and Z. Shen, Dual wavelet frames and Riesz bases in Sobolev spaces,, Constructive Approximation, 29 (2009), 369.  doi: 10.1007/s00365-008-9027-x.  Google Scholar

[20]

M. Langer, P. Cloetens and F. Peyrin, Regularization of phase retrieval with phase attenuation duality prior for 3-D holotomography,, IEEE Trans. Image Process, 19 (2010), 2425.  doi: 10.1109/TIP.2010.2048608.  Google Scholar

[21]

M. Langer, P. Cloetens, J. P. Guigay and F. Peyrin, Quantitative comparison of direct phase retrieval algorithms in in-line phase tomography,, Medical Physics, 35 (2008), 4556.   Google Scholar

[22]

A. Momose, T. Takeda, Y. Tai, A. Yoneyama and K. Hirano, Phase-contrast tomographic imaging using an X-ray interferometer,, J. Synchrotron. Rad., 5 (1998), 309.   Google Scholar

[23]

R. D. Nowak, S. J. Wright and M. A. T. Figueiredo, Sparse reconstruction by separable approximation,, IEEE Trans. Sig. Proc., 57 (2009), 2479.  doi: 10.1109/TSP.2009.2016892.  Google Scholar

[24]

K. A. Nugent, Coherent mehtods in the X-rays science,, Advances in Physics, 59 (2010), 1.   Google Scholar

[25]

S. Nuzzo, F. Peyrin, P. Cloetens, J. Baruchel and G. Boivin, Quantification of the degree of mineralization of bone in three dimensions using synchrotron radiation microtomography,, Med. Phys., 29 (2002), 2672.   Google Scholar

[26]

D. M. Paganin, "Coherent X-Ray Optics,", Oxford University Press, (2006).   Google Scholar

[27]

R. Ramlau, Morozov's discrepancy principle for Tikhonov-regularization of nonlinear operators,, J. Num. Funct. Anal. Opt., 23 (2002), 147.  doi: 10.1081/NFA-120003676.  Google Scholar

[28]

M. Salome, F. Peyrin, P. Cloetens, C. Odet, A. M. Laval-Jeantet, J. Baruchel and P. Spanne, A synchrotron radiation microtomography system for the analysis of trabecular bone samples,, Med. Phys., 26 (1999), 2194.   Google Scholar

[29]

O. Scherzer, M. Grasmair, H. Grossauer, M. Haltmeier and F. Lenzen, "Variational Methods in Imaging,", Applied Mathematical Sciences, 167 (2009).   Google Scholar

[30]

G. Teschke and C. Borries, Accelerated projected steepest descent method for nonlinear inverse problems with sparsity constraints,, Inverse Problems, 26 (2010).  doi: 10.1088/0266-5611/26/2/025007.  Google Scholar

[31]

G. Teschke and R. Ramlau, An iterative algorithm for nonlinear inverse problems with joint sparsity constraints in vector-valued regimes and an application to color image impainting,, Inverse Problems, 23 (2007), 1851.  doi: 10.1088/0266-5611/23/5/005.  Google Scholar

[32]

R. Ramlau and G. Teschke, A Tikhonov-based projection iteration for nonlinear ill-posed problems with sparsity constraints,, Numer. Math., 104 (2006), 177.  doi: 10.1007/s00211-006-0016-3.  Google Scholar

[33]

J. Tropp, Algorithm for simultaneous sparse approximation. Part II: Convex relaxation,, IEEE Transactions on Signal Processing, 86 (2006), 589.   Google Scholar

[34]

T. Weikamp, C. David, O. Bunk, J. Bruder, P. Cloetens and F. Pfeiffer, X-ray phase radiography and tomography of soft tissue using grating interferometry,, Eur. J. Radiol., 68 (2008).   Google Scholar

[35]

S. W. Wilkins, T. E. Gureyev, D. Gao, A. Pogany and A. W. Stevenson, Phase contrast imaging using polychromatic X-rays,, Nature, 384 (1996), 335.   Google Scholar

[36]

E. Zeidler, "Nonlinear Functional Analysis and its Applications. II/B. Nonlinear Monotone Operators,", Springer-Verlag, (1990).   Google Scholar

[1]

Thorsten Hohage, Mihaela Pricop. Nonlinear Tikhonov regularization in Hilbert scales for inverse boundary value problems with random noise. Inverse Problems & Imaging, 2008, 2 (2) : 271-290. doi: 10.3934/ipi.2008.2.271

[2]

Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems & Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017

[3]

Armin Lechleiter, Marcel Rennoch. Non-linear Tikhonov regularization in Banach spaces for inverse scattering from anisotropic penetrable media. Inverse Problems & Imaging, 2017, 11 (1) : 151-176. doi: 10.3934/ipi.2017008

[4]

Gabriel Peyré, Sébastien Bougleux, Laurent Cohen. Non-local regularization of inverse problems. Inverse Problems & Imaging, 2011, 5 (2) : 511-530. doi: 10.3934/ipi.2011.5.511

[5]

Frank Pörner, Daniel Wachsmuth. Tikhonov regularization of optimal control problems governed by semi-linear partial differential equations. Mathematical Control & Related Fields, 2018, 8 (1) : 315-335. doi: 10.3934/mcrf.2018013

[6]

Daniela Calvetti, Erkki Somersalo. Microlocal sequential regularization in imaging. Inverse Problems & Imaging, 2007, 1 (1) : 1-11. doi: 10.3934/ipi.2007.1.1

[7]

Pei Yean Lee, John B Moore. Gauss-Newton-on-manifold for pose estimation. Journal of Industrial & Management Optimization, 2005, 1 (4) : 565-587. doi: 10.3934/jimo.2005.1.565

[8]

Ke Chen, Yiqiu Dong, Michael Hintermüller. A nonlinear multigrid solver with line Gauss-Seidel-semismooth-Newton smoother for the Fenchel pre-dual in total variation based image restoration. Inverse Problems & Imaging, 2011, 5 (2) : 323-339. doi: 10.3934/ipi.2011.5.323

[9]

Stefan Kindermann, Andreas Neubauer. On the convergence of the quasioptimality criterion for (iterated) Tikhonov regularization. Inverse Problems & Imaging, 2008, 2 (2) : 291-299. doi: 10.3934/ipi.2008.2.291

[10]

Vinicius Albani, Adriano De Cezaro, Jorge P. Zubelli. On the choice of the Tikhonov regularization parameter and the discretization level: A discrepancy-based strategy. Inverse Problems & Imaging, 2016, 10 (1) : 1-25. doi: 10.3934/ipi.2016.10.1

[11]

Jun Lai, Ming Li, Peijun Li, Wei Li. A fast direct imaging method for the inverse obstacle scattering problem with nonlinear point scatterers. Inverse Problems & Imaging, 2018, 12 (3) : 635-665. doi: 10.3934/ipi.2018027

[12]

Bernadette N. Hahn. Dynamic linear inverse problems with moderate movements of the object: Ill-posedness and regularization. Inverse Problems & Imaging, 2015, 9 (2) : 395-413. doi: 10.3934/ipi.2015.9.395

[13]

Bruno Sixou, Cyril Mory. Kullback-Leibler residual and regularization for inverse problems with noisy data and noisy operator. Inverse Problems & Imaging, 2019, 13 (5) : 1113-1137. doi: 10.3934/ipi.2019050

[14]

Ji Li, Tie Zhou. Numerical optimization algorithms for wavefront phase retrieval from multiple measurements. Inverse Problems & Imaging, 2017, 11 (4) : 721-743. doi: 10.3934/ipi.2017034

[15]

Guozhi Dong, Bert Jüttler, Otmar Scherzer, Thomas Takacs. Convergence of Tikhonov regularization for solving ill-posed operator equations with solutions defined on surfaces. Inverse Problems & Imaging, 2017, 11 (2) : 221-246. doi: 10.3934/ipi.2017011

[16]

Vinicius Albani, Adriano De Cezaro. A connection between uniqueness of minimizers in Tikhonov-type regularization and Morozov-like discrepancy principles. Inverse Problems & Imaging, 2019, 13 (1) : 211-229. doi: 10.3934/ipi.2019012

[17]

Luca Rondi. On the regularization of the inverse conductivity problem with discontinuous conductivities. Inverse Problems & Imaging, 2008, 2 (3) : 397-409. doi: 10.3934/ipi.2008.2.397

[18]

Giuseppina di Blasio, Filomena Feo, Maria Rosaria Posteraro. Existence results for nonlinear elliptic equations related to Gauss measure in a limit case. Communications on Pure & Applied Analysis, 2008, 7 (6) : 1497-1506. doi: 10.3934/cpaa.2008.7.1497

[19]

Zhong-Qing Wang, Li-Lian Wang. A Legendre-Gauss collocation method for nonlinear delay differential equations. Discrete & Continuous Dynamical Systems - B, 2010, 13 (3) : 685-708. doi: 10.3934/dcdsb.2010.13.685

[20]

Giovanni Bozza, Massimo Brignone, Matteo Pastorino, Andrea Randazzo, Michele Piana. Imaging of unknown targets inside inhomogeneous backgrounds by means of qualitative inverse scattering. Inverse Problems & Imaging, 2009, 3 (2) : 231-241. doi: 10.3934/ipi.2009.3.231

2018 Impact Factor: 1.469

Metrics

  • PDF downloads (11)
  • HTML views (0)
  • Cited by (6)

[Back to Top]