Mining subsidence boundary identification method and device

Abstract

The invention discloses a method and device for identifying mining subsidence boundaries in a mining area, and relates to the technical field of mining subsidence monitoring. The method comprises obtaining N times before and after the mining start time. 1 The synthetic aperture radar image data is used as the first image data to obtain N images within the set time of the subsidence area. 2 The synthetic aperture radar image data is used as the second image data; the first image data and the second image data are processed by SARscape software to obtain an image map in ArcGIS data format. 1 and ImageMap 2 ; Map the image 1 and the image map 2 Input into ArcGIS to obtain 3D visualization map 1‑3D With Map 2‑3D ; 3D visualization map 1‑3D With Map 2‑3D The boundary of the public area is determined as the boundary of the stable area. The problem of determining the boundary of the stable area is solved.

Description

Technical Field

[0001] The invention relates to the technical field of coal mining subsidence monitoring in mining areas, and in particular to a method and device for identifying mining subsidence boundaries in mining areas. Background Art

[0002] Mining subsidence refers to the phenomenon that underground mineral resource mining causes rock movement and surface subsidence. After underground mineral resources are mined, the original mechanical equilibrium state of the rock mass around the mining area is destroyed, the rock layer moves and deforms, and spreads to the surface. Houses, buildings, railways, etc. located within the mining area will be deformed or damaged, threatening the lives and property of local residents. Therefore, it is necessary to monitor mining subsidence in mining areas.

[0003] At present, the short baseline set interferometry (SBAS-InSAR) technology is proposed by Berardino and Lanari et al. on the basis of synthetic aperture radar differential interferometry (D-InSAR). All radar (SAR) images are combined into several sets (baseline sets with short time and space baselines) according to the interference requirements, and the terrain deformation sequence of each small set is obtained by the least squares method. Then, the singular value decomposition method (SVD) is used to jointly solve the rank-deficient equation of multiple small baseline sets to obtain the deformation of ground subsidence over a period of time. Unlike the permanent scatterer interferometry (PS-InSAR), SBAS-InSAR does not require the same main image, and selects interferometric pairs with short baselines to combine into different sets, and then uses the singular value decomposition method to jointly solve. SBAS-InSAR not only limits the temporal and spatial incoherence and improves the utilization rate of images, but also enriches the observation information of images. It has been a research hotspot in the field of mining subsidence in recent years.

[0004] Coal mining subsidence areas are divided into unstable and stable areas. For unstable areas, the surface deformation is still changing, and conventional radar interferometry can accurately monitor the deformation value and subsidence boundary. However, for stable areas, the underground rock formation tends to be stable, and the surface deformation is small. At this time, if the selected time baseline is short, the surface deformation will not be monitored when using conventional radar interferometry technology to monitor it; if the selected time baseline is long, it will cause data time incoherence and cannot obtain accurate results.

[0005] In the related art, the invention patent application with application number 201510769347.7 discloses a method for multi-source data monitoring of the spatiotemporal characteristics of mining area deformation and identification of cross-border mining. First, the synthetic aperture radar (SAR) is used to obtain and optimize the processing data, and then the SAR data monitoring data of multiple satellite platforms are fused to obtain the settlement result map. In the fused settlement result map, the "large deformation area" or "decoherence area" in the measured mining area is supplemented with monitoring, and the coordinates of the multi-satellite platform data monitoring results are unified. Finally, the spatiotemporal characteristics of the surface deformation of the mining area are inverted, and finally the ArcGIS platform is used to identify the cross-border mining area. The method can effectively invert the spatiotemporal characteristics of the surface settlement of the mining area.

[0006] However, the subsidence boundary determined by this method is the boundary of the non-stable subsidence area after mining, and it does not provide a solution for identifying the boundary of the stable subsidence area formed after a period of time (one year, two years or longer) after mining. Summary of the invention

[0007] The technical problem to be solved by the present invention is how to improve the extraction accuracy of mining subsidence information in a stable area.

[0008] The present invention solves the above technical problems through the following technical means:

[0009] On the one hand, a method for identifying mining subsidence boundaries in a mining area is used, the method comprising:

[0010] Get the number of N before and after the mining start time 1 The synthetic aperture radar image data is used as the first image data, and N images within the set time of the subsidence area are obtained. 2 using the scene synthetic aperture radar image data as the second image data;

[0011] The first image data and the second image data are processed by SARscape software to obtain an image map in ArcGIS data format. 1 and ImageMap 2 ;

[0012] The image Map 1 and the image map 2 Input into ArcGIS to obtain 3D visualization map 1-3D With Map 2-3D ;

[0013] The 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area;

[0014] From the 3D visualization map 2-3DGet the boundary of the unstable area.

[0015] The present invention selects SAR data before and after mining in a mining area as comparison and includes SAR data within a set time of a subsidence area, thereby determining the mining range of a stable subsidence area in a mining area, and inverts an image in an ArcGIS data format by processing the two types of data, and inputs the subsidence map obtained by the SBAS technology into ArcGIS software for post-processing, so as to more intuitively illustrate the subsidence boundary; the present invention is particularly suitable for situations where the mining subsidence boundary cannot be determined due to the long mining time in a mining area, and can avoid the technical defects of traditional measurement methods such as being affected by weather, environment, etc. and the temporal incoherence and spatial incoherence of D-InSAR, thereby improving the work efficiency and the accuracy of terrain deformation monitoring, so as to better provide technical support for production and life in mining areas.

[0016] Furthermore, the first image data and the second image data are processed by SARscape software to obtain an image map in ArcGIS data format. 1 and ImageMap 2 ,include:

[0017] Interference image pairs are respectively performed on the first image data and the second image data to obtain a first interference image pair and a second interference image pair.

[0018] Interference processing is performed on the first interference image pair and the second interference image to obtain a first processed image and a second processed image.

[0019] Track refinement and re-flattening are performed on the first processed image and the second processed image respectively to obtain first track data and second track data after track refinement.

[0020] Inverting the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result;

[0021] After geocoding the first SBAS result and the second SBAS result, the image Map is converted into an ArcGIS data format. 1 and the image map 2 .

[0022] Further, the using SBAS to invert the first orbit data and the second orbit data to obtain a first SBAS result and a second SBAS result includes:

[0023] Inverting the first orbit data and the second orbit data using the SBAS to estimate corresponding deformation rate and residual topography;

[0024] Based on the corresponding deformation rate, the atmospheric filter is customized, the atmospheric phase is estimated and removed, and the displacement results on the corresponding time series are obtained as the first SBAS result and the second SBAS result.

[0025] Furthermore, the N 1 =N 2 ≥15.

[0026] In a second aspect, the present invention adopts a mining subsidence boundary identification device, the device comprising:

[0027] Acquisition module, used to obtain N before and after the mining start time 1 The synthetic aperture radar image data is used as the first image data, and N images within the set time of the subsidence area are obtained. 2 using the scene synthetic aperture radar image data as the second image data;

[0028] A processing module is used to process the first image data and the second image data respectively using SARscape software to obtain an image map in ArcGIS data format. 1 and ImageMap 2 ;

[0029] A visualization module is used to map the image 1 and the image map 2 Input into ArcGIS to obtain 3D visualization map 1-3D With Map 2-3D ;

[0030] The first boundary determination module is used to map the 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area;

[0031] The second boundary determination module is used to determine the boundary from the 3D visualization map 2-3D Get the boundary of the unstable area.

[0032] Furthermore, the processing module includes:

[0033] The pairing unit is used to pair the first image data with the second image data to obtain a first interference image pair and a second interference image pair.

[0034] The interference unit is used to perform interference processing on the first interference image pair and the second interference image respectively to obtain a first processed image and a second processed image.

[0035] The track refining unit is used to perform track refinement and re-flattening on the first processed image and the second processed image respectively to obtain the first track data and the second track data after track refinement.

[0036] an inversion unit, configured to invert the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result;

[0037] A coding conversion unit is used to convert the first SBAS result and the second SBAS result into the image Map in ArcGIS data format after geocoding. 1 and the image map 2 .

[0038] Furthermore, the inversion unit comprises:

[0039] an estimation subunit, configured to estimate corresponding deformation rate and residual topography by inverting the first orbit data and the second orbit data using the SBAS;

[0040] The displacement calculation subunit is used to customize the atmospheric filtering, estimate and remove the atmospheric phase based on the corresponding deformation rate, and obtain the displacement results on the corresponding time series as the first SBAS result and the second SBAS result.

[0041] Furthermore, the N 1 =N 2 ≥15.

[0042] In a third aspect, the present invention adopts a device, which includes a memory and a processor; wherein the processor runs a program corresponding to the executable program code by reading the executable program code stored in the memory, so as to implement the method as described above.

[0043] In a fourth aspect, the present invention adopts a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method as described above.

[0044] The advantages of the present invention are:

[0045] (1) The present invention selects SAR data before and after mining in a mining area as comparison and includes SAR data within a set time of the subsidence area, thereby determining the mining range of the stable subsidence area of ​​the mining area. By processing these two types of data, an image in ArcGIS data format is inverted, and the subsidence map obtained by SBAS technology is input into ArcGIS software for post-processing, so as to more intuitively illustrate the subsidence boundary. It is particularly suitable for situations where the mining subsidence boundary cannot be determined due to the long mining time in the mining area, and can avoid the technical defects of the traditional measurement method such as being affected by weather, environment, etc. and the temporal incoherence and spatial incoherence of D-InSAR, thereby improving the work efficiency and the accuracy of terrain deformation monitoring, so as to better provide technical support for production and life in the mining area.

[0046] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] Figure 1 It is a flow chart of a method for identifying mining subsidence boundaries in a mining area in Embodiment 1 of the present invention;

[0048] Figure 2 This is an overall flow chart of the method for identifying mining subsidence boundaries in a mining area in Embodiment 1 of the present invention;

[0049] Figure 3 This is a schematic diagram of mining subsidence in a mining area in Embodiment 1 of the present invention;

[0050] Figure 4 It is a structural diagram of a mining subsidence boundary identification device in a second embodiment of the present invention. DETAILED DESCRIPTION

[0051] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in combination with the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0052] like Figure 1 As shown, the first embodiment of the present invention proposes a method for identifying mining subsidence boundaries in a mining area, and the method comprises the following steps:

[0053] S10, obtain N before and after the mining start time 1 The synthetic aperture radar image data is used as the first image data, and N images within the set time of the subsidence area are obtained. 2 using the scene synthetic aperture radar image data as the second image data;

[0054] It should be noted that, in this embodiment, mining data of the mining area is collected to determine the mining start time of the mining area, and the total N including the coal seam before and after mining is obtained according to the mining start time. 1 Based on the mining status of the mining area, we can obtain N data including the subsidence area in the past two years. 2 Scene SAR image data.

[0055] It should be understood that the set time value of 2 years in this embodiment is only for example. Those skilled in the art may set the set time to one year, two years or longer according to actual conditions.

[0056] This embodiment obtains SAR data before mining in the mining area for comparison, so as to determine the mining range of the stable subsidence area of ​​the mining area, and obtains SAR data within a set time of the subsidence area. The set time refers to the time when the boundary of the stable subsidence area can be formed after mining in the mining area, while identifying the boundary of the stable subsidence area.

[0057] S20, using SARscape software to process the first image data and the second image data respectively to obtain an image map in ArcGIS data format 1 and ImageMap 2 ;

[0058] S30, Map the image 1 and the image map 2 Input into ArcGIS to obtain 3D visualization map 1-3D With Map 2-3D ;

[0059] S40, mapping the 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area;

[0060] S50, from the 3D visualization map 2-3D Get the boundary of the unstable area.

[0061] It should be noted that by selecting SAR data before and after mining in the mining area as a comparison and including SAR data within a set time in the subsidence area, the mining scope of the stable subsidence area of ​​the mining area is determined. By processing these two types of data, images in ArcGIS data format are inverted and input into ArcGIS software to post-process the subsidence map obtained by SBAS technology to more intuitively illustrate the subsidence boundary. It is especially suitable for situations where the mining subsidence boundary cannot be determined due to long mining time in the mining area. It can avoid the influence of weather, environment, etc. in traditional measurement methods and technical defects such as time incoherence and space incoherence of D-InSAR, thereby improving work efficiency and terrain deformation monitoring accuracy, so as to better provide technical support for production and life in the mining area.

[0062] It should be noted that this embodiment uses ArcGIS software to post-process the subsidence map obtained by SBAS technology, and uses the three-dimensional display function of ArcGIS software to more intuitively determine the real mining subsidence boundary. This is completely different from the related art that uses ArcGIS software to post-process the subsidence area map, which is mainly reflected in the spatial analysis of the subsidence rate and then the inversion of the subsidence boundary.

[0063] In one embodiment, if Figure 2 As shown, the step S20 specifically includes the following steps:

[0064] S21 , performing interference image pair matching on the first image data and the second image data to obtain a first interference image pair and a second interference image pair.

[0065] It should be noted that, in this embodiment, the input data is paired with interference images, and the input N (N 1 Or N 2 ) scene data, the maximum number of pairings that can be obtained is (N*(N-1)) / 2, and the connection diagram generation tool will select the optimal combination for pairing.

[0066] S22, performing interference processing on the first interference image pair and the second interference image respectively to obtain a first processed image and a second processed image;

[0067] It should be noted that in this embodiment, all paired interferometric image pairs are subjected to interference processing, and all data pairs are registered to the super master image from coherence generation, deflating, filtering and phase unwrapping.

[0068] S23, performing track refinement and re-flattening on the first processed image and the second processed image respectively to obtain first track data and second track data after track refinement;

[0069] It should be noted that this embodiment estimates and removes the residual constant phase and the phase ramp that still exists after unwrapping through trajectory refinement and re-flattening.

[0070] S24, inverting the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result;

[0071] S25, geocoding the first SBAS result and the second SBAS result, and converting them into the image Map in ArcGIS data format 1 and the image map 2 .

[0072] It should be noted that, in this embodiment, the SBAS result is geocoded and converted into an image map in ArcGIS data format. 1 With Map 2 ), which can then be reprocessed in ArcGIS.

[0073] This embodiment processes SAR data by using the SBAS process processing function of SARscape software to obtain a mining area subsidence map, and uses the three-dimensional display function of ArcGIS software to more intuitively determine the real mining subsidence boundary. The advantage of using SARscape software is that SAR data can be processed in a process-based manner, the operation is simple, and the results are accurate; the advantage of using SBAS technology is that it not only limits the temporal and spatial incoherence, improves the utilization rate of images, but also enriches the observation information of images.

[0074] In one embodiment, the step S24 includes the following steps:

[0075] Inverting the first orbit data and the second orbit data using the SBAS to estimate corresponding deformation rate and residual topography;

[0076] Based on the corresponding deformation rate, the atmospheric filter is customized, the atmospheric phase is estimated and removed, and the displacement results on the corresponding time series are obtained as the first SBAS result and the second SBAS result.

[0077] It should be noted that the inversion in this embodiment is divided into two steps: the first step is to estimate the deformation rate and residual terrain for the first time; the second step is to calculate the displacement in the time series. Based on the deformation rate obtained in the first step of SBAS inversion, a customized atmospheric filter is performed to estimate and remove the atmospheric phase to obtain a more pure final displacement result in the time series. The two options of atmospheric high pass and atmospheric low pass estimate the atmospheric impact, and finally subtract these atmospheric parts from the measured displacement at each time.

[0078] In one embodiment, the N1 =N 2 ≥15.

[0079] It should be noted that, in this embodiment, N 1 、N 2 The first image data and the second image data are data obtained by imaging the same surface observation scene by radar and completing image registration.

[0080] Specifically, enter Map in ArcGIS 1 With Map 2 , respectively in Map 1 With Map 2 Select the subsidence area, extract the subsidence area and overlay the mining plane and other data, generate contour lines and perform 3D visualization to obtain a 3D visualization map 1-3D With Map 2-3D . Figure 3 A schematic diagram of mining subsidence is given. Figure 3 In the figure, A indicates the point before the ground subsidence, B indicates the point before the ground subsidence, C indicates the point before the ground subsidence, O indicates the boundary of the stable subsidence area, A' indicates the point after the ground subsidence of point A, B' indicates the point after the ground subsidence of point B, Q indicates the boundary of the unstable subsidence area formed after mining for T1 time, B' indicates the position of the final subsidence of point B, C' indicates the point after the ground subsidence of point C, and Q' indicates the boundary of the unstable subsidence area formed after mining for T2 time; OQ indicates Map 1-3D Range, OQ' represents Map 2-3D Compare to Map 1-3D With Map 2-3D The boundary of the common area between the two is the boundary of the stable settlement area, and the boundary of the non-stable settlement area can be obtained from Map 2-3D Obtained from the figure.

[0081] It should be noted that because Map 2-3D The map contains the subsidence area map of the selected experiment end time. The area that still subsides during this time is the unstable subsidence area, so it can be directly obtained from the Map 2-3D The boundary of the non-stable area can be observed in the figure.

[0082] like Figure 4 As shown, the second embodiment of the present invention provides a mining subsidence boundary identification device, the device comprising:

[0083] Acquisition module 10, used to obtain N before and after the mining start time 1 The synthetic aperture radar image data is used as the first image data, and N images within the set time of the subsidence area are obtained. 2using the scene synthetic aperture radar image data as the second image data;

[0084] The processing module 20 is used to process the first image data and the second image data respectively using SARscape software to obtain an image map in ArcGIS data format. 1 and ImageMap 2 ;

[0085] The visualization module 30 is used to map the image 1 and the image map 2 Input into ArcGIS to obtain 3D visualization map 1-3D With Map 2-3D ;

[0086] The first boundary determination module 40 is used to determine the 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area;

[0087] The second boundary determination module 50 is used to determine the boundary from the 3D visualization map 2-3D Get the boundary of the unstable area.

[0088] In one embodiment, the processing module 20 includes:

[0089] The pairing unit is used to pair the first image data with the second image data to obtain a first interference image pair and a second interference image pair.

[0090] The interference unit is used to perform interference processing on the first interference image pair and the second interference image respectively to obtain a first processed image and a second processed image.

[0091] The track refining unit is used to perform track refinement and re-flattening on the first processed image and the second processed image respectively to obtain the first track data and the second track data after track refinement.

[0092] an inversion unit, configured to invert the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result;

[0093] A coding conversion unit is used to convert the first SBAS result and the second SBAS result into the image Map in ArcGIS data format after geocoding. 1 and the image map 2 .

[0094] In one embodiment, the inversion unit comprises:

[0095] an estimation subunit, configured to estimate corresponding deformation rate and residual topography by inverting the first orbit data and the second orbit data using the SBAS;

[0096] The displacement calculation subunit is used to customize the atmospheric filtering based on the corresponding deformation rate, estimate and remove the atmospheric phase, and obtain the displacement results on the corresponding time series as the first SBAS result and the second SBAS result.

[0097] In one embodiment, the N 1 =N 2 ≥15.

[0098] It should be noted that other embodiments of the mining subsidence boundary identification device of the present invention or implementation methods can refer to the above-mentioned method embodiments, which are not repeated here.

[0099] A third embodiment of the present invention provides a device, which includes a memory and a processor; wherein the processor runs a program corresponding to the executable program code by reading the executable program code stored in the memory, so as to implement the method described in the above embodiment.

[0100] A fourth embodiment of the present invention provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the method described in the above embodiment is implemented.

[0101] It should be noted that the logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be specifically implemented in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection portion with one or more wirings (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0102] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0103] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0104] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0105] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A method for identifying mining subsidence boundaries in a mining area, characterized in that: The method comprises: Acquire N1 scene synthetic aperture radar image data before and after the start time of mining as first image data, and acquire N2 scene synthetic aperture radar image data within a set time including the subsidence area as second image data; The first image data and the second image data are processed respectively by using SARscape software to obtain an image Map1 and an image Map2 in ArcGIS data format; The image Map1 and the image Map2 are input into ArcGIS to obtain 3D visualization maps Map 1-3D With Map 2-3D ; The 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area; From the 3D visualization map 2-3D Get the boundary of the unstable area.

2. The method for identifying mining subsidence boundaries according to claim 1, characterized in that: The method of using SARscape software to process the first image data and the second image data respectively to obtain an image Map1 and an image Map2 in ArcGIS data format includes: Perform interference image pair pairing on the first image data and the second image data to obtain a first interference image pair and a second interference image pair; performing interference processing on the first interference image pair and the second interference image respectively to obtain a first processed image and a second processed image; performing track refinement and re-flattening on the first processed image and the second processed image respectively to obtain first track data and second track data after track refinement; Inverting the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result; After geocoding, the first SBAS result and the second SBAS result are converted into the image Map1 and the image Map2 in ArcGIS data format.

3. The method for identifying mining subsidence boundaries according to claim 2, characterized in that: The method of using SBAS to invert the first orbit data and the second orbit data to obtain a first SBAS result and a second SBAS result includes: Inverting the first orbit data and the second orbit data using the SBAS to estimate corresponding deformation rate and residual topography; Based on the corresponding deformation rate, the atmospheric filter is customized, the atmospheric phase is estimated and removed, and the displacement results on the corresponding time series are obtained as the first SBAS result and the second SBAS result.

4. The method for identifying mining subsidence boundaries according to claim 1, characterized in that: The N1=N2≥15.

5. A mining subsidence boundary identification device, characterized in that: The device comprises: An acquisition module is used to acquire N1 scene synthetic aperture radar image data before and after the mining start time as the first image data, and to acquire N2 scene synthetic aperture radar image data within a set time including the subsidence area as the second image data; A processing module, used to process the first image data and the second image data respectively using SARscape software to obtain an image Map1 and an image Map2 in ArcGIS data format; A visualization module is used to input the image Map1 and the image Map2 into ArcGIS to obtain 3D visualization maps Map 1-3D With Map 2-3D ; The first boundary determination module is used to map the 3D visualization map 1-3D With Map 2-3D The boundary of the public area is determined as the boundary of the settled area; The second boundary determination module is used to determine the boundary from the 3D visualization map 2-3D Get the boundary of the unstable area.

6. The mining subsidence boundary identification device according to claim 5, characterized in that: The processing module comprises: A pairing unit, used for pairing the first image data and the second image data with interference images to obtain a first interference image pair and a second interference image pair; An interference unit, used for performing interference processing on the first interference image pair and the second interference image respectively to obtain a first processed image and a second processed image; A track refining unit, used for performing track refinement and re-flattening on the first processed image and the second processed image respectively, to obtain first track data and second track data after track refinement; an inversion unit, configured to invert the first orbit data and the second orbit data using SBAS to obtain a first SBAS result and a second SBAS result; The encoding conversion unit is used for converting the first SBAS result and the second SBAS result into the image Map1 and the image Map2 in ArcGIS data format after geocoding.

7. The mining subsidence boundary identification device according to claim 6, characterized in that: The inversion unit comprises: an estimation subunit, configured to estimate corresponding deformation rate and residual topography by inverting the first orbit data and the second orbit data using the SBAS; The displacement calculation subunit is used to customize the atmospheric filtering based on the corresponding deformation rate, estimate and remove the atmospheric phase, and obtain the displacement results on the corresponding time series as the first SBAS result and the second SBAS result.

8. The mining subsidence boundary identification device according to claim 5, characterized in that: The N1=N2≥15.

9. A device, characterized in that: The device includes a memory and a processor; wherein the processor runs a program corresponding to the executable program code by reading the executable program code stored in the memory, so as to implement the method according to any one of claims 1 to 4.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 4 is implemented.

Images