CN115585785B - Urban road sedimentation InSAR evaluation method - Google Patents

Abstract

The invention relates to an urban road sedimentation InSAR evaluation method, which comprises the steps of taking a time sequence InSAR result in a multi-period monitoring area, wherein the time sequence InSAR result comprises a monitoring point space position, a deformation rate and a gradual deformation amount; judging and screening deformation abnormal points in the deformation quantity; and (3) solving the spatial intersection of the road and the abnormal point by adopting a spatial analysis method, and selecting to obtain the abnormal point overlapped with the road. According to the urban road settlement InSAR evaluation method, inSAR software can identify and screen out points with mutation in deformation, then a spatial intersection of a road and the mutation points is obtained through a spatial analysis method, abnormal points overlapped with the road are selected, and road settlement conditions can be intuitively represented by using InSAR technology and used for evaluating the road state.

Description

InSAR evaluation method of urban road settlement

Technical field

The invention relates to the field of urban geological disaster prevention and control, particularly an InSAR evaluation method for urban road settlement.

Background technique

Synthetic Aperture Radar Interferometry (InSAR) is a developing microwave remote sensing technology with great potential. Synthetic aperture radar emits radar waves to illuminate the target object, receives reflected signals, and obtains information about the observed object. It is an active observation technology. InSAR technology is widely used in the monitoring of cities, important buildings, geological disasters, etc., and can provide deformation trend information of target monitoring objects. In recent years, Xu Qiang and others have proposed geological disaster identification methods based on the integration of space, ground and space, and based on the "three-shaped" information obtained by optical remote sensing, InSAR (Synthetic Aperture Radar Interference) and LiDAR (Lidar) multi-source remote sensing methods, disaster identification identification technology solutions, among which InSAR technology is an important basis for "situation" judgment. Time series InSAR uses multi-period images to remove interference information, which can greatly improve the accuracy of deformation monitoring. The two time series InSAR technologies currently widely used are SBAS-InSAR and PS-InSAR technology. Among them, SBAS-InSAR technology uses multi-period images to form interference pairs by combining them in pairs, selects interference pairs with high interference quality, and removes interference pairs with poor interference quality. Through multi-period interference, various types of noise phases within the observation time can be found. rules to remove them. PS-InSAR technology looks for so-called "permanent scatterer points" in multi-phase images, which are points that show small long-term changes in radar images, such as various structures, exposed rocks and other fixed objects. For these objects By conducting continuous observations, very high deformation monitoring accuracy can be obtained for permanent scatterers. Time series InSAR technology has been widely used in wide-area deformation monitoring work, such as goaf areas, urban buildings, and important infrastructure.

Time series InSAR technology can obtain two levels of information:

(1) Deformation rate. Figure 1 shows the deformation rate distribution map of Chengdu. The hierarchical coloring method can be used to reflect the spatial distribution of the deformation rate. The spatial distribution of the rate can be explained by mapping the deformation rate;

(2) Time series changes in deformation amount. InSAR technology can obtain the change in deformation amount in each period relative to the first period, that is, the change characteristics of deformation amount in the time dimension, and can also reflect the deformation characteristics of the monitoring target, as shown in Figure 2 .

Various applications mainly focus on points with abnormal deformation rates in InSAR monitoring results:

(1) The first type of point: a point with a large deformation rate, as shown in Figure 3;

(2) The second type of point: a point where the deformation suddenly undergoes large settlement after a certain point in time, as shown in Figure 4.

However, the existing technology only uses hierarchical color maps to map points with larger deformation rates (the first type of points) and clarify the spatial distribution characteristics of the second type of points. Since InSAR software mostly uses fitting deformation rates, The deformation rate of the second type of points is generally small and easy to be ignored. In practice, the sudden change of the deformation amount of the second type of points often corresponds to a sudden imbalance of surface stability, which may be caused by geological disasters such as road collapse and landslides. Early characterization.

Contents of the invention

The purpose of the present invention is to provide an InSAR evaluation method for urban road settlement in view of the problems existing in the existing technology.

In order to achieve the above objects, the technical solutions adopted by the present invention are:

An InSAR evaluation method for urban road settlement, including:

Step 1: Get the time series InSAR results in the monitoring area for multiple periods. The time series InSAR results include the spatial position of the monitoring point, the deformation rate, and the period-by-period deformation amount;

Step 2: There are a total of n periods of observation data. The time series InSAR data results of a single monitoring point include:

S1. Draw a diagonal line Line1 through the coordinates of the first period and the last period, and calculate the slope and intercept of Line1;

Assume that the coordinates of the deformation variable in the first and last periods are (x1, y1), (x2, y2), the 60% quantile of the deformation variable is f1, the 80% quantile is f2, and the 90% quantile is f3 ;

The slope of Line1: a=(y1-y2)/(x1-x2);

Intercept of Line1: b=(y2*x1-y1*x2)/(x1-x2);

S2. If a<0, set the single-period deformation variable coordinates as (x, y), calculate w=-a*x+y+b, if w>0, then the single-period deformation variable is on the right side of Line1, and the statistics are The total number of periods on the right side of Line1 is c;

S3. Set the distance from the point to Line1 as d, d=|(-a*x+yb)|/(a*a+b*b) ^1/2 , calculate the distance from all date points to Line1, and set the maximum distance is d _max , the maximum value corresponds to the mth period;

S4. Calculate the average deformation amount from the 1st period to the m-th period as u; calculate the average deformation amount from the m-th period to the last period as v;

S5. If c>0.8*n and u≤f1 and v<f2 and y2<f3, then the point is screened as an abnormal deformation point;

Step 3: Use the spatial analysis method to find the spatial intersection of the road and the abnormal points, and select the abnormal points that overlap with the road.

Adopting an urban road settlement InSAR evaluation method described in the present invention, the InSAR software can identify and screen out the points where the deformation value has sudden changes, and then use the spatial analysis method to find the spatial intersection of the road and the sudden change points, and select the points with the road Overlapping abnormal points can then be used to visually represent road settlement conditions using InSAR technology, which can be used to evaluate road conditions.

Preferably, the urban road settlement InSAR evaluation method also includes:

Step 4: Use field reconnaissance methods to look for signs of surface deformation at abnormal points.

Further preferably, the urban road settlement InSAR evaluation method also includes:

Step 5: Classify abnormal points according to urgency.

Preferably, the time series InSAR results are obtained through PSInSAR, SBASInSAR, StampsMTI or IPTA.

The present invention also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the urban road settlement InSAR evaluation method as described in any of the above items is implemented.

The invention also provides an electronic device, including:

A memory on which a computer program is stored;

A processor, configured to execute the program in the memory to implement the urban road settlement InSAR evaluation method as described in any one of the above.

In summary, due to the adoption of the above technical solutions, the beneficial effects of the present invention are:

The InSAR evaluation method of urban road settlement described in the present invention enables the InSAR software to identify and screen out points with sudden changes in deformation, and then uses the spatial analysis method to find the spatial intersection of the road and the sudden change points, and select the points that overlap with the road. The abnormal points can then be used to visually represent the road settlement situation using InSAR technology, which can be used to evaluate the road status.

Description of the drawings

Figure 1 is a schematic diagram of deformation rate graded coloring;

Figure 2 is a schematic diagram of the deformation variable time series curve;

Figure 3 shows the points with larger deformation rates in the InSAR monitoring results;

Figure 4 shows the point where the deformation variable suddenly changes in the InSAR monitoring results;

Figure 5 is a schematic flow chart of the InSAR evaluation method for urban road settlement;

Figure 6 is a schematic diagram of the time series InSAR data results of a single monitoring point;

Figure 7 is a schematic diagram of the selected abnormal deformation points.

Detailed ways

The present invention will be described in detail below with reference to the accompanying drawings.

In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention and are not intended to limit the present invention.

Example 1

As shown in Figures 5 to 7, an urban road settlement InSAR evaluation method according to the present invention includes the following steps:

Step 1: Get the time series InSAR (PSInSAR, SBASInSAR, StampsMTI, IPTA and other technologies) results in the monitoring area for 12 periods or more. The time series InSAR results mainly include the spatial position of the monitoring point, the deformation rate, and the period-by-period deformation amount.

Step 2: Screen out abnormal points of deformation

Obtain all the deformation rates in the monitoring area, and select points with a deformation rate less than 2% as abnormal points;

Screen the points where there is a mutation in the deformation variable;

Assuming n periods of observation data, as shown in Figure 6, the time series InSAR data results of a single monitoring point are as follows:

1. Make a diagonal line Line1 through the first and last period coordinates, and calculate the slope and intercept of Line1;

Let the coordinates of the deformation variable in the first and last periods be (x1, y1) and (x2, y2), and let the total number of periods be n; the 60% quantile of the deformation variable is f1, and the 80% quantile is f2. The 90% quantile is f3;

The slope of Line1: a=(y1-y2)/(x1-x2);

Intercept of Line1: b=(y2*x1-y1*x2)/(x1-x2);

2. If a<0, set the single-period deformation variable coordinates to (x, y), calculate w=-a*x+y+b, if w>0, then the single-period deformation variable is on the right side of Line1, and the statistics are The total number of periods on the right side of Line1 is c;

3. Set the distance from the point to Line1 as d, d=|(-a*x+yb)|/(a*a+b*b) ^1/2 , calculate the distance from all date points to Line1, and set the maximum distance is d _max , the maximum value corresponds to the mth period;

4. Calculate the average deformation amount from the 1st period to the m-th period as u; calculate the average deformation amount from the m-th period to the last period as v;

5. If c>0.8*n and u≤f1 and v<f2 and y2<f3, then the point is screened as an abnormal deformation point, as shown in Figure 7.

Step 3: Select the point that overlaps the road

The spatial analysis method is used to find the spatial intersection of roads and abnormal points, and the abnormal points that overlap with the roads are selected.

Step 4: Questionnaire about signs of deformation

Use field reconnaissance methods to look for signs of surface deformation at abnormal points and record them.

Step 5: Classify abnormal points

Classify the deformation variable abnormal points selected in step 2. The classification method is shown in Table 1.

Table 1. Abnormal point classification comparison table

The abnormal points obtained after preliminary screening are classified according to the principles in Table 2. According to the degree of urgency, they are divided into A: emergency; B, C: general. For points in category A, it is recommended to carry out on-site investigation and continue to pay attention to them. Category B and C: Points need to continue to pay attention to their deformation and development.

Table 2. Abnormal point classification comparison table

An urban road settlement InSAR evaluation method described in this embodiment is adopted, so that the InSAR software can identify and screen out the points where the deformation variable has a sudden change, and then use the spatial analysis method to find the spatial intersection of the road and the sudden change point, and select the The abnormal points of road overlap can then be used to visually represent the road settlement situation using InSAR technology, which can be used to evaluate the road status.

Example 2

The present invention is a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the urban road settlement InSAR evaluation method as described in Embodiment 1 is implemented.

Computer-readable storage media are used to store various types of data to support operations on the electronic device. These data may include, for example, instructions for any application program or method operating on the electronic device, as well as application-related data. ; Computer-readable storage media can be implemented by any type of volatile or non-volatile storage device or their combination, such as static random access memory (Static Random Access Memory, SRAM), electrically erasable programmable read-only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (Read-Only Memory, ROM), magnetic memory, flash memory, magnetic disk or optical disk.

Example 3

An electronic device according to the present invention includes:

A memory on which a computer program is stored;

A processor, configured to execute the program in the memory to implement the InSAR evaluation method for urban road settlement as described in Embodiment 1.

As a preferred solution of this embodiment, the electronic device may include: a processor and a memory. The electronic device may also include one or more of a multimedia component, an input/output (I/O) interface, and a communication component.

Wherein, the processor is used to control the overall operation of the electronic device to complete all or part of the steps in the above-mentioned urban road settlement InSAR evaluation method.

The memory is used to store various types of data to support operations on the electronic device. These data may include, for example, instructions for any application program or method operating on the electronic device, as well as application-related data; the memory may be composed of Any type of volatile or non-volatile storage device or their combination implementation, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (Electrically Erasable Programmable Read-Only) Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM) , magnetic memory, flash memory, magnetic disk or optical disk.

The multimedia component may include a screen and an audio component, where the screen may be a touch screen, for example, and the audio component may be used to output and/or input audio signals; for example, the audio component may include a microphone used to receive external audio signals, and the received audio signals It can be further stored in the memory or sent through the communication component; the audio component also includes at least one speaker for outputting the audio signal.

The I/O interface provides an interface between the processor and other interface modules. The other interface modules mentioned above can be keyboards, mice, buttons, etc.; these buttons can be virtual buttons or physical buttons.

The communication component is used for wired or wireless communication between the electronic device and other devices; wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G or 5G, or among them One or a combination of several, so the corresponding communication component may include: Wi-Fi module, Bluetooth module, NFC module, and mobile phone communication module.

As a preferred solution of this embodiment, the electronic device can be configured by one or more application specific integrated circuits (Application Specific Integrated Circuit, ASIC), digital signal processor (Digital Signal Processor, DSP), digital signal processing equipment (Digital Signal Processing) Device (DSPD), Programmable Logic Device (PLD), Field Programmable Gate Array (FPGA), controller, microcontroller, microprocessor or other electronic component implementation, used to perform the above InSAR evaluation method of urban road settlement.

In addition, the computer-readable storage medium provided by the embodiment of the present disclosure can be the above-mentioned memory including program instructions. The above-mentioned program instructions can be executed by the processor of the electronic device to complete the above-mentioned urban road settlement InSAR evaluation method.

The above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention shall be included in the protection of the present invention. within the range.

Claims (6)

1. An InSAR evaluation method for urban road settlement, which is characterized by including: Step 1: Get the time series InSAR results in the monitoring area for multiple periods. The time series InSAR results include the spatial position of the monitoring point, the deformation rate, and the period-by-period deformation amount; Step 2: There are a total of n periods of observation data. The time series InSAR data results of a single monitoring point include: S1. Draw a diagonal line Line1 through the coordinates of the first period and the last period, and calculate the slope and intercept of Line1; Assume that the coordinates of the deformation variable in the first and last periods are (x1, y1), (x2, y2), the 60% quantile of the deformation variable is f1, the 80% quantile is f2, and the 90% quantile is f3 ; The slope of Line1: a=(y1-y2)/(x1-x2); Intercept of Line1: b=(y2*x1-y1*x2)/(x1-x2); S2. If a<0, set the single-period deformation variable coordinates as (x, y), calculate w=-a*x+y+b, if w>0, then the single-period deformation variable is on the right side of Line1, and the statistics are The total number of periods on the right side of Line1 is c; S3. Set the distance from the point to Line1 as d, d=|(-a*x+yb)|/(a*a+b*b) ^1/2 , calculate the distance from all date points to Line1, and set the maximum distance is d _max , the maximum value corresponds to the mth period; S4. Calculate the average deformation amount from the 1st period to the m-th period as u; calculate the average deformation amount from the m-th period to the last period as v; S5. If c>0.8*n and u≤f1 and v<f2 and y2<f3, then the point is screened as an abnormal deformation point; Step 3: Use the spatial analysis method to find the spatial intersection of the road and the abnormal points, and select the abnormal points that overlap with the road. 2. The urban road settlement InSAR evaluation method according to claim 1, characterized in that it also includes: Step 4: Use field reconnaissance methods to look for signs of surface deformation at abnormal points. 3. The urban road settlement InSAR evaluation method according to claim 2, further comprising: Step 5: Classify abnormal points according to urgency. 4. The urban road settlement InSAR evaluation method according to any one of claims 1 to 3, characterized in that the time series InSAR results are obtained through PSInSAR, SBASInSAR, StampsMTI or IPTA. 5. A computer-readable storage medium with a computer program stored thereon, characterized in that when the program is executed by a processor, the urban road settlement InSAR evaluation method according to any one of claims 1-4 is implemented. 6. An electronic device, characterized in that it includes: A memory on which a computer program is stored; A processor, configured to execute the program in the memory to implement the urban road settlement InSAR evaluation method according to any one of claims 1 to 4.