Abstract
Easy to use, oriented half kernels are reliable in image analysis. These thin filters, rotated in all the desired directions are useful to detect edges, or extract precisely their orientations, even concerning highly noisy images. Usually, the filtering process corresponds to convolutions with Gaussians and their derivatives. Other filters exist and can be implemented in order to build half kernels. However, functions used for the smoothing and derivative parts have not been studied in depth. The goal of this paper is to evaluate different types of half filters as a function of the noise level. The studied kernels have the same spatial support, enabling easier comparisons. To address the robustness of the studied filters against noise, the image quality is gradually worsened. Then, their performances are compared through objective evaluations of both segmentation and gradient direction.
You have full access to this open access chapter, Download conference paper PDF
Similar content being viewed by others
Keywords
1 Introduction
Edge detection is a fundamental process and remains widely used in image analysis and computer vision applications. Moreover, gradient and edge orientations are essential information for the interpretation and exploitation of digital images [19, 26]. Classical methods use first-order fixed operators [4, 5, 21, 24] or the structure tensor [11] to detect gradient and its orientation. Gradient-based edge and orientation estimators are frequently utilized. Among these operators, DoB (Difference of Boxes) filters [21] often remain implemented in many applications as in robotics. DoB simply computes the difference between two mean filters of same sizes on both sides of the considered pixel (usually represented by a list of 1 or \(-1\)). However, DoB filters are sensitive to noise and orientations not reliable. The robustness of the detection is generally improved by smoothing the image by a low pass filter followed by a gradient detection using vertical and horizontal masks of type [\(-1\) 0 1]. Edge detection methods differ in the types of smoothing filters that are applied [4,5,6, 24], as equations presented in Table 3. The computations of the gradient magnitude \(|\nabla I|\) and its orientation \(\eta \) are presented in Table 1(left). Among all the edge detectors, box [21] and exponential [4, 24] filters do not delocalize contour points [12], whereas they are sensitive to noise and blur while contours near a corner are rounded by Gaussian and Deriche filters (cf. Fig. 2(a)). In the orientation and edge detection domains, steerable feature detectors represent popular and efficient tools [8,9,10]. This concept provides valuable orientations and edges for many image processing applications, especially when the filter is elongated (cf. Fig. 2(b)). The concept was generalized in [20] by decomposing a given filter kernel optimally in a set of basic filters approximating an Anisotropic Gaussian Kernel (AGK). Mathematically, Table 1(middle column) presents the orientation and the gradient magnitude computation using this type of filter, corresponding to the maximum energy in the maximal response direction. The AGK possesses a common shortcoming, as a matter of fact, only one \(\pi \)-periodic orientation is extracted efficiently [20], so the impossibility of these filters to estimate in a relevant way several coexisting orientations at the same pixel, as illustrated in Fig. 1(d). Otherwise, multi steerable filters are designed and adjusted to a pattern of interest or a region to determine more precisely the directions of edges with \(2\pi \)-periodic templates [13, 18, 23]. On the other hand, the asymmetric filter developed in [16] estimates the homogeneity in multiple local directions. This filter consists in a network of several parallel lines in which a homogeneity is computed and enables an estimation of edge directions modulo 2\(\pi \). Lastly, wedge [18, 23] or asymmetric oriented filters [15, 16] sound better suited to multiple edge directions detection or modeling a template. The next section details the edge detection process using \(2\pi \)-periodic filters.
2 Edge Detection Using Half Oriented Kernels
2.1 Advantages of Half Oriented Kernels
Edge detection techniques using elongated kernels are efficient to detect large linear structures correctly [9, 10, 20, 22]. The robustness against noise depends strongly on the smoothing parameters of the filter, i.e., the parameter of the filter elongation. If the filter length increases, the detection is less sensitive to noise, but small structures are considered as noise and their edges are not extracted. Consequently, the accuracy of the detected edge points decreases strongly at corner points and for non-straight object contour parts, as illustrated in Fig. 1(d). To bypass this undesirable effect, an anisotropic edge detection method is developed in [15, 17]. Indeed, the proposed technique is able to detect crossing edges and corners due to two elongated and oriented filters in two different directions. The simplest solution is to consider paths crossing each pixel in several directions. The idea proposed in [15, 17] is to “cut” the elongated kernel using a Heaviside function into two parts, i.e. two directions (see Fig. 2(c)). Then, the half kernel (HK) is rotated in several directions from 0 to \(2\pi \) (bilinear rotation) and computes a derivative information at each desired angle (illustrated in Fig. 3, for each \(\pi /36\) radian angle). Thus, in a contour direction, a half smoothing is performed, whereas in the perpendicular direction, a derivative information is computed, as illustrated in Fig. 3(a),(d)–(g). In order to better understand this technique to extract edge, the filter support of a HK on a straight contour is equivalent to 1/2 on both sides of the edge, as for elongated and circular or box filters, illustrated in Fig. 2. On the contrary, for a corner point with a \(\pi /2\) radian angle, the support of the half filter remains 1/2 on both sides of the edge, whereas it is around 1/4 and 3/4 concerning other filter supports. Such a manner to extract edges enables to compute a gradient of blurred contours, even corners as illustrated in Fig. 1(e), contrary to other filtering methods. Eventually, as detailed in Table 1(right), the gradient corresponds to the maximum value minus the minimum value of the convolution of the oriented half filter with the image. These two directions represent to the two main orientations of a contour and the gradient direction \(\eta \) corresponds to the bisector between these two directions (cf. Fig. 2(c)). These orientations are useful and efficient for image restoration via PDE [14], corner detection [3] or image descriptor [25].
2.2 Types of Implemented Filters
Edge detection using oriented half kernels is reliable, robust against noise, and, outperforms classical edge detectors [14, 15]. As pointed in [17] and illustrated in Figs. 5(l)–(q), the HK can be implemented using different equations:
-
a derivative part of the filter, represented by an equation in Table 2,
-
a half smoothing part, which corresponds to an equation in Table 3.
In order to create an elongated filter, the support of the half smoothing part must be higher than the derivative support. The parameters of the different filter equations can be chosen in Table 5 for the half smoothing part, where S represents the length of the support where 98% of the filtered information is computed, i.e., the width of the HK (cf. Fig. 3(c)). As for the derivative part, the parameters are available in Table 4, where D represents the length of the support where 98% of the filtered information is computed (i.e., width of the HK, cf. Fig. 3(b)). As an example, a half filter using Bourennane equation, of length \(S = 46\) pixels and \(D = 11\) pixels, the parameter corresponds to \(s = 3.22\) for the smoothing part and \(s = 0.33\) for the derivative part. The half filter can be generated using a half DoB, see first rows in Tables 2 and 3. Thus, S and D numbers, indicate directly the size of the half DoB, which is rotated from 0 to \(2\pi \). It is important to note that the parameters of the half smoothing parts in Table 5 correspond to parameters representing a pass bass filter centered at 0, and only the causal part is preserved here, so S corresponds to a smoothing support size of a half filter. As a result, the derivative part computes the oriented gradient information, whereas the half smoothing part enables a smooth along edges, preserving them and their directions. Finally, on the one hand, the length of the HK must be sufficiently large to be robust against noise. On the other hand, the width of the HK must be thin to ensure precise edge detection and orientations. As an example in Fig. 3(left), when S and D are small, all the signals are noisy, whereas when S increases, the edge directions highlight incrementally, even though Shen and Bourennane HK remain always sensitive to noise.
The Matlab code is available on the MathWorks website: https://fr.mathworks.com/matlabcentral/fileexchange/66853-edge-detection-methods-based-on-oriented-half-kernels?s_tid=srchtitle.
3 Evaluation and Results
The aim of the experiments is to obtain the best edge map in a supervised way. In the experiments, 7 edge detection methods based on filtering gradient computation are compared: isotropic Gaussian filter [5], AGK [9] and HK formed by the DoB, Shen [24], Bourennane [4], Deriche [6], and Gaussian filters [15, 17]. These HK possess the same spatial support (AGK has two times the spacial support, Fig. 5(m)). Filters must be thin and elongated in order to compute accurate contours and associated directions (cf. Fig. 3(d)–(g)). Thus, the derivative parameter concerning isotropic Gaussian filter [5] remains the same than others (\(D=7\)), see Fig. 5(l)–(q). HK and AGK are rotated each \(\pi \)/72 radian angles. Finally, after a non-maximum suppression [26], an objective assessment is performed by varying hysteresis thresholds on normalized thin edges until the Relative Distance Error (RDE, inspired by [7]) evaluation obtains the minimum score [2]:
where \(D_c\) represents the desired contour map, \(G_t\) the ground truth edge image (Fig. 5(b)), \(|\cdot |\) the number of contour points of a set, and, \(d_{A} ( p)\) is the minimal Euclidian distance between a pixel p and a set A [1, 7]. These scores are recorded and plotted as a function of the noise level in the original image, as presented in Fig. 5(j). Hence, a plotted curve increases monotonously with the noise level (Gaussian noise). The second evaluation measure concerns the gradient angle, \(\eta \) (cf. Table 1). Once \(D_c\) is created, considering \(\mathcal{C}_{D_c}\) the set of contour chains in \(D_c\) (i.e., at least 2 pixels per chain), the gradient evaluation as follows:
where \(d_k\) represents a contour pixel belonging to W, a \(3\times 3\) window centered on p, \(\overrightarrow{\eta _{d_k}}\) the gradient orientation of \(d_k\) and \(c_k\) the number of contour pixels in W, minus the central pixel. This evaluation linearly ranges from 0 for identical angles of \(\overrightarrow{\eta _{p}}\) and \(\overrightarrow{\eta _{d_k}}\) to 1 for angles that differs (Fig. 4). Note that angles of \(\overrightarrow{\eta _{d_k}}\) and \(\overrightarrow{\eta _{p}}\) belong to \([0; \, \pi [\) and when one direction approximates 0 and the other direction \(\pi \), the evaluation of these two directions remains close to 0 (see Fig. 4(c)).
Illustration of \({\eta _{p}}\) and \({\eta _{d_k}}\) directions, in a \(3\times 3\) window, where each white square represents a contour pixel of \(D_c\), whereas gray squares correspond to non-edge pixels. Example of angular single error computation \( e = 1 - {abs(\pi / 2 - abs( {\eta _{p}} - {\eta _{d_k}}))}/{\pi /2}\).
The presented segmentations in Fig. 5 correspond to the original image for a PSNR = 14 dB. Clearly, all the HK obtains better segmentations than the isotropic Gaussian. The HK using Gaussian outperforms six other compared filters, in terms of both RDE evaluation and visualization. AGK obtains good scores, but the contour image presents undesirable straight contours disturbing its interpretation. HK using DoB produces many false positive points, especially close to the true edges; it performs in presence of salt&pepper noise, but not concerning white Gaussian noise. Otherwise, HK using Shen and Bourennane filters obtain close the same results, however, they remain too sensitive to noise, and, especially, many edges are missing. HK with Deriche filter is more robust to noise but edges are still missing. Concerning \(\eta \) evaluation, \(E(\mathcal{C}_{D_c}, \!\eta )\) indicates that HK using Gaussian remains more reliable than other HK. This information means that HK using Gaussian is less disturbed by noise than other kernels. The gradient evaluation of AGK is close to 0 because almost all the edges obtained are straight; consequently, the gradient angle remains close the same for all the contour chain.
4 Conclusion
Oriented half kernels (HK) represent thin filters enable to describe the image information all around a considered pixel. As they are robust against noise and their edge directions are accurate, HK are utilized in the context of many image processing problems [3, 14, 25]. This study presents different manners to build HK devoted to edge detection in digital images. Five HK have been objectively compared and HK using Gaussian remains the most reliable. Especifically, this HK and the anisotropic Gaussian kernel obtain similar evaluation results kernel but the obtained edges are more usable with HK. For this purpose, we plan in a future study to implement multiscale and subpixel half Gaussian kernels.
References
Abdulrahman, H., Magnier, B., Montesinos, P.: From contours to ground truth: how to evaluate edge detectors by filtering. J. WSCG 25(2), 133–142 (2017)
Abdulrahman, H., Magnier, B., Montesinos, P.: A new objective supervised edge detection assessment using hysteresis thresholds. In: International Workshop on Brain-Inspired Computer Vision, held as part of ICIAP, pp. 3–14 (2017)
Abdulrahman, H., Magnier, B., Montesinos, P.: Oriented asymmetric kernels for corner detection. In: IEEE EUSIPCO, pp. 778–782 (2017)
Bourennane, E., Gouton, P., Paindavoine, M., Truchetet, F.: Generalization of Canny-Deriche filter for detection of noisy exponential edge. Sig. Process. 82(10), 1317–1328 (2002)
Canny, J.F.: A computational approach to edge detection. IEEE TPAMI 8(6), 679–698 (1986)
Deriche, R.: Using Canny’s criteria to derive a recursively implemented optimal edge detector. Int. J. Comput. Vis. 1, 167–187 (1987)
Dubuisson, M.P., Jain, A.K.: A modified Hausdorff distance for object matching, In: 12th IAPR International Conference on Pattern Recognition, vol. 1, pp. 566–568 (1994)
Freeman, W.T., Adelson, E.H.: The design and use of steerable filters. IEEE TPAMI 13, 89–906 (1991)
Geusebroek, J.-M., Smeulders, A.W.M., van de Weijer, J.: Fast anisotropic gauss filtering. In: Heyden, A., Sparr, G., Nielsen, M., Johansen, P. (eds.) ECCV 2002. LNCS, vol. 2350, pp. 99–112. Springer, Heidelberg (2002). https://doi.org/10.1007/3-540-47969-4_7
Jacob, M., Unser, M.: Design of steerable filters for feature detection using Canny-like criteria. IEEE TPAMI 26(8), 1007–1019 (2004)
Köthe, U.: Reliable low-level image analysis. Habilitation thesis (2007)
Laligant, O., Truchetet, F., Meriaudeau, F.: Regularization preserving localization of close edges. IEEE Sig. Process. Letters 14(3), 185–188 (2007)
Magnier, B., Montesinos, P., Diep, D.: Texture removal by pixel classification using a rotating filter. In: IEEE ICASSP, pp. 1097–1100 (2011)
Magnier, B., Montesinos, P.: Evolution of image regularization with PDEs toward a new anisotropic smoothing based on half kernels. In: IS&T/SPIE Electronic Imaging, International Society for Optics and Photonics, p. 86550M (2013)
Magnier, B., Montesinos, P., Diep, D.: Fast anisotropic edge detection using gamma correction in color images. In: IEEE ISPA, pp. 212–217 (2011)
Michelet, F., Da Costa, J.-P., Lavialle, O., Berthoumieu, Y., Baylou, P., Germain, C.: Estimating local multiple orientations. Sig. Process. 87(7), 1655–1669 (2007)
Montesinos, P., Magnier, B.: A new perceptual edge detector in color images. In: Blanc-Talon, J., Bone, D., Philips, W., Popescu, D., Scheunders, P. (eds.) ACIVS 2010. LNCS, vol. 6474, pp. 209–220. Springer, Heidelberg (2010). https://doi.org/10.1007/978-3-642-17688-3_21
Mühlich, M., Friedrich, D., Aach, T.: Design and implementation of multisteerable matched filters. IEEE TPAMI 34(2), 279–291 (2012)
Papari, G., Petkov, N.: Edge and line oriented contour detection: state of the art. Image Vis. Comput. 29(2), 79–103 (2011)
Perona, P.: Steerable-scalable kernels for edge detection and junction analysis. IVC 10(10), 663–672 (1992)
Prewitt, J.M.S.: Object enhancement and extraction. In: Lipkin, B., Rosenfeld, A. (eds.) Picture Processing and Psychopictorics. Academic Press, New York (1970)
Püspöki, Z., Storath, M., Sage, D., Unser, M.: Transforms and operators for directional bioimage analysis: a survey. In: De Vos, W.H., Munck, S., Timmermans, J.-P. (eds.) Focus on Bio-Image Informatics. AAECB, vol. 219, pp. 69–93. Springer, Cham (2016). https://doi.org/10.1007/978-3-319-28549-8_3
Simoncelli, E.P., Farid, H.: Steerable wedge filters for local orientation analysis. IEEE TIP 5(9), 1377–1382 (1996). Citeseer
Shen, J., Castan, S.: An optimal linear operator for edge detection. IEEE CVPR 86, 109–114 (1986)
Venkatrayappa, D., Montesinos, P., Diep, D., Magnier, B.: A novel image descriptor based on anisotropic filtering. In: Azzopardi, G., Petkov, N. (eds.) CAIP 2015. LNCS, vol. 9256, pp. 161–173. Springer, Cham (2015). https://doi.org/10.1007/978-3-319-23192-1_14
Ziou, D., Tabbone, S.: Edge detection techniques: an overview. Int. J. Pattern Recogn. Image Anal. 8(4), 537–559 (1998)
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2018 Springer International Publishing AG, part of Springer Nature
About this paper
Cite this paper
Magnier, B. (2018). An Objective Evaluation of Edge Detection Methods Based on Oriented Half Kernels. In: Mansouri, A., El Moataz, A., Nouboud, F., Mammass, D. (eds) Image and Signal Processing. ICISP 2018. Lecture Notes in Computer Science(), vol 10884. Springer, Cham. https://doi.org/10.1007/978-3-319-94211-7_10
Download citation
DOI: https://doi.org/10.1007/978-3-319-94211-7_10
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-94210-0
Online ISBN: 978-3-319-94211-7
eBook Packages: Computer ScienceComputer Science (R0)