Next Article in Journal
Factors Predictive for Immunomodulatory Therapy Response and Survival in Patients with Hypersensitivity Pneumonitis—Retrospective Cohort Analysis
Next Article in Special Issue
Machine Learning for Online Automatic Prediction of Common Disease Attributes Using Never-Ending Image Learner
Previous Article in Journal
IDH Mutations Are Potentially the Intrinsic Genetic Link among the Multiple Neoplastic Lesions in Ollier Disease and Maffucci Syndrome: A Clinicopathologic Analysis from a Single Institute in Shanghai, China
Previous Article in Special Issue
Current Status of 68Ga-Pentixafor in Solid Tumours
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Noise Suppression and Edge Preservation for Low-Dose COVID-19 CT Images Using NLM and Method Noise Thresholding in Shearlet Domain

1
Computer Science and Engineering Department, Graphic Era Deemed to be University, Dehradun 248007, India
2
School of Computer Science Engineering and Technology, Bennett University, Greater Noida 201310, India
3
Department of Basic Science, College of Science and Theoretical Studies, Saudi Electronic University, Riyadh-Male Campus, Riyadh 13316, Saudi Arabia
4
Department of Computer Science and Engineering, Amity School of Engineering and Technology, Amity University, Noida 201303, India
5
Center for Artificial Intelligence, Prince Mohammad Bin Fahd University, Khobar 34754, Saudi Arabia
6
Department of Computer Science and Information Engineering, National Cheng Kung University, Tainan 701, Taiwan
7
School of Computer Science, University of Petroleum and Energy Studies, Dehradun 248007, India
*
Author to whom correspondence should be addressed.
Diagnostics 2022, 12(11), 2766; https://doi.org/10.3390/diagnostics12112766
Submission received: 28 September 2022 / Revised: 9 November 2022 / Accepted: 9 November 2022 / Published: 12 November 2022

Abstract

:
In the COVID-19 era, it may be possible to detect COVID-19 by detecting lesions in scans, i.e., ground-glass opacity, consolidation, nodules, reticulation, or thickened interlobular septa, and lesion distribution, but it becomes difficult at the early stages due to embryonic lesion growth and the restricted use of high dose X-ray detection. Therefore, it may be possible for a patient who may or may not be infected with coronavirus to consider using high-dose X-rays, but it may cause more risks. Conclusively, using low-dose X-rays to produce CT scans and then adding a rigorous denoising algorithm to the scans is the best way to protect patients from side effects or a high dose X-ray when diagnosing coronavirus involvement early. Hence, this paper proposed a denoising scheme using an NLM filter and method noise thresholding concept in the shearlet domain for noisy COVID CT images. Low-dose COVID CT images can be further utilized. The results and comparative analysis showed that, in most cases, the proposed method gives better outcomes than existing ones.

1. Introduction

The novel coronavirus disease 2019 (COVID-19) emerged in early December 2019 in the Hubei province of the People’s Republic of China, and it was caused by the virus Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). It rapidly spread to various countries and developed into a global outbreak. The World Health Organization (WHO) declared the outbreak as a global pandemic on 11 March 2020 [1,2]. This also led to a massive strain on healthcare systems worldwide.
While CT scan imagery might be effective in coronavirus detection, it has certain limitations when it comes to continuous use and its subsequent hazardous side effects. The detection of COVID-19 can be performed by identifying lesions in scans, namely, ground-glass opacity, consolidation, nodules, reticulation, or thickened interlobular septa, and lesion distribution [3,4], but in the early stages, due to the nascent development of lesions and limited use of high dose X-rays, identification is difficult. These symptoms, when not detected early, multiply the spread not only for the victim, but also for others around the victim. The dilemma here is that, for early diagnoses of infection, doctors require either the continuous use of X-rays or the use of high-dose X-rays, both of which lead to hazardous side effects such as vomiting, bleeding [5,6,7,8,9], fainting, hair loss, and the loss of skin and hair and cancerous cells [9,10,11,12,13,14,15]. This leaves us to use low-dose X-rays for CT scans, which give us scans corroded by noise and disturbance. Thus, deliberating the above, the use of high-dose X-rays may create further complications for a patient who may or may not be infected with coronavirus [16,17,18]. Conclusively, the only way to prevent patients from side effects consistent with high-dose X-rays while diagnosing the presence of coronavirus early is using low-dose X-rays to generate CT scans, then implementing a robust denoising algorithm to the scans [19,20,21,22,23]. As explained above, the importance and need for denoising have developed more than ever. Moreover, image segmentation for the classification and identification of the type of lesion is also a substantial task in the diagnosis procedure. There are different methods for it, and various researchers have contributed to the cause [24,25,26,27,28,29].
Medical image denoising [30,31,32,33,34,35,36] is the process of removing corrosive and corrupting elements from the image in order to obtain an image of superior clarity and quality. While there are multiple techniques used to denoise images, they are graded into two major purviews, namely spatial techniques such as stationary multiplicative speckle model filters (SMSM), non-stationary multiplicative speckle model filters (NSMSM), among others, and transform techniques such as threshold-based techniques, Bayesian techniques, directional transform techniques, non-local means (NLM), block-matching 3D (BM3D), and others, each having variegated techniques with respective advantages and disadvantages [37]. Another robust technique in the field of image processing for handling tasks such as image denoising and image segmentation in shearlet transform consists of a multi-scale framework for the multi-dimensional representation of images. Recently, shearlet transform is being implemented in variegated domains; for example, the study [38] used optimized threshold shearlet transform for denoising microscopic data by eliminating random noise. Thresholding is performed for shearlet coefficients by shrinking the same the microscopic data random noise that is distributed unevenly. The novelty of the method lies in the use of adaptive thresholding coefficients, which are obtained by the use of the adjustment factor for each fundamental sub-band individually [39].
Another state-of-the-art denoising technique implements a shearlet transform that amalgamates the same, with an oriented second-order partial equation, to remove multiplicative speckle noise in electronic speckle pattern interferometry (ESPI) fringe patterns. The novel proposition of this study lies in the use of shearlet transform in ESPI fringe patterns with the aim to denoise. It leverages the concept that the spatial-domain filtering method SOOPDE and shearlet transform are symbiotic to each other on the grounds of mutual benefit [40]. This study leveraged the unique ability of shearlet transform to handle multi-modal data at a multi-scale level such as a wide range of orientations, complex geometric figures, etc. It aimed to gauge the performance of shearlet transform for denoising and segmentation on the grounds of space-time complexity for a set of GPUs by implementing 2D and 3D on standard GPU and proposed accelerated GPU. Preserving edges and the texture of images when denoising is a challenging task. Non-subsampled shearlet transform (NSST) can perform a tenable job, but amalgamating with a twin support vector machine can provide better results, as conducted in the presented study [41]. Primarily, it divides the noisy image into multiple sub-bands oriented to NSST, later using a twin support vector for feature detection and denoising [42].
For a similar study that implemented NSST for image denoising, its novelty lied in amalgamating sparse representation classification with NSST. This study followed a multi-step procedure. NSST coefficients are primarily classified into noise or edge coefficients using spars representation via variable splitting and augmented Lagrangian (SUnSAL) classifiers, posting which noise coefficients are denoised by the shrink method with adaptive Bayesian threshold. Finally, inverse NSST reverses the coefficients to obtain the denoised image [43,44,45,46,47,48,49,50]. A similar technique that implements shearlet transform via signal detecting operation by using a state-of-the-art proposed strategy that uses the correlation with energy distribution of tri-component microseismic signals shows substantial performance both on stimulated and real noisy data sets [51,52]. Optimizing shearlet transformation with fellow wavelet techniques is common, but it can be performed by manipulating the coefficients shown by the presented work [53]. This study shows a groundbreaking increase in the performance of shearlet transform by manipulating its constraints of residual coefficients.
This study proposes a novel low-dose CT image denoising technique by amalgamating two state-of-the-art methodologies, the NLM filter and Bivariant thresholding, that provides superior results to the classic NLM filter, which suppresses sharp noise and thresholding, which works on edge preservation. The presented research framework evaluated and analyzed the results obtained in imageries via qualitative and quantitative analysis and graphical and statistical means. The rest of the paper is organized as follows: A brief overview of methods such as NLM and shearlet transform are given in Section 2. In Section 3, the proposed methodology is described. In Section 4, results are analyzed and discussed. Finally, the conclusion is drawn in Section 5.

2. Preliminaries

In this section, some major concepts are discussed that were utilized in the proposed methodology.

2.1. NSST or Non-Subsampled Shearlet Transform

Through the use of decomposition and the preservation of translation invariance, NSST is able to accomplish multi-scale and multi-directional transformations while avoiding the Gibbs phenomenon. The pseudo-grid is processed using a two-dimensional Fourier transform before being passed through a 1D sub-band filter. As a result, essential sampling can be eliminated. Due to its ability to adapt to the diverse geometrical properties of multi-dimensional and multi-scale data, NSST has an edge over other transform methods.
N S L P j + 1 = A j f = A h j 1 k = 1 j 1 A h k 0 f
Primarily, the low-pass image and the high-pass image are obtained by decomposing the original image. The coefficient for the pseudo grid is obtained by calculation, which is then used to obtain frequency-domain coefficients by bandpass filtering. Finally, non-subsampled shearlet transform coefficients are obtained by implementing an inverse transformation of the fast Fourier transform (FFT) on the pseudo grid.
For p > 0, q ∈ R, t ∈ R2, the shearlets can be expressed as:
Ψ p . q . f = | d e f F p . q | 1 Δ Ψ F p q 1 x f
where   F p , q = B q A p = ( p p q 0 p ) , A = p 0 0 p , B q = 1 q 0 1 .
Only the shear matrix controls the shearlet’s direction. As a result, the shearlet transform depends on the scale (p), the orientation (q), and the position. Each f ∈ L2 (R2) can be recovered by:
f = R 2 + 0 + < f ,   Ψ p . q . f > Ψ p . q . f d p p 3 d q d t
Because of its mathematical basis, the discrete shearlet transform can be effectively represented in a variety of contexts. The discrete implementation of the shearlets transform offers the ability to take into account functions with multiple dimensions. With this consideration, p = 2 – 2J, q = −L with J, t = K ∈ Z2, and L ∈ Z. The discrete shearlet transform is expressed as:
Ψ p . q . f = | det   A 0 | J 2 Ψ B 0 L A 0 L x K
where A 0 = 4 0 0 2 , and B 0 = 1 1 0 1 .
For each L 2   R 2 , the formula can be reproduced with the characteristic Ψ , as shown below:
f = J , L     Z , K     Z 2   f ,   ψ J L . K ψ J L . K

2.2. Non-Local Mean (NLM)

NLM effectively preserves the variegated regions of multiple types of textures while cleaning the input image. While other filters deliberate the mean of neighboring pixels of the target image, i.e., pixels surround target pixels, NLM deliberates the mean of all pixels in the target image. However, the weights depend upon the similarity of pixels, or simply put, the higher the similarities higher the weights are.
N L v i = j I w i , j v j
This approach results in a lesser degree of distortion of the denoised pixel from the original when compared to other methods. The mean of all pixels across is taken as a whole using a Gaussian window, and the value of a particular pixel is evaluated based on its position in the image, up to a range of values determined by the mean pixel value of said image.

3. Proposed Methodology

With the assumption that CT images are noisy due to additive white Gaussian noise, a new method was proposed to remove the noise from CT images. Figure 1 shows the complete overview of the proposed framework. This noise predominantly attacks mainly on high frequencies, although its effect on low frequencies is non-trivial. In the proposed method, the noise removal process is included so that noise is suppressed effectively and the edges are preserved at the same time. In the proposed methodology, the characteristics of the NLM filter are utilized, and a novel thresholding function was proposed for noise suppression and edge preservation.

Proposed Thresholding Function

Let the CT image be noisy, which can be denoted as:
X(p,q) = Y(p,q) + n(p,q)
where X(p,q) is a noisy CT image, Y is a clean CT image, and n is an additive Gaussian noise.
This is used in the filtration of high-frequency bands because it is a derivative of bilateral thresholding. Taking the position of a coefficient of a wavelet function and substituting it with its parent of a higher scale forms the background of this function. Assuming that w1k is the child of w2k, where w1k is the kth complex wavelet coefficient at the same position as that of w1k, then we can express them as:
y 1 k = w 1 k + η 1 k y 2 k = w 2 k + η 2 k
where y1k and y2k are noisy complex wavelet coefficients, η1k and η2k are additive noise coefficients, and wk = (w1k, w2k), yk = (y1k, y2k), and ηk = (η1k, η2k). From this equation, after a bit of mathematical manipulation, the bivariate shrinkage function or thresholding can be expressed as:
w ˆ 1 k =   y 1 k 2 +   y 2 k 2     λ k y 1 k 2 +   y 2 k 2   y 1 k
The function is defined as:
g a * = a                                     i f   a > 0 0                               o t h e r w i s e
In the suggested thresholding equation, the multi-variate function is modified by the definition of process noise with the SURE-LET technique. Therefore, it can be expressed as:
z c = y c w ^ c
θ y c = d 1 w ^ c + d 2 z c
where the weight values that should be optimized to achieve the best outcomes are d1 and d2.
It can be expressed as:
θ y c = d 1 d 2 w ^ c + d 2 y c
Let a 1 = d 2 ,   a 2 = d 1 d 2 ,   W c 1 = w ^ c , and W c 2 = y c , and then the above equation can be expressed as:
θ y c = a 1 W k 1 + a 2 W k 2
where the parameter a c   1 c 2 is linear.
It is possible to re-design the equation above as:
θ y c = s = 1 2 a s φ s y c
A thresholding function is proposed for the above linear function, as follows:
θ y new _ c = 0 , if   θ y c < 0 θ y c   , Otherwise
A composition of different estimates of noise-unfastened coefficients is the suggested thresholding function. For the best denoising outcomes, the best ak weight values must be optimized. An unbiased mean-square-error (MSE) estimator is known to optimize weight values, followed by the SURE-LET method. It can be conveyed as:
ε = 1 M N c = 1 M N θ 2 y c 2 y c θ y c + 2 σ 2 θ y c + s = 1 M N w c 2
To estimate the parameter a , only ε = 1 M N ( θ 2 y 2 y θ y ) MSE for a s can be expressed as:
f = 1 4 P s . f a ζ Q s = 0
where P s . f = φ s y φ ζ y and Q s = y φ s y σ 2 φ s y . The above equation can be solved by the linear method: a ζ = P s . f 1 Q s
Figure 1 successfully illustrates the denoising scheme followed in the presented framework. It is evident that low-dose CT images have different levels of accuracy due to noise intensity and pattern that vary considerably; thus, two images were taken to generate a more accurate result. It must be noted that, although two images may require two CT scans, the total radiation of the two combined would still be considerably less than the radiation exposure from a normal CT scan.
The summary of the proposed work is shown in Algorithm 1:
Algorithm 1: CT image denoising.
Step 1: Firstly, read input noisy CT image.
Step 2: Non-subsampled shearlet transform is applied to both noisy images, which divides the image into two parts:
 a. Approximation part (H)
 b. Detailed part (D)
Step 3: Apply k and l directional circular shift to obtain n high-frequency sub-bands of both input images.
       H k , l s y = circular _ shift   H k , l s , k shift , l shift
Step 4: Perform NLM filter on both approximation part.
Step 5: Perform average operation on the outcomes of Step 4.
Step 6: For all levels in high-frequency sub-bands of both input images:
(a) Calculate the threshold value
(b) Apply shrinkage rule using Equation (16)
Step 7: To obtain an enhanced high-frequency sub-band, calculate the weighted average based on patch variance on the outcome of Step 6:
         H ^ k , l s = y = 1 n α s y . H k , l s y
where α s y = var 1 X w , s y n y = 1 n var 1 H w , s y n , var−1 (.) represents the inverse of threshold, and H(k,l)s is the final threshold value.
Step 8: To obtain the final output image, perform the inverse of the circular shift using the outcome of Step 5 and Step 7:
             H k , l sy = circular _ shift H ^ k , l s y , k shift , l shift

4. Results and Discussion

The experiments were carried out in Matlab 2018 running on a system comprising an i5 8250 with 8 GB of RAM running 64-bit Windows 10 as the operating system. All images used were of dimensions 512 × 512 pixels, as illustrated in Figure 2, which were taken from the Cornell open-access library dataset.
In the proposed method, a thresholding function was performed in detail parts of NSST domain. In the approximation part of NSST, the NLM filter was performed. The denoising process was performed using patch-wise circular shifting. To show the significance of the proposed method, performance metrics (PSNR, SSIM, DIV, and ED) were estimated without circular shifting and with circular shifting using different patches. The average values (87 images) of the performance metrics are shown in Table 1. Here, it can be clearly analyzed that the best values were given by the proposed method, with circular shifting using 5 × 5 patches. Hence, in our proposed method, 5 × 5 patches were used to evaluate the result analysis.

4.1. Comparative Analysis

To maintain consistency, all the Gaussian noise present in the images was artificially applied. Gaussian noise was added to all the images at different intensities, with sigma values between 30 and 5. The noisy images are shown in Figure 3. Figure 4, Figure 5, Figure 6 and Figure 7 show the results of the denoising process, as performed by the different competing algorithms that were compared. Figure 5 shows the denoised output from all eight different frameworks when applied on the first image from Figure 1, respectively. Similarly, Figure 6, Figure 7 and Figure 8 show the results of the denoising process for the case of images 2, 3, and 4 from Figure 2 as outputted by the frameworks [10,11,12,13,14,15,16], respectively.
As shown in Figure 4, Figure 4a, corresponding to [10], shows the least feasible result of the group. The edges were not preserved, the contrast was altered, and rogue pixels were introduced, creating noise. Although the result was better than the noisy image, the loss of detail and the excessive fuzziness take away from its ability to be useful in real-world applications. Figure 4b, corresponding to [11], shows slightly better results, although the overall quality was still much lower than the quality of the other techniques. This framework introduces less noise but even more blur than [10]. The details of [12] and [13] are very incorrigible, with visually noticeable noise generated from the region with intermediate grey values.
In Figure 5, the frameworks deal with images that need more contrast preservation and heterogeneous integrity. Figure 5a, corresponding to [10], gave the worst output, with no noticeable improvement from the noisy image. Figure 4b,c shows the results of [12] and [13]. The best results demonstrated by Figure 5g,h correspond to [16] and the proposed method. No doubt, the proposed system came out on top; however, [16] came pretty close in terms of contrast preservation. Nothing special or noteworthy was illustrated, showing a generally average performance with a slight loss in the intensity of black pixels in all of [11,14,15].
While all the frameworks performed very close to each other for this image, [10] showed a high level of fuzziness or blur introduced as denoising artifacts by the overcorrection of the image. The [11] framework showed good edge preservation, with the contrasts almost accurate to the original, but the presence of increased noise in the background degraded the overall quality of the image. Similarly, [12] sufficiently removed noise with proper clarity and no added noise; the detail preservation was actually inferior to [11], especially in heterogeneous noisy areas. The [13] framework changed the background and increased the overall pixel value of the image by giving pixels that were black a lower color value. Amongst all the methods tested, [14] had the most inferior edge preservation results, introducing noise at the edges of the images, particularly when there was a high contrast between the subject and the background, which further corroded the image quality. Both [15,16] performed very well, with almost negligible performance between each other and the proposed methodology.
The outcome comparison of the fourth and final image is shown in Figure 7. The quantity of blur pertaining to the white noise introduced in [10] made the image as bad as the noisy image. Almost no improvement is seen, which implies that scans with regions dense in muscle will absolutely not work with [10]. Ref. [11], shown in Figure 7b, had better contrast and slightly better edge detection, although it was still way off from the original or even its competing denoised counterparts. Figure 7c, corresponding to [13], performed a satisfactory job, but it can be observed that the contrast was reduced all throughout the image, with the bones getting darker and the background getting lighter. Figure 7e,f shows nothing noticeably special. Satisfactory edge preservation, contrast maintenance, and noise suppression were conducted as expected. Figure 7g, the outcome from [15], was deemed as the best among the older techniques, coming closest to the original after the proposed technique. The only difference between it and the proposed framework was the combination of increased sharpness and slightly better contrast that is clear from the pattern of the spine in the two images. Figure 7g is less sharp, and the individual vertebrae are harder to distinguish, while Figure 7h is both whiter and easier to discern.
Figure 8 shows a line used for intensity profiling in CT image for the analysis of all the frameworks. Figure 9 shows the intensity profiles of the original image against the [10,11,12,13,14,15,16] frameworks and the proposed framework, respectively. From Figure 8 and Figure 9, certain observations about the frameworks were made. While all the frameworks performed reasonably well, [10] showed a high level of fuzziness or blur introduced as denoising artifacts by the overcorrection of the image. The [11] framework showed good edge preservation, with the contrasts almost accurate to the original, but the presence of increased noise in the background degraded the overall quality of the image. Similarly, [12,13] performed sufficient removal of the noise with proper clarity and no added noise; the detail preservation was actually inferior to [11], especially in heterogeneous noisy areas. The [14] framework had the most inferior edge preservation results amongst all the methods tested, introducing noise at the edges of the images, particularly when there was a high contrast between the subject and the background, which further corroded the image quality. Both [15,16] performed very well, with almost negligible performance between each other and the proposed methodology.
The proposed framework was found to be supplanting all the other frameworks, scoring better in all aspects, including but not limited to edge preservation, noise suppression for both heterogeneous and contrasting areas of the input image, detail preservation, and overall similarity to the original image. The heterogeneous noise removal was good, with equally good contrast preservation to go with it. Thus, it can be confidently said that the proposed framework competed with and outperformed the other state-of-the-art frameworks in almost all aspects of speckle image denoising in medical CT scan images.

4.2. Performance Metrics

SSIM: The structural similarity index is a full reference metric used to quantify image degradation by comparing two windows, X and Y, to the original and processed image. This method works in a manner similar to humans and, hence, can be taken as a crude estimate close to what a human would grade the image. The original image is required for the use of this metric, the formula for which is given below:
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 μx, μy are the mean values of noisy and denoised images, respectively. C1 and C2 are the constant values; σ is the variance of the respective image.
PSNR: Peak signal-to-noise ratio, as is evident from the name, is the ratio between the total power, or original image in this case, to the power of the noise responsible for declining the image quality. It is also a full reference metric, and the formula is:
P S N R = 10 log 10 I M S M S E
where IMS is the image size in m*n pixels, and MSE is the mean square error calculated as: M S E = 1 N j = 0 N 1 X Y 2 .
Entropy Difference (ED): Entropy is the randomness or chaos in a system. In the image, it is the chaos present, i.e., the difference from the original source. The entropy difference is the mean of the Sharon entropy (SE) between the clean and denoised image.
E D = S E X S E Y
Difference in Variance (DIV): Performance measure by variance (Var) used in statistics.
D I V = 1 V A R X V A R Y
Table 2 shows the average of all the PSNR and SSIM values of all the images that were tested. The testing performed on a sample set of 90 images was enough to normalize and negate the impact of any bias that a technique might have for one particular image, whilst increasing the legitimacy of this technique. The Σ symbol corresponds to the value of sigma, which defines the amount of noise introduced in the images. The higher the value, the greater the noise introduced in the image and the lower the quality of the output image even after denoising. In this table, the first column shows the different frameworks, and the next four columns show the values of PSNR corresponding to the particular framework at a particular noise level. This is why the values kept decreasing as the value of Σ kept increasing. In almost all the cases, the values generated by the different frameworks were close to each other, with differences noticeable. Except for Σ = 30, the proposed framework was superior to all other frameworks by a noticeable margin of 0.4–0.9. This is greater than the average difference between results amongst the other frameworks, proving the efficacy of this technique. Overall, by comparing on the basis of SSIM and PSNR, it was observed that [15] was the best among the existing techniques, and although the proposed technique is better than [15], the difference was small. The [11] framework showed the most consistency, with a graceful decline in scores with the increase of noise levels. The other techniques, [10,12,13,14,16], performed reasonably without anything notable worth mentioning.
Similar to Table 2, Table 3 shows the ED and DIV values of the denoised images generated from the frameworks at different noise levels. Unlike the previous techniques, these values increased with the increase in noise and subsequently for the dissimilarity between images. In this case, [14] was the most volatile and shows the largest variation with a large range and huge jumps in scores with the increase in noise. The [16] framework scored relatively well, coming very close to the proposed framework at lower noises but rapidly lost its performance with an increase in noise. The [10,11] frameworks had the lowest scores overall, which is in support of the observations made in the qualitative analysis. The [13] framework had an average set of scores, neither too low to be comparable to the proposed framework nor too high to be the worst of the set. The [12] framework’s ED value came closest to the proposed framework at the lowest noise level but showed no significant score in the rows pertaining to the higher noise values. Overall, the proposed framework beat all the systems, similar to its superior performance seen in Table 1. The DIV score of 0.02 at Σ = 10 is an order of magnitude better than the other competitions.
A total of 112 experts pertaining to the domain of image processing were consulted for their opinion on the quality of the image and, therein, their judgment on the performance of the algorithms. Table 4 shows the average of the scores granted by the experts to the various techniques in various aspects or features used in tandem with image denoising quality measurement. The scores range from 1 to 3, with 3 denoting the best performance. As seen in the table, the proposed framework garnered high approval in most of the fields, according to the experts. This is a clear demonstration that adds weight to the validity of the framework.

5. Conclusions

This study presented a novel hybrid denoising framework for low-dose CT images that combined non-subsampled shearlet transform and bivariate thresholding to obtain pristine results that surpassed the state-of-the-art techniques on the grounds of edge preservation, feature and detail preservation, denoising contrast, and heterogeneous area of the image, and it prevented the creation of artifacts during the process of denoising. The results obtained were validated by comparing the proposed framework with the latest techniques in the domain via qualitative analysis and quantitative analysis. Reference image indexes, namely, PSNR, SSIM, entropy difference, and variance difference of the proposed and state-of-the-art-methods, were compared to support the quantitative analysis. Similarly, the qualitative analysis was supported by the comparison of image intensity profiles between the original and denoised images. The presented study addressed the dilemma of using high-dose X-rays that provide noise-free images but have multiple deleterious effects versus using low-dose CT images that are corrupted by noise but have comparatively very few or negligible side effects. The presented study provided noise-free CT images that aid COVID-19 diagnoses without creating any salutary hazards to the patient due to X-rays, thereby being applied widely in clinical medical image capturing, classification, and diagnosis of coronavirus.

Author Contributions

Conceptualization, M.D.; Data curation, M.D., P.S. and M.J.; Formal analysis, M.D., E.B. and V.R.; Funding acquisition, C.S. and V.R.; Investigation, P.S. and V.R.; Methodology, M.D., P.S. and E.B.; Project administration, M.J. and V.R.; Resources, E.B., M.J. and V.R.; Software, P.S. and E.B.; Supervision, V.R.; Visualization, M.J.; Writing–original draft, M.D.; Writing–review & editing, M.J., K.U.S. and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shi, F.; Wang, J.; Shi, J.; Wu, Z.; Wang, Q.; Tang, Z.; He, K.; Shi, Y.; Shen, D. Review of artificial intelligence techniques in imaging data acquisition, segmentation and diagnosis for covid-19. IEEE Rev. Biomed. Eng. 2021, 14, 4–15. [Google Scholar] [CrossRef] [Green Version]
  2. Dong, D.; Tang, Z.; Wang, S.; Hui, H.; Gong, L.; Lu, Y.; Xue, Z.; Liao, H.; Chen, F.; Yang, F.; et al. The role of imaging in the detection and management of COVID-19: A review. IEEE Rev. Biomed. Eng. 2021, 14, 16–29. [Google Scholar] [CrossRef]
  3. Ulhaq, A.; Khan, A.; Gomes, D.; Paul, M. Computer Vision for COVID-19 Control: A Survey. arXiv 2020, arXiv:2004.09420. [Google Scholar] [CrossRef]
  4. Li, K.; Wu, X.; Zhong, Y.; Qin, W.; Zhang, Z. Diagnostic performance of CT and its key signs for COVID-19: A systematic review and meta-analysis. medRxiv 2020. [Google Scholar] [CrossRef]
  5. Nair, A.; Rodrigues, J.C.L.; Hare, S.; Edey, A.; Devaraj, A.; Jacob, J.; Robinson, G. A British Society of Thoracic Imaging statement: Considerations in designing local imaging diagnostic algorithms for the COVID-19 pandemic. Clin. Radiol. 2020, 75, 329–334. [Google Scholar] [CrossRef]
  6. Huang, Q.; Li, J.; Lyu, S.; Liang, W.; Yang, R.; Zhang, R.; Chen, F. COVID-19 associated kidney impairment in adult: Qualitative and quantitative analyses with non-enhanced CT on admission. Eur. J. Radiol. 2020, 2020, 109240. [Google Scholar] [CrossRef]
  7. Hamilton, N.E.; Adam, G.H.; Ifan, D.L.; Lam, S.S.; Johnson, K.; Vedwan, K.A.G.; Abbas, A. Diagnostic utility of additional whole-chest CT as part of an acute abdominal pain CT imaging pathway during the COVID-19 pandemic. Clin. Radiol. 2020, 75, 592–598. [Google Scholar] [CrossRef]
  8. Jiang, M.; Chen, P.; Li, T.; Tang, Y.; Chen, X.; Chen, X.; Ruan, X. Chest CT imaging features and clinical outcome of coronavirus disease 2019 (COVID-19): A single-center case study in Ningbo, China. Clin. Imaging 2021, 69, 27–32. [Google Scholar] [CrossRef]
  9. Bao, C.; Liu, X.; Zhang, H.; Li, Y.; Liu, J. Coronavirus disease 2019 (COVID-19) CT findings: A systematic review and meta-analysis. J. Am. Coll. Radiol. 2020, 17, 701. [Google Scholar] [CrossRef]
  10. Diwakar, M.; Singh, P. CT image denoising using multivariate model and its method noise thresholding in non-subsampled shearlet domain. Biomed. Signal Processing Control 2020, 57, 101754. [Google Scholar] [CrossRef]
  11. Routray, S.; Malla, P.P.; Sharma, S.K.; Panda, S.K.; Palai, G. A new image denoising framework using bilateral filtering based non-subsampled shearlet transform. Optik 2020, 216, 164903. [Google Scholar] [CrossRef]
  12. Goyal, B.; Dogra, A.; Agrawal, S.; Sohi, B.S. Two-dimensional gray scale image denoising via morphological operations in NSST domain and bitonic filtering. Future Gener. Comput. Syst. 2018, 82, 158–175. [Google Scholar] [CrossRef]
  13. Panigrahi, S.K.; Gupta, S.; Sahu, P.K. Curvelet-based multiscale denoising using non-local means and guided image filter. IET Image Processing 2018, 12, 909–918. [Google Scholar] [CrossRef]
  14. Liu, X.; Tang, Y.; Yang, Y. Primal-dual algorithm to solve the constrained second-order total generalized variational model for image denoising. J. Electron. Imaging 2019, 28, 043017. [Google Scholar] [CrossRef]
  15. Dhabal, S.; Chakrabarti, R.; Mishra, N.S.; Venkateswaran, P. An improved image denoising technique using differential evolution-based salp swarm algorithm. Soft Computing 2020, 25, 1941–1961. [Google Scholar] [CrossRef]
  16. Liu, Y.; Li, S.; Zhang, H. Multibaseline Interferometric Phase Denoising Based on Kurtosis in the NSST Domain. Sensors 2020, 20, 551. [Google Scholar] [CrossRef] [Green Version]
  17. Al-Mefty, O.; Kersh, J.E.; Routh, A.; Smith, R.R. The long-term side effects of radiation therapy for benign brain tumors in adults. J. Neurosurg. 1990, 73, 502–512. [Google Scholar] [CrossRef]
  18. Kursun, S.; Öztas, B.; Atas, H.; Tastekin, M. Effects of X-rays and magnetic resonance imaging on mercury release from dental amalgam into artificial saliva. Oral Radiol. 2014, 30, 142–146. [Google Scholar] [CrossRef]
  19. Thanh, D.; Prasath, S. A review on CT and X-ray images denoising methods. Informatica 2019, 43, 151–159. [Google Scholar] [CrossRef] [Green Version]
  20. Fan, L.; Zhang, F.; Fan, H.; Zhang, C. Brief review of image denoising techniques. Vis. Comput. Ind. Biomed. Art 2019, 2, 7. [Google Scholar] [CrossRef]
  21. Tang, N.; Zhao, X.; Li, Y.; Zhu, D. Adaptive threshold shearlet transform for surface microseismic data denoising. J. Appl. Geophys. 2018, 153, 64–74. [Google Scholar] [CrossRef]
  22. Xu, W.; Tang, C.; Gu, F.; Cheng, J. Combination of oriented partial differential equation and shearlet transform for denoising in electronic speckle pattern interferometry fringe patterns. Appl. Opt. 2017, 56, 2843–2850. [Google Scholar] [CrossRef] [PubMed]
  23. Gibert, X.; Patel, V.M.; Labate, D.; Chellappa, R. Discrete shearlet transform on GPU with applications in anomaly detection and denoising. EURASIP J. Adv. Signal Processing 2014, 2014, 64. [Google Scholar] [CrossRef]
  24. Yang, H.Y.; Wang, X.Y.; Niu, P.P.; Liu, Y.C. Image denoising using nonsubsampled shearlet transform and twin support vector machines. Neural Netw. 2014, 57, 152–165. [Google Scholar] [CrossRef] [PubMed]
  25. Shahdoosti, H.R.; Khayat, O. Image denoising using sparse representation classification and non-subsampled shearlet transform. Signal Image Video Process. 2016, 10, 1081–1087. [Google Scholar] [CrossRef]
  26. Zhang, C.; van der Baan, M. Multicomponent microseismic data denoising by 3D shearlet transform. Geophysics 2018, 83, A45–A51. [Google Scholar] [CrossRef]
  27. Wang, J.; Bao, Y.; Wen, Y.; Lu, H.; Luo, H.; Xiang, Y.; Qian, D. Prior-Attention Residual Learning for More Discriminative COVID-19 Screening in CT Images. IEEE Trans. Med. Imaging 2020, 39, 2572–2583. [Google Scholar] [CrossRef]
  28. Hu, S.; Gao, Y.; Niu, Z.; Jiang, Y.; Li, L.; Xiao, X.; Yang, G. Weakly supervised deep learning for COVID-19 infection detection and classification from ct images. IEEE Access 2020, 8, 118869–118883. [Google Scholar] [CrossRef]
  29. Sun, L.; Mo, Z.; Yan, F.; Xia, L.; Shan, F.; Ding, Z.; Shen, D. Adaptive feature selection guided deep forest for COVID-19 classification with chest ct. IEEE J. Biomed. Health Inform. 2020, 24, 2798–2805. [Google Scholar] [CrossRef]
  30. Shiri, I.; Akhavanallaf, A.; Sanaat, A.; Salimi, Y.; Askari, D.; Mansouri, Z.; Zaidi, H. Ultra-low-dose chest CT imaging of COVID-19 patients using a deep residual neural network. Eur. Radiol. 2020, 31, 1420–1432. [Google Scholar] [CrossRef]
  31. Zhao, W.; Zhong, Z.; Xie, X.; Yu, Q.; Liu, J. Relation between chest CT findings and clinical conditions of coronavirus disease (COVID-19) pneumonia: A multicenter study. Am. J. Roentgenol. 2020, 214, 1072–1077. [Google Scholar] [CrossRef] [PubMed]
  32. Tan, Y.; Wang, X.; Yang, W.; Cheng, Z.; Cao, Q.; Pan, A.; Qin, L. COVID-19 patients with progressive and non-progressive CT manifestations. Radiol. Infect. Dis. 2020, 7, 97–105. [Google Scholar] [CrossRef] [PubMed]
  33. Abolyazid, S.; Alshareef, S.; Abdullah, N.; Khalil, A.; Hamza, N.; Salem, A. COVID-19 pneumonia identified by CT of the abdomen: A report of three emergency patients presenting with abdominal pain. Radiol. Case Rep. 2020, 15, 2098–2103. [Google Scholar] [CrossRef] [PubMed]
  34. Ambrosetti, M.C.; Battocchio, G.; Zamboni, G.A.; Fava, C.; Tacconelli, E.; Mansueto, G. Rapid onset of bronchiectasis in COVID-19 Pneumonia: Two cases studied with CT. Radiol. Case Rep. 2020, 15, 2098–2103. [Google Scholar] [CrossRef]
  35. Johnson, L.N.; Vesselle, H. COVID-19 in an asymptomatic patient undergoing FDG PET/CT. Radiol. Case Rep. 2020, 15, 1809–1812. [Google Scholar] [CrossRef] [PubMed]
  36. Ferrando-Castagnetto, F.; Wakfie-Corieh, C.G.; García, A.M.B.; García-Esquinas, M.G.; Caro, R.M.C.; Delgado, J.L.C. Incidental and simultaneous finding of pulmonary thrombus and COVID-19 pneumonia in a cancer patient derived to 18F-FDG PET/CT. New pathophysiological insights from hybrid imaging. Radiol. Case Rep. 2020, 15, 1803–1805. [Google Scholar] [CrossRef]
  37. Liu, H.; Liu, F.; Li, J.; Zhang, T.; Wang, D.; Lan, W. Clinical and CT imaging features of the COVID-19 pneumonia: Focus on pregnant women and children. J. Infect. 2020, 80, e7–e13. [Google Scholar] [CrossRef]
  38. Wang, K.; Kang, S.; Tian, R.; Zhang, X.; Wang, Y. Imaging manifestations and diagnostic value of chest CT of coronavirus disease 2019 (COVID-19) in the Xiaogan area. Clin. Radiol. 2020, 75, 341–347. [Google Scholar] [CrossRef] [Green Version]
  39. Brogna, B.; Bignardi, E.; Salvatore, P.; Alberigo, M.; Brogna, C.; Megliola, A.; Musto, L. Unusual presentations of COVID-19 pneumonia on CT scans with spontaneous pneumomediastinum and loculated pneumothorax: A report of two cases and a review of the literature. Heart Lung 2020, 49, 864–868. [Google Scholar] [CrossRef]
  40. Ardakani, A.A.; Kanafi, A.R.; Acharya, U.R.; Khadem, N.; Mohammadi, A. Application of deep learning technique to manage COVID-19 in routine clinical practice using CT images: Results of 10 convolutional neural networks. Comput. Biol. Med. 2020, 212, 103795. [Google Scholar] [CrossRef]
  41. Skalidis, I.; Nguyen, V.K.; Bothorel, H.; Poli, L.; Da Costa, R.R.; Younossian, A.B.; Kherad, O. Unenhanced computed tomography (CT) utility for triage at the emergency department during COVID-19 pandemic. Am. J. Emerg. Med. 2020, 46, 260–265. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, G.; Liu, X.; Li, C.; Xu, Z.; Ruan, J.; Zhu, H.; Zhang, S. A noise-robust framework for automatic segmentation of COVID-19 pneumonia lesions from ct images. IEEE Trans. Med. Imaging 2020, 39, 2653–2663. [Google Scholar] [CrossRef] [PubMed]
  43. Isola, P.; Zhu, J.Y.; Zhou, T.; Efros, A.A. Image-to-image translation with conditional adversarial networks. In Proceedings of the IEEE Conference on Computer vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 1125–1134. [Google Scholar]
  44. Bhonsle, D.; Chandra, V.K.; Sinha, G.R. De-noising of CT Images using Combined Bivariate Shrinkage and Enhanced Total Variation Technique. i-Manag. J. Electron. Eng. 2018, 8, 12. [Google Scholar] [CrossRef]
  45. Leal, L.; Castillo, M.; Juarez, F.; Ramirez, E.; Aspuac, M.; Letona, D. Convolutional-LSTM for Multi-Image to Single Output Medical Prediction. arXiv 2020, arXiv:2010.10004. [Google Scholar]
  46. Kallel, F.; Hamida, A.B. A new adaptive gamma correction based algorithm using DWT-SVD for non-contrast CT image enhancement. IEEE Trans. Nanobioscience 2017, 16, 666–675. [Google Scholar] [CrossRef] [PubMed]
  47. Sinha, G.R. An Optimized Framework Using Adaptive Wavelet Thresholding and Total Variation Technique for De-noising Medical Images. J. Adv. Rrsearch Dyn. Control. Syst. 2018, 10, 953–965. [Google Scholar]
  48. Adachi, T.; Takata, T. A three-dimensional cross-directional bilateral filter for edge-preserving noise reduction of low-dose computed tomography images. Comput. Biol. Med. 2019, 111, 103353. [Google Scholar] [CrossRef]
  49. Zhang, B.; Allebach, J.P. Adaptive Bilateral Filter for Sharpness Enhancement and Noise Removal. IEEE Trans. Image Processing 2008, 17, 664–678. [Google Scholar] [CrossRef] [PubMed]
  50. Deng, L.; Zhu, H.; Yang, Z.; Li, Y. Hessian matrix-based fourth-order anisotropic diffusion filter for image denoising. Opt. Laser Technol. 2019, 110, 184–190. [Google Scholar] [CrossRef]
  51. Pang, Z.F.; Zhou, Y.M.; Wu, T.; Li, D.J. Image denoising via a new anisotropic total-variation-based model. Signal Process. Image Commun. 2019, 74, 140–152. [Google Scholar] [CrossRef]
  52. Rabbouch, H.; Saadaoui, F. A wavelet-assisted subband denoising for tomographic image reconstruction. J. Vis. Commun. Image Represent. 2018, 55, 115–130. [Google Scholar] [CrossRef]
  53. Islam, M.Z.; Islam, M.M.; Asraf, A. A combined deep CNN-LSTM network for the detection of novel coronavirus (COVID-19) using X-ray images. Inform. Med. Unlocked 2020, 20, 100412. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Proposed algorithm.
Figure 1. Proposed algorithm.
Diagnostics 12 02766 g001
Figure 2. COVID-19 reference CT image dataset; (a) CT1 image; (b) CT2 image; (c) CT3 image; (d) CT4 image.
Figure 2. COVID-19 reference CT image dataset; (a) CT1 image; (b) CT2 image; (c) CT3 image; (d) CT4 image.
Diagnostics 12 02766 g002
Figure 3. COVID-19 noisy CT image dataset corresponding to Figure 2; (a) Noisy CT1 image; (b) Noisy CT2 image; (c) Noisy CT3 image; (d) Noisy CT4 image.
Figure 3. COVID-19 noisy CT image dataset corresponding to Figure 2; (a) Noisy CT1 image; (b) Noisy CT2 image; (c) Noisy CT3 image; (d) Noisy CT4 image.
Diagnostics 12 02766 g003
Figure 4. Denoised COVID-19 CT image 1 from the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Figure 4. Denoised COVID-19 CT image 1 from the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Diagnostics 12 02766 g004
Figure 5. Denoised COVID-19 CT image 2 for the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Figure 5. Denoised COVID-19 CT image 2 for the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Diagnostics 12 02766 g005
Figure 6. Denoised COVID-19 CT image 3 from the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Figure 6. Denoised COVID-19 CT image 3 from the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Diagnostics 12 02766 g006
Figure 7. Denoised COVID-19 CT image 4 for the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Figure 7. Denoised COVID-19 CT image 4 for the (a) [10], (b) [11], (c) [12], (d) [13], (e) [14], (f) [11], (g) [16], and (h) proposed frameworks, respectively.
Diagnostics 12 02766 g007
Figure 8. The line used for intensity profiling in Image 1 for all the frameworks.
Figure 8. The line used for intensity profiling in Image 1 for all the frameworks.
Diagnostics 12 02766 g008
Figure 9. Intensity profiles of the original image against the [10,11,12,13,14,15,16] framework and the proposed framework, respectively; (a) Intensity profile between reference image and denoised image [10]; (b) Intensity profile between reference image and denoised image [11]; (c) Intensity profile between reference image and denoised image [12]; (d) Intensity profile between reference image and denoised image [13]; (e) Intensity profile between reference image and denoised image [14]; (f) Intensity profile between reference image and denoised image [15]; (g) Intensity profile between reference image and denoised image [16]; (h) Intensity profile between reference image and denoised image of proposed method.
Figure 9. Intensity profiles of the original image against the [10,11,12,13,14,15,16] framework and the proposed framework, respectively; (a) Intensity profile between reference image and denoised image [10]; (b) Intensity profile between reference image and denoised image [11]; (c) Intensity profile between reference image and denoised image [12]; (d) Intensity profile between reference image and denoised image [13]; (e) Intensity profile between reference image and denoised image [14]; (f) Intensity profile between reference image and denoised image [15]; (g) Intensity profile between reference image and denoised image [16]; (h) Intensity profile between reference image and denoised image of proposed method.
Diagnostics 12 02766 g009aDiagnostics 12 02766 g009b
Table 1. Average results (87 images) of performance metrics (PSNR, SSIM, DIV, and ED).
Table 1. Average results (87 images) of performance metrics (PSNR, SSIM, DIV, and ED).
Proposed Method (without Circular Shifting)Proposed Method (with Patch-Wise Circular Shifting)
PSNRSSIMPSNRSSIM
Noise Variance 3 × 35 × 57 × 79 × 93 × 35 × 57 × 79 × 9
1028.010.740828.329.428.828.50.76330.78550.75340.7432
2026.130.676326.527.427.626.20.69330.70110.69226823
3024.210.602424.625.725.524.70.61240.63210.62110.6121
4022.030.573222.323.823.922.20.59110.60620.59820.5827
EDDIVEDDIV
100.81310.57410.42320.32310.31230.43410.44020.31320.47130.4928
201.93121.67421.75421.03411.12121.33121.23361.03421.18721.4342
303.62452.53223.45212.92313.43253.32252.12212.00222.32112.2232
405.78454.22155.14315.10425.25115.42254.19324.03154.19284.2051
Table 2. PSNR and SSIM of noiseless CT images.
Table 2. PSNR and SSIM of noiseless CT images.
InputGaussian Noisy COVID-19 CT Image Dataset 1 (Average Results on 90 Images)
PSNRSSIM
1020304010203040
[10]29.354427.435625.344622.34210.73440.66930.61330.5993
[11]29.443427.454725.342622.235510.78880.66740.61960.5977
[12]29.442527.323725.566722.34230.74660.64520.61040.5951
[13]29.644627.233425.246822.32240.73430.63570.61240.5918
[14]29.233627.344425.345822.74550.73540.64320.61240.5951
[15]29.454727.544825.743622.56350.75480.63470.61540.5918
[16]29.123727.455325.563322.45310.73550.68590.61940.5905
Proposed30.023828.133925.183923.22390.79720.69930.62070.6021
Table 3. ED and DIV of noiseless CT images.
Table 3. ED and DIV of noiseless CT images.
InputGaussian Noisy COVID-19 CT Image Dataset 1 (Average Results on 90 Images)
EDDIV
Σ1020304010203040
[10]0.58211.32932.34483.45330.35411.43532.34434.3422
[11]0.68381.34642.23463.23370.44381.45442.34224.2355
[12]0.38641.56722.45643.45510.44241.32322.56664.3422
[13]0.68281.12772.46643.32280.64481.23372.24644.3222
[14]0.78641.43722.56643.24510.23341.34422.34564.7454
[15]0.78281.56772.64343.34580.45481.54472.74334.5633
[16]0.48551.67592.23543.56250.12351.45592.56324.4535
Proposed0.29821.06942.05773.01210.02321.13342.18323.2236
Table 4. Visual analysis from the experts.
Table 4. Visual analysis from the experts.
[10][11][12][13][14][15][16]Proposed
Edge Preservation12233233
Noise Suppression12322323
Artifact Removal22123333
Heterogeneous Consistency12212323
Data retention21222223
Medical feature Extraction33232333
Clarity/Sharpness12123333
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Diwakar, M.; Singh, P.; Swarup, C.; Bajal, E.; Jindal, M.; Ravi, V.; Singh, K.U.; Singh, T. Noise Suppression and Edge Preservation for Low-Dose COVID-19 CT Images Using NLM and Method Noise Thresholding in Shearlet Domain. Diagnostics 2022, 12, 2766. https://doi.org/10.3390/diagnostics12112766

AMA Style

Diwakar M, Singh P, Swarup C, Bajal E, Jindal M, Ravi V, Singh KU, Singh T. Noise Suppression and Edge Preservation for Low-Dose COVID-19 CT Images Using NLM and Method Noise Thresholding in Shearlet Domain. Diagnostics. 2022; 12(11):2766. https://doi.org/10.3390/diagnostics12112766

Chicago/Turabian Style

Diwakar, Manoj, Prabhishek Singh, Chetan Swarup, Eshan Bajal, Muskan Jindal, Vinayakumar Ravi, Kamred Udham Singh, and Teekam Singh. 2022. "Noise Suppression and Edge Preservation for Low-Dose COVID-19 CT Images Using NLM and Method Noise Thresholding in Shearlet Domain" Diagnostics 12, no. 11: 2766. https://doi.org/10.3390/diagnostics12112766

APA Style

Diwakar, M., Singh, P., Swarup, C., Bajal, E., Jindal, M., Ravi, V., Singh, K. U., & Singh, T. (2022). Noise Suppression and Edge Preservation for Low-Dose COVID-19 CT Images Using NLM and Method Noise Thresholding in Shearlet Domain. Diagnostics, 12(11), 2766. https://doi.org/10.3390/diagnostics12112766

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