Next Article in Journal
Evaluation of Daytime Evaporative Fraction from MODIS TOA Radiances Using FLUXNET Observations
Next Article in Special Issue
Assessing Long-Term Changes in the Beach Width of Reef Islands Based on Temporally Fragmented Remote Sensing Data
Previous Article in Journal
Modifying SEBAL Model Based on the Trapezoidal Relationship between Land Surface Temperature and Vegetation Index for Actual Evapotranspiration Estimation
Previous Article in Special Issue
Comparing Two Photo-Reconstruction Methods to Produce High Density Point Clouds and DEMs in the Corral del Veleta Rock Glacier (Sierra Nevada, Spain)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 6
Issue 7
10.3390/rs6075938
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TecLines: A MATLAB-Based Toolbox for Tectonic Lineament Analysis from Satellite Images and DEMs, Part 1: Line Segment Detection and Extraction

by
Mehdi Rahnama
1,2,* and
Richard Gloaguen
1,2
1
Remote Sensing Group, Institute of Geology, TU Bergakademie Freiberg (TUBAF), Bernhard-von-Cotta-Str. 2, 09599 Freiberg (Sachsen), Germany
2
Remote Sensing Group, Helmholtz Institute Freiberg of Resource Technology, Halsbruecker Str. 34, 09599 Freiberg (Sachsen), Germany
*
Author to whom correspondence should be addressed.
Remote Sens. 2014, 6(7), 5938-5958; https://doi.org/10.3390/rs6075938
Submission received: 22 April 2014 / Revised: 11 June 2014 / Accepted: 11 June 2014 / Published: 25 June 2014
(This article belongs to the Special Issue Remote Sensing in Geomorphology)
Browse Figures
Versions Notes

Abstract

:
Geological structures, such as faults and fractures, appear as image discontinuities or lineaments in remote sensing data. Geologic lineament mapping is a very important issue in geo-engineering, especially for construction site selection, seismic, and risk assessment, mineral exploration and hydrogeological research. Classical methods of lineaments extraction are based on semi-automated (or visual) interpretation of optical data and digital elevation models. We developed a freely available Matlab based toolbox TecLines (Tectonic Lineament Analysis) for locating and quantifying lineament patterns using satellite data and digital elevation models. TecLines consists of a set of functions including frequency filtering, spatial filtering, tensor voting, Hough transformation, and polynomial fitting. Due to differences in the mathematical background of the edge detection and edge linking procedure as well as the breadth of the methods, we introduce the approach in two-parts. In this first study, we present the steps that lead to edge detection. We introduce the data pre-processing using selected filters in spatial and frequency domains. We then describe the application of the tensor-voting framework to improve position and length accuracies of the detected lineaments. We demonstrate the robustness of the approach in a complex area in the northeast of Afghanistan using a panchromatic QUICKBIRD-2 image with 1-meter resolution. Finally, we compare the results of TecLines with manual lineament extraction, and other lineament extraction algorithms, as well as a published fault map of the study area.

1. Introduction

Detection and extraction of lineaments are commonly used for construction site selection (dams, bridges, roads, etc.) [1–4], seismic and risk assessment [5–11], water resources and hydrogeological investigations [12–16], mineral exploration [17–22], and in the study of the structural or tectonic history of a region [23–25]. Most of the tectonic features are associated with straight linear elements (further denoted as “lineaments”) on satellite images. Several regions cannot be studied in the field for various reasons such as being impassable or unsafe area [26,27]. The use of aerial photographs and satellite images in the regional scale study of linear features such as faults, fracture zones, shears zones, and lithological contacts has greatly reduced the effects of these limitations [17,27].
Traditionally, lineament mapping is based on a visual or semi-automatic interpretation (photo-geology). Lineaments are often extracted manually by digitizing, which is subjective, time consuming, expensive and requires expertise, training and adequate scientific skills. In addition, it cannot produce results for large scale areas [28]. Lineament extraction could be more advantageous if the results were reproducible. Automatic lineament extraction is therefore needed.
Increasing spatial and radiometric resolution in satellite images favors the development of automatic, or criteria-based, lineament extraction algorithms [29,30]. However, algorithms are usually slow and commonly generate false features that are related to artificial linear features such as roads and power lines [31]. Linear feature enhancement and detection mostly have been done in the spatial domain (based on direct manipulation of pixels in an image) by using Sobel [32], Prewitt [32,33], and LOG filters [12,34–36], as well as morphological filters [37] and frequency domain (based on modifying the spectral transform of an image). The binary edge maps that have been derived in the edge detection step have been used as inputs for extraction algorithms, such as edge following (graph searching [29]), edge linking operators (standard Hough Transform [27,35,38–40], and edge tracing algorithms (STA, START, and ALERT algorithm [41]).
The success of automatic lineament extraction procedures depends on the reliability and accuracy of edge detection mechanism [28,42,43]. Until now it has been impossible to recognize or measure linear features of interest that have to be detected and to remove all remaining artifacts [43]. Selection of image preprocessing procedures and appropriate edge detection methods are very important in tectonic linear feature extraction, because they have significant effects on the accuracy of the final results.
The main goal of this study is to develop a new MATLAB based toolbox (TecLines) for automatic lineaments mapping from satellite images and digital elevation models (DEM). TecLines contains a set of functions for detecting and extracting potential edges with integration between image frequency and spatial filtering, Tensor voting framework, Hough transformation and polynomial fitting methodology.
The extraction of lineaments using the TecLines toolbox is performed in two main steps: edge detection and edge linking. Edge detection methods should ideally generate sets of pixels lying only on edges. In practice, these pixels seldom completely characterize edges because of unwanted noise and breaks. The general shape of the objects may be initially unknown, but in many cases they can be approximated by piecewise linear segments. It is not easily feasible to fit linear segments to all the edges in an image and discard false segments. These problems are addressed in the edge-linking step. The goal of edge linking is to describe an edge as a linear segment of specified shape and estimate the missing edge pixels from the assumed equation of the curvilinear segment. Due to differences in the mathematical background of the edge detection and edge linking procedure as well as the breadth of topics, we present the approach in a two-part paper. We chose MATLAB because it provides a multipurpose environment for mathematical processing based on a high-level programming language. Additionally, the functions provided in MATLAB are easy to modify and open for improvements.
The specific objective of the first part is to describe and assess the comprehensive edge detection procedure by integration of frequency domain filters, edge detection methods, and tensor voting framework. We introduce this part, with a focus on tectonic lineament extraction (binary edge maps with efficient thickness, length and pixel connectivity, and increased degree of accuracy of the edge detection). Finally, we validate this functionality using a synthetic image with known discontinuities and a high-resolution satellite image (QUICKBIRD-2) from an active tectonic area: the Andarab fault zone in NE Afghanistan.

2. Data

In this paper, we demonstrate the performance of the TecLines for edge detection, where validation has been performed on a synthetic and a real dataset.

2.1. Synthetic Dataset

In recent decades, several approaches for performance evaluation of edge detection methods have been proposed according to the presence or absence of ground truth data [44–50]. These approaches are based on the characteristic of the images (i.e., real images, synthetic images) [44,49–51]. Most of the edge detection methods that rely on ground truth use simple synthetic images [52–55], because it is easy to specify the ground truth edge locations [44,45]. In these cases, the edge detection can be quantitatively evaluated based on the known ideal detection considered to be the ground truth [55]. The synthetic DEM (Figure 1) used here is the result of landscape evolution algorithm created using set river incision and different uplift rates across tectonic faults. The drainage system adapts to the evolving surface conditions.

2.2. Real Dataset

Study Area and Data

We evaluated the performance of TecLines on satellite images of the active Andarab fault zone, northeastern part of Afghanistan (Figure 2a). The selected area has a long history of damaging earthquakes [56]. Furthermore, it encompasses regions interpreted as deformed by primarily transtensional forces in the Trans-Himalayan orogenic belt [57,58]. The Andarab fault is dextral and coincides with an approximately 150 km long, east-west-trending valley north of the intersection between the Paghman and Hari-Rud faults. In the high valley of Darya, many evidences of recent tectonic movements were observed along the fault trace [57]. In this study, we used panchromatic band of the Quickbird-2 (1 m spatial resolution) for 2 March 2006 (Figure 2b). This data is in UTM coordinate system, datum “WGS84” and zone 42 N.

3. Methodology

TecLines is a new MATLAB based framework that contains various functions for automatic detection and extraction of tectonic lineaments from satellite images and digital elevation models (DEM). Besides import and export functions that support the raster and vectors in standard file formats, TecLines provides functions for image filtering in the frequency and spatial domains to produce primary binary edge maps. Final binary edge maps in TecLines are produced by performing the computation of the Tensor voting framework. In addition, TecLines extract line segments from final binary edge maps by employing standard Hough transformation functions. A set of functions serves the grouping and merging line segments, which will be resulted in final lineaments maps. Comparing results with published/non-published lineament maps and lineament analysis are also possible in TecLines. An overview of TecLines divisions can be found in Figure 3.

3.1. Frequency Domain Filtering

Image enhancement methods can contribute significantly to the automated extraction of linear features by using noise reduction and edges enhancement in the input images. There are two main categories for image enhancement methods: (1) spatial domain methods, such as Mean, Median, and Mode filters, and (2) frequency domain methods, such as Gaussian and Butterworth high, low, and band pass frequency filters. The main difference between both categories is that the spatial domain methods are usually applied locally, while the frequency domain methods are usually done in the global context. Unfortunately, there is no common principle for determining what is “good” image enhancement. However, when image enhancement methods are used as a pre-processing phase for other image processing methods, then quantitative measures can determine which methods are most appropriate.
In this study, we used Butterworth band-pass filter in the frequency domain Equation (1) in order to enhance edges by suppressing low frequencies and filter noise by attenuating high frequencies [20]. The Butterworth band-pass method has several advantages with respect to other methods such as the smooth transfer function without any discontinuity or clear cut-off frequency to reduce the sharp truncation effect and the easy control of the range of frequencies passes (bandwidth) dependently on the order of the filter. In addition, the Butterworth band-pass method can improve the performance of the Ideal band-pass method while reducing the ring and blurring effect.
H B ( f ) = ( 1 1 + [ f / f L ] 2 n ) ⋅ ( 1 1 + [ f / f H ] 2 n )
where fL and fH are the low and high cut frequencies respectively, f is the distance from the origin, and n is the order of the filter. The Butterworth filter is available in MATLAB but it shows some limitations. According to MATLAB documentation, numerical problems can happen for frequency filter orders as low as 6 when the transfer function coefficient form ([b, a] form) is used. These problems are due to round-off errors. The reason is that higher order filters may need extremely precise tap coefficients to get the desired performance. Therefore, numerical error in the calculation of the tap coefficients can destroy the performance. In fact, the computed filter is unstable. Accordingly, we implemented a new MATLAB based code to apply the band-pass filter on the satellite images. This code allows edge enhancement by frequency domain filtering at the specified band-pass parameters such as the order of the filter and the upper and lower frequency cutoffs of the filter for a given bandwidth that satisfies the rate of transmission at the cutoffs while maintaining the desired performance.

3.2. Spatial Domain Filtering

The principal physical edges correspond to significant variation in the reflectance, illumination, orientation and depth of scene surface [59]. Conceptually, edge detection refers to the process of identifying and locating sharp discontinuities in an image in three steps: differentiation, smoothing and labeling [36,60–62]. Several types of methods are available for edge detection in the spatial domain. These methods are classified in: gradient based method (first order derivative), laplacian based methods (second order derivative) and optimal edge detection methods [59,63]. We successfully tested and implemented in TecLines three common methods: Sobel and Prewitt [32,33] (gradient based), LOG [36,64,65] (Laplacian based), and Canny [60] (optimum) edge detection methods. In the gradient-based edge detection methods, the magnitude of the gradient vector reflects the strength of the edge, or edge response, at any given point. The effect of the noise in the signal will appear on the small local maxima in magnitude (edge strength), thus the resulting map of local maxima is thresholding to distinguish convincing edges [66–68].
In the Laplacian based methods such as LOG method, the gaussian smoothing filter is used for decreasing the sensitivity to noise in the input image, in other words to slightly blur the image. Then, Laplacian is applied to detect regions of rapid intensity change. The disadvantages of these edge detectors are their sensitivity to noise also producing two-pixel thick edges [69].
The Canny edge detection method is an optimum method and provides a multi-stage algorithm to detect a wide-range of edges in images. Those should be at a minimum distance to the actual edge in the real image. In addition, the detected edges should have minimal response. In other words, the discovered edges should have only one response to a single edge and completely eliminate the possibility of multiple responses to an edge by using adaptive thresholds with hysteresis.

3.3. Morphological Image Processing

Binary edge maps produced by Sobel, LOG, and Canny edge detection methods can have small-scale edges and isolated (or island) edge pixels. A common approach for improving these results is to use mathematical morphology methods to eliminate extra edge pixels. In this study, we used opening morphological filter using bwareaopen command in MATLAB. Morphological opening is an erosion followed by a dilation. Erosion eliminates small-scale details by removing outlying pixels and isolated pixels. Dilation restores all remained edges to their original size. The opening method presents many advantages because it relies only on the relative ordering of the pixel value, not on their numerical values. The opening is anti-extensive. In addition, opening is idempotent operation because it can be applied multiple times without changing the result beyond the initial application.

3.4. Tensor Voting Framework

Tensor voting is a non-iterative methodology to the inference of statistically salient features from possibly sparse and noisy data. In the tensor voting framework [70], the input data is encoded as elementary second order, symmetric, non-negative definite tensors (position/orientation pairs), then support information (including proximity and smoothness of continuity) is propagated by vote casting According to these principles for each possible dimension, each input token, where the input token is an unstructured point cloud with no priori orientation, encodes a local potential orientation [71]. Therefore, every token is a location where an orientation is defined. The representation is in the form of a second order, symmetric, non-negative definite tensor which essentially indicates the saliency of each type of perceptual structure (curve, junction or region in 2-D) the token may belong to its preferred normal and tangent orientations [71,72]:
T = ( λ 1 − λ 2 ) e 1 e 1 T + λ 2 ( e 1 e 1 T + e 2 e 2 T )
The tensor can specify its preferred tangent, normal orientation, as well as saliency corresponding to its perceptual structures. In Figure 4, the major axis e1 is the preferred normal orientation of a potential curve segment. The magnitude of the stick component (λ1−λ2) indicates curve saliency.
Tensor voting was implemented using tensor fields to pre-compute and store the votes from both stick and ball voters in receivers at various distances and angles (Figure 5) [70–72]. Vote orientation corresponds to the smoothest local curve continuation from voter to the recipient, while vote strength decays exponentially with distance and curvature. The second order vote also is a stick tensor and has a normal lying along the circular arc from voter O to recipient P (Figure 5). According to the Gestalt principles, vote orientation corresponds to the smoothest local curve continuation between two points, O as a voter and P as a recipient, while vote strength decays exponentially with distance and curvature. The magnitude of the vote is described by a function of confidence in spherical coordinates that the voter O and the receiver P indeed belong to the same perceptual structure. This magnitude can be calculated according to [70]:
DF ( S , K , σ ) = e − ( S 2 + CK 2 σ 2 )
where S is the distance between voter (O) and receiver (P), K is the curvature, σ is the scale of voting, which determines the effective neighborhood size and C as a function of the scale controls the degree of decay with curvature (Figure 5). In the feature extraction process, a saliency map is produced for each feature type, which is assigning a second order symmetric tensor that estimates the structure of the feature type and the associated saliency. The characteristics of the elementary tensors in 2-D are shown in Table 1.
In TecLines, we wrote and implemented a set of MATLAB based functions for tensor voting to compute the gradient vector and the tensor matrix at the edge pixels of the binary image, which were extracted by edge detection methods. Then, voting tensor field generated by casting votes on all the edge pixels of the image were used to construct the stick saliency map. Finally, we extracted the distinguished edge pixels based on the stick saliency map.

3.5. Accuracy Measurements

There are specific statistical measures that we can use to assess erroneous measures [48,73,74]. The difference between edge pixels obtained by edge detection methods and a reference map (ground truth) can be inferred [75,76]. True positive is the number of correctly detected edge pixels, false positive is the number of pixels erroneously classified as edge pixels, and false negative is the amount of pixels that were not classified as edge pixels [77,78].
TP = | E IM ∩ E GTP | | P |
FP = | E IM ∩ − E GTP | | P |
FN = | − E IM ∩ E GTP | | P |
where TP is the total number of true positive edge. FN and FP are total number of false negative and false positive edges, respectively. EIM and EGTP are the total number of edge pixels in the image and reference dataset, respectively. P is the subset of positions in the image at which edges occur.
Accuracy requires that edges should be detected as close as possible to their correct positions. In a given image, the edge positions and numbers can be vary according to resolution and procedures [75]. In this study, by comparing the edges detected using TecLines with a reference map, the accuracy was computed as follows.
Ac = TP TP + FN + FP + TP E GTP 2 × 100
where Ac stands for accuracy. The value of Ac for accurate edge detection methods should be close to 100%. In order to evaluate the performance of the proposed edge detection procedure a reference dataset is required. The references dataset for both synthetic and QuickBird 2 images are determined based on manual extraction (Figure 6), because it is easy to distinguish geological and non-geological lineaments by visual interpretation [35,41,51,79,80].

4. Testing and Evaluating TecLines

4.1. Performance Evaluation of the Edge Detection Methods on a Synthetic DEM

Accuracy assessment is based on 21 known lineaments, which have been shown by 1180 edge pixels in the synthetic DEM (Figure 6b). The range of parameters to be used in edge detection procedure should be large enough to cover a wide range of detection results. We implemented 27 sets of parameter combinations of sigma, low and high thresholds for Canny edge detection method, where sigma ∈ (300, 500, 700), low ∈ (0.01, 0.05, 0.1), and high ∈ (1, 3, 5) In the opening morphological filter, minimum number of connected neighboring for each pixel was 4, 8, and 24 pixels. In the first step, we extract edges from original dataset without anterior Butterworth band-pass filtering (Figure 7a). The accuracy assessment shows that 30% of the known lineaments are detected as true positive edge pixels and 70% are detected as false negative edge pixels. The percentage of pixels that are false positives is about 45%. The tensor voting framework is applied to improve the result by considering the adjustment between increasing the information and decreasing the noise in the detection procedure. As seen in Figure 7b, the results after tensor voting framework demonstrate an increase in accuracy. 40% of the known lineaments are correctly detected. The percentage of pixels that are detected as false negative edge pixels is around 60%.
In the next step, we used a synthetic DEM preprocessed using a Butterworth band-pass filter in frequency domain. Due to local intensity changes in the image, edges are better detected in the west of the image (Figure 7c). The accuracy assessment of the Canny method shows that 78% of the known lineaments are detected as true positive edge pixels and 22% are detected as false negative edge pixel. Drainage patterns witness the critical relationship between streams, tectonic, and erosion changes [26,81,82]. Tectonic activity disrupts drainage networks. Studying the nature of this disruption can give clues about the magnitude and orientation of the original tectonic activity [83–85]. Thus, they are potential instruments for tectonic-geomorphology analysis [26,86]. The result of the visual analysis shows that nearly 60% of the streams are also detected as edge pixels. Due to noise, the results also include many fragmented edge pixels that lead to the extremely low accuracy.
As seen in Figure 7d, 93% of the known lineaments are correctly detected after performing tensor voting framework. The percentage of false negatives is less than 7%. However, these results consist of many discontinuous edges, instead of edge contours. Results show that employing a Butterworth band-pass filtering prior to Canny edge detector improves the its performance. Additionally, using the tensor-voting framework both improved binary edge extraction and merged neighboring edges with similar direction. The overall accuracy with the proposed approach is about 52% and is higher than the overall accuracy achieved by the other methods (Table 2). In particular, a comparison with true positive edge pixels directly validates the superiority of the implemented edge detection procedure based on frequency band-pass filtering, edge detection, and tensor-voting framework implemented in TecLines.

4.2. Performance Evaluation of the Edge Detection Methods on a Satellite Image

In the first step, we applied a Butterworth band-pass filter in the frequency domain. We used 0.2 as an upper-lower cutoffs frequency and 6 for order numbers of filter such that these preserved as many significant edges (tectonic linear features) present in the truth as possible (Figure 8). Frequencies above and below these cutoffs are in the stop-band. Therefore, TecLines cannot find it in the edge detection step. Clearly, using a wide bandwidth and high order number of filter can decrease this assumption and preserve small edges but also increase the number of the wrong edges resulting from noises. After filtering in frequency domain, we used Sobel, LOG and Canny methods to detect linear features (edges). We used a 3 × 3 pixel masking window size in order to ensure the detection of sufficiently small feature and preserve image details. We classified image pixels in two groups: edge points (marked as 0) and non-edge points (marked as 255). We used a fixed threshold of 60 for Sobel filter. We selected a sigma (standard deviation of the noise) of 1.2 and a threshold of 20 for LOG filter. We set the hysteresis thresholding to 0.25 and sigma to 0.08 for Canny filter. The number of edge points produced in thresholding procedure can be modified on a case-by-case basis. The preliminary binary edge maps produced by Sobel, LOG and Canny filters are shown in Figure 9.
We used opening morphological operation in order to thin the edges. This operation removes all pixels that are below a threshold of connected pixels. The minimum number of connected neighbors for each pixel was set to 4 pixels. The final binary edge maps are produced after removing all residual pixels.
We used the tensor-voting framework to improve edge detected map accuracy and merge neighboring edges with similar direction. This provides better results than other methods for detecting common edges. It should be noted that using tensor-voting framework could lead to loss of small edges surrounded by bigger ones. After the voting process, a saliency map is produced by assigning a second order symmetric tensor that estimates the structure of the feature type and the associated saliency (Figure 10). Then, edges are extracted based on that certain orientations of the features coexist at a given location (Figure 10). If we compare the extracted edges with and without the use of tensor voting framework it appears that the extracted edges become more apparent and continuous after tensor voting. This technique allows the production of longer and smoother edge lines than the common edge detection methods.

Accuracy Measurements

The accuracy was computed by comparing the edges detected using TecLines to the set of edges in the reference dataset (ground truth), The reference dataset is based on visual image interpretation and manual extraction. We analyzed the results before and after tensor voting framework. The results are shown in Table 3. This table shows that the tensor-voting framework significantly improves the true positive percentage in all of three edge detection methods. The number of false positive pixels is also drastically reduced. The highest overall accuracy was achieved by applying a combination of the Canny method and tensor voting framework and reaches about 74.5%. The accuracy of sobel and LOG methods were also improved after performing tensor-voting framework.

5. Conclusions

The aim of this work is to develop a MATLAB based toolbox (TecLines) in order to extract discontinuities from satellite images and digital elevation models. The high edge detection accuracy achieved by using TecLines shows that the proposed edge detection procedure can be used for geological purposes using high-resolution satellite images and DEM. The combined process of Butterworth band-pass filter in the frequency domain, edge detection methods in spatial domain and tensor voting framework slightly improves the edge detection accuracy. Consequently using the tensor-voting framework, the accuracy of edge detection is improved by 20% to 30%. We show that the combination of Sobel or LOG methods and tensor voting framework have similar overall accuracy for both the QuickBird and the synthetic images. However, the combination of the Canny and tensor-voting framework is considered more effective for edge detection. There are several important points that should be taken into account for further investigation. The percentage of achieved accuracy of edge detection methods largely depends on the gradient thresholds. The values of thresholds are still based on trial and error. For low thresholds, more edges will be detected but the resulting image will be increasingly susceptible to noise and trigger the detection of irrelevant features. A high threshold may miss subtle edges, or result in fragmented edges. Therefore, it is recommended to investigate further interactive thresholding in future edge detection methods. The outcome of this study is based on the panchromatic band of the QuickBird 2, which has a very high spatial resolution (1 m resolution). The proposed edge detection procedure is also adequate for medium resolution satellite images such as Landsat. The Landsat images have some advantages over QuickBird 2, such as free to availability, more coverage area, as well as shorter repeat intervals. While these results for edge detection procedures from satellite images are promising, the further mathematical approach presented in part 2 is still needed to extract linear segments from edges. In the second part of TecLines paper series, we describe the procedure of edge linking using the Hough transform and polynomial curve fitting using B-spline and Tavares methods. The quantitative representation of the binary edge maps can be exported as shape files, Geotiff and ENVI headers. Users can write their own MATLAB functions in order to expand the TecLines capabilities. The proposed toolbox is available from the TecLines website [87]. This provides new opportunities for education and research in tectonic lineament analysis.

Acknowledgments

High resolution QuickBird image was freely obtained from the GeoEye Foundation. We would like to thanks Louis Andreani for review and editing of this manuscript. Thanks to two anonymous reviewers whose input improved this paper.

Author Contributions

Mehdi Rahnama and Richard Gloaguen designed the research. Mehdi Rahnama wrote the manuscript and was responsible for the MATLAB programming, mathematical models and data analysis. Both authors contributed in editing and reviewing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Şener, Ş.; Sener, E.; Karagüzel, R. Solid waste disposal site selection with GIS and AHP methodology: A case study in Senirkent–Uluborlu (Isparta) Basin, Turkey. Environ. Monit. Assess 2011, 173, 533–554. [Google Scholar]
  2. Mukherjee, S. Microzonation of seismic and landslide prone areas for alternate highway alignment in a part of western coast of India using remote sensing techniques. J. Indian Soc. Remote Sens 1999, 27, 81–90. [Google Scholar]
  3. Charnpratheep, K.; Zhou, Q.; Garner, B. Preliminary landfill site screening using fuzzy geographical information systems. Waste Manag. Res 1997, 15, 197–215. [Google Scholar]
  4. Rangzan, K.; Charchi, A.; Abshirini, E.; Dinger, J. Remote sensing and GIS approach for water-well site selection, southwest Iran. Environ. Eng. Geosci 2008, 14, 315–326. [Google Scholar]
  5. Riedel, M.; Collett, T.S.; Kumar, P.; Sathe, A.V.; Cook, A. Seismic imaging of a fractured gas hydrate system in the Krishna–Godavari Basin offshore India. Mar. Pet. Geol 2010, 27, 1476–1493. [Google Scholar]
  6. Bartolome, R.; Gràcia, E.; Stich, D.; Martínez-Loriente, S.; Klaeschen, D.; de Lis Mancilla, F.; Iacono, C.L.; Dañobeitia, J.J.; Zitellini, N. Evidence for active strike-slip faulting along the Eurasia-Africa convergence zone: Implications for seismic hazard in the southwest Iberian margin. Geology 2012, 40, 495–498. [Google Scholar]
  7. Safari, H.O.; Pirasteh, S.; Pradhan, B.; Gharibvand, L.K. Use of remote sensing data and GIS tools for seismic hazard assessment for shallow oilfields and its impact on the settlements at Masjed-I-soleiman area, Zagros Mountains, Iran. Remote Sens 2010, 2, 1364–1377. [Google Scholar]
  8. Gomez, H.; Kavzoglu, T. Assessment of shallow landslide susceptibility using artificial neural networks In jabonosa River Basin, Venezuela. Eng. Geol 2005, 78, 11–27. [Google Scholar]
  9. Neuhäuser, B.; Terhorst, B. Landslide susceptibility assessment using “weights-of-evidence” applied to a study area at the Jurassic Escarpment (SW-Germany). Geomorphology 2007, 86, 12–24. [Google Scholar]
  10. Geiß, C.; Taubenböck, H. Remote sensing contributing to assess earthquake risk: From a literature review towards a roadmap. Nat. Hazards 2013, 68, 7–48. [Google Scholar]
  11. Pradhan, B.; Lee, S. Landslide susceptibility assessment and factor effect analysis: Backpropagation artificial neural networks and their comparison with frequency ratio and bivariate logistic regression modelling. Environ. Model. Softw 2010, 25, 747–759. [Google Scholar]
  12. Comparison of Edge Detection and Hough Transform Techniques for the Extraction of Geologic Features. Available online: http://www.cartesia.org/geodoc/isprs2004/comm3/papers/376.pdf (accessed on 19 June 2014).
  13. Mabee, S.B.; Hardcastle, K.C.; Wise, D.U. A method of collecting and analyzing lineaments for regional-scale fractured-bedrock aquifer studies. Ground Water 1994, 32, 884–894. [Google Scholar]
  14. Chowdary, V.M.; Ramakrishnan, D.; Srivastava, Y.K.; Chandran, V.; Jeyaram, A. Integrated water resource development plan for sustainable management of mayurakshi watershed, India using remote sensing and GIS. Water Resour. Manag 2009, 23, 1581–1602. [Google Scholar]
  15. Mabee, S.B.; Curry, P.J.; Hardcastle, K.C. Correlation of lineaments to ground water inflows in a bedrock tunnel. Groundwater 2002, 40, 37–43. [Google Scholar]
  16. Krishnamurthy, J.; Venkatesa Kumar, N.; Jayaraman, V.; Manivel, M. An approach to demarcate ground water potential zones through remote sensing and a geographical information system. Int. J. Remote Sens 1996, 17, 1867–1884. [Google Scholar]
  17. Rowan, L.C.; Lathram, E.H. Chapter 17 Mineral Exploration. In Remote Sensing in Geology; Siegal, B.S., Gillespie, A.R., Eds.; John Wiley and sons: New York, NY, USA, 1980; pp. 553–605. [Google Scholar]
  18. Lee, M.; Morris, W.; Harris, J.; Leblanc, G. An automatic network-extraction algorithm applied to magnetic survey data for the identification and extraction of geologic lineaments. Lead. Edge 2012, 31, 26–31. [Google Scholar]
  19. Agterberg, F.P.; Bonham-Carter, G.F.; Wright, D.F. Statistical Pattern Integration for Mineral Exploration. In Computer Applications in Resource Estimation Prediction and Assessment for Metals and Petroleum; Pergamon Press: Oxford, UK, 1990; pp. 1–21. [Google Scholar]
  20. Sabins, F.F. Remote sensing for mineral exploration. Ore Geol. Rev 1999, 14, 157–183. [Google Scholar]
  21. Halbouty, M.T. Application of Landsat imagery to petroleum and mineral exploration. AAPG Bull 1976, 60, 745–793. [Google Scholar]
  22. Cooper, S.M. Potential field investigation of the Liberia Basin, west Africa. J. Am. Sci 2010, 6, 199–207. [Google Scholar]
  23. Bartholomew, I.D.; Peters, J.M.; Powell, C.M. Regional Structural Evolution of the North Sea: Oblique Slip and the Reactivation of Basement Lineaments. In Geological Society, London, Petroleum Geology Conference Series, 1993; Geological Society of London: London, UK, 1993; pp. 1109–1122. [Google Scholar]
  24. Salvini, F.; Storti, F. Cenozoic tectonic lineaments of the Terra Nova Bay Region, Ross Embayment, Antarctica. Glob. Planet. Chang 1999, 23, 129–144. [Google Scholar]
  25. Zevenbergen, L.W.; Thorne, C.R. Quantitative analysis of land surface topography. Earth Surf. Process. Landf 1987, 12, 47–56. [Google Scholar]
  26. Shahzad, F.; Gloaguen, R. Tecdem: A matlab based toolbox for tectonic geomorphology, part 1: Drainage network preprocessing and stream profile analysis. Comput. Geosci 2011, 37, 250–260. [Google Scholar]
  27. Karnieli, A.; Meisels, A.; Fisher, L.; Arkin, Y. Automatic extraction and evaluation of geological linear features from digital remote sensing data using a hough transform. Photogramm. Eng. Remote Sens 1996, 62, 525–531. [Google Scholar]
  28. Wladis, D. Automatic lineament detection using digital elevation models with second derivative filters. Photogramm. Eng. Remote Sens 1999, 65, 453–458. [Google Scholar]
  29. Wang, J.; Howarth, P.J. Use of the hough transform in automated lineament. Geosci. Remote Sens. IEEE Trans 1990, 28, 561–567. [Google Scholar]
  30. Moore, G.K.; Waltz, F.A. Objective procedures for lineament enhancement and extraction (digital convolution enhanced images). Photogramm. Eng. Remote Sens 1983, 49, 641–647. [Google Scholar]
  31. Gustafsson, P. Spot satellite data for exploration of fractured aquifers in a semi-arid area in southeastern Botswana. Appl. Hydrogeol 1994, 2, 9–18. [Google Scholar]
  32. Duda, R.O.; Hart, P.E. Pattern Classification and Scene Analysis; Wiley: New York, NY, USA, 1973. [Google Scholar]
  33. Prewitt, J.M.S. Object Enhancement and Extraction; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  34. Morris, K. Using knowledge-base rules to map the three-dimensional nature of geological features. Photogramm. Eng. Remote Sens 1991, 57, 1209–1216. [Google Scholar]
  35. Suzen, M.L.; Toprak, V. Filtering of satellite images in geological lineament analyses: An application to a fault zone in central Turkey. Int. J. Remote Sens 1998, 19, 1101–1114. [Google Scholar]
  36. Marr, D.; Hildreth, E. Theory of edge detection. Proc. R. Soc. Lond. Ser. B Biol. Sci 1980, 207, 187–217. [Google Scholar]
  37. Tripathi, N.K.; Gokhale, K.; Siddiqui, M.U. Directional morphological image transforms for lineament extraction from remotely sensed images. Int. J. Remote Sens 2000, 21, 3281–3292. [Google Scholar]
  38. Fitton, N.C.; Cox, S.J.D. Optimising the application of the hough transform for automatic feature extraction from geoscientific images. Comput. Geosci 1998, 24, 933–951. [Google Scholar]
  39. Hough, P.V.C. Method and Means for Recognizing Complex Patterns, U.S. Patent 3069654. 18 December 1962.
  40. Duda, R.O.; Hart, P.E. Use of the hough transformation to detect lines and curves in pictures. Commun. ACM 1972, 15, 11–15. [Google Scholar]
  41. Koike, K.; Nagano, S.; Kawaba, K. Construction and analysis of interpreted fracture planes through combination of satellite-image derived lineaments and digital elevation model data. Comput. Geosci 1998, 24, 573–583. [Google Scholar]
  42. Argialas, D.P.; Harlow, C.A. Computational image interpretation models: An overview and a perspective. Photogramm. Eng. Remote Sens 1990, 56, 871–886. [Google Scholar]
  43. Karantzalos, K.; Argialas, D. Improving edge detection and watershed segmentation with anisotropic diffusion and morphological levellings. Int. J. Remote Sens 2006, 27, 5427–5434. [Google Scholar]
  44. Ziou, D.; Koukam, A. Knowledge-based assistant for the selection of edge detectors. Pattern Recognit 1998, 31, 587–596. [Google Scholar]
  45. Nguyen, T.B.; Ziou, D. Contextual and non-contextual performance evaluation of edge detectors. Pattern Recognit. Lett 2000, 21, 805–816. [Google Scholar]
  46. Heath, M.D.; Sarkar, S.; Sanocki, T.; Bowyer, K.W. A robust visual method for assessing the relative performance of edge-detection algorithms. IEEE Trans. Pattern Anal. Mach. Intell 1997, 19, 1338–1359. [Google Scholar]
  47. Cho, K.; Meer, P.; Cabrera, J. Performance assessment through bootstrap. Pattern Anal. Mach. Intell. IEEE Trans 1997, 19, 1185–1198. [Google Scholar]
  48. Lopez-Molina, C.; de Baets, B.; Bustince, H. Quantitative error measures for edge detection. Pattern Recognit 2013, 46, 1125–1139. [Google Scholar]
  49. Ji, Q.; Haralick, R.M. Efficient facet edge detection and quantitative performance evaluation. Pattern Recognit 2002, 35, 689–700. [Google Scholar]
  50. Román-Roldán, R.; Gómez-Lopera, J.F.; Atae-Allah, J; Martinez-Aroza, J.; Luque-Escamilla, P.L. A measure of quality for evaluating methods of segmentation and edge detection. Pattern Recognit 2001, 34, 969–980. [Google Scholar]
  51. Evaluation of Edge Detectors: Critics and Proposal. Available online: http://www.vision.auc.dk/~hic/perf-proc.html (accessed on 20 June 2014).
  52. Shin, M.C.; Goldgof, D.B.; Bowyer, K.W. Comparison of edge detector performance through use in an object recognition task. Comput. Vision Image Underst 2001, 84, 160–178. [Google Scholar]
  53. Heath, M.; Sarkar, S.; Sanocki, T.; Bowyer, K. Comparison of Edge Detectors: A Methodology and Initial Study, Computer Vision and Pattern Recognition, 1996. Proceedings of the 1996 IEEE Computer Society Conference on CVPR96, San Franciso, CA, USA, 18–20 June 1996; pp. 143–148.
  54. Pellegrino, F.A.; Vanzella, W.; Torre, V. Edge detection revisited. Syst. Man Cybern. Part B 2004, 34, 1500–1518. [Google Scholar]
  55. Yitzhaky, Y.; Peli, E. A method for objective edge detection evaluation and detector parameter selection. IEEE Trans. Pattern Anal. Mach. Intell 2003, 25, 1027–1033. [Google Scholar]
  56. Ambraseys, N.; Bilham, R. The Tectonic Setting of Bamiyan and the Seismicity in and Near Afghanistan for the Past 12 Centuries. In After the Destruction of Giant Buddha Statues in Bamiyan (Afghanistan) in 2001; Margottini, C., Ed.; Springer Berlin Heidelberg: Berlin, Germany, 2009; pp. 67–94. [Google Scholar]
  57. Map and Database of Probable and Possible Quaternary Faults in Afghanistan. Available online: http://pubs.usgs.gov/of/2007/1103/downloads/pdf/ofr07-1103plate.pdf (accessed on 20 June 2014).
  58. Wheeler, R.L.; Bufe, C.G.; Johnson, M.L.; Dart, R.L. Seismotectonic Map of Afghanistan, with Annotated Bibliography; US Department of the Interior, US Geological Survey: Reston, VA, USA, 2005. [Google Scholar]
  59. Ziou, D.; Tabbone, S. Edge detection techniques—An overview. Pattern Recognit. Image Anal 1998, 8, 537–559. [Google Scholar]
  60. Canny, J. A computational approach to edge detection. IEEE Trans. Pattern Anal. Mach. Intell 1986, PAMI-8, 679–698. [Google Scholar]
  61. Deriche, R. Using cannys criteria to derive a recursively implemented optimal edge detector. Int. J. Comput. Vision 1987, 1, 167–187. [Google Scholar]
  62. Ziou, D. Line detection using an optimal IIR filter. Pattern Recognition 1991, 24, 465–478. [Google Scholar]
  63. Maini, R.; Aggarwal, H. Study and comparison of various image edge detection techniques. Int. J. Image Process. (IJIP) 2009, 3, 1–11. [Google Scholar]
  64. Edge Detection Analysis. Available online: http://disp.ee.ntu.edu.tw/meeting/%E8%87%AA%E6%81%86/edge%20detection/edge_detection.pdf (accessed on 19 June 2014).
  65. Pratt, W.K. Digital Image Processing. In Image Enhancement; Wiley-Interscience: New York, NY, USA, 2007; pp. 247–307. [Google Scholar]
  66. Chai, H.Y.; Wee, L.K.; Supriyanto, E. Edge Detection in Ultrasound Images Using Speckle Reducing Anisotropic Diffusion in Canny Edge Detector Framework. Proceedings of the 15th WSEAS International Conference on SYSTEMS, Corfu Island, Greece, 14–16 July 2011; pp. 226–231.
  67. Gonzalez, R.C.; Woods, R.E.; Eddins, S.L. Digital Image Processing Using Matlab; Gatesmark Publishing: Knoxville, TN, USA, 2009. [Google Scholar]
  68. Han, C.; Kerwin, W.S.; Hatsukami, T.S.; Hwang, J.-N.; Yuan, C. Detecting objects in image sequences using rule-based control in an active contour model. IEEE Trans. Biomed. Eng 2003, 50, 705–710. [Google Scholar]
  69. Sonka, M.; Hlavac, V.; Boyle, R. Image Processing, Analysis, and Machine Vision; Springer: Berlin, Germany, 1999. [Google Scholar]
  70. Tensor Voting: Theory and Applications. Available online: http://www.sci.utah.edu/~gerig/CS7960-S2010/handouts/Medioni_tensor_voting.pdf (accessed on 19 June 2014).
  71. Mordohai, P.; Medioni, G. Tensor voting: A perceptual organization approach to computer vision and machine learning. Synth. Lect. Image Video Multimed. Process 2006, 2, 1–136. [Google Scholar]
  72. Loss, L.; Bebis, G.; Nicolescu, M.; Skurikhin, A. An iterative multi-scale tensor voting scheme for perceptual grouping of natural shapes in cluttered backgrounds. Comput. Vision Image Underst 2009, 113, 126–149. [Google Scholar]
  73. Abdou, I.E.; Pratt, W. Quantitative design and evaluation of enhancement/thresholding edge detectors. Proc. IEEE 1979, 67, 753–763. [Google Scholar]
  74. Kanungo, T.; Jaisimha, M.Y.; Palmer, J.; Haralick, R.M. A methodology for quantitative performance evaluation of detection algorithms. IEEE Trans. Image Proc 1995, 4, 1667–1674. [Google Scholar]
  75. Mokhtarian, F.; Mohanna, F. Performance evaluation of corner detectors using consistency and accuracy measures. Comput. Vision Image Underst 2006, 102, 81–94. [Google Scholar]
  76. Viola, P.; Jones, M. Rapid Object Detection Using a Boosted Cascade of Simple Features, CVPR 2001. Proceedings of the 2001 IEEE Computer Society Conference on Computer 2011 Vision and Pattern Recognition, Kauai, HI, USA, 8–14 December 2001; pp. I-511–I-518.
  77. Lu, D.-S.; Chen, C.-C. Edge detection improvement by ant colony optimization. Pattern Recognit. Lett 2008, 29, 416–425. [Google Scholar]
  78. Manish, T.I.; Murugan, D.; Kumar, G. Hybrid edge detection using canny and ant colony optimization. Commun. Inf. Sci. Manag. Eng. 2013, 3, 402–405. [Google Scholar]
  79. Koike, K.; Nagano, S.; Ohmi, M. Lineament analysis of satellite images using a Segment Tracing Algorithm (STA). Comput. Geosci 1995, 21, 1091–1104. [Google Scholar]
  80. Novak, I.D.; Soulakellis, N. Identifying geomorphic features using Landsat-5/TM data processing techniques on Lesvos, Greece. Geomorphology 2000, 34, 101–109. [Google Scholar]
  81. Mallast, U.; Gloaguen, R.; Geyer, S.; Rodiger, T.; Siebert, C. Derivation of groundwater flow-paths based on semi-automatic extraction of lineaments from remote sensing data. Hydrol. Earth Syst. Sci 2011, 15, 2665–2678. [Google Scholar] [Green Version]
  82. Kirby, E.; Whipple, K. Quantifying differential rock-uplift rates via stream profile analysis. Geology 2001, 29, 415–418. [Google Scholar]
  83. Gloaguen, R.; Kabner, A.; Wobbe, F.; Shahzad, F.; Mahmood, A. Remote Sensing Analysis of Crustal Deformation Using River Networks. In. Proceedings of the IEEE International Conference on Geoscience and Remote Sensing Symposium, IGARSS 2008, Boston, MA, USA, 7–11 July 2008; pp. IV-1–IV-4.
  84. Gloaguen, R.; Marpu, P.R.; Niemeyer, I.; Alonso, C.; Castellanos, M.T.; Redondo, J.M. Automatic extraction of faults and fractal analysis from remote sensing data. Nonlinear Process. Geophys 2007, 14, 131–138. [Google Scholar]
  85. Kurz, T.; Gloaguen, R.; Ebinger, C.; Casey, M.; Abebe, B. Deformation distribution and type in the Main Ethiopian Rift (MER): A remote sensing study. J. Afr. Earth Sci 2007, 48, 100–114. [Google Scholar]
  86. Shahzad, F.; Gloaguen, R. Tecdem: A matlab based toolbox for tectonic geomorphology, part 2: Surface dynamics and basin analysis. Comput. Geosci 2011, 37, 261–271. [Google Scholar]
  87. Tectonic Lineament Analysis using Remote Sensing Datasets Tollbox (TecLines). Available online: http://www.teclines.net (accessed on 19 June 2014).
Remotesensing 06 05938f1 1024
Figure 1. The synthetic DEM that is the result of landscape evolution algorithm created using set river incision and different uplift rates across tectonic faults. The drainage system adapts to the evolving surface conditions.
Figure 1. The synthetic DEM that is the result of landscape evolution algorithm created using set river incision and different uplift rates across tectonic faults. The drainage system adapts to the evolving surface conditions.
Remotesensing 06 05938f1
Remotesensing 06 05938f2 1024
Figure 2. (a) location of the study area in Northeast Afghanistan. (b) panchromatic band of the Quickbird-2 (1-m spatial resolution) for 2 March 2006, of the study area.
Figure 2. (a) location of the study area in Northeast Afghanistan. (b) panchromatic band of the Quickbird-2 (1-m spatial resolution) for 2 March 2006, of the study area.
Remotesensing 06 05938f2
Remotesensing 06 05938f3 1024
Figure 3. Overview of the essential components of lineament mapping using TecLines.
Figure 3. Overview of the essential components of lineament mapping using TecLines.
Remotesensing 06 05938f3
Remotesensing 06 05938f4 1024
Figure 4. Tensor decomposition.
Figure 4. Tensor decomposition.
Remotesensing 06 05938f4
Remotesensing 06 05938f5 1024
Figure 5. Second order vote casting by a stick tensor located at the origin [71]. Where S is distance between voter (O) and receiver (P), K is the curvature, σ is the scale of voting, which determines the effective neighborhood size and C as a function of the scale controls the degree of decay with curvature.
Figure 5. Second order vote casting by a stick tensor located at the origin [71]. Where S is distance between voter (O) and receiver (P), K is the curvature, σ is the scale of voting, which determines the effective neighborhood size and C as a function of the scale controls the degree of decay with curvature.
Remotesensing 06 05938f5
Remotesensing 06 05938f6 1024
Figure 6. (a) The reference lineament map for real dataset that is based on manual extraction from panchromatic band of QuickBird-2; (b) The reference map of the synthetic DEM consists in the digitized traces of the modeled faults (black line).
Figure 6. (a) The reference lineament map for real dataset that is based on manual extraction from panchromatic band of QuickBird-2; (b) The reference map of the synthetic DEM consists in the digitized traces of the modeled faults (black line).
Remotesensing 06 05938f6
Remotesensing 06 05938f7 1024
Figure 7. (a) and (b) preliminary binary edge map produced by Canny edge detection method and binary edge map obtained by tensor voting framework; (c) preliminary binary edge map produced by Canny edge detection method with prior Butterworth band-pass filter; (d) Final resulting binary edge map was obtained by tensor voting framework.
Figure 7. (a) and (b) preliminary binary edge map produced by Canny edge detection method and binary edge map obtained by tensor voting framework; (c) preliminary binary edge map produced by Canny edge detection method with prior Butterworth band-pass filter; (d) Final resulting binary edge map was obtained by tensor voting framework.
Remotesensing 06 05938f7
Remotesensing 06 05938f8 1024
Figure 8. (a) Quickbird-2 (1 m spatial resolution) for 2 March 2006 of the study area shown in pseudo color; (b) Filtered image after Butterworth band-pass filter. Colors represent the image brightness digital number (DN) values.
Figure 8. (a) Quickbird-2 (1 m spatial resolution) for 2 March 2006 of the study area shown in pseudo color; (b) Filtered image after Butterworth band-pass filter. Colors represent the image brightness digital number (DN) values.
Remotesensing 06 05938f8
Remotesensing 06 05938f9 1024
Figure 9. Preliminary binary edge maps resulting from Sobel (a), LOG (b) and Canny (c) methods and morphological filtering. We used a fixed threshold of 60 for Sobel, a sigma (standard deviation of the noise) of 1.2 and a threshold of 20 for LOG, and the hysteresis thresholding to 0.25 and sigma to 0.08 for Canny.
Figure 9. Preliminary binary edge maps resulting from Sobel (a), LOG (b) and Canny (c) methods and morphological filtering. We used a fixed threshold of 60 for Sobel, a sigma (standard deviation of the noise) of 1.2 and a threshold of 20 for LOG, and the hysteresis thresholding to 0.25 and sigma to 0.08 for Canny.
Remotesensing 06 05938f9
Remotesensing 06 05938f10 1024
Figure 10. (Left Pannel) Saliency density maps produced individually for (a) Sobel, (b) Log, and (c) Canny edge detection methods, respectively. The saliency density maps were constructed by using voting tensor fields that are generated from casting votes on all of the edge pixels. (Right Pannel) final extracted edge map by using tensor-voting framework after (a) Sobel, (b) Log, and (c) Canny filtering.
Figure 10. (Left Pannel) Saliency density maps produced individually for (a) Sobel, (b) Log, and (c) Canny edge detection methods, respectively. The saliency density maps were constructed by using voting tensor fields that are generated from casting votes on all of the edge pixels. (Right Pannel) final extracted edge map by using tensor-voting framework after (a) Sobel, (b) Log, and (c) Canny filtering.
Remotesensing 06 05938f10
Table
Table 1. Elementary tensors in 2-D, where n and t represent the normal and tangent vector respectively and the geometric features extracted after 2-D tensor voting [71,72].
Table 1. Elementary tensors in 2-D, where n and t represent the normal and tangent vector respectively and the geometric features extracted after 2-D tensor voting [71,72].
Featureλ1λ2e1e2TensorSaliencyNormalTangentNormalTensor
Point11Any orthonormal basisBallλ1 all2 > 1NoneAny orthonormal basisNone
Curve10n tStickλ1 − λ2 > λ2e1e2e1
Table
Table 2. Quantitative measures and overall accuracy obtained by TecLines for synthetic digital elevation models (DEM). True positive (TP) is the number of correctly detected edge pixels. False positive (FP) is the number of pixels erroneously classified as edge pixels. False negative (FN) is the amount of pixels that were not classified as edge pixels.
Table 2. Quantitative measures and overall accuracy obtained by TecLines for synthetic digital elevation models (DEM). True positive (TP) is the number of correctly detected edge pixels. False positive (FP) is the number of pixels erroneously classified as edge pixels. False negative (FN) is the amount of pixels that were not classified as edge pixels.
MethodTPFPFNOverall Accuracy (%)
Without Butterworth band-pass filteringCanny36026,37082015.9
Canny + TVF48017,38070021.6
With Butterworth band-pass filteringCanny92012,24026042.4
Canny + TVF109883398252.3
Table
Table 3. Quantitative measures obtained by TecLines for panchromatic band of QuickBird-2. True positive (TP) is the number of correctly detected edge pixels. False positive (FP) is the number of pixels erroneously classified as edge pixels. False negative (FN) is the amount of pixels that were not classified as edge pixels.
Table 3. Quantitative measures obtained by TecLines for panchromatic band of QuickBird-2. True positive (TP) is the number of correctly detected edge pixels. False positive (FP) is the number of pixels erroneously classified as edge pixels. False negative (FN) is the amount of pixels that were not classified as edge pixels.
MethodTPFPFNOverall Accuracy (%)
Sobel37268657246442.7
Sobel + TVF43331075185764.8
LOG38389721235243.1
LOG + TVF42711684191960.62
Canny41476293204346.5
Canny + TVF4891795129974.5

Share and Cite

MDPI and ACS Style

Rahnama, M.; Gloaguen, R. TecLines: A MATLAB-Based Toolbox for Tectonic Lineament Analysis from Satellite Images and DEMs, Part 1: Line Segment Detection and Extraction. Remote Sens. 2014, 6, 5938-5958. https://doi.org/10.3390/rs6075938

AMA Style

Rahnama M, Gloaguen R. TecLines: A MATLAB-Based Toolbox for Tectonic Lineament Analysis from Satellite Images and DEMs, Part 1: Line Segment Detection and Extraction. Remote Sensing. 2014; 6(7):5938-5958. https://doi.org/10.3390/rs6075938

Chicago/Turabian Style

Rahnama, Mehdi, and Richard Gloaguen. 2014. "TecLines: A MATLAB-Based Toolbox for Tectonic Lineament Analysis from Satellite Images and DEMs, Part 1: Line Segment Detection and Extraction" Remote Sensing 6, no. 7: 5938-5958. https://doi.org/10.3390/rs6075938

APA Style

Rahnama, M., & Gloaguen, R. (2014). TecLines: A MATLAB-Based Toolbox for Tectonic Lineament Analysis from Satellite Images and DEMs, Part 1: Line Segment Detection and Extraction. Remote Sensing, 6(7), 5938-5958. https://doi.org/10.3390/rs6075938

Article Metrics

Back to TopTop