Adaptive Gradient-Based and Anisotropic Diffusion Equation Filtering Algorithm for Microscopic Image Preprocessing ()
1. Introduction
1.1. Quality of Microscopic Images
The limitation of the quality of microscopic images can be attributed to image intensity levels, image resolution, image pixel size and dynamic ranges. The former two factors are related to the optical system and the latter two are related to the detector system.
1.2. Noise in Microscopic Images
Due to the influence of image quantization, adding quantization noises exist in microscopic images and the whole image background shows nonuniform distribution. In Figure 1, the gray level distributions of three horizontal areas were computed in three different positions Y = 30, 1215 and 1900. It is evident that the gray level distribution in image is not constant and nonuniform exists, which is smaller in background regions and larger in boundary regions. Another characteristic of microscopic image is that objects are always in motion. Due to the amplification effect of microscope, the image is sensitive to focal length and usually the images appear blurry near the boundaries and halo-type boundaries appear. In this case, the precise detection of boundary becomes very difficult. Figure 2 demonstrates three typical morphologies of cells in the image, which all have blurry and halotype boundaries. Because microscopic images are fluid images, apart from the objects, there exists some other objects like bubbles, sediments etc.
1.3. Microscopic Image Preprocessing Research Status
Given the complexity of factors introducing noise and intensity heterogeneity, there are two approaches for improving the quality of microscopic images. First, one could focus on improving specimen preparation and imaging conditions, for example improved fluorescent dyes. Second, one could attempt to correct pixel intensities after image acquisition, as it is the case for image restoration. Apart from some traditional filtering algorithm such as self-adaptive median filtering algorithm [1,2], hybrid filtering algorithm [3], anisotropic filtering algorithm [4], feature fusion filtering algorithm [5-7], wavelet based filtering algorithm [8-11] etc., some specific filtering and restoration algorithms were also proposed such as the empirical correction methods for intensity loss [12], constant thresholding [13], iterative correction methods [14], 2D histogram or estimations of intensity decay function [15]. Although these methods work well for certain category of microscopic image, most of them require the image meet some basic requirements such as the photo bleaching is a spatially homogeneous in a lateral plane, which greatly limit the universality and availability.
1.4. Quality Degration Model of Microscopic Images
Normally in image denoising techniques, noise can be
ideally divided into adding noise and multiplicative noise. Suppose the degraded image as g, the ideal image as f, and the noise as n. The adding model and the multiplicative model can be described in formula (1) and (2). Background noise in microscopic images can normally be regarded as adding Gaussian noise and low magnification lens tiny particles as adding pulse noise.
(1)
(2)
Besides noises, there are some other quality degradation factors like inappropriate focusing, uneven illumination etc. For better cell segmentation and recognition, the system should design a good preprocessing algorithm to overcome noise and other image quality degration factors.
The boundary of cell is very important for quantitative analysis of microscopic images. As described above, the qualities of microscopic image are always degraded by adding noise and halo-type boundaries. So it is very important to design an appropriate filtering algorithm that can filter out the noise and will not further degrade the image quality of boundaries. In this paper, an adaptive gradient-base and anisotropic diffusion equation filtering algorithm was proposed to deal with this problem.
2. Proposed Algorithm
In order to accelerate the filtering speed, a judgment should be made to each pixel in image and then the filtering manipulation is implemented on pixels that are judged as noises. For non-noise pixel, no filtering manipulation is implemented.
2.1. Gradient Based Background Noise Detection
Based on the description in introduction section, it is known that in background area, the gray scale changes less while in boundary area or noise it changes more. Therefore, we can compute the variation of gray scale in a sliding window while it slides through the whole image to judge the property of each pixel, boundary or noise. Then the filtering is implemented on possible noise candidate pixel. Suppose the gradient of each direction as, i denotes the all eight directions. The judgment rule is described as following equations:
(3)
(4)
After the judgment is made for all pixels, the filtering manipulation is implemented on the pixels that are judged as boundary or noise.
2.2. Anisotropic Filtering Algorithm
According to the concept of divergence field, anisotropic diffusion equation can be described in the following equation:
(5)
In Equation (5), is divergence operator, is gradient operator, I is the function of, is spatial scale function, namely the diffusion coefficient, which is a nonnegative monotonic decreasing function, Ii is the derivative of I with respect to time t, t is the time of thermal diffusion.
The selection of diffusion coefficient will directly influence the filtering effect and in general, diffusion coefficient is given a value of the norm of a vector, which has the following attributes:
1) In interior area.
2) In the boundary area, , is the unit vector of the gradient of point in the boundary area. is the strength difference at each side of the boundary.
We can determine as
(6)
The classical can be described as
(7)
Or
(8)
Here, k is the controlment coefficient of diffusion strength, x is the gradient of the diffusion point. According to the definition of diffusion coefficient, in different directions different coefficients are adopted. In addition, the monotonic decreasing function was taken in different directions as the diffusion coefficient. In background area or interior area of image, the gray level values are similar and gradient is very small, so the diffusion coefficient is large to implement smoothing. On the contrary, in the area of boundary, gray level value changes much and accordingly the gradient increases, diffusion coefficient is very small to preserve the boundary information. Here we consider the processing in eight directions. Figure 3 demonstrate the flow chart of the proposed algorithm.
The algorithm is achieved using iterations that are described in Table 1 as the discretization of the thermal diffusion process. The detailed algorithms are as follows.
3. Experimental Results
The filtering was implemented using the proposed algorithm of gradient-based background noise detection and anisotropic filtering. The diffusion coefficient is set as 10, iteration step size set as 0.21, iteration time as 5, 10, 15, 20 and thresholding value for boundary and noise as 30.
Table 1. The discretization process of the proposed algorithm.
Figure 3. Flow chart of proposed algorithm.
A sample microscopic image was used, which incorporates the typical features of microscopic images with adding background noise, pulse noise and tiny particle noise. The results of the algorithm are compared with that of median filtering. It is very clear that the filtering algorithm proposed above can obtain satisfactory results in removing background noise as well as preserving object boundaries. In addition, by using noise detection algorithm, the computational speed improved by 3 times compared with using anisotropic filtering algorithm alone.
Figure 4 demonstrates the comparison of experimental results between the proposed algorithm and the median filter. The top row shows the filtering results of the proposed algorithm using different iteration times. The bottom row shows the filtering results of traditional median filtering algorithm using different padding sizes. It is evident that our algorithm obtained better performance in removing noises while preserving boundary details.
To demonstrate the advantages of the proposed algorithm to other filtering algorithm, the performance of proposed algorithm (Aniso) was compared with other four widely used filtering algorithms in medical image processing including homomorphic filtering algorithm (Homo) [16], DCT-based filtering algorithm [17], wavelet based filtering algorithm (Wavelet) [18], isotropic diffusion based filtering algorithm (Iso) [19]. The results in Figure 5 demonstrate that the proposed algorithm outperforms other algorithms in filtering out noise while preserving boundaries.
Figure 6 demonstrates the surface map of the original image and the filtered image using proposed algorithm and other four algorithms. Surface map draws a wire frame mesh with color determined by the intensity value in the image and the color is proportional to the surface height. It is evident that the proposed anisotropic diffusion filtering algorithm can obtain better filtering performance in removing noise and preserving boundary