CN103424744A - Interference SAR layover area digital elevation model rebuilding method - Google Patents

Interference SAR layover area digital elevation model rebuilding method Download PDF

Info

Publication number
CN103424744A
CN103424744A CN2012101524772A CN201210152477A CN103424744A CN 103424744 A CN103424744 A CN 103424744A CN 2012101524772 A CN2012101524772 A CN 2012101524772A CN 201210152477 A CN201210152477 A CN 201210152477A CN 103424744 A CN103424744 A CN 103424744A
Authority
CN
China
Prior art keywords
folded
interference sar
zone
solution
covering
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
Application number
CN2012101524772A
Other languages
Chinese (zh)
Other versions
CN103424744B (en
Inventor
李芳芳
胡东辉
仇晓兰
丁赤飚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201210152477.2A priority Critical patent/CN103424744B/en
Publication of CN103424744A publication Critical patent/CN103424744A/en
Application granted granted Critical
Publication of CN103424744B publication Critical patent/CN103424744B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an interference SAR layover area digital elevation model rebuilding method. According to the method, the characteristic that the frequency of a layover area interferometric phase in the distance direction is opposite to that of areas on the two sides of a layover area is used, the frequency of the interferometric phases of the areas on the two sides of the layover area is estimated, the phase of a boundary of the layover area can be obtained, and accordingly an unwrapping phase inside the layover area is obtained, the DEM rebuilding problem of the layover area is solved, and an effective method is provided for topographic mapping achieved by means of an interference SAR in the layover area.

Description

The folded method of covering the reconstruction of zone digit elevation model of interference SAR
Technical field
The present invention relates to interference SAR topographic mapping field, relate in particular to a kind of folded method of covering zone digit elevation model (Digital Elevation Model, vehicle economy M) reconstruction of interference SAR.
Background technology
Interference synthetic aperture radar (Synthetic Aperture Radar, be called for short SAR) be to utilize the phase information of SAR complex data to obtain a technology of elevation information or the change information on earth's surface, there are round-the-clock, round-the-clock, high-precision characteristics, therefore in a plurality of fields, all be widely used, topographic mapping is that it mainly one of is applied.Yet, due to SAR side-looking imaging, making under precipitous landform or discontinuous surface conditions, the SAR image can form the folded phenomenon of covering.The folded image-region of covering phenomenon is appearring, the echo of a plurality of scattering sources of differing heights is added in a resolution element, thereby can cause error propagation at the folded facet phases of covering when discontinuous, cause traditional single base line interference SAR can't recover correct interferometric phase, thereby the precision that causes DEM to rebuild is lower.
Multi-baseline interference SAR is a kind of new technology proposed the nineties in 20th century, and the interference SAR altimetry precision can be provided.More document solves the folded regional mapping problem of covering with multi-baseline interference SAR, solution mainly can be divided into two classes: a class is that the method with spectrum estimation that the people such as F.Gini proposes is recovered to fold and covered regional phase place and elevation (Gini F, Lombardini F, Montanari M.Layover Solution in Multibaseline SAR Interferometry.IEEE Trans.on Aerospace and Electronic Systems.2002,38 (4): 1344-1356.); The method by statistics that the people such as the another kind of G.Ferraiuolo of being propose is carried out modeling to multi-baseline interference SAR signal, utilize the methods such as maximal possibility estimation, maximum a posteriori probability to estimate the folded regional elevation (G.Ferraiuolo that covers, V.Pascazio, G.Schirinzi.Maximum A Postseriori Estimation of Height Profiles in InSAR Imaging.IEEE Geoscience and Remote Sensing Letters.2004,1 (2): 66-70.).Yet, under complicated topographic condition, may there be a plurality of scattering sources stacks to form folded covering in practice, need more baseline number reach the purpose of differentiating a plurality of signals, therefore many Baselines are covered regional DEM and are difficult to realize in rebuilding actual interference SAR is folded.
Summary of the invention
(1) technical matters that will solve
For solving above-mentioned one or more problems, the invention provides a kind of folded method of covering the reconstruction of zone digit elevation model of interference SAR, cover folded the shortcoming that facet phases can cause error propagation when discontinuous to overcome the conventional interference phase unwrapping method, utilize single base line interference SAR to recover correct interferometric phase, realize that the folded zone digit elevation model of covering of interference SAR rebuilds.
(2) technical scheme
According to an aspect of the present invention, provide a kind of interference SAR the folded method of covering the reconstruction of zone digit elevation model, having comprised: the folded zone of covering of detecting the first interference SAR complex pattern and the second interference SAR complex pattern; The first interference SAR complex pattern and the second interference SAR complex pattern are carried out to conjugate multiplication, the result of conjugate multiplication is got to phase place, thereby obtain interferometric phase image; To fold and cover zone as mask, and on interferometric phase image, except folded, cover the extra-regional non-folded zone of covering and carry out the interferometric phase solution and twine, obtain non-ly foldedly covering regional solution and twining phase place; In each orientation upwards, by non-, foldedly cover regional solution and twine phase place, obtain respectively folded non-folded the first solution of covering zone of covering regional left border and twine phase place
Figure BDA00001645658700021
Twine phase place with non-folded the 3rd solution of covering zone of right side boundary
Figure BDA00001645658700022
In each orientation to, folded cover regional solution and twine phase place by non-, calculate folded cover regional left side along distance to front interferometric phase frequency f Before, and right side along the distance to rear interferometric phase frequency f AfterUpwards, by the first solution, twine phase place in each orientation
Figure BDA00001645658700023
Twine phase place with the 3rd solution
Figure BDA00001645658700024
And front phase frequency f BeforeWith rear phase frequency f After, calculate this orientation to the folded solution of covering regional left border twine phase place
Figure BDA00001645658700025
Twine phase place with the solution of right side boundary
Figure BDA00001645658700026
In each orientation upwards, twine phase place by the folded solution of covering regional left border
Figure BDA00001645658700027
Twine phase place with the solution of right side boundary Come this orientation of matching upwards to fold the solution of covering intra-zone and twine phase place
Figure BDA00001645658700029
In each orientation to each pixel, twine phase place by the solution of this pixel
Figure BDA000016456587000210
Digital elevation model h with this pixel of interference SAR calculation of parameter.
(3) beneficial effect
The present invention utilize folded cover regional interferometric phase along distance to the frequency characteristics contrary with its two side areas, by estimating the folded interferometric phase frequency of covering two side areas, can obtain the folded phase place of covering regional border, further matching obtains the folded solution of covering intra-zone and twines phase place, can cause the shortcoming of error propagation at the folded facet phases of covering when discontinuous thereby overcome the conventional interference phase unwrapping method, and only need to utilize the interference SAR data of single baseline can realize that the folded zone digit elevation model of covering of interference SAR rebuilds.
The accompanying drawing explanation
Fig. 1 is the folded process flow diagram of covering zone digit elevation model method for reconstructing of embodiment of the present invention interference SAR.
To be that SAR is folded cover regional imaging geometry figure to Fig. 2;
Fig. 3 be in Fig. 2 corresponding folded cover zone along distance to the interferometric phase schematic diagram;
The schematic diagram that Fig. 4 is the interference SAR geometric relationship;
Fig. 5 is the elevation schematic diagram for the treatment of simulating area;
The interferometric phase image that Fig. 6 simulates according to simulation parameter for the elevation schematic diagram by Fig. 5;
Fig. 7 is for folding and cover the resulting DEM reconstructed results of zone digit elevation model method for reconstructing schematic diagram according to embodiment of the present invention interference SAR;
Fig. 8 is along a certain orientation upwards, according to the folded resulting DEM of zone digit elevation model method for reconstructing and the desirable DEM correlation curve figure of covering of embodiment of the present invention interference SAR.
Embodiment
For making the purpose, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and, with reference to accompanying drawing, the present invention is described in more detail.
It should be noted that, in accompanying drawing or instructions description, similar or identical part is all used identical figure number.And in the accompanying drawings, to simplify or convenient the sign.Moreover the implementation that does not illustrate in accompanying drawing or describe, be form known to a person of ordinary skill in the art in affiliated technical field.In addition, although this paper can provide the demonstration of the parameter that comprises particular value, should be appreciated that, parameter is without definitely equaling corresponding value, but can in acceptable error margin or design constraint, be similar to corresponding value.
According to an aspect of the present invention, provide a kind of interference SAR the folded method of covering regional DEM reconstruction.Fig. 1 is the folded process flow diagram of covering zone digit elevation model method for reconstructing of embodiment of the present invention interference SAR.As shown in Figure 1, the present embodiment comprises the following steps:
Step S102: the first interference SAR complex pattern obtained by first passage and the second interference SAR complex pattern obtained by second channel are carried out to registration;
To the first interference SAR complex pattern of the first passage that is obtained from interference SAR system and second channel and the method that the second interference SAR complex pattern carries out registration, can be maximum coherence Y-factor method Y or maximum spectrum registration method.
Step S104: the first interference SAR complex pattern after the detection registration and the folded zone of covering of the second interference SAR complex pattern;
The folded zone of covering shows as obvious clear zone on magnitude image, and its corresponding gray-scale value is regional larger more than 3 times than other usually.According to these characteristics, just can detect the folded zone of covering in two width interference SAR complex patterns after registration.Due to two width interference SAR complex patterns registration, only need be detected the amplitude of any width interference SAR complex pattern wherein, the first interference SAR complex pattern of take is example, specifically comprises:
Sub-step S104a, ask for the averaged amplitude value M of whole pixels of this first interference SAR complex pattern a
Sub-step S104b, set a global threshold K a, make K a=cM a, wherein c is empirical value, value is 2≤c≤5 usually, preferably, c=3;
Sub-step S104c, be greater than K by amplitude in the first interference SAR complex pattern aPixel tentatively be judged as the folded zone of covering.
In addition, in order detecting more accurately to fold, to cover zone, can also to be revised preliminary judged result obtained above, can also comprise after above-mentioned steps S104c:
Sub-step S104d, calculate the coefficient of coherence of every part in the coefficient of coherence figure of the first interference SAR complex pattern and the second interference SAR complex pattern;
Wherein, the computing method of two width interference SAR complex pattern coefficient of coherence are as follows:
γ ( m , n ) = | Σ k = - P P Σ l = - Q Q s 1 | ( m + k , n + l ) s 2 * ( m + k , n + l ) [ Σ k = - P P Σ l = - Q Q | s 1 ( m + k , n + l ) | 2 ] [ Σ k = - P P Σ l = - Q Q | s 2 ( m + k , n + l ) | 2 ] - - - ( 1 )
Wherein, γ is the coefficient of coherence calculated, s 1And s 2Be two width interference SAR complex patterns, m, n be respectively orientation to distance while to ，Dui each position, position, calculating its coefficient of coherence all centered by this position, at image s 1And s 2In get (2P+1) * (2Q+1) window calculated.
Sub-step S104e, set the coefficient of coherence threshold k c
Sub-step S104f, folded the covering in zone coefficient of coherence higher than K gone out in Preliminary detection cPart remove, cover zone thereby obtain folding more accurately, this coefficient of coherence threshold k cValue be generally 0.9.
In addition, it should be noted that, detecting the folded zone of covering can also be judged according to the folded characteristics of covering regional interferometric phase.Certainly, those skilled in the art it will also be appreciated that other method, within should being included in protection scope of the present invention equally.In this step, by Detection and Extraction, fold and cover zone, just can when the subsequent phase solution twines, to it, carry out special processing.
Step S106: the first interference SAR complex pattern after registration and the second interference SAR complex pattern are carried out to conjugate multiplication, the result of conjugate multiplication is got to phase place, thereby obtain interferometric phase image;
Step S108: will fold and cover zone and do mask, and on interferometric phase image, except folded, cover the extra-regional non-folded zone of covering and carry out the interferometric phase solution and twine, obtain non-ly foldedly covering regional solution and twining phase place;
Wherein, the folded zone of covering is done after mask and done the purpose that the interferometric phase solution twines again and be: the phase error of having avoided the folded uncontinuity of covering regional phase place to cause is propagated.In addition, above-mentionedly to non-, folded cover zone to carry out the method that the interferometric phase solution twines can be Goldstein branch cutting method, or the path trace method that instructs of Quality Map, these methods are the common practise of this area, no longer repeat herein.
Step S110, non-ly folded covered regional solution and twine phase place by described, obtains in each orientation upwards, and folded non-folded the first solution of covering zone of covering regional left border twines phase place
Figure BDA00001645658700051
Folded non-folded the 3rd solution of covering zone of covering regional right side boundary twines phase place
Figure BDA00001645658700052
In each orientation to, get and folded cover adjacent non-folded of regional left border and cover regional solution and twine phase place and be Get and folded cover adjacent non-folded of regional right side boundary and cover regional solution and twine phase place and be
Step S112: in each orientation to, non-ly folded cover regional solution and twine phase place by described, calculate folded cover the regional left and right sides along distance to front interferometric phase frequency f BeforeWith rear interferometric phase frequency f After.
To an orientation wherein always saying, left side along distance to the interferometric phase frequency f BeforeConcrete method of estimation as follows:
Step S112a: to this orientation to, cover folded the estimation window that (2M+1) * (2N+1) chosen in regional left side, the phase model of this estimation window can be expressed as:
Figure BDA00001645658700055
Wherein,
Figure BDA00001645658700056
For window center point phase place, in window, the displacement of the relative window center of pixel is k, l, f x, f yBe respectively this window along the distance to orientation to the interferometric phase frequency.Here only need to estimate along distance to the interferometric phase frequency;
Step S112b: by maximization estimate the cost function of window phase model estimate this distance to the interferometric phase frequency:
Figure BDA00001645658700057
Wherein,
Figure BDA00001645658700058
To the interferometric phase frequency, estimate that the size of window is (2M+1) * (2N+1) for the distance estimating to obtain.Here
Figure BDA00001645658700059
Be the folded interferometric phase frequency f of covering the left side, zone Before.
With f BeforeMethod of estimation identical, S112a, to folded, covering the estimation window that (2M+1) * (2N+1) chosen on regional right side, press in each orientation, the method for S112b can estimate folds the interferometric phase frequency f of covering regional right side After.
In this step, the value of M and N is respectively 16,16.By this step, estimate the interferometric phase frequency, in order to twine phase place according to the folded folded solution of covering border of regional phase place characteristics calculating of covering.
Step S114: in each orientation upwards, according to folded the first solution of covering regional both sides, twine phase place
Figure BDA00001645658700061
Twine phase place with the 3rd solution
Figure BDA00001645658700062
And front interferometric phase frequency f BeforeWith rear interferometric phase frequency f After, calculate this orientation to the folded solution of covering regional left border and right side boundary twine phase place With
Figure BDA00001645658700064
To be that SAR is folded cover regional imaging geometry figure to Fig. 2.Zone A, B, C is added to from r in the SAR image E1To r E2Distance in resolution element, form the folded zone of covering.Fig. 3 be in Fig. 2 corresponding folded cover zone along distance to the interferometric phase schematic diagram.Part in Fig. 3 between two vertical dotted line is the folded zone of covering,
Figure BDA00001645658700065
Arrive
Figure BDA00001645658700066
Between the interferometric phase of phantom line segments corresponding region A,
Figure BDA00001645658700067
Arrive
Figure BDA00001645658700068
Between the interferometric phase of phantom line segments corresponding region C,
Figure BDA00001645658700069
Arrive
Figure BDA000016456587000610
Between the interferometric phase of real segment corresponding region B.Usually meet domatic scattering strength in actual SAR image larger, therefore fold and cover in zone, the interferometric phase of regional B is occupied an leading position.
As can be seen from Figure 3, folded cover the saltus step that has interferometric phase on regional border, calculate and foldedly cover the formula that the zone boundary solution twines phase place and be:
Figure BDA000016456587000611
Figure BDA000016456587000612
Here, suppose the image along continuous straight runs for distance to, vertical direction be orientation to, wherein,
Figure BDA000016456587000613
For folded cover regional left border non-folded obtained by step S110 covered regional solution and twines phase place,
Figure BDA000016456587000614
For folded cover regional right side boundary non-folded that step S110 obtains covered regional solution and twines phase place,
Figure BDA000016456587000615
For the folded true solution of covering regional right side boundary twines phase place,
Figure BDA000016456587000616
For the folded true solution of covering regional left border twines phase place.F BeforeFold for left side is non-and cover regional distance to interferometric phase frequency, f AfterFold for right side is non-and cover regional distance to interferometric phase frequency, N PixelFor folded cover zone along distance to pixel count.
Step S116: in each orientation upwards, according to the folded solution of covering zone boundary twine phase place come this orientation of matching upwards the folded solution of covering intra-zone twine phase place;
Each orientation obtained according to step S114 to foldedly cover regional facet phases With
Figure BDA000016456587000618
This orientation of linear fit twines phase place to the folded solution of covering intra-zone:
Figure BDA000016456587000619
Wherein,
Figure BDA000016456587000620
Figure BDA000016456587000621
For the folded facet phases of covering the regional left and right sides that step 4 obtains, N PixelFor folded cover zone along distance to pixel count, the pixel that n is phase place to be calculated is with respect to the folded pixel count of regional left margin along distance to skew of covering, n=0,1 ... N Pixel-1.
Step S118: according to each orientation upwards the folded solution of covering intra-zone twine phase place and first passage carrier aircraft height, first passage oblique distance and base length and baseline angle come inverting in each orientation to the DEM of each pixel.The schematic diagram that Fig. 4 is the interference SAR geometric relationship.In conjunction with Fig. 4, inverting is folded cover each orientation in zone to the computing method of DEM of each pixel as follows:
Figure BDA00001645658700071
θ = α + a sin ( R 1 2 + B 2 - R 2 2 2 R 1 B ) - - - ( 8 )
h＝H-R 1cos θ (9)
Wherein, R 1For the first passage oblique distance, the wavelength that λ is radar used, R 2For the second channel oblique distance, α is baseline angle, and B is base length, the visual angle that θ is first passage, and h is the DEM that is finally inversed by target, the height that H is first passage.
Below in conjunction with specific embodiment, and, with reference to accompanying drawing, the present invention is described in more detail.Confirmed the validity of the inventive method by emulation experiment.Fig. 5 is the elevation schematic diagram for the treatment of simulating area.The interferometric phase image that Fig. 6 simulates according to simulation parameter for the elevation schematic diagram by Fig. 5.The simulation parameter adopted is as shown in table 1:
Table 1 simulation parameter table
Wavelength Flying height Visual angle Base length Baseline angle
0.03m 3000m 53deg 2m 0.5deg
As shown in Figure 6, folded cover the stack of zone due to a plurality of scattering sources, interference fringe is unintelligible, but still that the interferometric phase that can find out regional B accounts for is leading, and it is regional contrary that the frequency of phase place and the left and right sides non-folded covered.Cover zone and processed folded by method of the present invention, finally obtain the DEM reconstructed results on oblique distance plane as shown in Figure 7, in order to verify that the inventive method covers to folded the correctness that regional DEM rebuilds, the actual oblique distance DEM obtained and desirable oblique distance DEM are contrasted.Fig. 8 is along a certain orientation upwards, according to the folded resulting DEM of zone digit elevation model method for reconstructing and the desirable DEM correlation curve figure of covering of embodiment of the present invention interference SAR.As shown in Figure 8, both are basically identical.The inventive method has been avoided the phase unwrapping error propagation that the folded phase hit of covering zone boundary causes as can be seen here, has recovered preferably the folded regional DEM that covers.
It should be noted that, the above-mentioned definition to each step is not limited in the various specific implementations of mentioning in embodiment, and those of ordinary skill in the art can replace simply to it with knowing, for example:
(1) detecting the folded zone of covering can also be judged according to the folded characteristics of covering regional interferometric phase;
(2) calculating the interferometric phase frequency can replace by methods such as maximal possibility estimation or MUSIC (Multi-Signal Classification) super-resolution estimations.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of making, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (10)

1. an interference SAR is folded and is covered the method that the zone digit elevation model is rebuild, and comprising:
The folded zone of covering of the first interference SAR complex pattern that detection is obtained by first passage and the second interference SAR complex pattern obtained by second channel;
The first interference SAR complex pattern and the second interference SAR complex pattern are carried out to conjugate multiplication, the result of conjugate multiplication is got to phase place, thereby obtain interferometric phase image;
To fold and cover zone as mask, and on interferometric phase image, except folded, cover the extra-regional non-folded zone of covering and carry out the interferometric phase solution and twine, obtain non-ly foldedly covering regional solution and twining phase place;
In each orientation upwards, by described, non-ly foldedly cover regional solution and twine phase place, obtain respectively folded non-folded the first solution of covering zone of covering regional left border and twine phase place
Figure FDA00001645658600011
Twine phase place with non-folded the 3rd solution of covering zone of right side boundary
Figure FDA00001645658600012
In each orientation to, non-ly folded cover regional solution and twine phase place by described, calculate folded cover regional left side along distance to front interferometric phase frequency f Before, and right side along the distance to rear interferometric phase frequency f After
Upwards, by the first solution, twine phase place in each orientation
Figure FDA00001645658600013
Twine phase place with the 3rd solution
Figure FDA00001645658600014
And front phase frequency f BeforeWith rear phase frequency f After, calculate this orientation to the folded solution of covering regional left border twine phase place
Figure FDA00001645658600015
Twine phase place with the solution of right side boundary
Figure FDA00001645658600016
In each orientation upwards, twine phase place by the folded solution of covering regional left border Twine phase place with the solution of right side boundary
Figure FDA00001645658600018
The solution of coming this orientation of matching upwards to fold each pixel of covering intra-zone twines phase place
Figure FDA00001645658600019
In each orientation to each pixel, twine phase place by the solution of this pixel Digital elevation model h with this pixel of interference SAR calculation of parameter.
2. interference SAR according to claim 1 is folded covers zone digit elevation model method for reconstructing, wherein, describedly by non-, foldedly covered regional solution and twines phase place, calculate folded cover regional left side along distance to front interferometric phase frequency f Before, and right side along the distance to rear interferometric phase frequency f AfterStep comprise:
Cover regional left side and the estimation window of (2M+1) * (2N+1) is respectively chosen on right side folded;
Corresponding estimate that the cost function of the phase model of window estimates front interferometric phase frequency f by maximizing BeforeWith rear interferometric phase frequency f After.
3. interference SAR according to claim 2 is folded covers zone digit elevation model method for reconstructing, wherein, and described M=16, N=16.
4. interference SAR according to claim 1 is folded covers zone digit elevation model method for reconstructing, wherein, according to following formula, calculate this orientation to the folded solution of covering regional left border twine phase place
Figure FDA00001645658600021
Twine phase place with the solution of right side boundary
Figure FDA00001645658600023
Figure FDA00001645658600024
Wherein, N PixelFor folded cover zone along distance to pixel count.
5. interference SAR according to claim 4 is folded covers zone digit elevation model method for reconstructing, wherein, according to following formula, by the folded solution of covering regional left border, twines phase place Twine phase place with the solution of right side boundary
Figure FDA00001645658600026
Come this orientation of matching upwards to fold the solution of covering intra-zone and twine phase place
Figure FDA00001645658600027
Figure FDA00001645658600028
Wherein, the pixel that n is phase place to be calculated is covered regional left margin along distance to the pixel count be offset with respect to folded, n=0, and 1 ... N Pixel-1.
6. interference SAR according to claim 1 is folded covers zone digit elevation model method for reconstructing, wherein, according to following formula, by the solution of pixel, twines phase place With the digital elevation model h of this pixel of interference SAR calculation of parameter:
Figure FDA000016456586000210
θ = α + a sin ( R 1 2 + B 2 - R 2 2 2 R 1 B )
h＝H-R 1cosθ
Wherein, R 1For the oblique distance of first passage, the wavelength that λ is radar used, R 2For the oblique distance of second channel, α is baseline angle, and B is base length, the visual angle that θ is first passage, and h is the DEM that is finally inversed by target, the height that H is first passage.
7. cover zone digit elevation model method for reconstructing according to the described interference SAR of any one in claim 1 to 6 is folded, wherein, described detection the first interference SAR complex pattern and the second interference SAR complex pattern folded also comprises before covering regional step:
The first interference SAR complex pattern and the second interference SAR complex pattern are carried out to registration.
8. interference SAR according to claim 7 is folded and is covered zone digit elevation model method for reconstructing, wherein, take following mode one of them, the first interference SAR complex pattern and the second interference SAR complex pattern are carried out to registration: maximum coherence Y-factor method Y or maximum spectrum registration method.
9. interference SAR according to claim 1 is folded covers zone digit elevation model method for reconstructing, and wherein, described detection the first interference SAR complex pattern and the second interference SAR complex pattern folded covered regional step and comprise:
Ask for the averaged amplitude value M of whole pixels of interference SAR complex pattern a
Set a global threshold K a, make K a=cM a, 2≤c≤5 wherein;
Amplitude in described interference SAR complex pattern is greater than to K aPixel as the folded zone of covering;
Wherein, described interference SAR complex pattern is the first interference SAR complex pattern or the second interference SAR complex pattern.
10. interference SAR according to claim 9 is folded covers zone digit elevation model method for reconstructing, wherein, described amplitude in the interference SAR complex pattern is greater than to K aPixel as folded, also comprise after covering regional step:
Set the coefficient of coherence threshold k c
Using described folded folded the cover zone of zone as preliminary judgement of covering;
Calculate the coefficient of coherence of every part in the coefficient of coherence figure of the first interference SAR complex pattern and the second interference SAR complex pattern;
In preliminary judgement folded covered zone by coefficient of coherence higher than K cPart remove, thereby obtain the final folded zone of covering.
CN201210152477.2A 2012-05-16 2012-05-16 Interference SAR layover area digital elevation model rebuilding method Active CN103424744B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210152477.2A CN103424744B (en) 2012-05-16 2012-05-16 Interference SAR layover area digital elevation model rebuilding method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210152477.2A CN103424744B (en) 2012-05-16 2012-05-16 Interference SAR layover area digital elevation model rebuilding method

Publications (2)

Publication Number Publication Date
CN103424744A true CN103424744A (en) 2013-12-04
CN103424744B CN103424744B (en) 2015-03-25

Family

ID=49649763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210152477.2A Active CN103424744B (en) 2012-05-16 2012-05-16 Interference SAR layover area digital elevation model rebuilding method

Country Status (1)

Country Link
CN (1) CN103424744B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249234A (en) * 2016-07-08 2016-12-21 中国科学院电子学研究所 A kind of InSAR water body region interferometric phase unwrapping method
CN108132468A (en) * 2017-12-25 2018-06-08 中南大学 A kind of more baseline polarimetric SAR interferometry depth of building extracting methods
CN108761458A (en) * 2018-08-15 2018-11-06 中国科学院电子学研究所 Interference SAR water body digital elevation model modification method based on morphologic thinning
CN108957454A (en) * 2018-08-01 2018-12-07 中国科学院电子学研究所 The interferometric phase emulation mode of airborne Interference synthetic aperture radar
CN109212522A (en) * 2018-05-28 2019-01-15 中国科学院电子学研究所 A kind of method and apparatus obtaining numerical map
CN111696207A (en) * 2020-05-28 2020-09-22 电子科技大学 Multi-baseline DEM fusion method based on guided filtering
CN113640758A (en) * 2021-08-23 2021-11-12 中国科学院空天信息创新研究院 SAR image scaler placement method and system under urban complex environment
WO2021227423A1 (en) * 2020-05-13 2021-11-18 深圳大学 Insar digital elevation model construction method and system based on dynamic baseline

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369019A (en) * 2008-10-10 2009-02-18 清华大学 Polarization interference synthetic aperture radar three-dimensional imaging method based on polarization data amalgamation
CN101551455A (en) * 2009-05-13 2009-10-07 西安电子科技大学 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369019A (en) * 2008-10-10 2009-02-18 清华大学 Polarization interference synthetic aperture radar three-dimensional imaging method based on polarization data amalgamation
CN101551455A (en) * 2009-05-13 2009-10-07 西安电子科技大学 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱俊杰 等: "基于高分辨率SAR图像的建筑物高度提取", 《现代雷达》 *
李芳芳 等: "一种基于质量指导的InSAR相位解缠快速实现方法", 《雷达学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249234A (en) * 2016-07-08 2016-12-21 中国科学院电子学研究所 A kind of InSAR water body region interferometric phase unwrapping method
CN106249234B (en) * 2016-07-08 2019-06-04 中国科学院电子学研究所 A kind of InSAR water area interferometric phase unwrapping method
CN108132468A (en) * 2017-12-25 2018-06-08 中南大学 A kind of more baseline polarimetric SAR interferometry depth of building extracting methods
CN109212522A (en) * 2018-05-28 2019-01-15 中国科学院电子学研究所 A kind of method and apparatus obtaining numerical map
CN108957454A (en) * 2018-08-01 2018-12-07 中国科学院电子学研究所 The interferometric phase emulation mode of airborne Interference synthetic aperture radar
CN108761458A (en) * 2018-08-15 2018-11-06 中国科学院电子学研究所 Interference SAR water body digital elevation model modification method based on morphologic thinning
CN108761458B (en) * 2018-08-15 2021-06-29 中国科学院电子学研究所 Morphological refinement-based interference SAR water body digital elevation model correction method
WO2021227423A1 (en) * 2020-05-13 2021-11-18 深圳大学 Insar digital elevation model construction method and system based on dynamic baseline
CN111696207A (en) * 2020-05-28 2020-09-22 电子科技大学 Multi-baseline DEM fusion method based on guided filtering
CN111696207B (en) * 2020-05-28 2022-10-11 电子科技大学 Multi-baseline DEM fusion method based on guided filtering
CN113640758A (en) * 2021-08-23 2021-11-12 中国科学院空天信息创新研究院 SAR image scaler placement method and system under urban complex environment

Also Published As

Publication number Publication date
CN103424744B (en) 2015-03-25

Similar Documents

Publication Publication Date Title
CN103424744A (en) Interference SAR layover area digital elevation model rebuilding method
CN104459633A (en) Wavelet domain InSAR interferometric phase filtering method combined with local frequency estimation
CN104111456B (en) A kind of line of high-speed railway Ground Deformation high-resolution InSAR monitoring methods
CN103217674B (en) A kind of target three-dimensional scattering center reconstructing method of inverse synthetic aperture radar (ISAR)
CN103454636B (en) Differential interferometric phase estimation method based on multi-pixel covariance matrixes
CN101650439B (en) Method for detecting change of remote sensing image based on difference edge and joint probability consistency
CN102999759A (en) Light stream based vehicle motion state estimating method
CN106772342A (en) A kind of Timing Difference radar interference method suitable for big gradient surface subsidence monitoring
CN103576149B (en) A kind of foundation interference radar three-dimensional deformation extraction method based on amplitude information
CN102629380B (en) Remote sensing image change detection method based on multi-group filtering and dimension reduction
CN104007439B (en) Interferential circular SAR elevation estimation processing method
CN103439708B (en) Polarized InSAR interferogram estimation method based on generalized scattering vector
CN102645652A (en) Method for detecting single channel synthetic aperture radar (SAR) ground motion target
CN103630900B (en) The method of 3-D SAR wavenumber domain fast imaging
CN103616686B (en) A kind of optimum angle method of estimation of the complete polarization interference synthetic aperture radar image based on mixed mode
CN103823219B (en) The non local interfering synthetic aperture radar interferometric phase filtering method of adaptive iteration
CN102024268B (en) Model simplification method for keeping appearance characteristics
CN103777188A (en) InISAR imaging method based on double-frequency conjugation processing
CN102800113B (en) Digital image analysis method based on fractal dimension
CN104515980A (en) Method and device for ground moving target indication based on InSAR (interferometric synthetic aperture radar) formation
CN105975957A (en) Lane-line-edge-based road plane detection method
CN103809180A (en) Azimuth pre-filtering processing method for Interferometric Synthetic Aperture Radar (InSAR) topographic survey
CN102903078B (en) A kind of motion blur image method for parameter estimation based on multiresolution Fourier analysis theory
CN104200458A (en) MeanShift based high-resolution remote sensing image segmentation distance measurement optimization method
CN103226194A (en) Interferometric synthetic aperture radar (InSAR) interferometric phase filtering method based on empirical mode decomposition

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