CN110058234A - A method of it resolving mining area surface and settles three-dimensional deformation - Google Patents
A method of it resolving mining area surface and settles three-dimensional deformation Download PDFInfo
- Publication number
- CN110058234A CN110058234A CN201910418824.3A CN201910418824A CN110058234A CN 110058234 A CN110058234 A CN 110058234A CN 201910418824 A CN201910418824 A CN 201910418824A CN 110058234 A CN110058234 A CN 110058234A
- Authority
- CN
- China
- Prior art keywords
- deformation
- data
- width
- time
- sentinel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
Abstract
A kind of method that the present invention discloses resolving mining area surface sedimentation three-dimensional deformation, includes the following steps: the N+1 scape Sentinel-1A data of S1, the research time for obtaining survey region;S2, the N+1 scape Sentinel-1A data for studying the time for cutting survey region;System parameter and evaluation time and Space Baseline needed for S3, selection data processing;S4, extract survey region the research time sight to deformation map;S5, the orientation deformation map for studying the time for extracting survey region;S6, calculate area-of-interest sight to accumulation deformation quantity;S7, the orientation for calculating area-of-interest accumulate deformation quantity;S8, the three-dimensional shaped variable field for resolving each accumulation deformation quantity;A kind of data source is used only in the present invention, and not against any external data, the accurate location and its settling amount that earth's surface three-dimensional sedimentation information in survey region can be provided realize quick, efficient, cheap, the accurate quantificationization monitoring of a wide range of mining area surface three-dimensional sedimentation information based on radar satellite remotely-sensed data.
Description
Technical field
The present invention relates to resolving mining area surface to settle three-dimensional deformation field, radar satellite remote sensing is utilized more particularly to a kind of
Data are quick, efficient, cheap, a wide range of mining area subsidence three-dimensional deformation method of precise quantification.
Background technique
Pit mining is one of the main reason for causing mining area surface to settle, and the accurate mining area surface that extracts settles three-dimensional deformation
It is all of great significance to protecting and restoring periphery ecological environment, humanistic environment and geological environment.
Had at present using more method: based on the multi-platform or multi-track SAR number for covering the identical regional and same time
According to resolving three-dimensional deformation;Based on identical area and the ascending, descending rail SAR data of same time is covered, deviated in conjunction with DInSAR, pixel
One of amount (pixel offset tracking), multiple aperture InSAR (MAI), two or three of technology extract three-dimensional jointly
Deformation;InSAR (DInSAR, PS-InSAR and SBAS-InSAR) technology combination external GPS data, Levelling data or combination
Crustal deformation model resolves three-dimensional deformation information etc..
Since requirement multi-platform, that multi-track data are to data is harsher, lead to traditional method in practical applications
Limited promise;Since the precision of pixel offset tracking is lower than MAI technology, it is applied to monitoring earthquake, glacier mostly
The big magnitude deformation of movement, volcanic eruption etc., so being not particularly suited for the monitoring of mining area three-dimensional deformation;In addition, not all need
Wanting area to be monitored, all there is exactly drop rail data, to be unable to satisfy at least three kinds of imaging geometry data for resolving three-dimensional deformation
Condition at least obtains sight to the condition with orientation deformation by two orbital datas;Although being monitoring shape in conjunction with external data
Become the higher method of precision, but the manpower as needed for geodesic survey, financial resources, material resources expend huge leads to its spatial resolution
It is low, and select influence of the suitable crustal deformation model for calculation accuracy significant, it is not suitable for mining area in this way
Large-area three-dimensional deformation monitoring.Therefore there are data acquisition difficulties, the external priori of dependence for traditional three-dimensional deformation extracting method
The shortcomings that data, required cost are high and are not applied for generally existing a wide range of ground settlement.In current National mining area
Geological disaster takes place frequently and under the serious background of ecological environment destruction, needs to invent a kind of quick, efficient, cheap, accurate standard
Method carrys out Quantitative Monitoring ground settlement three-dimensional deformation information.
Summary of the invention
It is above-mentioned of the existing technology to solve the object of the present invention is to provide a kind of method that resolving mining area three dimensional type becomes
Problem.
To achieve the above object, the present invention provides following schemes:
S1, the N+1 scape Sentinel-1A data for studying the time for obtaining survey region
The IW Mode S entinel-1A SAR data for downloading the research time of survey region, decompresses compressed data,
And Sentinel-1A data are imported into SARscape software, by Sentinel-1A data conversion at standard data format;
S2, the N+1 scape Sentinel-1A data for studying the time for cutting survey region
On the intensity data figure of the Sentinel-1A normal data derived from step S1, chosen by vector boundary interested
Region cuts area-of-interest, to improve data processing speed;
System parameter and evaluation time and Space Baseline needed for S3, selection data processing
The Sentinel-1A data system parameter SENTINEL-TOPSAR for being suitable for IW mode is chosen, and estimates N number of shadow
The time of picture pair and Space Baseline, according to parallax range, time interval, seasonal variations and vegetation growth to image to coherence
Influence, choose the time reference line of image pair, Space Baseline;
S4, extract survey region the research time sight to deformation map
According to the area-of-interest SAR data that step S2 is cut out, the N width sight in research time is obtained to deformation map, is tied
4 data of SRTM-3Version are closed successively to carry out interference pattern generation to the SLC SAR image that polarization mode is VV, filtering, be concerned with
Property calculating, phase unwrapping, control point selection, track is concise and goes to put down again, after phase turns deformation and geocoding operation, output
Sight is to deformation map;
S5, the orientation deformation map for studying the time for extracting survey region
According to the area-of-interest SAR data that step S2 is cut out, the N width orientation deformation map in research time, knot are obtained
4 data of SRTM-3Version are closed successively to carry out interference pattern generation to the SLC SAR image that polarization mode is VV, filtering, be concerned with
Property calculating, phase unwrapping, control point selection, track is concise and goes to put down again, after phase turns deformation and geocoding operation, output
Orientation deformation map;
S6, calculate area-of-interest sight to accumulation deformation quantity
The 1st width sight obtained by step S4 is to based on deformation map, being sequentially overlapped the 2nd adjacent width, the 2nd, 3 width ...,
2nd, 3 ..., N-1, N width deformation map, get the N-1 width sight of opposite 1st width deformation map variation to accumulation deformation map;
S7, the orientation for calculating area-of-interest accumulate deformation quantity
Based on the 1st width orientation deformation map that above-mentioned steps S5 is obtained, it is sequentially overlapped the 2nd adjacent width, the 2nd, 3
Width ..., the 2nd, 3 ..., N-1, N width deformation map, get the N-1 width orientation accumulation deformation of opposite 1st width deformation map variation
Figure;
S8, the three-dimensional shaped variable field for resolving each accumulation deformation quantity
Sight is to deformation according to actual vertically to, East and West direction and north-south deformation vector superposition composition, orientation deformation
It is constituted according to actual East and West direction and the superposition of north-south deformation vector, based on sight to specific between deformation, orientation deformation
Functional relation can calculate the practical three-dimensional deformation of mining area surface sedimentation by adjustment means.
The invention discloses following technical effects: being difficult to obtain simultaneously with multi-platform, multi-track data in the prior art
It compares, a kind of SAR data, which is used only, in the present invention can calculate three-dimensional deformation;With pixel offset tracking technology
Monitoring accuracy is low to be compared, and the present invention is in the wrong using DInSAR and MAI technical office and pays a return visit period SAR data, and it is big can to reduce monitoring
There are the probability of large error for gradient deformation, and can obtain high-precision deformation results;With need combine external geodetic surveying data
It compares, the present invention has sufficiently excavated the application potential of radar remote sensing satellite data, can not only round-the-clock, round-the-clock to earth's surface
High precision monitor is carried out, high quality monitoring result can also be got in a wide range of interior do not kept off by cloud, canopy;With hull shape in combination
Varying model is compared, and the present invention is not bound by the influence for thinking inverted parameters, can more really reflect practical deformation situation.Therefore, with biography
The mining area surface sedimentation three-dimensional deformation calculation method of system is compared, and a kind of data source is used only in the present invention, and not against any outside
Data, so that it may which the accurate location and its settling amount for providing earth's surface three-dimensional sedimentation information in survey region are realized and defended based on radar
Quick, efficient, cheap, the accurate quantificationization monitoring of a wide range of mining area surface three-dimensional sedimentation information of star remotely-sensed data.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is flow diagram of the invention;
Fig. 2 is pavilion hilllock mining area three-dimensional deformation cylindricality schematic diagram of the invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
- 2 present invention provide a kind of method for resolving mining area surface three-dimensional deformation referring to Fig.1, include the following steps:
S1, the N+1 scape Sentinel-1A data for studying the time for obtaining survey region
Download Sentinel-1A of the survey region from totally 21 scape IW modes on June 6,9 days to 2018 October in 2017
SAR data decompresses compressed data, and Sentinel-1A data are imported SARscape software, by Sentinel-1A
Data conversion is at standard data format;The azimuth resolution and range resolution of downloaded Sentinel-1A image are distinguished
For 5 meters and 20 meters, after multiple look processing, azimuth resolution and range resolution are respectively 13.9018 meters and 17.674 meters;
S2, the N+1 scape Sentinel-1A data for studying the time for cutting survey region
On the intensity data figure of the Sentinel-1A normal data derived from step S1, by the vector boundary in pavilion hilllock mining area
Area-of-interest is chosen, area-of-interest is cut out, to improve data processing speed;
System parameter and evaluation time and Space Baseline needed for S3, selection data processing
The Sentinel-1A data system parameter SENTINEL-TOPSAR for being suitable for IW mode is chosen, and estimates N number of shadow
The time of picture pair and Space Baseline, according to parallax range, time interval, seasonal variations and vegetation growth to image to coherence
Influence, the time reference line of selected image pair is 12 days, and Space Baseline is less than 120 meters;
S4, extract survey region the research time sight to deformation map
According to the 21 scape SAR datas in the covering pavilion hilllock mining area cut out step S2, June 69 days to 2018 October in 2017
Day totally 20 width sights to deformation map, in conjunction with 4 data of SRTM-3Version to the SLC SAR image that polarization mode is VV successively into
The generation of row interference pattern, filtering, coherence calculation, phase unwrapping, control point selection, track is concise and goes to put down again, phase turns deformation
And after geocoding operation, sight is exported to deformation map;
S5, the orientation deformation map for studying the time for extracting survey region
According to the 21 scape SAR datas in the covering pavilion hilllock mining area cut out step S2,9 days to 2018 October in 2017 is obtained
June 6 totally 20 width orientation deformation map, in conjunction with 4 data of SRTM-3Version to the SLC SAR image that polarization mode is VV according to
Secondary progress interference pattern generation, filtering, coherence calculation, phase unwrapping, control point selection, track is concise and goes to put down again, phase turns
After deformation and geocoding operation, orientation deformation map is exported;
S6, calculate area-of-interest sight to accumulation deformation quantity
The 20 width sights obtained by step S4 to based on deformation map, respectively with the 1st width sight to deformation map be starting when
Between point, be sequentially overlapped the 2nd adjacent width, the 2nd, 3 width ..., the 2nd, 3 ..., 19,20 width deformation maps, get opposite 1st width shape
Change schemes the 20 width sights changed to accumulation deformation map;
S7, the orientation for calculating area-of-interest accumulate deformation quantity
Based on the 20 width orientation deformation maps that above-mentioned steps S5 is obtained, it is with the 1st width orientation deformation map respectively
Begin time point, be sequentially overlapped the 2nd adjacent width, the 2nd, 3 width ..., the 2nd, 3 ..., 19,20 width deformation maps, get the opposite 1st
20 width orientations of width deformation map variation accumulate deformation map;
S8, the three-dimensional shaped variable field for resolving each accumulation deformation quantity
According to the geometry imaging relations of radar, sight can be from actual vertically to, East and West direction and north-south shape to deformation
Become vector superposed composition, orientation deformation can be made of actual East and West direction and the superposition of north-south deformation vector, such as following formula institute
Show:
In formula 1, θ is radar wave incidence angle, and α is satellite course angle.DLOSAnd DAZORespectively acquired sight to
Orientation deformation;Dv, De, DnIt is respectively true vertically to East and West direction and north-south deformation vector;Write formula 1 as matrix form
Are as follows:
Based on formula 2, unique generalized inverse of coefficient matrix B is asked, and then the only optimal solution for calculating unknown three-dimensional deformation is public
Formula are as follows:
According to formula 3, the three-dimensional deformation information of mining area surface sedimentation can be calculated.
In the description of the present invention, it is to be understood that, term " longitudinal direction ", " transverse direction ", "upper", "lower", "front", "rear",
The orientation or positional relationship of the instructions such as "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outside" is based on attached drawing institute
The orientation or positional relationship shown is merely for convenience of the description present invention, rather than the device or element of indication or suggestion meaning must
There must be specific orientation, be constructed and operated in a specific orientation, therefore be not considered as limiting the invention.
Embodiment described above is only that preferred embodiment of the invention is described, and is not carried out to the scope of the present invention
It limits, without departing from the spirit of the design of the present invention, those of ordinary skill in the art make technical solution of the present invention
Various changes and improvements, should all fall into claims of the present invention determine protection scope in.
Claims (1)
1. a kind of method for resolving mining area surface three-dimensional deformation, characterized by the following steps:
S1, the N+1 scape Sentinel-1A data for studying the time for obtaining survey region
The IW Mode S entinel-1A SAR data for downloading the research time of survey region, decompresses compressed data, and will
Sentinel-1A data import SARscape software, by Sentinel-1A data conversion at standard data format;
S2, the N+1 scape Sentinel-1A data for studying the time for cutting survey region
On the intensity data figure of the Sentinel-1A normal data derived from step S1, area-of-interest is chosen by vector boundary,
Area-of-interest is cut, to improve data processing speed;
System parameter and evaluation time and Space Baseline needed for S3, selection data processing
The Sentinel-1A data system parameter SENTINEL-TOPSAR for being suitable for IW mode is chosen, and estimates N number of image pair
Time and Space Baseline, according to parallax range, time interval, seasonal variations and vegetation growth to image to the shadow of coherence
It rings, chooses time reference line, the Space Baseline of image pair;
S4, extract survey region the research time sight to deformation map
The area-of-interest SAR data cut out according to step S2, obtains the N width sight in research time to deformation map, in conjunction with
4 data of SRTM-3Version successively carry out interference pattern generation, filtering, coherence to the SLC SAR image that polarization mode is VV
Calculating, phase unwrapping, control point selection, track is concise and goes to put down again, after phase turns deformation and geocoding operation, output view
Line is to deformation map;
S5, the orientation deformation map for studying the time for extracting survey region
According to the area-of-interest SAR data that step S2 is cut out, the N width orientation deformation map in research time is obtained, in conjunction with
4 data of SRTM-3Version successively carry out interference pattern generation, filtering, coherence to the SLC SAR image that polarization mode is VV
Calculating, phase unwrapping, control point selection, track is concise and goes to put down again, after phase turns deformation and geocoding operation, output side
Position is to deformation map;
S6, calculate area-of-interest sight to accumulation deformation quantity
The 1st width sight obtained by step S4 is to based on deformation map, being sequentially overlapped the 2nd adjacent width, the 2nd, 3 width ..., the 2nd,
3 ..., N-1, N width deformation map gets the N-1 width sight of opposite 1st width deformation map variation to accumulation deformation map;
S7, the orientation for calculating area-of-interest accumulate deformation quantity
Based on the 1st width orientation deformation map that above-mentioned steps S5 is obtained, it is sequentially overlapped the 2nd adjacent width, the 2nd, 3 width ...,
2nd, 3 ..., N-1, N width deformation map, get the N-1 width orientation accumulation deformation map of opposite 1st width deformation map variation;
S8, the three-dimensional shaped variable field for resolving each accumulation deformation quantity
Sight to deformation according to actual vertically to, East and West direction and north-south deformation vector superposition composition, orientation deformation according to
Actual East and West direction and the superposition of north-south deformation vector are constituted, based on sight to the specific function between deformation, orientation deformation
Relational expression can calculate the practical three-dimensional deformation of mining area surface sedimentation by adjustment means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418824.3A CN110058234A (en) | 2019-05-20 | 2019-05-20 | A method of it resolving mining area surface and settles three-dimensional deformation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418824.3A CN110058234A (en) | 2019-05-20 | 2019-05-20 | A method of it resolving mining area surface and settles three-dimensional deformation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110058234A true CN110058234A (en) | 2019-07-26 |
Family
ID=67323715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910418824.3A Pending CN110058234A (en) | 2019-05-20 | 2019-05-20 | A method of it resolving mining area surface and settles three-dimensional deformation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110058234A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111076704A (en) * | 2019-12-23 | 2020-04-28 | 煤炭科学技术研究院有限公司 | Method for accurately calculating ground surface subsidence of coal mining subsidence area by using INSAR |
CN111562557A (en) * | 2020-05-12 | 2020-08-21 | 太原理工大学 | Remote sensing prediction method for ground settlement of loess plateau Xuangang mining area |
CN111650579A (en) * | 2020-06-12 | 2020-09-11 | 中南大学 | InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium |
CN111650587A (en) * | 2020-06-24 | 2020-09-11 | 中南大学 | Mining area ground surface three-dimensional dynamic deformation estimation method and device considering movement law and storage medium |
CN111859689A (en) * | 2020-07-27 | 2020-10-30 | 太原理工大学 | Method for determining loess plateau mining area ground settlement estimation model parameters |
CN111895903A (en) * | 2020-07-21 | 2020-11-06 | 太原理工大学 | Remote sensing estimation method for snow depth in northern area of Xinjiang |
CN112241712A (en) * | 2020-10-22 | 2021-01-19 | 山东省地质矿产勘查开发局第一地质大队 | Mineral resource acquisition and monitoring system |
CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
CN113900117A (en) * | 2021-09-13 | 2022-01-07 | 东华理工大学 | Underground non-evidence exploitation identification method integrating PS-InSAR and optical remote sensing |
CN114236541A (en) * | 2021-12-08 | 2022-03-25 | 电子科技大学 | Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090284758A1 (en) * | 2008-05-19 | 2009-11-19 | Toyonaka Kenkyusho Co., Ltd. | Displacement measuring method, displacement measuring apparatus and displacement measuring target |
CN106526590A (en) * | 2016-11-04 | 2017-03-22 | 山东科技大学 | Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image |
CN107102332A (en) * | 2017-05-11 | 2017-08-29 | 中南大学 | The three-dimensional earth's surface deformation monitoring methods of InSAR based on variance components estimate and strees strain model |
CN107271998A (en) * | 2017-07-01 | 2017-10-20 | 东华理工大学 | A kind of integrated D InSAR and the underground illegal mining recognition methods of GIS technology and system |
CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
CN109029344A (en) * | 2018-07-10 | 2018-12-18 | 湖南中科星图信息技术有限公司 | A kind of dykes and dams Monitoring method of the subsidence based on high score image and lift rail InSAR |
CN109471104A (en) * | 2018-11-15 | 2019-03-15 | 首都师范大学 | A method of obtaining earth's surface three-dimensional amount of movement from two parallel orbit SAR datas |
CN109696152A (en) * | 2019-02-13 | 2019-04-30 | 太原理工大学 | A kind of low coherence regional land subsidence amount evaluation method |
CN109738892A (en) * | 2019-01-24 | 2019-05-10 | 中南大学 | A kind of mining area surface high-spatial and temporal resolution three-dimensional deformation estimation method |
-
2019
- 2019-05-20 CN CN201910418824.3A patent/CN110058234A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090284758A1 (en) * | 2008-05-19 | 2009-11-19 | Toyonaka Kenkyusho Co., Ltd. | Displacement measuring method, displacement measuring apparatus and displacement measuring target |
CN106526590A (en) * | 2016-11-04 | 2017-03-22 | 山东科技大学 | Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image |
CN107102332A (en) * | 2017-05-11 | 2017-08-29 | 中南大学 | The three-dimensional earth's surface deformation monitoring methods of InSAR based on variance components estimate and strees strain model |
CN107271998A (en) * | 2017-07-01 | 2017-10-20 | 东华理工大学 | A kind of integrated D InSAR and the underground illegal mining recognition methods of GIS technology and system |
CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
CN109029344A (en) * | 2018-07-10 | 2018-12-18 | 湖南中科星图信息技术有限公司 | A kind of dykes and dams Monitoring method of the subsidence based on high score image and lift rail InSAR |
CN109471104A (en) * | 2018-11-15 | 2019-03-15 | 首都师范大学 | A method of obtaining earth's surface three-dimensional amount of movement from two parallel orbit SAR datas |
CN109738892A (en) * | 2019-01-24 | 2019-05-10 | 中南大学 | A kind of mining area surface high-spatial and temporal resolution three-dimensional deformation estimation method |
CN109696152A (en) * | 2019-02-13 | 2019-04-30 | 太原理工大学 | A kind of low coherence regional land subsidence amount evaluation method |
Non-Patent Citations (2)
Title |
---|
BEI ZHANG ET AL.: "Using Phase Unwrapping Methods to Apply D-InSAR in Mining Areas", 《CANADIAN JOURNAL OF REMOTE SENSING》 * |
尹宏杰 等: "基于SBAS的矿区形变监测研究", 《测绘学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111076704B (en) * | 2019-12-23 | 2022-05-20 | 煤炭科学技术研究院有限公司 | Method for accurately resolving ground surface subsidence of coal mining subsidence area by using INSAR |
CN111076704A (en) * | 2019-12-23 | 2020-04-28 | 煤炭科学技术研究院有限公司 | Method for accurately calculating ground surface subsidence of coal mining subsidence area by using INSAR |
CN111562557A (en) * | 2020-05-12 | 2020-08-21 | 太原理工大学 | Remote sensing prediction method for ground settlement of loess plateau Xuangang mining area |
CN111650579A (en) * | 2020-06-12 | 2020-09-11 | 中南大学 | InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium |
CN111650579B (en) * | 2020-06-12 | 2022-09-30 | 中南大学 | InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium |
CN111650587A (en) * | 2020-06-24 | 2020-09-11 | 中南大学 | Mining area ground surface three-dimensional dynamic deformation estimation method and device considering movement law and storage medium |
CN111895903A (en) * | 2020-07-21 | 2020-11-06 | 太原理工大学 | Remote sensing estimation method for snow depth in northern area of Xinjiang |
CN111895903B (en) * | 2020-07-21 | 2021-06-01 | 太原理工大学 | Remote sensing estimation method for snow depth of detection area |
CN111859689A (en) * | 2020-07-27 | 2020-10-30 | 太原理工大学 | Method for determining loess plateau mining area ground settlement estimation model parameters |
CN111859689B (en) * | 2020-07-27 | 2023-06-16 | 太原理工大学 | Determination method for loess plateau mining area ground subsidence estimation model parameters |
CN112241712A (en) * | 2020-10-22 | 2021-01-19 | 山东省地质矿产勘查开发局第一地质大队 | Mineral resource acquisition and monitoring system |
CN112241712B (en) * | 2020-10-22 | 2023-04-07 | 山东省地质矿产勘查开发局第一地质大队 | Mineral resource acquisition and monitoring system |
CN113091598B (en) * | 2021-04-06 | 2022-02-08 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
CN113900117A (en) * | 2021-09-13 | 2022-01-07 | 东华理工大学 | Underground non-evidence exploitation identification method integrating PS-InSAR and optical remote sensing |
CN114236541A (en) * | 2021-12-08 | 2022-03-25 | 电子科技大学 | Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image |
CN114236541B (en) * | 2021-12-08 | 2023-05-16 | 电子科技大学 | Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110058234A (en) | A method of it resolving mining area surface and settles three-dimensional deformation | |
Strozzi et al. | JERS SAR interferometry for land subsidence monitoring | |
Delacourt et al. | Velocity field of the “La Clapière” landslide measured by the correlation of aerial and QuickBird satellite images | |
CN102927934B (en) | Method for obtaining mining area earth surface three-dimensional deformation fields through single interferometric synthetic aperture radar (InSAR) interference pair | |
Long et al. | Acquisition and analysis of terrestrial gravity data | |
Guo et al. | Land subsidence in Tianjin for 2015 to 2016 revealed by the analysis of Sentinel-1A with SBAS-InSAR | |
Haase et al. | Coordinates and maps of the Apollo 17 landing site | |
Wolf | Surveying and mapping: History, current status, and future projections | |
Sam et al. | Estimation of dune celerity and sand flux in part of West Rajasthan, Gadra area of the Thar Desert using temporal remote sensing data | |
Fueten et al. | Stratigraphy and structure of interior layered deposits in west Candor Chasma, Mars, from High Resolution Stereo Camera (HRSC) stereo imagery and derived elevations | |
Leberl | Imaging radar applications to mapping and charting | |
Xie et al. | Determination of the height of Mount Everest using the shallow layer method | |
Gerstbach et al. | A small CCD zenith camera (ZC-G1)-developed for rapid geoid monitoring in difficult projects | |
Mondal | Mapping India since 1767: transformation from colonial to postcolonial image | |
Hargitai et al. | Methods in planetary topographic mapping: a Review | |
Li | Coastline mapping and change detection using one-meter resolution satellite imagery | |
Xu | Three-dimensional virtual geology: photorealistic outcrops, and their acquisition, visualization and analysis | |
Ata | A test of the validity of morphometric analysis in determining tectonic activity from aster derived DEMS in the Jordan-dead sea transform zone | |
Ping | Digital elevation model extraction from ASTER in support of the coal fire and environmental research project, China | |
Karachevtseva et al. | Cartography of the Soviet Lunokhods’ Routes on the Moon | |
Alberti | GIS analysis of geological surfaces orientations: the qgSurf plugin for QGIS | |
Al-Siaede | Tectonic Interpretation of the Relation between Rivers Courses and The Subsurface Structures in Central Maysan Petroleum Province SE Iraq | |
Riedel | Airborne-based geophysical investigation in Dronning Maud Land, Antarctica | |
Dickerson et al. | Geological methods | |
PICU | CONTRIBUTIONS TO THE FURTHER DEVELOPMENT OF SMART CITIES |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190726 |
|
RJ01 | Rejection of invention patent application after publication |