Hierarchical Transmission Tower Detection from High-Resolution SAR Image

: The small scale of transmission towers and the environmental diversity around their situations make their detection in Synthetic Aperture Radar (SAR) images a challenging task. This paper presents a new hierarchical detection algorithm for transmission towers. First, Signal-to-Clutter Ratios (SCRs) of pixels are calculated. Afterwards, a SCR threshold is set. Since transmission towers possess strong scattering characteristics, pixels with SCRs above the threshold are considered as potential transmission tower pixels. Second, spatial densities of potential transmission tower pixels are calculated. According to the aggregation characteristics of transmission tower pixels, some potential transmission tower pixels with small spatial densities are removed. The remained potential transmission tower pixels are considered as candidate transmission tower pixels. The candidate transmission tower pixels are grouped by the nearest neighbour scheme such that in each group the distance between pixels is under a given threshold. Thus, each of the groups is viewed as a quasi-transmission tower. Convex-hulls of quasi-transmission towers are built, and then Minimum Bounding Rectangle (MBR) for each convex-hull is generated. According to the rectangle aspect ratios of MBRs, the real transmission towers are extracted. C-band HH-polarization GaoFen-3 (GF-3) amplitude images are used for experiments and four of the most popular transmission tower detection algorithms are selected as comparing algorithms to validate the proposed algorithms. The detection performance of transmission towers is evaluated with detection rate and quality factor. Experimental results verify that the proposed algorithm can efﬁciently and accurately detect transmission towers while maintaining the transmission tower geometry to a certain extent, which indicates that the proposed algorithm is efﬁcient and promising.


Introduction
As protection and support equipment for overhead transmission lines, transmission towers are important infrastructure for the power industry [1]. They are situated in a wide variety of natural environments such as the suburbs, farmland, forests and other field areas, and are very vulnerable to various natural disasters such as floods, landslide, wildfire, wind damage and so on [2]. To ensure the uninterrupted and safe transmission of electricity, effective detection and maintenance of transmission towers are essential [3,4]. Initially, transmission towers are mostly monitored by unmanned aerial vehicles (UAV) [5,6] and helicopter patrols to collect information about the status of transmission towers. Obviously, these methods consume a lot of manpower, material and financial resources and are prone to measurement errors. Therefore, the transmission tower detection based on remote sensing images is developed. Transmission towers are symmetrical in structure, unchanged for a long time and numerous in number. However, they are vertical structures that occupy a small area, so it is often difficult to detect transmission towers in remote sensing images. Synthetic Aperture Radar (SAR) system [7,8] can effectively avoid the influence of cloud, fog, rain, snow and other climatic conditions, and can continuously observe the ground surface for a long time [9,10], so it is widely applied in various fields such as target detection [11,12] and disaster monitoring [13,14]. However, SAR images have strong speckle noise [15][16][17], which increases the difficulty of the object detection task. With the continuous improvement of SAR system resolution, the information of transmission towers in the images becomes more and more abundant, so the detection of transmission towers from SAR images has gradually received the attention of many researchers [18].
Over the past decades, numerous transmission tower detection algorithms have been proposed. Depending on the data types, they are mainly divided into airborne or spaceborne SAR-based. In the early transmission tower detection, airborne SAR images are mainly used. For example, Woods et al. [19] proposed to construct Digital Surface Models (DSMs) using airborne InSAR [20] images and to detect vertical obstacles such as various towers based on the DSMs by modeling their geometries. However, this algorithm did not fully utilize the spectral information within the SAR images. Yang et al. [21] proposed an algorithm for automatic detection of transmission towers based on airborne fully polarized SAR images. This algorithm analyzes the four-component scattering model [22] and the spiral scattering component separately, extracts the point targets, and then connects all the isolated points as the final transmission tower lines. Currently, it is becoming a trend that spaceborne SAR images are used to detect transmission towers. Liu et al. [23] used the TerraSAR-X SSC images to detect and estimate the height of transmission towers. From the detection results of these algorithms, there still exist large false alarms.
Depending on the cardinal principles used, transmission tower detection methods can be divided into Constant False Alarm Rate (CFAR) [24,25] methods and deep learning methods [4,24,26]. For transmission tower detection in SAR images, the CFAR algorithm is one of the most commonly used algorithms, which has been used and studied for many years due to its computational simplicity and stable performance. For example, He et al. [27] proposed a local CFAR-based detection model to reconstruct transmission towers by local CFAR model, extract scale parameters, and eliminate coherent scattering noise to implement the detection algorithm. Therefore, Zhang et al. [28] combined CFAR with Extended Fractal (EF) to detect transmission towers as point targets in SAR intensity images. However, EF features have limited ability to describe geometric features, resulting in less accurate detection results. Subsequently, based on the characteristics of transmission towers. Guo et al. [29] proposed an algorithm to detect transmission towers using contrario theory. However, the algorithm is limited and needs to satisfy the equidistant linear arrangement, and the wrong detection is more serious. In addition, as an important research method of target detection, deep learning is also sought after by many scholars. Such as, Region-based CNNs (R-CNN) [30], Faster-RCNN [31], YOLO (You Only Look Once) [32] and SSD (Single Shot Multibox) [33] have been developed in the recent years for object detection tasks [26]. However, as the depth of the network increases, the convolutional neural network often encounters some negative effects, such as overfitting and gradient explosion, which affect the performance of the network itself. To this end, Gao et al. [34] proposed a deep cascaded network for change detection in SAR images. In addition, it is well known that speckle noise in SAR images has a great impact on detection. Zhang et al. [35] proposed a two-stage target-based deep learning method for multi-temporal SAR image change detection. Although this method reduced the false alarm rate to a certain extent, it can be found from the results that there was a loss of spatial details, over-reliance on the characteristics of deep learning, can only handle specialized tasks, and was not applicable to other complex tasks. Apparently, deep learning methods are mainly used for transmission detection by training a large number of transmission tower samples. Zeng et al. [36] used Polar Coordinate Semi-Variogram (PCSV) to extract the geometric features of Region of Interest (ROI) and transmission tower detection is implemented by using a three-layer neural network. However, the algorithm requires extensive preprocessing before the experiment, resulting in the detection of a large number of false alarm targets with similar geometric features to the transmission towers. Gao et al. [4] proposed an improved SSD-based transmission tower detection method to speed up the training by reducing the categories in the dataset and changing the scale and aspect ratio of default boxes in SSD, which improves the detection speed but tends to reduce the detection accuracy. In summary, from the results and basic ideas of these proposed algorithms, for large-scale images, all of them extract transmission towers as point targets, and the detection results have high false alarm rates or improve the detection rate through the idea of image enhancement and noise suppression [37]. Additionally, the algorithms that consider the geometric shape features of transmission towers ignore the amplitude information of transmission towers in SAR images.
This paper proposes a hierarchical algorithm from high-resolution SAR images. First, the Signal-to-Clutter Ratio (SCR) [38] features of the transmission towers in SAR images are used to extract the potential transmission tower pixels, which is different from the above literature that constructs an a priori model to detect the transmission towers. Secondly, based on the aggregation characteristics of transmission tower pixels, some potential transmission tower pixels with small spatial density are eliminated, and the remaining potential transmission tower pixels are considered as candidate transmission tower pixels. Finally, the convex-hulls [39] of quasi-transmission towers are generated by calculating the minimum distance between the candidate transmission tower pixel pairs, grouping the quasi-transmission tower pixel groups, and generating their minimum bounding rectangles (MBRs) [40]. The MBR aspect ratio extraction criterion is set according to the actual size of the transmission tower and its imaging morphology in the SAR images to complete the accurate extraction of transmission tower. Fast, robust and effective detection of transmission towers is realized to meet the detection needs of transmission towers.
The remainder of this paper is organized as follows. Section 2 describes the proposed algorithm in detail. Section 3 presents the experimental results of the proposed algorithm. Section 4 presents the comparison and analysis of the proposed algorithm and the comparative experimental results. Finally, Section 5 concludes this paper.

Preliminary
Since transmission towers are made of steel metal materials with stronger reflective properties than their surrounding background, and appear as strongly bright small elongated targets in high-resolution SAR images. Meanwhile, according to the geographical conditions on which transmission towers are situated, it can be known that the background around the transmission towers can be usually viewed as distributed targets, such as land and forest. Therefore, the transmission tower pixels have a stronger correlation each other compared to the background pixels in spectral and spatial. In conclusion, there are the following characteristics of transmission towers.
(1) The backscattering intensity of background clutter is required to be much lower than the peak signals of transmission towers, so it is reasonable to select potential transmission tower pixels by the SCRs.
(2) The spatial density of potential transmission tower pixels detected on a transmission tower area is greater than for ones on background area.
(3) The transmission towers in the SAR images appear as simple and recognizable long strips.

Potential Transmission Tower Pixels Detection
For a given preprocessed and cropped amplitude SAR image z = {z i (x i , y i ); (x i , y i ) ∈ D, i ∈ N}, where D is the image domain, i is the pixel index, (x i , y i ) is the grid position coordinate of pixel i, N = {1, . . . , n} is the pixel index set, n is the total number of pixels, and z i is the spectral reflection amplitude of pixel i. The idea behind potential transmission tower pixels detection corresponds to the transmission tower characteristics (1) in Section 2.1. That is, using transmission towers possessing strong scattering characteristics can distinguish potential transmission tower pixels with strong scattering characteristics from background pixels with weak scattering characteristics. In order to judge the scattering ability of pixels more accurately, it is necessary to convert the original SAR image to its SCR image, as the SCR image is more representative of the scattering ability by calculating the spectral values of central pixel and its neighboring pixels.
The definition of SCR is motivated by the spatial distribution of potential transmission tower pixels and the sliding detection window, as shown in Figure 1. Each transmission tower presents a number of aggregated pixels. The SCR value of the pixels within each window can be obtained by sliding the window. That is, first, an appropriate detection window W is established on the SAR image to be detected. By sliding the window, the pixels in the window are arranged in ascending order, the pixel with the largest amplitude is regarded as the target pixel, and the average value of the other smallest 90% of the pixels is regarded as the background pixel, thereby calculating the signal-to-noise ratio of the target pixel. The schematic diagram is shown in Figure 1. Let W be a detection window and its size be set according to the number of pixels occupied by the transmission tower and the SAR image resolution to minimize the detection of pseudo-pixels while ensuring that the transmission tower pixels are not missed. The detection window with too-large size will miss the transmission tower pixels, otherwise, it will result in too many pseudo pixels. Accordingly, for pixel i, its SCR is calculated as follows, where SCR i is the SCR value at pixel i, z imax is the maximum spectral reflection amplitude of pixels within the detection window central on pixel i, i.e., z i max = max{z i , i ∈ N i }, N i is the index set of pixel i and its neighboring pixels in the detection window, N i90% is the index set of pixels of the lowest amplitude 90% in the detection window, n i90% is the number of pixels in N i90% . Therefore, the generated SCR image can be expressed as To obtain the potential transmission tower pixels, the SCR threshold of the SCR image is set to T scr , i.e., the potential transmission tower pixels can be detected as SCR i ≥ T scr . Due to the different background complexity within the image, the statistical distribution of the pixels' amplitudes is prone to different situations such as single-peak distribution or multi-peak distribution, and the corresponding SCR images also show a variety of distributions, as shown in Figure 2a,b. In order to accurately describe the SCR image statistical distribution, a Gaussian mixture model is used, where f (·) is the probability density function (pdf) of the Gaussian mixture model, Φ k is the weight, µ k and σ k are the mean and variance of kth class, Φ k , µ k and σ k can be calculated by Expectation-Maximization (EM) algorithm. The threshold is chosen according to the Gaussian mixture model distribution at the wave crest corresponding to the maximum and next largest µ k . Without loss of generality, µ k 's , k = 1, . . . , K, are arranged in ascending order. Thus, it can be said that the SCR threshold is between K-1 and K (K > 2) means. To be accurate, SCR threshold can be calculated as follows, As a result, the set of indexes of potential transmission tower pixels can be obtained as

Candidate Transmission Tower Pixels Detection
According to the transmission tower characteristic (2) in Section 2.1, the spatial density of potential transmission tower pixels is higher than the spatial density of pixels in the background area. Therefore, candidate transmission tower pixels can be further filtered by calculating spatial density of the potential transmission tower pixel and based on a density threshold. Constructing a square window W i centered at potential transmission tower pixel i, and the spatial density of the W i can be expressed as for all potential transmission tower pixels, their spatial density set is where # denotes the operation of counting the number of elements in a set. As shown in Figure 3, the number of pixels contained in different windows varies. By calculating the pixel spatial density, the windows with higher density are detected by combining the predefined density threshold T ∆ to obtain the set of candidate transmission tower pixels. Taking the threshold T ∆ = 16 as an example, Figure 3 shows the calculation of the pixel spatial density for candidate transmission tower pixels, in which the interfering pixels in the two green square windows will be filtered out. Therefore, the index set of candidate transmission tower pixels can be expressed as C = {i, ∆ i ≥ T ∆ , i ∈ P}. By statistical analysis of the image, T ∆ can be obtained as follows: where · is the downward rounding symbol, n W i represents the total number of pixels in the window W i . The reason is that if they are transmission tower pixels, they should dominate within the n W i and correspond to a higher number of pixels, otherwise, smaller than T ∆ are background pixels.

Transmission Tower Detection
In order to remove false alarm targets in detected candidate transmission tower pixels, such as ground background, road edge lines and surrounding artificial targets, the geometric feature-based filtering algorithm is designed. Based on the imaging characteristics of transmission towers in SAR images, it is found that the constituent pixels of the same transmission tower will inevitably break and thus be divided into several subsets. Surrounding them are also interspersed with independent interfering pixels that are not removed by spatial density detection. Therefore, the candidate transmission tower pixel sets are first grouped by calculating the pixel distances and the interfering pixels will be removed as well.
The distance between the pixels of the candidate transmission towers is as, Setting the distance threshold as T G , candidate transmission tower pixels that satisfy the threshold range will be divided into one group, and each group of pixels is considered as a quasi-transmission tower. Therefore, the quasi-transmission tower pixels can be expressed as G = {i, ∀i ∈ C, ∃d ii < T G , i ∈ C}. The distance threshold T G can be calculated by the following: where · is the upward rounding symbol, β denotes the maximum break length between subsets of pixels in one transmission tower, which is determined by counting the transmission tower pixels within the image. L W i is the length of window W i . Furthermore, if the distances to any of the other pixels are larger than T G , it will be judged as a background interference pixel and filtered out. The schematic diagram of candidate transmission tower pixels grouped into several quasi-transmission towers is shown in Figure 4. Observing Figure 4, which contains two sets of candidate transmission tower pixels, A and B. The single pixel marked as #1 in A is considered as interference background pixels, so is #3 in B. Therefore, the two groups of pixels labeled #2 and #4 in A and B constitute two quasi-transmission towers, respectively. The proposed method can effectively filter out the independent interference pixels around the transmission towers.
Therefore, the quasi-transmission towers can be denoted as G = G j , j = 1, . . . , M , where G j is the set of quasi-transmission tower pixels corresponding to transmission tower j, satisfying U i=1 G j = G; G j ⊂ C, M is the number of transmission towers. Then, the convex-hull is generated for each quasi-transmission tower G j which can be denoted as H = H j , j = 1, . . . , M .
After that, a MBR is constructed for each convex-hull which is used to define the geometric shape of the transmission tower. The MBRs of convex-hulls can be denoted as R = R j ; j = 1, . . . , M . Figure 5a-c show the MBR corresponding to quasi-transmission tower, the MBR superimposed onto the original image and the sketch of the convex-hull.
Finally, the aspect ratio is determined in accordance with transmission tower geometry to further filter out quasi-transmission towers that do not satisfy the geometric parameters of transmission towers. That is, let is the coordinate of the center position of the rectangle, l i is the length of the rectangle, w i is the width of the rectangle, and θ i is the angle between the long side of the rectangle and the horizontal coordinate axis. Therefore, the ratio of length to width r j is, The aspect ratio thresholds are chosen based on the size of the transmission towers and the resolution of the images, and the MBRs that do not satisfy the aspect ratio are further filtered to obtain the exact MBRs. Additionally, the real transmission towers results can be denoted as T = T j , r j ≥ T r , j = 1, . . . , M , T r is threshold of r for real transmission tower.
The detailed procedures of the proposed Algorithm 1 are described as follows.

Algorithm 1. The detailed procedures of the proposed algorithm.
Input: Amplitude SAR image z. Output: Detection results T.
Step 1. Setting the window size for detecting pixels as W.
Step 2. Calculating the SCR by using Formula (1) and SCR threshold to obtain potential transmission tower pixels P.
Step 3. Calculating spatial density by using Formula (4) to obtain candidate transmission tower pixels C.
Step 4. Calculating pixels distances by Formula (6) and distance threshold to obtain quasi-transmission towers G.
Step 5. Generating convex-hulls H of G j , constructing MBRs R of the convex-hulls and set aspect ratio of MBRs by Formula (8).
Step 6. Obtaining the detection results T. Therefore, the overall flowchart of the proposed method is shown in Figure 6.

Data Used
In the experiment, C-band HH-polarization GF-3 amplitude images are employed to validate the proposed algorithm. Figure 7

Experimental Results and Analysis
To validate the feasibility and effectiveness of the proposed algorithm, transmission towers detection experiments in real SAR remote sensing images are performed on a PC with Intel (R) Core (TM) i5-4460, 3.20 GHz and 8 GB memory using MATLAB R2020b. Figure 8 shows the normalized SCR image results, with white pixels the high SCR region and black pixels the low SCR region. As seen in Figure 8, the detail information within the SCR image is stronger than that of the original image. Although the transmission tower is hidden in a large amount of background clutter, it can be seen using Formula (1) that the proposed algorithm fully considers the relationship between transmission tower pixels and neighboring pixels, and has a stronger SCR value compared to other white background pixel clutter.   Figure 8a-d, respectively. It can be seen from Figure 10a-d, for each transmission tower, there are always part of its pixels being detected, although there are still some false alarm pixels in the potential transmission tower pixels. Figure 11a-h show the results using gray threshold method with different threshold in Figure 7a-d, respectively. Comparing the detection results of the two methods, it is found that the detection using SCR images is more advantageous compared to the original images. A large number of background interference pixels are removed while ensuring that transmission towers are not missed, and the geometry of transmission tower pixels is preserved. Taking the transmission tower in Figure 10a as an example, a red box is used to mark the corresponding location of the transmission tower and it is zoomed in to show that the pixel distribution and geometry of the transmission tower can be clearly observed without losing important component pixels. However, when the threshold is increased, most of the transmission tower pixels are filtered out, leaving only a few main pixels, as shown in Figure 11e-h. Among them, by observing Figure 11f, it can be found that the geometry of the transmission tower cannot be identified, and a transmission tower is missed in Figure 11h. Therefore, a comparison with Figure 11 indicates that the proposed algorithm is not simply and irregularly filtered by the amplitude values. By calculating the spatial densities of potential transmission tower pixels, candidate transmission tower pixels can be obtained by thresholding the densities. Observing the detection results in Figure 10, the transmission tower pixels have obvious spatial aggregation. According to Formula (5), the density threshold ranges of Figure 10a-d are set to 18-66. Therefore, Figure 12a-d show the candidate transmission tower pixels by thresholding the densities of potential transmission tower pixels in Figure 10a-d, respectively. Observing Figure 12, it can be noticed that most of the false alarm pixels have been filtered out and only a few pixels are detected as candidates for the number of transmission tower pixels. Additionally, all the transmission tower main pixels within the four images have also been retained, which verifies the effectiveness of the density thresholding algorithm. The results of grouping the candidate transmission tower pixels by the nearest neighbor scheme are shown in Figure 13a-d. By calculating the distance between pixel pairs, several sets of pixels are determined and each set of pixels is considered as a quasitransmission tower. Figure 13a-d contain 7, 3, 6 and 11 quasi-transmission towers, respectively. Figure 14a Table 1. (It is worth noting that in this paper, based on the imaging characteristics of transmission towers, the default horizontal axis is the length and the vertical axis is the width.) Among them, the aspect ratio thresholds of Figure 14a-d are set to 1.0, 3.5, 1.2, 3.5, respectively. Therefore, the non-transmission towers are further filtered to obtain the real transmission tower detection results, as shown in Figure 15a-d, and the real and the false alarm transmission towers are marked with red rectangles and yellow triangles, respectively. Observing Figure 15a-d, most of the false alarm targets are eliminated, and it can be seen that the detection results do not show any missed detection, only Figure 15c shows one false alarm. Overall, as expected, the proposed algorithm can detect the transmission towers well, so as to better improve the detection rates. Figure 16 shows the effect of superimposing the final transmission tower detection results onto the original image, which can better show the detection effect of the proposed algorithm.  Figure 14a,c for each quasi-transmission tower aspect ratio.

Discussion
For comparison, some known and proven detection algorithm are utilized. In this paper, Log-normal distribution of CFAR detection algorithm, the combined CFAR and EF algorithm [28], the combined CFAR and Contrario algorithm [29] and deep learning method [35] are used as comparison algorithms to confirm the detection performance of the proposed algorithm. Figure 17 shows results from the comparison algorithms for transmission towers detection in Figure 7a-d. Among them, Figure 17(a1)-(d1) show the results of Log-normal distribution of CFAR detection algorithm. Figure 17(a2)-(d2) show the results of the combined CFAR and EF algorithm. Figure 17(a3)-(d3) show the results of and the combined CFAR and Contrario algorithm. Additionally, Figure 17(a4)-(d4) show the results of the proposed algorithm. In Figure 17, the red and yellow rectangle marks the real and false alarm transmission towers, respectively, and the green oval marks the missed detection transmission towers.
As seen from the first row of Figure 17, the log-normal distribution of CFAR detection algorithm is simple and can describe the asymmetrical distribution of tailing phenomenon, it is not suitable for the identification of small targets under complex scenes. For example, in experimental images shown in Figure 7a,c with relatively complex background, there are a large number of false alarm transmission towers in the detection results, resulting in wrong detection. Compared with the log-normal distribution of CFAR detection algorithm, the false alarm transmission towers detected by the combined CFAR and EF algorithm are significantly reduced, but the EF feature has a weak ability to describe the geometric features of targets. As a result, the detection results cannot effectively distinguish transmission towers from non-transmission towers. The combined CFAR and Contrario algorithm is based on the characteristic of being linearly aligned to detect transmission towers, and it can be seen from Figure 17(a3)-(d3) that the detection results have a large number of false alarm targets, and even the alignments of false alarm targets. Obviously, the algorithm first filters out some candidate points and then detects the transmission towers by using the feature that they present a linear arrangement. However, many false alarm targets can easily be detected by the linear arrangement, and the arrangement of real transmission towers is very demanding, for example, they must be equally spaced and arranged in a straight line, and will be missed if they are turned around.
The proposed algorithm not only accurately identifies the transmission towers, but also preserves the geometric characteristics of the transmission towers. Among them, only one false alarm transmission tower appears in Figure 17(c4). Moreover, it can be observed that the size and geometry of the false alarm target is similar to that of the transmission tower, which is the reason why the false alarm target is mistakenly detected as a transmission tower.
In addition, in order to verify the superiority of the proposed algorithm, the currently popular deep learning method is used for experimental comparison. As shown in Figure 18a-d, are the detection results of YOLO v5, and Figure 18e-h show the results of the proposed algorithm superimposed on the original images, so as to show the comparison of the results of the two algorithms more intuitively. Observing Figure 18a-d, it can be clearly found that the detection results of the deep learning method are poor, and there are a large number of missed detections and false detections in the four images. The reason is that the detection results are not ideal due to the lack of enough transmission tower training samples and public data sets. Moreover, the small size of the transmission tower makes it difficult to judge its accurate geometry in SAR images with large speckle noise, which is also the reason for the poor results of deep learning methods.
To quantitatively evaluate the detection performance of the five algorithms, this paper uses the detection rate β Rate , quality factor δ Quality and the number of false detections n false to evaluate them, β Rate = n Correct n Real (9) δ Quality = n Correct n Real + n false (10) n false = n missing + n false-alarm (11) where n Correct is the correct number of detections, n Real is the real number of transmission towers, n missing is the number of missed detections, and n false-alarm is the false alarm target. The higher the detection rate, the better the ability to detect really transmission towers, and the higher the quality factor, the more sensitive the detection algorithm is in discriminating between transmission towers and non-transmission towers. Table 2 presents the detection results of the five algorithms. Comparing the indicators in Table 2, the following conclusions can be drawn: (1) The other three algorithms have high detection rates , except for the combined CFAR and Contrario algorithm and YOLO v5 method. In the experiment, it is found that the combined CFAR and Contrario algorithm is more suitable for detection when the background is simple and the number of transmission towers is large. If the number of transmission towers is small, it is easy to cause other false alarm targets that are linearly arranged to be misdetected. In addition, transmission towers are usually erected based on the local geological and environmental conditions, so it is not guaranteed that all transmission towers are kept evenly distributed and equally spaced. This case leads to the low detection rate of the combined CFAR and Contrario algorithm. Furthermore, the deep learning method needs to train a large number of samples, which is time-consuming and labor-intensive. In the case of lack of sufficient samples and few targets in the image, it is easy to lead to low detection accuracy. Moreover, the small size of the transmission tower brings great difficulties to the training results. (2) Among the five algorithms, only the proposed algorithm has a higher accuracy of quality factor. Due to the existence of a large number of false alarm targets, the complexity of the scene increases, such as Figure 7a,c, especially the log-normal distribution of CFAR detection algorithm and the combined CFAR and Contrario algorithm, so the number of false detections increases and subsequently the quality factor decreases to an unsatisfiable accuracy range. In contrast, the proposed algorithm has a small number of false detections, so the quality factor accuracy remains high in all four experimental images. This is not the case for one false alarm in Figure 17(c4), which causes the quality factor to drop to 83.33%. After comparison, it is found that proposed algorithm has the strongest detection ability compared with the other four algorithms, and the identification results for transmission and non-transmission towers are more accurate.
Computational complexities are measured by the algorithms' CPU time. Table 3 presents the running time of different algorithms, as shown in Figure 7a. For the comparison algorithms, the Log-normal distribution of CFAR detection algorithm has the shortest running time for the same image size, followed by the proposed algorithm and deep learning method, while the combined CFAR and EF algorithm and the combined CFAR and Contrario algorithm have a relatively longer running time. As the EF algorithm needs to consider both the intensity and geometric features of the target, it takes a longer time to run. Moreover, the Contrario algorithm needs to filter all point targets detected by CFAR detection algorithm, which takes more time. In addition, the combined CFAR and Contrario algorithm has a much greater computational complexity than other algorithms. On the other hand, the hierarchical detection process leads to a longer time for the proposed algorithm to perform the stepped thresholding. However, the computational complexity is much lower than the other four algorithms, and does not need to train a large number of samples.

Conclusions
The main contribution of this paper is to provide a new method for hierarchical detection of transmission towers based on SAR images.
The method is simple, fast, it has high precision, and fully combines various characteristics of the transmission tower to achieve high-precision detection results. Compared with the other known and proven techniques, the proposed algorithm has the following advantages. (1) The idea of combining high SCR with geometric features is fully utilized to hierarchical detection, without relying too much on any one step. (2) There is no need to do much preprocessing work and to consider a priori knowledge. (3) The proposed algorithm can better demonstrate the imaging morphology of transmission towers within SAR images. In addition, comparative experiments are conducted to verify the superiority of the proposed algorithm. Although the proposed algorithm is effective in distinguishing transmission and non-transmission towers, it has some shortcomings. For example, a large number of small targets with similar geometric and scattering characteristics to transmission towers within the image are not accurately detected and classified. In future work, the study will focus on detecting targets with characteristics similar to transmission towers and perform more accurate analysis.