CN105842691A - Automatic adjustment and inlay method of large area InSAR deformation result and device thereof - Google Patents
Automatic adjustment and inlay method of large area InSAR deformation result and device thereof Download PDFInfo
- Publication number
- CN105842691A CN105842691A CN201610154331.XA CN201610154331A CN105842691A CN 105842691 A CN105842691 A CN 105842691A CN 201610154331 A CN201610154331 A CN 201610154331A CN 105842691 A CN105842691 A CN 105842691A
- Authority
- CN
- China
- Prior art keywords
- adjustment
- deformation
- big region
- insar
- equations
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
Abstract
An embodiment of the invention provides an automatic adjustment and inlay method of a large area InSAR deformation result and a device thereof. The method comprises the following steps of automatically identifying a homonymy high coherence point on a large area InSAR adjacent scene SAR image; using the automatically-identified homonymy high coherence point to establish an adjustment model of a large area deformation result and generating an adjustment equation; resolving the adjustment equation and acquiring an adjustment correction amount of a deformation rate of each single scene coverage area in the large area; and using the resolved adjustment correction amounts of the deformation rates of the plurality of single scene coverage areas in the large area to carry out automatic inlay on a large area deformation result. By using the above technical scheme, the scheme can be automatically realized; and large-area coordinated and consistent surface deformation information can be rapidly obtained so that large area surface deformation can be well monitored.
Description
Technical field
The present invention relates to InSAR (Synthetic Aperture Radar Interferometry, synthetic aperture radar interferometry)
Technical field, particularly relates to automatization's adjustment of a kind of big region InSAR deformation results and method for embedding and device.
Background technology
Difference synthetic aperture radar interferometry (D-InSAR, Differential Interferometric based on satellite repeat track
SAR) technology is the new technique of the Ground Deformation fields of measurement that the eighties in last century occurs, not only can obtain the precision of Centimeter Level,
And the information on energy acquisition face, compared to layering mark, GPS (Global Positioning System, global positioning system)
And the surface subsidence technology such as the measurement of the level has the biggest advantage.Traditional D-InSAR is limited by that space dephasing is dry, time dephasing
The factors such as dry and atmospheric interference, it is difficult to ensure its stability applied.In the time series that late 1990s occurs
InSAR technology breaches these limiting factors the most completely, is the once innovation to D-InSAR technology.Time series InSAR
The feature of technology is: process based on time series SAR image；The object processed is not whole pixels of view picture image,
But wherein there is stable scattering properties or keep high relevant pixel subset (to call in the following text " high relevant in the long period is spaced
Point ").Time series InSAR technology may be summarized to be two classes generally: interferes (Permanent Scatterer with Permanent scatterers
Or Persistent Scatterer Interferometry or PS-InSAR) it is single main image time series InSAR of representative
Technology and be that representative is many with Small Baseline Subset technology (Small baseline subset interferometry, or SBAS InSAR)
Main image time series InSAR technology.Time series InSAR technology in volcano, earthquake, landslide, surface subsidence
Widely applied etc. in the deformation monitoring caused.
Surface subsidence is under nature and anthropic factor effect, causes regional ground elevation due to the compression of the veneer of crust soil body
A kind of geological phenomenon reduced, is a kind of special Ground Deformation, extensively occurs in China.Early 1990s, Shanghai,
16 provinces (district, city) the surface subsidence areas such as Tianjin, Beijing, Jiangsu, Zhejiang, Hebei are about 48700km2, by 2003
Reach 93855km2, define 4 surface subsidence such as Yangtze River Delta, the North China Plain, the Delta of the Pearl River and Fen-Wei basin
The serious district of disaster.Visible, surface subsidence has become as a kind of regional disasters phenomenon.The time carried out the most both at home and abroad
Sequence InSAR deformation monitoring is studied, and the overwhelming majority is only limitted to the relatively small region that a scape SAR image is covered.Pin
The InSAR deformation monitoring needing the big region covered by many scapes SAR image the most rails (track) SAR image is ground
Study carefully seldom.
How on the basis of single scape time series InSAR deformation data extracts, carry out many scapes SAR image deformation results
Adjustment and inlaying, to realize the monitoring of big regional deformation, this is the key issue that those skilled in the art are urgently to be resolved hurrily.
Summary of the invention
The embodiment of the present invention provides automatization's adjustment of a kind of big region InSAR deformation results and method for embedding and device, to solve
The certainly auto_mosaicking problem of many scapes SAR image deformation results, thus preferably realize the Ground Deformation monitoring of big region.
On the one hand, embodiments provide automatization's adjustment and the method for embedding of a kind of big region InSAR deformation results,
Described method includes:
High coherent point of the same name on big region InSAR adjacent scape SAR image is identified automatically；
Utilize the of the same name high coherent point automatically identified, set up the adjustment Models of described big regional deformation result, generate adjustment side
Journey；
Described adjustment Equations is resolved, obtains the adjustment correction of the rate of deformation of each single scape area of coverage in described big region
Amount；
Utilize the adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region calculated, to described big region
Deformation results carries out auto_mosaicking.
On the other hand, embodiments provide automatization's adjustment of a kind of big region InSAR deformation results and inlay dress
Putting, described device includes:
Automatic identification unit, for automatically identifying high coherent point of the same name on big region InSAR adjacent scape SAR image；
Adjustment Models sets up unit, for utilizing the of the same name high coherent point automatically identified, sets up described big regional deformation result
Adjustment Models, generate adjustment Equations；
Adjustment Equations solving unit, for resolving described adjustment Equations, obtains each single scape in described big region and covers
The adjustment correction amount of the rate of deformation in district；
Auto_mosaicking unit, for utilizing the adjustment of the rate of deformation of multiple single scape areas of coverage in the described big region calculated to change
Positive quantity, carries out auto_mosaicking to described big regional deformation result.
Technique scheme has the advantages that and can automatically realize, and can quickly obtain the ground that big regional coordination is consistent
Table deformation data, thus preferably realize the Ground Deformation monitoring of big region.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing skill
In art description, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only the present invention
Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to root
Other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 is automatization's adjustment and the method for embedding flow chart of a kind of big region of embodiment of the present invention InSAR deformation results；
Fig. 2 is that the 4 scape SAR images that the embodiment of the present invention uses cover (the filling such as oblique line, point, black of big area schematic
District represents two scape images, three scape images, the overlapping region of 4 four scape images respectively)；
Fig. 3 is Radarsat-2 image coverage condition and the bench mark distribution of embodiment of the present invention use；
Fig. 4 is the Jing-jin-ji region big region ground settlement speed mosaic map mosaic that the embodiment of the present invention obtains；
Fig. 5 is the Jing-jin-ji region being conventionally treated big region ground settlement speed mosaic map mosaic；
Fig. 6 is the Nogata of the difference of the rate of settling that InSAR and the measurement of the level obtain on 142 check points of application example of the present invention
Figure；
Fig. 7 is automatization's adjustment and the inlaying device structural representation of a kind of big region of embodiment of the present invention InSAR deformation results
Figure.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clearly and completely
Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on this
Embodiment in invention, the every other reality that those of ordinary skill in the art are obtained under not making creative work premise
Execute example, broadly fall into the scope of protection of the invention.
As it is shown in figure 1, be automatization's adjustment and the method for embedding of a kind of big region of embodiment of the present invention InSAR deformation results
Flow chart, described method includes:
101, high coherent point of the same name on big region InSAR adjacent scape SAR image is identified automatically；
Preferably, described on big region InSAR adjacent scape SAR image, high coherent point of the same name identifies automatically before,
Described method includes: the SAR image of described big region InSAR is used same type of Law of DEM Data DEM
Carry out geocoding.
The identification of high coherent point of the same name is adjustment and the basis inlayed.Although utilizing Image Matching technology can pass through adjacent S AR
The coupling of amplitude image completes the automatic identification of high coherent point of the same name, but the most oversize.It is true that it is same after utilizing geocoding
The high coherent point of name has this principle of identical geodetic coordinates can quickly solve this problem.Should be noted that all SAR images
Geocoding should use same type of DEM, such as SRTM DEM, so can guarantee that adjacent scape SAR image ground
Identical plane precision is had after reason coding.Fig. 2 showing, 4 scape SAR images cover the schematic diagram in big region, has more shadow
As cover larger area situation can the rest may be inferred.Rule of judgment based on geocoding secondary homonym high coherent point is:
| P1 (x, y)-P2 (x, y) | ＜ (Lx, Ly) (1)
Wherein (x, y) represents the geodetic coordinates after geocoding, and (x, y), (x y) represents adjacent S AR image upper to P2 to P1
In ((x y) represents the geodetic coordinates after geocoding, L for x, high coherent point y)x, LyRepresent thing after geocoding respectively
To the output Pixel-space with north-south.The high coherent point meeting condition (1) is considered as same place.
102, utilize the of the same name high coherent point automatically identified, set up the adjustment Models of described big regional deformation result, generate
Adjustment Equations；
Preferably, described adjustment Equations is made up of two class equations: control point adjustment Equations and high coherent point adjustment Equations of the same name；
When not having control point in described big region InSAR, the most described adjustment Equations is made up of high coherent point adjustment Equations of the same name.
Each single scape area of coverage tried to achieve by time series InSAR technology be relative to this scape in certain point relative shape
Variable Rate.In the case of not specifying the datum mark (reference point) of rate of deformation, time series InSAR technology is utilized to try to achieve
Linear deformation rate and actual value between exist a system deviation.This system deviation can be in big regional deformation result
Try to achieve during adjustment.The situation that above-mentioned 4 scape images are covered, it is assumed that the rate of deformation of single scape image is that (x y), puts down Vi
Difference reduction is Δ i, i=1 ..., 4.Equation in adjustment Models is made up of two classes, and a class is control point adjustment Equations, separately
One class is high coherent point adjustment Equations of the same name.Under outside ground control point measured data is supported, set up control point adjustment variance,
That is:
Vi (x, y)+Δ i=GCp_V (x, y) (2)
Wherein (x y) represents in control point (x, y) the deformation speed that place obtains via level or other means measurements GCp_V
Rate value.They are the reference datas of the rate of deformation that InSAR obtains.
For the of the same name high coherent point of adjacent scape image, set up high coherent point adjustment Equations of the same name, it may be assumed that
Vi (x, y)+Δ i=Vj (x, y)+Δ j, i, j ∈ { 1,2,3,4} (3)
After this equation represents adjustment, same place on different images should have consistent rate of deformation.Respectively with tiltedly in Fig. 2
Line, point, black isotype marked two scape images, three scape images, the overlapping region of four scape images respectively.When point (x, y) position
In two scapes, three scapes, four scape images overlay region in time, the shape such as adjustment Equations of (3) formula has 1,2,3 respectively.
103, described adjustment Equations is resolved, obtain the flat of the rate of deformation of each single scape area of coverage in described big region
Difference reduction；
Preferably, described described adjustment Equations is resolved, obtain the deformation speed of each single scape area of coverage in described big region
The adjustment correction amount of rate, including: according to the principle of least square, described adjustment Equations is resolved, obtain in described big region
The adjustment correction amount of the rate of deformation of each single scape area of coverage.
Above-mentioned two class equations are united, constitute the system of linear equations about 4 unknown numbers (Δ 1, Δ 2 Δ 3, Δ 4).
Because same place number is a lot, the number of observational equation, far more than unknown number number, can solve according to the principle of least square
Unknown number, completes adjustment and the absolute calibration of the rate of deformation of every scape area of coverage the most simultaneously.Do not account for control point herein
Weight.There will be different control point in reality application and have the situation of different credibility.Now, can according to control point can
Reliability arranges different weights in formula (2), pays attention to, to improve calculation accuracy in least square resolves.
When monitoring section does not has control point, the most only Equations of The Second Kind equation.Adjustment now can be with the rate of deformation of arbitrary area of coverage
On the basis of value, having one degree of freedom in unknown number, therefore number is kept to 3.Rate of deformation after so adjustment completes is with true
A systematic bias is still there is between real-valued.
104, the adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region calculated is utilized, to described
Big regional deformation result carries out auto_mosaicking.
Preferably, the adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region that described utilization calculates,
Described big regional deformation result is carried out auto_mosaicking, including: utilize multiple single scape areas of coverage in the described big region calculated
The adjustment correction amount of rate of deformation, according to the geodetic coordinates after geocoding, and the overlapping number of times of same place, to described greatly
Regional deformation result carries out auto_mosaicking.
After completing the deformation results adjustment of multiple single scape area of coverage, according to the geodetic coordinates after geocoding, can complete
The auto_mosaicking of big regional deformation result.In locus, (x, y) deformation results at place is
In formula, N represents at position (x, y) the overlapping number of times of place's high coherent point.Corresponding in Fig. 2, at oblique line, point-like, black
In territory, zone, the value of N is 2,3,4 respectively.And in other regions, there is no an overlap, only one of which measured value, N's
Value is 1.Rate of deformation value in whole high coherent point in thus obtaining big region.
In order to effectiveness and the superiority of technical solution of the present invention are better described, existing to application example of the present invention and conventional treatment
Method carries out following relative analysis: as it is shown on figure 3, be the Radarsat-2 image coverage condition that uses of the embodiment of the present invention and
Bench mark scattergram, monitoring section is Beijing-tianjin-hebei Region, area about 100,000 square kilometres, totally 8 scape Radarsat-2 wide mode shadow
As covering, the overlay area of single scape image is about 150km*150km.Time series SAR on these 8 single scape areas of coverage
Image number and initial acquisition time are as shown in table 1.While utilizing InSAR to carry out surface subsidence monitoring, from Beijing
And the Mapping departments of Tianjin obtain some measurement of the level data.The measurement of the level measurement work annual 10 of Beijing and Tianjin
Moon first arrival completes by the end of November, obtains the three phase leveling datas of 2012 to 2014 altogether.Please note the effective of bench mark
Period is not completely superposed with the effectual time of SAR image.The position of these bench marks is shown in Fig. 3 with black circle,
Wherein the point of 1/3 represents the bench mark participating in compensating computation, and totally 70, the point of 2/3 represents the bench mark for precision test,
Totally 142.
Radarsat-2 image information used by table 1 Jing-jin-ji region surface subsidence monitoring
Utilize time series InSAR technology obtain radar line of sight on every scape area of coverage to averaged deformation speed, then change
Become the ground settlement speed of vertical direction.Wherein negative value represents surface subsidence, on the occasion of representing earth's surface lifting.Then by described flat
Difference method processes, and obtains the average sinking rate figure that calibrated big region is inlayed, and as shown in Figure 4, different single scapes cover
The rate of settling between cover region has good concordance, and this shows that the inventive method is effective and reliable.
Fig. 5 is the Jing-jin-ji region big regional deformation speed mosaic map mosaic that conventional treatment method obtains, due to only Beijing and Tianjin
The measurement of the level data in territory, are therefore only capable of the three scape rates of deformation to being numbered 1,2,3,4 in Fig. 3 and directly calibrate,
Remaining 4 scape rates of deformation can only rely on the same place between the most scaled data to carry out transmission and realize calibration indirectly.?
After when generating the mosaic map mosaic in big region due to the calibration precision of different scapes inconsistent will and occur significantly " splicing seams ".
Can see that, compared with Fig. 5 with Fig. 4 result, there is obvious error.
Utilize 142 bench marks that the big regional subsidence rate result shown in Fig. 3 has been carried out precision test.The two difference straight
Side's figure as shown in Figure 6, the mean square deviation of difference be 4.5 millimeters annual, difference maximum be 12.1 millimeters annual, minima is
-11.2 millimeters annual.Wherein the difference of the rate of settling on 90 points is positioned at [-3mm/year, 3mm/year], and this also illustrates this
The Jing-jin-ji region big regional subsidence rate result that invention application example obtains has the highest precision.Certainly, SAR image is initial
Time is not completely superposed with the initial time of bench mark, and the difference of the two calculated rate of settling also can be caused one by this
Fixing sound.
As fully visible, automatization's adjustment of application example of the present invention big region InSAR deformation results and method for embedding have two big
Advantage, the first can automatically realize, and can quickly obtain the Ground Deformation information that big regional coordination is consistent, its two be by with
Name point and external control point combine consideration, synchronize to realize adjustment and the absolute calibration of InSAR deformation monitoring result.The method exists
On a large scale Ground Deformation monitoring has important using value.
As it is shown in fig. 7, be automatization's adjustment and the inlaying device of a kind of big region of embodiment of the present invention InSAR deformation results
Structural representation, described device includes:
Automatic identification unit 71, for automatically knowing high coherent point of the same name on big region InSAR adjacent scape SAR image
Not；
Adjustment Models sets up unit 72, for utilizing the of the same name high coherent point automatically identified, sets up described big regional deformation
The adjustment Models of result, generates adjustment Equations；
Adjustment Equations solving unit 73, for resolving described adjustment Equations, obtains each single scape in described big region
The adjustment correction amount of the rate of deformation of the area of coverage；
Auto_mosaicking unit 74, in the described big region that utilization calculates, the rate of deformation of multiple single scape areas of coverage is flat
Difference reduction, carries out auto_mosaicking to described big regional deformation result.
Preferably, described device also includes: geocoding unit 70, for the SAR image to described big region InSAR
Same type of Law of DEM Data DEM is used to carry out geocoding.
Preferably, described adjustment Equations is made up of two class equations: control point adjustment Equations and high coherent point adjustment Equations of the same name；
When not having control point in described big region InSAR, the most described adjustment Equations is made up of high coherent point adjustment Equations of the same name.
Preferably, described adjustment Equations solving unit 73, further specifically for according to the principle of least square to described adjustment
Equation resolves, and obtains the adjustment correction amount of the rate of deformation of each single scape area of coverage in described big region.
Preferably, described auto_mosaicking unit 74, further specifically for utilizing multiple lists in the described big region calculated
The adjustment correction amount of the rate of deformation of the scape area of coverage, according to the geodetic coordinates after geocoding, and the overlapping number of times of same place,
Described big regional deformation result is carried out auto_mosaicking.
Technique scheme has the advantages that and can automatically realize, and can quickly obtain the ground that big regional coordination is consistent
Table deformation data, thus preferably realize the Ground Deformation monitoring of big region.
The particular order of the step during disclosed in should be understood that or level are the examples of illustrative methods.Based on design partially
Good, it should be appreciated that during the particular order of step or level can be in the case of without departing from the protection domain of the disclosure
Rearranged.Appended claim to a method gives the key element of various step with exemplary order, and is not meant to
It is limited to described particular order or level.
In above-mentioned detailed description, various features combine together in single embodiment, to simplify the disclosure.Should not
This open method is construed to reflect such intention by this, i.e. the embodiment of theme required for protection needs ratio clear
The more feature of feature stated in each claim to Chu.On the contrary, reflected such as appending claims that
Sample, the present invention is in the state fewer than whole features of disclosed single embodiment.Therefore, appending claims
Being hereby expressly incorporated in detailed description, wherein each claim is alone as the single preferred embodiment of the present invention.
For making any technical staff in this area be capable of or use the present invention, above disclosed embodiment is carried out
Describe.To those skilled in the art；The various amendment modes of these embodiments will be apparent from, and herein
The General Principle of definition can also be applicable to other embodiments on the basis of without departing from the spirit and scope of the disclosure.Cause
This, the disclosure is not limited to embodiments set forth herein, but with principle disclosed in the present application and the widest model of novel features
Enclose consistent.
Described above includes the citing of one or more embodiment.Certainly, describe to describe above-described embodiment parts or
The all possible combination of method is impossible, but it will be appreciated by one of ordinary skill in the art that each embodiment can
To do further combinations and permutations.Therefore, embodiment described herein is intended to fall into the guarantor of appended claims
All such changes, modifications and variations in the range of protecting.Additionally, with regard to the term " bag used in description or claims
Contain ", the mode that contains of this word is similar to term and " includes ", is used as link word in the claims just as " including, "
As being explained.Additionally, use any one term in the description of claims " or " be intended to represent " non-
Exclusive or ".
Those skilled in the art are it will also be appreciated that the various illustrative components, blocks (illustrative that list of the embodiment of the present invention
Logical block), unit, and step can pass through electronic hardware, computer software, or both combinations realize.For
Clearly show that the replaceability (interchangeability) of hardware and software, above-mentioned various illustrative components (illustrative
Components), unit and step the most universally describe their function.Such function is by hardware or soft
Part realizes depending on specifically applying the design requirement with whole system.Those skilled in the art can be specific for every kind
Application, it is possible to use the function described in the realization of various methods, but this realization is understood not to protect beyond the embodiment of the present invention
The scope protected.
Various illustrative logical block described in the embodiment of the present invention, or unit can pass through general processor, numeral
Signal processor, special IC (ASIC), field programmable gate array or other programmable logic device, discrete gate
Or transistor logic, discrete hardware components, or the design of any of the above described combination realize or operate described function.General
Processor can be microprocessor, and alternatively, this general processor can also be any traditional processor, controller, micro-
Controller or state machine.Processor can also realize by calculating the combination of device, such as digital signal processor and Wei Chu
Reason device, multi-microprocessor, one or more microprocessors one Digital Signal Processor Core of associating, or any other like
Configuration realize.
Method or the step of algorithm described in the embodiment of the present invention can be directly embedded into hardware, the software mould of processor execution
Block or the combination of both.Software module can be stored in RAM memory, flash memory, ROM memory, EPROM
In memorizer, eeprom memory, depositor, hard disk, moveable magnetic disc, CD-ROM or this area, other is any
In the storage medium of form.Exemplarily, storage medium can be connected with processor, so that processor can be from storage matchmaker
Jie reads information, it is possible to deposit write information to storage medium.Alternatively, storage medium can also be integrated in processor.
Processor and storage medium can be arranged in ASIC, and ASIC can be arranged in user terminal.Alternatively, processor and
Storage medium can also be arranged in the different parts in user terminal.
In one or more exemplary designs, the above-mentioned functions described by the embodiment of the present invention can hardware, software,
The combination in any of firmware or this three realizes.If realized in software, these functions can store the matchmaker with computer-readable
On Jie, or it is transmitted on the medium of computer-readable with one or more instructions or code form.Computer readable medium includes computer
Store medium and be easy to so that allowing computer program transfer to the telecommunication media in other place from a place.Storage medium can be
Any general or special computer can be with the useable medium of access.Such as, such computer readable media can include but not
It is limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device,
Other any may be used for carrying or storage with instruction or data structure and other can by general or special computer or general or
Special processor reads the medium of the program code of form.Additionally, any connection can be properly termed computer-readable
Medium, such as, if software is by coaxial cable, a light from a web-site, server or other remote resource
Fiber-optic cable, twisted-pair feeder, Digital Subscriber Line (DSL) or with the wireless way for transmittings such as the most infrared, wireless and microwave also by
It is included in defined computer readable medium.Described video disc (disk) and disk (disc) include Zip disk, radium-shine
Dish, CD, DVD, floppy disk and Blu-ray Disc, disk is generally with magnetic duplication data, and video disc generally carries out light with laser
Learn and replicate data.Combinations of the above can also be included in computer readable medium.
Above-described detailed description of the invention, has been carried out the purpose of the present invention, technical scheme and beneficial effect the most in detail
Illustrate, be it should be understood that the detailed description of the invention that the foregoing is only the present invention, be not intended to limit the present invention
Protection domain, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, all should wrap
Within being contained in protection scope of the present invention.
Claims (10)
1. automatization's adjustment of a big region InSAR deformation results and method for embedding, it is characterised in that described method bag
Include:
High coherent point of the same name on big region InSAR adjacent scape SAR image is identified automatically；
Utilize the of the same name high coherent point automatically identified, set up the adjustment Models of described big regional deformation result, generate adjustment side
Journey；
Described adjustment Equations is resolved, obtains the adjustment correction of the rate of deformation of each single scape area of coverage in described big region
Amount；
Utilize the adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region calculated, to described big region
Deformation results carries out auto_mosaicking.
2. as claimed in claim 1 automatization's adjustment of big region InSAR deformation results and method for embedding, it is characterised in that
Described on big region InSAR adjacent scape SAR image, high coherent point of the same name identifies automatically before, described method includes:
Same type of Law of DEM Data DEM is used to carry out geography the SAR image of described big region InSAR
Coding.
3. as claimed in claim 1 automatization's adjustment of big region InSAR deformation results and method for embedding, it is characterised in that
Described adjustment Equations is made up of two class equations: control point adjustment Equations and high coherent point adjustment Equations of the same name；When described big region
When not having control point in InSAR, the most described adjustment Equations is made up of high coherent point adjustment Equations of the same name.
4. as claimed in claim 1 automatization's adjustment of big region InSAR deformation results and method for embedding, it is characterised in that
Described described adjustment Equations is resolved, obtain the adjustment correction of the rate of deformation of each single scape area of coverage in described big region
Amount, including:
According to the principle of least square, described adjustment Equations is resolved, obtain the shape of each single scape area of coverage in described big region
The adjustment correction amount of variable Rate.
5. as claimed in claim 1 automatization's adjustment of big region InSAR deformation results and method for embedding, it is characterised in that
The adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region that described utilization calculates, to described big region
Deformation results carries out auto_mosaicking, including:
Utilize the adjustment correction amount of the rate of deformation of multiple single scape areas of coverage in the described big region calculated, according to geocoding
After geodetic coordinates, and the overlapping number of times of same place, described big regional deformation result is carried out auto_mosaicking.
6. automatization's adjustment of a big region InSAR deformation results and inlaying device, it is characterised in that described device bag
Include:
Automatic identification unit, for automatically identifying high coherent point of the same name on big region InSAR adjacent scape SAR image；
Adjustment Models sets up unit, for utilizing the of the same name high coherent point automatically identified, sets up described big regional deformation result
Adjustment Models, generate adjustment Equations；
Adjustment Equations solving unit, for resolving described adjustment Equations, obtains each single scape in described big region and covers
The adjustment correction amount of the rate of deformation in district；
Auto_mosaicking unit, for utilizing the adjustment of the rate of deformation of multiple single scape areas of coverage in the described big region calculated to change
Positive quantity, carries out auto_mosaicking to described big regional deformation result.
7. as claimed in claim 6 automatization's adjustment of big region InSAR deformation results and inlaying device, it is characterised in that
Described device also includes:
Geocoding unit, for using same type of digital elevation model to the SAR image of described big region InSAR
Data DEM carry out geocoding.
8. as claimed in claim 6 automatization's adjustment of big region InSAR deformation results and inlaying device, it is characterised in that
Described adjustment Equations is made up of two class equations: control point adjustment Equations and high coherent point adjustment Equations of the same name；When described big region
When not having control point in InSAR, the most described adjustment Equations is made up of high coherent point adjustment Equations of the same name.
9. as claimed in claim 6 automatization's adjustment of big region InSAR deformation results and inlaying device, it is characterised in that
Described adjustment Equations solving unit, further specifically for solving described adjustment Equations according to the principle of least square
Calculate, obtain the adjustment correction amount of the rate of deformation of each single scape area of coverage in described big region.
10. as claimed in claim 6 automatization's adjustment of big region InSAR deformation results and inlaying device, its feature exists
In,
Described auto_mosaicking unit, further specifically for utilizing the shape of multiple single scape areas of coverage in the described big region calculated
The adjustment correction amount of variable Rate, according to the geodetic coordinates after geocoding, and the overlapping number of times of same place, to described big region
Deformation results carries out auto_mosaicking.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610154331.XA CN105842691B (en) | 2016-03-17 | 2016-03-17 | The automation adjustment and method for embedding and device of big region InSAR deformation results |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610154331.XA CN105842691B (en) | 2016-03-17 | 2016-03-17 | The automation adjustment and method for embedding and device of big region InSAR deformation results |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105842691A true CN105842691A (en) | 2016-08-10 |
| CN105842691B CN105842691B (en) | 2018-01-12 |
Family
ID=56588348
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610154331.XA Expired - Fee Related CN105842691B (en) | 2016-03-17 | 2016-03-17 | The automation adjustment and method for embedding and device of big region InSAR deformation results |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105842691B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107610161A (en) * | 2017-10-10 | 2018-01-19 | 电子科技大学 | A kind of InSAR fast image registration methods based on Quadtree Partition |
| CN108564656A (en) * | 2017-12-11 | 2018-09-21 | 高力 | A kind of DSM data treating method and apparatus |
| CN114046774A (en) * | 2022-01-05 | 2022-02-15 | 中国测绘科学研究院 | Ground deformation continuous monitoring method integrating CORS network and multi-source data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090237297A1 (en) * | 2008-02-06 | 2009-09-24 | Halliburton Energy Services, Inc. | Geodesy Via GPS and INSAR Integration |
| KR20110048925A (en) * | 2009-11-03 | 2011-05-12 | 연세대학교 산학협력단 | The method for three-dimensional deformation measurement and the apparatus thereof |
| CN102520407A (en) * | 2011-12-29 | 2012-06-27 | 北京理工大学 | Polarimetric SAR (Synthetic Aperture Radar) interferometry image registration method based on optimal correlation coefficient linearity |
| CN103489176A (en) * | 2012-06-13 | 2014-01-01 | 中国科学院电子学研究所 | Method for extracting TPs from SAR image of serious geometric distortion |
| CN104239419A (en) * | 2014-08-20 | 2014-12-24 | 中国测绘科学研究院 | Discontinuous subnet connecting method and device oriented to time sequence InSAR |
-
2016
- 2016-03-17 CN CN201610154331.XA patent/CN105842691B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090237297A1 (en) * | 2008-02-06 | 2009-09-24 | Halliburton Energy Services, Inc. | Geodesy Via GPS and INSAR Integration |
| KR20110048925A (en) * | 2009-11-03 | 2011-05-12 | 연세대학교 산학협력단 | The method for three-dimensional deformation measurement and the apparatus thereof |
| CN102520407A (en) * | 2011-12-29 | 2012-06-27 | 北京理工大学 | Polarimetric SAR (Synthetic Aperture Radar) interferometry image registration method based on optimal correlation coefficient linearity |
| CN103489176A (en) * | 2012-06-13 | 2014-01-01 | 中国科学院电子学研究所 | Method for extracting TPs from SAR image of serious geometric distortion |
| CN104239419A (en) * | 2014-08-20 | 2014-12-24 | 中国测绘科学研究院 | Discontinuous subnet connecting method and device oriented to time sequence InSAR |
Non-Patent Citations (4)
| Title |
|---|
| LING ZHANG等: ""Seasonal subsidence retrieval with coherent point Target SAR Interferometry：A case study in Dezhou City"", 《GEOSCIENCE AND REMOTE SENSING SYMPOSIUM,2011 IEEE INTERNATIONAL》 * |
| M.FEKIR 等: ""Coherent targets selection based on small baseline interferogram subsets"", 《RECENT ADVANCES IN SPACE TECHNOLOGIES,2013 6TH INTERNATIONAL CONFERENCE ON》 * |
| 张永红 等: ""时间序列INSAR技术中的形变模型研究"", 《测绘学报》 * |
| 李莹莹 等: ""高分辨率SAR和可见图像同名点自动匹配技术"", 《测绘通报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107610161A (en) * | 2017-10-10 | 2018-01-19 | 电子科技大学 | A kind of InSAR fast image registration methods based on Quadtree Partition |
| CN108564656A (en) * | 2017-12-11 | 2018-09-21 | 高力 | A kind of DSM data treating method and apparatus |
| CN108564656B (en) * | 2017-12-11 | 2022-03-22 | 高力 | DSM data processing method and device |
| CN114046774A (en) * | 2022-01-05 | 2022-02-15 | 中国测绘科学研究院 | Ground deformation continuous monitoring method integrating CORS network and multi-source data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105842691B (en) | 2018-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| DeLong et al. | Multitemporal ALSM change detection, sediment delivery, and process mapping at an active earthflow | |
| Hsieh et al. | Using differential SAR interferometry to map land subsidence: a case study in the Pingtung Plain of SW Taiwan | |
| CN107144265B (en) | A kind of measuring method based on the graphic three-dimensional property right body space boundary mark of projection | |
| CN105842691A (en) | Automatic adjustment and inlay method of large area InSAR deformation result and device thereof | |
| CN110133655A (en) | Discharge of river based on multi-source radar remote sensing technology monitors inversion method | |
| Grigoriadis et al. | Gravity data inversion for Moho depth modeling in the Hellenic area | |
| Shi et al. | Early soil consolidation from magnetic extensometers and full resolution SAR interferometry over highly decorrelated reclaimed lands | |
| Baiocchi et al. | Accuracy of 1908 high to medium scale cartography of Rome and its surroundings and related georeferencing problems | |
| Fryer et al. | On the accuracy of heighting from aerial photographs and maps: implications to process modellers | |
| Yao et al. | Research on Surface Deformation of Ordos Coal Mining Area by Integrating Multitemporal D‐InSAR and Offset Tracking Technology | |
| CN113312695B (en) | Ancient landform restoration method based on trend surface analysis | |
| Kuok et al. | Flood map development by coupling satellite maps and three-dimensional drafting software: Case study of the Sarawak River Basin | |
| Ghomsi et al. | Crustal thickness beneath Atlas region from gravity, topographic, sediment and seismic data | |
| Widyaningrum et al. | Accuracy comparison of vhr systematic-ortho satellite imageries against vhr orthorectified imageries using GCP | |
| Parker | The integration of bathymetry, topography and shoreline and the vertical datum transformations behind it | |
| CN103759705B (en) | Three-point closed laying measurement method for composite water level | |
| CN108038086A (en) | A kind of dem data error assessment and bearing calibration based on pixel dimension | |
| Rifai et al. | Continuous monitoring of fluvial dike breaching by a laser profilometry technique | |
| Lee et al. | Transformation of vertical datum surface in the coastal area using hybrid geoid models | |
| Mo et al. | Evaluation of AW3D30 Elevation Accuracy in China. | |
| Yilmaz et al. | A research of consistencies and progresses of geoid models in Turkey | |
| Martin et al. | Lateral spread deformations from the 2011 Christchurch, New Zealand earthquake measured from Satellite images and optical image correlation | |
| Arabov | Improving the accuracy of the open access digital elevation models using GNSS: case Fergana valley | |
| Razumović et al. | Impact of height datum parameters on height reference system realizations: case study in Croatia | |
| Xu et al. | Implicit expansion rate of the Earth in global plate motions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180112 Termination date: 20190317 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |