The Improvement of the Region Homogeneity Measure Operator and Its Noise Disposed

Based on the analysis of the characteristic of the image texture, we extend the region homogeneity measure operator (RHMO) to a more reasonable case, thus avoiding some omission situation. Moreover, we compare the RHMO with the convolution function to find out the more rational range of region homogeneity measure (RHM). At the end of this paper, we use the improved RHM method for some images in a new Partial Differential Equation (PDE); the results indicate that the improved concept has achieved a good result.


Introduction
In the recent thirty years, various digital image processing methods based on PDE [1][2][3] have been put forward and studied widely.These methods will overcome some defects of the classical image processing in some aspects, so they have attracted the attention of numerous researchers and have been applied to image processing widely.
In the field of classical image processing [4], Fourier method can be regarded as the main mathematical method.In the algorithm of the classical image processing, the smooth method has a serious defect: while doing the image smooth processing, it will also do the same processing on the edge of the image, which is a main characteristic of an image.So after smoothness, some characteristics of edge have been eliminated.Generally speaking, the traditional methods in image processing are difficult to achieve the ideal effect.But the digital image processing method which is based on PDE has overcome the difficulty described above to some degree.
Different from the traditional method, the RHMO [5] does not use the gray gradient which is defined among several pixel points.Instead, the gradient is replaced by the average gray gradient which includes all the pixel points in the neighborhood.This simple replacement can achieve a similar effect in image processing as the traditional method.RHM is defined as follows.
In (1), the noise's characteristics of single pixel and region have been considered.As to the single pixel, we are very approximate to regard it as a solid circle (as Figure 1).While dividing the area, if the diameter ℓ passes through  which is the center of a circle, the value of   is zero.On the other hand, because the distribution of the noise is random, in all whirling directions of ℓ, the distribution and the intensity of the noise are all the same on probability (as Figure 2).The influence of the noise is counteracted.So the value of    will not change obviously-while its value can reflect conformity of area's characteristic.At the same time, the value can reflect the possibility of the image's edge.In the subsequent part of this paper, we improve the RHMO by the analysis of the characteristic about the noise and the texture, thus preventing some situations that RHMO will ignore in the processing of the image.On the other hand, we provide another improvement of RHMO through the analysis of the convolution term and the RHMO term and that lead to a more reasonable range of the value of RHMO.

The Improvement of the Model
2.1.The New Division for the Region.During the process of calculating RHM, choosing the semidiameter  depends on the characteristic of the image.If the distribution of the edge is dispersed, the value of  can be a bigger one; for the opposite situation (intensive situation), the value of  should be a smaller one.But the treatment may be insensitive to some particular texture.So we may regard the corresponding texture as the noise to deal with.We consider such an area (as Figures 3 and 4 show, the distribution of texture has axial symmetry property) according to the definition of RHMO at this moment; we regard the texture as the noise to deal with.This kind of situation may appear, for example, in the picture of fingerprint, and an area of 20×20 pixels may distribute 2 or 3 strips of fingerprint.Of course, we can break the symmetry of Figures 3 and 4 by dwindling the value of ; thus the texture of the image can be searched out.Because we cannot distinguish the noise and true texture in a small region, the value of  cannot be too small (especially in this situation, the semidiameter of the noise is bigger than the normal size).In fact, we can use a new definition about the RHM which does not calculate the difference between two areas divided by the diameter.We can make a series of diameters through the center of a circle; the angles between the two adjoint semidiameters are equal to each other.Then these small partitions (for the square region, these partitions are triangles) are divided into two regions where the areas are equal to each other.At the end, we can get the minimum absolute value of these two pieces' difference.For example, we can cross the center of the region to make two vertical diameters (such as Figures 3 and 4) and get the valve by the steps mentioned above.From Figures 3 and 4, by the improvement method, we can calculate the RHM's value: just the sum of two triangular regions from head to foot deducts the sum of the other two triangular regions.But according to the algorithm mentioned first, the value of RHM will be zero, so the texture of the picture will be ignored as noise.
Of course, dividing the region with several parallel straight lines will also be a working method (as Figure 5).We can do this step first and then make up them into two regions to get the RHM's value.But each small region for summing cannot be equal to each other, because in each small region the area of the noise must be equal to each other.Compared with the previous method, this division is more difficult to achieve.Hence, we adopt the first kind of method to define RHM.
Therefore, we provide the following definition of the improved RHM.Definition 2 (the improvement RHM).Let (, ) be a circle in which the center is  = ( 0 ,  0 ) and the semidiameter is .Set /2 as integer.The diameters  1 ,  2 ,  3 , . . .,   divide (, ) into 2 part ( 1 ,  2 ,  3 , . . .,  2 ), and the minimum angles which are divided by these diameters are equal to each other, and these 2 partitions make up two equal parts in area without repeated  1 and  2 .The angle between -axis and  1 is , and we write (, , )            , then, we define    as the improved RHM of  and   as the improved RHMO of .

The Adjustment for the RHM.
In the field of digital image processing, the anisotropic diffusion model which is set up by Perona and Malik [1] is a typical one and used widely.
The model replaces the gauss smoothing with coefficient directionality distribution, so it arouses the attention of numerous researchers.A kind of nonlinear parabolic PDE used for image processing has been put forward [3] by on the basis of the forefathers: (, , 0) =  0 (, ) . ( 0 (, ) is the grey gradient of the primitive image.(, , ) is the image when it is at time . is a smooth kernel (e.g., gauss smooth kernel). *  is the local estimation of the gradient ||.() is the nonincreasing real function, while | ⋅ |  → ∞, () tends to 0. This degenerative model treats the image with the alternative smooth method.In this processing, the edge of the image is maintained while the other place of it becomes smooth.
In this paper, we couple the Alvarez model locally with a two-order PDE.Furthermore, we combine the HRMO with this coupling processing; it leads to a new PDE model of anisotropic diffusion.Exactly, we substitute the |∇  * | term with the RHM term (  ).
Through rough estimation, as to a 21 × 21 pixels region (such as a region of a bmp image of 256 grade), the range of    is from 0 to 256 × 21 × 10.This range of the value may cause a relatively large margin in (  ) term in the PDE.But the function of (  ) is used for strengthening the edge of the image.In fact, the (  ) term has controlled the pace spread in the image processing; the value of (  ) is in inverse proportion to rate of diffusion.This may lead to the fact that the image disposed may show a very different effect in each region after the diffusion of the image (in the keeping of the edge).So we consider to adjust the value field of   .
The method we can use is through comparing the convolution function with    term in the value.As to the convolution term |∇  * |, the value field is far little than the range of    (the detailed reason can be seen in the following analysis).But it is difficult to define the convolution term |∇  * |.So we can consider the method used in [2], that is, using the |∇  (; , , , )| to replace |∇  * |.
To the arbitrary (, , ) ∈ ([0, ];  2 (Ω)), we have As to  0 (),   (),  mentioned above, we have the following Fourier expansion: By combining ( 5) with (4), we can get   .Then we can calculate the value of |∇  (; , , , )| 2 .So we can provide the restrain coefficient; that is, In fact, we can provide a roughness estimation about the value of .By calculating a two-dimensional discrete convolution  =  ×  in a 3 × 3 pixels area, where ,  represent the concrete pixel of a picture; , , respectively, represent the convolution kernel and image.As to convolution kernel, we can let for the matrix the convolution of them As to a 3×3 pixels region of bmp image which is 256 grade, the absolute value of the sum of the matrix's element which is given by the convolution kernel is 0 ∼ 256 × 9 (consider the sum between the negative and positive elements of the matrix).So we can get the coefficient of the restraining; that is,  = (256 × 21 × 10)/(256 × 9) ≈ 23.

Setting-Up of the Model
From the above-mentioned analysis, combining with (2), we will propose the improvement model which comes from the Alvarez PDE model (3).It can overcome some shortcomings of the original model and get better result while handling the noise problem.
Let Ω be an open set of R 2 ; let Γ be its smooth border.Consider the following initial-boundary value problem: (, , 0) =  0 .
Similar to Alvarez model,  0 in (12) is the initial signal,  :  + →  + is a decreasing function, and (0) = 1, lim  → +∞ () = 0. Let () = 1/(1 +  2 ); then  → () is smooth.The difference is as follows: in the model,    is the RHM and (, , ) is the unique solution of the following linear questions (achieving the same effect as the convolution function [6]): Through choosing the appropriate threshold value  and solving the equation above to get the solution (, , ), we can get the sharp image and the clear border of the image.

Numerical Scheme for the Model
Suppose that Ω is the Δℎ * Δℎ rectangle area where Δℎ is the unit length.For example, the side of the rectangle neighborhood is  = 20Δℎ.We can write (, , )  .(14) For the situation of double diameter, when the diameter whirl in  ∈ [0 ∘ , 180 ∘ ] one time and ℎ ∈ [0, ], we have furthermore For the situation of many diameters, we can also provide a similar definition by the formula given above.Consider the difficulty of realizing of the process; we can use the method which multiplies with matrix.For example, as to the four diameters situation, we can get the sum of the pixel of every triangle region through the multiplication of the matrix.Then calculate the RHM by the way of making up every triangle region at random (there are  3 7 = 35 kinds of combination altogether).
The pixel's sum of shadow region in Figure 6 (a 5×5 pixels area) is (the sum of other regions can be similarly given) A(i, j) (i − 1, j) (i + 1, j + 1) Figure 8 where For the main equation of (11), we can use the 3 × 3 (as Figure 7) difference scheme.
When we calculate the diffusion operator || div(/||), we must do approximation on the mesh.Suppose that  in which the coordinate is (, ) is the cross point of the mesh.Then we can define the four directions, respectively, as follows: (20 ) to approximate || div(/||), and  is the angle between each direction and ; then   can be defined as follows: where  0 is a constant.Using Δ  to stand for Laplace operator, the two-order differential coefficient about the four directions is  2   / 2  . is step of the time.The difference equation about ( 11) is As to the (13), for obtaining better data precision, we can use the following implicit scheme-the solution of (13) meets the following Douglas-Rachford difference equation: where  = 2, . . .,  − 1;  = 2, . . .,  − 1;  = Δ/4Δℎ 2 and where  ∈ {1, . . ., },  ∈ {1, . . ., }.

The Choice of the Parameters.
In the experiment, we try to make the adjustment to algorithm time () and length of stride (step), diligently making the balance in the computation load and the picture quality.Here, the processed objects which are brain's MRI images are all for 128 × 128 pixels bmp images.The original image and image added to the gauss noise are, respectively, as the pictures below (Figures 8  and 9).The left parts of Figure 10∼Figure 15 that come from Figure 9 are given by the algorithm in this paper.The left part of Figure 10 is the denoising image of Figure 9, while the right side image is edge extraction of the left side image (using the canny operator).The edge extraction image may reflect the quality of denoising image-reflecting the retention of the edge after the denoising.Figure 11∼Figure 15 are also presented in the same way.Following the concrete result, we analyze the following.
(1) The  of Figure 10∼Figure 13 is 1 s, while the step of each other, respectively, is 0.25 s, 0.05 s, 0.02 s, and  (2) The  of Figure 14∼Figure 15 all is 0.005 s, while the step of each, respectively, is 0.0001 s and 0.00005 s; the two images eliminate background noise well; moreover, two images achieve a more ideal effect in the details maintained; comparatively, Figure 14 is clearer than Figure 15 (this can be seen by comparing it with the original image in Figure 8); moreover, the step times of Figure 14 are 50, while of Figure 15 are 100; therefore, considering the denoising effect and the algorithmic load, Figure 14 surpasses obviously Figure 15.
Considering the above two kinds of situations ( and step), in order to get the high image denoising quality, we may use the  and step which Figure 14 uses; but in the situations in which image quality request is not high and computation real-time request is high, the scheme which Figure 11 uses can be considered.

Comparison of the Models.
We will present some numerical results of the different models by the computer simulation.Because we cannot distinguish the noise and true texture in a small region, the value of  cannot be too small.We choose a semidiameter  = 10 during the process of calculating the picture of fingerprint.In this setting, the RHM will not deal with the texture of the fingerprint as the noise.So the most valuable texture will be retained.Figures 23 and 24 are the noising fingerprint and the contours of it.
In order to verify the effectiveness of our model, three different methods which include P-M model, Alvarez model, and improved RHM model are compared with experiments.In Figure 25∼Figure 27,  = 0.004, while the step of each other, respectively, is 0.0001 s.The level set which comes from the diffusion processing of P-M model can be seen in Figure 25.In Figures 26 and 27, we deal in Figure 23 with the Alvarez model and improved RHMO method, respectively.From Figure 25∼Figure 27, we can find out some subtle differences in the contours.In fact, we can find that the P-M model and Alvarez method form some false lines structures in local fingerprint information, while improved RHM method can achieve a good enhanced image, without generating false lines structure.This shows that our model can deal with more complex situations.

Conclusions
In this paper, we do some further researches about RHM based on the analysis of the RHM operator and the property of the image texture.We establish the new model by coupling the Alvarez model locally with a two-order PDE.In the new model, we substitute the |∇  * | term with the RHM term (  ).This model is a new PDE model of anisotropic diffusion.From the convolution operator, we deduce a rational range of the coefficient of RHM in our new PDE model.We test our model by the numerical method.From the experiment of brain's MRI images, we select the different algorithm time and length of stride to achieve a more reasonable diffusion effect of the picture.Comparing the result with some classic models and the original RHM model, we are convinced that our model achieves a good result in the diffusion of subtle texture.

Figure 17 The
Figure 17

Figure 20 see
Figure 20