Next Article in Journal
Image Features Based on Characteristic Curves and Local Binary Patterns for Automated HER2 Scoring
Previous Article in Journal
Partition and Inclusion Hierarchies of Images: A Comprehensive Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm

Faculty of Engineering, University of Rijeka, Vukovarska 58, 51000 Rijeka, Croatia
*
Author to whom correspondence should be addressed.
J. Imaging 2018, 4(2), 34; https://doi.org/10.3390/jimaging4020034
Submission received: 5 December 2017 / Revised: 1 February 2018 / Accepted: 2 February 2018 / Published: 5 February 2018

Abstract

:
Diagnostics and treatments of numerous diseases are highly dependent on the quality of captured medical images. However, noise (during both acquisition and transmission) is one of the main factors that reduce their quality. This paper proposes an adaptive image denoising algorithm applied to enhance X-ray images. The algorithm is based on the modification of the intersection of confidence intervals (ICI) rule, called relative intersection of confidence intervals (RICI) rule. For each image pixel apart, a 2D mask of adaptive size and shape is calculated and used in designing the 2D local polynomial approximation (LPA) filters for noise removal. One of the advantages of the proposed method is the fact that the estimation of the noise free pixel is performed independently for each image pixel and thus, the method is applicable for easy parallelization in order to improve its computational efficiency. The proposed method was compared to the Gaussian smoothing filters, total variation denoising and fixed size median filtering and was shown to outperform them both visually and in terms of the peak signal-to-noise ratio (PSNR) by up to 7.99 dB.

1. Introduction

Denoising (both of one-dimensional and multi-dimensional signals) is one of the most important fields in signal processing and hence, the pursuit for novel and more efficient denoising methods is constant [1]. Numerous methods developed over the last few decades found their applications in various fields, including medical image processing and analysis leading to enhancements in examining the interior of the human body without surgeries. However, since both revealing and treatment of large number of diseases today is highly dependent on medical images (such as computed tomography (CT), ultrasound, magnetic resonance imaging (MRI), X-rays etc.), improving their quality is an essential precondition for their analysis. One of the main problems with captured medical images is the presence of noise (introduced during both acquisition and transmission processes), which complicates their visual inspection and pathology pattern recognition. Thus, noise has to be suppressed in order to decrease the probability of possible misinterpretations and incorrect diagnoses [2].
Due to their cost and a decrease in the dose of ionizing radiation, X-rays are still the most frequently used medical imaging technique. However, lower doses of radiation require more efficient noise reduction methods such that the important image features (object contours, edges, textures, etc.) are preserved while the noise is smoothed.
There are numerous techniques for reducing noise from medical images, which can be classified into two main categories. First are the methods for image processing in the spatial domain and the second are the transform domain filtering approaches [3]. Spatial filtering methods may be divided into linear (such as, for example, mean and Wiener based filtering) and non-linear methods (for example, median and weighted median based filtering) [3]. However, the linear spatial domain filtering fails in case of signal dependent noise and, unfortunately, it also introduces blurring artefacts to denoised images [4,5]. On the other hand, transform domain image processing approaches may be divided into data adaptive methods (for example, independent component analysis) and non-data adaptive methods (such as wavelet domain and spatial-frequency domain processing algorithms) [3]. Unlike the spatial methods, the transform domain-based methods transform images from spatial into transform domain and noise is removed in the transform domain [4]. Next, the estimate of the noise-free image is obtained by applying the inverse transform.
The method proposed in the paper introduces an algorithm for medical image processing based on the adaptive, data-adaptive 2D spatial filters designed using the local polynomial approximation (LPA) combined with the modification of the intersection of confidence intervals (ICI) rule, called the relative intersection of confidence intervals (RICI) rule (where RICI rule determines the size of the LPA estimators).
The 2D LPA-RICI based algorithm has been shown to be an efficient noise-free image estimator, since the RICI based algorithm calculates the near optimal estimator size and its 2D shape for each considered pixel [6,7]. Unlike the approaches based on the fixed size filters, the given method (due to its couture and edge preserving property) smoothes the noise locally in the vicinity of the considered pixel and hence avoids blurring artefacts in the resulting denoised images.
Once calculated, adaptive 2D regions are used as binary masks for segmentation and extraction of the regions of interest from the original image. The noise-free estimate of the considered pixel is found by applying the LPA weighted averaging to the extracted region. The described procedure is repeated independently for each image pixel separately. Thus, one of the important advantages of the proposed 2D LPA-RICI based image denoising algorithm is that it may be easily parallelized in order to improve the method’s computational efficiency.
The proposed adaptive 2D LPA-RICI method was applied to denoising of real-life medical X-ray images, outperforming competitive methods (Gaussian smoothing filters, total variation denoising and fixed size median filtering) in most cases.
The paper is organized as follows. Section 2 introduces the original ICI rule and its modification (the RICI rule) extended to 2D image processing. Section 3 gives an elaboration on the achieved results. The conclusion is found in Section 4.

2. The ICI Rule and Its Modification

2.1. The ICI Algorithm

The LPA-ICI method utilizes the LPA for designing estimators, the size of which is data-driven and, hence, chosen adaptively using the ICI based rule. This section presents original one-dimensional ICI rule and its modification (the RICI rule) extended to denoising 2D medical images. A software for the original 1D ICI and 1D RICI signal denoising is made available by the author in [8].
Therefore, before we propose the 2D version of the RICI algorithm, let us introduce the original 1D signal filtering procedure.
Firstly, in order to obtain the estimate of the noise free signal, the ICI algorithm introduce a range of K estimators with increasing widths hi [6,7]:
H = {h1 < h2 < ⋯ < hK},
and the corresponding confidence intervals D h k ( n ) , 1 ≤ k ≤ K, defined for each signal sample n with interval limits (upper U h k ( n ) and lower L h k ( n )) defined as [6,7] (Katkovnik et al. 2005, 2002):
U h k ( n ) = y ^ h k ( n ) + Г σ h k ( n ) ,
L h k ( n ) = y ^ h k ( n ) Г σ h k ( n ) ,
where Г defines the confidence level, σ h k ( n ) represents the standard deviation of the estimation error and y ^ h k ( n ) is calculated as the LPA weighted average of k samples neighboring the considered n -th sample [6,7].
The procedure is repeated as long as all subsequent k intervals are overlapping (or until the edges of the signal are reached). Namely, the ICI algorithm tracks the values of the smallest upper U ¯ h k ( n ) and the largest lower limits L ¯ h k ( n ) of the confidence interval D h k ( n ) defined as [6,7]:
U ¯ h k ( n ) = min i = 1 , , k U h i ( n ) ,
L ¯ h k ( n ) = max i = 1 , , k L h i ( n ) ,
as long as the lower interval limit is equal to or smaller than the upper interval limit (meaning that all subsequent intervals are overlapping) [6,7]:
L ¯ h k ( n ) U ¯ h k ( n ) .
Finally, the optimal filter width is determined as the largest h k for which all previous k confidence intervals, including the k -th interval, are overlapping [6,7].
As shown in [9], small hk values increase the estimate error variance and decrease its bias. On the other hand, large hk values decrease the estimate error variance and increase its bias (dependent on the unknown high-order signal derivatives) [9]. The adaptive ICI based algorithm calculates the maximal hk, value which ensures the optimal trade-off between the estimation bias and the variance. In other words, proper Г value results in the optimal smoothing effect since it ensures selecting the largest vicinity neighboring the considered point such that the LPA fits well to data [9]. Namely, small hk values are chosen for points close to rapid changes in the signal, while otherwise larger hk values are selected [9].
In order to justify previous claims, let us consider the absolute estimation error e h k ( n ) , calculated as:
e h k ( n ) = | y ( n ) y ^ h k ( n ) | ,
where y ( n ) is a noise-free signal and y ^ h k ( n ) its LPA estimate. The estimation error can be written as a sum of the bias b h k ( n ) and the zero-mean random error e h k 0 ( n ) [9]:
e h k ( n ) = b h k ( n ) + e h k 0 ( n ) .
Hence, the following inequity stands true [9]:
e h k ( n ) | b ¯ h k ( n ) | + | e h k 0 ( n ) | ,
where | b ¯ h k ( n ) | is the maximal value of | b h k ( n ) | and e h k 0 ( n ) is, in case of the Gaussian noise, zero-mean estimation error with the standard deviation σ h k ( n ) . Furthermore, the following inequity holds true with the probability p = 1 β :
| e h k 0 ( n ) | χ 1 β / 2 · σ h k ( n ) ,
where χ 1 β / 2 is ( 1 β / 2 )-th qantile of the normal distribution N ( 0 , 1 ) [9]. In other words, the estimation error e h k 0 ( n ) is found inside the interval [ χ 1 β / 2 , + χ 1 β / 2 ] with the probability p . From Equations (9) and (10) it follows that [9]:
| e h k ( n ) | | b ¯ h k ( n ) | + χ 1 β / 2 · σ h k ( n )
holds true with the same probability p [9]. It was shown that there exists an optimal h k , denoted as h , providing an optimal estimation error bias to the variance trade-off, such that for h k h [9]:
| e h k ( n ) | Г · σ h k ( n ) .
Thus, the inequity (12) can be written as:
| y ( n ) y ^ h k ( n ) | Г · σ h k ( n ) ,
or as:
y ^ h k ( n ) Г · σ h k ( n ) y ( n ) y ^ h k ( n ) + Г · σ h k ( n ) ,
where y ^ h k ( n ) Г · σ h k ( n ) is lower and y ^ h k ( n ) + Г · σ h k ( n ) upper confidence interval limit, as given in Equations (2) and (3). Since both inequities (11) and (14) hold true with the probability p , we can claim with the same probability that y ( n ) is inside the interval y ^ h k ( n ) ± Г · σ h k ( n ) for h k h [9]. In other words, a nonempty intersection of all subsequent confidence intervals ensures h k h . On the other hand, an empty intersection of all subsequent confidence intervals means that h k > h [9]. Thus, the largest h k , as close as possible to h , for which the intersection of all subsequent confidence intervals is nonempty, is to be chosen as the one giving the near optimal balance between the estimation error bias and the variance, as shown in Equation (6).
An illustration of tracking of confidence intervals by the ICI algorithm is shown in Figure 1. As it can be seen, the first three confidence intervals ( D 1 , D 2   i   D 3 ) have common points, i.e., their intersection is nonempty. On the other hand, D 4 is overlapping with D 2 and D 3 but since it does not overlap with all previous confidence intervals (it is not overlapping with D 1 ), the ICI rule results in ℎ3 being chosen as the optimal estimator width.

2.2. The RICI Algorithm

As shown in Equations (2) and (3), the ICI algorithm is highly affected by the selected value of the Г parameter. Namely, too large Г values result in oversmoothing (i.e., noise is removed with object edges being blurred) and too small Г values cause signal undersmoothing (i.e., significant amount of noise is left in the image) [6,7].
Thus, we propose a modification of the ICI rule (further extended to 2D for medical images denoising) shown to be significantly more robust to the suboptimal Г selections [10,11]. The method is based on tracking the value of the ratio (and thus is called the RICI rule) of the overlapping of the confidence intervals and the length of the current confidence interval, defined as [10,11]:
R h k ( n ) = U ¯ h k ( n ) L ¯ h k ( n ) U h k ( n ) L h k ( n ) .
Next, R h k ( n ) , given in Equation (15), is then used as an additional criterion for selecting the adaptive filter width [10,11]:
R h k ( n ) R c ,
where R c stands for the preset data-driven threshold value [10,11].
The RICI rule allows selecting larger Г values (ensuring efficient noise removal), while additional criterion Equation (16) protects from blurring artefacts. Hence, the RICI algorithm smoothes the noise and, at the same time, preserves object contours and edges in the denoised image.

2.3. Extending the ICI and RICI Algorithm to Image Processing

The above described ICI and RICI based algorithms for 1D signal processing were upgraded for detecting two-dimensional regions of interest in the vicinity of the considered image pixel.
Namely, the lower and upper interval limits defined in Equations (2) and (3) were extended to two dimensions (resulting in calculating U h k ( i ,   j ) and L h k ( i ,   j ) ), where i and j define pixel indexes.
The size and the shape of the 2D regions, used in the estimation of the noise-free image pixel, depends on image content and varies from pixel to pixel. Once detected, denoised pixel value is estimated using the 2D LPA based weighted averaging of the pixels in the detected region. The procedure is repeated for each image pixel separately. Namely, in this paper the regions were calculated using two, four and eight lines intersecting in the considered pixel forming quadrilateral, octagonal and hexadecagonal regions, respectively. Examples of adaptive quadrilateral, octagonal and hexadecagonal regions are shown in Figure 2.
The proposed 2D LPA-RICI image denoising method is planned to be published in an extensive virtual instrument together with other methods for image denoising as a part of our future work (as it was case for the ICI and RICI based 1D signal denoising methods made available in [8]).
The next section elaborates on the obtained results in X-ray image denoising.

3. Results and Discussion

This section presents denoising results achieved using the proposed 2D LPA-RICI method applied to three test X-ray images (namely, X-ray scans of shoulder, chest and ankle). The method was implemented in the Matlab 2015b and denoising performances were measured in terms of the peak signal-to-noise ratio (PSNR).
Shoulder X-ray images (with resolution 1024 × 1018) are presented in Figure 3. Namely, the noise-free image is given in Figure 3a. The noisy image is shown in Figure 3b (corrupted by additive white Gaussian noise (AWGN) the standard deviation of which is σ = 25). Figure 3c–e present images denoised using the RICI based method. Namely, Figure 3c is denoised using the quadrilateral regions, Figure 3d is denoised using the octagonal regions and Figure 3e is denoised using the hexadecagonal regions.
Denosing results, in terms of the PSNR, for the shoulder X-ray image are given in Table 1 (for various Г and Rc values). The first column gives regions (octagonal, octagonal and hexadecagonal, respectively). The second and the third column present parameter values for the 2D LPA-RICI denoising method. The fourth column presents the standard deviation of the AWGN. The fifth column provides the noisy image PSNR followed by the column giving the PSNR for image denoised using the proposed RICI based method. The last three columns of the Table 1 show the PSNRs of the images denoised using fixed size filtering (2D median filter), Gaussian smoothing filters and total variation denoising, respectively. For the Gaussian smoothing filtering a Matlab function imgaussfilt was used with standard deviation values chosen same as the standard deviations of the AWGN (σ = 20, 25 and 30). For the total variation denoising [12] Matlab function provided at the MathWorks was used with parameter λ chosen as λ = 0.2, λ = 0.25 and λ = 0.3 for the three tested AWGN levels (σ = 20, 25 and 30), respectively.
Table 2 gives PSNR improvements obtained using the 2D LPA-RICI method vs. noisy image and images denoised using the fixed size filtering, Gaussian smoothing filters and total variation method, respectively. As shown in Table 2, the proposed adaptive RICI based method outperforms the fixed size filters in all cases. Furthermore, it outperforms Gaussian smoothing filters and total variation denoising when quadrilateral and octagonal regions are used (hexadecagonal regions do not outperform Gaussian smoothing and total variation denoising methods).
Namely, the 2D LPA-RICI with quadrilateral regions increased the PSNR of the denoised image by up to 9.13 dB, with octagonal regions increased the PSNR by up to 10.50 dB and with hexadecagonal regions increased the PSNR by up to 4.69 dB (when compared to the noisy image). Also, the proposed 2D LPA-RICI method outperformed the fixed size 2D median filtering by up to 7.99 dB, the Gaussian smoothing filters by up to 6.08 dB and the total variation denoising by up to 4.35 dB for shoulder X-ray scan.
Denoising time for the fixed size filtering was up to 1 s and up to 3 s for the Gaussian smoothing filter and up to 7 min for the total variation denoising. However, the denoising time for the 2D LPA-RICI method without parallelization was up to 40 min for shoulder X-ray scan, depending on the used region. Namely, denoising time for quadrilateral region was up to 24 min, for octagonal region up to 34 min and for hexadecagonal region up to 40 min (greater number of polygonal angles leads to longer execution time).
Denoising of all test images was performed on a personal computer Dell Inspiron 15R 5521 with Intel processor i5 d generation, 6 GB RAM DDR3, 1 TB HDD SATA 1.5 Gbps and Windows 8.1 (64-bit) operating system.
The results for the second test image (chest X-ray image with the resolution 2431 × 1782) are given in Figure 4. Namely, Figure 4a,b shows the noise-free and the noisy image (also corrupted with AWGN with σ = 25), respectively. Figure 4c–e present denoised images obtained using the proposed 2D LPA-RICI algorithm with quadrilateral, octagonal and hexadecagonal regions, respectively.
The PSNR results for denoised chest X-ray images are found in Table 3. As it can be seen from Table 3, the proposed 2D LPA-RICI method (as it was the case for the shoulder X-ray image) outperformed fixed size 2D filtering in all cases. Furthermore, the RICI based denoising outperformed Gaussian smoothing filters and total variation denoising when quadrilateral and octagonal regions were used.
Table 4 gives PSNR improvements for chest X-ray scan denoised using the 2D LAP-RICI method compared to noisy images and images denoised using the fixed size filtering, Gaussian smoothing filters and the total variation method, respectively. Namely, the LPA-RICI method enhanced denoised image quality, both visually and in terms of the PSNR, by up to 9.40 dB when quadrilateral regions were used, by up to 10.23 dB when octagonal regions were used and by up to 5.25 dB when hexadecagonal regions were used (when compared to the noisy image PSNR). It also outperforms fixed size 2D median filtering in all cases, increasing the PSNR by up to 7.51 dB. In addition, the 2D LPA-RICI method outperformed the Gaussian smoothing filters by up to 5.35 dB when quadrilateral and octagonal regions were used. Furthermore, it also outperformed the total variation denoising increasing the PSNR by up to 3.30 dB (in case of the adaptive 2D quadrilateral and octagonal regions).
Denoising time for the chest X-ray scans for the fixed size filtering was up to 1.5 s and for the Gaussian smoothing filter up to 6 s. Furthermore, the denoising time for the total variation method was up to 10.5 min and for the 2D LPA-RICI without parallelization up to 60 min, depending on the used regions (up to 36 min for quadrilateral regions, up to 51 min for octagonal regions and up to 60 min for hexadecagonal regions).
Denoised ankle X-ray images (with resolution 600 × 682) are given in Figure 5. Figure 5a,b show noise-free and noisy images, respectively. Figure 5c–e shows the denoised images obtained using the 2D LPA-RICI method with quadrilateral, octagonal and hexadecagonal regions, respectively.
Table 5 gives the PSNRs of the denoised ankle X-ray images. The proposed 2D LPA-RICI method, as it was the case for other two test images, outperformed fixed size 2D filtering in all cases. Furthermore, it also increased the denoised image PSNR when compared to the Gaussian smoothing filters and the total variation denoising (when quadrilateral or octagonal regions were used).
As it can be seen from Table 6, ankle X-ray scan denoised by the 2D LPA-RICI method with quadrilateral regions resulted in an increased PSNR (compared to the noisy image) by up to 8.09 dB, by up to 7.46 dB with octagonal regions and by up to 3.93 dB with hexadecagonal regions. Once again, the proposed 2D LPA-RICI method was shown to significantly outperform the fixed size 2D median filtering. In addition, it also outperformed Gaussian smoothing filter method for all tested regions (including hexadecagonal regions, which was not the case in the previous examples) by up to 5.46 dB. Furthermore, it increased the denoised PSNR when compared to the total variation denoising method by up to 2.80 dB (when quadrilateral and octagonal regions were used).
Denoising of the ankle X-ray scan took by up to 0.5 s for fixed size filtering, by up to 2 s for the Gaussian smoothing filtering and by up to 5 min for the total variation denoising. Denoising using the LPA-RICI method without parallelization took up to 23 min (for quadrilateral region up to 15.5 min, for octagonal region up to 26 min and for hexadecagonal region up to 32 min).
The provided study showed that the proposed 2D LPA-RICI method performs competitively in terms of medical image denoising due to its adaptivity to image content and calculation of the 2D regions neighbouring the considered pixel for each pixel independently. Furthermore, it was shown that the proposed method outperforms fixed size filters for all tested cases (as shown in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6). Also, the proposed 2D LPA-RICI method outperforms Gaussian smoothing filters and total variation denoising when adaptive quadrilateral and octagonal regions were used. Furthermore, the obtained results suggest that increasing the number of polygonal angles does not necessarily result in improved estimation accuracy (quadrilateral and octagonal regions performed better than the hexadecagonal regions). However, a larger number of polygonal angles increases computational burden of the 2D LPA-RICI method, thus the study suggests using the 2D LPA-RICI with quadrilateral and octagonal regions in medical imaging applications.
Note also that greater improvements may be obtained in cases of higher resolution images (as it is case for images shown in Figure 3 and Figure 4 compared to the ankle X-ray image in Figure 5).
In addition, since denoising of each image pixel is done independently, the method is applicable for easy parallelization which would reduce its computational expenses.

4. Conclusions

The paper introduces an adaptive denoising method applied to medical image denoising. The proposed method is based on a modification of the ICI algorithm, called the RICI algorithm, extended to local 2D LPA based filtering. The paper shows that the proposed 2D LPA-RICI method performs competitively in terms of medical image denoising due to its adaptivity to image content by readjusting the size and shape of 2D regions from pixel to pixel independently. Thus, the proposed method is easily parallelized for reducing its computational expenses.
The 2D LPA-RICI with quadrilateral, octagonal and hexadecagonal regions was tested in the paper, showing that increasing the number of polygonal angles does not necessarily result in an improved denoising accuracy. On the other hand, a larger number of polynomial angles increases the method’s computational burden.
The method was tested on real-life X-ray images and compared to the fixed size 2D median filtering, Gaussian smoothing filters and the total variation denoising. The 2D LPA-RICI method increased the image PSNR by up to 10.50 dB when compared to the noisy image PSNR. Furthermore, it outperformed the fixed size filtering by up to 7.99 dB, as well as the Gaussian smoothing filters when quadrilateral and octagonal regions were used by up to 6.08 dB. Furthermore, the 2D LPA-RICI method with quadrilateral and octagonal regions outperformed the total variation denoising by up to 4.35 dB.

Acknowledgments

This work has been supported in part by the University of Rijeka under the project number 16.09.2.2.04.

Author Contributions

Ivica Mandić has implemented 2D LPA-RICI method in Matlab and applied it to test images. He contributed in writing draft of the paper as well as in comparing the results obtained by the proposed method to the total variation denoising. Hajdi Peić has implemented adaptive regions (quadrilateral, octagonal and hexadecagonal) in Matlab and analysed PSNRs with respect to tested regions. She also contributed in writing draft of the paper and comparing the results obtained using the proposed method to Gaussian and fixed size filtering. Jonatan Lerga has proposed extending his LPA-RICI method to the 2D LPA-RICI method. He rewrote numerous sections in order to improve the paper’s readability. Furthermore, as a corresponding author, he wrote reply to reviewers’ comments. Ivan Štajduhar contributed in finalizing the revised version of the paper modified according to the reviewers’ suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Buades, A.; Coll, B.; Morel, J.M. On Image Denoising Methods. CMLA Preprint. 2004, Volume 5. Available online: http://www3.cs.stonybrook.edu/~cse577/image.smoothing.morel.pdf (accessed on 4 February 2018).
  2. Bansal, M.; Devi, M.; Jain, N.; Kukreja, C. A proposed approach for biomedical image denoising using PCA-NLM. Int. J. Bio-Sci. Bio-Technol. 2014, 6, 13–20. [Google Scholar]
  3. Samundeeswari, E.S.; Saranya, P.K.; Manavalan, R. M2 filter for speckle noise suppression in breast ultrasound images. ICTACT J. Image Video Process. 2015, 6, 1137–1144. [Google Scholar]
  4. Vijay, M.; Suhha, S.V.; Karthik, K. Adaptive spatial and wavelet multiscale products thresholding method for medical image denoising. In Proceedings of the International Conference on Computing, Electronics and Electrical Technologies, Kumaracoil, India, 21–22 March 2002; pp. 1109–1114. [Google Scholar]
  5. Ali, H.M. MRI medical image denoising by fundamental filters. SCIREA J. Comput. 2017, 12, 2–26. [Google Scholar]
  6. Katkovnik, V.; Egiazarian, K.; Astola, J. A spatially adaptive nonparametric regression image deblurring. IEEE Trans. Image Process. 2005, 14, 1469–1478. [Google Scholar] [CrossRef] [PubMed]
  7. Katkovnik, V.; Egiazarian, K.; Astola, J. Adaptive window size image de-noising based on intersection of confidence intervals (ICI) rule. J. Math. Imaging Vis. 2002, 16, 223–235. [Google Scholar] [CrossRef]
  8. Segon, G.; Lerga, J.; Sucic, V. Improved LPA-ICI-based estimators embedded in a signal denoising virtual instrument. Signal Image Video Process. 2017, 11, 211–218. [Google Scholar] [CrossRef]
  9. Katkovnik, V.; Egiazarian, K.; Astola, J. Local Approximation Techniques in Signal and Image Processing; SPIE Publications—The International Society for Optical Engineering: Washington, WA, USA, 2006. [Google Scholar]
  10. Lerga, J.; Vrankic, M.; Sucic, V. A signal denoising method based on the improved ICI rule. IEEE Signal Process. Lett. 2008, 15, 601–604. [Google Scholar] [CrossRef]
  11. Lerga, J.; Sucic, V.; Grbac, E. An adaptive method for video denoising based on the ICI rule. Eng. Rev. 2012, 32, 33–40. [Google Scholar]
  12. Rudin, L.I.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D Nonlinear Phenom. 1992, 60, 259–268. [Google Scholar] [CrossRef]
Figure 1. Illustration of the ICI rule tracking the intersection of the confidence intervals.
Figure 1. Illustration of the ICI rule tracking the intersection of the confidence intervals.
Jimaging 04 00034 g001
Figure 2. Examples of the detected regions obtained using the proposed method. (a) Quadrilateral region; (b) Octagonal region; (c) Hexadecagonal region.
Figure 2. Examples of the detected regions obtained using the proposed method. (a) Quadrilateral region; (b) Octagonal region; (c) Hexadecagonal region.
Jimaging 04 00034 g002
Figure 3. Shoulder X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Figure 3. Shoulder X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Jimaging 04 00034 g003aJimaging 04 00034 g003b
Figure 4. Chest X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Figure 4. Chest X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Jimaging 04 00034 g004aJimaging 04 00034 g004b
Figure 5. Ankle X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Figure 5. Ankle X-ray scan. (a) Original noise-free image; (b) Noisy image (AWGN with σ = 25); (c) Image denoised using the 2D LPA-RICI method (quadrilateral region, Г = 1.8, Rc = 0.8); (d) Image denoised using the 2D LPA-RICI method (octagonal region, Г = 1.8, Rc = 0.8); (e) Image denoised using the 2D LPA-RICI method (hexadecagonal region, Г = 1.8, Rc = 0.8); (f) Image denoised using Gaussian smoothing filters; (g) Image denoised using total variation denoising.
Jimaging 04 00034 g005
Table 1. PSNRs for shoulder X-ray images denoised using the proposed 2D LPA-RICI method and competitive methods.
Table 1. PSNRs for shoulder X-ray images denoised using the proposed 2D LPA-RICI method and competitive methods.
Parameters PSNR (dB)
RegionГRcσNoisy ImageRICIFixed Size FiltersGaussian SmoothingTotal Variation
Quadrilateral1.80.82029.8238.7333.434.6935.02
Octagonal1.80.82029.8239.3733.434.6935.02
Hexadecagonal1.80.82029.8233.6133.434.6935.02
Quadrilateral1.70.92029.8238.6333.434.6935.02
Octagonal1.70.92029.8239.133.434.6935.02
Hexadecagonal1.70.92029.8233.4233.434.6935.02
Quadrilateral1.60.62029.8238.7233.434.6935.02
Octagonal1.60.62029.8239.1933.434.6935.02
Hexadecagonal1.60.62029.8233.533.434.6935.02
Quadrilateral1.80.82529.3238.4532.2433.9936.69
Octagonal1.80.82529.3238.9632.3433.9936.69
Hexadecagonal1.80.82529.3233.5132.3433.9936.69
Quadrilateral1.70.92529.3238.1732.2433.9936.69
Octagonal1.70.92529.3238.8232.2433.9936.69
Hexadecagonal1.70.92529.3233.532.2433.9936.69
Quadrilateral1.60.62529.3238.2732.2433.9936.69
Octagonal1.60.62529.3238.8432.2433.9936.69
Hexadecagonal1.60.62529.3233.5732.2433.9936.69
Quadrilateral1.80.83028.9938.0931.533.4135.62
Octagonal1.80.83028.9938.7631.533.4135.62
Hexadecagonal1.80.83028.9933.6831.533.4135.62
Quadrilateral1.70.93028.9937.7831.533.4135.62
Octagonal1.70.93028.9938.5831.533.4135.62
Hexadecagonal1.70.93028.9933.431.533.4135.62
Quadrilateral1.60.63028.9937.8231.533.4135.62
Octagonal1.60.63028.9939.4931.533.4135.62
Hexadecagonal1.60.63028.9933.6631.533.4135.62
Table 2. PSNR improvements for shoulder X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
Table 2. PSNR improvements for shoulder X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
Parameters PSNR (dB)
RegionГRcσRICI vs. NoisyRICI vs. FixedRICI vs. GAUSSIANRICI vs. Total Variation
Quadrilateral1.80.8208.915.334.043.71
Octagonal1.80.8209.555.974.684.35
Hexadecagonal1.80.8203.790.21−1.08−1.41
Quadrilateral1.70.9208.815.233.943.61
Octagonal1.70.9209.285.74.414.08
Hexadecagonal1.70.9203.60.02−1.27−1.6
Quadrilateral1.60.6208.95.324.033.7
Octagonal1.60.6209.375.794.54.17
Hexadecagonal1.60.6203.680.1−1.19−1.52
Quadrilateral1.80.8259.136.214.461.76
Octagonal1.80.8259.646.624.972.27
Hexadecagonal1.80.8254.191.17−0.48−3.18
Quadrilateral1.70.9258.855.934.181.48
Octagonal1.70.9259.56.584.832.13
Hexadecagonal1.70.9254.181.26−0.49−3.19
Quadrilateral1.60.6258.956.034.281.58
Octagonal1.60.6259.526.64.852.15
Hexadecagonal1.60.6254.251.33−0.42−3.12
Quadrilateral1.80.8309.16.594.682.47
Octagonal1.80.8309.777.265.353.14
Hexadecagonal1.80.8304.692.180.27−1.94
Quadrilateral1.70.9308.796.284.372.16
Octagonal1.70.9309.597.085.172.96
Hexadecagonal1.70.9304.411.9−0.01−2.22
Quadrilateral1.60.6308.836.324.412.2
Octagonal1.60.63010.57.996.083.87
Hexadecagonal1.60.6304.672.160.25−1.96
Table 3. PSNRs for chest X-ray images denoised using the proposed 2D LPA-RICI method and the competitive methods.
Table 3. PSNRs for chest X-ray images denoised using the proposed 2D LPA-RICI method and the competitive methods.
Parameters PSNR (dB)
RegionГRcσNoisy ImageRICIFixed Size FiltersGaussian SmoothingTotal Variation
Quadrilateral1.80.82029.3636.833.2334.1835.34
Octagonal1.80.82029.3638.0933.2334.1835.34
Hexadecagonal1.80.82029.3633.7633.2334.1835.34
Quadrilateral1.70.92029.3636.7433.2334.1835.34
Octagonal1.70.92029.3638.2333.2334.1835.34
Hexadecagonal1.70.92029.3633.8133.2334.1835.34
Quadrilateral1.60.62029.3636.633.2334.1835.34
Octagonal1.60.62029.3638.0633.2334.1835.34
Hexadecagonal1.60.62029.3633.333.2334.1835.34
Quadrilateral1.80.82528.2637.2831.8633.7636.43
Octagonal1.80.82528.2638.4931.8633.7636.43
Hexadecagonal1.80.82528.2633.5131.8633.7636.43
Quadrilateral1.70.92528.2637.6631.8633.7636.43
Octagonal1.70.92528.2638.4731.8633.7636.43
Hexadecagonal1.70.92528.2632.5831.8633.7636.43
Quadrilateral1.60.62528.2637.2631.8633.7636.43
Octagonal1.60.62528.2638.4831.8633.7636.43
Hexadecagonal1.60.62528.2633.0931.8633.7636.43
Quadrilateral1.80.83028.8737.7431.0133.1735.22
Octagonal1.80.83028.8738.4731.0133.1735.22
Hexadecagonal1.80.83028.8733.0131.0133.1735.22
Quadrilateral1.70.93028.8736.9631.0133.1735.22
Octagonal1.70.93028.8738.5231.0133.1735.22
Hexadecagonal1.70.93028.8733.1431.0133.1735.22
Quadrilateral1.60.63028.8737.6231.0133.1735.22
Octagonal1.60.63028.8737.9431.0133.1735.22
Hexadecagonal1.60.63028.8732.9231.0133.1735.22
Table 4. PSNR improvements for chest X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
Table 4. PSNR improvements for chest X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
ParametersPSNR (dB)
RegionГRcσRICI vs. NoisyRICI vs. FixedRICI vs. GaussianRICI vs. Total Variation
Quadrilateral1.80.8207.443.572.621.46
Octagonal1.80.8208.734.863.912.75
Hexadecagonal 1.80.8204.40.53−0.42−1.58
Quadrilateral1.70.9207.383.512.561.4
Octagonal1.70.9208.8754.052.89
Hexadecagonal 1.70.9204.450.58−0.37−1.53
Quadrilateral1.60.6207.243.372.421.26
Octagonal1.60.6208.74.833.882.72
Hexadecagonal 1.60.6203.940.07−0.88−2.04
Quadrilateral1.80.8259.025.423.520.85
Octagonal1.80.82510.236.634.732.06
Hexadecagonal 1.80.8255.251.65−0.25−2.92
Quadrilateral1.70.9259.45.83.91.23
Octagonal1.70.92510.216.614.712.04
Hexadecagonal 1.70.9254.320.72−1.18−3.85
Quadrilateral1.60.62595.43.50.83
Octagonal1.60.62510.226.624.722.05
Hexadecagonal 1.60.6254.831.23−0.67−3.34
Quadrilateral1.80.8308.876.734.572.52
Octagonal1.80.8309.67.465.33.25
Hexadecagonal 1.80.8304.142−0.16−2.21
Quadrilateral1.70.9308.095.953.791.74
Octagonal1.70.9309.657.515.353.3
Hexadecagonal 1.70.9304.272.13−0.03−2.08
Quadrilateral1.60.6308.756.614.452.4
Octagonal1.60.6309.076.934.772.72
Hexadecagonal 1.60.6304.051.91−0.25−2.3
Table 5. PSNRs for ankle X-ray images denoised using the proposed 2D LPA-RICI method and the competitive methods.
Table 5. PSNRs for ankle X-ray images denoised using the proposed 2D LPA-RICI method and the competitive methods.
Parameters PSNR (dB)
RegionГRcσNoisy ImageRICIFixed Size FiltersGaussian SmoothingTotal Variation
Quadrilateral1.80.82029.8536.6832.7932.2234.45
Octagonal1.80.82029.8536.1632.7932.2234.45
Hexadecagonal 1.80.82029.8532.9232.7932.2234.45
Quadrilateral1.70.92029.8536.9732.7932.2234.45
Octagonal1.70.92029.8536.2832.7932.2234.45
Hexadecagonal 1.70.92029.8532.8332.7932.2234.45
Quadrilateral1.60.62029.8536.9332.7932.2234.45
Octagonal1.60.62029.8536.4132.7932.2234.45
Hexadecagonal 1.60.62029.8532.9532.7932.2234.45
Quadrilateral1.80.82529.3737.2331.9631.9335.61
Octagonal1.80.82529.3736.5731.9631.9335.61
Hexadecagonal 1.80.82529.3732.9431.9631.9335.61
Quadrilateral1.70.92529.3737.0731.9631.9335.61
Octagonal1.70.92529.3736.6731.9631.9335.61
Hexadecagonal 1.70.92529.3732.9831.9631.9335.61
Quadrilateral1.60.62529.3737.0331.9631.9335.61
Octagonal1.60.62529.3736.8331.9631.9335.61
Hexadecagonal 1.60.62529.3732.9831.9631.9335.61
Quadrilateral1.80.83029.0637.1531.1531.6934.35
Octagonal1.80.83029.0636.3531.1531.6934.35
Hexadecagonal 1.80.83029.0632.6331.1531.6934.35
Quadrilateral1.70.93029.0636.9731.1531.6934.35
Octagonal1.70.93029.0636.2831.1531.6934.35
Hexadecagonal 1.70.93029.0632.8331.1531.6934.35
Quadrilateral1.60.63029.0637.1131.1531.6934.35
Octagonal1.60.63029.0636.5131.1531.6934.35
Hexadecagonal 1.60.63029.0632.9931.1531.6934.35
Table 6. PSNR improvements for ankle X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
Table 6. PSNR improvements for ankle X-ray images denoised using the proposed 2D LPA-RICI method vs. the competitive methods.
Parameters PSNR (dB)
RegionГRcσRICI vs. NoisyRICI vs. FixedRICI vs. GaussianRICI vs. Total Variation
Quadrilateral1.80.8206.833.894.462.23
Octagonal1.80.8206.313.373.941.71
Hexadecagonal 1.80.8203.070.130.70−1.53
Quadrilateral1.70.9207.124.184.752.52
Octagonal1.70.9206.433.494.061.83
Hexadecagonal 1.70.9202.980.040.61−1.62
Quadrilateral1.60.6207.084.144.712.48
Octagonal1.60.6206.563.624.191.96
Hexadecagonal 1.60.6203.100.160.73−1.50
Quadrilateral1.80.8257.865.275.301.62
Octagonal1.80.8257.204.614.640.96
Hexadecagonal 1.80.8253.570.981.01−2.67
Quadrilateral1.70.9257.705.115.141.46
Octagonal1.70.9257.304.714.741.06
Hexadecagonal 1.70.9253.611.021.05−2.63
Quadrilateral1.60.6257.665.075.101.42
Octagonal1.60.6257.464.874.901.22
Hexadecagonal 1.60.6253.611.021.05−2.63
Quadrilateral1.80.8308.096.005.462.80
Octagonal1.80.8307.295.204.662.00
Hexadecagonal 1.80.8303.571.480.94−1.72
Quadrilateral1.70.9307.915.825.282.62
Octagonal1.70.9307.225.134.591.93
Hexadecagonal 1.70.9303.771.681.14−1.52
Quadrilateral1.60.6308.055.965.422.76
Octagonal1.60.6307.455.364.822.16
Hexadecagonal 1.60.6303.931.841.30−1.36

Share and Cite

MDPI and ACS Style

Mandić, I.; Peić, H.; Lerga, J.; Štajduhar, I. Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm. J. Imaging 2018, 4, 34. https://doi.org/10.3390/jimaging4020034

AMA Style

Mandić I, Peić H, Lerga J, Štajduhar I. Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm. Journal of Imaging. 2018; 4(2):34. https://doi.org/10.3390/jimaging4020034

Chicago/Turabian Style

Mandić, Ivica, Hajdi Peić, Jonatan Lerga, and Ivan Štajduhar. 2018. "Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm" Journal of Imaging 4, no. 2: 34. https://doi.org/10.3390/jimaging4020034

APA Style

Mandić, I., Peić, H., Lerga, J., & Štajduhar, I. (2018). Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm. Journal of Imaging, 4(2), 34. https://doi.org/10.3390/jimaging4020034

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