Characterizing, measuring, and utilizing the resolution of CT imagery for improved quantification of fine-scale features

Abstract Quantitative results extracted from computed tomographic (CT) data sets should be the same across resolutions and between different instruments and laboratory groups. Despite the proliferation of scanners and data processing methods and tools, and scientific studies utilizing them, relatively little emphasis has been given to ensuring that these results are comparable or reproducible. This issue is particularly pertinent when the features being imaged and measured are of the same order size as data voxels, as is often the case with fracture apertures, pore throats, and cell walls. We have created a tool that facilitates quantification of the spatial resolution of CT data via its point-spread function (PSF), in which the user draws a traverse across a sharp interface between two materials and a Gaussian PSF is fitted to the blurring across that interface. Geometric corrections account for voxel shape and the angle of the traverse to the interface, which does not need to be orthogonal. We use the tool to investigate a series of grid phantoms scanned at varying conditions and observe how the PSF varies within and between slices. The PSF increases with increasing radial distance within slices, and can increase tangentially with increasing radial distance in CT data sets acquired with relatively few projections. The PSF between CT slices is similar to that within slices when a 2-D detector is used, but is much sharper when the data are acquired one slice at a time with a collimated linear detector array. The capability described here can be used not only to calibrate processing algorithms that use deconvolution operations, but it can also help evaluate scans on a routine basis within and between CT research groups, and with respect to the features within the imagery that are being measured.

[1]  T. Ryan,et al.  Preliminary observations on the calcaneal trabecular microarchitecture of extant large-bodied hominoids. , 2006, American journal of physical anthropology.

[2]  Richard A. Ketcham Accurate Three‐dimensional Measurements of Features in Geological Materials from X‐ray Computed Tomography Data , 2010 .

[3]  Manuel Dierick,et al.  Bronnikov-aided correction for x-ray computed tomography. , 2009, Journal of the Optical Society of America. A, Optics, image science, and vision.

[4]  R. Ketcham,et al.  Acquisition, optimization and interpretation of X-ray computed tomographic imagery: applications to the geosciences , 2001 .

[5]  Jack Bresenham,et al.  Algorithm for computer control of a digital plotter , 1965, IBM Syst. J..

[6]  R. A. Johns,et al.  Nondestructive measurements of fracture aperture in crystalline rock cores using X ray computed tomography , 1993 .

[7]  Richard A. Ketcham,et al.  Computational methods for quantitative analysis of three-dimensional features in geological specimens , 2005 .

[8]  Maarten Van Geet,et al.  Quantitative 3D‐fracture analysis by means of microfocus X‐Ray Computer Tomography (µCT): An example from coal , 2001 .

[9]  P. Suetens,et al.  Metal streak artifacts in X-ray computed tomography: a simulation study , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[10]  P. Cloetens,et al.  Phase objects in synchrotron radiation hard x-ray imaging , 1996 .

[11]  R. Ketcham,et al.  Three-dimensional measurement of fractures in heterogeneous materials using high-resolution X-ray computed tomography , 2010 .

[12]  N. Otsu A threshold selection method from gray level histograms , 1979 .