Compensation of ring artefacts in synchrotron tomographic images

Tomographic images are often superimposed by so called ring artefacts. Ring artefacts are concentric rings in the image s round the center of rotation of the tomographic setup caused e.g. by differen ces in the individual pixel response of the detector. They complicate the p ost rocessing of the data, i.e. the segmentation of individual image infor mation. Hence, for a quantitative analysis of the tomographic images a sign ifica t reduction of these artefacts is essential. In this paper, a simple but e fficient method to eliminate such artefacts during the reconstruction is prop osed. © 2006 Optical Society of America OCIS codes: (100.6950) Tomographic image processing References and links 1. G. T. Herman, Image Reconstruction from Projections The Fundamentals of Computed Tomography, (Academic Press New York, 1980). 2. R. Davis, J. C. Elliott, “X-ray microtomography scanner usi ng time-delay integration for elimination of ring artefacts in the reconstructed image,” Nucl. Instrum. Method s Phys. Res. A394, 157-162, (1997). 3. W. Görner, M. P. Hentschel, B. R. M̈ uller, H. Riesemeier, M. Krumrey, G. Ulm, W. Diete, U. Klein, R. F rahm, “BAMline: the first hard X-ray beamline at BESSY II,” Nucl. Ins trum. Methods Phys. Res. A 467-468, 703-706, (2001). 4. J. Sijbers, A. Postnov, “Reduction of ring artifacts in hi gh resolution micro-CT reconstructions,” Phys. Med. Biol. 49, 247-253, (2004). 5. C. Raven, “Numerical removal of ring artifacts in microtomogr aphy,” Rev. Sci. Instrum. 69, 2978-2980, (1998). 6. C. Antoine, P. Nyg̊ ard, Ø. Gregersen, R. Holmstad, T. Weitkamp, C. Rau, “3D images o f paper obtained by phase-contrast X-ray microtomography: image quality and bina risation,” Nucl. Instrum. Methods Phys. Res. A 490, 392-402, (2002). 7. A. Haibel, C. Scheuerlein, “Synchrotron tomography for th e study of void formation in internal tin Nb3Sn superconductors,” submitted to IEEE Trans. Appl. Supercond. 8. A. Goedeke, “Performance Boundaries on Nb 3Sn Superconductors,” Ph.D. thesis, University of Twente, E nschede, Netherlands, 2005. 9. PyHST, http://ftp.esrf.fr/pub/scisoft/ESRF sw/doc/PyHST/GettingPyHST.html, last visited Aug. 2006. 10. M. Rivers, “Tutorial Introduction to X-ray Computed Micr otomography Data Processing,” http://wwwfp.mcs.anl.gov/xray-cmt/rivers/tutorial.html, Universit y of Chicago, 1998. 11. W. Smith, “The Scientist and Engineer’s Guide to Digital S ignal Processing,” http://www.dspguide.com/, last visited Aug. 2006.


Introduction
In synchrotron tomography, ring artefacts are caused by imperfect detector pixel elements as well as by defects or impurities on the scintillator crystals.They appear on synchrotron tomographic images as a number of dark concentric rings superimposed on the image structures.As the grey levels in the reconstructed images are influenced by these ring artefacts, a quantitative analysis of the measured data is difficult.Post processing such as noise reduction, binarization, or segmentation of image information is significantly complicated by the presence of such artefacts.
A common method to reduce ring artefacts is known as flatfield correction [1].Thereby images of the background without sample before and after the data acquisition and in certain intervals during the data acquisition are measured.
The resulting flatfields include the non-uniform sensitivity of the CCD camera pixels, the non-uniform response of the scintillator screen, as well as the inhomogeneities in the incident X-ray beam.The measurement data are corrected by dividing the normalized radiographic images of the samples and the flatfields (I corr = I measurement /I f lat f ield ).However, ring artefacts are not removed completely by this method when different camera elements have intensity dependent, nonlinear response functions or the incident beam has time dependent non-uniformities.Additionally, the effect of the non-uniform sensitivity of different detector pixel elements can be avoided by moving the sample or the detector system during the acquisition in defined horizontal and vertical steps [2,3].The accurately defined shifts are counted back after the measurement.Thus the characteristic of all detector elements is averaged, which leads to significantly reduced ring artefacts.However, this requires very precise translation motor components, especially for measurements with spatial resolutions of about 1µm or less.This is absolutely nessecary to avoid motion artefacts.
A third way to compensate ring artefacts is the post-processing of the radiographic or tomographic images.They can be reduced by different image processing algorithms, e.g.size and shape filtering after the reconstruction or sinogram processing during the reconstruction [4,5,6].
In this paper, a simple but efficient method is presented to eliminate ring artefacts by sinogram processing during the reconstruction.As an example for this artefact a reconstructed slice of a tomographic image superimposed by strong ring artefacts is shown in Fig. 1.The image represents a cross section of a multi-filament wire consisting of superconducting Nb 3 Sn filaments covered by a copper shield [7,8].The diameter of this wire is about 0.7 mm.The reconstruction algorithm, that we used here, is the Filtered Back Projection (see e.g.[1]).It is implemented in the reconstruction program PyHST [9] which was developed at the ESRF, Grenoble.
The image in Fig. 1 contains strong concentric rings that stretch across the whole image.As described the information of such artefacts is stored in all radiographic projections that are captured during the measurement.One line in a radiographic image represents the projection of a slice of the specimen under a certain angle.Arranged among each other the projections of all angles between 0 and 180 degrees produce the so called sinogram which contains the information to reconstruct one tomographic slice.

Filtering
In the case of ring artefacts a sinogram shows dark, thin, vertical lines (see Fig. 2) which correspond to the ring shaped structures in the reconstructed tomographic image.The vertical lines in the sinogram are more easy to detect than the ring shaped structures in the tomographic slice.Hence, it is obvious to filter out the rings by modifying the grey values of these lines in the sinogram in such a way, that they disappear in their environment.If this is achieved, the rings in the tomographic slice will disappear equally.The use of common filtering methods, like unsharp masking and plain low pass filtering of the sinogram, do not show useable results [6].Besides, the filtering parameters of these methods require a lot of adjustment for each ring and over all slices.Due to the automation of the artefact reduction during the reconstruction these filtering methods are out of the question.At first the sum y(i) of all grey values of each column i of a sinogram is calculated over each projection p n .
The graph of these sums shows local minima (peaks) at the same positions that represent the ring artefacts (see red curve in Fig. 2).The next step is to fit a function with the same shape but without the peaks.To align the grey values at the peak positions with the information of their neighbourhood usually high frequency corrections are carried out [5,10].This is by one-dimensional low-pass filtering in fourier space.Figure 3 shows the result of this correction with comparatively good success.Although the artefacts are widely removed the fourier filter does not delete all rings completely.Moreover, the fourier transformation can cause overshoots at jump discontinuities, known as the Gibbs effect [11], which can lead to new ring artefacts.Instead, we concentrate on developing a compensation algorithm in real space.Here, a more simple but efficient method with only one adjustable variable N is presented.With the help of the Moving Average Filter [11] it is possible to delete the peaks (see blue curve in Fig. 2).The applied filtering method replaces a value by averaging the grey values of its neighbourhood with a span width of 2N + 1.The result from Eq. ( 2), whose shape depends on the span factor N, is the smoothed function y s of y.The larger the value of N the more smooth the resulting function will look like.
To filter out the rings in the reconstructed image it is nessecary to find the best value for N. Choosing a small value for N leads to weak filtering and is therefore closer to an unfiltered reconstruction.Otherwise, a very large value results in a blurring of the whole tomographic image.A compromise should be used for an ideal filtering.A clue for that is the number of pixels of the peak width.With both functions, the sum y(i) and the fitted graph y s (i) (see Fig. 2) we are now able to normalize each projection in the sinogram by the ratio of both.
The values of y and y s differ from each other strongly in the area of the peaks, but not outside of them.Thus, these large grey value deviations are smoothed out whereas the rest of the projection is almost unchanged by this normalization.

Results and Conclusion
Ideally, only the grey values at the peak positions should be changed.The rest of them should be unmodified and the sinogram is freed from the vertical lines corresponding to the ring artefacts.Although, in most cases the filtering will affect the whole image, the effect of smoothing for image areas which are outside the ring artefacts is very weak.Only the grey value distribution around the peaks in the function y(i) is smoothed but not each projection in the sinogram itself.The reconstruction of those modified projections shows that the rings disappeared (Fig. 4), which enables the analysis of the intrinsic sample structures, e.g. the segmentation of the image information, without artefacts.The presented algorithm shows a successful way to compensate ring artefacts from tomographic images.The strength is adjustable by the span width of the Moving Average Filter.As a further advantage, the process could be implemented independently from the reconstruction method and is applicable in automation tasks.We implemented this compensation method in the reconstruction program PyHST as a plugin.It is now possible to filter ring artefacts in a tomographic image automatically during its reconstruction.

Fig. 1 .
Fig. 1. 2D slice of a tomographic image: The cross section of a superconductor strand (binary Powder-in-Tube strand of Nb 3 Sn and copper) superimposed by strong ring artefacts.

Fig. 2 .
Fig. 2. The background of the image shows a sinogram, which consists of all measured projections p n (i) of a slice for angles n between 0 and 180 degrees that are arranged among each other.Here, the sinogram reveals vertical lines that are responsible for ring artefacts, the inset in the top right corner emphasises this.The diagram shows the sum y(i) of all grey values for each pixel column in the sinogram and the fitted graph y s (i) of this function.

Fig. 3 .
Fig. 3. Tomographic slice after the high frequency correction [5, 10] of the sinogram.In comparison to the original image the slice contains less rings.But there are still some ring artefacts left.

Fig. 4 .
Fig. 4. Tomographic slice after using the new sinogram correction algorithm without ring artefacts.Here, the algorithm was used with the span factor N = 20.