A wide-area deformation area InSAR monitoring precision evaluation method

Abstract

The application discloses a kind of wide-area deformation area InSAR monitoring precision evaluation method, it is applicable to deformation monitoring field.In wide-area deformation area, level route SZ is laid out, n level observation points are made in situ along SZ, and its plane position coordinate is measured;Surface deformation information is extracted from m SAR images covering wide-area deformation area;The level measurement method with designed measurement precision δ s Is used to observe SZ m times, obtain the elevation of each level observation point, calculate the elevation deformation value of level observation point between two level measurement interval periods;The deformation monitoring value of the corresponding pixel of level observation point is extracted;And the cumulative deformation value of each level observation point and the cumulative deformation value monitored by InSAR are calculated, the deformation difference value between adjacent level points is obtained, to calculate the evaluation value of InSAR monitoring precision.It is simple in steps, and the deformation difference value of adjacent deformation points is used for precision evaluation, and the effect of InSAR precision evaluation for wide-area deformation area is good.

Description

Technical Field

[0001] This invention relates to a method for assessing monitoring accuracy, and is particularly applicable to a method for assessing the accuracy of InSAR monitoring in a wide-area deformation region, belonging to the field of InSAR monitoring accuracy assessment in wide-area deformation regions. Background Technology

[0002] InSAR deformation monitoring belongs to the field of image geodesy. It refers to the extraction of movement and deformation information of the Earth's surface using multi-temporal SAR images covering the same geographical area. This technology has been widely used in fields such as coal mining subsidence monitoring, dam deformation monitoring, earthquake disaster monitoring, and landslide disaster monitoring.

[0003] Unlike traditional methods of monitoring Earth's surface deformation, such as leveling, total station surveying, and Global Navigation Satellite Systems (GNSS) surveying, which can only observe deformation information at discrete points, InSAR technology can acquire planar deformation information of the Earth's surface in a non-contact sensing manner. Furthermore, synthetic aperture radar (SAR) uses side-looking imaging, enabling it to penetrate clouds, rain, fog, and haze, and acquire ground information in all weather conditions. Methods for obtaining Earth's surface deformation information from SAR imagery have undergone long-term development, with early methods primarily using D-InSAR, which could extract ground deformation information using only two SAR images. However, conventional D-InSAR technology is affected by adverse factors such as spatiotemporal decorrelation and atmospheric effects. With the launch of more SAR satellites and the accumulation of historical data, to address the shortcomings of conventional D-InSAR technology, scholars both domestically and internationally have proposed a series of time-series InSAR methods. These methods process long-term SAR images covering the same area to extract deformation information from highly coherent point targets with strong backscattering characteristics. Representative temporal InSAR methods include Persistent Scatterer InSAR (PS-InSAR), Small Baseline Subset InSAR (SBAS-InSAR), and Distributed Scatterer InSAR (DS-InSAR). In addition, some researchers have studied the use of pixel-tracking techniques to extract the Earth's surface deformation field under large deformation gradients from SAR images.

[0004] Although InSAR technology, which extracts surface deformation from SAR imagery, has been extensively studied and applied in many fields, the accuracy of surface deformation results obtained using InSAR has remained a challenge for researchers and engineers engaged in its research and application. Existing research typically uses deformation values ​​from a subset of observation points monitored using leveling or GNSS techniques as reference values ​​to estimate the accuracy of InSAR monitoring results. However, this approach suffers from two problems: first, it ignores the inherent errors in surface deformation values ​​monitored by leveling or GNSS techniques; second, when surface deformation occurs over a wide area, and the control points used for leveling or GNSS monitoring are located within the deformation zone, the surface deformation results obtained using these techniques are inherently unreliable due to the influence of the widespread deformation. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this paper presents a theoretically rigorous, highly accurate, and effective method for evaluating the accuracy of InSAR monitoring in wide-area deformation regions.

[0006] To achieve the above objectives, the present invention provides a method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region, comprising the following steps: a) designing a leveling route SZ in the wide-area deformation region, wherein the leveling route SZ contains n leveling points with equal spacing between each leveling point; and constructing n leveling observation points S according to the leveling point locations designed in the leveling route SZ. i (i = 1, 2, ..., n);

[0007] b. Measure and record the S values ​​at each leveling observation point. i The planar position coordinates (X) of (i = 1, 2, ..., n) i ,Y i ), i = 1, 2, ..., n, where X i It is the i-th leveling point S i The ordinate in a Cartesian coordinate system, Y i It is the i-th leveling point S i The x-coordinate in a Cartesian coordinate system;

[0008] c. Acquire m SAR images covering a wide deformation area. t (t=1,2,...,m), and register, geometrically correct and geocode m SAR images, where the number of SAR images m≥2;

[0009] d. Using InSAR technology, from m SAR images I t Image CI, which records the surface deformation value, was extracted from (t=1,2,...,m).t (t=1,2,...,m-1), where CI t Indicates from SAR image I t and SAR image I t+1 Deformation image of wide-area deformation region A extracted from it;

[0010] e employs a design measurement accuracy of δ s The leveling measurement method for each leveling observation point S i (i = 1, 2, ..., n) Perform synchronous leveling to acquire leveling data. The measurement time is the same as the acquisition time of SAR image I. t The time intervals (t = 1, 2, ..., p) are kept synchronized; a total of m synchronized leveling measurements are performed. Based on the elevations of the leveling control points and the leveling measurement data, the elevations h(t, i) (t = 1, 2, ..., m; i = 1, 2, ..., n) of each leveling observation point are obtained. The elevation deformation value CS of the leveling observation point during the interval between two leveling measurements is calculated. t,i (t = 1, 2, ..., m-1; i = 1, 2, ..., n), where h(t, i) represents the i-th leveling point S. i In acquiring the t-th image I t The elevation values ​​obtained from leveling measurements on the same day, CS t,i =h(t+1,i)-h(t,i) represents the i-th leveling point S i In acquiring SAR image I t and SAR image I t+1 The reference values ​​for surface deformation within the time interval, wherein the leveling control points and their number required for leveling meet the leveling measurement specifications;

[0011] Using formula Calculate the measurement accuracy δ s ,

[0012] Where δ 站 This represents the theoretical mean square error of the elevation difference observation at each station. Based on leveling measurement specifications, δ is designed and determined by comprehensively considering the accuracy of the leveling instruments used, the measuring rod, and the specific observation requirements. 站 q represents the number of stations between two adjacent leveling points;

[0013] f is based on the measured leveling observation points S i Planar position coordinates (X i ,Y i ), in the deformation image CI obtained using InSAR technology t Extract the deformation value CI of the image pixels corresponding to the leveling observation points from (t=1,2,...,m-1). t,i(t=1,2,...,m-1; i=1,2,...,n), where CI t,i S represents the i-th leveling point. i In acquiring SAR image I t and SAR image I t+1 Surface deformation monitoring values ​​at time intervals;

[0014] g using formula Calculate the S of each leveling point i (i = 1, 2, ..., n) are the cumulative deformation reference values ​​CSL obtained by leveling. i (i = 1, 2, ..., n);

[0015] Using formula Calculate S at each leveling observation point i (i = 1, 2, ..., n) Cumulative deformation monitoring values ​​(CIL) obtained using InSAR monitoring technology i (i = 1, 2, ..., n), where ∑ is the summation symbol commonly used in mathematics;

[0016] h uses the formula CSLN i,i+1 =CSL i -CSL i+1 (i = 1, 3, ..., n-1) Calculate the deformation difference reference value between adjacent leveling points obtained by leveling method; use the formula CILN i,i+1 =CIL i -CIL i+1 (i = 1, 3, ..., n-1) Calculate the deformation difference monitoring values ​​between adjacent leveling observation points obtained using InSAR monitoring technology; using the deformation difference reference value and deformation difference monitoring value between adjacent leveling observation points, the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in the wide-area deformation region can be calculated. L And the evaluation value δ1 for the accuracy of InSAR cumulative deformation monitoring in a wide-area deformation region.

[0017] Furthermore, the characteristics of the wide-area deformation region in step a are as follows: the ground surface undergoes large-scale deformation. When the leveling method is used for deformation monitoring, the leveling control points are located in the wide-area deformation region. The leveling control points are affected by the large-scale deformation of the ground surface and thus deform, which leads to inaccurate leveling point deformation values ​​calculated based on the leveling control points.

[0018] Furthermore, in step a, the number n of leveling points laid out on the leveling route SZ is an even number. In flat terrain areas, the distance between two adjacent leveling points is basically equal, and in other areas, the number of leveling stations between two adjacent observation points is the same.

[0019] Furthermore, in step d, InSAR technology refers to a method for extracting surface deformation from SAR images, including any one of the following methods: SBAS-InSAR, PS-InSAR, DS-InSAR, D-InSAR, and Pixel-Tracking.

[0020] Furthermore, in the leveling measurement described in step e:

[0021] If, when conducting leveling surveys in flat areas, the theoretical value of the mean square error δ per kilometer of leveling measurement is determined in the design... k Then, based on the number of stations z that need to be measured per kilometer, first apply the formula... Calculate δ 站 Then according to the formula To calculate the measurement accuracy, the distance between the two leveling observation points is 1 kilometer, and the formula is used. Calculate measurement accuracy;

[0022] If the leveling plan involves forward and backward measurements, then the formula is used. The measured value δ is calculated from the measured data to determine the measurement accuracy. s ′, where h i往 h is the observed value of the cross-sectional elevation difference between the i-th leveling point and the (i+1)-th leveling point. i返 It is the observed value of the back-measured elevation difference between the i-th leveling point and the (i+1)-th leveling point.

[0023] Furthermore, in step h, the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in the wide-area deformation region is calculated. L The formula for evaluating the accuracy of a single InSAR deformation monitoring in a wide-area deformation region, δ1, is as follows:

[0024]

[0025]

[0026] In the formula: [·] is the Gaussian summation symbol. It is the symbol for taking the square root in mathematics.

[0027] Furthermore, when it is necessary to further reduce the impact of leveling measurement errors on accuracy estimation to meet higher accuracy requirements, the following formula is used to calculate the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in a wide-area deformation region. L The evaluation value δ1 for the accuracy of InSAR deformation monitoring in a wide-area deformation region:

[0028]

[0029]

[0030] When the leveling survey scheme involves forward and backward measurements, the measured value δ of the measurement accuracy can be used in the formula. s ′ replace δ s Perform accuracy estimation.

[0031] Beneficial Effects: This invention addresses the problem of accuracy estimation for surface deformation monitoring results based on InSAR technology in wide-area deformation regions. When leveling control points are located within the deformation region and are subject to deformation, it effectively estimates the accuracy of InSAR deformation monitoring in wide-area deformation regions by constructing deformation difference values ​​unaffected by the deformation of the leveling control points. Furthermore, considering the measurement errors inherent in leveling results, it further obtains an estimate of the InSAR monitoring accuracy in wide-area deformation regions unaffected by leveling measurement errors. Therefore, this invention provides a new approach for estimating the accuracy of deformation monitoring in wide-area deformation regions using InSAR technology. Attached Figure Description

[0032] Figure 1 This is a flowchart of a method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to the present invention.

[0033] Figure 2 This is a schematic diagram of the deformation area and leveling point layout for the InSAR monitoring accuracy assessment method for wide-area deformation areas according to the present invention. Detailed Implementation

[0034] The present invention will now be described in further detail with reference to the accompanying drawings and embodiments;

[0035] like Figure 1 As shown, the present invention provides a method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region, comprising the following steps:

[0036] Step a: In Figure 2 As shown, a wide-area deformation region A is arranged with a line as shown in the figure. Figure 2 The leveling route SZ shown contains n=18 leveling points. The actual leveling process along SZ is as follows: Figure 2 The blue dots indicate the n leveling observation points S. i (i = 1, 2, ..., n); where the wide-area deformation region A refers to a large-scale deformation of the earth's surface. When using leveling methods for deformation monitoring, such as... Figure 2 The leveling control points indicated by the red dots are located in a wide-area deformation zone. These leveling control points are affected by large-scale surface deformation, resulting in inaccurate leveling deformation values ​​calculated based on them. The number of leveling points n=18 along the leveling route SZ satisfies the requirement of being an even number. The distance between adjacent leveling points is approximately equal in flat terrain and requires the same number of leveling stations between adjacent observation points in other areas.

[0037] Step b: Measure and record the levels S. i The planar position coordinates (X) of (i = 1, 2, ..., n) i ,Y i ), i = 1, 2, ..., n, where X i It is the i-th leveling point S i The ordinate in a Cartesian coordinate system, Y i It is the i-th leveling point S i The x-coordinate in a Cartesian coordinate system;

[0038] Step c: Acquire m SAR images I covering the wide-area deformation region A, after registration, geometric correction, and geocoding. t (t=1,2,...,m), where m=16, satisfying the requirement that the number of SAR images m≥2;

[0039] Step d: Use InSAR technology to extract m SAR images I t Image CI, which records the surface deformation value, was extracted from (t=1,2,...,m). t (t=1,2,...,m-1), where CI t Indicates from SAR image I t and SAR image I t+1 The deformation image of region A extracted from SAR images; where InSAR technology refers to the method of extracting surface deformation from SAR images, and there is no particular limitation on which method among SBAS-InSAR, PS-InSAR, DS-InSAR, D-InSAR, Pixel-Tracking, etc. is used. In this embodiment, PS-InSAR is used to extract surface deformation from SAR images.

[0040] Step e: Acquiring SAR image I t On the same day (t=1,2,...,m), a measurement with a designed accuracy of δ was used. s The leveling method is used for each leveling point S on SZ. i Synchronous leveling is conducted for (i = 1, 2, ..., n), requiring a total of m leveling measurements. Based on the elevations of the leveling control points and the leveling measurement data, the elevations h(t, i) (t = 1, 2, ..., m; i = 1, 2, ..., n) of each leveling observation point are obtained. The elevation deformation CS of the leveling points during the interval between two leveling measurements is calculated. t,i (t = 1, 2, ..., m-1; i = 1, 2, ..., n), where h(t, i) represents the i-th leveling point S. i In acquiring the t-th image I t The elevation values ​​obtained from leveling measurements on the same day, CSt,i =h(t+1,i)-h(t,i) represents the i-th leveling point S i In acquiring SAR image I t and SAR image I t+1 The time interval is the reference value of surface deformation; among which, the leveling control points and their number required for leveling meet the leveling measurement specifications, and the measurement accuracy δ s Formulas can be used Calculate, where δ 站 This refers to the theoretical value of the mean square error of the elevation difference observation at each station. It can be determined by referring to the leveling measurement specifications and comprehensively considering the accuracy of the leveling instruments used, the measuring rod, and the specific observation requirements. 站 q is the number of stations between two adjacent leveling points; if the theoretical error δ per kilometer of leveling measurement is determined when conducting leveling measurements in a flat area... k Then, based on the number of stations z that need to be measured per kilometer, first apply the formula... Calculate δ 站 Then according to the formula Calculate the measurement accuracy; however, it is recommended to use the formula based on the distance l kilometers between the two benchmarks. Calculate the measurement accuracy; if the leveling plan includes forward and backward measurements, the formula can be used. The measured value δ is calculated from the measured data to determine the measurement accuracy. s ′, where h i往 h is the observed value of the cross-sectional elevation difference between the i-th leveling point and the (i+1)-th leveling point. i返 It is the observed value of the back-measured elevation difference between the i-th leveling point and the (i+1)-th leveling point;

[0041] Step f: Based on the measurements of each leveling point S i Planar position coordinates (X i ,Y i ), in the deformation image CI obtained using InSAR technology t Extract the deformation value CI of the pixel corresponding to the leveling point position from (t=1,2,...,m-1). t,i (t=1,2,...,m-1; i=1,2,...,n), where CI t,i S represents the i-th leveling point. i In acquiring SAR image I t and SAR image I t+1 Surface deformation monitoring values ​​at time intervals;

[0042] Step g: Using the formula Calculate the S of each leveling point i (i = 1, 2, ..., n) are the cumulative deformation reference values ​​CSL obtained by leveling.i (i = 1, 2, ..., n); using the formula Calculate the S of each leveling point i (i = 1, 2, ..., n) Cumulative deformation monitoring values ​​(CIL) obtained using InSAR monitoring technology i (i = 1, 2, ..., n), where ∑ is the summation symbol commonly used in mathematics;

[0043] Step h: Using the formula CSLN i,i+1 =CSL i -CSL i+1 (i = 1, 3, ..., n-1) Calculate the deformation difference reference value between adjacent leveling points obtained by leveling method; use the formula CILN i,i+1 =CIL i -CIL i+1 (i = 1, 3, ..., n-1) Calculate the differential deformation monitoring values ​​between adjacent leveling points obtained using InSAR monitoring technology; based on this, calculate the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in the wide-area deformation region. L The evaluation value δ1 for the accuracy of single-time monitoring of InSAR cumulative deformation in a wide-area deformation region; the evaluation value δ used to calculate the accuracy of InSAR cumulative deformation monitoring in a wide-area deformation region. L The formula for evaluating the accuracy of a single InSAR deformation monitoring in a wide-area deformation region, δ1, is as follows:

[0044]

[0045]

[0046] In the formula: [·] is the Gaussian summation symbol. The square root symbol is used in mathematics. When it is necessary to further reduce the impact of leveling errors on accuracy estimation, the following formula can be used to calculate the assessment value δ of the InSAR cumulative deformation monitoring accuracy in a wide-area deformation region. L The evaluation value δ1 for the accuracy of InSAR deformation monitoring in a wide-area deformation region:

[0047]

[0048]

[0049] When the leveling survey scheme involves forward and backward measurements, the measured value δ of the measurement accuracy can be used in the formula. s ′ replace δ s Perform accuracy estimation.

[0050] It should be understood that any parts not described in detail in this specification belong to the prior art.

[0051] It should be understood that the above description of the implementation examples is quite detailed, but it should not be considered as a limitation on the scope of protection of this invention. Those skilled in the art, under the guidance of this invention, can make substitutions or modifications without departing from the scope of protection of the claims of this invention, and all such substitutions or modifications fall within the scope of protection of this invention. The scope of protection of this invention should be determined by the appended claims.

Claims

1. A method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region, characterized in that... Includes the following steps: In a wide-area deformation region, a leveling route SZ is designed, containing n leveling points with equal spacing between them. Based on the locations of the leveling points designed in the leveling route SZ, n leveling observation points S are constructed. i (i = 1, 2, ..., n); b. Measure and record the S values ​​at each leveling observation point. i The planar position coordinates (X) of (i = 1, 2, ..., n) i ,Y i ), i = 1, 2, ..., n, where X i It is the i-th leveling point S i The ordinate in a Cartesian coordinate system, Y i It is the i-th leveling point S i The x-coordinate in a Cartesian coordinate system; c. Acquire m SAR images covering a wide deformation area. t (t=1,2,...,m), and register, geometrically correct and geocode m SAR images, where the number of SAR images m≥2; d. Using InSAR technology, from m SAR images I t Image CI, which records surface deformation values, was extracted from (t=1,2,...,m). t (t=1,2,...,m-1), where CI t Indicates from SAR image I t and SAR image I t+1 Deformation image of wide-area deformation region A extracted from it; e employs a design measurement accuracy of δ s The leveling measurement method for each leveling observation point S i (i = 1, 2, ..., n) Perform synchronous leveling to acquire leveling data. The measurement time is the same as the acquisition time of SAR image I. t The time intervals (t = 1, 2, ..., p) are kept synchronized; a total of m synchronized leveling measurements are performed. Based on the elevations of the leveling control points and the leveling measurement data, the elevations h(t, i) (t = 1, 2, ..., m; i = 1, 2, ..., n) of each leveling observation point are obtained. The elevation deformation value CS of the leveling observation point during the interval between two leveling measurements is calculated. t,i (t = 1, 2, ..., m-1; i = 1, 2, ..., n), where h(t, i) represents the i-th leveling point S. i In acquiring the t-th image I t The elevation values ​​obtained from leveling measurements on the same day, CS t,i =h(t+1,i)-h(t,i) represents the i-th leveling point S i In acquiring SAR image I t and SAR image I t+1 The reference values ​​for surface deformation within the time interval, wherein the leveling control points and their number required for leveling meet the leveling measurement specifications; Using formula Calculate the measurement accuracy δ s , Where δ 站 This represents the theoretical mean square error of the elevation difference observation at each station. Based on leveling measurement specifications, δ is designed and determined by comprehensively considering the accuracy of the leveling instruments used, the measuring rod, and the specific observation requirements. 站 q represents the number of stations between two adjacent leveling points; f is based on the measured leveling observation points S i Planar position coordinates (X i ,Y i ), in the deformation image CI obtained using InSAR technology t Extract the deformation value CI of the image pixels corresponding to the leveling observation points from (t=1,2,...,m-1). t,i (t=1,2,...,m-1; i=1,2,...,n), where CI t,i S represents the i-th leveling point. i In acquiring SAR image I t and SAR image I t+1 Surface deformation monitoring values ​​at time intervals; g using formula Calculate the S of each leveling point i (i = 1, 2, ..., n) are the cumulative deformation reference values ​​CSL obtained by leveling. i (i = 1, 2, ..., n); Using formula Calculate S at each leveling observation point i (i = 1, 2, ..., n) Cumulative deformation monitoring values ​​(CIL) obtained using InSAR monitoring technology i (i = 1, 2, ..., n), where ∑ is the summation symbol commonly used in mathematics; h uses the formula CSLN i,i+1 =CSL i -CSL i+1 (i = 1, 3, ..., n-1) Calculate the deformation difference reference value between adjacent leveling points obtained by leveling method; use the formula CILN i,i+1 =CIL i -CIL i+1 (i = 1, 3, ..., n-1) Calculate the deformation difference monitoring values ​​between adjacent leveling observation points obtained using InSAR monitoring technology; using the deformation difference reference value and deformation difference monitoring value between adjacent leveling observation points, the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in the wide-area deformation region can be calculated. L And the evaluation value δ1 for the accuracy of InSAR cumulative deformation monitoring in a wide-area deformation region.

2. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 1, characterized in that: The characteristics of the wide-area deformation region in step a are as follows: the surface undergoes large-scale deformation. When the leveling method is used for deformation monitoring, the leveling control points are located in the wide-area deformation region. The leveling control points are affected by the large-scale deformation of the surface and thus deform, which leads to inaccurate leveling point deformation values ​​calculated based on the leveling control points.

3. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 1, characterized in that: In step a, the number of leveling points n laid out on the leveling route SZ is an even number. In flat terrain, the distance between two adjacent leveling points is basically equal, while in other areas, the number of leveling stations between two adjacent observation points is the same.

4. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 1, characterized in that: In step d, InSAR technology refers to the method of extracting surface deformation from SAR images, including any one of the following methods: SBAS-InSAR, PS-InSAR, DS-InSAR, D-InSAR, and Pixel-Tracking.

5. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 1, characterized in that: In the leveling measurement described in step e: If, when conducting leveling surveys in flat areas, the theoretical value of the mean square error δ per kilometer of leveling measurement is determined in the design... k Then, based on the number of stations z that need to be measured per kilometer, first apply the formula... Calculate δ 站 Then according to the formula To calculate the measurement accuracy, the distance between the two leveling observation points is l kilometers, and the formula is used. Calculate measurement accuracy; If the leveling plan involves forward and backward measurements, then the formula is used. The measured value δ is calculated from the measured data to determine the measurement accuracy. s ′, where h i往 h is the observed value of the cross-sectional elevation difference between the i-th leveling point and the (i+1)-th leveling point. i返 It is the observed value of the back-measured elevation difference between the i-th leveling point and the (i+1)-th leveling point.

6. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 1, characterized in that: In step h, the evaluation value δ for the InSAR cumulative deformation monitoring accuracy in the wide-area deformation region is calculated. L The formula for evaluating the accuracy of a single InSAR deformation monitoring in a wide-area deformation region, δ1, is as follows: In the formula: [·] is the Gaussian summation symbol. It is the symbol for taking the square root in mathematics.

7. The method for evaluating the accuracy of InSAR monitoring in a wide-area deformation region according to claim 6, characterized in that: When it is necessary to further reduce the impact of leveling measurement errors on accuracy estimation to meet higher accuracy requirements, the following formula is used to calculate the evaluation value δ of the InSAR cumulative deformation monitoring accuracy in a wide-area deformation region. L The evaluation value δ1 for the accuracy of InSAR deformation monitoring in a wide-area deformation region: When the leveling survey scheme involves forward and backward measurements, the measured value δ′ of the measurement accuracy can be used in the formula. s Replace δ s Perform accuracy estimation.

Images