Transmission tower deformation monitoring method based on frame interval time compression

Abstract

The invention discloses a transmission tower deformation monitoring method based on frame interval time compression, which adopts a multiple-input multiple-output synthetic aperture radar (MIMOSAR) to carry out deformation monitoring, and changes the form that a real interferometric phase is obtained by unwrapping data after deformation monitoring in a traditional method. The radar monitoring frame time without phase winding is obtained through finite element analysis, a longer-time monitoring result is obtained through accumulation of multi-frame short-time deformation phases, and winding is avoided from the precondition; during PS point selection, a six-step screening method for automatic identification of a power transmission tower cluster is adopted, the PS point cluster which belongs to a power transmission tower and is high and stable in amplitude and good in space connectivity can be automatically identified, and two kinds of PS point sets which can be used for atmospheric interference phase linear fitting and power transmission tower deformation analysis are obtained. According to the method, a block processing mode is adopted during data processing, computing resources are saved, and rapid analysis of deformation of the target point of the power transmission tower is finally achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of synthetic aperture radar, in particular to a MIMOSAR monitoring method for monitoring deformation of a power transmission tower. BACKGROUND

[0002] As one of the most important energies in modern society, electricity plays an irreplaceable role in all aspects of society. The stability of power grids and power critical facilities is not only the basis for the operation of modern society, but also directly related to the improvement of people's living quality and the sustainability of the overall high-quality development of the country. However, due to the influence of environmental and climatic factors, the structure of power critical facilities will inevitably produce a certain degree of deformation. When the deformation value exceeds the safety threshold, the power facilities will collapse, and the stable operation of the power grid will also be destroyed, causing huge economic and property losses to the country, society and people, and even threatening the lives and safety of relevant personnel. Therefore, it is necessary to monitor the deformation of power critical facilities, which is of great significance to the safe and stable operation of the power grid.

[0003] In recent years, SAR aperture synthesis technology has been widely used in the deformation monitoring of some large buildings and has achieved good monitoring results. SAR technology has the advantages of high precision, non-contact monitoring, less influence of bad weather, and the ability to monitor in low visibility environment, which is very suitable for deformation monitoring of power transmission towers under the action of wind. However, these technologies usually perform unwrapping processing on the data obtained after monitoring when processing data, and the power transmission tower is a truss structure with poor spatial continuity. When the deformation speed is too fast or the deformation span is large, it is easy to cause failure of unwrapping of the interference phase or large unwrapping error, thereby causing large errors in the final deformation results. SUMMARY

[0004] Based on the above, the present application provides a power transmission tower deformation monitoring method based on frame interval time compression, which can automatically identify the power transmission tower target and avoid phase wrapping to obtain real deformation data, so as to realize accurate unwrapping even in the case of poor spatial continuity of the monitoring target, and select high-quality deformation analysis points.

[0005] To achieve the above method, the present application adopts the following technical scheme: The power transmission tower deformation monitoring method based on frame interval time compression has the following steps: Step 1, set the radar parameters and receive the radar echo data to obtain the RA imaging data matrix set of the first frame to the Nth frame of the power transmission tower , and based on the RA imaging data matrix set obtain the interference phase matrix set between adjacent frames , wherein, Indicates the first Frame RA imaging data matrix, express and the Frame RA imaging data matrix Interference phase between them; Step 2: Use a six-step screening method to select from the RA imaging data matrix set Select permanent scattering point set from and permanent scattering point set of transmission towers ,in, express The Middle A permanent scattering point, express The Middle A permanent scattering point of a power transmission tower; Step 3, according to ,extract Linear fitting of the atmospheric disturbance phase model is performed on the points in the model to obtain the atmospheric disturbance phase model, and then the disturbance is removed. The atmospheric interference phase contained in the interference phase of the points in the middle is obtained. Deformation phase set corresponding to all points inside ; Step 4: Based on inter-frame deformation phase set Accumulate and acquire PS point sets Set of all point deformation values ;in express Inner The deformation of each point in the radar line of sight during the monitoring period.

[0006] The deformation monitoring method for transmission towers based on frame interval time compression described in this invention is characterized in that step 1 includes the following steps: Step 1.1: Based on the environmental factors of the transmission tower and the physical characteristics of the tower itself, use the finite element analysis method to determine the deformation value of the transmission tower that does not exceed [the specified value]. The time t1 is 4 / 4, where... Wavelength; Step 1.2: Set the number of radar data acquisition frames to N, the frame interval time to t = t1 / 2, and the sampling rate to N. f The number of single-pulse sampling points is m And the frequency modulation slope is k ; Step 1.3: Obtain the radar signal bandwidth using equations (1), (2), and (3) respectively. B Distance resolution Maximum detection range : (1) (2) (3) Step 1.4: Set the radar at its maximum detection range in front of the power transmission tower. Inside, radar echo data of the power transmission tower is obtained; Step 1.5: Based on radar echo data, generate the RA imaging data matrix set of the first to Nth frames of the transmission tower. ,in, This represents the RA imaging data matrix of the j-th frame. ; Step 1.6, based on Data in row a, column b and Data in row a, column b Equation (4) is used to obtain the interference phase of the a-th row and b-th column between adjacent frames. Thus obtain : (4) In equation (4), express The conjugate; This indicates the calculation of the phase function.

[0007] Furthermore, step 2 includes the following steps: Step 2.1: Use the correlation coefficient method to... Candidate scattering points are obtained by screening based on correlation coefficients; by Data in row a, column b Centered on, in Define a window of dimension m×n, and then use equation (5) to obtain and Correlation coefficient between Then, using equation (6) we obtain Mean correlation coefficient of data in row a, column b : (5) (6) when , then it means The data in row a and column b are candidate scattering points; among them, This represents the threshold value for the correlation coefficient of the data in row a and column b. Step 2.2: Perform amplitude threshold screening on the candidate scattering points to obtain quasi-scattering points; Calculate the average amplitude corresponding to each RA imaging data matrix in s, and denote it as the average amplitude set. ,in, express The corresponding average amplitude; from The minimum average amplitude is selected and denoted as . And used as the amplitude threshold, iterating through Is the amplitude of each candidate scattering point in the corresponding RA imaging data matrix greater than... If the value is greater than 0, the corresponding candidate scattering point is taken as the quasi-scattering point; otherwise, the corresponding candidate scattering point is discarded. Step 2.3: Use the amplitude dispersion index method to filter the scattering points using the amplitude dispersion index to obtain the permanent scattering point set. : Get The set of amplitude values ​​of the data in the a-th row and b-th column of each RA imaging data matrix ,in, express The magnitude of the data in row a, column b; Using equations (7), (8), and (9) respectively, obtain The average amplitude of the data in row a and column b Amplitude variance Amplitude dispersion index : (7) (8) (9) when When this happens, the quasi-scattering point is taken as the permanent scattering point and added. In, among them, The threshold value for the amplitude dispersion index of the data in row a and column b; Step 2.4, for The spatial connectivity of permanent scattering points within the tower is used to filter candidate transmission tower scattering points. Build Delaunay triangulation of all permanent scattering points within the region, statistical The number of triangle sides connected to any permanent scattering point within the area The number of permanent scattering points h within a certain length of its neighborhood; when and Then the corresponding permanent scattering point set As candidate scattering points for power transmission towers, among which... The threshold representing the number of triangle sides. The threshold for the number of permanent scattering points in the neighborhood; Step 2.5: Use the DBSCAN algorithm to cluster the scattering points of the candidate transmission towers to obtain PS point clusters. Then, based on the target physical characteristics and scattering characteristics of the transmission towers, identify the PS point clusters to obtain the PS point clusters of the transmission towers. Step 2.6: After sparsifying the PS point cluster of the transmission tower, the permanent scattering point set of the transmission tower is obtained. .

[0008] Furthermore, step 3 includes the following steps: Step 3.1: Based on the interferometric phase matrix set between adjacent frames middle Extract permanent scattering point set The set of inter-frame interference phases corresponding to each point in the data ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The interference phase generated between them; Step 3.2, based on The atmospheric delay coefficient is obtained by linearly fitting the traditional atmospheric disturbance linear phase model. and constant Thus, by using equation (10), we can obtain middle Atmospheric interference phase at position a, b : (10) In equation (10), for The line-of-sight distance of the target corresponding to the data in row a, column b; Step 3.3, based on middle Interference phase at row a, column b Phase with atmospheric interference Using equation (11) to obtain Deformation phase at row a, column b Thus obtain Deformation phase set corresponding to all points inside ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The deformation phase generated between; (11) Furthermore, step 4 includes the following steps: Step 4.1: After removing atmospheric interference phase errors, when the radar frame interval time setting meets the requirements... At that time, based on inter-frame deformation phase The set of interferometric phase matrices within the monitoring period is obtained using equation (12). Accumulated deformation phase at row a, column b : (12) Step 4.2: Based on cumulative deformation phase The RA imaging data matrix set within the monitoring period is obtained using equation (13). The cumulative deformation generated at position a, b from frame 1 to frame N at row a, column b Thus obtain : (13) The present invention provides an electronic device, including a memory and a processor, characterized in that the memory is used to store a program that supports the processor in executing the transmission tower deformation monitoring method, and the processor is configured to execute the program stored in the memory.

[0009] The present invention discloses a computer-readable storage medium storing a computer program, characterized in that the computer program, when executed by a processor, performs the steps of the transmission tower deformation monitoring method.

[0010] Compared with the prior art, the beneficial effects of the present invention are as follows: 1. When observing the deformation of transmission towers, this invention changes the traditional method of untangling the deformation monitoring data to obtain the true interference phase. Instead, it uses finite element analysis to obtain radar monitoring frame times that do not involve phase entanglement. By accumulating multiple short-time deformation phases, it obtains monitoring results over a longer period of time, thus avoiding entanglement under preconditions. Even if the spatial continuity of the monitored target is poor or the deformation is abrupt, this method can still achieve accurate untangling.

[0011] 2. In selecting PS points, this invention employs a six-step screening method for automatic identification of transmission tower clusters. This method integrates the correlation and amplitude information of pixels to obtain a set of PS points suitable for linear fitting of atmospheric interference phase. Considering the overall characteristics of transmission tower deformation, a Delaunay triangulation is constructed. Isolated points or points with poor spatial connectivity are removed based on the number of triangle edges connected by pixels and the number of PS points pre-screened in the neighborhood. The PS point density is reasonably adjusted. Finally, combined with the DBSCAN algorithm and the physical characteristics of the tower, PS point clustering and automatic extraction of PS point clusters that conform to the characteristics of transmission towers are achieved. This identifies a valid set of PS points that can be used for deformation result analysis, thereby identifying higher-quality valid target points and avoiding the selection of poor-quality target points or non-target points for atmospheric interference removal or deformation state analysis.

[0012] 3. Since this invention monitors by compressing frame time intervals, the amount of data is large. Data processing adopts a block processing method. This method has lower requirements for data processing equipment because the amount of data processed each time is much smaller than the amount of data in batch processing, which can greatly speed up radar data processing and provide deformation monitoring results for multiple target points. Attached Figure Description

[0013] Figure 1 This is a flowchart of the transmission tower deformation monitoring method based on frame interval time compression according to the present invention; Figure 2 This is a model diagram of the power transmission tower for this invention; Figure 3 This is a wind load time history diagram generated at a certain wind speed according to the present invention; Figure 4 This is a flowchart of the data block reading and processing of the present invention. Detailed Implementation

[0014] In this embodiment, a method for monitoring the deformation of transmission towers based on frame interval time compression is described, such as... Figure 1 As shown, it includes the following steps: Step 1: Determine the frame interval time parameter t.

[0015] Step 1.1: Before monitoring begins, determine the actual physical characteristics of the selected transmission tower, including its dimensions, material, damping, and connection relationships.

[0016] Step 1.2: After obtaining the actual physical characteristics of the transmission tower, perform physical modeling on it using finite element analysis software, such as... Figure 2 As shown.

[0017] Step 1.3, Calculation of wind load time history curve: The wind speed time history curve consists of mean wind and fluctuating wind. The mean wind can be obtained by converting the mean wind profile relationship using the exponential law, as shown in equation (1): (1) In equation (1), The wind speed at a height of 10m. The height that needs to be calculated. This is a coefficient related to ground roughness.

[0018] Pulsating wind is a random variable with zero mean. By using the traditional harmonic synthesis method for time history simulation, the time history curve of pulsating wind speed can be obtained. The wind speed time history curve is then obtained by superimposing it with the average wind at the corresponding location. According to Bernoulli's equation, the relationship between wind speed and wind pressure is shown in equation (2): (2) In equation (2), For wind load, This represents air density. Based on the above principles, a wind load time history curve at a given wind speed is obtained through conversion, as shown below. Figure 3 The figure shows the wind load time history curve generated at an average wind speed of 10 m / s.

[0019] The wind load at this time cannot be directly applied to the transmission tower. It is necessary to perform a certain conversion on the wind load applied to the transmission tower according to GB 50009 Code for Structural Loads of Buildings, as shown in Equation (3): (3) In equation (3), This refers to the wind load experienced by the transmission tower at a corresponding height. To obtain the wind load time history curve above, The tower shape coefficient can be found or calculated in GB50009 Code for Structural Loads of Buildings or the Technical Specification for Structural Design of Overhead Transmission Line Towers (DL / T 5154-2012). This refers to the projected area of ​​the components on the windward side of the power transmission tower.

[0020] Obtain the wind load at the corresponding height of the transmission tower. Then, the corresponding wind loads were applied to the transmission towers in the finite element analysis software. Analysis to obtain the deformation value of the transmission tower not exceeding The time t1 is 4 / 4, where... Given the wavelength, and to allow for a sufficient threshold, the frame interval is set to t = t1 / 2. The frame interval needs to be small enough, ranging from a few milliseconds to tens of milliseconds, so that the phase caused by atmospheric interference in the interference phase between adjacent frames is small enough or negligible.

[0021] Step 2: Set the radar parameters and receive radar echo data to obtain the RA imaging data matrix set of the power transmission tower from frame 1 to frame N. , RA imaging is an imaging algorithm that utilizes FFT to compress pulses from FMCW radar data to obtain high-resolution range images. It combines this with a digital beamforming (DBF) algorithm to calculate the steering vector and perform spatially matched filtering with the array data to obtain high-resolution angular images. This is based on a RA imaging data matrix set. Obtain the set of interferometric phase matrices between adjacent frames ,in, Indicates the first Frame RA imaging data matrix, express and the Frame RA imaging data matrix The interference phase between them.

[0022] Step 2.1: Set the number of radar data acquisition frames to N, the frame interval time to t = t1 / 2, and the sampling rate to N. f The number of single-pulse sampling points is m And the frequency modulation slope is k ; Step 2.2: To avoid exceeding the radar's monitoring range, the monitoring distance needs to be calculated based on the parameters set by the radar. The radar's signal bandwidth is obtained using equations (4), (5), and (6) respectively. B Distance resolution Maximum detection range : (4) (5) (6) Step 2.3: Set the radar to its maximum detection range in front of the power transmission tower. Inside, radar echo data of the power transmission tower can be obtained.

[0023] Step 2.4: Based on radar echo data, a data block processing method is adopted, reading only L frames of data at a time. Data is processed frame by frame in the L-frame data. Once processing is complete, the next L-frame data is processed. Figure 4 As shown; the final generated RA imaging data matrix set of the transmission tower from frame 1 to frame N is obtained. ,in, This represents the RA imaging data matrix of the j-th frame. .

[0024] Step 2.5, based on Data in row a, column b and Data in row a, column b The interference phase in the a-th row and b-th column between adjacent frames is obtained using equation (7). Thus obtain : (7) In equation (7), express The conjugate; This indicates the calculation of the phase function.

[0025] Step 3: Use the six-step screening method to select from the RA imaging data matrix set Select permanent scattering point set from and permanent scattering point set of transmission towers ,in, express The Middle A permanent scattering point, express The Middle A permanent scattering point of a power transmission tower; Step 3.1: Use the correlation coefficient method to... Candidate scattering points are obtained by screening based on correlation coefficients; by Data in row a, column b Centered on, in Define a window of dimension m×n, and then use equation (8) to obtain and Correlation coefficient between Then, using equation (9) we obtain Mean correlation coefficient of data in row a, column b : (8) (9) when , then it means The data in row a and column b are candidate scattering points; among them, This represents the correlation coefficient threshold of the data in row a and column b.

[0026] Step 3.2: Perform amplitude threshold screening on the candidate scattering points to obtain quasi-scattering points; Calculate the average amplitude corresponding to each RA imaging data matrix in s, and denote it as the average amplitude set. ,in, express The corresponding average amplitude; from The minimum average amplitude is selected and denoted as . And used as the amplitude threshold, iterating through Is the amplitude of each candidate scattering point in the corresponding RA imaging data matrix greater than... If the value is greater than 0, the corresponding candidate scattering point is taken as the quasi-scattering point; otherwise, the corresponding candidate scattering point is discarded.

[0027] Step 3.3: Use the amplitude dispersion index method to filter the scattering points by amplitude dispersion index to obtain the permanent scattering point set. : Get The set of amplitude values ​​of the data in the a-th row and b-th column of each RA imaging data matrix ,in, express The magnitude of the data in row a, column b; Using equations (10), (11), and (12) respectively, obtain The average amplitude of the data in row a and column b Amplitude variance Amplitude dispersion index : (10) (11) (12) when When this happens, the quasi-scattering point is taken as the permanent scattering point and added. In, among them, The threshold value for the amplitude dispersion index of the data in row a and column b.

[0028] Step 3.4, for The spatial connectivity of permanent scattering points within the tower is used to filter candidate transmission tower scattering points. Build Delaunay triangulation of all permanent scattering points within the region, statistical The number of triangle sides connected to any permanent scattering point within the area The number of permanent scattering points h within a certain length of its neighborhood; when and Then the corresponding permanent scattering point set As candidate scattering points for power transmission towers, among which... The threshold representing the number of triangle sides. This is a threshold for the number of permanent scattering points in the neighborhood.

[0029] Step 3.5: Use the DBSCAN algorithm to cluster the scattering points of the candidate transmission towers to obtain PS point clusters. At this time, the obtained point clusters include not only the transmission tower point clusters, but also point clusters formed by other non-monitored targets in the environment, such as trees and buildings. Therefore, it is necessary to further identify the PS point clusters based on the physical characteristics and scattering characteristics of the transmission tower targets to obtain the PS point clusters of the transmission towers. Step 3.6: Sparsify the points within the PS (Permanent Scattering Point) cluster of the transmission tower. The sparsification process is as follows: a PS point is only designated as a PS point if its amplitude is greater than the average amplitude of all PS points within a 3x3 region centered on it, or a multiple thereof. After sparsification, the permanent scattering point set of the transmission tower is obtained. .

[0030] Step 4, according to ,extract The interference phase at the points in the model is linearly fitted to the atmospheric interference phase model to obtain the atmospheric interference phase model, and then the interference is removed. The atmospheric interference phase contained in the interference phase of the point is obtained. Deformation phase set corresponding to all points inside ; Step 4.1: Based on the interferometric phase matrix set between adjacent frames middle Extract permanent scattering point set The set of inter-frame interference phases corresponding to each point in the data ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The interference phase generated between them.

[0031] Step 4.2, based on The atmospheric delay coefficient is obtained by linearly fitting the traditional atmospheric disturbance linear phase model. and constant Thus, by using equation (13), we can obtain middle Atmospheric interference phase at position a, b : (13) In equation (13), for The line-of-sight distance of the target corresponding to the data in row a and column b.

[0032] Step 4.3, based on middle Interference phase at row a, column b Phase with atmospheric interference Using equation (14) to obtain Deformation phase at row a, column b Thus obtain Deformation phase set corresponding to all points inside ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The deformation phase generated between; (14) Step 5: Based on inter-frame deformation phase set Accumulate PS point set Set of all point deformation values ;in express Inner The deformation of each point in the radar line of sight during the monitoring period.

[0033] Step 5.1: After removing atmospheric interference phase errors, when the radar frame interval time setting meets the requirements... At that time, based on inter-frame deformation phase The set of interferometric phase matrices within the monitoring period is obtained using equation (15). Accumulated deformation phase at row a, column b : (15) Step 5.2: Based on cumulative deformation phase The RA imaging data matrix set within the monitoring period is obtained using equation (16). The cumulative deformation generated at position a, b from frame 1 to frame N at row a, column b Thus, the point set is obtained. Set of all point deformation values : (16).

[0034] In this embodiment, an electronic device includes a memory and a processor. The memory stores a program that supports the processor in executing the above-described method, and the processor is configured to execute the program stored in the memory.

[0035] In this embodiment, a computer-readable storage medium stores a computer program, which is executed by a processor to perform the steps of the above method.

Claims

1. A method for monitoring the deformation of transmission towers based on frame interval time compression, characterized in that, Includes the following steps: Step 1: Set radar parameters and receive radar echo data to obtain the RA imaging data matrix set of the power transmission tower from frame 1 to frame N. , And based on the RA imaging data matrix set Obtain the set of interferometric phase matrices between adjacent frames ,in, Indicates the first Frame RA imaging data matrix, express and the Frame RA imaging data matrix Interference phase between them; Step 2: Use a six-step screening method to select from the RA imaging data matrix set Select permanent scattering point set from and permanent scattering point set of transmission towers ,in, express The Middle A permanent scattering point, express The Middle A permanent scattering point of a power transmission tower; Step 3, according to ,extract Linear fitting of the atmospheric disturbance phase model is performed on the points in the model to obtain the atmospheric disturbance phase model, and then the disturbance is removed. The atmospheric interference phase contained in the interference phase of the points in the middle is obtained. Deformation phase set corresponding to all points inside ; Step 4: Based on inter-frame deformation phase set Accumulate and acquire PS point sets Set of all point deformation values ;in express Inner The deformation of each point in the radar line of sight during the monitoring period.

2. The method for monitoring the deformation of transmission towers based on frame interval time compression according to claim 1, characterized in that, Step 1 includes the following steps: Step 1.1: Based on the environmental factors of the transmission tower and the physical characteristics of the tower itself, use the finite element analysis method to determine the deformation value of the transmission tower that does not exceed [the specified value]. The time t1 is 4 / 4, where... Wavelength; Step 1.2: Set the number of radar data acquisition frames to N, the frame interval time to t = t1 / 2, and the sampling rate to N. f The number of single-pulse sampling points is m And the frequency modulation slope is k ; Step 1.3: Obtain the radar signal bandwidth using equations (1), (2), and (3) respectively. B Distance resolution Maximum detection range : (1) (2) (3) Step 1.4: Set the radar at its maximum detection range in front of the power transmission tower. Inside, radar echo data of the power transmission tower is obtained; Step 1.5: Based on radar echo data, generate the RA imaging data matrix set of the first to Nth frames of the transmission tower. ,in, This represents the RA imaging data matrix of the j-th frame. ; Step 1.6, based on Data in row a, column b and Data in row a, column b Equation (4) is used to obtain the interference phase of the a-th row and b-th column between adjacent frames. Thus obtain : (4) In equation (4), express The conjugate; This indicates the calculation of the phase function.

3. The method for monitoring the deformation of transmission towers based on frame interval time compression according to claim 2, characterized in that, Step 2 includes the following steps: Step 2.1: Use the correlation coefficient method to... Candidate scattering points are obtained by screening based on correlation coefficients; by Data in row a, column b Centered on, in Define a window of dimension m×n, and then use equation (5) to obtain and Correlation coefficient between Then, using equation (6) we obtain Mean correlation coefficient of data in row a, column b : (5) (6) when , then it means The data in row a and column b are candidate scattering points; among them, This represents the threshold value for the correlation coefficient of the data in row a and column b. Step 2.2: Perform amplitude threshold screening on the candidate scattering points to obtain quasi-scattering points; Calculate the average amplitude corresponding to each RA imaging data matrix in s, and denote it as the average amplitude set. ,in, express The corresponding average amplitude; from The minimum average amplitude is selected and denoted as . And used as the amplitude threshold, iterating through Is the amplitude of each candidate scattering point in the corresponding RA imaging data matrix greater than... If the value is greater than 0, the corresponding candidate scattering point is taken as the quasi-scattering point; otherwise, the corresponding candidate scattering point is discarded. Step 2.3: Use the amplitude dispersion index method to filter the scattering points using the amplitude dispersion index to obtain the permanent scattering point set. : Get The set of amplitude values ​​of the data in the a-th row and b-th column of each RA imaging data matrix ,in, express The magnitude of the data in row a, column b; Using equations (7), (8), and (9) respectively, obtain The average amplitude of the data in row a and column b Amplitude variance Amplitude dispersion index : (7) (8) (9) when When this happens, the quasi-scattering point is taken as the permanent scattering point and added. In, among them, The threshold value for the amplitude dispersion index of the data in row a and column b; Step 2.4, for The spatial connectivity of permanent scattering points within the tower is used to filter candidate transmission tower scattering points. Build Delaunay triangulation of all permanent scattering points within the region, statistical The number of triangle sides connected to any permanent scattering point within the area The number of permanent scattering points h within a certain length of its neighborhood; when and Then the corresponding permanent scattering point set As candidate scattering points for power transmission towers, among which... The threshold representing the number of triangle sides. The threshold for the number of permanent scattering points in the neighborhood; Step 2.5: Use the DBSCAN algorithm to cluster the scattering points of the candidate transmission towers to obtain PS point clusters. Then, based on the target physical characteristics and scattering characteristics of the transmission towers, identify the PS point clusters to obtain the PS point clusters of the transmission towers. Step 2.6: After sparsifying the PS point cluster of the transmission tower, the permanent scattering point set of the transmission tower is obtained. .

4. The method for monitoring the deformation of transmission towers based on frame interval time compression according to claim 3, characterized in that, Step 3 includes the following steps: Step 3.1: Based on the interferometric phase matrix set between adjacent frames middle Extract permanent scattering point set The set of inter-frame interference phases corresponding to each point in the data ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The interference phase generated between them; Step 3.2, based on The atmospheric delay coefficient is obtained by linearly fitting the traditional atmospheric disturbance linear phase model. and constant Thus, by using equation (10), we can obtain middle Atmospheric interference phase at position a, b : (10) In equation (10), for The line-of-sight distance of the target corresponding to the data in row a, column b; Step 3.3, based on middle Interference phase at row a, column b Phase with atmospheric interference Using equation (11) to obtain Deformation phase at row a, column b Thus obtain Deformation phase set corresponding to all points inside ,in, express Inner The data at point j in the j-th frame With the data of frame j+1 The deformation phase generated between; (11)。 5. The method for monitoring the deformation of transmission towers based on frame interval time compression according to claim 4, characterized in that, Step 4 includes the following steps: Step 4.1: After removing atmospheric interference phase errors, when the radar frame interval time setting meets the requirements... At that time, based on inter-frame deformation phase The set of interferometric phase matrices within the monitoring period is obtained using equation (12). Accumulated deformation phase at row a, column b : (12) Step 4.2: Based on cumulative deformation phase The RA imaging data matrix set within the monitoring period is obtained using equation (13). The cumulative deformation generated at position a, b from frame 1 to frame N at row a, column b Thus obtain : (13)。 6. An electronic device, comprising a memory and a processor, characterized in that, The memory is used to store a program that supports the processor in executing the transmission tower deformation monitoring method according to any one of claims 1-5, and the processor is configured to execute the program stored in the memory.

7. A computer-readable storage medium storing a computer program thereon, characterized in that, When the computer program is run by the processor, it executes the steps of the transmission tower deformation monitoring method according to any one of claims 1-5.

Images