Adaptive weighted total variation regularized phase retrieval in differential phase-contrast imaging

Abstract. Phase retrieval with unidirectional differential phase-contrast image requires integration with noisy data, which is an illposed inverse problem. Conventional direct integration method would result in severe streak artifacts. Total variation (TV) regularization-based method would reduce the streak artifacts, but the edges parallel with phase-contrast sensitivity direction are likely to be over smoothed. We propose an improved weighted TV regularization phase retrieval method by introducing a weighting factor to the conventional TV term. When applied to simulation and experimental data, this method shows an advantage of preserving the sharpness of the edges while preserving the ability of reducing streak artifacts compared with conventional TV-regularization method.

[1]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization. , 2010, Physics in medicine and biology.

[2]  D. Donoho.  For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[3]  Miles N Wernick,et al.  Computation of mass-density images from x-ray refraction-angle images. , 2006, Physics in medicine and biology.

[4]  Peter Modregger,et al.  Non-linear regularized phase retrieval for unidirectional X-ray differential phase contrast radiography. , 2011, Optics express.

[5]  P. Munro,et al.  Noninterferometric phase-contrast images obtained with incoherent x-ray sources. , 2011, Applied optics.

[6]  Richard J. Fitzgerald,et al.  Phase‐Sensitive X‐Ray Imaging , 2000 .

[7]  E. Pisano,et al.  Diffraction enhanced x-ray imaging. , 1997, Physics in medicine and biology.

[8]  Peter M Edic,et al.  A Fourier-domain algorithm for total-variation regularized phase retrieval in differential X-ray phase contrast imaging. , 2014, Optics express.

[9]  Gerhard Martens,et al.  Contrast-to-noise in X-ray differential phase contrast imaging , 2011 .

[10]  S. Wilkins,et al.  Phase-contrast imaging of weakly absorbing materials using hard X-rays , 1995, Nature.

[11]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization , 2010, Physics in medicine and biology.

[12]  Tadashi Hattori,et al.  X-ray phase imaging: from synchrotron to hospital , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[13]  O. Bunk,et al.  A two-directional approach for grating based differential phase contrast imaging using hard x-rays. , 2007, Optics express.

[14]  O. Bunk,et al.  Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources , 2006 .

[15]  Michael Unser,et al.  Improved phase retrieval for X-ray differential phase-contrast radiography , 2014 .

[16]  L D Chapman,et al.  Common characteristics shared by different differential phase contrast imaging methods. , 2014, Applied optics.

[17]  Michael Unser,et al.  Spline based iterative phase retrieval algorithm for X-ray differential phase contrast radiography. , 2015, Optics express.

[18]  Timm Weitkamp,et al.  Two-dimensional x-ray grating interferometer. , 2010, Physical review letters.

[19]  Jorge L Flores,et al.  Phase retrieval from one partial derivative. , 2013, Optics letters.

[20]  Wei Huang,et al.  Computed tomography algorithm based on diffraction-enhanced imaging setup , 2005 .

[21]  Atsushi Momose,et al.  Demonstration of X-Ray Talbot Interferometry , 2003 .