Next Article in Journal
Containment Control-Guided Boundary Information for Semantic Segmentation
Next Article in Special Issue
A Method for All-Weather Unstructured Road Drivable Area Detection Based on Improved Lite-Mobilenetv2
Previous Article in Journal
Easy Fabrication of Ultrafiltration Membrane via Polyethersulfone-Fumed Silica
Previous Article in Special Issue
Personalized Video Summarization: A Comprehensive Survey of Methods and Datasets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methods for Reducing Ring Artifacts in Tomographic Images Using Wavelet Decomposition and Averaging Techniques

by
Paweł Lipowicz
1,*,
Marta Borowska
1,* and
Agnieszka Dardzińska-Głębocka
2
1
Institute of Biomedical Engineering, Faculty of Mechanical Engineering, Bialystok University of Technology, Wiejska 45C, 15-351 Bialystok, Poland
2
Department of Mechanics and Applied Computer Science, Faculty of Mechanical Engineering, Bialystok University of Technology, Wiejska 45C, 15-351 Bialystok, Poland
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 7292; https://doi.org/10.3390/app14167292
Submission received: 26 June 2024 / Revised: 19 July 2024 / Accepted: 14 August 2024 / Published: 19 August 2024
(This article belongs to the Special Issue Novel Research on Image and Video Processing Technology)

Abstract

:
Computed tomography (CT) is one of the fundamental imaging modalities used in medicine, allowing for the acquisition of accurate cross-sectional images of internal body tissues. However, during the acquisition and reconstruction process, various artifacts can arise, and one of them is ring artifacts. These artifacts result from the inherent limitations of CT scanner components and the properties of the scanned material, such as detector defects, non-uniform distribution of radiation from the source, or the presence of metallic elements within the scanning region. The purpose of this study was to identify and reduce ring artifacts in tomographic images using image decomposition and average filtering methods. In this study, tests were conducted on the effectiveness of identifying ring artifacts using wavelet decomposition methods for images. The test was performed on a Shepp–Logan phantom with implemented artifacts of different intensity levels. The analysis was performed using different wavelet families, and linear approximation methods were used to filter the image in the identified areas. Additional filtering was performed using moving average methods and empirical mode decomposition (EMD) techniques. Image comparison methods, i.e., RMSE, SSIM and MS-SSIM, were used to evaluate performance. The results of this study showed a significant improvement in the quality of tomographic phantom images. The authors obtained more than 50% improvement in image quality with reference to the image without any filtration. The different wavelet families had different efficiencies with relation to the identification of the induction regions of ring artifacts. The Haar wavelet and Coiflet 1 showed the best performance in identifying artifact induction regions, with comparative RMSE values for these wavelets of 0.1477 for Haar and 0.1469 for Coiflet 1. The applied additional moving average filtering and EMD permitted us to improve image quality, which is confirmed by the results of the image comparison. The obtained results allow us to assess how the used methods affect the reduction in ring artifacts in phantom images with induced artifacts.

1. Introduction

Computed tomography is an increasingly popular imaging technique with applications in the fields of medicine, diagnostics, scientific research and industry. A variety of artifacts arise during tomographic imaging. They are caused by various physical phenomena during the emission of radiation. Ring artifacts are one of the most common artifacts. They arise as a result of imperfections in the detector’s pixel elements or fiber optic matrix, or defects on scintillator crystals [1]. Ring artifacts can also result from contamination of any of the previously mentioned elements, as well as various types of drift, for example, thermal drift or drift in the X-ray-gene spectrum [2]. Scintillators are the element located in front of the detector optical fiber matrix and have the function of converting ionizing radiation into visible or near-visible emission. There are a number of materials that exhibit scintillation properties that can exist in crystalline, liquid or gaseous form.
Optical fibers behind the scintillator direct radiation in the visible spectrum emitted by the scintillator directly to individual pixels located on the CCD or CMOS array of the detector (Figure 1). The attenuation or extinguishing of individual pixels during acquisition causes the formation of artifacts in the image [3].
Ring artifacts visible on tomograms appear in the form of circles spreading from the center of the image to its edges. There may be several to dozens of them. The visible circles can range in width from single pixels up to several pixels. They can have different grayscale brightness. They significantly affect the qualitative and quantitative analysis of the measured data. They can also distort important morphometric information in the tomographic image. Artifacts hinder such processes as noise reduction, binarization and segmentation in the image [4,5,6,7].
Methods for removing these artifacts can be divided into three types: pre-acquisition, pre-reconstruction and post-reconstruction methods. The commonly used pre-acquisition method of reducing ring artifacts resulting from unbalanced detector sensitivity is called flat-field correction. Flat-field correction involves pixel correction so that all pixels have an identical gray value when they are exposed to a beam with specific parameters. These parameters are identical to those used during the acquisition of the research object [8]. In this method, measurement of the bright- and dark-field images is performed. The bright-field image is collected during homogeneous unconstrained radiation of the detector with the X-ray beam, and the dark-field image is the image formed while the radiation is disabled. The correction is calculated as follows for a pixel in the (x, y) position: The intensity of the pixel for the bright field is I B ( x , y ) , and for the dark field is I D ( x , y ) , and the average intensity for both fields is μ B and μ D . Based on this information, the specific gain of the pixel G ( x , y ) and the displacement of the pixel O ( x , y ) can be counted according to the following formula:
G ( x , y ) = μ B μ D I B x , y I D ( x , y )
O x , y = μ B G ( x , y ) I B ( x , y )
The resulting data can be applied to correct the raw signal I U ( x , y ) obtained during acquisition, the result being a corrected pixel value I C ( x , y ) :
I C x , y = I u x , y G x , y + O ( x , y )
The presented method significantly improves image acquisition; however, it does not provide complete removal of ring artifacts [9]. This is due to the nonlinearity of the response function of the photosensitive components of the detector, and to the nonlinearity in time of the uniformity of the radiation beam.
Pre-reconstruction and post-reconstruction methods involve post-acquisition pixel correction. They implement image filtering methods involving entire images or specific areas of artifact occurrence. Pre-reconstruction methods involve translating projections into sinograms. The damaged area in the sinogram is represented by a darker strip running across all or part of the sinogram. The simplest methods are based on determining the positions of the stripes based on the sum of the pixel values in the image columns. The positions with a sum under a specific threshold are the location of the stripes which then are processed. Such a proposal for artifact reduction is presented in the work of Mark Rivers [10]. This method is effective for very significant artifacts, for example, in the case of damaged pixels. More advanced methods are based on the processing of the obtained sinograms to visualize striations of variable intensity. Well-developed methods use frequency processing using Fourier transform or wavelet decomposition. Such a method is presented in the work of Munch et al., which uses wavelet decomposition and frequency filtering followed by reconstruction to the original sinogram [11]. This method is more efficient; however, it significantly impacts the pixel values within the entire sinogram. The work of Tang et al. presents a hybrid approach to reducing artifacts in the pre-reconstruction process. It involves using wavelet decomposition to identify the positions of artifact induction and then transforming the image within the raw sinogram at predefined positions [12].
The process of reducing artifacts after reconstruction involves transforming tomograms to polar coordinates, identifying artifacts, and removing them using various image processing techniques. In their paper, Sijbers et al. presented such a method using local statistical evaluation of pixel values to identify stripes in a tomogram [13].
In Table 1, the authors present various papers on ring artifact reduction as well. They show on which data the cited methods were tested and present some relevant results of ring artifact reduction if they were presented in the paper.The authors of the paper present a method for pre-reconstruction artifact reduction that uses the identification of artifact induction sites using wavelet decomposition. They apply sinogram processing in two stages to the identified locations and within the entire sinogram. The authors of the current study have developed an original model for reducing ring artifacts using discrete wavelet decomposition supplemented by filtering using empirical mode decomposition (EMD). The authors have not encountered in the literature the use of wavelet decomposition combined with EMD in the filtering of medical tomographic images.

2. Materials and Methods

In this paper, ring artifact reduction was performed on a Shepp–Logan (Yu-Ye-Wang type) tomographic phantom [16]. The phantom was generated at a 512 × 512 × 512 px resolution. All operations, reconstruction and calculations were performed in the Python 3.8 environment with the Tomographic Iterative GPU-based Reconstruction Toolbox (TIGRE toolbox v.2.6) [17]. The python environment uses the libraries NumPy v. 1.21.5 (https://numpy.org, accessed on 20 April 2024), OpenCV v. 4.6.0 (https://opencv.org, accessed on 20 April 2024), Matplotlib v. 3.6.2 [18], PyWavelets v. 1.4.1 [19] and EMD v. 0.5.5 [20]. The development environment was run on a Windows 11 Pro computer equipped with an i7-12700KF 3.60 GHz processor, an NVIDIA GeForce RTX 3070Ti OC 8 GB (New Taipei, Taiwan) graphics card with a 6144 CUDA core, and 64 GB of DDR4 RAM.
After performing ring artifact reduction, the phantom was reconstructed using the FDK algorithm with Shepp–Logan filtering. The process of performing the reconstruction and comparing the results of different reconstruction algorithms with their description was presented in the authors’ earlier work [21].
Before performing ring artifact reduction, artifacts were introduced into the generated phantom. Translation of the phantom to sinograms was performed, and 18 stripes with a width of one pixel were introduced at an equal distance from the center of the sinogram (Figure 2). Tests were performed using stripes with values ranging from 2 to 20% larger than the initial grayscale pixel value. Subsequently, all sinograms S k ( n , j ) were subjected to discrete two-dimensional first-level wavelet decomposition using wavelets from the “Haar”, “Daubechies”, “Symlets”, “Coiflets”, “Biorthogonal” and “Discrete FIR approximation of Mayer wavelet” families. Further analysis was performed on the image after high-pass filtering horizontally scaled and low-pass filtering vertically scaled to the dimensions of the original sinogram [22,23] (Figure 3a).
For each sinogram after decomposition T k ( n , j ) , the lower b L and upper b H levels of the correct pixel values were determined. They were calculated based on the average value of the normalized image of the sinogram before decomposition w s = M e a n ( S k n n , j ) , the average standard deviation of the sinogram after decomposition m 1 = S D ( T k x , y ) , the average value of the pixels in the image after decomposition m 2 = M e a n T k ( x , y ) , and a window width factor k selected experimentally. The calculation of the levels of pixel values and binarization follows the following formulas:
w 1 = k w s w s 4
w 2 = k ( 0.2 w s ) w s 4
b L = m 1 w 2 m 2
b H = m 1 + w 1 m 2
i f   T k n , j b H   a n d   T k n , j b L ,   T k n , j = 0 ;
i f   T k n , j b L   a n d   T k n , j b H ,   T k n , j = 1 .
As a result of these operations, we obtained a matrix of size for image decomposition T k n , j of zeros and ones B k n , j (Figure 3b). The matrix B k n , j was summed with respect to the columns, and the values were divided by the number of rows in the matrix. We obtained a single row matrix N k j = n   B k n , j n . The column positions in the image j of stripe occurrences were identified when N k j T , where T is the thresholding level determined experimentally for a particular data set (Figure 4). The positions of artifact occurrences were subjected to local filtering using line interpolation in a window of 0.5% of the image width rounded to the nearest whole odd number, but not less than 3 pixels (pascal). The operation result was a corrected image S k n , j .
The use of a certain thresholding level resulted in the fact that the striations with very low intensity below the thresholding level were not identified. This occurred especially in situations where tomography was performed on an element with a very different density. In such a situation, it was proposed to perform additional filtering of the entire sinogram. For this purpose, filtering was performed with moving average [24] methods and using empirical mode decomposition [25,26].
In order to perform the filtering of the moving average, the grayscale values of the pixels of the columns of the corrected sinogram were summed, obtaining the curve y j = n S k n , j . Then, the curve was filtered according to the formula in Equation (10), and the whole sinogram was multiplied by the correction factor (Equation (11)):
y s j = 1 2 N + 1 k = N N y ( j + k ) .
S k n , j = S k ( n , j ) y s ( j ) y ( j ) .
where N is the span factor.
Global filtering with mask shift EMD was also performed. The exact process and description are presented in the work of Deering et al. [26]. As a result of decomposition to the 4th level, the harmonics of the signal were extracted. The filtered curve y s j was the curve with the highest harmonics removed. The correction process then proceeded according to Equation (11).
The images resulting from the presented processes were subjected to further analysis at a depth of 32 bits. The corrected and reconstructed images were compared with the reference image using image comparison methods. The Root-Mean-Square Error (RMSE) (Equation (12)) [27], Structural Similarity Index (SSIM) (Equations (13) and (14) [28]), and Multi-scale Structural Similarity Index (MS-SSIM) (Equation (15)) [29] methods were used for comparison.
R M S E = 1 M × N i = 0 ,       j = 0 M 1 ,       N 1 [ k i , j k i , j ] 2
where k   is the reference image and k is the corrected image.
S S I M x , y = l x , y · c x , y β · s x , y γ
S S I M x , y = ( 2 μ x μ y + C 1 ) ( 2 σ x y + C 2 ) ( μ x 2 μ y 2 + C 1 ) ( σ x 2 σ y 2 + C 2 )
where C 1 and C 2 are constants using two scalar constants K 1 and K 2 and the dynamic range of the image. In this paper, K 1 = 0.01 and K 2 = 0.03 . The SSIM indexing algorithm uses a sliding window method to evaluate the image quality. SSIM takes values from −1 to 1; S S I M x , y = 1   if and only if x = y . The window is moved pixel by pixel in the whole image. In this paper, the window size is 8 × 8 px.
M S S S I M x , y = l M x , y M · j = 1 M c j x , y β j · s j x , y γ j
In this paper K 1 = 0.01 , K 2 = 0.03 and the window size is 11 × 11 px.

3. Results and Discussion

In the first stage, experimental determination of the value of the k -factor was performed. For this purpose, the noisy Shepp–Logan phantom was subjected to wavelet decomposition using a Haar wavelet. A ring artifact reduction process was carried out by changing the value of k in the range from 5 to 35, and the thresholding level was fixed at T = 0.4 . Global filtering of the whole sinogram was not performed. The reconstruction result was compared to the reference image using the RMSE and MSE methods. The results of the comparison are presented in the graph in Figure 5. The results presented in the graph show the effect of the k parameter on the effectiveness of reducing ring artifacts in the Shepp–Logan phantom image. From the graph, it can be seen that the optimal parameter for this phantom is values above 21. A further increase in the value has no effect on reducing the errors between the reference image and the one obtained as a result of processing and reconstruction.
After determining the value of the k -factor, the sensitivity of the presented method of reducing ring artifacts to striations of different intensities was checked. For this purpose, the percentage of introduced noise was increased in the range of 2–20%. Reconstruction was performed and the obtained images were compared to the reference image.
Table 2 shows the results of the image comparison for reconstructed images without ring artifact reduction and reconstructed images with artifact reduction. The filtering process used a first-level Haar wavelet, a constant of k = 21 , a thresholding level of T = 0.4 and a span factor of N = 3 (Figure 6).The last row of the Table 2 shows the result of comparing the image of the generated phantom and the phantom after reconstruction without introduced artifacts. The analysis of the results presented in the table shows that the reduction in ring artifacts is effective and significantly improves image quality. The results of the MSE method suggest that image noise below 5% as a result of filtering no longer improves image quality. This is due to global filtering, in which, at such a low noise level, the changes resulting from filtering the entire image make more difference than the noise by itself. However, the SSIM and MS-SSIM methods suggest that the filtered image is of similar or even slightly better quality than without filtering at noise levels of 2% AND 5%. With high-intensity stripes, the image has two or even three times the similarity to the reference image.The quantitative results are worse than in the other works cited [2,14], however, it should be noted that in this work the authors use a virtually generated phantom rather than a phantom reconstruction image as a reference image. The filtered back projection (FBP) reconstruction process further introduces some imperfections and simplifications. The work of author Selim M. [14] uses reconstruction and iterative filtering which gives much better results. However, it requires many iterations is much time-consuming and has much higher hardware requirements. The method proposed in this work is much easier to implement and has more practical application today.
Wavelet decomposition was also performed using wavelet families other than the Haar wavelet. This operation was aimed at checking which of the wavelet families used for decomposition would most precisely identify the vertical stripes and separate them from the non-sine data in the sinogram of the reconstructed object. A phantom with a noise level of 10% was subjected to a series of decompositions. After decomposition, it was subjected to a ring artifact reduction process and compared with the reference image. The results for each phantom are recorded in Table 2 and in the graph (Figure 6).
When analyzing the results, it should be noted that with the RMSE method, the similarity coefficient tends to 0, so the lower the value, the more similar the obtained image is to the reference image. In the case of the SSIM and MS-SSIM methods, the values tend to 1, that is, the closer the value is to one, the greater the similarity of the images. In the results shown in the Table 3 and Figure 7, you can see four wavelets that obtained significantly higher results than the others: Haar, Coiflet 1, Biorthogonal 2.2 and 3.1 and Daubechies 2. The values of the coefficients are included in Table 3, and the four mentioned wavelets are in bold.
It can be noted that the more complex the wavelet, the lower the similarity to the reference image obtained. The lowest efficiency was obtained by performing artifact reduction with the “Discrete” Meyers wavelet. The use of global filtering of the entire image may also have affected the values obtained. A slight reduction in the value of all pixels relative to the reference image is visible in the results as a lower similarity even though the morphometric information of the image was not lost.

4. Conclusions

The results presented in this paper show a method of reducing ring artifacts that effectively improves the image in which ring artifacts of different intensities occur. Depending on the intensity of the ring artifacts, the results obtained were up to 50% better than the image not subjected to any filtering. The comparison uses a computer-generated cross-sectional image as a reference, rather than a noise-free reconstructed phantom image. This causes the comparison results to be lower, due to the fact that the reconstruction itself introduces other errors and imperfections in the image, not only ring artifacts. The FBP method is not perfect and introduces some approximations to the calculation. The authors made a comparison between the phantom and various reconstruction methods in their earlier paper [21]. The authors intend in this way to show the level of imperfection of the image relative to the real object. An important aspect to pay attention to when performing reduction is the selection of the k parameter. Its imprecise selection can significantly reduce the effectiveness of the method. The need for experimental selection of the parameter k is a drawback of the model; therefore, further research based on varied real data should parametrize this parameter for different materials. A second filtering step based on EMD introduces changes to all pixels in the image, as shown by the comparative results for very slight artifacts (less than 5%). The authors want to conduct further research to develop this method to provide it with more selective performance. Performing a test of different wavelet families in wavelet decomposition showed that the complexity of the wavelet has a negative effect on the effectiveness of artifact reduction. The best results were obtained using the following wavelets: Haar, Coiflet 1, Biorthogonal 2.2 and 3.1 and Daubechies 2. The lowest efficiency was obtained by performing artifact reduction with the “Discrete” Meyers wavelet. Greater variation in the results for individual wavelets is expected when reducing ring artifacts in real images due to the greater complexity of such images. This type of work is planned in further work on artifact reduction in tomographic images. The authors are in the process of testing the method on real data obtained from the micro-CT Skyscan 1172, and the results of this work will be presented in the authors’ next article.

Author Contributions

Conceptualization, P.L., M.B. and A.D.-G.; methodology, P.L. and M.B.; software, P.L.; validation, P.L. and M.B.; formal analysis, M.B. and A.D.-G.; investigation, P.L.; resources, P.L.; data curation, P.L.; writing—original draft preparation, P.L.; writing—review and editing, M.B. and A.D.-G.; visualization, P.L.; supervision, M.B.; project administration, P.L.; funding acquisition, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Science and Higher Education of Poland under research project No. WI/WM-IIB/8/2023.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This paper did not use external data; phantom images and artifacts were generated in software. Source codes can be made available by email after contacting the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kyriakou, Y.; Prell, D.; Kalender, W.A. Ring artifact correction for high-resolution micro CT. Phys. Med. Biol. 2009, 54, N385–N391. [Google Scholar] [CrossRef] [PubMed]
  2. Anas, E.M.; Kim, J.G.; Lee, S.Y.; Hasan, K. Comparison of ring artifact removal methods using flat panel detector based CT images. Biomed. Eng. Online 2011, 10, 72. [Google Scholar] [CrossRef] [PubMed]
  3. Kalender, W.A.; Kyriakou, Y. Flat-detector computed tomography (FD-CT). Eur. Radiol. 2007, 17, 2767–2779. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, Y.; Baek, J.; Hwang, D. Ring artifact correction using detector line-ratios in computed tomography. Opt. Express 2014, 22, 13380–13392. [Google Scholar] [CrossRef]
  5. Çimen, T.; Düzgün, S.; Akyüz, İ.E.; Topçuoğlu, H.S. The effect of cone beam computerized tomography voxel size and the presence of root filling on the assessment of middle mesial canals in mandibular molar teeth. Clin. Oral Investig. 2024, 28, 394. [Google Scholar] [CrossRef] [PubMed]
  6. Hiller, J.; Hornberger, P. Measurement accuracy in X-ray computed tomography metrology: Toward a systematic analysis of interference effects in tomographic imaging. Precis. Eng. 2016, 45, 18–32. [Google Scholar] [CrossRef]
  7. Pauwels, R.; Stamatakis, H.; Bosmans, H.; Bogaerts, R.; Jacobs, R.; Horner, K.; Tsiklakis, K. Quantification of metal artifacts on cone beam computed tomography images. Clin. Oral Implant. Res. 2013, 24, 94–99. [Google Scholar] [CrossRef] [PubMed]
  8. Kwan, A.L.; Seibert, J.A.; Boone, J.M. An improved method for flat-field correction of flat panel X-ray detector. Med. Phys. 2006, 33, 391–393. [Google Scholar] [CrossRef] [PubMed]
  9. Lifton, J.; Liu, T. Ring artefact reduction via multi-point piecewise linear flat field correction for X-ray computed tomography. Opt. Express 2019, 27, 3217–3228. [Google Scholar] [CrossRef]
  10. Rivers, M. Tutorial Introduction to X-ray Computed Microtomography Data Processing. Available online: https://www.mcs.anl.gov/research/projects/X-ray-cmt/rivers/tutorial.html (accessed on 7 December 2023).
  11. Munch, B.; Trtik, P.; Marone, F.; Stampanoni, M. Stripe and ring artifact removal with combined wavelet—Fourier filtering. Opt. Express 2009, 17, 8567–8591. [Google Scholar] [CrossRef]
  12. Tang, X.; Ning, R.; Yu, R.; Conover, D. Cone beam volume CT image artifacts caused by defective cells in x-ray flat panel imagers and the artifact removal using a wavelet-analysis-based algorithm. Med. Phys. 2001, 28, 812–825. [Google Scholar] [CrossRef] [PubMed]
  13. Sijbers, J.; Postnov, A. Reduction of ring artefacts in high resolution micro-CT reconstructions. Phys. Med. Biol. 2004, 49, N247. [Google Scholar] [CrossRef] [PubMed]
  14. Selim, M.; Rashed, E.A.; Atiea, M.A.; Kudo, H. Sparsity-based method for ring artifact elimination in computed tomography. PLoS ONE 2022, 17, e0268410. [Google Scholar] [CrossRef] [PubMed]
  15. Anas, E.M.A.; Lee, S.Y.; Hasan, M.K. Removal of ring artifacts in CT imaging through detection and correction of stripes in the sinogram. Phys. Med. Biol. 2010, 55, 6911. [Google Scholar] [CrossRef] [PubMed]
  16. Yu, H.; Ye, Y.; Wang, G. Katsevich-type algorithims for variable radius spiral cone-beam CT. In Proceedings of the Developments in X-Ray Tomography IV, Bellingham, WA, USA, 26 October 2004; pp. 550–557. [Google Scholar]
  17. Biguri, A.; Dosanjh, M.; Hancock, S.; Soleimani, M. TIGRE: A MATLAB-GPU toolbox for CBCT image reconstruction. Biomed. Phys. Eng. Express 2016, 2, 055010. [Google Scholar] [CrossRef]
  18. Hunter, J.D. Matplotlib: A 2D Graphics Environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  19. Lee, G.; Gommers, R.; Waselewski, F.; Wohlfahrt, K.; O’Leary, A. PyWavelets: A Python package for wavelet analysis. J. Open Source Softw. 2019, 4, 1237. [Google Scholar] [CrossRef]
  20. Quinn, A.; Lopes-Dos-Santos, V.; Dupret, D.; Nobre, A.; Woolrich, M. EMD: Empirical Mode Decomposition and Hilbert-Huang Spectral Analyses in Python. J. Open Source Softw. 2021, 6, 2977. [Google Scholar] [CrossRef]
  21. Lipowicz, P.; Dardzińska-Głębocka, A.; Borowska, M.; Biguri, A. Comparison of Analytical and Iterative Algorithms for Reconstruction of Microtomographic Phantom Images and Rat Mandibular Scans. In Information Technology in Biomedicine, Proceedings of the 9th International Conference, Kamień Śląski, Poland, 20–22 June 2022; Springer Nature: Cham, Switzerland, 2022; pp. 107–118. [Google Scholar]
  22. Mallat, S.; Hwang, W.L. Singularity detection and processing with wavelets. IEEE Trans. Inf. Theory 1992, 38, 617–643. [Google Scholar] [CrossRef]
  23. Toufik, B.; Mokhtar, N. The wavelet transform for image processing applications. Adv. Wavelet Theory Their Appl. Eng. Phys. Technol. 2012, 17, 395–422. [Google Scholar]
  24. Boin, M.; Haibel, A. Compensation of ring artefacts in synchrotron tomographic images. Opt. Express 2006, 14, 12071–12075. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  26. Deering, R.; Kaiser, J.F. The use of a masking signal to improve empirical mode decomposition. In Proceedings of the (ICASSP ‘05) IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005, Philadelphia, PA, USA, 23–23 March 2005; Volume 484, pp. iv/485–iv/488. [Google Scholar]
  27. Asamoah, D.; Ofori, E.; Opoku, S.; Danso, J. Measuring the performance of image contrast enhancement technique. Int. J. Comput. Appl. 2018, 181, 6–13. [Google Scholar] [CrossRef]
  28. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE Trans. Image Process 2004, 13, 600–612. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, Z.; Simoncelli, E.P.; Bovik, A.C. Multiscale structural similarity for image quality assessment. In Proceedings of the Thrity-Seventh Asilomar Conference on Signals, Systems & Computers, 2003, Pacific Grove, CA, USA, 9–12 November 2003; pp. 1398–1402. [Google Scholar]
Figure 1. Schematic representation of the structure of the microtomograph’s CCD detector, with the X-ray and visible-spectrum radiation areas marked.
Figure 1. Schematic representation of the structure of the microtomograph’s CCD detector, with the X-ray and visible-spectrum radiation areas marked.
Applsci 14 07292 g001
Figure 2. Reference image of the Shepp–Logan phantom (a) and image showing the Shepp–Logan phantom with ring artifacts inserted (b).
Figure 2. Reference image of the Shepp–Logan phantom (a) and image showing the Shepp–Logan phantom with ring artifacts inserted (b).
Applsci 14 07292 g002
Figure 3. The image of the sinogram after wavelet decomposition using the Haar wavelet (a) and the result of the binarization process (b).
Figure 3. The image of the sinogram after wavelet decomposition using the Haar wavelet (a) and the result of the binarization process (b).
Applsci 14 07292 g003
Figure 4. A graph showing the N k j matrix after summing the B k n , j matrix relative to the columns. The T line represents the thresholding level to identify ring artifact stripes.
Figure 4. A graph showing the N k j matrix after summing the B k n , j matrix relative to the columns. The T line represents the thresholding level to identify ring artifact stripes.
Applsci 14 07292 g004
Figure 5. A graph showing the result of comparing images of the reference Shepp–Logan phantom and the reconstructed noisy and filtered image of the Shepp–Logan phantom using the RMSE (red line) and MSE (black line) methods for changing parameter k.
Figure 5. A graph showing the result of comparing images of the reference Shepp–Logan phantom and the reconstructed noisy and filtered image of the Shepp–Logan phantom using the RMSE (red line) and MSE (black line) methods for changing parameter k.
Applsci 14 07292 g005
Figure 6. Images showing Shepp–Logan phantom reconstructions (a) without ring artifact reduction (stripe level 10%) and (b) after performing ring artifact reduction with constant k = 21, thresholding level T = 0.4, and span factor N = 3.
Figure 6. Images showing Shepp–Logan phantom reconstructions (a) without ring artifact reduction (stripe level 10%) and (b) after performing ring artifact reduction with constant k = 21, thresholding level T = 0.4, and span factor N = 3.
Applsci 14 07292 g006
Figure 7. Graph showing results of comparative analysis of images after reconstruction using wavelets of different families and reference image using RMSE, SSIM and MS-SSIM methods.
Figure 7. Graph showing results of comparative analysis of images after reconstruction using wavelets of different families and reference image using RMSE, SSIM and MS-SSIM methods.
Applsci 14 07292 g007
Table 1. A table showing papers on methods for reducing ring artifacts, along with the type of data on which the method was tested and example results.
Table 1. A table showing papers on methods for reducing ring artifacts, along with the type of data on which the method was tested and example results.
Authors (year)MethodsDataResults
Selim et al. (2022) [14]Sparsity-based method, during iterative reconstruction algorithmsReal data and simulations0.9682 (SSIM)
Sijbers et al. (2004) [13]Post-reconstruction filtration method Real dataOnly visulal results, no quantitative data
Tang et al. (2001) [12]Wawelet analysis pre-reconstruction method Callibration cellsOnly visulal results, no quantitative data
Kim et al. (2014) [4] Pre-reconstruction method of calculating the ratio of adjacent detector elements in the projection data, termed the line ratioSheep–Logan phantom and human organs models 7.13 × 10−7 (MSE) for Sheep-Logan phantom
Abu et al. (2011) [2]Comparison five ring artifact removal methods: Modified Wavelet (MWPN) [12], Wavelet-Fourier (WF) method [11], Ring corrections using homogeneity test (RCHT) [13], Ring correction in polar coordinate (RCP) [1], and Strength based ring correction (SBRC) method [15]Head Phantom and real dataquantitative data for head Phantom (MSSIM):
MWPN—0.97
WF—0.81
RCHT—0.82
RCP—0.77
SBRC—0.98
Table 2. Values of calculated similarity coefficients of the post-filtered and noisy unfiltered images of the Sheep-Logan phantom to the reference image.
Table 2. Values of calculated similarity coefficients of the post-filtered and noisy unfiltered images of the Sheep-Logan phantom to the reference image.
PSNR (dB)MSESSIMMS-SSIM
No Filtration 20% noise21.992880.070740.376110.60218
Filtered 20% noise27.091040.021870.800900.82214
No Filtration 10% noise27.087070.021890.476120.76688
Filtered 10% noise29.353440.012990.802260.82219
No Filtration 5% noise34.590040.003890.611200.80143
Filtered 5% noise27.085090.021900.802650.82261
No Filtration 2% noise37.747950.001880.794280.88367
Filtered 2% noise27.093030.021860.805340.82340
Noise-free reconstruction 38.614330.001540.921640.98685
Table 3. Table showing results of comparative analysis of reconstructed images using wavelets of different families and reference image using RMSE, SSIM and MS-SSIM methods. The results with the highest similarity values to the reference image are shown in bold.
Table 3. Table showing results of comparative analysis of reconstructed images using wavelets of different families and reference image using RMSE, SSIM and MS-SSIM methods. The results with the highest similarity values to the reference image are shown in bold.
Wavelets TypePSNR (dB)RMSESSIMMS-SSIM
Haar29.04710.14770.80310.8231
“Discrete” Meyer28.97290.15020.61400.7211
Coiflets 129.07100.14690.78980.8155
Coiflets 529.02790.14830.63370.7244
Coiflets 1029.,04200.14790.65330.7484
Coiflets 1529.03200.14820.64340.7320
Biorthogonal 1.529.07430.14680.67080.7485
Biorthogonal 2.229.06600.14700.76310.8052
Biorthogonal 2.829.07130.14690.65310.7320
Biorthogonal 3.129.05130.14750.78330.8125
Biorthogonal 3.929.08860.14630.66010.7470
Biorthogonal 4.429.10730.14570.71110.7823
Biorthogonal 5.529.04920.14760.65430.7451
Biorthogonal 6.829.03640.14810.62950.7083
Daubechies 229.04150.14790.77000.8080
Daubechies 1029.03530.14810.66330.7551
Daubechies 2029.01820.14870.62930.7366
Daubechies 3028.99430.14950.66040.7428
Symlets 429.08080.14650.72760.7685
Symlets 829.06520.14710.67020.7565
Symlets 1229.01400.14880.61830.7119
Symlets 1629.03170.14820.62430.7255
Symlets 2029.06330.14710.63090.7397
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lipowicz, P.; Borowska, M.; Dardzińska-Głębocka, A. Methods for Reducing Ring Artifacts in Tomographic Images Using Wavelet Decomposition and Averaging Techniques. Appl. Sci. 2024, 14, 7292. https://doi.org/10.3390/app14167292

AMA Style

Lipowicz P, Borowska M, Dardzińska-Głębocka A. Methods for Reducing Ring Artifacts in Tomographic Images Using Wavelet Decomposition and Averaging Techniques. Applied Sciences. 2024; 14(16):7292. https://doi.org/10.3390/app14167292

Chicago/Turabian Style

Lipowicz, Paweł, Marta Borowska, and Agnieszka Dardzińska-Głębocka. 2024. "Methods for Reducing Ring Artifacts in Tomographic Images Using Wavelet Decomposition and Averaging Techniques" Applied Sciences 14, no. 16: 7292. https://doi.org/10.3390/app14167292

APA Style

Lipowicz, P., Borowska, M., & Dardzińska-Głębocka, A. (2024). Methods for Reducing Ring Artifacts in Tomographic Images Using Wavelet Decomposition and Averaging Techniques. Applied Sciences, 14(16), 7292. https://doi.org/10.3390/app14167292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop