Improving spatial-resolution in high cone-angle micro-CT by source deblurring

Micro scale computed tomography (CT) can resolve many features in cellular structures, bone formations, minerals properties and composite materials not seen at lower spatial-resolution. Those features enable us to build a more comprehensive model for the object of interest. CT resolution is limited by a fundamental trade off between source size and signal-to-noise ratio (SNR) for a given acquisition time. There is a limit on the X-ray flux that can be emitted from a certain source size, and fewer photons cause a lower SNR. A large source size creates penumbral blurring in the radiograph, limiting the effective spatial-resolution in the reconstruction. High cone-angle CT improves SNR by increasing the X-ray solid angle that passes through the sample. In the high cone-angle regime current source deblurring methods break down due to incomplete modelling of the physical process. This paper presents high cone-angle source de-blurring models. We implement these models using a novel multi-slice Richardson-Lucy (M-RL) and 3D Conjugate Gradient deconvolution on experimental high cone-angle data to improve the spatial-resolution of the reconstructed volume. In M-RL, we slice the back projection volume into subsets which can be considered to have a relative uniform convolution kernel. We compare these results to those obtained from standard reconstruction techniques and current source deblurring methods (i.e. 2D Richardson-Lucy in the radiograph and the volume respectively).

[1]  P.J. La Riviere,et al.  Correction for Resolution Non-Uniformities Caused by Anode Angulation in Computed Tomography , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[2]  L. Canioni,et al.  Propagation beam consideration for 3D THz computed tomography. , 2012, Optics express.

[3]  Mario Bertero,et al.  A simple method for the reduction of boundary effects in the Richardson-Lucy approach to image deconvolution , 2005 .

[4]  Martin J Blunt,et al.  Pore-network extraction from micro-computerized-tomography images. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[6]  T Varslot,et al.  High-resolution helical cone-beam micro-CT with theoretically-exact reconstruction from experimental data. , 2011, Medical physics.

[7]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[8]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[9]  Aggelos K. Katsaggelos,et al.  A regularized iterative image restoration algorithm , 1991, IEEE Trans. Signal Process..

[10]  L. Lucy An iterative technique for the rectification of observed distributions , 1974 .

[11]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[12]  A. Kingston,et al.  Extending reference scan drift correction to high-magnification high-cone-angle tomography. , 2011, Optics letters.

[13]  T Varslot,et al.  Reliable automatic alignment of tomographic projection data by passive auto-focus. , 2011, Medical physics.

[14]  Patrick J. La Rivière,et al.  Correction for Resolution Nonuniformities Caused by Anode Angulation in Computed Tomography , 2008, IEEE Transactions on Medical Imaging.

[15]  X. Zhang,et al.  Pulsed terahertz tomography , 2004 .

[16]  C. Braithwaite Carbonate Rocks , 1967, Nature.

[17]  David M. Paganin,et al.  Coherent X-Ray Optics , 2006 .