CN103698749B - A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers - Google Patents
A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers Download PDFInfo
- Publication number
- CN103698749B CN103698749B CN201310751219.0A CN201310751219A CN103698749B CN 103698749 B CN103698749 B CN 103698749B CN 201310751219 A CN201310751219 A CN 201310751219A CN 103698749 B CN103698749 B CN 103698749B
- Authority
- CN
- China
- Prior art keywords
- pixel
- sar image
- noise ratio
- domain signal
- amplitude
- 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.)
- Active
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
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 invention provides a kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers.Technical scheme utilizes the sequence comprising 8 to 15 width SAR image, first carries out the pre-service of SAR image sequence, then carries out spatial domain signal to noise ratio and estimate and the estimation of time domain signal to noise ratio, judge, obtain PS and extract result finally by PS.The invention has the beneficial effects as follows: realization is simple, accuracy is high, applied widely.
Description
Technical field
The present invention relates to the interleaving techniques field of microwave remote sensing and signal transacting, particularly one utilizes small data set SAR(Synthetic Aperture Radar, synthetic-aperture radar) method of image zooming-out Permanent scatterers.
Background technology
PS(Permanent Scatterers, Permanent scatterers) technology is a kind of typical long-term sequence differential interferometry technology, by process Spaceborne SAR System not in the same time to the same area obtain some width SAR image, millimetre-sized Ground Deformation information can be obtained, therefore be widely used in Natural calamity monitoring, Ground Deformation extraction, city prevention and cure of settlement etc., there is important using value.
PS technical finesse flow process comprises PS extraction, PS network struction, PS network parameter solve and optimize, the step such as Surface Deformation Field extraction, wherein, the processing accuracy of the follow-up each step of accuracy impact that PS extracts and sane degree, finally affect the extraction accuracy of Surface Deformation Field.
According to Italian scholar A.Ferretti to the definition of PS, PS refers to that, in longer time range, scattering phase place keeps stable point target.Existing PS extracting method has amplitude departure method and coefficient of coherence method etc.The basic thought of amplitude departure method is, due under given conditions, scattering amplitude and the scattering phase place of pixel have close degree of stability, and the amplitude stability by pixel extracts PS.This method needs the SAR image of more than process tens width to obtain effectively to extract result usually.Coefficient of coherence method utilizes the average coherence coefficient of pixel on whole image sequence as the evaluation index of its scattering stability, but the defect of this algorithm is the block clustering phenomena of PS easily occurred.By the restriction of SAR system data retrieval capabilities, in real world applications, often need the SAR image sequential extraction procedures regional deformation field utilizing small data set (picture number is 8 to 15 width).Be difficult to obtain gratifying extraction result because existing PS extracting method is applied to small data set SAR image sequence, at this moment we will adopt new, effective PS extracting method.
Summary of the invention
The technical problem to be solved in the present invention is, provides a kind of that be applicable to small data set, that effective PS extracts method.This method make use of the time-space domain information of SAR image data, thus decreases the requirement to SAR image number, extracts PS realize breaking through utilizing small data set SAR image.
The basic ideas of technical solution of the present invention are, because pixel signal to noise ratio and scattering phase stability exist positive correlation, first utilize the time-space domain information of SAR image data accurately to estimate the signal to noise ratio of each pixel in SAR image, then judge each pixel whether as PS according to the size of signal to noise ratio.
Technical scheme of the present invention is, based on Spaceborne SAR System not in the same time to M (8≤M≤15) the width SAR image sequence S that the same area obtains
i(1≤i≤M), carry out following process:
Step one: SAR image sequence pre-service.
The pre-service of SAR image sequence comprises SAR image sequence registration, resampling and radiant correction three steps, concrete grammar is with reference to Ph.D. Dissertation's " research based on Permanent scatterers radar differential interferometry search coverage Ground Deformation " (Southwest Jiaotong University, Chen Qiang, 2006.8,40th page to the 66th page) process.
Step 2: signal to noise ratio is estimated.This step comprises spatial domain signal to noise ratio and estimates and the estimation of time domain signal to noise ratio.
Wherein, spatial domain signal to noise ratio is estimated to comprise following three steps:
(1) step, the filtering of SAR image sequence time domain average.
To pretreated M width SAR image sequence S
i, utilize following formula to carry out time domain average filtering, obtain average amplitude figure
(2) step, estimates spatial domain signal to noise ratio.To average amplitude figure
following formula is utilized to estimate the noise intensity of a jth pixel
j represents average amplitude figure
pixel number.
Wherein, Ω
jfor the clutter estimation region of a jth pixel, use spatial domain windowing method selected (concrete grammar can with reference to " synthetic aperture radar target etection theory, algorithm and application " the 122 to 125 page).A
jkfor average amplitude figure
a jth pixel clutter estimation region Ω
jthe amplitude of an interior kth pixel, K is clutter estimation region Ω
jpixel number.
Suppose that a jth pixel is at average amplitude figure
amplitude be P
j, then its spatial domain signal to noise ratio
for:
In addition, time domain signal to noise ratio is estimated to comprise following two steps:
(1) step, estimates the amplitude deviation value of pixel.
Based on the SAR image sequence obtained after pre-service, following formula is utilized to estimate the amplitude deviation value of a jth pixel, j=1,2 ..., J:
Wherein, a
ijfor a jth pixel is at the SAR image amplitude in i moment.
(2) step, utilizes
estimate the time domain signal to noise ratio of a jth pixel
be shown below:
Step 3: PS judges.
Setting signal to noise ratio threshold value is 8db(decibel, decibel), according to the spatial domain signal to noise ratio of a jth pixel
with time domain signal to noise ratio
judge it whether as PS:
As
with
all be greater than threshold value, then judge that a jth pixel is as PS; Otherwise, judge that this pixel is as non-PS.After so all pixels being judged, final PS can be obtained and extract result.
The invention has the beneficial effects as follows: realization is simple, accuracy is high, applied widely.Owing to make use of spatial domain and the time-domain information of SAR image data, the present invention on extraction accuracy higher than existing extracting method.The present invention relaxes the requirement to SAR image number, even if utilize the SAR image sequence of small data set, also can obtain and extract result more accurately, and therefore realization is simple, usable range is wide, and this is that existing method is difficult to realize.
Accompanying drawing explanation
Fig. 1 is principle process schematic diagram of the present invention;
Fig. 2 utilizes the embodiment of the present invention to carry out the result figure tested;
Fig. 3 is for carrying out the conclusion of PS Accuracy evaluation to experimental result shown in Fig. 2.
Embodiment
Fig. 1 is principle process schematic diagram of the present invention.Technical scheme comprises carries out the pre-service of SAR image sequence to based on satellite-borne SAR image sequence, then carries out spatial domain signal to noise ratio and estimates and the estimation of time domain signal to noise ratio, judged, obtain PS and extract result by PS.
Fig. 2 utilizes the embodiment of the present invention to carry out the result figure tested.Wherein, SAR image corresponding region is the Anaheim city of the U.S., horizontal ordinate be orientation to (km), ordinate be distance to (km), the PS of extraction indicates with white round dot.Wherein, the SAR image sequence utilized comprises 12 width images, i.e. M=12, hollow windowing method is used to select clutter estimation region, the total size of rectangular window is 12 pixel × 12 pixels, window hollow area size 5 pixel × 5 pixel, i.e. pixel-5 pixel × 5, pixel number K=119=12 pixel × 12 pixel of clutter estimation region.
Fig. 3 is for carrying out the conclusion of PS Accuracy evaluation to experimental result shown in Fig. 2.
The PS accuracy index adopted is that (circular is with reference to international conference paper " Testing And Evaluation Of Permanent Scatterers Candidates Selection Methods " for PS network optimization rate, the National University of Defense Technology, Chen Shao strength, 2013.05).
SAR image sequence in utilizing Fig. 2 to test, adopt existing amplitude departure method and coefficient of coherence method (with reference to National University of Defense Technology grandson Chiron Ph.D. Dissertation " SAR chromatography and difference chromatography imaging technique are studied " the 30th page to the 40th page) and this method to extract PS respectively, then calculate the PS network optimization rate index of each method.Result of calculation shows, when process 12 scape SAR image, the PS accuracy index of this method is far above existing amplitude departure method and coefficient of coherence method, and this has absolutely proved the superiority that this method is utilizing small data set SAR image to extract PS.
Claims (2)
1. utilize a method for small data set SAR image sequential extraction procedures Permanent scatterers, based on Spaceborne SAR System not in the same time to the M width SAR image sequence S that the same area obtains
i, wherein SAR represents synthetic-aperture radar, 8≤M≤15, and 1≤i≤M, is characterized in that, comprises the steps:
Step one: SAR image sequence pre-service:
The pre-service of SAR image sequence comprises SAR image sequence registration, resampling and radiant correction three steps;
Step 2: signal to noise ratio is estimated:
This step comprises spatial domain signal to noise ratio and estimates and the estimation of time domain signal to noise ratio;
Wherein, spatial domain signal to noise ratio is estimated to comprise following two steps:
(1) step, the filtering of SAR image sequence time domain average:
To pretreated M width SAR image sequence S
i, utilize following formula to carry out time domain average filtering, obtain average amplitude figure
(2) step, estimate spatial domain signal to noise ratio:
To average amplitude figure
following formula is utilized to estimate the noise intensity of a jth pixel
j represents average amplitude figure
pixel number;
Wherein, Ω
jfor the clutter estimation region of a jth pixel, windowing method in spatial domain is used to select; a
jkfor average amplitude figure
a jth pixel clutter estimation region Ω
jthe amplitude of an interior kth pixel, K is clutter estimation region Ω
jpixel number;
Suppose that a jth pixel is at average amplitude figure
amplitude be P
j, then its spatial domain signal to noise ratio
for:
Wherein, time domain signal to noise ratio is estimated to comprise following two steps:
(1) step, estimate the amplitude deviation value of pixel:
Based on the SAR image sequence obtained after pre-service, following formula is utilized to estimate the amplitude deviation value of a jth pixel, j=1,2 ..., J:
Wherein, a
ijfor a jth pixel is at the SAR image amplitude in i moment;
(2) step, utilizes
estimate the time domain signal to noise ratio of a jth pixel
be shown below:
Step 3: PS judges;
Setting signal to noise ratio threshold value is 8 decibels, according to the spatial domain signal to noise ratio of a jth pixel
with time domain signal to noise ratio
judge it whether as PS:
As
with
all be greater than threshold value, then judge that a jth pixel is as PS; Otherwise, judge that this pixel is as non-PS; After so all pixels being judged, final PS can be obtained and extract result; Wherein PS represents Permanent scatterers.
2. the method utilizing small data set SAR image sequential extraction procedures Permanent scatterers according to claim 1, it is characterized in that, when using hollow windowing method to select clutter estimation region, the total size of rectangular window is 12 pixel × 12 pixels, window hollow area size 5 pixel × 5 pixel.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310751219.0A CN103698749B (en) | 2013-12-31 | 2013-12-31 | A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310751219.0A CN103698749B (en) | 2013-12-31 | 2013-12-31 | A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698749A CN103698749A (en) | 2014-04-02 |
CN103698749B true CN103698749B (en) | 2015-10-21 |
Family
ID=50360343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310751219.0A Active CN103698749B (en) | 2013-12-31 | 2013-12-31 | A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698749B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833971B (en) * | 2015-04-15 | 2017-06-20 | 北京理工大学 | Based on the Bistatic Radar System image PS point correlating methods for sliding scattering center |
CN107144213A (en) * | 2017-06-29 | 2017-09-08 | 中南大学 | The big magnitude three-D sequential deformation method of estimation in mining area and device based on SAR intensity images |
CN112797886B (en) * | 2021-01-27 | 2022-04-22 | 中南大学 | Winding phase oriented InSAR time sequence three-dimensional deformation monitoring method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE9804417L (en) * | 1998-12-18 | 1999-12-20 | Foersvarets Forskningsanstalt | A SAR radar system |
CN102144174B (en) * | 2008-07-04 | 2015-04-15 | 电视广播有限公司 | Identification and analysis of persistent scatterers in series of sar images |
-
2013
- 2013-12-31 CN CN201310751219.0A patent/CN103698749B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN103698749A (en) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bruni et al. | On the sensitivity of urban hydrodynamic modelling to rainfall spatial and temporal resolution | |
CN104360332B (en) | Atmospheric phase screen extraction method based on ground-based SAR (synthetic aperture radar) interference | |
Fan et al. | A new method of three‐dimensional location for low‐frequency electric field detection array | |
Rueben et al. | Optical measurements of tsunami inundation through an urban waterfront modeled in a large-scale laboratory basin | |
CN103698764B (en) | Interference synthetic aperture radar formation method under a kind of sparse sampling condition | |
CN104076404B (en) | Use multichannel to be concerned with and suppress the magnetic anomaly detection method of earth magnetism background noise | |
CN104316936B (en) | A kind of comprehensive DME pulse interference suppression method | |
CN103698749B (en) | A kind of method utilizing small data set SAR image sequential extraction procedures Permanent scatterers | |
WO2021051863A1 (en) | Water level pre-warning and monitoring method and apparatus, and storage medium and electronic device | |
CN104224188B (en) | For the faint breath signal detection method of life detection radar | |
CN106446919A (en) | Quick detection method for ground penetrating radar hyperbolic curve target | |
CN109190510A (en) | Underground cavity based on Ground Penetrating Radar quantifies recognition methods | |
CN105321163A (en) | Method and apparatus for detecting variation region of fully polarimetric SAR (Synthetic Aperture Radar) image | |
CN103675783A (en) | A broadband multiband imaging coherent processing method | |
Wang et al. | Subway tunnels identification through Cosmo-SkyMed PSInSAR analysis in Shanghai | |
Li et al. | Effectiveness evaluation of DS-InSAR method fused PS points in surface deformation monitoring: a case study of Hongta District, Yuxi City, China | |
Li et al. | Generalized Hough transform and ANN for subsurface cylindrical object location and parameters inversion from GPR data | |
CN104954298A (en) | Frequency estimation method for signals under impact noises and with data loss | |
Li et al. | Research on methods of high coherent target extraction in urban area based on PSINSAR technology | |
US11079461B2 (en) | Consistent arrival time measurement and determination of discharge polarity | |
CN203502585U (en) | Wireless monitoring system for mountain landslide monitoring | |
Yang et al. | Adaptive 4D PSI-based change detection | |
KR20170104100A (en) | System and Method For Ground Clutter Removing of WindProfiler | |
Wu et al. | Automatic Detection of Widely Distributed Local-Scale Subsidence Bowls in Rapidly Urbanizing Metropolitan Region Using Time-Series InSAR and Deep Learning Methods | |
Testi et al. | Successive-MFCW modulation for ultra-fast narrowband radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |