CN107884772A - Method and device of the optimal interference images of time series InSAR to selection - Google Patents
Method and device of the optimal interference images of time series InSAR to selection Download PDFInfo
- Publication number
- CN107884772A CN107884772A CN201711012284.6A CN201711012284A CN107884772A CN 107884772 A CN107884772 A CN 107884772A CN 201711012284 A CN201711012284 A CN 201711012284A CN 107884772 A CN107884772 A CN 107884772A
- Authority
- CN
- China
- Prior art keywords
- width
- interference
- interference pattern
- pixel
- point target
- 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.)
- Withdrawn
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides a kind of time series InSAR optimal interference image to system of selection and device, is related to InSAR synthetic aperture radar interferometry technical fields.Method is：The selected characteristic region in monitored area；Each SAR striographs of characteristic area are generated into multi-frame interferometry figure；According to the signal to noise ratio of each pixel in every width SAR image features region, candidate's point target is selected；The common factor of candidate's point target of each SAR images is taken to be characterized the point target in region；The average coherence coefficient of the every width interference pattern of quantitative comparison, chooses optimal interference image pair.The present invention is by characteristic area coherence's quantitative comparison, overcoming and selecting the problem of interference image is to blindness in the prior art, avoid the leakage choosing of interference image pair and wrong choosing, realize the complete extraction of optimal interference image pair.
Description
Technical field
The present invention relates to InSAR (Synthetic Aperture Radar Interferometry, synthetic aperture radar
Interferometry) optimal interference image in technical field, more particularly to time series InSAR technologies is to system of selection and device.
Background technology
Synthetic aperture radar (SAR, Synthetic Aperture Radar) is to grow up most the 1950s
Important earth observation technology.By interferometry technology and SAR technologies with reference to and the synthetic aperture radar interferometry skill that is formed
Art (InSAR, Synthetic Aperture Radar Interferometry), which provides, obtains the brand-new of Three Dimensional Ground information
Method, it is observed or by common antenna parallel observation twice simultaneously by two slave antennas, obtains the multiple figure of the same landscape in ground
As right, according to phase difference of the ground each point in two width complex patterns, the path length difference of each point microwave in being imaged twice is drawn, so as to
Obtain the elevation information or deformation data of ground target.
Due to the InSAR technologies based on repeat track easily by space dephasing is dry, time dephasing is dry and atmospheric interference etc. because
The influence of element, it is difficult to monitored applied to stable Ground Deformation.In the prior art, in order to overcome conventional InSAR technology not
Foot, from the end of the nineties in last century, time series InSAR treatment technologies are suggested.Time series InSAR technologies are generally divided into two
Class：Interfere (Permanent Scatterer or Persistent Scatterer with Permanent scatterers
Interferometry, or PS-InSAR) for the single main image time series InSAR technologies of representative and with Small Baseline Subset technology
(Small baseline subset interferometry, or SBAS InSAR) is more main image time serieses of representative
InSAR technologies.
, it is necessary to select a number of interference pattern to carry out deformation parameter calculating in time series InSAR technologies.Where is selection
A little interference patterns are handled, and are had a significant impact for the precision of distortion measurement.
In process of the present invention is realized, inventor has found that at least there are the following problems in the prior art：
In PS-InSAR technologies, it is public main image that the problem of choosing which interference pattern, which changes into and selects which width image,
The interference map number that this method obtains is constant, i.e., M width image can only obtain M-1 width interference patterns.Interference pattern quantity is few, sternly
Ghost image rings the measurement accuracy of Ground Deformation.SBAS methods are in the prior art based on most important in time series InSAR technologies
Interference image is to system of selection.SBAS methods do not require to fix main image, but use space-time baseline hard threshold method selection interference image
It is right, hereafter referred to as space-time baseline hard threshold method.It overcomes PS-InSAR technologies and obtains the shortcomings that interference pattern quantity is few, it
By setting less vertical parallax threshold value and time reference line threshold value, all time reference lines and vertical parallax are less than given threshold
Interference pattern choose, so as to obtain more interference pattern.But space-time baseline hard threshold method selection interference image is blind to existing
Mesh, obtained interference pattern are of low quality.Reason is as follows：Sometimes, the interference pattern of time reference line very little is because the original of Changes in weather
Cause, coherence is very poor, but is less than given threshold and is selected；On the other hand, a small number of Space Baselines are smaller but time reference line is larger
The high quality interference pattern of (or Space Baseline is larger but time reference line is smaller) can not be chosen because exceeding baseline threshold.How to subtract
Wrong choosing and the leakage choosing of few interference image pair, the low-quality interference image pair of automatic rejection, to ensure that institute's generation interference pattern is high-quality does
Figure is related to, this is key issue urgently to be resolved hurrily in time series InSAR processing.
The content of the invention
The embodiment of the present invention provides one kind towards the optimal interference images of time series InSAR to automation system of selection and dress
Put.Its purpose is to the blindness for solving to select picture pair in the prior art, and the wrong choosing and leakage for reducing interference image pair are selected, so as to real
The complete extraction of existing high-quality interference image pair, beneficial to the monitoring of high-precision Ground Deformation.
On the one hand, the embodiments of the invention provide a kind of optimal interference image to include to automating system of selection, methods described：
The selected characteristic region in monitored area；
Each synthetic aperture radar SAR striographs of the characteristic area are generated into multi-frame interferometry figure；
According to the signal to noise ratio of each pixel in every width SAR image features region, the candidate point mesh of selection every width SAR images
Mark；
The common factor of candidate's point target of each SAR images is the point target of the characteristic area；
According to the point target of the characteristic area, the average coherence coefficient of every width interference pattern is calculated；
According to the average coherence coefficient of every width interference pattern, optimal interference image pair is chosen.
On the other hand, the embodiments of the invention provide a kind of optimal interference image to automating the device of selection, described device
Including：
First chooses unit, for choosing characteristic area in monitored area；
Into figure unit, in selected characteristic area, synthetic aperture radar SAR striographs to be generated into multi-frame interferometry figure；
Second chooses unit, for the signal to noise ratio according to each pixel in every width SAR image features region, chooses every width
Candidate's point target of SAR images；
Determining unit, for taking the common factor of candidate's point target of each SAR images to be defined as the point target of the characteristic area；
Computing unit, for the point target according to the characteristic area, calculate the average coherence coefficient of every width interference pattern；
3rd chooses unit, for the average coherence coefficient according to every width interference pattern, chooses optimal interference image pair.
One technical scheme of above-mentioned technical proposal has the advantages that：Because by each SAR images of characteristic area
Figure generation multi-frame interferometry figure, to choose optimal interference image to providing candidate's interference pattern as much as possible；Because according to every width SAR
The signal to noise ratio of each pixel in image feature region, chooses candidate's point target of every width SAR images, and takes the candidate of each SAR images
The common factor of point target is characterized the point target in region, avoids noise spot and enters candidate's point target, improves characteristic area point mesh
Target quality；Because according to the characteristic area point target of high quality, the average coherence coefficient of calculation interferogram；Quantitative comparison is interfered
The average coherence coefficient of figure chooses optimal interference pattern, and optimal interference image pair is determined by interference pattern, overcomes prior art selection
Interference pattern does not account for the blindness of the quality of actual interference pattern, reduces wrong choosing and the leakage choosing of interference image pair, and then
Optimal interference image improves the precision that Ground Deformation monitors to the Stability and veracity of selection.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is flow chart of the optimal interference image of the embodiment of the present invention to system of selection；
Fig. 2 is apparatus structure schematic diagram of the optimal interference image of the embodiment of the present invention to system of selection；
Fig. 3 is average amplitude figure of the embodiment of the present invention；
Fig. 4 is the interference image of selection of the embodiment of the present invention to baseline time-space distribution graph；
Fig. 5 is 78 interference images of prior art selection to time-space distribution graph；
Fig. 6 is different interference images that the embodiment of the present invention and prior art are chosen to coherence's distribution map；
Fig. 7 is 167 interference images of prior art selection to time-space distribution graph；
Fig. 8 is the interference pattern phase of two groups of selected interference images pair of the embodiment of the present invention；Wherein：Fig. 8 (a) is interference image pair
2 interference pattern；Fig. 8 (b) be interference image to 11 interference pattern；
Fig. 9 is the interference pattern phase of two groups of selected interference images pair of prior art；Wherein：Fig. 9 (a) is interference image to 31
Interference pattern；Fig. 9 (b) be interference image to 33 interference pattern；
Embodiment
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 clear, complete
Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under the premise of creative work is not made
Embodiment, belong to the scope of protection of the invention.
As shown in figure 1, be the optimal interference image of the embodiment of the present invention to automate system of selection flow chart.
101st, the selected characteristic region in monitored area；
102nd, each synthetic aperture radar SAR striographs of characteristic area are generated into multi-frame interferometry figure；
103rd, according to the signal to noise ratio of each pixel in every width SAR image features region, the candidate point mesh of every width SAR images is chosen
Mark；
104th, the common factor of candidate's point target of each SAR images is taken as the point target of the characteristic area；
105th, according to the point target of the characteristic area, the average coherence coefficient of every width interference pattern is calculated；
106th, according to the average coherence coefficient of every width interference pattern, optimal interference image pair is chosen.
Preferably, the characteristic area, it is representative high relevant characteristic area in monitored area.The model of characteristic area
Enclosing should be smaller, and the scope typically chosen is 600 × 600 pixels to 1000 × 1000 pixels, it is further preferred that selection 600 ×
600 pixels, 700 × 700 pixels, it is easy to the quick processing of data.It is further preferred that characteristic area include artificial works compared with
More Urban Areas, it is easy to choose more point target.
Preferably, each synthetic aperture radar SAR striographs by characteristic area generate multi-frame interferometry figure, specific bag
Include：
By any two width combination producing interference image pair of SAR striographs all in characteristic area；
Any two width of all M width SAR images combines to form N group interference images pair in the characteristic area：
By the N groups interference image of generation to all generating corresponding N width interference patterns.
Preferably, the signal to noise ratio according to each pixel in every width SAR image features region, every width SAR shadows are chosen
Candidate's point target of picture, is specifically included：
Obtain each pixel in every width SAR image features region；
Calculate the signal to noise ratio of each pixel in every width SAR image features region；
By the signal to noise ratio of pixel described in quantitative comparison, candidate's point target of every width SAR images is chosen；
Preferably, the signal to noise ratio of the pixel is the ratio of amplitude average and amplitude standards difference in the neighborhood centered on pixel
Value；
Wherein mp, σpAmplitude average and standard deviation in Square Neighborhood respectively centered on P.
The contiguous range is 5 × 5 pixels to 7 × 7 pixels；Preferably, Size of Neighborhood is 5 × 5 pixels or 7 × 7 pixels.
It is further preferred that the signal to noise ratio of quantitative comparison pixel, is specifically included：
The default snr threshold of setting；
The noise of every each pixel of width SAR images is contrasted one by one with predetermined threshold value；
When the signal to noise ratio of pixel is more than predetermined threshold value, the pixel is candidate's point target.
If local SNR SCRPMore than given threshold value, then P is defined as candidate's point target.It is special to choose every width SAR images
Pixel all in region is levied, chooses local SNR SCRPCandidate's point target of the SAR images is obtained more than given threshold value.
The common factor of candidate's point target in all M width SAR image features regions is taken, just obtains the point target in characteristic area.
Preferably, according to the point target of the characteristic area, the average coherence coefficient of every width interference pattern is calculated, it is specific to wrap
Include：
Each point target using each point target of characteristic area as every width interference pattern；
Calculate the coherence factor of each point target in characteristic area；
According to the coherence factor of each point target, the average coherence coefficient of every width interference pattern is calculated；
Preferably, the average coherence coefficient according to every width interference pattern, chooses optimal interference image pair, specifically includes：
The average coherence coefficient of each width interference pattern of quantitative comparison, chooses optimal interference pattern,
Interference image corresponding to the optimal interference pattern is to for optimal interference image pair；
It is further preferred that the average coherence coefficient of each width interference pattern of quantitative comparison, is specifically included：
The default coherence factor threshold value of setting；
The average coherence coefficient of every width interference pattern is contrasted one by one with predetermined threshold value；
When the average coherence coefficient of interference pattern is more than predetermined threshold value, the interference pattern is optimal interference pattern.
As shown in Fig. 2 for the optimal interference image of the embodiment of the present invention to automate system of selection apparatus structure schematic diagram.It is described
Device includes：
First chooses unit 21, for choosing characteristic area in monitored area；
Into figure unit 22, in selected characteristic area, synthetic aperture radar SAR striographs to be generated into multi-frame interferometry
Figure；
Second chooses unit 23, for the signal to noise ratio according to each pixel in every width SAR image features region, chooses described every
Candidate's point target of width SAR images；
Determining unit 24, for taking the common factor of candidate's point target of each SAR images to be defined as the point mesh of the characteristic area
Mark；
Computing unit 25, for the point target according to the characteristic area, calculate the average coherence coefficient of every width interference pattern；
3rd chooses unit 26, for the average coherence coefficient according to every width interference pattern, chooses optimal interference image pair.
Preferably, the characteristic area is representative high relevant characteristic area；The scope of characteristic area should be smaller, institute
The scope for stating characteristic area is 600 × 600 pixels to 1000 × 1000 pixels.It is further preferred that 600 × 600 pixels of selection,
700 × 700 pixels, it is easy to the quick processing of data.It is further preferred that characteristic area includes the more cities and towns of artificial works
Area, it is easy to choose more point target.
Preferably, it is described into figure unit 22, it is specifically used for：
By any two width combination producing interference image pair of SAR striographs all in characteristic area；
By all interference images of generation to all generating corresponding interference pattern.
Any two width of all M width SAR images combines to form N group interference images pair in the characteristic area：
By the N groups interference image of generation to all generating corresponding N width interference patterns.
Preferably, described second unit 23 is chosen, specifically included：
Acquisition module, for obtaining each pixel in every width SAR image features region；
Computing module, for calculating the signal to noise ratio of each pixel in every width SAR image features region；
First comparison module, for the signal to noise ratio of pixel described in quantitative comparison, choose the candidate point mesh of every width SAR images
Mark；
It is further preferred that first comparison module, is specifically used for：
Set snr threshold；
The signal to noise ratio of every each pixel of width SAR images is contrasted one by one with predetermined threshold value；
Choose the pixel signal to noise ratio be more than predetermined threshold value pixel be candidate's point target.To every width SAR image features
Pixel P in region, define local SNR
Wherein mp, σpAmplitude average and standard deviation in Square Neighborhood respectively centered on P, Size of Neighborhood are recommended to make
With 5 × 5 pixels or 7 × 7 pixels.If local SNR SCRPMore than given threshold value, then P is defined as candidate's point target.Traversal
Pixel in characteristic area, so that it may obtain candidate's point target of width SAR images.To being waited in whole M width SAR image features regions
Reconnaissance target takes common factor, then obtains the point target in characteristic area.
Preferably, the computing unit 25, is specifically used for：
Obtain each point target of each point target respectively as every width interference pattern of characteristic area；
Calculate the coherence factor of each point target in characteristic area；
Calculate the average coherence coefficient of every width interference pattern；
Preferably, the described 3rd unit 26 is chosen, specifically included：
Second comparison module, for the average coherence coefficient of each width interference pattern of quantitative comparison, choose optimal interference pattern；
Selecting module, for interference image corresponding to the optimal interference pattern will to be selected to for optimal interference image pair；
It is further preferred that second comparison module, is specifically used for：
Set coherence factor threshold value；
The average coherence coefficient of every width interference pattern is contrasted one by one with predetermined threshold value；
The interference pattern that the average coherence coefficient of selected interference pattern is more than predetermined threshold value is optimal interference pattern.
Above-mentioned technical proposal has the advantages that：By choose in whole monitored area it is representational it is high coherently
As characteristic area, the characteristic area of selection is 600 × 600 pixels, and scope is smaller, is easy to the quick processing of data, is chosen in area
Characteristic area include more artificial works, to ensure enough candidate's point targets.The characteristic area is all
Any two width combination producing whole interference pattern of SAR images, by the way of combination, obtained more interference image pair, so as to protect
The interference pattern quantity for having demonstrate,proved candidate is enough.According to the signal to noise ratio of each pixel in every width SAR image features region, every width is chosen
Candidate's point target of SAR images, and take the common factor of candidate's point target of each SAR images to be characterized the point target in region, avoid
Noise spot enters candidate's point target, improves the quality of characteristic area point target.Because according to the characteristic area of high quality point mesh
Mark, the average coherence coefficient of calculation interferogram；The average coherence coefficient of quantitative comparison interference pattern chooses optimal interference pattern, by dry
Relate to figure and determine optimal interference image pair, overcome the blindness that prior art selection interference pattern does not account for the quality of actual interference pattern
Property, reduce wrong choosing and the leakage choosing of interference image pair, and then Stability and veracity of the optimal interference image to selection, improve
The precision of Ground Deformation monitoring.
Above-mentioned technical proposal of the embodiment of the present invention is described in detail below in conjunction with application example：
In order to which the superiority of technical solution of the present invention is better described, application example is lifted below and is illustrated：Use this hair
Bright method and prior art SBAS methods (space-time baseline hard threshold method) select interference image to comparing respectively.
Fig. 3 is refer to, the scape COSMO-SkyMed of the Tianjin middle part 29 drop rail images that application example of the present invention uses
Average amplitude figure, after orientation 3 regards depending on × distance to 3, image size is 10000 pixels (orientation) × 5733 pixel
(distance to), resolution ratio is about 5 meters, area coverage about 1500km2。
First, optimal interference image pair is chosen using the method for the present invention.Choose the representational high coherence area in Tianjin city
The pixel of domain 600 × 600 is as characteristic area.Optional two width combination producing, 406 groups of interference images pair of 29 scape SAR images, generation
406 groups of interference images are to generating 406 width interference patterns.Calculated using the square sliding window of contiguous range 5 × 5 in each characteristic area
Local SNR, given threshold 2.0, the local SNR of each pixel of quantitative comparison；If the local SNR of pixel is big
In 2.0, then the pixel is candidate's point target；The common factor of candidate's point target in 29 scape SAR image features regions is taken as feature
The point target in region, 7276 altogether.According to the coherence factor of 7276 characteristic area point targets, the flat of every width interference pattern is calculated
Equal coherence factor, the average coherence coefficient of 406 width interference patterns is calculated altogether；The threshold value of average coherence coefficient is set as 0.7, is selected
Average coherence coefficient is more than 0.7 78 interference patterns, selects 78 groups of interference images altogether to for optimal interference image pair.It refer to Fig. 4,
Fig. 4 is the interference image of selection of the embodiment of the present invention to baseline time-space distribution graph.Fig. 4 abscissa is time reference line away from unit
" day ", " 01/01/12 " is on January 1st, 2012 in Fig. 4, and ordinate is spatial vertical baseline distance, unit " rice ", as being pressed to numbering
Coherence factor height sorts.
Then, optimal interference image pair is chosen using prior art.In order to choose the interference image with the bright method of we as much
It is right, the threshold value of time reference line is set to 300 days, vertical parallax threshold value is set to 450 meters.Table one is refer to, table one provides Tianjin
The time of the scape COSMO-SkyMed images of city middle part 29 and vertical parallax parameter：
Table one
The first row is refer to, the first scape is radix, and vertical parallax is away from, away from being 0, selection time baseline is small with time reference line
In 300 days and interference image pair of the vertical parallax less than 450 meters, 78 groups of interference images pair are obtained.Fig. 5 is refer to, Fig. 5 is existing skill
The interference image that art is chosen to time-space distribution graph, Fig. 5 abscissa for time reference line away from, ordinate is spatial vertical baseline distance, as
Numbering is in chronological sequence sorted.
Contrast it can be found that the optimal interference image of two methods selection is to the widely different of baseline spatial and temporal distributions, two kinds of sides
78 groups of interference images pair that method determines, wherein identical interference image is to there is 47 groups, and different interference images is to there is 31 groups, different interference
As being 39.7% to proportion.It refer to Fig. 6, Fig. 6 compared for the embodiment of the present invention and two methods of prior art obtain
31 groups of difference interference images are to coherence's distribution map, and left side is the 31 groups of interference images pair chosen using the method for prior art in Fig. 6,
Average coherence coefficient of the interference image to the point target of the interference pattern of generation, as seen from the figure, the average of its coherence factor is
0.6974；31 groups of interference images pair that right side is chosen for the inventive method in Fig. 6, point mesh of the interference image to the interference pattern of generation
Target average coherence coefficient, as seen from the figure, the average of its coherence factor is 0.7576；It can be seen that the inventive method is chosen as to obtaining
The coherence of interference pattern is selected as to the coherence of obtained interference pattern, absolutely proving, percentage is high apparently higher than prior art
Up to 39.7% different interference images pair, the interference image that the inventive method is chosen is to more excellent.
Below, we change the quality that a kind of angle contrast analyzes the interference image pair that two methods determine.Using the present invention's
Method choose 78 groups of interference images pair, its corresponding maximum time baseline be 896 days, maximum perpendicular baseline be 478 meters, therefore by when
Between the threshold value of baseline be set to 896 days, vertical parallax threshold value is set to 478 meters.Selection was less than 896 days and vertical parallax is less than 478 meters
Interference image pair, 167 interference images pair are obtained, refer to Fig. 7, Fig. 7 show prior art using 896 days+478 meters as when
167 interference images that space base line hard threshold method is chosen are to time-space distribution graph.As seen from the figure, can not possibly be chosen using prior art
The optimal interference image pair of the inventive method generation.
Next, we analyze interference image that two methods obtain to quality.First, prior art, which can leak, selects some interference
Quality is preferable, but baseline it is larger picture it is right.The interference image for the part high quality that the inventive method determines is adopted to not appearing in
The interference image centering chosen with prior art space-time, if the picture in Fig. 4 is to 2 (20140224-20150211) and as to 11
(20111030-20120930), please searched in Fig. 4 as to 2 and as the position to 11, its interference pattern according to Fig. 4 abscissa
Phase according to Fig. 8 as shown in figure 8, understand as preferable to 11 interference quality to 2 and picture, but as to 2 and as not having to 11
The interference image centering that present prior art is chosen.Second, prior art wrong can select the relatively low picture of some interference quality right.Such as Fig. 5
In picture to 31 (20130205-20130309) and as to 33 (20130205-20130917), its interference pattern phase such as Fig. 9 institutes
Show, as shown in Figure 9, the interference of the two interference images pair is second-rate.As to 31 and as there is not the inventive method choosing to 33
The interference image centering taken, illustrate that similar low quality interference image can exclude to the inventive method.
To find out its cause, prior art space-time baseline hard threshold method is true according to time reference line and Space Baseline dual thresholds
Determine interference image pair, there is blindness, selected interference image can not be ensured to high quality interference pattern can be generated.Under normal circumstances, do
Relate to as it is all smaller to time reference line, Space Baseline when, it is interfered, and plot quality is all higher, but this is not absolute.When major-minor image
When imaging moment weather meteorological condition differs greatly：Such as different vegetation growing way, rainfall snowfall, interference pattern coherence will be compared with
It is low.In this application example, by inquiring about Tianjin weather history data, there is rain and snow mixed Tianjin on January 31st, 2013,
There is within 2 3rd, 2013 light to moderate snow, have slight snow within 2 5th, 2013, continuous snowfall causes 20130205 (2 months 2013 5
Day) obtained image accumulated snow is more, interference pattern 31 (20130205-20130309) and interference pattern 33 in Fig. 5 formed therewith
(20130205-20130917) quality is relatively low.It follows that only could inherently it be evaluated by analyzing average coherence coefficient
The quality of interference pattern, and then filter out the interference preferable interference image pair of quality.
In summary, compared with prior art, the inventive method can avoid time dephasing is dry to be done with space dephasing completely
Bringing etc. influences, and avoids the leakage of optimal interference image pair from selecting and wrong choosing, extracts whole optimal interference images pair, ensure that interference image pair
The Stability and veracity of selection, be advantageous to the monitoring of high-precision Ground Deformation.
Those skilled in the art will also be appreciated that the various illustrative components, blocks that the embodiment of the present invention is listed
(illustrative logical block), unit, and step can pass through the knot of electronic hardware, computer software, or both
Conjunction is realized.To clearly show that the replaceability of hardware and software (interchangeability), above-mentioned various explanations
Property part (illustrative components), unit and step universally describe their function.Such work(
Can be that specific application and the design requirement of whole system are depended on to realize by hardware or software.Those skilled in the art
Various methods can be used to realize described function, but this realization is understood not to for every kind of specific application
Beyond the scope of protection of the embodiment of the present invention.
Various illustrative logical blocks described in the embodiment of the present invention, or unit can by general processor,
Digital signal processor, application specific integrated circuit (ASIC), field programmable gate array or other programmable logic devices, discrete gate
Or the design of transistor logic, discrete hardware components, or any of the above described combination is come the function described by realizing or operate.General place
It can be microprocessor to manage device, and alternatively, the general processor can also be any traditional processor, controller, microcontroller
Device or state machine.Processor can also be realized by the combination of computing device, such as digital signal processor and microprocessor,
Multi-microprocessor, one or more microprocessors combine a Digital Signal Processor Core, or any other like configuration
To realize.
The step of method or algorithm described in the embodiment of the present invention can be directly embedded into hardware, computing device it is soft
Part module or the combination of both.Software module can be stored in RAM memory, flash memory, ROM memory, EPROM storages
Other any form of storaging mediums in device, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area
In.Exemplarily, storaging medium can be connected with processor, to allow processor to read information from storaging medium, and
Write information can be deposited to storaging medium.Alternatively, storaging medium can also be integrated into processor.Processor and storaging medium can
To be arranged in ASIC, ASIC can be arranged in user terminal.Alternatively, processor and storaging medium can also be arranged at use
In different parts in the terminal of family.
In one or more exemplary designs, above-mentioned function described by the embodiment of the present invention can be in hardware, soft
Part, firmware or any combination of this three are realized.If realized in software, these functions can store and computer-readable
On medium, or with one or more instruction or code form be transmitted on the medium of computer-readable.Computer readable medium includes electricity
Brain storaging medium and it is easy to so that allowing computer program to be transferred to other local telecommunication medias from a place.Storaging medium can be with
It is that any general or special computer can be with the useable medium of access.For example, such computer readable media can include but
It is not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage devices, or other
What can be used for carrying or store with instruct or data structure and it is other can be by general or special computer or general or specially treated
The medium of the program code of device reading form.In addition, any connection can be properly termed computer readable medium, example
Such as, if software is to pass through a coaxial cable, fiber optic cables, double from a web-site, server or other remote resources
Twisted wire, Digital Subscriber Line (DSL) or with defined in being also contained in of the wireless way for transmitting such as infrared, wireless and microwave
In computer readable medium.Described disk (disk) and disk (disc) include Zip disk, radium-shine disk, CD, DVD, floppy disk
And Blu-ray Disc, disk is generally with magnetic duplication data, and disk generally carries out optical reproduction data with laser.Combinations of the above
It can also be included in computer readable medium.
Above-described embodiment, the purpose of the present invention, technical scheme and beneficial effect are carried out further
Describe in detail, should be understood that the embodiment that the foregoing is only the present invention, be not intended to limit the present invention
Protection domain, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., all should include
Within protection scope of the present invention.
Claims (10)
1. a kind of optimal interference images of time series InSAR are to system of selection, it is characterised in that methods described includes：
The selected characteristic region in monitored area；
Each synthetic aperture radar SAR striographs of the characteristic area are generated into multi-frame interferometry figure；
According to the signal to noise ratio of each pixel in every width SAR image features region, candidate's point target of selection every width SAR images；
The common factor of candidate's point target of each SAR images is taken as the point target of the characteristic area；
According to the point target of the characteristic area, the average coherence coefficient of every width interference pattern is calculated；
According to the average coherence coefficient of every width interference pattern, optimal interference image pair is chosen.
2. the optimal interference images of time series InSAR as claimed in claim 1 are to system of selection, it is characterised in that the feature
Region, it is representative, high relevant characteristic area；The scope of each characteristic area is 600 × 600 pixels to 1000 × 1000
Pixel.
3. the optimal interference images of time series InSAR as claimed in claim 1 are to system of selection, it is characterised in that by characteristic area
Each synthetic aperture radar SAR striographs generation multi-frame interferometry figure in domain, is specifically included：
By any two width combination producing interference image pair of SAR striographs all in characteristic area；
By all interference images of generation to all generating corresponding interference pattern.
4. the optimal interference images of time series InSAR as claimed in claim 1 are to system of selection, it is characterised in that the basis
The signal to noise ratio of each pixel in every width SAR image features region, candidate's point target of every width SAR images is chosen, is specifically included：
Obtain each pixel in every width SAR image features region；
Calculate the signal to noise ratio of each pixel in every width SAR image features region；
By the signal to noise ratio of pixel described in quantitative comparison, candidate's point target of every width SAR images is chosen；
The signal to noise ratio of the pixel is the ratio poor with amplitude standards of amplitude average in the neighborhood centered on pixel；
The contiguous range is 5 × 5 pixels to 7 × 7 pixels；
The signal to noise ratio by pixel described in quantitative comparison, specifically include：
The default snr threshold of setting；
The signal to noise ratio of every each pixel of width SAR images is contrasted one by one with predetermined threshold value；
When the signal to noise ratio of pixel is more than predetermined threshold value, the pixel is candidate's point target.
5. the optimal interference images of time series InSAR as claimed in claim 1 are to system of selection, it is characterised in that according to described
The point target of characteristic area, the average coherence coefficient of every width interference pattern is calculated, is specifically included：
Each point target using each point target of characteristic area as every width interference pattern；
Calculate the coherence factor of each point target in characteristic area；
According to the coherence factor of each point target, the average coherence coefficient of every width interference pattern is calculated；
The average coherence coefficient according to every width interference pattern, optimal interference image pair is chosen, is specifically included：
The average coherence coefficient of each width interference pattern of quantitative comparison, chooses optimal interference pattern,
Interference image corresponding to the optimal interference pattern is to for optimal interference image pair；
The average coherence coefficient of each width interference pattern of quantitative comparison, is specifically included：
The default coherence factor threshold value of setting；
The average coherence coefficient of every width interference pattern is contrasted one by one with predetermined threshold value；
When the average coherence coefficient of interference pattern is more than predetermined threshold value, the interference pattern is optimal interference pattern.
6. a kind of optimal interference images of time series InSAR are to the device of selection, it is characterised in that described device includes：
First chooses unit, for choosing characteristic area in monitored area；
Into figure unit, in selected characteristic area, synthetic aperture radar SAR striographs to be generated into multi-frame interferometry figure；
Second chooses unit, for the signal to noise ratio according to each pixel in every width SAR image features region, chooses every width SAR
Candidate's point target of image；
Determining unit, for taking the common factor of candidate's point target of each SAR images to be defined as the point target of the characteristic area；
Computing unit, for the point target according to the characteristic area, calculate the average coherence coefficient of every width interference pattern；
3rd chooses unit, for the average coherence coefficient according to every width interference pattern, chooses optimal interference image pair.
7. the optimal interference images of time series InSAR as claimed in claim 6 are to the device of selection, it is characterised in that the spy
Region is levied as representative high relevant characteristic area；The scope of the characteristic area is 600 × 600 pixels to 1000 × 1000
Pixel.
8. the optimal interference images of time series InSAR as claimed in claim 6 are to the device of selection, it is characterised in that it is described into
Figure unit, is specifically used for：
By any two width combination producing interference image pair of SAR striographs all in characteristic area；
By all interference images of generation to all generating corresponding interference pattern.
9. the optimal interference images of time series InSAR as claimed in claim 6 are to the device of selection, it is characterised in that described
Two choose unit, specifically include：
Acquisition module, for obtaining each pixel in every width SAR image features region；
Computing module, for calculating the signal to noise ratio of each pixel in every width SAR image features region；
First comparison module, for the signal to noise ratio of pixel described in quantitative comparison, choose candidate's point target of every width SAR images；
The signal to noise ratio of each pixel is the ratio of amplitude average and amplitude standards difference in neighborhood centered on pixel；
The contiguous range is 5 × 5 pixels to 7 × 7 pixels；
First comparison module, is specifically used for：
Set snr threshold；
The signal to noise ratio of every each pixel of width SAR images is contrasted one by one with predetermined threshold value；
Choose the pixel signal to noise ratio be more than predetermined threshold value pixel be candidate's point target.
10. the optimal interference images of time series InSAR as claimed in claim 6 are to the device of selection, it is characterised in that the meter
Unit is calculated, is specifically used for：
Obtain each point target of each point target respectively as every width interference pattern of characteristic area；
Calculate the coherence factor of each point target in characteristic area；
Calculate the average coherence coefficient of every width interference pattern；
Described 3rd chooses unit, specifically includes：
Second comparison module, for the average coherence coefficient of each width interference pattern of quantitative comparison, choose optimal interference pattern；
Selecting module, for interference image corresponding to the optimal interference pattern will to be selected to for optimal interference image pair；
Second comparison module, is specifically used for：
Set coherence factor threshold value；
The average coherence coefficient of every width interference pattern is contrasted one by one with predetermined threshold value；
The interference pattern that the average coherence coefficient of selected interference pattern is more than predetermined threshold value is optimal interference pattern.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711012284.6A CN107884772A (en) | 2017-10-26 | 2017-10-26 | Method and device of the optimal interference images of time series InSAR to selection |
| CN201810417761.5A CN108802729B (en) | 2017-10-26 | 2018-05-04 | Method and device for selecting time sequence InSAR optimal interference image pair |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711012284.6A CN107884772A (en) | 2017-10-26 | 2017-10-26 | Method and device of the optimal interference images of time series InSAR to selection |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107884772A true CN107884772A (en) | 2018-04-06 |
Family
ID=61782349
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711012284.6A Withdrawn CN107884772A (en) | 2017-10-26 | 2017-10-26 | Method and device of the optimal interference images of time series InSAR to selection |
| CN201810417761.5A Expired - Fee Related CN108802729B (en) | 2017-10-26 | 2018-05-04 | Method and device for selecting time sequence InSAR optimal interference image pair |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810417761.5A Expired - Fee Related CN108802729B (en) | 2017-10-26 | 2018-05-04 | Method and device for selecting time sequence InSAR optimal interference image pair |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN107884772A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108549080A (en) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | A kind of transmission tower position extracting method and system |
| CN109443532A (en) * | 2018-10-23 | 2019-03-08 | 北京工业大学 | A kind of mutual coherence factor test device of laser light field |
| CN117609533A (en) * | 2024-01-24 | 2024-02-27 | 航天宏图信息技术股份有限公司 | Automatic screening method, device and equipment for synthetic aperture radar interference image pair |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111220980B (en) * | 2020-01-19 | 2022-03-08 | 北京理工大学 | Ground-based SAR nonlinear atmospheric phase compensation method |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102034239B (en) * | 2010-12-07 | 2012-08-22 | 北京理工大学 | Local gray abrupt change-based infrared small target detection method |
| CN102608584B (en) * | 2012-03-19 | 2014-04-16 | 中国测绘科学研究院 | Time sequence InSAR (Interferometric Synthetic Aperture Radar) deformation monitoring method and device based on polynomial inversion model |
| DE102012204977B4 (en) * | 2012-03-28 | 2017-05-18 | Siemens Healthcare Gmbh | Method for iterative image reconstruction for bi-modal CT data and corresponding computation unit, computed tomography system, computer program and data carrier |
| US9148186B1 (en) * | 2014-04-08 | 2015-09-29 | Broadcom Corporation | Highly linear receiver front-end with thermal and phase noise cancellation |
| CN104035131A (en) * | 2014-05-23 | 2014-09-10 | 中国科学院电子学研究所 | Crust migration estimating method and device |
| CN104331899B (en) * | 2014-11-24 | 2018-06-19 | 中国科学院电子学研究所 | A kind of SAR image registration method and device |
-
2017
- 2017-10-26 CN CN201711012284.6A patent/CN107884772A/en not_active Withdrawn
-
2018
- 2018-05-04 CN CN201810417761.5A patent/CN108802729B/en not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108549080A (en) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | A kind of transmission tower position extracting method and system |
| CN109443532A (en) * | 2018-10-23 | 2019-03-08 | 北京工业大学 | A kind of mutual coherence factor test device of laser light field |
| CN117609533A (en) * | 2024-01-24 | 2024-02-27 | 航天宏图信息技术股份有限公司 | Automatic screening method, device and equipment for synthetic aperture radar interference image pair |
| CN117609533B (en) * | 2024-01-24 | 2024-04-05 | 航天宏图信息技术股份有限公司 | Automatic screening method, device and equipment for synthetic aperture radar interference image pair |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108802729A (en) | 2018-11-13 |
| CN108802729B (en) | 2021-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107884772A (en) | Method and device of the optimal interference images of time series InSAR to selection | |
| CN102608584B (en) | Time sequence InSAR (Interferometric Synthetic Aperture Radar) deformation monitoring method and device based on polynomial inversion model | |
| Hui et al. | Modelling spatial‐temporal change of Poyang Lake using multitemporal Landsat imagery | |
| CN109840553B (en) | Extraction method and system of cultivated land crop type, storage medium and electronic equipment | |
| CN103679202B (en) | Method and device suitable for classifying vegetation through optical remote sensing satellite images | |
| Yue et al. | A new accuracy evaluation method for water body extraction | |
| CN101470179A (en) | Method and apparatus for distortion calibration in magnetic resonance imaging | |
| CN104523274B (en) | A kind of MR imaging method utilizing steady state free precession | |
| CN110988909A (en) | TLS-based vegetation coverage determination method for sandy land vegetation in alpine and fragile areas | |
| CN111982822B (en) | Long-time sequence high-precision vegetation index improvement algorithm | |
| CN103489171A (en) | Wide-range remote-sensing image automatic dodging and color uniformizing method based on standard color library | |
| CN104239419B (en) | A kind of discontiguous subnetwork connection method and device towards time series InSAR | |
| CN101738172B (en) | Method for three-dimensional measurement of high sampling density color structured light based on green stripe segmentation | |
| CN112132946B (en) | Data extraction and display method for three-dimensional ground penetrating radar | |
| CN112254656A (en) | Stereoscopic vision three-dimensional displacement measurement method based on structural surface point characteristics | |
| CN110992366A (en) | Image semantic segmentation method and device and storage medium | |
| CN106991702A (en) | A kind of projector calibrating method and device | |
| Zhang et al. | A study on coastline extraction and its trend based on remote sensing image data mining | |
| CN108761458A (en) | Interference SAR water body digital elevation model modification method based on morphologic thinning | |
| CN104166112B (en) | MR image reconstruction method and device | |
| CN105842691B (en) | The automation adjustment and method for embedding and device of big region InSAR deformation results | |
| CN111539297B (en) | Semi-automatic extraction method for road information of high-resolution remote sensing image | |
| Bescoby | Detecting Roman land boundaries in aerial photographs using Radon transforms | |
| CN105044773A (en) | Converted-wave well-seismic calibration method and apparatus | |
| CN106908797A (en) | X-band areal rainfall radar data extracting method based on GIS |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180406 |