1. Introduction
In recent years, fractional differential theory has gradually attracted great attention from engineering and technical personnel, and it has been applied to electrochemistry [
1], material mechanics [
2], fluid field theory [
3], control theory [
4,
5], electromagnetic field theory [
6], biomedicine [
7], and signal processing [
8]. Fractional differential theory has achieved fruitful research results. It was found that fractional calculus theory is suitable for studying uncertain signals with nonlinear, non-causal, and non-stationary characteristics. Therefore, many classical algorithms in fractional calculus theory are very suitable for modern signal analysis and processing. In this paper, the fractional differential model was applied to the image enhancement of the on-line structured light vision measurement of strongly reflective surfaces.
In the process of light stripe image enhancement, the key is to enhance the speckle noise highlighted by the sense of particles while preserving the effective fringe image area with continuous gray scale. The commonly used image enhancement methods, both in the spatial domain and the frequency domain, have the contradiction of enhancing the high-frequency image and weakening the low-frequency part of the image [
9]. For example, the histogram equalization method enhances the high-frequency component of the image, and the image histogram obeys the uniform distribution; however, the edge of the image is excessively cut and sharpened, which affects the identification of image information [
10]. Later, the theory of Retinex proposed by Land [
11] and others was widely used in image enhancement. The most classical algorithms are single-scale Retinex (SSR) [
12] and multi-scale Retinex (MSR) [
13]. They all use the Gaussian filter as a low-pass filter. Although the algorithm is simple, the image halo phenomenon is serious after enhancement, and the algorithm is based on the time-domain enhancement method. Additionally, the details are lost when processing the image. Andrivanov, Nikita et al., based on the deep double stochastic Gaussian model combined with nonlinear double stochastic filtering, achieved that the error of the recovered image was only 9% when only 50% of the image information was used [
14]. Integer order differential gradient operators (first-order differential operators such as Sobel, Prewitt, and Canny operators; second-order differential operators such as Laplace operators) are usually effective in enhancing high-frequency details, but at the same time, low-frequency components are lost [
15,
16].
2. Fractional Differential Image Enhancement Model
Differential operation can enhance the high-frequency component of the signal. As shown in
Figure 1, the second-order differential (
) increases the amplitude of the high-frequency part (
), but the amplitude of the low-frequency part decreases significantly. The main difference between speckle noise on a strongly reflective surface and an effective continuous light stripe lies in the gray dispersion of image, that is, the difference between granular and continuous stripes. If the differential enhancement image processing method is used to highlight the graining characteristics of noise, then the high-frequency variation characteristics of the light stripe edge can be highlighted, but at the same time, it is hoped to effectively protect the information of low-frequency continuous fringes to extract the effective light stripe center by using the conventional gray center of the gravity method after removing the speckle noise of a strongly reflective surface. If the differential order is
, although the rise degree of the curve of the high-frequency part (
) is lower than that of the integer order differential, then the amplitude of the low-frequency part is significantly higher than that of the integer order. Using fractional differentiation to process the image with a strong reflection speckle can highlight the edge of the image and enhance the texture information at the same time.
Under the Euclid measurement, the Grumwald-Letnikov fractional calculus is defined as [
17,
18,
19]:
∀∈R. Rounding
is
. If signal
has a continuous
derivative, when
, m takes at least the integer as
, then the
-order fractional derivative of
can be found in the following equation:
In the formula, stands for step size, represents the signal value, is the degree of continuous derivative, represents the time, and and represent integer order and fractional order, respectively.
To make
the limit of approaches zero,
when
and
. Namely,
. For image signals
, the fractional partial differential is obtained for two independent variables
and
, respectively, and the corresponding backward difference expression is obtained by the following [
20,
21,
22,
23],
Similarly, fractional differential expressions in six other directions such as negative direction, negative direction, and left-right diagonal direction can be obtained.
In Equations (3) and (4), the corresponding coefficients of the first
term are equal, and the coefficients of the first
term are as follows:
Generally, for the image signal
of
, M, N represents the height and width of the image. the image processing in the form of
is expressed as:
where
,
and
are strong reflection image filter mask operators. In order to obtain a complete enhanced image, all pixels whose distance between the center point of the mask and the edge of the original image
is within
and
need to be fractional differential operation in eight directions one by one, as shown in
Figure 2.
The same-sex filter in each direction constructed by Equation (5) can obtain the
fractional differential operator in eight symmetrical directions, as shown in
Figure 3, and the elements not in the eight symmetrical operation directions are assigned a value of 0.
The fractional differential mask operator shown in
Figure 3 has rotational isotropy and symmetry. In the application, by comprehensively analyzing the complexity of the differential operator and the effect of image enhancement, the operation scale is selected. In this paper, a
operator as shown in
Figure 4 was used as the strong reflection image filtering operator
.
3. Fractional Differential Image Enhancement Model and Performance Evaluation
In order to make the best use of the autocorrelation of adjacent pixels and to make the weight of all pixels with a distance of 2 half of the original non-zero weight, the pixels with a weight of 0 were not omitted. At the same time, in order to make the pixel value of the enhanced image change within a certain range of the original pixel value, each item of the mask operator in
Figure 4 was divided by
for normalization. W
0, W
1, and W
2 were used to represent the values of 8a
0, a
1, and a
2 divided by
in
Figure 4, respectively. The result is shown in Formula (7). The value in the
mask was replaced with W. The final
mask operator is shown in
Figure 5.
The light bar images of three types of diffuse reflection samples were selected: PCB board, metal roller, and hard disk substrate. The image resolution was 1000 dpi and the size was 1200 pixels. The Sobel algorithm, Laplace algorithm, and fractional differential algorithm (order
v = 0.5,
v = 0.7 and
v = 0.8) were selected for enhancement. The enhancement effects of different algorithms are shown in
Figure 6,
Figure 7 and
Figure 8.
Table 1 shows the entropy and gradient values of the three images processed by different enhancement algorithms. The entropy and gradient values of the image after fractional differential enhancement were significantly higher than those of the other two algorithms, indicating that the fractional differential method retained more texture details after image enhancement. However, with the change of fractional order, it was found that when v = 0.7, the texture was significantly richer than that when v = 0.5, while when v = 0.8, there was excessive sharpening and enhancement. Therefore, v = 0.7 was selected in this paper.
4. Experiment and Results of Structured Light Vision Measurement
After the light stripe image was enhanced by the fractional differential of the structured light vision measurement, the continuous light stripe and speckle noise were separated by connecting area statistics, and the speckle interference was removed to obtain the light stripe image. The effective fringe center was extracted using the gray center of gravity. The algorithm flow is shown in
Figure 9. The center coordinates of the light bar image were obtained by the gray center of gravity method:
where
G(
) is the gray value of row i and column j in the light bar image, and (
,
) is the gray value of the extracted light bar image coordinates of the heart.
In order to verify the superiority of the method in this paper, a comparative experiment was set up to compare the algorithm in this paper with Laplacian denoising.
Figure 10 shows the processing of the hard disk baseband highlight area by different algorithms. After denoising using the Laplace method, the center of light stripe was obtained by the gray center of gravity method. As shown in
Figure 10a, the red ellipse marks the strong reflection area, and the extraction of the center of light stripe failed. The center of the light stripe was extracted by the fractional differentiation algorithm, as shown in
Figure 10b. The image under the dark background dark light bar was taken and the light bar center was extracted, as shown in
Figure 10c. The light bar center in the strong reflection area was extracted successfully. The performance comparison of extracting the light stripe center after the two methods is shown in
Figure 10d. It can be seen from the figure that the error of the average extraction result of the fractional differential algorithm in the center of the hard disk baseband was reduced by 0.53 pixels compared to the Laplace algorithm. This proves that the algorithm in this paper had better performance when dealing with bright surfaces.
Images of the PCB board and roller were used for comparative experimental verification.
Figure 11 shows the light stripe image in the strong reflection area of the PCB board. The center of the light stripe was extracted after the Steger method processing, as shown in
Figure 11a; the center of the light stripe was extracted after fractional differential processing, as shown in
Figure 11b; the center of the light stripe was extracted from the image under the dark background, as shown in
Figure 11c; and the performance comparison is shown in
Figure 11d. The average error of the extraction center of fractional differential processing was reduced by 0.61 pixels.
Figure 12 shows the light stripe image in the strong reflection area of the roller. The center of the light stripe was extracted after the Steger method processing, as shown in
Figure 12a; the center of the light stripe was extracted after fractional differential processing; as shown in
Figure 12b; the center of the light stripe was extracted from the image under dark background, as shown in
Figure 12c; and the performance comparison is shown in
Figure 12d. The average error of the extraction center of fractional differential processing was reduced by 0.9 pixels.
After the strong reflection area was enhanced by fractional differentiation, the image outline was clearer, as shown in
Figure 13.
5. Comparative Experiment
In order to verify the effectiveness of the algorithm in this paper, five algorithms in recent years were selected, and a comparative experiment was performed by processing a PCB and metal roller. Ref. [
24] firstly Fourier transformed the pixels in the spatial domain and then performed inverse Fourier transform to obtain a synthesized image. Ref. [
25] selected the maximum value in each pixel channel to initialize the image light map, then refined the initial light map by strongly adding a structural prior, and finally synthesized the enhanced map according to the Retinex theory. Ref. [
26] showed the pixel difference through a two-dimensional histogram, and used the image enhancement mapping function to process the histogram to export the image. Ref. [
27] performed filtering based on the Laplacian operator, added a patch difference model, used the generalized Gamma distribution for parameter estimation, and determined the image enhancement algorithm according to the self-guided and external-guided characteristics of the filter. Ref. [
28] adopted an improved AGC algorithm to achieve the CE of bright images with a new negative image strategy and used truncated CDF-modulated gamma correction to enhance dark images, which can alleviate local over-enhancement and structural distortion. The comparison results of PCB and metal roller processing are shown in
Figure 14 and
Figure 15.
6. Image Quality Assessment
Using NIQE (Natural Image Quality Evaluator), the design idea of this model was based on constructing a series of features for measuring image quality, and these features were used to fit a multivariate Gaussian model. These features were extracted from some simple and highly regular natural landscapes. This model actually measured the difference in the multivariate distribution of a test image, which was constructed with these features extracted from a series of normal natural images.
The NIQE quality evaluation model does not require the main performance evaluation score of the original image. It extracts image features from the original image library and then uses a multivariate Gaussian (MYG) model for modeling. Generally, the distance between the NSS feature model and the MVG extracted from the distorted image features is used to represent the quality of the distorted image, as shown in the Formula (9):
represent the mean vector and covariance matrix of the natural MVG model and the distorted image MVG model, respectively. A smaller result for this value indicates better image quality. The picture NIQE indices of the PCB and metal roller processed by different algorithms are shown in
Table 2 and
Table 3.
It can be seen from the above table that, by comparing the NIQE values of the PCB and metal roller images processed by several algorithms, the fractional differential algorithm proposed in this paper has the minimum NIQE value in the case of dark light source; that is, it is the best match with the original image. It can be seen that the fractional differential algorithm in this paper has a good effect in processing high reflection images.
7. Conclusions
In this paper, a fractional differential image enhancement model was applied to structured light vision measurement to enhance texture information while highlighting image edges. This paper selected three types of diffuse reflection samples of PCB, metal roller, and hard disk substrate and selected the Sobel algorithm, the Laplace algorithm, and the fractional differentiation algorithm (the order is v = 0.5, v = 0.7, and v = 0.8, respectively) for enhancement. The image entropy and gradient after order difference enhancement were obviously higher than those of the other two algorithms. When fractional order v = 0.7, the effect was more detailed. The center of the light band was obtained by the Gray barycentric method, and the center of light band was extracted by fractional differential algorithm. The experimental results show that the average extraction error of the fractional differential processing center was reduced by 0.53 pixels. Finally, for the strip image in the strong reflection area of the PCB, the Steger method was used to extract the strip center, and fractional differential processing was used to extract the strip center. By comparing the image results of the algorithm in this paper with several existing algorithms, as well as the quality evaluation index NIQE, it was found that the average error of the extraction center of fractional differential processing was relatively small, and the image outline was clearer. Because the algorithm can highlight the surface of the object more clearly, the algorithm can be used in factory production for the quality of workpiece classification, and it can be used in industrial production. Compared to the deep learning algorithm, this algorithm has some defects. It cannot be detected in real time and is greatly affected by environmental factors. In the future, this algorithm is considered to be extended to the deep learning framework, which is expected to improve the detection accuracy and efficiency.