Abstract

It is difficult to give a fractional global threshold (FGT) that works well on all images as the image contents are totally different. This paper describes an interesting use of fractional calculus in the field of digital image processing. In the proposed method, the fractional global threshold-based edge detector (FGTED) is established using the Riemann–Liouville fractional integral operator. FGTED is used to find the microedges in minimum time for any input digital images. The results demonstrate that the FGTED outperforms conventional techniques for detecting microtype edges. The image with a higher entropy was produced by the FGT value-based approach. Tables and images are used to summarize the output performance analysis of various images using structural similarity index measure, F-score (F-measure), precision and recall, signal-to-noise ratio, peak signal-to-noise ratio, and computational time. The FGTED can be used to detect very thin or microtype edges more accurately in minimum time without training or prior knowledge.

1. Introduction

Fractional calculus is concerned with the mathematical investigation and application of integrals and derivatives of arbitrary order [1]. Given the usefulness and significance of fractional calculus in the analysis of real-world problems, the properties and applications of a variety of fractional integral operators have been studied by a large number of researchers such as McBride [2], Kilbas and Sebastian [3], Kiryakova [4], Baleanu et al. [5], Gaur Sanjay et al. [6, 7], etc. It has been recently found that the application of fractional calculus can be useful for image processing, particularly to enhance the quality of edge detection for the images [810]. Edge detection is one of the fundamental tasks in digital image processing. It plays a prominent role in image analysis, pattern recognition, and other deep-level processing [11, 12].

Edges are a collection of pixels where the intensity of each pixel fluctuates dramatically [13]. To identify, segment and recognize an item in an image, edge detection is a very useful and difficult job in computer vision. It is also used in a system for feature identification, compression, and picture retrieval [1416]. The process of edge detection consists of four main steps: localization, detection, enhancement, and smoothing. The quality of edges diminishes in noisy environments and, in certain cases, after the noise has been suppressed, blurring the image’s major transitions. Edge detection necessitates some sort of smoothing since it relies on differentiating the image function, which amplifies all high-frequency components. Low pass filters are the most commonly employed for this [15].

During the last four decades, the topic of fractional calculus, which deals with the study of arbitrary order integrals and derivatives, has grown in relevance and interest. It is mostly owing to its wide range of possible applications in science, engineering, and technology. Various authors have investigated various expansions of several fractional integration operators and their applications in recent years.

Pu et al. [17] investigated a fractional differential approach to the digital processing field for enhancement. The authors have proposed six fractional multiscale differential masks for the enhancement of images. Traditional integral differential-based algorithms perform worse than the suggested fractional differential-based method. They discovered it to be ideal for real-time HDTV, legibility of subsequent bionic medical images (for example, X images, cell pictures, mammography, PET, MRI, CT, and plutonic images), improved bank ticket quality, and improving remote sensing image, among other applications.

For the numerical solution of fractional differential equations, Diethelm et al. [18] described an Adams-type predictor–corrector approach. They have analyzed the results for both linear and nonlinear problems. Finally, it was recommended that their approach can be extended to multiterm equations using more than one differential operator.

Sarikaya and Yaldiz [19] developed a novel Hermite–Hadamard type fractional integral inequalities utilizing a wide class of fractional integral operators. They also applied fractional integral inequalities to develop inequalities of the Hermite–Hadamard type, employing the well-known Riemann–Liouville fractional integral (RLFI) operators.

Ntouyas et al. [20] studied Polya–Szego type integral inequalities using the RLFI operator and utilized them to show various Chebyshev-type fractional integral inequalities affecting the integral of the product of two functions and the product of two integrals. Khmag et al. [21], Khmag [22], and Khmag and Kamarudin [23] studied clustering-based natural image denoising using a dictionary learning approach in the wavelet domain and semisoft thresholding approach and natural image deblurring using a neural network.

Martin et al. [24] have proposed a Berkeley Segmentation Dataset 500 (BSDS500), which is regarded as the standard benchmark for the detection of the edges. Initially, the given database was designed for investigating natural edge detection, which includes object contours, its interior, and boundaries. This dataset includes five hundred natural-type images with regressive annotated boundaries collected from different users. The abovementioned dataset is used in this paper to compare the performance of the fractional global threshold-based edge detector (FGTED) method with the existing prevailed edge detectors.

The Riemann–Liouville fractional calculus integral operator as left-sided fractional integrals is explained in Section 2.

2. RLFI Operator

The fractional order integral operators are efficient tools to describe the complex phenomenon in much better ways due to their nonlocal nature and having a tendency to preserve hereditary properties related to the image. Since the gray-level values have a high correlation with the neighboring pixels so the fractional order operators are more efficient to describe the hereditary relation and give better neighboring gray-level information. There are several forms of the fractional integrals for a continuous function f(t) have been given by various authors [25]. One of which, known as the Riemann–Liouville fractional calculus operator, has been studied extensively for applications in the area of image processing [811].

Further, the left-sided Riemann–Liouville fractional calculus operator, defined by Kilbas et al. [1] for ς ∈ C, ℜ(ς) > 0, is applied as follows:where Γ (ς) is the Gamma function, the above integrals are called left-sided fractional integrals. When ς = n N, Equation (1) can be solved by the well-known nth integral formula by Cauchy as given below:

Riemann–Liouville fractional integration of the power functions (t − a)ρ−1, (ρ) > 0 is then used to generate power functions of the same form, as given below:

This paper emphasizes upon the need to detect the edge in digital images by introducing a fractional global threshold (FGT) value. Further, the investigation demonstrated an interesting application of the RLFI operator [1] for edge detection in image processing. The edge detection based on noninteger integration of the general statistical average improves the detection selectivity of an edge detector.

3. Edge Detection by Global Thresholding

The proposed methodology for edge detection consists of three steps: preprocessing and normalization of an image, computation of global threshold by RLFI operator, binarisation, and removal of isolated pixels. A block diagram of the working is shown in Figure 1.

After the camera captures the image (Input Image), it is routed via the preprocessing unit, which processes the raw image. Preprocessing focuses on two fundamental operations: grayscale conversion and contrast stretching, which allows the system to perform consistently well for all images.

Furthermore, the edge identification procedure is aided if the RGB color input picture is transformed to a 256 grayscale image using the YIQ color space used by the NTSC color system [5], as shown in the following equation:

Here R(i, j), G(i, j), and B(i, j) are the red, green, and blue intensities of the (i, j)th pixel of the color image, respectively.

In the next step, the filtered image P(i, j) is normalized in the range 0–1 using the following equation:

Further, the left-sided Riemann–Liouville fractional calculus operator for a continuous function f (x) is defined by Kilbas et al. [1] for ς ∈ , ℜ(ς) > 0 is applied as follows:

Riemann–Liouville fractional integration of the power function (x − a)ρ − 1, (ρ) > 0 is then used to generate power functions of the same form as given below:

Now, if f (t) represents the general statistical average, then by virtue of Equation (6), it can be written as follows:

The result of Equation (8) is applied to evaluate FGT for image t(i, j) as follows:

Now, the FGT is used to convert normalized t(i, j) into binary image b(i, j) with the help of computed FGT as follows:

Now Algorithm 1 removes the unwanted pixels, leaving only boundary pixels, and cleans the uncoordinated or broken edges from the binary images (consider n = ∞, to remove the pixels that do not belong to the boundary so that it shrinks to a significantly connected stroke, and other portion with holes shrinks to a ring halfway between the outer boundary and the hole). Finally, the edges of the image are mixed by combining detected edges in the high-frequency as well as low-frequency portion.

Input: Binary image: b of size m × n; Set of neighbor pixels N = 4, ςl= 0.08 and ςh= 0.35.
Output: Optimized edges: oe of size m × n; Set of operations to be executed oe.
1.  if ς = ςl
2.   for all {i and j ∈ N}
3.    if b4 (i, j) = 1
4.     then oel(i, j) = 0;
5.    else oel(i, j) = 1;
6.   end if
7.  end for
8. else if ς = ςh
9.  for all {i and j ∈ N}
10.   if b4 (i, j) = 1
11.    then oeh (i, j) = 0;
12.   else oeh (i, j) = 1;
13.  end if
14.  end for
15. end if
16. then oe = oeloeh

The smaller value of ςl would result in detecting more fine edges in the higher changes in intensity portion of the image. On the other hand, a higher value ςh could detect edges in the lower change of intensity. A tradeoff could be achieved by choosing two different values of ς, which is similar to the feature’s synthesis step in the Canny edge detector. Thus, we have two different values of ς, ςh > ςl (through the experimental results for different values of ς, we found that ςl= 0.08 and ςh= 0.35). Algorithm 1 merges the detected edges for higher and lower values of ς; depending on the satisfaction of connectivity criteria, the edges are linked, which results in the direction of the majority of edges in the image.

4. Results Analysis

The FGTED technique is performed in Matlab-2023(a) with Intel(R) Core(TM) i3-3220 CPU @ 3.30 GHz, and it outperforms traditional approaches. The performance of the FGTED was carried out on the Berkeley Segmentation Dataset and Benchmark (BSDS 500), which is composed of 500 images that have manually annotated ground truth contours for each image. Edge detection accuracy is evaluated using six standard measures: structural similarity index measure (SSIM), F-score (F-measure), precision and recall (PR), signal-to-noise ratio (SNR), peak signal-to-noise ratio (PSNR), and computational time (CT).

The results are shown in Figure 2 for the images, which are randomly selected as experimental samples to compare it with the existing traditional approaches. Figure 2(a)–2(e) illustrates the comparative outcomes produced by several conventional techniques and the proposed FGTED on randomly selected images from BSD500. Figure 2(a) shows the ground truth image, Figure 2(b) shows the edges detected by the Prewitt edge detector, Figure 2(c) shows edges detected by the Sobel edge detector, Figure 2(d) shows the edges detected by Canny edge detector, and finally, Figure 2(e) shows edges detected by the FGTED.

From visual perceptions, Figure 2(e) presents better results than Figure 2(a)–2(d) for both the input images. The FGTED is capable of detecting very thin or microtype edges. The FGT value-based approach yielded a higher entropy image in comparison with the previously existing methods. The results of the statistical performance parameter are summarized and compared in Tables 16.

The visual comparison in Figure 2 also shows that the images depicted by the proposed edge detector are much better for the outer shape of the object than the other images obtained by the different existing edge detectors. The values of different parameters are shown in Tables 16, which also support the inferred results.

The output performance analysis of various images was also done using SSIM, F-score (F-measure), PR, SNR, PSNR, and CT, as shown in Tables 16. The achieved results are as follows:(1)The SSIM, F-measure, SNR, and PSNR are better for FGTED for all the images shown in Tables 1, 2, 4, and 5. The result shows that FGTED performs better compared to all three edge detectors, viz. Prewitt, Sobel, and Canny.(2)Results in Table 3 show the PR score, which is comparatively better than Prewitt and Sobel but lacking for some images compared to Canny, which can be further improved by adaptively changing the threshold as per need of the local region.(3)Table 6 compares the total computation time of all the edge detectors, which also shows that the proposed edge detector is very fast compared to its other rivals.

This shows that the proposed FGTED detects the edges very fast, with significant improvement in the performance parameters compared to the other rival edge detectors.

5. Conclusion

In this article, a novel FGTED is presented for computing microedges in minimum time without training or prior knowledge. It is concluded that Prewitt, Sobel, and Canny edge detectors eliminate very thin or microtype edges, while the proposed FGTED technique takes it into consideration. Figure 2(a)–2(d) shows the edges obtained by different classical edge detection techniques, whereas Figure 2(e) is the edge images detected by the FGTED. The FGTED technique exhibits better performance parameters for texture based images than the other prevailed methods. By the visual perceptions, the above results clearly indicate the superiority of the FGTED technique over the existing techniques. Hence, it is concluded that the FGTED works well for the detection of strong edges in comparison with the classical techniques. However, it is observed that the region near the hat-brim for BSD image no. 189,080 has broken edges, whereas for 135,069, it detected false edges around the birds. The scope of improvement is possible by applying local or adaptive thresholds using a fractional approach.

This technique gives improved results in minimum time, if intensity variations are correlated. Due to the fast speed for computing edges as well as fine outer boundary detection, this approach can be used in real-time applications where fast processing is required in object detection. The authors also suggest that the presented FGTED is a pathway for future innovations in FGT to get better edge detection for more digital images.

Data Availability

No data were used to support this study.

Conflicts of Interest

There are no conflicts of interest regarding the publication of this article.