CN108153979A - Deformation information extraction method based on InSAR, terminal and storage medium - Google Patents
Deformation information extraction method based on InSAR, terminal and storage medium Download PDFInfo
- Publication number
- CN108153979A CN108153979A CN201711430532.9A CN201711430532A CN108153979A CN 108153979 A CN108153979 A CN 108153979A CN 201711430532 A CN201711430532 A CN 201711430532A CN 108153979 A CN108153979 A CN 108153979A
- Authority
- CN
- China
- Prior art keywords
- points
- initial
- deformation
- point
- insar
- 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.)
- Granted
Links
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a deformation information extraction method based on InSAR, a terminal and a storage medium, wherein the method comprises the following steps: A. extracting a plurality of initial PS points of a target building; B. performing clustering analysis processing on the plurality of initial PS points to obtain a plurality of aggregation points, and recording the number of the initial PS points corresponding to each aggregation point; C. calculating the distance between two adjacent aggregation points, and calculating the quantity difference of the initial PS points corresponding to the two adjacent aggregation points; judging whether two adjacent aggregation points meeting a preset singular aggregation point judgment condition exist or not, wherein the preset singular aggregation point judgment condition is that the distance is smaller than a preset distance threshold value, and the quantity difference value is larger than a preset quantity threshold value; D. and if so, rejecting aggregation points with less initial PS points in two adjacent aggregation points, and repeating the step C until no aggregation points exist. E. And if not, extracting deformation information according to the remaining aggregation points. The invention not only improves the processing rate, but also improves the analysis accuracy.
Description
Technical field
The present invention relates to synthetic aperture radar and big data analysis technical field more particularly to a kind of shape based on InSAR
Become information extracting method, terminal and storage medium.
Background technology
Surface subsidence is also known as land subsidence or earth subsidence, it is under the influence of economic activities of human beings, since underground is loose
Formation consolidation compresses, a kind of local descending motion (or engineeriing geological phenomenon) that earth's crust surface level is caused to reduce.From result
On see, surface subsidence influences the foundation structure of building, and then reduces its security performance.Especially differential settlement, it may
The local dip of building is caused, changes the stress state of its foundation structure, the exterior wall of building is likely to occur crack when serious
And crackle, if not taking reinforcement measure, with the further accumulation of settling amount, in some instances it may even be possible to lead to the safe thing of building
Therefore.
According to technological means of today, it is hardly possible to avoid Earth Subsidence Hazards completely.Therefore, the management in city
Persons are wished through applied science technology, the risk of accurate evaluation city owned building, and then it is effectively integrated
Management reduces loss caused by Earth Subsidence Hazards as far as possible.
Spaceborne InSAR (InSAR) technology belongs to Global observation by remote sensing, and there are other technologies
Not available advantage.It is not illuminated by the light the limitation with weather condition, and earth's surface letter can be obtained under round-the-clock, all weather conditions
Breath.By taking Italian COSMO-SkyMed systems as an example, the coverage area observed every time is up to thousands of square kilometres, spatial resolution
Highest can reach 1 meter, and distortion measurement precision is up to mm magnitudes.It is defended simultaneously as the system includes four SAR with rail operation
Star, the temporal resolution of distortion measurement reach as high as 4 times/16 days.InSAR remote sensing technologies have the ability to supervise for urban ground subsidence
It surveys and building risk assessment application provides data and supports.
By the development of twenties years, InSAR distortion measurement technologies were widely used in engineering practice.For example:
InSAR distortion measurements technology is applied to the extraction of building deformation data and the deformation data according to extraction is analyzed and processed.
But the building deformation data for being directly based upon InSAR extractions, which carries out analyzing and processing, there is problems:
For building target, if some a part or whole part deforms upon, then it is hidden that this solitary building has safety
The possibility of trouble is larger.But if only small range deforms upon, the possibility to damage to building is little.Therefore,
If InSAR distortion measurement technologies extract the small-scale deformation data, and the small-scale deformation data is carried out at analysis
Reason, then can reduce precision of analysis.
Invention content
The purpose of the present invention is to provide a kind of deformation data extracting method, terminal and storage mediums based on InSAR, pick
It is analyzed except the deformation data for not influencing safe condition, and then based on the deformation data after rejecting, it is accurate so as to improve analysis
True rate.
To solve the above-mentioned problems, the present invention provides a kind of deformation data extracting method based on InSAR, including:
A, multiple initial p S points of target construction are extracted;
B, cluster analysis processing is carried out to multiple initial p S points, obtains multiple congruent points, record each congruent point correspondence
Initial p S points quantity;
C, the distance between two adjacent polymeric points are calculated, and calculate the number of the corresponding initial p S points of two adjacent polymeric points
Measure difference;It judges whether to meet two adjacent polymeric points for presetting unusual congruent point decision condition, presets unusual congruent point
Decision condition is that distance is less than pre-determined distance threshold value, and number differences are more than preset quantity threshold value;
If D, in the presence of, reject two adjacent polymeric points in correspond to the few congruent point of initial p S points, repeat step C,
Until there is no two adjacent polymeric points for meeting default unusual congruent point decision condition;
If E, being not present, deformation data is extracted according to remaining congruent point.
As a further improvement on the present invention, step A, including:
Under preset coordinate system, multiple PS points in target zone are acquired;
Obtain space range information corresponding with target construction and PS point three-dimensional space position information;
Spatial data analysis treatment is carried out to multiple PS points according to space range information and PS point three-dimensional space position information,
To extract multiple initial p S points corresponding with target construction.
As a further improvement on the present invention, step B, including:
B1, two adjacent initial p S points form a PS point pair, and analyse whether there is the PS for meeting default fusion conditions
Point pair;
If B2, in the presence of to meeting the PS points of default fusion conditions to carrying out fusion treatment, to generate new PS points, newly
PS points correspond to two initial p S points, step B1 is repeated, until there is no meet the PS points pair of default fusion conditions;
If B3, being not present, the corresponding all initial p S points of each remaining PS point are obtained, and to all initial p S points
Polymerization processing is carried out, obtains a congruent point.
As a further improvement on the present invention, analyse whether there is the PS points pair for meeting default fusion conditions in step B1,
Including:
The relevant parameter of PS points pair is calculated, relevant parameter includes space length and deformation deviation standard deviation;
Judge whether that space length is less than pre-set space distance threshold, and deformation deviation standard deviation is less than preset standard
The PS points pair of poor threshold value;
To be less than pre-set space distance threshold, and deformation deviation standard deviation is less than preset standard difference threshold value when there are space lengths
PS points pair when, then PS points are to meeting default fusion conditions.
As a further improvement on the present invention, to meeting the PS points of default fusion conditions to carrying out at fusion in step B2
Reason, to generate new PS points, including:
It obtains with PS points to the location information and deformation history curve of corresponding two initial p S points;
Mean place information is calculated according to the location information of two initial p S points and deformation history curve and averaged deformation is gone through
History curve;
New PS points are generated according to mean place information and averaged deformation history curve.
As a further improvement on the present invention, in step B3 all initial p S points are carried out with polymerization processing, one is obtained and gathers
Chalaza, including:
Obtain the location information and deformation history curve of all initial p S points;
Mean place information is calculated according to the location information of all initial p S points and deformation history curve and averaged deformation is gone through
History curve;
According to mean place information and averaged deformation history curve generation congruent point.
As a further improvement on the present invention, step E includes later:
Output shows congruent point, and different congruent points corresponds to different output modes, and output mode includes color and/or shape
Shape.
As a further improvement on the present invention, step E includes later:
The congruent point without time point is obtained, and linear fit processing is carried out to the congruent point of different time points, to obtain
Deformation history curve.
To solve the above-mentioned problems, the present invention also provides a kind of terminal, including processor, memory and display,
Processor couples memory, display, and the computer program that can be run on a processor is stored on memory;
When processor performs computer program, the above-mentioned deformation data extracting method based on InSAR is realized.
To solve the above-mentioned problems, the present invention also provides a kind of storage mediums, are stored thereon with computer program, calculate
When machine program is executed by processor, the above-mentioned deformation data extracting method based on InSAR is realized.
The present invention is handled by cluster analysis, and multiple similar PS points are polymerized to a congruent point, and will be in congruent point
Unusual congruent point is (i.e.:Deformation is significantly unusual) it rejects, to obtain the congruent point for including crucial deformation data, in this case, one
Aspect is handled by cluster analysis, reduces data processing quantity, so as to improve processing speed.On the other hand, based on rejecting
Congruent point after unusual congruent point is analyzed and processed, and improves analysis accuracy rate.
Description of the drawings
Fig. 1 is the flow diagram of deformation data extracting method one embodiment the present invention is based on InSAR;
Fig. 2 is that the present invention is based on the flows that PS points in the deformation data extracting method of InSAR extract flow one embodiment
Schematic diagram;
Fig. 3 is the flow the present invention is based on cluster analysis flow one embodiment in the deformation data extracting method of InSAR
Schematic diagram;
Fig. 4, which is that the present invention is based on fusion conditions in the deformation data extracting method of InSAR, to judge flow one embodiment
Flow diagram;
Fig. 5 is that the present invention is based on the flow signals that flow one embodiment is merged in the deformation data extracting method of InSAR
Figure;
Fig. 6 is that the present invention is based on the flow signals of polymerization process one embodiment in the deformation data extracting method of InSAR
Figure;
Fig. 7 is the flow diagram of second embodiment of deformation data extracting method the present invention is based on InSAR;
Fig. 8 is the showing interface schematic diagram of terminal one embodiment of the present invention;
Fig. 9 is the flow diagram of the deformation data extracting method third embodiment the present invention is based on InSAR;
Figure 10 is the showing interface schematic diagram of second embodiment of terminal of the present invention;
Figure 11 is the showing interface schematic diagram of terminal third embodiment of the present invention;
Figure 12 is the showing interface schematic diagram of 4th embodiment of terminal of the present invention;
Figure 13 is the circuit theory schematic diagram of terminal one embodiment of the present invention;
Figure 14 is the high-level schematic functional block diagram of terminal one embodiment of the present invention;
Figure 15 is the high-level schematic functional block diagram of extraction module one embodiment of terminal of the present invention;
Figure 16 is the high-level schematic functional block diagram of cluster analysis processing module one embodiment of terminal of the present invention;
Figure 17 is the high-level schematic functional block diagram of same alanysis submodule one embodiment of terminal of the present invention;
Figure 18 is the high-level schematic functional block diagram of fusion submodule one embodiment of terminal of the present invention;
Figure 19 is the high-level schematic functional block diagram of polymerization submodule one embodiment of terminal of the present invention;
Figure 20 is the high-level schematic functional block diagram of second embodiment of terminal of the present invention;
Figure 21 is the high-level schematic functional block diagram of terminal third embodiment of the present invention;.
Specific embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, it is right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used to limit the present invention.
Fig. 1 illustrates one embodiment of the deformation data extracting method the present invention is based on InSAR.In the present embodiment,
The deformation data extracting method based on InSAR includes:
Step S1 extracts multiple initial p S points of target construction.
It should be noted that the target point of high correlation can be kept in selection Long time scale in InSAR images, and will
It is defined as PS (Permament Scatterers, permanent to scatter) point.
On the basis of above-described embodiment, in other embodiment, referring to Fig. 2, step S1, including:
Step S10 under preset coordinate system, acquires multiple PS points in target zone.
It should be noted that the preset coordinate system in the present embodiment includes space polar coordinate system, WGS-84 coordinate systems
The spatial coordinate systems such as system.Preferably WGS-84 coordinate systems.
In step slo, which can be a region, not only have the data of building in the region, also have
The data of other targets such as road, bridge.Therefore, there is the PS points for belonging to non-building in collected PS points.
Step S11 obtains space range information corresponding with target construction and PS point three-dimensional space position information.
In step s 11, each building has corresponding space range information and PS point three-dimensional space positions
Information, and in advance store the space range information and PS point three-dimensional space position information.Specifically, the target model is searched for
Interior target construction is enclosed, and obtains space range information corresponding with the target construction and PS points three-dimensional space position letter
Breath.
Step S12 carries out spatial data according to space range information and PS point three-dimensional space position information to multiple PS points
Analyzing and processing, to extract multiple initial p S points corresponding with target construction.
The present embodiment is automatically by the corresponding space range information of target construction and PS point three-dimensional space position information pair
Collected PS points carry out Screening Treatment, have both avoided the influence of the PS points of non-building, and have decreased PS point quantity, thus into
One step scapegoat data processing rate.
Step S2 carries out cluster analysis processing to multiple initial p S points, obtains multiple congruent points, record each congruent point
The quantity of corresponding initial p S points.
On the basis of above-described embodiment, in other embodiment, referring to Fig. 3, step S2 includes:
Step S20, two adjacent initial p S points form a PS point pair, and analyse whether to exist and meet default fusion conditions
PS points pair.If there is the PS points pair for meeting default fusion conditions, step S21 is performed.If there is no meet default fusion item
The PS points pair of part, then perform step S22.
It should be noted that this step mainly carries out similar PS points operation.Two adjacent initial p S in the present embodiment
In adjacent definition can be:Any one initial p S point B in the range of the pre-set space of one initial p S points A, the two are formed
Two adjacent initial p S points.
On the basis of above-described embodiment, in other embodiment, referring to Fig. 4, analyse whether to exist satisfaction in step S20 in advance
If the PS points of fusion conditions are to including:
Step S200, calculates the relevant parameter of PS points pair, and relevant parameter includes space length and deformation deviation standard deviation.
In step s 200, it is assumed that it is PS point A and PS points B to form the PS points pair.Wherein, the coordinate information of A is A (xA,
yA, zA), the coordinate information of B is B (xB, yB, zB)。
Specifically, space length DAB, wherein, the calculation formula of space length is as follows:
Further, deformation deviation average value is MAB, wherein, the calculation formula of deformation deviation average value is as follows:
Wherein, N represents time sampling quantity, and i represents time sampling ordinal number, dA, iWith
dB, iDeformation quantities of the A and B at the i-th moment is represented respectively.
Further, deformation deviation standard deviation is SAB, wherein, the calculation formula of deformation deviation standard deviation is as follows:
Wherein, N represents time sampling quantity, and i represents time sampling ordinal number, dA, iAnd dB, iRepresent A and B at i-th respectively
The deformation quantity at quarter.
Step S201 judges whether that space length is less than pre-set space distance threshold, and deformation deviation standard deviation is small
In the PS points pair of preset standard difference threshold value.If there are space lengths to be less than pre-set space distance threshold, and deformation deviation standard deviation
Less than the PS points pair of the poor threshold value of preset standard, then step S202 is performed.
It should be noted that pre-set space distance threshold and preset standard difference threshold value in the present embodiment, it both can be advance
Setting, can also carry out sets itself according to the demand of user.Preferably, pre-set space distance threshold may be configured as 20m, preset
Standard deviation threshold method could be provided as 2mm.
Step S202, PS point are to meeting default fusion conditions.
Space length is less than pre-set space distance threshold by the present embodiment, and deformation deviation standard deviation is poor less than preset standard
The PS points pair of threshold value carry out fusion treatment, further reduce the processing quantity of PS points, so as to further improve processing speed
Rate.
Step S21 presets the PS points of fusion conditions to meeting to carrying out fusion treatment, to generate new PS points, new PS
The corresponding two initial p S points of point, repeat step S20, until there is no the PS points pair for meeting default fusion conditions.
On the basis of above-described embodiment, in other embodiment, referring to Fig. 5, to meeting default fusion conditions in step S21
PS points to carry out fusion treatment, included with the PS points for generating new:
Step S210 is obtained with PS points to the location information and deformation history curve of corresponding two initial p S points.
In step S210, it is assumed that for PS points to corresponding PS point A and PS point B, then the coordinate information for obtaining A is A (xA, yA,
zA), the coordinate information of B is B (xB, yB, zB).Further, it is assumed that the previous moment of PS points A is A0, and the latter moment is
A1 then connects A0, and the curve of A, A1 are inflection curves A.Similarly, the previous moment of PS points B is B0, and the latter moment is B1,
B0 is then connected, the curve of B, B1 are inflection curves B.
Step S211, according to the location information and deformation history curve calculating mean place information peace of two initial p S points
Equal deformation history curve.
Step S212 generates new PS points according to mean place information and averaged deformation history curve.
In step S212, it is assumed that new PS points are C, and the coordinate of C is C (xc, yc, zc)。
ThenSimilarly, averaged deformation history curve is also deformation
Curve A and inflection curves B synthesize to obtain.
Step S22 obtains the corresponding all initial p S points of each remaining PS point, and all initial p S points is gathered
Conjunction is handled, and obtains a congruent point.
On the basis of above-described embodiment, in other embodiment, referring to Fig. 6, all initial p S points are carried out in step S22
Polymerization is handled, and is obtained a congruent point and is included:
Step S220 obtains the location information and deformation history curve of all initial p S points.
In step S220, it is assumed that formd after PS point C, PS point C and PS points D fusions are formed after PS point A and PS points B fusions
PS points E.Finally, PS points E is remaining PS points, then the corresponding initial p S points of PS points E are:PS point A, PS point B and PS points D.Assuming that
The previous moment of PS points A is A0, and the latter moment is A1, then connects A0, the curve of A, A1 are inflection curves A.Similarly, PS points
The previous moment of B is B0, and the latter moment is B1, then connects B0, the curve of B, B1 are inflection curves B.PS points D's is previous
Moment is D0, and the latter moment is D1, then connects D0, the curve of B, D1 are inflection curves D.
Step S221 calculates mean place information peace according to the location information of all initial p S points and deformation history curve
Equal deformation history curve.
Step S222 generates congruent point according to mean place information and averaged deformation history curve.
Assuming that congruent point is H, the coordinate of H is H (xh, yh, zh)。
ThenSimilarly, averaged deformation history
Curve is also that inflection curves A, inflection curves B and inflection curves D synthesize to obtain.
The present embodiment is handled by cluster analysis, and multiple similar PS points are polymerized to a congruent point, reduce PS points
Processing quantity, so as to improve processing speed.Further, it is subsequently analyzed and processed, is further reduced based on congruent point
Data processing quantity, so as to further improve data processing rate.
Step S3 calculates the distance between two adjacent polymeric points, and calculates the corresponding initial p S of two adjacent polymeric points
The number differences of point;It judges whether to meet two adjacent polymeric points for presetting unusual congruent point decision condition, preset unusual
Congruent point decision condition is that distance is less than pre-determined distance threshold value, and number differences are more than preset quantity threshold value.Meet in advance if existing
If two adjacent polymeric points of unusual congruent point decision condition, then perform step S4.If it does not preset unusual congruent point in satisfaction to sentence
Two adjacent polymeric points of fixed condition, then perform step S5.
It should be noted that this step mainly carries out unusual congruent point analysis operation.Two in the present embodiment are adjacent
Adjacent definition can be in congruent point:Any one congruent point B in the range of the pre-set space of one congruent point A, the two structure
Into two adjacent polymeric points.
In step s3, it is assumed that there are meeting point F and congruent point G, wherein, the coordinate information of F is F (xF, yF, zF), G's
Coordinate information is G (xG, yG, zG), the quantity of the corresponding initial p S points of F is a, and the quantity of the corresponding initial p S points of G is b, then quantity
Difference c=| a-b |.
Specifically, distance is DFG, wherein, the calculation formula of space length is as follows:
It should be noted that pre-determined distance threshold value and preset quantity threshold value in the present embodiment, can both preset,
Sets itself can be carried out according to the demand of user.
Step S4 rejects in two adjacent polymeric points and corresponds to the few congruent point of initial p S points, repeats step S3, until
There is no two adjacent polymeric points for meeting default unusual congruent point decision condition.
In step s3, in distance DFGLess than pre-determined distance threshold value, and number differences c>During preset quantity threshold value, if a >
B then rejects congruent point G.
Step S5 then extracts deformation data according to remaining congruent point.
The present embodiment is handled by cluster analysis, and multiple similar PS points are polymerized to a congruent point, and by congruent point
In unusual congruent point (i.e.:Deformation is significantly unusual) it rejects, to obtain the congruent point for including crucial deformation data, in this case,
On the one hand, it is handled by cluster analysis, reduces data processing quantity, so as to improve processing speed.On the other hand, based on picking
Except the congruent point after unusual congruent point is analyzed and processed, analysis accuracy rate is improved.
The deformation data extracting method based on InSAR of the present invention is applied to during the use of terminal, needs difference
The congruent point of classification uses different output modes, therefore, on the basis of above-described embodiment, in other embodiment, referring to Fig. 7,
After step S5, further include:
Step S30, output show congruent point, and different congruent points corresponds to different output modes, and output mode includes face
Color and/or shape.
Specifically, referring to Fig. 8, different shape represents different types of PS point sets, and larger point represents such PS point set
Close corresponding congruent point.
The present embodiment represents different types of PS point sets by different shape and/or different colours, and using size into
Row congruent point and PS initial ps S points intuitively, quickly know the cluster analysis result of PS points convenient for user, so as to further improve
User experience.
The deformation data extracting method based on InSAR of the present invention is applied to during the use of terminal, needs to know
Therefore the deformation history of PS points, on the basis of above-described embodiment, in other embodiment, referring to Fig. 9, after step S5, is also wrapped
It includes:
Step S40 obtains the congruent point without time point, and the congruent point of different time points is carried out at linear fit
Reason, to obtain deformation history curve.
Measured data parameter information used in this example is as shown:
Satellite type:COSMO-SkyMed
Screening-mode:StripMap
Spatial resolution:3m*3m
Polarization mode:HH
Center incident angle:29.376°
Image quantity:44
Monitor initial time:2012.12.28
Monitor the end time:2016.08.16
Specifically, referring to Figure 10, the deformation history curve of the PS points in the similar PS point groups in exhibition building object north side, Yi Jiju
The deformation history curve of chalaza.Wherein, filament is the deformation history curve of PS points, and thick line is the deformation history curve of congruent point.
Referring to Figure 11, the deformation history curves of the PS points in the similar PS point groups in exhibition building object southern side and congruent point
Deformation history curve.Wherein, filament is the deformation history curve of PS points, and thick line is the deformation history curve of congruent point.
Referring to Figure 12, the deformation history curves of the PS points in the middle part of exhibition building object in similar PS point groups and congruent point
Deformation history curve.Wherein, filament is the deformation history curve of PS points, and thick line is the deformation history curve of congruent point.
The present embodiment by deformation history curve, can in order to user quick, the intuitive deformation situation that know building, from
And further improve user experience.
Figure 13 be terminal device of the present invention one embodiment, in the present embodiment, the terminal device including memory 2,
Processor 1 and display 3, processor 1 are respectively coupled to memory 2, display 3.
Wherein, processor 1 can be a general central processor, microprocessor, application-specific integrated circuit or one
Or multiple integrated circuits for application scheme program to be controlled to perform.
In addition, memory 2 can be read-only memory, it is the static storage device that static information and instruction can be stored, random
Access memory or can store information and instruction dynamic memory or Electrically Erasable Programmable Read-Only Memory,
CD-ROM or other optical disc storages, optical disc storage, magnetic disk storage medium or other magnetic storage apparatus.Memory 2 and processing
Device 1 can be connected by communication bus, can also be integrated with processor 1.
Above-mentioned memory 2 can be used for storage to perform the computer program of application scheme, and processor 1 is available for performing storage
The computer program stored in device 2, to realize the deformation data extracting method based on InSAR of above-described embodiment description.
Specifically, referring to Figure 14, in the present embodiment, which includes initial information extraction module 10, at cluster analysis
It manages module 11, singularity analysis module 12, reject processing module 13 and deformation data extraction module 14.
Wherein, initial information extraction module 10, for extracting multiple initial p S points of target construction;Cluster analysis is handled
Module 11 for carrying out cluster analysis processing to multiple initial p S points, obtains multiple congruent points, records each congruent point correspondence
Initial p S points quantity;Singularity analysis module 12 for calculating the distance between two adjacent polymeric points, and calculates two
The number differences of the corresponding initial p S points of adjacent polymeric point;It judges whether to meet and presets the two of unusual congruent point decision condition
A adjacent polymeric point, it is that distance is less than pre-determined distance threshold value, and number differences are more than default to preset unusual congruent point decision condition
Amount threshold;Processing module 13 is rejected, if meeting two adjacent polymeric points for presetting unusual congruent point decision condition for existing,
It then rejects in two adjacent polymeric points and corresponds to the few congruent point of initial p S points, repeat singularity analysis module 12, until not
In the presence of two adjacent polymeric points for meeting default unusual congruent point decision condition.Deformation data extraction module 14, if for not depositing
In two adjacent polymeric points for meeting default unusual congruent point decision condition, then deformation data is extracted according to remaining congruent point.
On the basis of above-described embodiment, in other embodiment, referring to Figure 15, initial information extraction module 10 includes acquisition
Unit 100, building related information acquiring unit 101 and extraction unit 102.
As a further improvement on the present invention, multiple initial p S points of target construction are extracted, including:
Wherein, collecting unit 100, under preset coordinate system, acquiring multiple PS points in target zone;Building
Related information acquiring unit 101, for obtaining space range information corresponding with target construction and PS point three-dimensional space positions
Information;Extraction unit 102, for carrying out space to multiple PS points according to space range information and PS point three-dimensional space position information
Data Analysis Services, to extract multiple initial p S points corresponding with target construction.
On the basis of above-described embodiment, in other embodiment, referring to Figure 16, cluster analysis processing module 11 includes similar
Analyze submodule 110, fusion submodule 111 and polymerization submodule 112.
Wherein, with alanysis submodule 110, a PS point pair is formed, and analyse whether for two adjacent initial p S points
In the presence of the PS points pair for meeting default fusion conditions;Submodule 111 is merged, if for there is the PS points for meeting default fusion conditions
It is right, then the PS points of fusion conditions are preset to meeting to carrying out fusion treatment, to generate new PS points, new PS points are corresponded at the beginning of two
Beginning PS points repeat same alanysis submodule 110, until there is no the PS points pair for meeting default fusion conditions;It polymerize submodule
For if there is no the PS points pair for meeting default fusion conditions, it is corresponding all first to obtain each remaining PS point for block 112
Beginning PS points, and all initial p S points are carried out with polymerization processing, obtain a congruent point.
On the basis of above-described embodiment, in other embodiment, referring to Figure 17, this includes meter with alanysis submodule 110
Calculate unit 1101, judging unit 1102 and confirmation unit 1103.
Wherein, computing unit 1101, for calculating the relevant parameter of PS points pair, relevant parameter includes space length and deformation
Deviation standard is poor;Judging unit 1102, for judging whether that space length is less than pre-set space distance threshold, and deformation is inclined
Poor standard deviation is less than the PS points pair of preset standard difference threshold value;Confirmation unit 1103, for pre- when being less than there are the space length
If word space threshold, and the deformation deviation standard deviation be less than preset standard difference threshold value PS points pair when, then PS points to meet
Default fusion conditions.
On the basis of above-described embodiment, in other embodiment, referring to Figure 18, which obtains including first
Take unit 1110, the first average calculation unit 1111 and new PS points generation unit 1112.
Wherein, first acquisition unit 1110 obtain and location information of the PS points to corresponding two initial p S points and deformation
History curve;First average calculation unit 1111, for the location information and deformation history curve according to two initial p S points
Calculate mean place information and averaged deformation history curve;New PS points generation unit 1112, for according to mean place information
New PS points are generated with averaged deformation history curve.
On the basis of above-described embodiment, in other embodiment, referring to Figure 19, which obtains including second
Take unit 1120, the second average calculation unit 1121 and congruent point generation unit 1122.
Wherein, second acquisition unit 1120, for obtaining the location information and deformation history curve of all initial p S points;The
Two average calculation units 1121 calculate mean place for the location information according to all initial p S points and deformation history curve
Information and averaged deformation history curve;Congruent point generation unit 1122, for according to mean place information and averaged deformation history
Curve generates congruent point.
On the basis of above-described embodiment, in other embodiment, referring to Figure 20, which further includes output module 20.
Wherein, output module 20 show congruent point for exporting, and different congruent points corresponds to different output modes, defeated
Go out pattern and include color and/or shape.
On the basis of above-described embodiment, in other embodiment, referring to Figure 21, which further includes the life of deformation history curve
Into module 30.
Wherein, deformation history curve generation module 30, for obtaining the congruent point without time point, and to different time points
Congruent point carry out linear fit processing, to obtain deformation history curve.
The other details of technical solution are realized about module each in above-mentioned seven embodiment terminal devices, reference can be made to above-mentioned reality
The description in the deformation data extracting method based on InSAR in example is applied, details are not described herein again.
It should be noted that each embodiment in this specification is described by the way of progressive, each embodiment weight
Point explanation is all difference from other examples, and just to refer each other for identical similar part between each embodiment.
For device class embodiment, since it is basicly similar to embodiment of the method, so description is fairly simple, related part is joined
See the part explanation of embodiment of the method.
The embodiment of the present application additionally provides a kind of storage medium, and for storing computer program, it includes for performing sheet
Apply for the computer program designed by the above-mentioned deformation data extracting method embodiment based on InSAR.It is situated between by performing the storage
The computer program stored in matter can realize the deformation data extracting method based on InSAR that the application provides.
The specific embodiment of invention is described in detail above, but it is only used as example, the present invention is not intended to limit
With specific embodiments described above.For a person skilled in the art, any equivalent modifications carried out to the invention
Or substitute also all among scope of the invention, therefore, the equalization made in the case where not departing from the spirit and principles in the present invention range
Transformation and modification, improvement etc., all should be contained within the scope of the invention.
Claims (10)
1. a kind of deformation data extracting method based on InSAR, which is characterized in that it includes:
A, multiple initial p S points of target construction are extracted;
B, cluster analysis processing is carried out to the multiple initial p S points, obtains multiple congruent points, record each congruent point correspondence
Initial p S points quantity;
C, the distance between two adjacent polymeric points are calculated, and calculate the number of the corresponding initial p S points of described two adjacent polymeric points
Measure difference;It judges whether to meet two adjacent polymeric points for presetting unusual congruent point decision condition, it is described default unusual poly-
Chalaza decision condition is less than pre-determined distance threshold value for the distance, and the number differences are more than preset quantity threshold value;
If D, in the presence of rejecting and the few congruent point of initial p S points corresponded in described two adjacent polymeric points, repeat the step
Rapid C, until there is no meet two adjacent polymeric points for presetting unusual congruent point decision condition;
If E, being not present, deformation data is extracted according to remaining congruent point.
2. the deformation data extracting method according to claim 1 based on InSAR, which is characterized in that step A, including:
Under preset coordinate system, multiple PS points in target zone are acquired;
Obtain space range information corresponding with the target construction and PS point three-dimensional space position information;
Spatial data is carried out to the multiple PS points according to the space range information and PS point three-dimensional space position information
Analyzing and processing, to extract multiple initial p S points corresponding with the target construction.
3. the deformation data extracting method according to claim 1 based on InSAR, which is characterized in that step B, including:
B1, two adjacent initial p S points form a PS point pair, and analyse whether there is the PS points pair for meeting default fusion conditions;
If B2, in the presence of described new to generate new PS points to meeting the PS points of default fusion conditions to carrying out fusion treatment
PS points correspond to two initial p S points, the step B1 is repeated, until there is no meet the PS of the default fusion conditions
Point pair;
If B3, being not present, the corresponding all initial p S points of each remaining PS point are obtained, and to all initial p S points
Polymerization processing is carried out, obtains a congruent point.
4. the deformation data extracting method according to claim 3 based on InSAR, which is characterized in that analyzed in step B1
The PS points pair of fusion conditions are preset with the presence or absence of satisfaction, including:
The relevant parameter of the PS points pair is calculated, the relevant parameter includes space length and deformation deviation standard deviation;
Judge whether that the space length is less than pre-set space distance threshold, and the deformation deviation standard deviation is less than default
The PS points pair of standard deviation threshold method;
To be less than pre-set space distance threshold, and the deformation deviation standard deviation is poor less than preset standard when there are the space lengths
During the PS points pair of threshold value, then PS points are to meeting the default fusion conditions.
5. the deformation data extracting method according to claim 3 based on InSAR, which is characterized in that full in step B2
The PS points of the default fusion conditions of foot are to carrying out fusion treatment, to generate new PS points, including:
It obtains with PS points to the location information and deformation history curve of corresponding two initial p S points;
Mean place information is calculated according to the location information of described two initial p S points and deformation history curve and averaged deformation is gone through
History curve;
New PS points are generated according to the mean place information and the averaged deformation history curve.
6. the deformation data extracting method according to claim 3 based on InSAR, which is characterized in that institute in step B3
It states all initial p S points and carries out polymerization processing, obtain a congruent point, including:
Obtain the location information and deformation history curve of all initial p S points;
Mean place information is calculated according to the location information of all initial p S points and deformation history curve and averaged deformation is gone through
History curve;
According to the mean place information and averaged deformation history curve generation congruent point.
7. the deformation data extracting method according to claim 1 based on InSAR, which is characterized in that step E is wrapped later
It includes:
Output shows the congruent point, and different congruent points corresponds to different output modes, the output mode include color with/
Or shape.
8. the deformation data extracting method according to claim 1 based on InSAR, which is characterized in that step E is wrapped later
It includes:
The congruent point without time point is obtained, and linear fit processing is carried out to the congruent point of the different time points, to obtain
Deformation history curve.
9. a kind of terminal, which is characterized in that it includes processor, memory and display, and the processor couples the storage
Device, the display are stored with the computer program that can be run on the processor on the memory;
When the processor performs the computer program, the shape based on InSAR described in any one of claim 1-8 is realized
Become information extracting method.
10. a kind of storage medium, is stored thereon with computer program, which is characterized in that the computer program is held by processor
During row, the deformation data extracting method based on InSAR described in any one of claim 1-8 is realized.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711430532.9A CN108153979B (en) | 2017-12-26 | 2017-12-26 | Deformation information extraction method based on InSAR, terminal and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711430532.9A CN108153979B (en) | 2017-12-26 | 2017-12-26 | Deformation information extraction method based on InSAR, terminal and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108153979A true CN108153979A (en) | 2018-06-12 |
CN108153979B CN108153979B (en) | 2021-05-04 |
Family
ID=62462282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711430532.9A Active CN108153979B (en) | 2017-12-26 | 2017-12-26 | Deformation information extraction method based on InSAR, terminal and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108153979B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019126972A1 (en) * | 2017-12-26 | 2019-07-04 | 深圳市城市公共安全技术研究院有限公司 | Deformation information extraction method using insar, terminal, and storage medium |
CN109991601A (en) * | 2018-11-29 | 2019-07-09 | 北京东方至远科技股份有限公司 | A kind of house methods of risk assessment based on PS-InSAR technology |
WO2019220574A1 (en) * | 2018-05-16 | 2019-11-21 | 日本電気株式会社 | Synthetic aperture radar signal analysis device, synthetic aperture radar signal analysis method, and synthetic aperture radar signal analysis program |
CN111308468A (en) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | Method for automatically identifying deformation risk area based on In SAR technology |
CN112923901A (en) * | 2021-02-04 | 2021-06-08 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | RS and GIS technology-based engineering settlement calculation method and device |
CN114809389A (en) * | 2022-04-26 | 2022-07-29 | 中铁二十局集团第三工程有限公司 | Collapse prevention construction control method, equipment, system and medium for steel structure |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551455A (en) * | 2009-05-13 | 2009-10-07 | 西安电子科技大学 | 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof |
US7768441B2 (en) * | 2008-02-06 | 2010-08-03 | Halliburton Energy Services, Inc. | Geodesy via GPS and INSAR integration |
CN103308925A (en) * | 2013-05-31 | 2013-09-18 | 中国科学院合肥物质科学研究院 | Integral three-dimensional color laser radar data point cloud generating method and device thereof |
CN103822598A (en) * | 2014-02-26 | 2014-05-28 | 北京理工大学 | Deformation monitoring method of foundation SAR in time decorrelation serious area |
CN106679620A (en) * | 2017-01-13 | 2017-05-17 | 首都师范大学 | Regional land subsidence space-time prediction method |
CN106772377A (en) * | 2017-01-18 | 2017-05-31 | 深圳市路桥建设集团有限公司 | A kind of building deformation monitoring method based on InSAR |
CN106940443A (en) * | 2017-01-16 | 2017-07-11 | 洪都天顺(深圳)科技有限公司 | Complicated city infrastructure PSInSAR deformation methods of estimation under the conditions of cloud-prone and raining |
-
2017
- 2017-12-26 CN CN201711430532.9A patent/CN108153979B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7768441B2 (en) * | 2008-02-06 | 2010-08-03 | Halliburton Energy Services, Inc. | Geodesy via GPS and INSAR integration |
CN101551455A (en) * | 2009-05-13 | 2009-10-07 | 西安电子科技大学 | 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof |
CN103308925A (en) * | 2013-05-31 | 2013-09-18 | 中国科学院合肥物质科学研究院 | Integral three-dimensional color laser radar data point cloud generating method and device thereof |
CN103822598A (en) * | 2014-02-26 | 2014-05-28 | 北京理工大学 | Deformation monitoring method of foundation SAR in time decorrelation serious area |
CN106679620A (en) * | 2017-01-13 | 2017-05-17 | 首都师范大学 | Regional land subsidence space-time prediction method |
CN106940443A (en) * | 2017-01-16 | 2017-07-11 | 洪都天顺(深圳)科技有限公司 | Complicated city infrastructure PSInSAR deformation methods of estimation under the conditions of cloud-prone and raining |
CN106772377A (en) * | 2017-01-18 | 2017-05-31 | 深圳市路桥建设集团有限公司 | A kind of building deformation monitoring method based on InSAR |
Non-Patent Citations (4)
Title |
---|
GERNHARDT S等: "Deformation monitoring of single buildings using meter - resolution SAR data in PSI", 《SPRS JOURNAL OF PHOTOGRAMMETRY AND REMOTE SENSING 》 * |
M. COSTANTINI等: "Persistent scatterer pair interferometry: approach and application to COSMO-SkyMed SAR data", <IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING> * |
王淦: "PSInSAR 技术中 PS 目标点提取方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
秦晓琼等: "应用高分辨率 PS-InSAR 技术监测上海动迁房歪斜形变", 《测绘通报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019126972A1 (en) * | 2017-12-26 | 2019-07-04 | 深圳市城市公共安全技术研究院有限公司 | Deformation information extraction method using insar, terminal, and storage medium |
WO2019220574A1 (en) * | 2018-05-16 | 2019-11-21 | 日本電気株式会社 | Synthetic aperture radar signal analysis device, synthetic aperture radar signal analysis method, and synthetic aperture radar signal analysis program |
JPWO2019220574A1 (en) * | 2018-05-16 | 2021-05-13 | 日本電気株式会社 | Synthetic aperture radar signal analyzer, synthetic aperture radar signal analysis method and synthetic aperture radar signal analysis program |
JP7006781B2 (en) | 2018-05-16 | 2022-01-24 | 日本電気株式会社 | Synthetic Aperture Radar Signal Analysis Device, Synthetic Aperture Radar Signal Analysis Method and Synthetic Aperture Radar Signal Analysis Program |
US11835619B2 (en) | 2018-05-16 | 2023-12-05 | Nec Corporation | Synthetic aperture radar signal analysis device, synthetic aperture radar signal analysis method, and synthetic aperture radar signal analysis program |
CN109991601A (en) * | 2018-11-29 | 2019-07-09 | 北京东方至远科技股份有限公司 | A kind of house methods of risk assessment based on PS-InSAR technology |
CN111308468A (en) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | Method for automatically identifying deformation risk area based on In SAR technology |
CN112923901A (en) * | 2021-02-04 | 2021-06-08 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | RS and GIS technology-based engineering settlement calculation method and device |
CN112923901B (en) * | 2021-02-04 | 2022-04-15 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | RS and GIS technology-based engineering settlement calculation method and device |
CN114809389A (en) * | 2022-04-26 | 2022-07-29 | 中铁二十局集团第三工程有限公司 | Collapse prevention construction control method, equipment, system and medium for steel structure |
Also Published As
Publication number | Publication date |
---|---|
CN108153979B (en) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153979A (en) | Deformation information extraction method based on InSAR, terminal and storage medium | |
CN111623722B (en) | Multi-sensor-based slope deformation three-dimensional monitoring system and method | |
Liu et al. | Automated matching of multiple terrestrial laser scans for stem mapping without the use of artificial references | |
CN105550199B (en) | A kind of point polymerization and device based on multi-source map | |
CN111985494B (en) | Road disease detection and classification method based on lightweight vehicle-mounted terminal | |
WO2019126972A1 (en) | Deformation information extraction method using insar, terminal, and storage medium | |
CN109991601A (en) | A kind of house methods of risk assessment based on PS-InSAR technology | |
CN103605978A (en) | Urban illegal building identification system and method based on three-dimensional live-action data | |
CN111583244B (en) | Bridge deformation detection method and system | |
CN109238227B (en) | Method for representing ground settlement space-time evolution | |
Chmelina et al. | A 3-d laser scanning system and scan data processing method for the monitoring of tunnel deformations | |
CN113240735B (en) | Slope displacement activity monitoring method | |
CN108830317B (en) | Rapid and fine evaluation method for joint attitude of surface mine slope rock mass based on digital photogrammetry | |
Mayr et al. | Multitemporal analysis of objects in 3D point clouds for landslide monitoring | |
CN109579827A (en) | A kind of magnetic target detection and localization method based on arcuate array | |
Zhou et al. | Vision-based modal parameter identification for bridges using a novel holographic visual sensor | |
CN117076928A (en) | Bridge health state monitoring method, device and system and electronic equipment | |
Xin et al. | Landslide surface horizontal displacement monitoring based on image recognition technology and computer vision | |
CN102254185A (en) | Background clutter quantizing method based on contrast ratio function | |
Chen et al. | A critical review of automated extraction of rock mass parameters using 3D point cloud data | |
CN116561536B (en) | Landslide hidden danger identification method, terminal equipment and medium | |
CN110969601B (en) | Structure rotation response non-contact identification method based on visual characteristic tracking algorithm | |
CN113156456A (en) | Pavement and tunnel integrated detection method and detection equipment and vehicle | |
Liu et al. | Displacement field reconstruction in landslide physical modeling by using a terrain laser scanner–Part 1: Methodology, error analysis and validation | |
CN117433440A (en) | Roadway deformation real-time detection method and device based on laser point cloud and storage medium |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |