Use of the Wavelet Transform for Digital Terrain Model Edge Detection (Special Issue—Wavelet Analysis)

The purpose of this work is to analyze the feasibility of using the wavelet transform in the edge detection of digital terrain models (DTM) obtained by Laser Scanner. The Haar wavelet transform and the edge detection method called Wavelet Transform Modulus Maxima (WTMM), both implemented in Matlab language, were used. In order to validate and verify the efficiency of WTMM, the edge detection of the same DTM was performed by the Roberts, Sobel-Feldman and Canny methods, chosen due to the wide use in the scientific community in the area of Image Processing and Remote Sensing. The comparison of the results showed superior performance of WTMM in terms of processing time.


Edge detection of images and digital terrain models (DTM) obtained by Remote
Sensing is an important research problem for mapping, especially in large scales (1:10,000 and larger), due to the need for greater detailing of the features and the large amount of digital data of the inputs used for these scales.
The wavelet transform (WT) had its origin in 1909 with the work of Alfred Haar [1] [2] [3] [4], who described the first orthonormal function system and used it for multiscale analysis in his doctoral thesis.However, the formalization of the theory was consolidated only from the 1980s, with the works of Morlet, Mallat, Meyer and Daubechies, among others [1] [5] [6] [7].The wavelet transform has been widely applied in the areas of knowledge that require analysis and visualization at different scales, which makes its proper use in image processing, digital terrain modeling, cartographic generalization and geoprocessing.
The use of WT for edge detection appears in this context as a tool with great potential due to the characteristics of ease of implementation, simplicity of algorithms and speed of processing.In this paper, we used the edge detection method called Wavelet Transform Modulus Maxima (WTMM), described by [8].
Among the traditional edge detection methods, Canny, Roberts and Sobel-Feldman are widely used in the Remote Sensing and Geoprocessing areas and for this reason they were used as reference for the validation, in this work, of the results obtained by edge detection by WTMM.

Fundamentals
In this section are made a brief study of the theoretical foundations and the main works carried out in the context of the edge detection with the operators of Canny, Roberts, Sobel-Feldman and WTMM.Several studies, such as those of [9] [10] [11] [12] [13], compared the edge detection operators of Roberts, Sobel-Feldman and Canny, among others.The general conclusion of these studies is that the Canny operator presented better results than those of Roberts and Sobel-Feldman, but it is a more complex algorithm with a higher computational cost.

Roberts Edge Detector
The Roberts edge detector was created by Lawrence Gilman Roberts in his Ph.D. thesis at MIT in 1963 [14].It is a classic edge detector that uses the directional derivative of the first order.Its implementation consists of a pair of 2 × 2 dimension convolution masks (matrices) as shown in Figure 1 for the detection of gradients.
These matrices provide maximum responses to the edges with a 45˚ direction in relation to the pixel grid and can be applied separately to the input grid to produce separate measurements of the gradient component in each orientation (called G x and G y ).The results can then be agreed to find the absolute magnitude of the gradient at each point and the orientation of that gradient.The magnitude of the gradient is given by: And the direction of the gradient is given by: Figure 1.Roberts operator template.
Roberts edge detector has the advantage of simplicity of the algorithm and disadvantage of noise susceptibility.

Sobel-Feldman Edge Detector
The Sobel edge detector operator, more properly called the Sobel-Feldman operator [15] uses the first derivative, similarly to the Roberts operator, for the detection of gradients and their orientation.The Sobel-Feldman operator consists of a pair of convolution masks (3 × 3) as shown in Figure 2, which produces responses to the 0˚ and 90˚ direction edges in relation to the pixel grid.Like the Roberts operator, the matrices can be applied separately to the input grid to produce separate measurements of the gradient components in each orientation (G x and G y ) and then agreed.By using a larger matrix, the Sobel-Feldman operator has the advantage of reducing errors due to the noise's effects [16].
The magnitude of the gradient is given by: And the direction of the gradient is given by:

3. Canny Edge Detector
The Canny edge detector [17] was created by John Canny in 1986 and is also known as Canny's Great Edge Detector.It performs better than Roberts and Sobel-Feldman, but requires more computational effort.In his work, Canny considered three criteria [18]: 1) Low probability of failure: the ideal detector must find all edges and not find false edges.
2) Location Criteria: The location of the detected edges must match the actual edges.2) Obtaining the gradient with the use of convolution matrices similar to those of Roberts or Sobel-Feldman.In this step, four convolution matrices are generally applied to obtain the edges in 0˚, 45˚, 90˚ and 135˚ directions.
3) Computation of the direction of the gradient.
4) Non-maximum suppression, which consists of thinning edges by suppressing pixels that are not local maxima.
5) Threshold with hysteresis, which consists in the use of thresholds for suppression of unconnected edges.
The Canny operator has more accurate edge detection, but its implementation is more complex and its computational cost is higher.

WTMM: Wavelet Transform Modulus Maxima for Edge Detector
According to Mallat and Hwang [8], if a wavelet function ( ) x ψ has a null moment, it can be taken as the first derivative of a smoothing function ( )

Results
Figures 5-8 are shown the results obtained in the edge and in Table 1 are shown the respective processing times.From the analysis of the figures and the table we can conclude that the Canny detector presented a better result, mainly in the detection of circular edges (northwest sector of the study region), but at a much higher computational cost.On the other hand, WTMM presented a similar result of Sobel-Feldman and better than Roberts, but with a much lower computational cost, although a user-made computer program is usually slower than libraries.The processing time of WTMM was about three times less than Roberts and Sobel-Feldman and 18 times less than Canny.
It was also observed that the WTMM has thicker edges than the other detectors.

Conclusion
The experiments carried out in this research had the objective of verifying the feasibility of using the wavelet transform in the edge detection of the digital terrain model obtained by Laser Scanner.The results show that edge detection by using WTMM can be seen as an alternative to edge detection due to its characte- How to cite this paper: Gaboardi, C. (2018) Use of the Wavelet Transform for Digital Terrain Model Edge Detection (Special Issue-Wavelet Analysis).Journal of Applied Mathematics and Physics, 6, 1997-2005.https://doi.org/10.4236/jamp.2018.610171

3 )
Simple Response Criteria: The detector should provide a single response to a single edge.The Canny detector implementation consists of 5 steps: 1) Application of a filter for smoothing the image, to reduce noise.

Figure 5 .
Figure 5. Edges obtained by the Roberts detector.

Figure 7 .
Figure 7. Edges obtained by the Canny detector.