A New Adaptive Spatial Filtering Method in the Wavelet Domain for Medical Images

Although there are many methods in the literature to eliminate noise from images, finding new methods remains a challenge in the field and, despite the complexity of existing methods, many of the methods do not reach a sufficient level of applicability, most often due to the relatively high calculation time. In addition, most existing methods perform well when the processed image is adapted to the algorithm, but otherwise fail or results in significant artifacts. The context of eliminating noise from images is similar to that of improving images and for this reason some notions necessary to understand the proposed method will be repeated. An adaptive spatial filter in the wavelet domain is proposed by soft truncation of the wavelet coefficients with threshold value adapted to the local statistics of the image and correction based on the hierarchical correlation map. The filter exploits, in a new way, both the inter-band and the bandwidth dependence of the wavelet coefficients, considering the minimization of computational resources.


State of the Art
Several statistical models for images have been proposed in the wavelet specialized literature. The simplest and most widespread model considers wavelet coefficients to be independent and identically distributed according to a generalized Gaussian distribution [1][2][3][4], a model successfully used in noise removal and image restoration applications. This model approximates the first-order statistics of wavelet coefficients, but does not take into account higher-order statistics, which leads to some limitations. The existing dependencies between the wavelet coefficients have been studied in order to implement efficient compression algorithms and they were explicitly formulated in [5] or implicitly in [6]. Most statistical models in the wavelet domain of images can be grouped into two categories: models that exploit inter-scale dependencies and models that exploit intra-scale dependencies.
Despite the processing complications, it has been shown [7][8][9][10] that using both types of dependencies leads to better results in white noise removal, compared to the cases where only one of the two dependencies are operated. importance in the wavelet transforms, especially for de-noising signals. Alfaouri and Daqrouq propose [49] a method of choosing the optimal threshold value for de-noising an electrocardiogram (ECG) signal.
X-ray imaging can also benefit from the advantages of wavelet domain. In [50], Purisha et al., use a controlled orthonormal wavelet domain for the reconstruction of a tomographic image.
Unwanted noise in clinical images and operations can be measured using several techniques [51,52]. These measurements are useful in assessing the potential denoising power when applying different denoising methods.
Dabov et al., extended the applications of denoising algorithms to color images through the use of a collaborative Wiener filter [53], based on block-matching fragments of the image into 3-dimensional group arrays (BM3D), thus obtaining very good PSNR (peak signal to noise ratio) results. A fast method for denoising images is the Non-Local Means (NLM) strategy, which was presented in [54], taking into consideration three aspects: noise, quality and MSE (Mean Square Error). This algorithm performed best in the case of textured images.
The methods proposed in [53,54] are applied for noise reduction in images. The results obtained in [53,54] are good, with applicability for the analysis of macroscopic images, the noise, especially the white Gaussian noise, is significantly reduced, but the contour is degraded. The algorithm proposed in this article represents a first stage in medical image processing, and the results obtained will be used in the development of a sub-pixelation algorithm with neural networks, which will be developed in the second stage.
Currently, medical imaging analyzes tumors on a macroscopic level. The interest in medical imaging is the transition from macroscopic to microscopic analysis by using image sub-pixelation algorithms so that the exact tumor contours in medical images can be determined.
The algorithm proposed by us achieves a good ratio between noise reduction and conservation of contours in medical images, taking into account the conservation of Hounsfield coefficients (radiodensity) used in medical imaging. The results obtained with this algorithm will be used in the future to develop a method of sub-pixelation of tumor regions to develop a new algorithm called re-segmentation using neural networks.
Applying the described algorithm on some medical images obtained by CT and magnetic resonance, it is found that the noise was reduced by up to 80%, and the maintenance of the details in the images was at least 80%. The proposed filter is superior to the Wiener filter proposed in [53] and [54], due to the sensitivity to contours and other important features in the images, retaining some high frequency data in the image. The direct spatial correlation of the undecided wavelet transform on several adjacent scales obviously improves the important contours in the image, while the noise and non-essential (small) features in the image are suppressed.
The compromise that this filter makes is between suppressing noise and retaining smaller details in the image. Those characteristics that induce coefficients higher than the noise are easily retained, as they are not degraded. Characteristics that induce coefficients of the same magnitude as noise are suppressed because they are not distinct from noise.

Scope and Objectives
Imaging techniques for the human body have grown increasingly popular in recent times, with the development of equipment in this field, as well as with the development of hardware and software systems in the field, capable of managing the volumes of data involved in medical imaging. This has allowed today more and more advanced methods in the field of MR (Magnetic Resonance) imaging, CT (Computer Tomography) or PET (Positron Emission Tomography). All of these methods require advanced image processing in order to extract useful diagnostic information from the available volume of data. For example, 3D rendering techniques allow better visualization of three-dimensional volumes, feature extraction techniques from images allow better assisted diagnosis, image segmentation techniques allow the separation of different areas or tissues of interest, and noise elimination techniques from images allow the improvement of image quality [55]. Some of these image processing techniques can be used in a processing chain. For example, in order to segment certain parts of the human body and obtain good segmentation results, the noise from the images must be removed [56][57][58].
In this article we propose and characterize by simulations a novel algorithm used to reduce the white Gaussian additive noise in the images. The algorithm exploits both types of dependencies between wavelet coefficients. We characterized the algorithm through simulations. The proposed algorithm exploits the intra-scale dependence of wavelet coefficients by determining the local dispersion of the signal, and the inter-scale dependence by determining the contours based on the hierarchical correlation map. The results obtained by applying this algorithm show that the reduction of noise in images is achieved when the quality of the contours is maintained or even improved. At the same time, by choosing a parameter, it is possible to establish the compromise between the amount by which the noise is reduced and the quality of preserving the details in the image [59].
This article is divided in four sections, which include Introduction, Materials and Methods, Results and Discussions and Conclusions. In the Introduction, we introduced the subject of our research and presented the state of the art related to the topic. The second part is Materials and Methods where all the steps of the filtering technique are presented: Adaptive spatial filtering based on inter-band and intra-band correlation of wavelet coefficients and Adaptive spatial soft truncation based on the hierarchical correlation map and with the establishment of the threshold value based on local statistics. In the third part we presented our Results and Discussions in which we conducted studies both on usual images and also on CT images. Finally, in the last part, we presented the Conclusions of the obtained results.

Materials and Methods
This study was carried out using documentation based on papers and articles presented in journals or at conferences, scientific literature related to adaptive spatial filtering and wavelet domain topics, scientific on-line databases such as Google Scholar, Google Academic, MDPI, Scopus, Science Direct and other web pages or research platforms. Table 1 presents the types of images considered as inputs for the filter, together with their peak signal to noise ratios (PSNR) and the noise dispersion (sigma):  This article makes use and conducts an in-depth review on a great number of papers, official documents, information related to de-noising methods in the wavelet domain and also findings from other research institutes relevant to our topic [60][61][62][63], such as the National Research and Development Institute for Cryogenic and Isotopic Technologies ICSI, Ramnicu Valcea, Romania, Automotive Department.

Theoretical Background
Based on both inter-band and intra-band correlation, J. Liu and P. Moulin [64] proposed an algorithm for reducing white Gaussian additive noise in images, whose results are compared by the authors with those obtained by the following methods: -Hard truncation with universal threshold; -A posteriori maximum probability estimator (MAP), also proposed by the authors in [65]; -A maximum probability estimator modeled by hidden Markov chains, proposed by J. Romberg in [66]; -Adaptive Wiener filtering in the wavelet domain, proposed in [67,68].
The authors compared their proposed method with the maximum probability estimator and adaptive Wiener filtering in the wavelet domain, showing that, depending on the processed image, significant improvements were obtained both in terms of MSE and visually, appreciating them (in terms of MSE) between 8% and 26% for the highest probability estimator, and between 1% and 7% for the Wiener filter. The authors also compared the proposed method with the spatially adaptive method proposed by Chang [69] and presented in the previous paragraph, in which case the results are lower than this (by almost 15% in MSE terms).
In this method, the inter-band correlation is used to obtain a map containing the position of those wavelet coefficients that have a high informational content (contours and texture), called the significance map. The wavelet coefficients thus classified are processed by soft truncation using threshold values adapted to each sub-band, obtained by modeling their distribution by a Laplace distribution.
The intra-band correlation is used to estimate the local signal scatter in the case of coefficients that do not belong to the significance map, which are processed with a maximum probability estimator.
Considering a three-level orthogonal wavelet decomposition, the algorithm proposed by J. Liu and P. Moulin [70] can be described as follows: -As in many other noise reduction methods using dyadic wavelet transformation, the lowest resolution sub-band is not processed (represented in blue in Figure 1); -Each of the three low-resolution detail sub-bands (corresponding to wavelet decomposition level 3, red in Figure 1) is processed by soft truncation; the threshold value for each sub-band is determined by modeling the distribution of its coefficients by a Laplace distribution of zero mean and dispersion: where: is the Gaussian noise dispersion, M being the number of coefficients in the considered sub-band and w the wavelet coefficients of the noise-degraded image.
In this case, the threshold value corresponding to sub-band j will be: in the case of the other detail sub-bands (corresponding to wavelet decomposition levels 2 and 1 and represented by the colors green and white in Figure 1), the significance map is drawn up based on the absolute value of the parent coefficient corresponding to each coefficient belonging to these sub-bands (in Figure 1, the relation descending coefficient-parent coefficient is also represented); thus, if the value of the parent coefficient is higher than a threshold value T, the coefficient considered is classified as significant, otherwise it is classified as insignificant; the two classes of coefficients, having different informational significance, have different statistical properties, therefore they are processed differently: o the significant classified coefficients are soft truncated, the threshold value used is determined based on the modeling of their distribution with a Laplace distribution of zero mean and dispersion: where: sig representing the significance map, and M s the number of coefficients in this map. o insignificant classified coefficients have small values and represent smooth regions; for each such coefficient the dispersion is estimated using a window of dimensions 5 × 5, but from which the significant coefficients are eliminated; ifσ 2 i represents the thus estimated dispersion of the signal, then the a posteriori maximum probability estimator will be: As it can be seen, the algorithm, just like other algorithms that determine the contours, introduces as a parameter a T-threshold value.
As it can be seen, the algorithm, just like other algorithms that determine the contours, introduces as a parameter a T-threshold value.  Wavelet decomposition (blue-lowest resolution sub-band is not processed; red-low resolution sub-bands, corresponding to wavelet decomposition level 3, are processed by soft truncation; green, white-detail sub-bands, corresponding to wavelet decomposition level 2 and 1, are processed based on the significance map).

Algorithm Proposed in This Work
In the article we propose a spatially adaptive soft Truncation algorithm based on the hierarchical correlation map and with the establishment of the threshold value based on local statistics. The proposed algorithm exploits the intra-scale dependence of the wavelet coefficients by determining the local dispersion of the signal, and the inter-scale dependence is exploited by determining the contours based on the hierarchical correlation map. Results obtained by applying this algorithm show that the reduction of noise in images is achieved and the quality of the contours is maintained or even improved. At the same time, by choosing a parameter, it is possible to establish the compromise between the amount by which the noise is reduced and the quality of preserving the details in the image.
The filtering technique proposed by us is based on the fact that the abrupt variations of the pixel intensity in the image determine wavelet coefficients that have high absolute values in several scales, while the noise induces wavelet coefficients that decrease in absolute value with increasing scale.
We define the direct spatial correlation of the undecided wavelet transformation on several adjacent scales to accurately determine the locations of contours or other important features in the image: where l is the number of scales considered in direct multiplication, M being the total number of scales, N is the number of samples in the signal, and w( j, n) are the wavelet coefficients of the degraded noise signal corresponding to scale j. The absence of contours or other significant features in a particular region will allow background noise to be removed from that region. The direct spatial correlation on several adjacent scales of the data will lead to the precise location of the important contours in the image.
If we consider an undecided wavelet decomposition on three levels (so we have three scales), then the direct spatial correlation will be: Data representing Corr 2 (1, n) are rescaled as a function of w(1, n) .
Appl. Sci. 2020, 10, 5693 9 of 29 The most important contours are identified by comparing the absolute values of Corr 2 (1, n) and w(1, n). It is assumed that the position n belongs to a contour if: The w(1, n) coefficient corresponding to position n, identified as belonging to the contour, is stored in a new matrix w new (1, n), the contour position is also stored and the value of w(1, n) is cancelled in the initial matrix. Finally, all determined contours are extracted from Corr 2 (1, n) and w(1, n) by cancelling their values in the locations identified as belonging to contours. We thus obtain two new data sets Corr 2 (1, n) and w (1, n) which form the basis of the next round of contour determination. We proceed further to a new data scaling Corr 2 (1, n) and a new contour identification.
This procedure, consisting of normalization, identification of contour locations, retention and cancellation, is performed iteratively several times until the power of the remaining coefficients in w (1, n) is less than a given value.
All contour information contained in the original data that is extracted from w(1, n) during the iterative process is stored in a new w new (1, n) vector. By replacing w(1, n) with w new (1, n) we will have the spatially filtered version of the first scale of the undecided wavelet transformation, a scale that actually contains most of the noise. This type of filter can be seen as a spatially dependent filter.
The block diagram of the proposed algorithm is presented further in Figure 2 and is discussed in more detail in the subsequent paragraphs. The whole procedure can be described with the following algorithm: A. For each scale, from = 1. . . ; B. We determine the data set ( , ), with the following itterative process: B1. We compute ( , ) and ( , ) B2. We then compute the correlation and wavelet factors: B3. We perform the scaling: B4. We identify the contour location by verifying if the following condition is met: The whole procedure can be described with the following algorithm: A. For each scale, from l = 1 . . . M; B.
We determine the data set w new (l, n), with the following itterative process:
We then compute the correlation and wavelet factors: PW(l) = n w 2 (l, n) B3. We perform the scaling: B4.
We identify the contour location by verifying if the following condition is met: Additionally, for each location identified as a contour, we execute: We then test the stopping condition of the iterative process: In case of non-fulfillment of this condition the process is resumed from step B2, and in case of fulfillment of the condition it is passed to the next scale (step B).
At the end of this procedure, in w new (m, n), we have the filtered values of the wavelet coefficients, corresponding to the undecimal wavelet transform of the estimated signal.
We tested the proposed procedure for noise reduction from images obtained from medical applications, using a three-level wavelet development (l = 2), which corresponds to 8 sub-bands. The results obtained by this method were compared with results obtained by Wiener filtration. It was thus observed that the images filtered by the proposed method are affected by less noise, better contoured and do not show artifacts generated by the Gibbs effect or other artifacts and have a better visual quality compared to images filtered with Wiener filter.
Maintaining the idea of punctually establishing the threshold value and exploiting inter-band dependence by determining a map of accounts using the concept of hierarchical, we propose a new spatial adaptive filter.
The proposed algorithm has 3 stages. In the first stage, the threshold values for each coefficient, necessary for the soft truncation of the wavelet coefficients, are determined based on the local statistics. The difference from the method set out in the previous paragraph is primarily in simplifying the determination of local signal dispersion. For this we discard the selection of the coefficients used to determine the signal dispersion based on the random variable z, taking into account all wavelet coefficients contained in a given window, but which has much smaller dimensions and we discard the use of the parent coefficient used to evaluate local dispersion.
In this case, the local dispersion corresponding to the w(i, j) coefficient will be given by: where: B(i, j) being the set of indices of the coefficients in the current window, and M is the number of coefficients contained in the window.

of 29
The current threshold value will be given by: Furthermore, in this stage, the contour map is determined based on the hierarchical correlation corresponding to the degraded image.
In the second stage, all wavelet coefficients of the detail sub bands are subjected to a soft truncation operation, determining the coefficients: η s being the threshold soft truncation operator.
In the third stage, a correction is made on the w s (i, j) coefficients obtained in the second stage, depending on the contour map. Thus, the coefficients that were established as belonging to a contour are estimated by the arithmetic mean between the initial value and the value obtained in the second processing stage, the rest of the coefficients remaining unchanged: We call this correction switched correction based on the hierarchical correlation map.
Another proposed possibility of correction on the w s (i, j) coefficients according to the hierarchical correlation map, this time a continuous one, is given by: where ρ(k, l) represents the normed value of the hierarchical correlation corresponding to the coordinate coefficient (k, l). Table 2 compares the results obtained in the case of processing the details sub bands corresponding to the first scale of a wavelet transformation by the method described by Chang and Vetterli [71] (adaptive soft truncation with contextual modeling), which was presented in the previous paragraph, and applying the first step of the proposed method.

Results and Discussions
As it can be seen from the data presented in Table 3, even in the case of applying only the first two steps of the proposed algorithm, better results are obtained both in terms of PSNR and in terms of contour conservation, given by the C coefficient. This coefficient is computed based on the reduction factor and represents how much of the contours were preserved from the original image.
Regarding the proposed method, it can be found that, as the window size increases, the results are weaker both in terms of PSNR and in terms of the C coefficient. This can be justified by the fact that those signal-induced wavelet coefficients are grouped, and with the increase in the window, their influence on the local dispersion decreases, approaching the noise dispersion. A first parameter of the proposed algorithm is the length of the window. The effect of this parameter, as well as of the number of wavelet decomposition levels used on the image quality obtained after the second processing step is presented in Table 4, compared to the results obtained by soft truncation with optimal threshold values. A third parameter that influences the performance of the algorithm is the threshold value used to obtain the contour map from the hierarchical correlation map, whose effect on the whole filtering process, is presented in Appendix A.
It can be observed that in all the considered cases, better results are obtained in PSNR terms than in the case of soft truncation with optimal threshold values, in many cases obtaining even higher C coefficients than the initial ones, hence an improvement of the contours.
Considering an ideal contour map, determined based on the image which is not degraded with noise, the simulations show (Appendix A) that, by the proposed method, an improvement of the contours can be obtained, expressed by higher values of the C coefficient, in the conditions in which the noise reduction is not affected, which is expressed by constant PSNR values. In the real case, only the noise-degraded image is available, the contour map being determined based on it. Table 2. Comparative results obtained by adaptive soft truncation with contextual modeling and applying the first two steps of our proposed method, in the case of processing only the details sub bands corresponding to the first scale of a wavelet transformation. The side of the analysis window is 2L + 1. As it can be seen in Figure 3, false contour points appear in this case, resulting in a smaller amount of noise removed. Therefore, on the hierarchical correlation map, before determining the contours, we proceed to process it with a median-hybrid filter, described in [72][73][74]. As it can be seen in Figure 3 and from the data in the Appendix A, this filtering operation leads to images characterized by PSNR values closer to those of the images obtained after the first two processing steps of the proposed algorithm, in the conditions under which an improvement of the contours takes place.

By Adaptive Soft Truncation with
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 29 As it can be seen in Figure 3, false contour points appear in this case, resulting in a smaller amount of noise removed. Therefore, on the hierarchical correlation map, before determining the contours, we proceed to process it with a median-hybrid filter, described in [72][73][74]. As it can be seen in Figure 3 and from the data in the Appendix A, this filtering operation leads to images characterized by PSNR values closer to those of the images obtained after the first two processing steps of the proposed algorithm, in the conditions under which an improvement of the contours takes place.  Figure 4.
The analysis of the data presented in the Appendix A shows that the effect of the threshold value used to obtain the contour map from the hierarchical correlation map does not depend on the processed image (as in the case of optimal threshold values in the case of soft truncation), but only on dispersion of the disturbing noise. Thus, for a noise dispersion = 0.05, the optimal value is between 0.1 and 0.15, and for = 0.1, between 0.15 and 0.175.  Table 3. Dependence of the results obtained in the first stage of the proposed algorithm, as a function of Table 2 (L + 1). In all cases, better results are obtained when using local statistics compared to optimal threshold, for both PSNR and C (improvement of the contours).  The proposed algorithm was tested on different images with dimensions of 256 × 256 pixels and 512 × 512 pixels, under white Gaussian additive dispersion noise degradation conditions: σ = 0.05 and σ = 0.1. The results are presented in Figure 4.
The analysis of the data presented in the Appendix A shows that the effect of the threshold value used to obtain the contour map from the hierarchical correlation map does not depend on the processed image (as in the case of optimal threshold values in the case of soft truncation), but only on dispersion of the disturbing noise. Thus, for a noise dispersion σ = 0.05, the optimal value is between 0.1 and 0.15, and for σ = 0.1, between 0.15 and 0.175. The proposed algorithm was tested on different images with dimensions of 256 × 256 pixels and 512 × 512 pixels, under white Gaussian additive dispersion noise degradation conditions: = 0.05 and = 0.1. The results are presented in Figure 4.
The analysis of the data presented in the Appendix A shows that the effect of the threshold value used to obtain the contour map from the hierarchical correlation map does not depend on the processed image (as in the case of optimal threshold values in the case of soft truncation), but only on dispersion of the disturbing noise. Thus, for a noise dispersion = 0.05, the optimal value is between 0.1 and 0.15, and for = 0.1, between 0.15 and 0.175.  The interest in such filters is determined by the fact that much of the existing information in the images is provided by contours, and most of the noise removal methods cause them to fade. Such filters can be designed starting from the finding that near the contours the noise is less noticeable than in the smooth regions, thus being possible to change the filter parameters depending on the level of activity of the processed region. In this case, at the cost of removing a smaller amount of noise in the regions containing contours, they can be better preserved. There is thus a compromise between the amount of noise removed and the quality of preserving the contours in the image.
In Table 4, we exemplify the results obtained for the filtration method proposed in two working variants. The first variant proceeds to segment the LL (low level) sub-band of the lowest resolution from a dyadic wavelet decomposition, and depending on the result of the segmentation operation, it chooses a threshold value necessary for the soft truncation of the wavelet coefficients. The method does not require a large computational effort, but the performance obtained by using it is not special. The second method also uses soft truncation of the coefficients, but this time the threshold value is calculated for each coefficient based on local statistics in the wavelet field. On the wavelet coefficients thus obtained, a correction is made according to a contour map determined on the basis of a hierarchical correlation map.
Two variants of correction are proposed: a switched correction, in which case either the truncated coefficient is chosen, or the average between the truncated and the unprocessed coefficient, depending on whether or not the respective coefficient belongs to a contour and a continuous correction, in which case a weighted sum of the two wavelet coefficients is performed, the values of the weights being calculated based on the hierarchical correlation coefficient. In the second case, images are obtained with higher PSNR values, especially in conditions of higher noise level, but which have more blurred contours. We approached this algorithm in invariant translation form in two ways: (A) By integrally applying the proposed algorithm in translation invariant form, mediating the results obtained by applying all three processing steps on the cyclic displacements of the input image; (B) By partially applying the proposed algorithm (only steps 1 and 2) in the translation invariant form and performing the correction according to the contour map only on the result thus obtained.
It can be observed from the data presented in Table 5, that the invariant approach to translation of the proposed algorithm leads to good results both in terms of PSNR and in terms of contour conservation. Only in one of the cases presented in Table 5 does the initial image have a higher C coefficient than the processed image, otherwise the processed images have a higher C coefficient. This means that with the removal of the noise, an improvement of the contours is obtained.
It is also found that variant A leads to better overall results than in variant B, but the calculation effort is higher in this case. The proposed algorithm leads both in the invariant approach to translation and in the simple one to good results both in terms of PSNR and C coefficient, as well as visually, as it can be seen in Figure 5 (Comparative Performance in Standard Images). Table 6 compares the two correction methods based on the hierarchical correlation map. It can be seen that in the case of continuous type correction, better results are obtained in terms of PSNR than in the case of switched type correction, especially for higher noise levels. However, this is done at the cost of poor preservation of the contours in the image, which is highlighted by lower values of the C coefficient.
Currently, imaging is present in all branches of science, and the acquisition of images that are not disturbed by external factors arouses interest in the field of scientific research. Digital images unaffected by disturbing elements, such as noise, contrast and uneven lighting, allow a correct understanding of its content.
Following the research of the specialized literature and the implemented applications, in this article we will confer a benefit in the extraction of useful information from medical images.
The objectives of the research are to model, process and extract useful information from digital medical images using image segmentation and fusion methods. These techniques generate higher quality images than the original ones.   The objectives of the research are to model, process and extract useful information from digital medical images using image segmentation and fusion methods. These techniques generate higher quality images than the original ones.

Image Initial Switched Type Correction Continuous Type Correction PSNR (dB) C (%) PSNR (dB) C (%) PSNR (dB) C (%)
To analyze the useful informational content in medical images, using hemodynamic modeling and imaging processing to improve diagnosis, we propose and will achieve the following objectives:  To analyze the useful informational content in medical images, using hemodynamic modeling and imaging processing to improve diagnosis, we propose and will achieve the following objectives: • Improving image quality by selectively eliminating disturbing information, such as noise, and eliminating other defects caused by the acquisition device by using adaptive filters.

•
Highlighting areas of interest by adjusting light intensity and contrast and accentuating contours and textures.

•
Improving the ability to detect edges by using the image fusion method based on the wavelet algorithm.

•
Developing new methods for evaluating segmentation.

•
Evaluating the performance of edge detection methods using operators that analyze the structural similarity of images.

•
The use of high-performance algorithms of Artificial Intelligence to perform sub-pixelation calculations and the transition from macroscopic analysis currently used to microscopic analysis for the analysis of tumors generated by cancer by analyzing medical images.
The innovative nature of the research will be based on the development of a high-performance algorithm to improve the results of medical image filtering using Gaussian functions, multidimensional wavelet decompositions and image fusion. We will also implement a high-performance algorithm by using a matched filtering that correlates a known signal, or a pattern, with an unknown signal to detect the presence of the pattern in the unknown signal, in order to improve the analyzed images. We will focus on using the method of merging wavelet images and neural networks, as it exceeds the limits of classical edge detection methods and provides a much more accurate map of decompositions at the sub-pixel level. Last but not least, we will develop an algorithm for the automatic selection of the detection threshold for the classification of brain medical images obtained from CT and MRI, using clustering algorithms.
From the research carried out in this article we propose an algorithm that responds to the first three points mentioned above.
The proposed image enhancement method using spatial filtering algorithm can be successfully applied to low contrast CT (Computer Tomography) images. Figure 6 shows the results obtained for such images. Evaluating the performance of edge detection methods using operators that analyze the structural similarity of images. • The use of high-performance algorithms of Artificial Intelligence to perform sub-pixelation calculations and the transition from macroscopic analysis currently used to microscopic analysis for the analysis of tumors generated by cancer by analyzing medical images.
The innovative nature of the research will be based on the development of a high-performance algorithm to improve the results of medical image filtering using Gaussian functions, multidimensional wavelet decompositions and image fusion. We will also implement a highperformance algorithm by using a matched filtering that correlates a known signal, or a pattern, with an unknown signal to detect the presence of the pattern in the unknown signal, in order to improve the analyzed images. We will focus on using the method of merging wavelet images and neural networks, as it exceeds the limits of classical edge detection methods and provides a much more accurate map of decompositions at the sub-pixel level. Last but not least, we will develop an algorithm for the automatic selection of the detection threshold for the classification of brain medical images obtained from CT and MRI, using clustering algorithms.
From the research carried out in this article we propose an algorithm that responds to the first three points mentioned above.
The proposed image enhancement method using spatial filtering algorithm can be successfully applied to low contrast CT (Computer Tomography) images. Figure 6 shows the results obtained for such images. As it can be seen, the proposed method provides robust results in eliminating noise from images [75], as well as improving image contrast and contour detection. Obtaining images in which the contours are best preserved is done primarily in conditions of increasing computational complexity, with a decrease in the amount of noise removed from the image taking place. Figures 7-11 show the noise elimination results obtained in the case of real CT images, using Daubechies scale filter order 2, order 4 and Haar wavelet transform [76]: As it can be seen, the proposed method provides robust results in eliminating noise from images [75], as well as improving image contrast and contour detection. Obtaining images in which the contours are best preserved is done primarily in conditions of increasing computational complexity, with a decrease in the amount of noise removed from the image taking place.          . Noise reduction from a CT image, using Daubechies scale filter order 4 (db4) [76]. Figure 9. Noise reduction from a CT image, using Daubechies scale filter order 4 (db4) [76]. Appl. Sci. 2020, 10, x FOR PEER REVIEW 19 of 29 Figure 10. Noise reduction from a CT image, using Haar wavelet [76]. Figure 11. Noise reduction from a CT image, using Daubechies scale filter order 4 (db4) [76].
As it can be seen, the proposed method and the determined operators offer a better efficiency in terms of eliminating noise and preserving the contours of an image, allowing the implementation of the proposed method on existing hardware structures.

Conclusions
The interest in wavelet filters is determined by the fact that much of the existing information in the images is provided by contours, and most of the noise removal methods cause them to fade. Such filters can be designed starting from the finding that near the contours the noise is less noticeable than in the smooth regions, thus being possible to change the filter parameters depending on the level of activity of the processed region. In this case, at the cost of removing a smaller amount of noise in the regions containing contours, they can be better preserved. There is a trade-off between the amount of noise removed and the quality of preserving the contours in the image.
The article classifies the methods of selective filtering in the wavelet domain according to the statistical modeling of the wavelet coefficients on which these methods are based. Thus, selective filtering methods in the wavelet domain are divided into three groups: methods based on the intraband dependence of the wavelet coefficients, methods based on the inter-band dependence of the wavelet coefficients and methods based on both the intra-band dependence of the wavelet coefficients and the inter-band one.
The proposed method uses soft truncation of coefficients, where the threshold value is calculated for each coefficient separately based on local statistics in the wavelet domain. On the wavelet coefficients thus obtained, a correction is made according to a contour map determined on the basis of a hierarchical correlation map.
The proposed algorithm uses noise dispersion as a parameter, but the switched version needs another parameter, namely the threshold value to obtain the contour map from the hierarchical correlation map. Figure 10. Noise reduction from a CT image, using Haar wavelet [76].
As it can be seen, the proposed method and the determined operators offer a better efficiency in terms of eliminating noise and preserving the contours of an image, allowing the implementation of the proposed method on existing hardware structures.

Conclusions
The interest in wavelet filters is determined by the fact that much of the existing information in the images is provided by contours, and most of the noise removal methods cause them to fade. Such filters can be designed starting from the finding that near the contours the noise is less noticeable than in the smooth regions, thus being possible to change the filter parameters depending on the level of activity of the processed region. In this case, at the cost of removing a smaller amount of noise in the regions containing contours, they can be better preserved. There is a trade-off between the amount of noise removed and the quality of preserving the contours in the image.
The article classifies the methods of selective filtering in the wavelet domain according to the statistical modeling of the wavelet coefficients on which these methods are based. Thus, selective filtering methods in the wavelet domain are divided into three groups: methods based on the intraband dependence of the wavelet coefficients, methods based on the inter-band dependence of the wavelet coefficients and methods based on both the intra-band dependence of the wavelet coefficients and the inter-band one.
The proposed method uses soft truncation of coefficients, where the threshold value is calculated for each coefficient separately based on local statistics in the wavelet domain. On the wavelet coefficients thus obtained, a correction is made according to a contour map determined on the basis of a hierarchical correlation map.
The proposed algorithm uses noise dispersion as a parameter, but the switched version needs another parameter, namely the threshold value to obtain the contour map from the hierarchical correlation map. Figure 11. Noise reduction from a CT image, using Daubechies scale filter order 4 (db4) [76].
As it can be seen, the proposed method and the determined operators offer a better efficiency in terms of eliminating noise and preserving the contours of an image, allowing the implementation of the proposed method on existing hardware structures.

Conclusions
The interest in wavelet filters is determined by the fact that much of the existing information in the images is provided by contours, and most of the noise removal methods cause them to fade. Such filters can be designed starting from the finding that near the contours the noise is less noticeable than in the smooth regions, thus being possible to change the filter parameters depending on the level of activity of the processed region. In this case, at the cost of removing a smaller amount of noise in the regions containing contours, they can be better preserved. There is a trade-off between the amount of noise removed and the quality of preserving the contours in the image.
The article classifies the methods of selective filtering in the wavelet domain according to the statistical modeling of the wavelet coefficients on which these methods are based. Thus, selective filtering methods in the wavelet domain are divided into three groups: methods based on the intra-band dependence of the wavelet coefficients, methods based on the inter-band dependence of the wavelet coefficients and methods based on both the intra-band dependence of the wavelet coefficients and the inter-band one.
The proposed method uses soft truncation of coefficients, where the threshold value is calculated for each coefficient separately based on local statistics in the wavelet domain. On the wavelet coefficients thus obtained, a correction is made according to a contour map determined on the basis of a hierarchical correlation map.
The proposed algorithm uses noise dispersion as a parameter, but the switched version needs another parameter, namely the threshold value to obtain the contour map from the hierarchical correlation map.
Unlike the methods presented in the literature, which in the case of very rich in details images leads to weaker results than soft truncation with optimal threshold values, the proposed method, by soft truncation of wavelet coefficients with threshold value adapted based on local statistics and correction based on the hierarchical correlation map, always leads to better results both in terms of PSNR and in terms of contour conservation.
Filtration methods mentioned in the literature prove to be very effective in terms of the amount of noise removed from the image, generally leading to high values of PSNR, but this is done at the cost of significant contour degradation. The proposed method leads to close values, although lower in terms of PSNR than other methods, but ensures a much better preservation of contours.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Dependence on the threshold value used to obtain the contour map from the hierarchical correlation map, in the case of the proposed method. Results obtained by using the ideal, real and filtered hierarchical correlation map are presented.