InSAR measurement method, device and storage medium based on temporary coherent scatterer

Abstract

An InSAR measurement method and device based on temporary coherent scatterers, and a storage medium. The invention discloses an InSAR measurement method and device based on temporary coherent scatterers. The invention obtains a coherent phase matrix by interfering every two SAR images in a SAR image sequence. The temporary coherent pixels in the SAR image are identified based on the coherent phase matrix. The phase unwrapping of the identified temporary coherent pixels is performed to obtain a deformation sequence of the temporary coherent pixels, thereby realizing the monitoring of the deformation amount of the temporary coherent pixels and obtaining the InSAR measurement result of the temporary coherent scatterer.

Description

Technical Field

[0001] The present invention relates to the field of InSAR technology, and in particular to an InSAR measurement method and device based on temporary coherent scatterers, and a storage medium. Background Art

[0002] How to discover and effectively identify the potential risks of major geological disasters in advance and take proactive measures to prevent and control them has become the focus and difficulty of the recent geological disaster prevention and control field. The subsidence of reclamation areas also threatens the safety and stability of buildings, public facilities and other structures.

[0003] Among them, the measurement of deformation in geological areas is the core technology for early discovery and effective identification of major geological disasters. Commonly used deformation measurement methods include total station observation method, level settlement observation method, etc. These measurement methods usually require on-site setting of measurement marks and regular manual on-site operations, which require a lot of manpower and material resources. The Beidou / GNSS observation technology, which has gradually matured, also requires on-site deployment and is relatively expensive.

[0004] Due to its advantages of large-scale, high-precision, and all-day and all-weather measurement, InSAR technology has gradually matured and been applied to urban ground subsidence observation, landslide deformation monitoring and other fields. However, vegetation growth or deformation gradient on the surface of the landslide body and construction on the surface of the reclamation area can cause the InSAR interferogram to lose coherence within a certain period of time. The existing permanent scatterer and identically distributed scatterer interferometry techniques mainly use pixel points that maintain high coherence throughout the entire observation period or maintain high coherence and continuity within a short observation period. For scatterers that only maintain coherence in part of the time period, effective monitoring results cannot be obtained. Summary of the invention

[0005] The main technical problem solved by the present invention is how to obtain the InSAR measurement result of the temporary coherent scatterer.

[0006] According to the first aspect, an embodiment provides an InSAR measurement method based on temporary coherent scatterers, comprising:

[0007] Acquire a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time;

[0008] Perform interference processing on every two SAR images in the SAR image sequence to obtain multiple interference graphs;

[0009] Determining temporary coherent pixel points and permanent coherent pixel points in the SAR image according to the plurality of interference patterns;

[0010] Performing phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points;

[0011] Converting the phase unwrapping sequence into a deformed sequence of temporary coherent pixels;

[0012] Based on the deformation sequence of the temporary coherent pixel points, the cumulative deformation sequence of the temporary coherent pixel points is fitted and acquired, and the cumulative deformation sequence of the temporary coherent pixel points is used to monitor the cumulative deformation amount of the temporary coherent pixel points.

[0013] According to the second aspect, an embodiment provides an InSAR measurement device based on temporary coherent scatterers, comprising:

[0014] An acquisition module, used for acquiring a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time;

[0015] An interference processing module is used to perform interference processing on every two SAR images in the SAR image sequence to obtain multiple interference graphs;

[0016] A temporary coherent pixel extraction module, used to determine temporary coherent pixel points and permanent coherent pixel points in the SAR image according to the plurality of interference patterns;

[0017] A phase unwrapping module, used to perform phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points;

[0018] A deformation sequence determination module, used for converting the phase unwrapping sequence into a deformation sequence of temporary coherent pixels;

[0019] The cumulative deformation sequence determination module is used to fit and obtain the cumulative deformation sequence of the temporary coherent pixel points based on the deformation sequence of the temporary coherent pixel points, and the cumulative deformation sequence of the temporary coherent pixel points is used to monitor the cumulative deformation amount of the temporary coherent pixel points.

[0020] According to a third aspect, an embodiment provides a computer-readable storage medium, including a program, wherein the program can be executed by a processor to implement the method described in the above embodiment.

[0021] According to the InSAR measurement method and device based on temporary coherent scatterers of the above-mentioned embodiment, by interfering every two SAR images in the SAR image sequence, a coherent phase matrix is ​​obtained, and temporary coherent pixels in the SAR image are identified based on the coherent phase matrix. Phase unwrapping is performed on the identified temporary coherent pixels to obtain a deformation sequence of the temporary coherent pixels, thereby realizing the monitoring of the deformation amount of the temporary coherent pixels to obtain the InSAR measurement results of the temporary coherent scatterers. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] Figure 1 A flowchart of an InSAR measurement method based on temporary coherent scatterers according to an embodiment;

[0023] Figure 2 A schematic structural diagram of an InSAR measurement device based on temporary coherent scatterers according to an embodiment;

[0024] Figure 3 This is the InSAR measurement result diagram using the existing method;

[0025] Figure 4 This is a diagram of InSAR measurement results using an embodiment of the present invention. DETAILED DESCRIPTION

[0026] The present invention is further described in detail below by specific embodiments in conjunction with the accompanying drawings. Wherein similar elements in different embodiments adopt associated similar element numbers. In the following embodiments, many detailed descriptions are for making the present application better understood. However, those skilled in the art can easily recognize that some features can be omitted in different situations, or can be replaced by other elements, materials, methods. In some cases, some operations related to the present application are not shown or described in the specification, this is to avoid the core part of the present application being overwhelmed by too much description, and for those skilled in the art, it is not necessary to describe these related operations in detail, and they can fully understand the related operations according to the description in the specification and the general technical knowledge in the art.

[0027] In addition, the features, operations or characteristics described in the specification can be combined in any appropriate manner to form various implementations. At the same time, the steps or actions in the method description can also be interchanged or adjusted in a manner that is obvious to those skilled in the art. Therefore, the various sequences in the specification and the drawings are only for the purpose of clearly describing a certain embodiment and are not meant to be a required sequence, unless otherwise specified that a certain sequence must be followed.

[0028] The serial numbers of the components in this document, such as "first", "second", etc., are only used to distinguish the objects described and do not have any order or technical meaning. The "connection" and "coupling" mentioned in this application, unless otherwise specified, include direct and indirect connections (couplings).

[0029] Please refer to Figure 1 , Figure 1 This is a flow chart of an InSAR measurement method based on temporary coherent scatterers according to an embodiment, which includes steps S10 to S50, which are described in detail below.

[0030] Step S10, obtaining a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time.

[0031] Step S20, performing interference processing on every two SAR images in the SAR image sequence to obtain a plurality of interference graphs.

[0032] Step S30: determining temporary coherent pixel points and permanent coherent pixel points in the SAR image according to the plurality of interference patterns.

[0033] Step S40, performing phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points.

[0034] Step S50: converting the phase unwrapping sequence into the deformation sequence of the temporary coherent pixel points, wherein the deformation sequence of the temporary coherent pixel points is used to represent the deformation of the temporary coherent pixel points in the SAR image.

[0035] Step S60, based on the deformation sequence of the temporary coherent pixel points, fitting and acquiring the cumulative deformation sequence of the temporary coherent pixel points, wherein the cumulative deformation sequence of the temporary coherent pixel points is used to monitor the cumulative deformation amount of the temporary coherent pixel points.

[0036] In this embodiment, it is assumed that the SAR image sequence of the area to be monitored includes N SAR images with continuous monitoring time, and an InSAR interference pattern sequence is obtained by interfering every two SAR images in the N SAR images, that is, N*(N-1) / 2 interference patterns are obtained. Moreover, all SAR images in the SAR image sequence are images taken in the same area, and have the same pixel points.

[0037] According to the N*(N-1) / 2 interference pattern obtained above, the temporary coherent pixels and permanent coherent pixels in the SAR image are determined, including:

[0038] Determine the interference phase and coherence coefficient of each coherent pixel in the interference pattern.

[0039] According to the interference phase and coherence coefficient of each coherent pixel point, a coherent phase matrix C corresponding to each coherent pixel point is formed, wherein the coherent phase matrix

[0040]

[0041] Among them, p and q are the serial numbers of the SAR image, and the values ​​of p and q are 1, 2...N. represents the coherence coefficient of SAR image p and SAR image q, φ pq is the interference phase of SAR image p and SAR image q.

[0042] In other words, all pixels in the SAR image sequence can obtain a coherent phase matrix after interference processing, that is, each pixel corresponds to a coherent phase matrix.

[0043] Perform the following traversal on all coherent pixels in the interference pattern:

[0044] If the coherent phase matrix C corresponding to the current coherent pixel point has an invalid interference phase during part of the monitoring time, the current coherent pixel point is a temporary coherent pixel point;

[0045] Among them, if there is a coherence coefficient less than the set threshold, the interference phase is an invalid interference phase, otherwise the interference phase is a valid interference phase. The invalid interference phase and the coherence coefficient are removed from the coherence phase matrix C, and then the effective monitoring period after removing the invalid interference phase is determined by connectivity judgment, wherein the connectivity judgment can adopt undirected connected graph judgment methods such as DFS, BFS, and Tarian. In order to increase the reliability of the effective monitoring time period judgment, this embodiment requires that the number of valid observations in the effective monitoring time period is greater than the number of SAR images in the time period, and the number of SAR images is ≥3, that is, the corresponding sub-coherence phase matrix Cs is obtained; otherwise, the interference phase in the monitoring period is an invalid interference phase.

[0046] Where s is the serial number of the sub-monitoring time period, and its value is [1, S]. Assuming that the m-th to n-th images are the acquired sub-monitoring time periods, the coherent phase matrix of the sub-monitoring time period is

[0047]

[0048] If the coherent phase matrix corresponding to the current coherent pixel has an effective interference phase in any monitoring time period, the current coherent pixel is a permanent coherent pixel.

[0049] In one embodiment, performing phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points includes:

[0050] Perform the following traversal on the temporary coherent pixels:

[0051] The phase sequence difference between the current temporary coherent pixel point and the pre-acquired reference pixel point is calculated to obtain the arc phase sequence.

[0052] Specifically, it includes: selecting permanent scatterer pixels or co-distributed coherent pixels around the current temporary coherent pixel as pre-acquired reference pixels. In this embodiment, the reference pixels are selected based on a dual threshold of spatial distance and quality of the reference pixels, and the quality of the reference pixels is mainly the amplitude deviation value of the permanent scatterer pixels or the coherence coefficient of the co-distributed target points.

[0053] The phase sequence difference between the current temporary coherent pixel and the pre-acquired reference pixel is calculated by the following formula:

[0054]

[0055] in, represents the phase sequence difference between the current temporary coherent pixel i and the reference pixel j at the kth monitoring time, represents the interference phase corresponding to the reference pixel j at the kth observation time, represents the interference phase corresponding to the current temporary coherent pixel point i at the k-th observation time; wherein the k-th observation time represents the k-th SAR image in the SAR image sequence.

[0056] The arc segment phase sequence is phase unwrapped in the time domain to obtain the unwrapped arc segment timing phase.

[0057] Specifically include: estimating the elevation error difference ΔE corresponding to the arc phase sequence j and the linear deformation rate difference Δv j .

[0058] This embodiment uses the temporal coherence coefficient as the quality standard for selecting valid arc segments. By maximizing the temporal coherence coefficient γ, the estimated value of the corresponding elevation error difference can be calculated. and an estimate of the difference in linear deformation rate If the coefficient γ exceeds the set threshold, the arc segment is considered valid. The temporal coherence coefficient can be expressed as:

[0059]

[0060] in, is the terrain error phase, is the phase value corresponding to the linear deformation rate difference, λ is the wavelength of the radar satellite electromagnetic wave, R is the distance from the radar satellite to the target point, θ is the incident angle of the electromagnetic wave, is the vertical baseline length of the interference pattern, Tsk is the time interval of the interference pattern.

[0061] The elevation error difference and linear deformation rate difference of the estimated arc phase sequence are added to the elevation error and deformation rate of the reference pixel point obtained in advance to obtain the elevation error of the current temporary coherent pixel point. and the linear deformation rate Among them, parameters such as the elevation error and deformation rate of the reference pixel point are calculated using existing methods.

[0062] If there are multiple valid arc segments, multiple elevation errors and linear deformation rates corresponding to the current temporary coherent pixel point can be obtained. This embodiment uses the median method to obtain the elevation error and linear deformation rate corresponding to the current temporary coherent pixel point.

[0063] Assuming that the number of valid arcs corresponding to the current temporary coherent pixel is H, the elevation error of the current temporary coherent pixel is and the linear deformation rate It can be expressed as:

[0064]

[0065]

[0066] Subtract the elevation error and linear deformation rate of the current temporary coherent pixel from the elevation error and deformation rate of the pre-acquired reference pixel to obtain the arc segment elevation error. and arc segment linear deformation rate

[0067] The phases corresponding to the arc segment elevation error and the arc segment linear deformation rate are removed from the arc segment phase sequence to obtain the residual arc segment timing phase.

[0068] In this embodiment, the phases corresponding to the arc segment elevation error and the arc segment linear deformation rate are removed from the arc segment phase sequence, and the residual arc segment timing phase is obtained. It is mainly phase noise, atmospheric phase residual term, etc., which can be considered to be less than π. Finally, the unwrapped arc segment timing phase can be expressed as:

[0069]

[0070] in, Indicates the arc segment elevation error The corresponding phase, Indicates the linear deformation rate of the arc segment The corresponding phase.

[0071] In one embodiment, performing phase unwrapping in the spatial domain on the unwrapped arc segment timing phase to obtain a phase unwrapping sequence of the current temporary coherent pixel point further includes:

[0072] If there are multiple reference pixels acquired in advance, multiple phase unwrapping sequences corresponding to the current temporary coherent pixel can be obtained; based on the multiple reference pixels acquired in advance, this embodiment can obtain phase unwrapping sequences corresponding to H temporary coherent pixels.

[0073] If all H phase unwrapping sequences are equal, the current temporary coherent pixel point is a credible monitoring point;

[0074] If the H phase unwrapping sequences are partially equal, and the number of equal phase unwrapping sequences is greater than the number of unequal phase unwrapping sequences, then the current temporary coherent pixel point is a credible monitoring point;

[0075] Otherwise, the current temporary coherent pixel is an unreliable monitoring point.

[0076] In this embodiment, the unwrapped phase sequence of the credible monitoring point can be converted into the deformation sequence d according to the wavelength of the incident wave of the radar satellite.

[0077] Since the deformation sequence corresponding to the temporary coherent pixels is discontinuous, the deformation cannot be transferred between the temporary coherent time periods, resulting in the inability to obtain the accumulated deformation.

[0078] This embodiment first calculates the deformation sequence difference Δd (arc deformation sequence) between the temporary coherent pixel point and the nearest reference pixel point in the surrounding area, and then performs linear or exponential function fitting on the arc deformation sequence in the first effective monitoring time period to calculate the function polynomial coefficient. For landslide monitoring, the present invention adopts a linear function fitting method. For reclamation settlement monitoring, according to the consolidation characteristics of soft soil, the present invention adopts an exponential function fitting method. According to the above fitting polynomial, the deformation at the start time of the second effective monitoring time period is predicted. Then, fitting and prediction are performed in sequence. Finally, the cumulative deformation sequence of the arc segment is obtained. The cumulative deformation sequence of the arc segment is added to the deformation sequence of the nearest reference point in the surrounding area to obtain the cumulative deformation sequence corresponding to the temporary coherent pixel point. According to its cumulative deformation sequence, the cumulative variable of the temporary coherent pixel point can be observed.

[0079] Please refer to Figure 2 , Figure 2 The present invention is a structural schematic diagram of an InSAR measurement device based on temporary coherent scatterers according to an embodiment, comprising: an acquisition module 10, an interference processing module 20, a temporary coherent pixel point extraction module 30, a phase unwrapping module 40, a deformation sequence determination module 50 and a cumulative deformation sequence determination module 60.

[0080] The acquisition module 10 is used to acquire a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time.

[0081] The interference processing module 20 is used to perform interference processing on every two SAR images in the SAR image sequence to obtain a plurality of interference graphs.

[0082] The temporary coherent pixel point extraction module 30 is used to determine the temporary coherent pixel points and the permanent coherent pixel points in the SAR image according to the plurality of interference patterns.

[0083] The phase unwrapping module 40 is used to perform phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points.

[0084] The deformation sequence determination module 50 is used to convert the phase unwrapping sequence into the deformation sequence of the temporary coherent pixel points, and the deformation sequence of the temporary coherent pixel points is used to express the deformation of the temporary coherent pixel points in the SAR image.

[0085] The accumulated deformation sequence determining module 60 is used to fit and obtain the accumulated deformation sequence of the temporary coherent pixel points based on the deformation sequence of the temporary coherent pixel points, wherein the accumulated deformation sequence of the temporary coherent pixel points is used to monitor the accumulated deformation amount of the temporary coherent pixel points.

[0086] It should be noted that the device module in this embodiment corresponds to the method steps in the above embodiment, and its specific implementation has been specifically described in the above embodiment, which will not be repeated here.

[0087] This example uses a landslide as an example for experiment. After data analysis and investigation, the slope continues to deform after the reservoir is filled with water, and tends to increase over time. Sentinel-1 image data is used as test data, with an image spatial resolution of 20m in azimuth and 5m in distance, the operating band is C band, the image shooting time range is from March 8, 2018 to July 1, 2019, and the image acquisition frequency is 1 scene every 12 days.

[0088] Firstly, the traditional time-series InSAR analysis method is used to obtain the deformation monitoring results of the slope surface. Figure 3 As shown in the figure, the types of monitoring points extracted include permanent scatterers and identically distributed target scatterers. Analysis of the monitoring results shows that based on Sentinel-1 satellite images, a certain number of deformation monitoring points can be obtained on the surface of the Guobu slope, and the deformation rate of the surrounding monitoring points is relatively large, with the maximum deformation rate of -282.2 mm / yr, located at the top of the slope. Analysis found that the interference pattern in the middle part of the Guobu slope was decoherent due to the excessive deformation rate, and no effective deformation monitoring points could be obtained.

[0089] This embodiment uses temporary coherent scatterer interferometry technology to process the Sentinel-1 satellite image to obtain the deformation monitoring results of the landslide body. Figure 4 As shown, Figure 4 is the deformation rate of the monitoring point obtained by the InSAR measurement method using temporary coherent pixels. Figure 3 Compared with the traditional data processing results in the above data, the embodiment of the present invention can obtain more monitoring points in the middle part of the slope where the deformation rate is relatively large, and the number of monitoring points on the surface of the slope can be increased by 323%. The deformation rate of the monitoring points on the slope surface is relatively large, with a maximum deformation rate of -372.2 mm / yr, located in the middle part of the slope. Analysis shows that for areas with excessive deformation rate and decoherence, the embodiment of the present invention can obtain more effective monitoring points, which is helpful for the identification of potential danger points of rapidly deforming slopes.

[0090] Those skilled in the art will appreciate that all or part of the functions of the various methods in the above-mentioned embodiments can be implemented by hardware or by computer programs. When all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can be stored in a computer-readable storage medium, and the storage medium can include: read-only memory, random access memory, disk, optical disk, hard disk, etc., and the program is executed by a computer to implement the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor, all or part of the above-mentioned functions can be implemented. In addition, when all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can also be stored in a storage medium such as a server, another computer, disk, optical disk, flash disk or mobile hard disk, and can be downloaded or copied and saved in the memory of the local device, or the system of the local device is updated, and when the program in the memory is executed by the processor, all or part of the functions in the above-mentioned embodiments can be implemented.

[0091] The above specific examples are used to illustrate the present invention, which is only used to help understand the present invention and is not intended to limit the present invention. For those skilled in the art, according to the concept of the present invention, some simple deductions, modifications or substitutions can be made.

Claims

1. An InSAR measurement method based on temporary coherent scatterers, characterized in that: include: Acquire a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time; Perform interference processing on every two SAR images in the SAR image sequence to obtain multiple interference graphs; According to the plurality of interference patterns, temporary coherent pixels and permanent coherent pixels in the SAR image are determined, wherein the interference phase and the coherence coefficient of each coherent pixel in the interference pattern are determined; according to the interference phase and the coherence coefficient of each coherent pixel, a coherent phase matrix corresponding to each coherent pixel is constructed; and the following traversal is performed on all coherent pixels in the interference pattern: if the interference phase of a part of the monitoring time period in the coherent phase matrix corresponding to the current coherent pixel is invalid, the current coherent pixel is a temporary coherent pixel; otherwise, the current coherent pixel is a permanent coherent pixel; Performing phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points; Converting the phase unwrapping sequence into a deformed sequence of temporary coherent pixels; Based on the deformation sequence of the temporary coherent pixel points, the cumulative deformation sequence of the temporary coherent pixel points is fitted and acquired, and the cumulative deformation sequence of the temporary coherent pixel points is used to monitor the cumulative deformation amount of the temporary coherent pixel points.

2. The InSAR measurement method based on temporary coherent scatterers according to claim 1, characterized in that: After performing phase unwrapping on the determined temporary coherent pixel points and before converting the phase unwrapping sequence into the deformed sequence of the temporary coherent pixel points, the method further includes: Based on the phase unwrapping sequence of the temporary coherent pixel point, determine whether the temporary coherent pixel point is a credible monitoring point, and output a determination result; If the judgment result is a credible monitoring point, the phase unwrapping sequence of the temporary coherent pixel points is converted into a deformation sequence of the temporary coherent pixel points.

3. The InSAR measurement method based on temporary coherent scatterers according to claim 1, characterized in that: If the interference phase of some monitoring time periods in the coherent phase matrix corresponding to the current coherent pixel is invalid, it includes: The interferometric phase invalidity of some monitoring time periods is determined by: If the coherence coefficient is less than the set threshold, the interference phase corresponding to the monitoring time period is an invalid interference phase.

4. The InSAR measurement method based on temporary coherent scatterers according to claim 1, characterized in that: Performing phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points includes: The temporary coherent pixel points are traversed as follows: Calculate the phase sequence difference between the current temporary coherent pixel point and the pre-acquired reference pixel point, that is, obtain the arc phase sequence; Performing phase unwrapping in the time domain on the arc segment phase sequence to obtain an unwrapped arc segment timing phase; The unwrapped arc segment temporal phase is phase unwrapped in the spatial domain to obtain the phase unwrapped sequence of the current temporary coherent pixel point.

5. The InSAR measurement method based on temporary coherent scatterers according to claim 4, characterized in that: The calculating of the phase sequence difference between the current temporary coherent pixel point and the pre-acquired reference pixel point includes: Selecting permanent scatterer pixels or co-distributed coherent pixels around the current temporary coherent pixel as pre-acquired reference pixels; The phase sequence difference between the current temporary coherent pixel and the pre-acquired reference pixel is calculated by the following formula: in, represents the phase sequence difference between the current temporary coherent pixel i and the reference pixel j at the kth observation time, represents the interference phase corresponding to the reference pixel j at the kth observation time, represents the interference phase corresponding to the current temporary coherent pixel point i at the k-th observation time; wherein the k-th observation time represents the k-th SAR image in the SAR image sequence.

6. The InSAR measurement method based on temporary coherent scatterers according to claim 4, characterized in that: Performing phase unwrapping in the time domain on the arc segment phase sequence to obtain an unwrapped arc segment timing phase, including: estimating the elevation error difference and the linear deformation rate difference corresponding to the arc segment phase sequence; Adding the estimated elevation error difference and linear deformation rate difference corresponding to the arc phase sequence to the elevation error and deformation rate of the reference pixel point acquired in advance to obtain the elevation error and linear deformation rate of the current temporary coherent pixel point; If there are multiple valid arc segments, multiple elevation errors and linear deformation rates corresponding to the current temporary coherent pixel points are obtained; Subtract the elevation error and linear deformation rate of the current temporary coherent pixel from the elevation error and deformation rate of the reference pixel obtained in advance to obtain the arc segment elevation error and arc segment linear deformation rate; The phases corresponding to the arc segment elevation error and the arc segment linear deformation rate are removed from the arc segment phase sequence to obtain a residual arc segment timing phase; Adding the phase corresponding to the arc segment elevation error and the arc segment linear deformation rate to the residual arc segment timing phase to obtain the unwrapped arc segment timing phase; The unwrapped arc segment temporal phase is added to the phase unwrapping sequence of the reference pixel point acquired in advance to obtain the phase unwrapping sequence of the temporary coherent pixel point.

7. The InSAR measurement method based on temporary coherent scatterers according to claim 4, characterized in that: The phase unwrapping of the unwrapped arc segment time sequence phase in the spatial domain to obtain the phase unwrapping sequence of the current temporary coherent pixel point also includes: If the number of the pre-acquired reference pixel points is multiple, multiple phase unwrapping sequences corresponding to the current temporary coherent pixel point can be obtained; If the multiple phase unwrapping sequences are all equal, the current temporary coherent pixel point is a credible monitoring point; If the plurality of phase unwrapping sequences are partially equal, and the number of equal phase unwrapping sequences is greater than the number of unequal phase unwrapping sequences, then the current temporary coherent pixel point is a credible monitoring point; Otherwise, the current temporary coherent pixel point is an unreliable monitoring point.

8. An InSAR measurement device based on temporary coherent scatterers, characterized in that: include: An acquisition module, used for acquiring a SAR image sequence of the area to be monitored, wherein the SAR image sequence includes multiple SAR images of continuous monitoring time; An interference processing module is used to perform interference processing on every two SAR images in the SAR image sequence to obtain multiple interference graphs; A temporary coherent pixel extraction module is used to determine temporary coherent pixels and permanent coherent pixels in the SAR image according to the plurality of interference patterns, wherein the interference phase and the coherence coefficient of each coherent pixel in the interference pattern are determined; a coherent phase matrix corresponding to each coherent pixel is constructed according to the interference phase and the coherence coefficient of each coherent pixel; and the following traversal is performed on all coherent pixels in the interference pattern: if the interference phase of a part of the monitoring time period in the coherent phase matrix corresponding to the current coherent pixel is invalid, the current coherent pixel is a temporary coherent pixel; otherwise, the current coherent pixel is a permanent coherent pixel; A phase unwrapping module, used to perform phase unwrapping on the determined temporary coherent pixel points to obtain a phase unwrapping sequence of the temporary coherent pixel points; A deformation sequence determination module, used for converting the phase unwrapping sequence into a deformation sequence of temporary coherent pixels; The cumulative deformation sequence determination module is used to fit and obtain the cumulative deformation sequence of the temporary coherent pixel points based on the deformation sequence of the temporary coherent pixel points, and the cumulative deformation sequence of the temporary coherent pixel points is used to monitor the cumulative deformation amount of the temporary coherent pixel points.

9. A computer-readable storage medium, characterized in that: The method comprises a program which can be executed by a processor to implement the method according to any one of claims 1 to 7.

Images