1. Introduction
In the field of high-resolution X-ray computed tomography (CT) with a micron and submicron spatial resolution, reconstructed CT data are often affected by severe ring artifacts. They appear as concentric ring-like features superimposed on the imaged scene and are centered on the object’s center of rotation creating either full rings (full scan over 360°) or half rings (half scan over 180°) [
1]. Ring artifacts are mainly caused by imperfect detector pixels, where a perfect pixel’s response should be linearly proportional to the amount of photons incident on the detector. There are many different underlying causes for individual pixels to have imperfect responses. These include defects in the scintillator, the detector itself, and the readout electronics [
2]. Moreover, the detector responses may vary due to numerous time-dependent drifts, such as thermal drifts, and also due to changes in the X-ray spectrum [
2]. No matter the cause, ring artifacts degrade the resulting image quality. Therefore, it is desirable to remove or at least significantly reduce the presence of ring artifacts in CT data.
Ring artifacts reduction methods can be divided into three approaches. The first is based on a flat-field correction of a detector [
3]. The proper flat-field correction should ideally remove all the detector sensitivity variations [
4]. However, ring artifacts may persist after this correction due to the detector pixels intensity dependencies and non-linear response functions, or due to time-dependent non-uniformities of the incident beam [
5]. To overcome these issues, advanced flat-field correction approaches were proposed in several recent works [
2,
4,
6,
7]. However, such sophisticated flat-field methods are not easily applicable in practice, because they require specific CT acquisition scenarios, and precise knowledge of used detection system is needed.
Second approach is the hardware-based ring artifacts reduction method. This method is based on moving the detector system in defined horizontal steps during the CT acquisition so that the object is projected on different regions of the detector during a CT scan [
8]. Through this, the effect of non-uniform detector responses is suppressed. Although the practical functionality of this method was reported in [
9] and [
10], this procedure’s disadvantage is reducing the spatial resolution of the CT data if the detector shifts are not accurately known [
2] or the movement precision is worse than the used detector pixel size. In general, this method is hardly applicable in nano CT systems due to such demanding requirements on the movement precision.
The third approach for the ring artifacts reduction are the image-based processing methods. These methods can be further divided, based on the domain of processed data, to sinogram-based (sinogram pre-processing) and tomogram-based (CT data post-processing) methods [
11]. Sinogram-based methods work directly with the sinogram data, where the ring artifacts appear as straight lines in a vertical direction and are therefore easier to detect and to process. Some of these methods assume the presence of a specific high-frequency component that is directly related to the ring artifacts. Therefore, they aim to filter out the artifacts using low pass filters [
12,
13,
14,
15]. Most of these, however, fail to remove the strong artifacts related to dead detector elements or damaged areas on the scintillator, in which case they create an extra band around the original ring [
10]. To overcome this, other methods first detect the ring artifacts elements and then correct them using various approaches: image inpainting [
1,
16,
17,
18] moving average and weighted moving average filters [
19,
20,
21], sensitivity equalization [
22]. However, even these methods have their limitations. Most of these methods are only suitable for suppressing a certain type of stripe. Moreover, they are generally difficult to use in practice due to many parameters needing to be adjusted when a wrong selection of parameters significantly affects the resulting quality. However, the work of Anas [
16] can be pointed out because it introduced a novel idea for classifying rings based on their statistical properties and for addressing each type separately.
On the contrary, the tomogram post-processing methods work with CT data after the tomographic reconstruction. These methods often use a conversion of the data from Cartesian to the polar coordinate system. After this conversion, the ring artifacts appear as stripes that can be further processed using similar assumptions and strategies as for sinogram-based methods [
23,
24,
25,
26,
27]. One method from these can be pointed out, Liang [
26] proposed a novel ring artifacts reduction approach integrating benefits of an efficient iterative framework together with relative total variations (RTV) algorithm for the texture extraction. However, this method uses a simple mean values analysis to detect and correct the artifacts, which is insufficient in the case of dead detector elements or damaged areas on the scintillator. Moreover, tomogram-based methods are, in general, strongly dependent on the quality of the used tomographic reconstruction because some extra artifacts might be created [
20]. Therefore, a novel class of methods lying between sinogram-based and tomogram-based approaches has been recently developed. The ring artifacts reduction is addressed directly during the reconstruction process using specific forms of regularizations (e.g., [
11,
28,
29]). Such regularizations can, however, be highly computationally demanding, which limits the practical applicability of those methods.
In this article, we present a new ring artifacts reduction procedure that combines several selected ideas from image-based processing methods into one complex sinogram-based method with a goal to overcome all previously mentioned limitations. The ring artifacts are classified into two types based on their cause and actual appearance in the CT data. We prefer to separate the detection and correction schemes for each type of artifact for their effective removal. We propose a two-step iterative correction scheme that deals with all the artifact types in the sinogram domain. Consequently, a significant influence of tomographic reconstruction on the efficiency of artifacts reduction is avoided. The reduction strategy was optimized for each artifact type separately to preserve the spatial resolution and sample’s structural information, which are the most important factors in the field of nano-tomography. Practical functionality of the prosed method was verified on both synthetic data and real CT data. It shows a high efficiency of ring artifacts removal, and a robustness to character of input data and used detection system in context of other tested ring artifacts correction techniques.
4. Discussion
The practical testing demonstrated that the proposed ring artifacts reduction procedure, compared to other methods, can achieve superior results in the following criteria: artifacts detection accuracy, overall performance, robustness to detection system, and the ability to preserve the spatial resolution. First, the method was tested in terms of HRA detection accuracy. It was found that for all the tested data, the proposed HRA detection scheme achieved a precision higher than 95% (see
Table 1), even for the increasing noise level present in the data (see
Figure 7). However, for all data, a certain amount of artifacts’ positions was not detected. Moreover, a direct proportion was found between the number of HRA positions that were not detected and the noise level. However, this amount was found to be negligible in terms of the total number of artifacts’ positions, as both the TPR and DSC metrics scored above 80% even for cases with a severe noise presence.
Although the proposed method did not detect all the HRA positions, the overall quantitative and qualitative results were superior to other tested methods. This was achieved by the proposed two-step correction scheme, when the HRA reduction algorithm and LRA reduction algorithm are working in tandem reducing all the artifacts effectively. A high ring artifacts presence in the case of the synthetic data made their correction very challenging, which was reflected by poor results of the M1, M2 and M3 methods (see
Table 2). As for the M1 method, residual rings persisted after the correction for all the tested data resulting in unusable data for further analysis. However, the worst overall results were achieved with the M2 method, as it failed to reduce all the artifacts and preserve the structural information of the input data. Specifically, in the case of ball phantom, the method failed to distinguish the artifacts and sample structure leading to an almost complete suppression of structural information (see
Figure 9), which is further represented by a negative SSIM value (
Table 2). The M3 method achieved visually acceptable results, but the quantitative evaluation on synthetic data revealed that the method had led to a poor preservation of structural information and intensity range of the input data (see
Figure 10). All these findings were further confirmed by testing on the real CT data. The acquired results corresponded to those from testing on the synthetic data. Compared to all other methods, the proposed method obtained the best results, as all the artifacts were reduced, and the sample structures were fully preserved.
The lastly considered criterion focused on the robustness of the proposed ring artifacts reduction procedure to the used detection system and its effect on the spatial resolution. In this analysis, the proposed method proved itself to preserve the spatial resolution within the accuracy limit of the regular standard for both detection systems (see
Table 3). Moreover, the proposed method was functional regardless of the detection system without the need for any settings’ optimization.
All the beforehand discussed aspects restrict the application of the M1, M2 and M3 methods in nano-tomography, where preservation of quality and structural information of input data are of key importance. Moreover, these methods were found highly dependent on the used parameters and the character of input data. Even following the original authors’ recommendations and optimization, the methods did not achieve acceptable results with one setting for all the tested data. This showed a limited robustness and applicability of these methods in practice. Only the proposed method achieved acceptable results for all the conducted tests and showed a high robustness to the character of input data in terms of structure complexity and also the used detection system.
5. Conclusions
The small size of detector pixels used in nano CT devices does not enable an application of any hardware-based method for removing ring artifacts, leaving image-based processing methods as the most promising way for an effective ring artifacts removal. Several approaches from this class exist but each with some deficiencies, such as the degradation of data quality and spatial resolution, which is inconsistent with the core purpose of nano-tomography. The procedure presented in this paper is based on a smart implementation of several ideas from existing methods and utilization of their advantages.
The ring artifacts are classified into two types based on their cause and actual appearance in the CT data. Each artifact class is then handled separately since it is impossible for a single approach to remove all of them. In our procedure, we applied novel iterative RTV-based algorithms in the sinogram domain to avoid any negative influence of tomographic reconstruction. The proposed procedure was optimized and tested on different types of data, cameras, and samples as well.
In confrontation with other advanced ring artifacts reduction methods, it proved its supremacy during practical tests, being robust regarding the character of input data and used detection system. Moreover, the method was able to fully preserve the input data, structural information and spatial resolution. Such features show a high potential of the proposed procedure for practical use in the field of synchrotron- or lab-based nano CT systems.