Data processing and optimization method of InSAR deformation monitoring in small area

Abstract

The present invention provides a small area InSAR deformation monitoring data processing and optimization method, which belongs to the field of measurement data technology, and specifically comprises: dividing the deformation monitoring area into a plurality of sub-areas according to the position of the deformation monitoring area, determining the associated sub-areas in the sub-areas according to the deformation association between the historical interferometric measurement data of different sub-areas and the distribution position area, determining the deformation risk factors of different associated sub-areas and the concerned sub-areas according to the analysis results of the historical interferometric measurement data of the associated sub-areas, determining the analysis processing frequency of the interferometric measurement data based on the distribution data of the associated sub-areas and the deformation risk factor, and performing analysis processing on the interferometric measurement data of the distribution position area and the concerned sub-area according to the distribution processing frequency, thereby improving the reliability of the monitoring data processing.

Description

Technical Field

[0001] The invention belongs to the technical field of measurement data, and in particular relates to a small-area InSAR deformation monitoring data processing and optimization method. Background Art

[0002] Interferometric radar refers to synthetic aperture radar (InSAR) that uses interferometric measurement technology. It is a recently developed space earth observation technology. It is the product of the combination of traditional SAR remote sensing technology and radio astronomy interferometry technology. It can calculate the topography, landforms and slight changes on the surface of the target area. It is often used in the establishment of digital elevation models and crustal deformation detection. At the same time, how to realize the processing of InSAR deformation monitoring data has become a technical problem that needs to be solved urgently.

[0003] In view of the above technical problems, the invention patent CN202411124200.8 "System and method for monitoring subsidence of artificial islands at sea" has realized full coverage subsidence monitoring from point to line to surface by integrating GNSS PPP, leveling and InSAR technology, overcoming the limitations of a single monitoring technology and improving the comprehensiveness, continuity and high precision of monitoring. However, the above technical solution has the following technical problems:

[0004] In a small area centered on a nuclear power plant, due to the large amount of interferometric radar data, if real-time analysis and monitoring are performed based on the interferometric radar data, it is impossible to improve the reliability and accuracy of deformation monitoring processing while reducing the difficulty of analyzing the interferometric radar data.

[0005] In view of the above technical problems, the present invention proposes a small area InSAR deformation monitoring data processing and optimization method. Summary of the invention

[0006] To achieve the purpose of the present invention, the present invention adopts the following technical solutions:

[0007] In order to solve the above technical problems, the present invention provides the following technical solutions to achieve the purpose of the present invention:

[0008] According to one aspect of the present invention, a small area InSAR deformation monitoring data processing and optimization method is provided.

[0009] The small area InSAR deformation monitoring data processing and optimization method specifically includes:

[0010] S1 determines the historical interferometric measurement data of the interferometric radar in the distribution location area based on the distribution location area of ​​the nuclear power plant, and determines the monitoring area threshold of the distribution location area in combination with the change of the historical interferometric measurement data;

[0011] S2 determines the deformation monitoring area of ​​the distribution location area based on the monitoring area threshold, divides the deformation monitoring area into a plurality of sub-areas according to the location of the deformation monitoring area, and determines the associated sub-areas in the sub-areas according to the historical interferometric measurement data of different sub-areas and the deformation association of the distribution location area;

[0012] S3, determining deformation risk factors of different associated sub-regions and concerned sub-regions based on the analysis results of the historical interferometric measurement data of the associated sub-regions;

[0013] S4 determines the analysis and processing frequency of the interference measurement data based on the distribution data of the associated sub-region and the deformation risk factor, and uses the analysis and processing frequency to perform analysis and processing on the interference measurement data of the distribution position region and the focus sub-region.

[0014] The beneficial effects of the present invention are:

[0015] The monitoring area threshold of the distribution location area is determined by the historical interferometric measurement data of the interferometric radar in the distribution location area and the changes in the historical interferometric measurement data, thereby avoiding the technical problem of inaccurate deformation analysis and processing results caused by considering only the interferometric measurement data of the interferometric radar in the distribution location area. By determining the monitoring area threshold, a comprehensive analysis of the deformation data of the associated sub-areas of the distribution location area can be achieved, thereby improving the reliability of the analysis and processing results of the interferometric measurement data.

[0016] The analysis and processing frequency of the interferometric measurement data is determined based on the distribution data of the associated sub-areas and the deformation risk factors. This not only takes into account the differences in deformation risks of the distribution position areas due to differences in the sizes of the deformation risk factors of the associated sub-areas, but also takes into account the discreteness of the distribution of the associated sub-areas, thereby achieving a comprehensive consideration of the difficulty of evaluating the deformation risks of the position distribution areas. This improves the reliability of the deformation risk evaluation and processing while reducing the difficulty of analyzing and processing the deformation data.

[0017] A further technical solution is that the distribution location area of ​​the nuclear power plant is determined according to the distribution areas of the nuclear island equipment, conventional island equipment and auxiliary systems of the nuclear power plant.

[0018] A further technical solution is that the historical interference measurement data includes measurement times corresponding to the historical interference measurement data of the interference radar in the distribution location area.

[0019] A further technical solution is that the change of the historical interference measurement data includes the change amount of the historical interference measurement data between different measurement moments.

[0020] A further technical solution is that the method for determining the analysis and processing frequency of the interference measurement data is:

[0021] Determine the distribution position of the associated sub-region based on the distribution data of the associated sub-region, determine the interval distance between different associated sub-regions in combination with the distribution position of the associated sub-region, and determine the distribution dispersion coefficient of the associated sub-region by the interval distance between different associated sub-regions;

[0022] Determining a comprehensive regional risk factor of the associated sub-region based on the number of the associated sub-regions and the deformation risk factors of different associated sub-regions;

[0023] The monitoring risk coefficient is determined by the comprehensive regional risk factor and the distribution dispersion coefficient of the associated sub-region, and the analysis and processing frequency of the interferometric measurement data is determined based on the monitoring risk coefficient.

[0024] A further technical solution is to determine the analysis and processing frequency of the interferometric measurement data based on the monitoring risk coefficient, specifically:

[0025] Based on the monitoring risk coefficient, the risk coefficient interval of the interference measurement data is determined, and the analysis and processing frequency of the interference measurement data is determined through a preset processing frequency corresponding to the risk coefficient interval.

[0026] Other features and advantages will be described in the following description, and partly become apparent from the description, or understood by practicing the invention. The purpose and other advantages of the invention are realized and obtained by the structures particularly pointed out in the description and the drawings.

[0027] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The above and other features and advantages of the present invention will become more apparent by describing in detail exemplary embodiments thereof with reference to the attached drawings.

[0029] Figure 1 It is a flow chart of the small area InSAR deformation monitoring data processing and optimization method;

[0030] Figure 2 is a flow chart of a method for determining a monitoring area threshold of a distribution location area;

[0031] Figure 3 is a flow chart of a method for determining associated sub-regions in sub-regions;

[0032] Figure 4is a flowchart of a method for determining a sub-region of interest. DETAILED DESCRIPTION

[0033] Example embodiments will now be described more fully with reference to the accompanying drawings. However, example embodiments can be implemented in a variety of forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided so that the present invention will be comprehensive and complete and fully convey the concepts of the example embodiments to those skilled in the art. The same reference numerals in the figures represent the same or similar structures, and thus their detailed description will be omitted.

[0034] The terms "a", "an", "the", and "said" are used to indicate the presence of one or more elements / components / etc.; the terms "comprising" and "having" are used to express an open-ended inclusive meaning and mean that additional elements / components / etc. may be present in addition to the listed elements / components / etc.

[0035] To solve the above problems, according to one aspect of the present invention, Figure 1 As shown, according to one aspect of the present invention, a small area InSAR deformation monitoring data processing and optimization method is provided, which specifically includes:

[0036] S1 determines the historical interferometric measurement data of the interferometric radar in the distribution location area based on the distribution location area of ​​the nuclear power plant, and determines the monitoring area threshold of the distribution location area in combination with the change of the historical interferometric measurement data;

[0037] Furthermore, the distribution location area of ​​the nuclear power plant is determined based on the distribution areas of the nuclear island equipment, conventional island equipment, and auxiliary systems of the nuclear power plant.

[0038] It can be understood that the historical interference measurement data includes the measurement time corresponding to the historical interference measurement data of the interference radar in the distribution location area.

[0039] It should be noted that the change of the historical interference measurement data includes the change amount of the historical interference measurement data between different measurement moments.

[0040] Specifically, Figure 2 As shown, the method for determining the monitoring area threshold of the distribution location area is:

[0041] Determining the deformation state of the distribution location area based on the analysis results of the historical interferometric measurement data of the interferometric radar of the distribution location area;

[0042] Based on the changes in the historical interferometric measurement data, determine the changes in the deformation state of the distribution location area within a preset period of time, and determine the duration of the abnormal deformation state within the preset period of time based on the changes in the deformation state within the preset period of time;

[0043] The comprehensive risk coefficient of the distribution location area is determined by the duration of the abnormal deformation state, and the monitoring area threshold of the distribution location area is determined based on the comprehensive risk coefficient.

[0044] Further, determining a monitoring area threshold of the distribution location area based on the comprehensive risk coefficient specifically includes:

[0045] Based on the comprehensive risk coefficient, a preset monitoring area threshold of the distribution location area under the comprehensive risk coefficient is determined, and the monitoring area threshold of the distribution location area is determined based on the preset monitoring area threshold.

[0046] It should also be noted that the comprehensive risk coefficient is determined based on the ratio of the duration of abnormal deformation state within a preset period of time to the preset duration.

[0047] Optionally, the method for determining the monitoring area threshold of the distribution location area includes steps S11-S13, specifically:

[0048] S11 determines the deformation state of the distribution location area based on the analysis result of the historical interferometric measurement data of the interferometric radar of the distribution location area, and determines the duration of the abnormal deformation state within the preset time period based on the analysis result of the deformation state, and determines the regional deformation risk coefficient of the distribution location area according to the duration;

[0049] S12 determines the change of the deformation state of the distribution location area within a preset period of time based on the change of the historical interference measurement data, and determines the change amount of the deformation state between different adjacent moments within the preset period of time based on the change of the deformation state within the preset period of time, determines the change moment based on the change amount of the deformation state between different adjacent moments, and determines the change risk coefficient of the distribution location area based on the proportion of the number of change moments;

[0050] S13 determines the comprehensive risk coefficient of the distribution location area through the change risk coefficient and the regional deformation risk coefficient, and determines the monitoring area threshold of the distribution location area based on the comprehensive risk coefficient.

[0051] Optionally, the above step S11 includes steps S111-S114, which are specifically:

[0052] S111 determines the deformation state of the distribution location area based on the analysis result of the historical interference measurement data of the interference radar of the distribution location area. When the deformation state of the distribution location area is not abnormal, the monitoring area threshold of the distribution location area is directly determined by using the preset area. When the deformation state of the distribution location area is abnormal, the process proceeds to step S112.

[0053] S112: determining the abnormal coefficient of the deformation state based on the analysis result of the deformation state; when the abnormal coefficient of the deformation state is not within the preset abnormal coefficient interval, directly using the specified area threshold to determine the monitoring area threshold of the distribution location area; when the abnormal coefficient of the deformation state is within the preset abnormal coefficient interval, proceeding to step S113;

[0054] S113 determines the duration of abnormal deformation state within the preset period based on the analysis result of the deformation state. When the duration of abnormal deformation state within the preset period is greater than the preset duration threshold, the monitoring area threshold of the distribution location area is directly determined by using the specified area threshold. When the duration of abnormal deformation state within the preset period is not greater than the preset duration threshold, the process proceeds to step S114.

[0055] S114 determines the regional deformation risk coefficient of the distribution location area according to the duration, and proceeds to step S12.

[0056] Optionally, the above step S12 includes steps S121-S123, specifically:

[0057] S121 determines the change of the deformation state of the distribution location area within the preset period based on the change of the historical interference measurement data, and determines the deformation state change amount between different adjacent moments within the preset period based on the change of the deformation state within the preset period. When there is no change moment within the preset period, the monitoring area threshold of the distribution location area is determined based on the regional deformation risk coefficient of the distribution location area. When there is a change moment within the preset period, the process proceeds to step S122.

[0058] S122 determines the change time by the deformation state change amount between different adjacent moments. When the number of change moments is greater than the preset number of moments, the monitoring area threshold of the distribution location area is directly determined by using the specified area threshold. When the number of change moments is not greater than the preset number of moments, the process proceeds to step S123.

[0059] S123 determines the change risk coefficient of the distribution location area based on the quantity proportion of the change time, and proceeds to step S13.

[0060] S2 determines the deformation monitoring area of ​​the distribution location area based on the monitoring area threshold, divides the deformation monitoring area into a plurality of sub-areas according to the location of the deformation monitoring area, and determines the associated sub-areas in the sub-areas according to the historical interferometric measurement data of different sub-areas and the deformation association of the distribution location area;

[0061] It should be noted that the deformation monitoring area is a monitoring area with the distribution position area as the center and the monitoring area threshold as the basis.

[0062] Furthermore, the deformation monitoring area is divided into a plurality of sub-areas, specifically including:

[0063] The deformation monitoring area is divided into a plurality of sub-areas according to the position of the deformation monitoring area in the distribution position area.

[0064] It will be appreciated that the locations are based on southeast, northeast, northwest and southwest.

[0065] Further, such as Figure 3 As shown, the method for determining the associated sub-regions in the sub-regions is:

[0066] Determining a data change time when the historical interferometry data of the sub-region changes based on an analysis result of the historical interferometry data of the sub-region;

[0067] Determine the data change time when the distribution position area changes by the change of the historical interference measurement data of the distribution position area at the data change time, and use it as the associated change time;

[0068] The comprehensive correlation coefficient between the sub-region and the distribution location region is determined according to the proportion of the correlation change moments in the data change moments, and whether the sub-region is a correlation sub-region is determined according to the comprehensive correlation coefficient.

[0069] Optionally, a method for determining an associated sub-region in the sub-region is:

[0070] S21, based on the analysis result of the historical interferometry data of the sub-region, determining the data change time when the historical interferometry data of the sub-region changes;

[0071] S22 determines the data change correlation coefficient at the data change time by the change of the historical interferometric measurement data of the distribution position area and the change of the historical interferometric measurement data of the sub-area at the data change time;

[0072] S23 determines a comprehensive correlation coefficient between the sub-region and the distribution location region using data change correlation coefficients at different data change moments, and determines whether the sub-region is a correlated sub-region using the comprehensive correlation coefficient.

[0073] Furthermore, the data change correlation coefficient at the data change moment is determined according to the ratio of the change proportions of the historical interference measurement data of the distribution location area and the sub-area.

[0074] It should also be noted that when the comprehensive correlation coefficient of the sub-region is greater than a preset correlation threshold, the sub-region is determined to be a correlation sub-region.

[0075] S3, determining deformation risk factors of different associated sub-regions and concerned sub-regions based on the analysis results of the historical interferometric measurement data of the associated sub-regions;

[0076] Specifically, Figure 4 As shown, the method for determining the focus sub-region is:

[0077] Determine the change of the historical interferometry data of the associated sub-region within a preset period of time based on the analysis result of the historical interferometry data of the associated sub-region;

[0078] Based on the change of the historical interferometric measurement data of the associated sub-area within the preset time period, the change amount of the historical interferometric measurement data at different times is determined, and the risk change time of the associated sub-area within the preset time period is determined by the change amount of the historical interferometric measurement data at different times;

[0079] The deformation risk factor of the associated sub-region is determined by the proportion of the number of risk change moments of the associated sub-region within a preset time period, and whether the associated sub-region is a concerned sub-region is determined by the deformation risk factor.

[0080] Furthermore, the risk change moment is the moment when the change amount of the historical interference measurement data is greater than a preset change amount threshold.

[0081] Optionally, the method for determining the focus sub-region includes steps S31-S33, specifically:

[0082] S31, based on the analysis result of the historical interferometry data of the associated sub-region, determining the change of the historical interferometry data of the associated sub-region within a preset period of time;

[0083] Optionally, before entering step S32, it is also necessary to determine based on the changes in the historical interference measurement data of the associated sub-area within the preset time period that if there is no moment in the preset time period when the change is not within the preset change range, then it is determined that the associated sub-area does not belong to the focus sub-area. When there is a moment in the preset time period when the change is within the preset change range, proceed to step S32.

[0084] S32: determining the change amount of the historical interferometry data at different times based on the change of the historical interferometry data of the associated sub-area within the preset time period, and determining the risk change time by the change amount of the historical interferometry data at different times;

[0085] Optionally, the above step S32 includes steps S321-S323, specifically:

[0086] S321 determines the change amount of the historical interference measurement data at different moments based on the change of the historical interference measurement data of the associated sub-area within the preset time period, and determines the risk change moment by the change amount of the historical interference measurement data at different moments. When the number of risk change moments is less than the preset number of risk moments, proceed to step S322; when the number of risk change moments is not less than the preset number of risk moments, proceed to step S33;

[0087] S322: When there is no moment when the change amount does not meet the requirement, it is determined that the associated sub-region does not belong to the focus sub-region. When there is a moment when the change amount does not meet the requirement, the process proceeds to step S33.

[0088] S33 divides the preset time period into different unit time periods, determines the deformation risk factor of the associated sub-area by the number of risk change moments in different unit time periods and the change amount of different risk change moments, and determines whether the associated sub-area is a focus sub-area by the deformation risk factor.

[0089] Optionally, the above step S33 includes steps S331-S333, specifically:

[0090] S331 divides the preset time period into different unit time periods. When there is a unit time period whose number ratio of risk change moments does not meet the requirement, the process proceeds to step S332. When there is no unit time period whose number ratio of risk change moments does not meet the requirement, it is determined that the associated sub-area does not belong to the concerned sub-area.

[0091] S332 determines the period deformation risk factors of different unit periods according to the number of risk change moments in different unit periods and the change amounts of different risk change moments. When there is no unit period whose period deformation risk factor does not meet the requirement, it is determined that the associated sub-area does not belong to the attention sub-area. When there is a unit period whose period deformation risk factor does not meet the requirement, the process proceeds to step S333.

[0092] S333: Determine a deformation risk factor of the associated sub-region based on the period deformation risk factors of the different unit time periods, and determine whether the associated sub-region is a focus sub-region according to the deformation risk factor.

[0093] S4 determines the analysis and processing frequency of the interference measurement data based on the distribution data of the associated sub-region and the deformation risk factor, and uses the analysis and processing frequency to perform analysis and processing on the interference measurement data of the distribution position region and the focus sub-region.

[0094] Specifically, the method for determining the analysis and processing frequency of the interference measurement data is:

[0095] Determine the distribution position of the associated sub-region based on the distribution data of the associated sub-region, determine the interval distance between different associated sub-regions in combination with the distribution position of the associated sub-region, and determine the distribution dispersion coefficient of the associated sub-region by the interval distance between different associated sub-regions;

[0096] Determining a comprehensive regional risk factor of the associated sub-region based on the number of the associated sub-regions and the deformation risk factors of different associated sub-regions;

[0097] The monitoring risk coefficient is determined by the comprehensive regional risk factor and the distribution dispersion coefficient of the associated sub-region, and the analysis and processing frequency of the interferometric measurement data is determined based on the monitoring risk coefficient.

[0098] Further, the analysis and processing frequency of the interferometric measurement data is determined based on the monitoring risk coefficient, specifically:

[0099] Based on the monitoring risk coefficient, the risk coefficient interval of the interference measurement data is determined, and the analysis and processing frequency of the interference measurement data is determined through a preset processing frequency corresponding to the risk coefficient interval.

[0100] In another possible embodiment:

[0101] Optionally, the above step S22 includes steps S221-S223, specifically:

[0102] S221 determines the data change time when the distribution position area changes according to the change of the historical interferometric measurement data of the distribution position area at the data change time, and uses it as the associated change time. When there is no associated change time, it is determined that the sub-area does not belong to the associated sub-area. When there is an associated change time, it goes to step S222;

[0103] S222: taking the ratio of the number of the associated change moments to the number of the data change moments as the ratio of the number of change moments; when the ratio of the number of change moments is greater than the preset ratio, determining that the sub-region belongs to the associated sub-region; and when the ratio of the number of change moments is not greater than the preset ratio, proceeding to step S223;

[0104] S223 determines the data change correlation coefficient at the data change moment through the changes in the historical interference measurement data of the distribution position area and the changes in the historical interference measurement data of the sub-area at the data change moment, and proceeds to step S23.

[0105] Optionally, the above step S23 includes steps S231-S233, specifically:

[0106] S231 determines that the sub-region does not belong to the associated sub-region when there is a data change correlation coefficient greater than the preset correlation coefficient threshold value at different data change moments based on the data change correlation coefficients at different data change moments. When there is a data change moment with a data change correlation coefficient greater than the preset correlation coefficient threshold value, the process proceeds to step S232;

[0107] S232: when the data change correlation coefficient is greater than the preset correlation coefficient threshold, the proportion of the number of data change moments in the data change moments is greater than the preset correlation moment proportion, it is determined that the sub-region belongs to the associated sub-region; when the data change correlation coefficient is greater than the preset correlation coefficient threshold, the proportion of the number of data change moments in the data change moments is not greater than the preset correlation moment proportion, the process proceeds to step S233;

[0108] S233 determines the comprehensive correlation coefficient between the sub-region and the distribution location region by using the data change correlation coefficients at different data change moments, and determines whether the sub-region is a correlated sub-region by using the comprehensive correlation coefficient.

[0109] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device, equipment, and non-volatile computer storage medium embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0110] The above is a description of a specific embodiment of the specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0111] The above description is only one or more embodiments of this specification and is not intended to limit this specification. For those skilled in the art, one or more embodiments of this specification may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of one or more embodiments of this specification shall be included in the scope of the claims of this specification.

Claims

1. A small area InSAR deformation monitoring data processing and optimization method, characterized in that: Specifically include: Based on the distribution location area of ​​the nuclear power plant, determine the historical interference measurement data of the interference radar in the distribution location area, and determine the monitoring area threshold of the distribution location area in combination with the change of the historical interference measurement data; Based on the monitoring area threshold, a deformation monitoring area of ​​the distribution location area is determined, and the deformation monitoring area is divided into a plurality of sub-areas according to the location of the deformation monitoring area, and the associated sub-areas in the sub-areas are determined according to the historical interferometric measurement data of different sub-areas and the deformation association of the distribution location area; Determine deformation risk factors of different associated sub-regions and concerned sub-regions based on analysis results of historical interferometric measurement data of the associated sub-regions; The analysis and processing frequency of the interference measurement data is determined based on the distribution data of the associated sub-region and the deformation risk factor, and the analysis and processing frequency is used to perform analysis and processing on the interference measurement data of the distribution position region and the focus sub-region.

2. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The distribution location area of ​​the nuclear power plant is determined according to the distribution areas of the nuclear island equipment, conventional island equipment, and auxiliary systems of the nuclear power plant.

3. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The historical interference measurement data includes the measurement time corresponding to the historical interference measurement data of the interference radar in the distribution location area.

4. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The change of the historical interference measurement data includes the change amount of the historical interference measurement data between different measurement moments.

5. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The method for determining the monitoring area threshold of the distribution location area is: Determining the deformation state of the distribution location area based on the analysis results of the historical interferometric measurement data of the interferometric radar of the distribution location area; Based on the changes in the historical interferometric measurement data, determine the changes in the deformation state of the distribution location area within a preset period of time, and determine the duration of the abnormal deformation state within the preset period of time based on the changes in the deformation state within the preset period of time; The comprehensive risk coefficient of the distribution location area is determined by the duration of the abnormal deformation state, and the monitoring area threshold of the distribution location area is determined based on the comprehensive risk coefficient.

6. The small area InSAR deformation monitoring data processing and optimization method according to claim 5, characterized in that: Determining the monitoring area threshold of the distribution location area based on the comprehensive risk coefficient specifically includes: Based on the comprehensive risk coefficient, a preset monitoring area threshold of the distribution location area under the comprehensive risk coefficient is determined, and the monitoring area threshold of the distribution location area is determined based on the preset monitoring area threshold.

7. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The deformation monitoring area is a monitoring area with the distribution position area as the center and a monitoring area threshold as the basis.

8. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The deformation monitoring area is divided into a plurality of sub-areas, specifically including: The deformation monitoring area is divided into a plurality of sub-areas according to the position of the deformation monitoring area in the distribution position area.

9. The small area InSAR deformation monitoring data processing and optimization method according to claim 1, characterized in that: The method for determining the analysis and processing frequency of the interferometric measurement data is: Determine the distribution position of the associated sub-region based on the distribution data of the associated sub-region, determine the interval distance between different associated sub-regions in combination with the distribution position of the associated sub-region, and determine the distribution dispersion coefficient of the associated sub-region by the interval distance between different associated sub-regions; Determining a comprehensive regional risk factor of the associated sub-region based on the number of the associated sub-regions and the deformation risk factors of different associated sub-regions; The monitoring risk coefficient is determined by the comprehensive regional risk factor and the distribution dispersion coefficient of the associated sub-region, and the analysis and processing frequency of the interferometric measurement data is determined based on the monitoring risk coefficient.

10. The small area InSAR deformation monitoring data processing and optimization method according to claim 9, characterized in that: The frequency of analyzing and processing the interferometric measurement data is determined based on the monitoring risk coefficient, specifically: Based on the monitoring risk coefficient, the risk coefficient interval of the interference measurement data is determined, and the analysis and processing frequency of the interference measurement data is determined through a preset processing frequency corresponding to the risk coefficient interval.

Images