{'fraction_non_alphanumeric': 0.031123723742249917, 'fraction_numerical': 0.0123384435053518, 'mean_word_length': 4.789285714285715, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
COVID-19 CT image denoising algorithm based on adaptive threshold and optimized weighted median filter

CT image of COVID-19 is disturbed by impulse noise during transmission and acquisition. Aiming at the problem that the early lesions of COVID-19 are not obvious and the density is low, which is easy to confuse with noise. A median filtering algorithm based on adaptive two-stage threshold is proposed to improve the accuracy for noise detection. In the advanced stage of ground-glass lesion, the density is uneven and the boundary is unclear. It has similar gray value to the CT images of suspected COVID-19 cases such as adenovirus pneumonia and mycoplasma pneumonia (reticular shadow and strip shadow). Aiming at the problem that the traditional weighted median filter has low contrast and fuzzy boundary, an adaptive weighted median filter image denoising method based on hybrid genetic algorithm is proposed. The weighted denoising parameters can adaptively change according to the detailed information of lung lobes and ground-glass lesions, and it can adaptively match the cross and mutation probability of genetic combined with the steady-state regional population density, so as to obtain a more accurate COVID-19 denoised image with relatively few iterations. The simulation results show that the improved algorithm under different density of impulse noise is significantly better than other algorithms in peak signal-to-noise ratio (PSNR), image enhancement factor (IEF) and mean absolute error (MSE). While protecting the details of lesions, it enhances the ability of image denoising.


Background & problem domain
Since December 2019, COVID-19 has spread widely around the world. According to the data released by the Johns Hopkins University, as of September 18, 2020, there were 30065728 confirmed cases in the world. The disease is highly infectious and it can cause severe acute dyspnea. Studies have shown that computed tomography (CT), as a noninvasive imaging method, is of great value in detecting lung lesions in patients with COVID-19. Therefore, CT can be used for detection and diagnosis of COVID-19. However, CT images will produce high-dose radiation in the process of acquisition, which will pose a great threat to the patient's health. Therefore, at present, the damage to the patient's body is usually reduced by decreasing the CT dose. But CT images taken with low dose usually have noise. The generation of noise will affect the quality of CT images, which will seriously affect the doctor's diagnosis for the patient's condition. Lung CT images are inevitably disturbed by various noises in the process of acquisition, transmission and storage, which will cause some random, discrete and isolated pixels on the image, that is, image noise will bring bad interference to lesion image segmentation and doctors' judgment of patients' condition.

Review of literature
At present, nucleic acid detection is the most common method to diagnose COVID-19 [1]. This method combines RNA reverse transcription and polymerase chain reaction (RT-PCR) to detect the viral RNA fragments [2,3], and the diagnosis can be confirmed by the positive nucleic acid test. However, the screening of RT-PCR has the problem of low sensitivity. Even if the RT-PCR result of suspected patient is negative, the possibility of SARS-Cov-2 infection cannot be completely ruled out [4]. In addition, nucleic acid testing has the disadvantage of timeconsuming and it requires a special test kits. Therefore, it is necessary to further accelerate the detection speed and reduce the cost [5]. Through the comparative experimental analysis of CT and RT-PCR, Ai [6] pointed out that the diagnosis based on chest CT image is faster and more efficient than RT-PCR. Generally speaking, the lungs of patients with COVID-19 have typical imaging features [7], including groundglass opacities (GGO), pulmonary sclerosis, pulmonary fibrosis and multiple lesions. In addition, related studies [8] have shown that the imaging information can play a vital role in the diagnosis of COVID-19. The autonomous diagnosis of COVID-19 based on artificial intelligence can minimize the contact between doctors and the virus, which reduces the chance of infection [9,10]. However, the manual analysis and diagnosis process based on CT images are highly dependent on professional knowledge [11], and the analysis of CT image features is timeconsuming, which is difficult to observe hidden lesions in the early stage, and it is hard to distinguish other viral pneumonia and bacterial pneumonia [12]. The automatic diagnosis of COVID-19 can convert the visual image into the deep-level feature information [13,14], which can be quantified. It helps to reduce manual operation and improve the efficiency of accurate quantitative analysis. Standard median filter can preferably filter impulse noise. However, the traditional median filter takes the median value of the neighborhood pixels, the information points in the image may be filtered out and replaced by noise points, which will lead to the loss of some details in the image. Therefore, many scholars have proposed some improved algorithms. Verma [15] proposed adaptive median filter (AMF) algorithm, which shows good filtering performance in low-intensity noise, but when the noise density is too high, a large number of noises needs to be filtered, and the filtering window will expand. The current pixel and the replaced median may lack correlation, resulting in information loss and image blur. Sheela [16] proposed the switching median filter (SMF) method. This method determines whether each pixel is noisy. Due to the threshold needs to be defined in advance when judging the noise, the threshold remains unchanged in the processing process, the final filtering effect is not very ideal. Gupta [17] proposed a discrimination-based dual threshold median filter (DBF) method to judge whether the pixel is noisy by comparing the current pixel with the dual threshold. However, when dealing with high-intensity noise, this method is easy to produce "tailing effect", and the filtering performance is greatly reduced. Ahmed [18] proposed an improved discriminant-based asymmetric clipping median filter algorithm. This method can obtain good filtering effect under lowintensity noise, but under high-intensity, the filtering effect is very poor and some details are lost. There are some state-of-the-art methods about denoising algorithm reported. Zhao [19] obtains the residual spectrum according to logarithmic amplitude spectrum of the image, and then the inverse Fourier transform is used to obtain the saliency map. Although this method can suppress the image noise, the subtraction operation is used in the process of calculating the log amplitude. It not only suppresses the noise, but also restrains the detail information, which is not conducive to the denoising of the COVID-19 image. The CNN denoising method [20] can effectively use the global features of the image to significantly improve the denoising effect. The lower the number of CNN layers, the more primitive the extracted low-level features, such as color and edge lines. However, with the deepening of the network, the features extracted by the CNN method become more advanced, and a large number of low-level features will be lost. Yan [21] proposed a KSVD image denoising algorithm based on K-means and singular value decomposition (SVD). This method is based on the bayesian reconstruction theory, and it uses the K-means clustering to obtain an optimal signal sparse representation dictionary set. Through the linear combination of the dictionary elements, various signals are constructed to remove the noise in the CT image. However, dictionary learning brings high computational cost, resulting in too long time for image denoising. In addition, the KSVD method is also insufficient for denoising of highproportion noisy images. Yahya [22] proposed block-matching and 3D filtering (BM3D) algorithm, which has achieved good results in various image denoising processes, but it is easy to cause tissue contour blur when denoising low-dose CT image. Moreover, there are a lot of parameters in the BM3D algorithm that need to be set manually, which severely limits the effect of the practical application. Gu [23] proposed a weighted kernel norm minimization (WNNM) denoising algorithm. This method can characterize the difference of image effect according to the matrix singular value, given different weights, and it highlights the important information for the image. However, this algorithm needs time-consuming SVD calculation and iterative approximation, which is slow and inefficient.

Gaps identified from review
The above improved algorithms have great improvements in filtering performance, but there are still some deficiencies in dealing with image detail blur. The lesions of COVID-19 mainly show various forms of ground-glass or consolidation shadow. In the early imaging of COVID-19, the change of lesions is not obvious, and the density is low, which is easy to confuse the early lesions of COVID-19 with noise. On the other hand, the lesions of ground-glass in the advanced stage are uneven and it have similar gray values to the CT images of suspected COVID-19 cases such as adenovirus pneumonia and mycoplasma pneumonia. The traditional weighted median filter denoising algorithm has the problems of low contrast and fuzzy boundary.

Highlights
The highlights of this work can be summarized as the following: (i) the adaptive two-stage threshold method is proposed, which reduces the misjudgment of noise signal points and improves the accuracy of detection for COVID-19; (ii) we develop the weighted median filter image denoising method based on adaptive hybrid genetic. The weighted denoising parameters can adaptively change according to the detailed information of lung lobes and ground-glass lesions, it can enhance the ability of image denoising while better protecting the details of the lesions; (iii) the cross and mutation probability of genetic algorithm based on the steady-state regional population density is advanced, which obtains a more accurate COVID-19 denoised image with relatively few iterations; (iv) in order to reduces the CPU time and improve the efficiency of the optimal solution, we investigate the doppler effect temperature improved simulated annealing algorithm.

Paper structure
The rest of this paper is organized as follows. Section 2 describes the multi-level threshold algorithm and its mathematical properties, and we present the algorithm of adaptive center weighted median filtering. In Section 3, we propose the adaptive hybrid genetic method and the advanced simulated annealing algorithm according to these properties. In Section 4, we show the simulation experiments for the different kinds of COVID-19 CT images. The paper is structured as shown in Fig.1.

Multi-level threshold center weighted adaptive median filtering algorithm
The CT image of early COVID-19 is easily confused by noise. The traditional extreme noise points judgment [24][25][26][27] can easily cause misdiagnosis of COVID-19 lesion signals. For the single threshold filtering method, the pixel whose value exceeds the preset threshold will be judged as noise [28][29][30]. Therefore, using a single threshold method may increase the probability of error detection for noise. Hence, a median filtering method based on adaptive multi threshold discrimination is proposed. In this method, the noise pixels will be discriminated in a relatively narrow interval. Thus, this will improve the accuracy of noise detection.

Multi-level threshold calculation method
Calculation of threshold T 0 : this threshold is the first step of noise point detection. The purpose is to find suspected noise points. It can calculate the gray value of each pixel, the median gray value of the filter window, the gray average value of the four rectangular neighborhood windows on the top, bottom, left and right of the point (i, j). The calculation of threshold T 0 is: The variable n and m represent the length and the width of the filter window respectively(the length and width are equal), and f(i,j) represents the gray value of the current pixel.
Calculation of the second threshold T 1 : this threshold is the second step for noise point detection. When the point is determined to be a suspected noise point, the average gray value of all pixels in the filter window is calculated. If the gray value of the point is significantly different from the average value, it can be determined that the current point is a noise point and needs to be replaced by the median value. Through the further accurate judgment for the suspected noise point by the second threshold T 1 , the calculation of T 1 is: The variable φ is the regulatory factor of threshold T 1 ,φ = F 1 + F 2 , the expression of f 1 and f 2 is: p is the image noise density, the expression is p = Npn n×n ; N pn represents the number of noise points in the window. The threshold T 1 increases with the increase of noise density, and the threshold value is related to the noise density, so as to realize the adaptive selection of the threshold in the same image under different noise density.

Adaptive weighted median filter algorithm
The lesions of ground-glass in the advanced COVID-19 show nonuniformity, which is similar to the CT images of COVID-19 suspected cases such as adenovirus pneumonia and mycoplasma pneumonia. Aiming at the problems of low contrast and fuzzy boundary in the traditional median filter denoising algorithm, a weighted median filter image denoising method based on adaptive hybrid genetic algorithm is proposed, so that the weighted denoising parameters can adapt to the changes according to the details of lung lobes and ground-glass lesions. Based on the steady-state regional population density, the cross and the mutation probability of genetic is adaptively matched to obtain a more accurate denoised CT image of COVID-19 with relatively few iterations. The calculation of the center weighted filtering [31][32][33][34] is: represents the gray value when the weight of f(i,j) is w(i,j), f(i, j) CWM is the center weighted median gray value, f(i,j) represents the gray value of pixel p(i,j), N o p(i,j) represents the hollow neighborhood of pixel p(i,j), and f(r,s) represents the gray value corresponding to all pixels in N o p(i,j) . The improved weight of the center weighted filtering is: α 1 and α 2 are weight parameters, which is optimized by the adaptive is used to characterize the similarity between the current pixel and the median pixel f med (representing uncontaminated point). The smaller the similarity value, the greater the correlation. The calculation of the weighted sub coefficient w 1 (i, j) is: T mean represents the average value for the square of the difference between each pixel value and its average value, which is calculated as following: N r is the total number of window pixels and f mean is the average value of all pixels in the current window. A weighted sub coefficient w 2 (i, j) is set in the weighted coefficient w(i,j), and w 2 (i, j) is used to characterize the similarity between the current pixel and the nearest noise point. According to the pixel space distance value, a distance correlation function Q(i,j) is designed to characterize the correlation between the current pixel and the noise points. The weighted coefficient w 2 (i, j) should meet the following conditions:①w 2 (i, j) is a non-decreasing function with Q(i, j), that is, the greater the Q(i,j), the greater the weight value to be added;②w 2 (i, j)When Q(i,j) is small, w(i,j) should increase slightly, and when Q(i,j) is large, w(i,j) should increase greatly, that is, w(i,j) is positively correlated with Q(i,j). Setting the current pixel as f(i,j), and f (s,t) is the noise point pixel nearest to the current pixel. Combined with the weighted distance and its normalization, the following weighted function is obtained, which meets positive correlation and the requirement of human vision. Setting the weighted function is:    Where, c is the distance constant, which is related to the size of the filter window, c = 1 (n×n+1)/2 . Therefore, the c value of 5x5 detection window used in this paper is c = 1 (5×5+1)/2 = 1 13 .Through the optimal combination for (α 1 ,α 2 ), the relatively optimal denoising effect in the spatial domain can be obtained. The pixel weighted coefficient is represented as shown in Fig.3.

Improvement steps
Step 1: find out the maximum gray value f max and minimum gray value f min in the filter window, which is compared with the current pixel , replace the weighted median value with the current point and repeat the above steps. The algorithm flow chart is shown in Fig.4.

Parameters optimization of weighted median filter based on adaptive hybrid genetic algorithm
The principle of the adaptive hybrid genetic optimization is as following: (i) Initializing the search space of X i (α 1 , α 2 ), leading the genetic operators to search with rmX i = X j +φ(Δx) to obtain the new solution rmV i (α 1 , α 2 ); (ii) Using the parameter factors α 1 and α 2 to carry   out the optimal weighted median filter; (iii) Calculating the value of PSNR as the fitness of the genetic algorithm; (iv) Using simulated annealing to cover the solution with high fitness over the low fitness as the current optimal solution. In this way, the cycle stops until the termination condition is established to obtain the optimal solution rmX k (α 1 , α 2 ) corresponding to the highest fitness.

Improved operator strategy of adaptive hybrid genetic
In this paper, an adaptive hybrid genetic operator is proposed. In addition to retaining the "excellent" individuals in the parents to the next generation, individuals with close to average fitness in the parents should also be considered. If the fitness value of the offspring generation extracted in the "steady-state region" is higher than the average and close to the upper bound, it indicates that such parent individuals contain excellent gene patterns, it should be retained to the next gen- The adaptive hybrid mutation probability can be expressed as: k 1 and k 2 are coefficients and k 1 , k 1 ⩽1.The adaptive hybrid cross operation is shown in Fig.5, the adaptive hybrid mutation operation is shown in Fig.6.

Improved simulated annealing algorithm
Derek [35] proposed metropolis simulated annealing algorithm criterion. This paper further improves the cooling function and designs an expression similar to the doppler effect temperature decline curve: The traditional exponential simulated annealing is difficult to reach the low temperature state when the number of iterations is small; the rapid simulated annealing drops to the low temperature state prematurely, making the subsequent iterative solution unchanged. The "doppler" curve takes the advantages of the two algorithms and eliminates their defects. It not only tends to the low temperature unurgency, but also  continuously "tempering and heating up" in the second half of annealing, so that the algorithm has multiple opportunities to jump out of the local optimum during the optimization process, making it easier to find out the global optimal solution. K is the total number of iterations, and k is the current number of iterations. The flow chart of the adaptive hybrid genetic combined with the doppler effect simulated annealing algorithm is shown in Fig.7.

Evaluation indexs
In order to objectively test the denoising effect of each method, mean square error (MSE), peak signal-to-noise ratio (PSNR) and image enhancement factor (IEF) are used to calculate the denoised image. The lower the MSE, the better the denoised image quality; the higher the peak signal-to-noise ratio and image enhancement factor, the better the denoising effect. The expression of MSE is: The expression of PSNR is: The expression of IEF is: f(i,j) represents the denoised pixel, f (i, j) represents the input pixel with impulse noise, x(i,j) represents the original pixel, M and N represents the length and width of the image respectively, the picture size is 512*512.

Denoising simulation experiments
Using the paper method and the comparison methods to do denoising simulations on the CT image of early COVID-19, and further comparing and analyzing the MSE, PSNR and IEF under different denoising methods. The comparative CT image denoising simulations of early COVID-19 are shown in Fig.8. The evaluation index values are shown in Table 1.
From the change trend of the data in the Table 1, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 5.5 dB, the value of MSE has been greatly reduced, the value of IEF has been increased by about 42 dB. It can be seen that the denoising effect for the CT image of early COVID-19 under TWMF is the best. There is no obvious change in the CT image of early COVID-19, the number of the lesion is small and the density is low. The traditional CT image denoising algorithms have the higher value of MSE and the lower value of IEF. It is easy to cause the CT image of early COVID-19 to be confused with noise, and it is easy to cause mistake diagnosis for patients with early COVID-19. The TWMF method improves the denoising accuracy for the CT image of early COVID-19. The CT image of early COVID-19 is simulated by denoising tests with different density (10%-70%) of impulse noise. The comparison for the MSE and IEF among different denoising methods by different density of the impulse noise is shown in Fig.13(a)(b). The evaluation index values are shown in Table 6.
When the value of MSE is smaller and the value of IEF is larger, it means that the denoised image is closer to the original image. It can be seen from the Fig.13(a)(b) that with the increase of the impulse noisy density, the PSNR is gradually reduced and the MES is gradually increased under different denoising methods. Compared with the different traditional denoising methods, the IEF of TWMF is still the largest, and the MSE of TWMF is still the smallest. In conclusion, compared with the traditional denoising methods, the improved median filter denoising algorithm in this paper increases the IEF and reduces the MSE for CT image of early COVID-19 to a certain extent with different noise density. It significantly improves the denoised effect. Using the TWMF and the comparison methods to do denoising simulations on the CT image of advanced COVID-19. The comparative CT image denoising simulations of advanced COVID-19 are shown in Fig.9. The evaluation index values are shown in Table 2.
From the change trend of the data in the Table 2, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 5.7 dB, the value of MSE has been greatly reduced, the value of IEF has been increased by about 45 dB. It can be seen that the denoising effect for the CT image of advanced COVID-19 under TWMF is best. The advanced stage of COVID-19 is characterized by increased number of lesions and increased density. The TWMF method improves the denoising accuracy for the CT image of advanced COVID-19 and it reduces the misdiagnosis rate for the advanced lesions. The CT image of advanced COVID-19 is simulated by denoising tests with different density(10%-70%) of impulse noise. The comparison for the MSE and IEF among different denoising methods with different density of the impulse noise is shown in Fig.13(c)(d).The evaluation index values are shown in Table 7.
It can be seen from the Fig.13(c)(d) and Table 7 that with the increase of the impulse noisy density, the IEF is gradually reduced and the MES is gradually increased under different denoising methods. Compared with the different traditional denoising methods, the IEF of TWMF is still the largest, and the MSE of TWMF is still the smallest. In conclusion, the TWMF significantly improves the denoised effect for CT image of advanced COVID-19 to a certain extent with different noisy density. Using the TWMF and the comparison methods to do denoising simulations on the CT image of adenovirus pneumonia (suspected cases of . The comparative CT image denoising simulations of adenovirus pneumonia are shown in Fig.10. The evaluation index values are shown in Table 3. From the change trend of the data in the Table 3, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 5.4 dB, the value of MSE has been greatly reduced, the value of IEF has been increased by about 41 dB. It can be seen that the denoising effect for the CT image of adenovirus pneumonia under TWMF is the best. The CT image of adenovirus pneumonia is patchy, small nodular ground-glass shadow, and then developing into consolidation shadow, bronchial wall thickening, visible mesh or strip shadow. The TWMF improves the denoising accuracy for the CT image of adenovirus pneumonia and reduces the misdiagnosis rate for the lesions. The CT image of adenovirus pneumonia are simulated by denoising tests with different density (10%-70%) of impulse noise. The comparison for the MSE and IEF among different methods with different density of the impulse noise is shown in Fig.13(e)(f). The evaluation index values are shown in Table 8.
It can be seen from the Fig.13(e)(f) and Table 8 that with the increase of the impulse noisy density, the IEF is gradually reduced and the MES is gradually increased under different denoising methods. Compared with the different traditional denoising methods, the IEF of TWMF is still the largest, and the MSE of TWMF is still the smallest. In conclusion, the TWMF method significantly improves the denoised effect for CT image of adenovirus pneumonia to a certain extent with different noisy density. Using the TWMF and the comparison methods to do denoising simulations on the CT image with non-symptom of early COVID-19. The comparative CT image denoising simulations are shown in Fig.11. The evaluation index values are shown in Table 4.
From the change trend of the data in the Table 4, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 5 dB, the value of MSE has been greatly reduced, the value of IEF has been increased by about 47 dB. It can be seen that the denoising effect for the CT image of early COVID-19 with non-symptom under TWMF is the best. The CT image of early COVID-19 with non-symptom are simulated by denoising test with different density (10%-70%) of impulse noise. The comparison for the MSE and IEF among different denoising methods with different density of the impulse noise is shown in Fig.13(g)(h). The evaluation index values are shown in Table 9.
It can be seen from the Fig.13(g)(h) and Table 9 that with the increase of the impulse noisy density, the IEF is gradually reduced and the MES is gradually increased under different denoising methods. Compared with the different traditional denoising methods, the IEF of TWMF is still the largest, and the MSE of TWMF is still the smallest. In conclusion, the TWMF significantly improves the denoised effect for CT image of early COVID-19 with non-symptom to a certain extent by different noisy density. Using the TWMF and the comparison methods to do denoising simulations on the CT image of advanced COVID-19 with non-symptom. The comparative CT image denoising simulations of advanced COVID-19 with non-symptom are shown in Fig.12. The evaluation index values are shown in Table 5.
From the change trend of the data in the Table 5, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 5.0 dB, the value of MSE has been greatly reduced, the value of IEF has been increased by about 60 dB. It can be seen that the denoising effect for the CT image of advanced COVID-19 with non-symptom under TWMF is the best. The CT image of advanced COVID-19 with nonsymptom is simulated by denoising test with different density (10%-70%) of impulse noise. The comparison for the MSE and IEF among different denoising methods with different density of the impulse noise is shown in Fig.13(i)(j).The evaluation index values are shown in Table 10.
It can be seen from the Fig.13(i)(j) and Table 10 that with the increase of the impulse noisy density, the IEF is gradually reduced and the MES is gradually increased under different denoising methods.
Compared with the different traditional denoising methods, the IEF of TWMF is still the largest, and the MSE of TWMF is still the smallest. In conclusion, the paper method significantly improves the denoised effect for CT image of advanced COVID-19 with non-symptom to a certain extent by different noisy density. Aiming at the different kinds of COVID-19 CT imagesadaptive hybrid genetic combined with exponential simulated annealing [41] (AHG + ESA) orthogonal adaptive genetic [42] combined with rapid simulated annealing [43] (OAG + RSA) traditional adaptive genetic [44] combined with exponential simulated annealing (TAG + ESA) adaptive hybrid genetic combined with doppler effect simulated annealing (AHG + DESA, paper method) are used for simulation comparison tests of denoising parameters optimization. The fitness evolution curve for different COVID-19 CT images are shown in Fig.14. The evaluation index values of time processing under different methods are shown in Table 11.
It can be seen from the Fig.14 and Table 11 that with the increase of the number of iterations, the standard deviation of the objective function value is gradually reduced and the processing time is gradually increased under different parameters optimization methods. Compared with the different comparison methods, the number of iterations of "HG + DESA" is still the smallest, and the processing time of "AHG + DESA" is still the fastest. In conclusion, the paper method (AHG + DESA) significantly improves the denoised effect for different kinds of COVID-19 CT images of parameters optimization.

Comparison to state-of-the-art methods
We compared the TWMF with 9 state-of-the-art approaches including TWSC [45], KSVD [21], activity-tuned image filtering(ATIF) [19], WNNM [23], BM3D [22], SAINT [46], DAAM [47], TNRD [48], DnCNN [20] to do denoising simulations on the CT image of early COVID-19, and further comparing and analyzing the MSE, PSNR and IEF under different denoising methods. The comparative CT image denoising simulations of early COVID-19 are shown in Fig.15. The evaluation index values are shown in Table 12. From the change trend of the data in the Table 12, compared with the comparative denoising methods, TWMF has increased the value of PSNR by about 1.2 dB, the value of MSE has been reduced by about 27 dB, the value of IEF has been increased by about 21 dB. It can be seen that the denoising effect for the CT image of early COVID-19 under TWMF is the best. The early stage of COVID-19 is characterized by less number of lesions and less density. The TWMF method improves the denoising accuracy for the CT image of early COVID-19 and reduces the missed diagnosis for the early lesions.

Conclusion
In this paper, a median filtering method based on multi-level threshold and parameters optimization is proposed. In this method, the pixels in the horizontal and vertical directions of the pixels to be processed, the multi-level threshold for detecting impulse noise is selected by comprehensive estimation, because the threshold changes with the movement of the window, which reflects the adaptability of threshold selection. In this paper, the algorithm gives different weight values to each pixel according to the gray distance between the pixel and the central pixel and the spatial distance between the pixel and the noise point to filter out the noise. The adaptive hybrid genetic algorithm combined with the doppler effect simulated annealing algorithm is used to determine the weight value of the improved median filter. The simulation results show that this method has greater advantages in denoising ability than other filtering methods of the iteration and the processing time, and the performance of detail protection is also greatly improved.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.