Iterative Spectral Molliﬁcation and Conjugation for Successive Edge Detection

Fourier concentration method is a well-established method for finding edges given a set of truncated Fourier data. However, it is hard to distinguish weak jumps from the Gibbs oscillations with the method unless the concentration is highly refined locally. In this note, we show that it is possible to find all edges whether strong or weak by using the method iteratively. The iterative method first finds strong edges followed by an iterative procedure in which previously found edges are smoothed by local adaptive mollifier. The iterative method is advantageous when the ratio of the largest to smallest jump magnitudes is significantly large.


Introduction
Detection of edges is a fundamental problem in a variety of applications, including image processing and the numerical solution of differential equations. In applications such as magnetic resonance imaging (MRI), it is required to construct images from Fourier data. Let {f k |k = 0, ±1, ±2, · · · } be the set of Fourier coefficients of f (x) ∈ L 2 [−π, π] given byf k =

2π
π −π f (x)e −ikx dx, and let f N be the Fourier partial sum f N = N k=−Nf k e ikx . When the underlying function is smooth and periodic, the Fourier reconstruction f N is accurate to spectral accuracy, but when edges are present, the reconstruction is plagued by the Gibbs phenomenon, also known as the Gibbs ringing in MRI applications.
Various methods have been proposed to address these issues and those methods consist of edge detection followed by reconstruction. Thus, the determination of edge locations is critical. Fourier concentration method has emerged over the past decade as a robust method for edge detection in a variety of circumstances and applications [5,6]. Essentially, a certain Fourier partial sum converges to the jump function as the number of Fourier coefficients increases and this convergence can be accelerated by what is known as concentration factors (functions). Use of different types of concentration factors tends to impart trade-offs between oscillations near jump discontinuities and significant non-zero concentration away from them [2]. Several methods have been devised to address this issue, as well as to treat special circumstances such as incomplete Fourier data and the presence of noise [1,4,13,15].
Thanks to the convergence property of the Fourier concentration to the jump function, the concentration method detects edges with large concentrations. Where the function is smooth, the concentration vanishes as the jump function vanishes as N → ∞. In practice, the concentration method is designed to detect edges with magnitudes larger than some given threshold, with the value of the used threshold being problem dependent. The value of the threshold cannot be arbitrarily small; otherwise, too many false edges can be detected. If the magnitude of weak edges is much smaller than other edges, those edges are considered insignificant, but for some cases weak edges are more important than strong edges. For example, it was shown that in the segmentation of MRI of the knee, the cartilage is better characterized by weak edges rather than strong edges for the separation from the tibia and femur [11].
This note shows that an iterative approach based on the successive conjugation and adaptive mollification can detect all edges without any prior threshold. This approach is similar to the iterative method in the context of the radial basis function method [3,9,10]. The iterative method is as follows: at each iteration step, all previously found edges are smoothed by a local mollification and new corresponding Fourier coefficients are computed. By applying conjugation and mollification successively, one can distinguish real edges from fake edges. This approach is useful and effective particularly for problems where the weak jump can significantly affect the global solution of differential equations or images where the interesting structure is represented by the weak edges [11].
In Sect. 2, a brief explanation of the Fourier concentration method is given. In Sect. 3, the proposed iterative method is explained based on the adaptive filtering method. The stopping criteria is also explained. Numerical examples with remarks are given. In Sect. 4, a brief concluding remark is provided.

Edge Detection Using Fourier Concentration Method
where the superscripts + and − denote the limits taken from the right and left, respectively. Given a finite set of Fourier coefficients, {f k } |k|≤N , the Fourier concentration method, developed in [5,6], computes the concentration as a sum of the form where the σ (·) are known as concentration factors and sgn(k) is the sign function. Given certain admissibility conditions [6], the sum converges to the jump function: where d(x) denotes the distance to the nearest edge and s depends on the concentration factor. Here we note that Eq. (2) shows that the concentration function S σ N [f ](x) recovers the jump function of f (x) as N → ∞ and the convergence may be slow. Equation (2) also implies that the absolute maximum value of the concentration function S σ N [f ](x) converges to the maximum jump. Accordingly we observe that strong jumps are relatively easier to detect than weak jumps. The common types of concentration factors satisfying the admissibility conditions are polynomial concentrations where p is a positive integer and η = |k|/N and exponential concentration functions where α > 0 is an order and C is a normalization constant. Cutoffs for edge detection, τ ∈ (0, 1], are with respect to the normalized concentration and the edge set, E, is defined as Several approaches have been developed for improving the concentration method. We refer readers some to [1,2,4,6,7,12,13,15,16]. All these methods are basically utilizing the edge map. Figure 1 is by the Fourier concentration method for f 1 (x)

Iterative Concentration Method
As clearly seen in Fig. 1, the Fourier concentration method may fail to detect the weak edge when the concentration of the weak edge is too small compared to the strong edge. To find all edges, we propose to apply the Fourier concentration method iteratively based on the local mollification using the local adaptive filtering method.

Local Adaptive Mollification
The local adaptive mollification is a key step for the iterative algorithm. Consider a smooth function φ ∈ C ∞ 0 [0, 2π ] which is compactly supported such that where With these properties, the limit property is given by where ( * ) operation denotes convolution. The parameter is free and it localizes the convolution and is known as the localization factor. The parameter is a fixed value for every x. Thus a global smoothing occurs everywhere including both the nonsmooth and smooth areas. However, we only want to apply the mollification locally to minimize the Gibbs oscillations near the jump. In order to achieve this, we use a two-parameter family of the spectral mollifier introduced by Gottlieb and Tadmor [8]. Consider the convolution of the Fourier partial sum f N (x) and the mollifier φ. Then by the definition of f N (x) and φ we have where D N is the Dirichlet kernel of degree N . The idea proposed in [8] is that one changes the degree of the Dirichlet kernel with the localization parameter so that the two-parameter family of the new mollifier is defined by where D p is the Dirichlet kernel of degree p. Then for all s, the error is given by [8] |φ where || · || L ∞ loc = sup (x− π,x+ π) | · |. The first term in the right hand side of the above inequality is the truncation error and the second term is the regularization error. As we see, the optimization of the error is determined by how the localization parameter and the degree p are balanced. In [8] those parameters were chosen such that = d π , where d is the distance to the nearest jump from the current position. The order of the Dirichlet kernel is chosen such that spectral convergence is achieved, say, p ≈ √ N . Here we note that a modification of the two-parameter mollifier for the enhancement of the convergence was proposed in [14], which was designed to reduce the Gibbs oscillations while it provides a sharp reconstruction up to the edge. Note that the adaptive mollifier was used to sharpen the concentration mapŜ in [2].
Our proposed iterative method is that once the edge is identified, the edge region is first localized using the value of so that f N in the region away from the detected edge is not affected by the mollification. This helps the next available edge to be preserved through the mollification of f N if existent. Thus as in [8], the localization factor is a function of the distance from the edge, d, i.e.
Then we adaptively mollify f N so that a heavy mollification using p is applied to reduce the Gibbs oscillations near the edge. The limit property of p is given as p → 0 if d → 0 and p → ∞ if d → 2π . In this work, we use the local adaptive filtering for the mollification.

Almost Automatic Stopping of the Iteration
To see the proposed method stops almost automatically, consider f (x) = x, x ∈ [−π, π] with the Fourier coefficientŝ There are two edges (x = ±π ) and conjugation and local adaptive filtering have the most effect at ±π . Therefore, by considering the local behavior near ±π , we assume a constant order of filtering p and of conjugation q, with functional forms of exp(− M (1 − |k/N|) p ) and exp(− M |k/N| q ), respectively. By letting φ p, and C σ N be the corresponding kernels and letting S and F denote conjugation and filtering, and after some simplification, we have This has Fourier coefficients From the sharp localization and heavy filtering near ±π , we choose p, q −→ 0. Then setting y = |k/N| yields which approaches 0 exponentially. Thus after all the edges are found through iteration, the concentration decays exponentially small. Thus if the stopping criteria η below is chosen small enough, e.g. η ∼ 10 −10 , the stopping of the iteration is guaranteed

Numerical Examples
We consider the case that the magnitude of the weak edge is highly small for the function f 1 (x) in Eq. (6) m w = 0.01.  Figure 3 shows similar result as in Fig. 2. As shown in the figure, the method is highly accurate and finds all edges including the highly weak edges.
As an application to the solution of PDEs, namely the shock-density wave interaction equation, we consider finding shocks in the density profile at t = 2 with the total number of grid points N = 300 computed with the WENO-Z method used in [9]. The left two figures of Fig. 4 show the edges (shocks) found by the Fourier concentration method while the right figure shows the edges (shocks) found by the iterative method. As shown in the figure, the iterative method find all the physical shocks accurately while the Fourier concentration method misses some of shocks.
For two-dimensional examples, we consider a Shepp-Logan image with a faint box added to comprise additional weak edges, and a brain image. To detect edges in two dimensions, edges are detected slicewise in the x and y directions. The x and y coordinates have a range of [−π, π]. For a 2N x +1×2N y +1 image, slices of f (x, y) are taken at evenly spaced x and y with x = 2π/(2N x + 1) and y = 2π/(2N y + 1), with −π included and π excluded. Within each slice, Fourier coefficients are computed by partial Fourier expansion and the iterative method is applied to find strong and weak edges. Calculation parameters for the two-dimensional calculations were similar to those for the one-dimensional calculations. An edge with a concentration magnitude at or above a fraction τ = 0.1 of the maximum magnitude concentration was considered strong. To detect strong edges, trigonometric concentration factors with α = π were used. Figure 5 shows the edges found by the proposed method for the Shepp-Logan image. As in the figure, the weak edges (square box with magnitude of 0.01) are successfully found by the method.
Remarks First, the proposed method is affected by noise as the original Fourier concentration. Consider Eq. (6). Letf k (m w = 0) be the Fourier coefficients with m w = 0. Thenf k (m w ) = m w /8 for k = 0 andf k (m w ) = m w 2πk sin(7kπ/4) − sin(3kπ/2) for k = 0. The weak edge translatesf k . Thus we expect that unless SNR is high enough, |f k (m w = 0) −f k (m w )| becomes easily smaller than the noise as m w decays. Figure 6 shows the concentration with m w = 0 (left), the concentration with SNR = 20 (middle) and with SNR = 10. As in the figure, the weak edges are easily indistinguishable as SNR decreases. As the main objective of this research is finding the weak edges, a proper noise reduction suitable for the proposed method should be investigated in our future research. Second, the proposed method is to find  Figure 7 shows the edge detection in the physical domain, with the forward difference, generated from the Fourier data (the left figure). The right figure shows the plot in x ∈ 3π 2 ≤ x ≤ 7π 4 where the weak edges exist. As in the figures, the weak edges are still hard to distinguish in the physical domain. Once the strong edges are removed in the Fourier domain and switching back and forth from the Fourier to physical domains, the weak edges are eventually found with the proposed method.

Summary
The following is the summary of the proposed iterative concentration method. The procedure stops eventually with a non-zero value of η > 0 in Eq. (10).

Conclusion
We showed that the iterative approach of the Fourier concentration method can detect all edges, which is not the case if the weak edges are too small. We showed that the proposed method is able to detect weak edges 3000 times weaker than the strongest edge, as long as the weak edges are well-separated from the stronger edges without noise and that the proposed method find all weak edges in a PDE application, namely the WENO calculation for the shock-density wave interaction. The iterative method also shows that it stops almost automatically after all the edges are found. Thus the proposed method is accurate and efficient.
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence and indicate if changes were made. The images or other third party material in this chapter are included in the chapter's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.