CN107422324A - Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information - Google Patents
Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information Download PDFInfo
- Publication number
- CN107422324A CN107422324A CN201710672671.6A CN201710672671A CN107422324A CN 107422324 A CN107422324 A CN 107422324A CN 201710672671 A CN201710672671 A CN 201710672671A CN 107422324 A CN107422324 A CN 107422324A
- Authority
- CN
- China
- Prior art keywords
- matrix
- polarization
- synthetic aperture
- aperture radar
- formula
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- 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/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
-
- 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/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information, methods described comprises the following steps:Obtain non-nominal data matrix [O];Non- nominal data is switched into covariance matrix, initialization phasor U, V, W, Z are arranged to 0,It is arranged to 1;Using initialization phasor U, V, W, Z,, adjustment variance matrix [C];Polarization orientation angle excursion matrix P corresponding to pattern-band is got by two waveband registration;According to excursion matrix P, U, V, W, Z are calculated;Calculate again;It is all only to meet that hypsography is flat and flight attitude stable condition to solve traditional polarimetric calibration method so that polarimetric calibration can not solve distributed object and be unsatisfactory for reflection symmetry to assume problem.
Description
Technical field
The present invention relates to SAR polarization radar calibrating methods, and in particular to reflects the synthetic aperture radar pole of asymmetrical information
Change calibrating method.
Background technology
It is related to for the processing method few people for calibrating front-reflection asymmetrical information, Kimura is excluded anti-using performance number
Remaining using area participates in calibration (Kimuraetal., 2004) after penetrating asymmetrical construction zone, and Freeman assumes uncertain
Systematic error before mark is relatively low to propose a kind of method of estimation of Faraday rotation (Freeman, 2004) with this.Domestic scholars have
Flight attitude angle is substituted into the trial (Zhang Lintao etc., 2011) of Whitt algorithms.Under normal conditions, the flight angle of pitch will not surpass
10 ° are crossed, and ground line gradient is many times, many regions all run far deeper than this degree, so airborne platform calibration needs are examined
It is more introducing and the parametric solution method that asymmetrical information is reflected under the influence of orographic condition to consider.But if landform is flat,
And flight attitude angle is when changing greatly, the introducing reflection asymmetrical information of this patent introduction is determined calibration method and is readily applicable to
This.
Traditional polarimetric calibration method is all only to meet that hypsography is flat and flight attitude stable condition so that polarization is fixed
Mark can not solve the hypothesis that distributed object is unsatisfactory for reflection symmetry.Using multi-source data, hypsography can be provided in advance
Region and the information of flight attitude unstable (polarization orientation angle).Using this prior information, propose that a kind of reflection is asymmetric known
In the case of introduce its information, the polarimetric calibration method unstable available for hypsography region and flight attitude.
The content of the invention
It is an object of the invention to provide one kind, and it is all only to meet that hypsography is flat to solve traditional polarimetric calibration method
With flight attitude stable condition so that polarimetric calibration can not solve distributed object and be unsatisfactory for the hypothesis of reflection symmetry to ask
Topic.
To solve above-mentioned technical problem, the present invention uses following technical scheme:Reflect the synthetic aperture of asymmetrical information
Polarization radar calibrating method, methods described comprise the following steps:
Step 1:Obtain non-nominal data matrix [O];
Step 2:Non- nominal data is switched into covariance matrix, initialization phasor U, V, W, Z are arranged to 0,It is arranged to 1;
Step 3:Using initialization phasor U, V, W, Z,Adjust variance matrix [C];
Step 4:Polarization orientation angle excursion matrix P corresponding to pattern-band is got by two waveband registration;
Step 5:According to excursion matrix P, U, V, W, Z are calculated;
Step 6:Calculate again
Step 7:Using corner reflector, judgeWhether restrain, if convergence, remove same polarization channel imbalance, perform step
Rapid 8;Initialization phasor U, V, W, Z are used if not restraining,Variance matrix [C] is adjusted, P is got by two waveband registration
Polarization orientation angle excursion matrix P corresponding to wave band;According to excursion matrix P, U, V, W, Z are calculated;Calculate againRepeat this step
Suddenly;
Step 8:Same polarization channel imbalance is removed, using U, V, W, Z,K corrects data.
Further technical scheme is that the data matrix of synthetic aperture radar described in step 1 includes synthetic aperture radar
VV passages echo data, the VH passages echo data of synthetic aperture radar, synthetic aperture radar HH passages echo data,
The HV passage echo datas of synthetic aperture radar.
Further technical scheme is, in step 3, the variance matrix [C]=[O] [O]H。
Further technical scheme is that in step 4, the two waveband registration gets polarization side corresponding to pattern-band
Parallactic angle excursion matrix P sub-step includes:
4a:DEM is generated using the interference data of single track double antenna X-band;
4b:FormulaPolarization orientation angle offset is calculated.The orientation gradient
Calculated with distance to the gradient by dem data.Wherein, incidence angle information comes from metadata;
4c:Main photography of the pattern-band HH polarization amplitudes image as registering image centering.Correspondingly, X-band HH polarization amplitudes
Image is considered as registering subpictures (target image to be matched), it should be pointed out that X-band HH polarization amplitudes image come from double antenna
Primary antenna part in interference data.By artificial reconnaissance and matching obtain two assemble to ground control point collection.Then use
Secondary geometrical model obtains the variation relation of two width images.
4d:Because dem data have with X-band primary antenna image as coordinate system, polarization orientation angle skew and X
Transformation relation between wave band image.After registration parameter conversion and resampling, polarised direction angular variation can obtain.
Further technical scheme is polarised direction angle offset and radar visual angle and the geometrical relationship of terrain slope
It is as follows
In formula, tan ω are the direction gradients, and tan γ are distances to the gradient, φ radars visual angle.
Further technical scheme is, if described in step 8Convergence, same polarization channel imbalance is removed, using U,
V, W, Z,K corrects data, specifically includes following sub-step:
6a:C=OOH=YMPS (SPMY)H=(YXFKP) E (PKFXY)H
In formulaIn order to build
Calibration process scaling factor Y is expressed to obtain.It is expressed as below:
E∝(PM)-1C((MP)H)-1=(PKFX)-1C((XFKP)H)-1
In formula, H represents conjugate transposition, and X represents crosstalk matrix, and E represents preferable covariance matrix.F and K is to intersect respectively
Polarization and the distortion matrix of same polarization passage.
6b:Assume preferable covariance matrix be:
* is conjugate of symbol in formula, and G represents that HV and VH is equal, and 80 values represent reflection symmetry in matrix E.
6c:Crosstalk and the unbalanced estimation of cross polarization
6d:Same polarization difference between diversity channels degree is estimated;
After crosstalk and the degree of balance are removed, reflector is laid in scene to estimate same polarization degree of unbalancedness.
Further technical scheme is that step 5 calculates U, V, W, Z according to excursion matrix P;Including following sub-step:
5a:According to basic model VanZyl, 1990, reflecting asymmetrical basic polarization distortion model has:
In O=RST+N formulas,
Above formula can be equivalent to by O, S, R, T, N definition
5b:Z in formula is represented in the case of no rotation, equal to S in formula.Only rotated, but not by other
The polarization scattering matrix that factor influences.If collision matrix belongs to reflection asymmetric target, then it with preferable collision matrix it
Between just have the relation of the individual anglec of rotation, can be expressed as following formula:
5c:Now, Z is expressed as rotating the polarization scattering matrix after θ angles.What θ was represented is polarization orientation angle skew joint
Above formula
5d:Two collision matrixes (O, S) are converted into vector form by Quegan and Ainsworth.Above formula can be changed to
5e:Above formula can be expressed as:
O=YMPS
Wherein, u, v, w, z are four crosstalks, k andIt is same polarization and cross-pole channel imbalance respectively
Degree, Y is the absolute scatter factor, and M is a day line distortion matrix, and P is polarization orientation angle spin matrix;In formula
Compared with prior art, the beneficial effects of the invention are as follows:
Excursion matrix P is added, reflection asymmetrical information is observed object on radar plane of incidence minute surface symmetric points
Offset information, that is, the information of hypsography region and flight attitude unstable (polarization orientation angle) is calculated, so as to landform
The synthetic aperture radar of undulating region and flight attitude unstable (polarization orientation angle) carries out polarimetric calibration.
Brief description of the drawings
Fig. 1 is the synthetic aperture radar polarimetric calibration method geometrical relationship figure of present invention reflection asymmetrical information.
Fig. 2 is that asymmetric letter is reflected in the introducing of the synthetic aperture radar polarimetric calibration method of present invention reflection asymmetrical information
The polarimetric calibration techniqueflow of breath.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
Embodiment 1
Fig. 1 to Fig. 2 is shown:The synthetic aperture radar polarimetric calibration method of asymmetrical information is reflected, methods described includes
Following steps:
Step 1:Obtain non-nominal data matrix [O];
Step 2:Non- nominal data is switched into covariance matrix, initialization phasor U, V, W, Z are arranged to 0,It is arranged to 1;
Step 3:Using initialization phasor U, V, W, Z,Adjust variance matrix [C];
Step 4:Polarization orientation angle excursion matrix P corresponding to pattern-band is got by two waveband registration;
Step 5:According to excursion matrix P, U, V, W, Z are calculated;
Step 6:Calculate again
Step 7:Using corner reflector, judgeWhether restrain, if convergence, remove same polarization channel imbalance, perform step
Rapid 8;Initialization phasor U, V, W, Z are used if not restraining,Variance matrix [C] is adjusted, P is got by two waveband registration
Polarization orientation angle excursion matrix P corresponding to wave band;According to excursion matrix P, U, V, W, Z are calculated;Calculate againRepeat this step
Suddenly;
Step 8:Same polarization channel imbalance is removed, using U, V, W, Z,K corrects data.
The present embodiment, polarization orientation angle corresponding to pattern-band is got by two waveband registration and offset.Then, using point
The reciprocity and reflection symmetry of cloth target are assumed by way of iteration, try to achieve crosstalk and cross aisle degree of unbalancedness.Most
Eventually, the scaling parameter obtained is calibrated using the last corner reflector and gets rid of same polarization channel imbalance.
Embodiment 2
On the basis of embodiment 1, the data matrix of synthetic aperture radar described in step 1 includes synthetic aperture radar
VV passages echo data, the VH passages echo data of synthetic aperture radar, the HH passages echo data of synthetic aperture radar, synthesis
The HV passage echo datas of aperture radar.
The present embodiment, describes the composition of the matrix parameter of synthetic aperture radar in detail, i.e. O phasors areZ phasors
For
Embodiment 3
On the basis of embodiment 1, in step 3, the variance matrix [C]=[O] [O]H。
The present embodiment, detail variance matrix [C] solution procedure.
Embodiment 4
On the basis of embodiment 1, in step 4, the two waveband registration gets polarization orientation corresponding to pattern-band
Angular variation matrix P sub-step includes:
4a:DEM is generated using the interference data of single track double antenna X-band;
4b:FormulaPolarization orientation angle offset is calculated.The orientation gradient
Calculated with distance to the gradient by dem data.Wherein, incidence angle information comes from metadata;
4c:Main photography of the pattern-band HH polarization amplitudes image as registering image centering.Correspondingly, X-band HH polarization amplitudes
Image is considered as registering subpictures (target image to be matched), it should be pointed out that X-band HH polarization amplitudes image come from double antenna
Primary antenna part in interference data.By artificial reconnaissance and matching obtain two assemble to ground control point collection.Then use
Secondary geometrical model obtains the variation relation of two width images.
4d:Because dem data have with X-band primary antenna image as coordinate system, polarization orientation angle skew and X
Transformation relation between wave band image.After registration parameter conversion and resampling, polarised direction angular variation can obtain.
The present embodiment, details excursion matrix P solution procedure, and traditional polarimetric calibration method is all only to meet landform
Rise and fall flat and flight attitude stable condition so that polarimetric calibration can not solve the vacation that distributed object is unsatisfactory for reflection symmetry
If.Using multi-source data, excursion matrix P is obtained, hypsography region and unstable (the polarization orientation of flight attitude can be provided in advance
Angle) information.
Embodiment 5
On the basis of embodiment 3, polarised direction angle offset and radar visual angle and the geometrical relationship of terrain slope
It is as follows
In formula, tan ω are the direction gradients, and tan γ are distances to the gradient, φ radars visual angle.
The present embodiment calculates, terrain slope angle.
Embodiment 6
On the basis of embodiment 1, if described in step 8Convergence, same polarization channel imbalance is removed, using U,
V, W, Z,K corrects data, specifically includes following sub-step:
6a:C=OOH=YMPS (SPMY)H=(YXFKP) E (PKFXY)H
In formulaIn order to build
Calibration process scaling factor Y is expressed to obtain.It is expressed as below:
E∝(PM)-1C((MP)H)-1=(PKFX)-1C((XFKP)H)-1
In formula, H represents conjugate transposition, and X represents crosstalk matrix, and E represents preferable covariance matrix.F and K is to intersect respectively
Polarization and the distortion matrix of same polarization passage.
6b:Assume preferable covariance matrix be:
* is conjugate of symbol in formula, and G represents that HV and VH is equal, and 80 values represent reflection symmetry in matrix E.
6c:Crosstalk and the unbalanced estimation of cross polarization
6d:Same polarization difference between diversity channels degree is estimated;
After crosstalk and the degree of balance are removed, reflector is laid in scene to estimate same polarization degree of unbalancedness.
The present embodiment, estimate crosstalk and degree of unbalancedness.
Embodiment 7
On the basis of any one described in embodiment 1 to 6, step 5 calculates U, V, W according to excursion matrix P,
Z;Including following sub-step:
5a:According to basic model VanZyl, 1990, reflecting asymmetrical basic polarization distortion model has:
In O=RST+N formulas,
Above formula can be equivalent to by O, S, R, T, N definition
5b:Z in formula is represented in the case of no rotation, equal to S in formula.Only rotated, but not by other
The polarization scattering matrix that factor influences.If collision matrix belongs to reflection asymmetric target, then it with preferable collision matrix it
Between just have the relation of the individual anglec of rotation, can be expressed as following formula:
5c:Now, Z is expressed as rotating the polarization scattering matrix after θ angles.What θ was represented is polarization orientation angle skew joint
Above formula
5d:Two collision matrixes (O, S) are converted into vector form by Quegan and Ainsworth.Above formula can be changed to
5e:Above formula can be expressed as:
O=YMPS
Wherein, u, v, w, z are four crosstalks, k andIt is same polarization and cross-pole channel imbalance degree respectively,
Y is the absolute scatter factor, and M is a day line distortion matrix, and P is polarization orientation angle spin matrix;In formula
The present embodiment, detail and solve u, v, w, z calculating process.
" one embodiment ", " another embodiment ", " embodiment ", " preferred embodiment " spoken of in this manual
Deng, refer to combine the embodiment description specific features, structure or feature be included in the application generality description at least
In one embodiment.It is not necessarily to refer to same embodiment that statement of the same race, which occur, in multiple places in the description.Further
For, when describing a specific features, structure or feature with reference to any embodiment, what is advocated is to combine other embodiment
To realize that this feature, structure or feature are also fallen within the scope of the present invention.
Although reference be made herein to invention has been described for multiple explanatory embodiments of the invention, however, it is to be understood that
Those skilled in the art can be designed that a lot of other modifications and embodiment, and these modifications and embodiment will fall in this Shen
Please be within disclosed spirit and spirit.More specifically, can in the range of disclosure, drawings and claims
A variety of variations and modifications are carried out with the building block to theme combination layout and/or layout.Except to building block and/or layout
Outside the modification and improvement of progress, to those skilled in the art, other purposes also will be apparent.
Claims (7)
1. reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information, it is characterised in that:Methods described includes following step
Suddenly:
Step 1:Obtain non-nominal data matrix [O];
Step 2:Non- nominal data is switched into covariance matrix, initialization phasor U, V, W, Z are arranged to 0,It is arranged to 1;
Step 3:Using initialization phasor U, V, W, Z,, adjustment variance matrix [C];
Step 4:Polarization orientation angle excursion matrix P corresponding to pattern-band is got by two waveband registration;
Step 5:According to excursion matrix P, U, V, W, Z are calculated;
Step 6:Calculate again;
Step 7:Using corner reflector, judgeWhether restrain, if convergence, remove same polarization channel imbalance, perform step 8;
Initialization phasor U, V, W, Z are used if not restraining,, adjustment variance matrix [C], pattern-band is got by two waveband registration
Corresponding polarization orientation angle excursion matrix P;According to excursion matrix P, U, V, W, Z are calculated;Calculate again, repeat this step;
Step 8:Same polarization channel imbalance is removed, using U, V, W, Z,, k correction data.
2. the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information according to claim 1, it is characterised in that:
The data matrix of synthetic aperture radar described in step 1 includes the VV passages echo data of synthetic aperture radar, synthetic aperture radar
VH passages echo data, HH passages echo data, the HV passage echo datas of synthetic aperture radar of synthetic aperture radar.
3. the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information according to claim 1, it is characterised in that:
In step 3, the variance matrix [C]=[O] [O]H。
4. the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information according to claim 1, it is characterised in that:
In step 4, the two waveband registration gets the sub-step of polarization orientation angle excursion matrix P corresponding to pattern-band and included:
4a:DEM is generated using the interference data of single track double antenna X-band;
4b:FormulaPolarization orientation angle offset is calculated.
The orientation gradient and distance are calculated to the gradient by dem data.
Wherein, incidence angle information comes from metadata;
4c:Main photography of the pattern-band HH polarization amplitudes image as registering image centering.
Correspondingly, X-band HH polarization amplitudes image is considered as registering subpictures(Target image to be matched), it should be pointed out that X ripples
Primary antenna part of the section HH polarization amplitudes image in double antenna interference data.
By artificial reconnaissance and matching obtain two assemble to ground control point collection.
The variation relation of two width images is then obtained using secondary geometrical model.
4d:Because dem data have with X-band primary antenna image as coordinate system, polarization orientation angle skew and X-band
Transformation relation between image.
After registration parameter conversion and resampling, polarised direction angular variation can obtain.
5. the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information according to claim 3, it is characterised in that:
Polarised direction angle offset and the geometrical relationship at radar visual angle and terrain slope are as follows
In formula,It is the direction gradient,It is distance to the gradient,Radar visual angle.
6. the synthetic aperture radar polarimetric calibration method of reflection asymmetrical information according to claim 1, it is characterised in that:
If described in step 8Convergence, same polarization channel imbalance is removed, using U, V, W, Z,, k corrects data, specific to wrap
Include following sub-step:
6a:C=OOH=YMPS(SPMY)H=(YXFKP)E(PKFXY)H
In formula,,In order to which construction expression is calibrated
Process scaling factor Y is obtained.
It is expressed as below:
In formula, H represents conjugate transposition, and X represents crosstalk matrix, and E represents preferable covariance matrix.
F and K is the distortion matrix of cross polarization and same polarization passage respectively.
6b:Assume preferable covariance matrix be:
* is conjugate of symbol in formula, and G represents that HV and VH is equal, and 80 values represent reflection symmetry in matrix E.
6c:Crosstalk and the unbalanced estimation of cross polarization
;
6d:Same polarization difference between diversity channels degree is estimated;
After crosstalk and the degree of balance are removed, reflector is laid in scene to estimate same polarization degree of unbalancedness.
7. the synthetic aperture thunder of the reflection asymmetrical information of any one in the claim according to claim 1 to 5
Up to polarimetric calibration method, it is characterised in that:Step 5 calculates U, V, W, Z according to excursion matrix P;Including following sub-step:
5a:According to basic model Van Zyl, 1990, reflecting asymmetrical basic polarization distortion model has:
In O=RST+N formulas,
Above formula can be equivalent to by O, S, R, T, N definition
5b:Z in formula is represented in the case of no rotation, equal to S in formula.
Only rotated, but the polarization scattering matrix not influenceed by other factors.
If collision matrix belongs to reflection asymmetric target, then it just has the pass of the individual anglec of rotation between preferable collision matrix
System, can be expressed as following formula:
5c:Now, Z is expressed as rotating the polarization scattering matrix after θ angles.
What θ was represented is polarization orientation angle skew joint above formula
=
5d:By Quegan and Ainsworth by two collision matrixes(O, S)Be converted to vector form.
Above formula can be changed to
5e:Above formula can be expressed as:
Wherein, u, v, w, z are four crosstalks, k andIt is same polarization and cross-pole channel imbalance degree respectively, Y is
The absolute scatter factor, M are a day line distortion matrixes, and P is polarization orientation angle spin matrix;In formula。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710672671.6A CN107422324A (en) | 2017-08-08 | 2017-08-08 | Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710672671.6A CN107422324A (en) | 2017-08-08 | 2017-08-08 | Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107422324A true CN107422324A (en) | 2017-12-01 |
Family
ID=60437560
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710672671.6A Pending CN107422324A (en) | 2017-08-08 | 2017-08-08 | Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107422324A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375770A (en) * | 2018-03-15 | 2018-08-07 | 中国科学院遥感与数字地球研究所 | A kind of polarimetric calibration method of the full-polarization SAR data based on rotational symmetry |
CN110261853A (en) * | 2019-06-12 | 2019-09-20 | 中国测绘科学研究院 | A kind of calibration processing method and processing device of polarimetric synthetic aperture radar |
CN113960545A (en) * | 2021-10-19 | 2022-01-21 | 安徽大学 | Satellite-borne SAR field-free geometric calibration method and system based on symmetric geometric configuration constraint |
CN115047416A (en) * | 2022-08-15 | 2022-09-13 | 中国科学院空天信息创新研究院 | Full-polarization SAR calibration method based on symmetric region surface distortion component |
CN115291218A (en) * | 2022-10-10 | 2022-11-04 | 中国电子科技集团公司第十四研究所 | Homologous common-view multiband interference SAR test system |
CN115718284A (en) * | 2023-01-10 | 2023-02-28 | 中国科学院空天信息创新研究院 | Common polarization channel unbalance calibration method based on circular polarization reflection symmetry |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393513A (en) * | 2011-08-31 | 2012-03-28 | 中国电子科技集团公司第三十八研究所 | Polarimetric calibration technique based on natural distribution scenes and rare calibrator |
CN105445708A (en) * | 2015-11-11 | 2016-03-30 | 西安电子科技大学 | Calibration method for polarization synthetic-aperture radar |
-
2017
- 2017-08-08 CN CN201710672671.6A patent/CN107422324A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393513A (en) * | 2011-08-31 | 2012-03-28 | 中国电子科技集团公司第三十八研究所 | Polarimetric calibration technique based on natural distribution scenes and rare calibrator |
CN105445708A (en) * | 2015-11-11 | 2016-03-30 | 西安电子科技大学 | Calibration method for polarization synthetic-aperture radar |
Non-Patent Citations (1)
Title |
---|
廖露: "机载SAR系统极化定标方法的优化研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375770A (en) * | 2018-03-15 | 2018-08-07 | 中国科学院遥感与数字地球研究所 | A kind of polarimetric calibration method of the full-polarization SAR data based on rotational symmetry |
CN108375770B (en) * | 2018-03-15 | 2020-07-17 | 中国科学院遥感与数字地球研究所 | Polarization calibration method of full-polarization SAR data based on rotational symmetry |
CN110261853A (en) * | 2019-06-12 | 2019-09-20 | 中国测绘科学研究院 | A kind of calibration processing method and processing device of polarimetric synthetic aperture radar |
CN110261853B (en) * | 2019-06-12 | 2021-01-05 | 中国测绘科学研究院 | Calibration processing method and device for polarized synthetic aperture radar |
CN113960545A (en) * | 2021-10-19 | 2022-01-21 | 安徽大学 | Satellite-borne SAR field-free geometric calibration method and system based on symmetric geometric configuration constraint |
CN115047416A (en) * | 2022-08-15 | 2022-09-13 | 中国科学院空天信息创新研究院 | Full-polarization SAR calibration method based on symmetric region surface distortion component |
CN115047416B (en) * | 2022-08-15 | 2022-10-25 | 中国科学院空天信息创新研究院 | Full-polarization SAR calibration method based on symmetric region surface distortion component |
CN115291218A (en) * | 2022-10-10 | 2022-11-04 | 中国电子科技集团公司第十四研究所 | Homologous common-view multiband interference SAR test system |
CN115291218B (en) * | 2022-10-10 | 2022-12-09 | 中国电子科技集团公司第十四研究所 | Homologous common-view multiband interference SAR test system |
CN115718284A (en) * | 2023-01-10 | 2023-02-28 | 中国科学院空天信息创新研究院 | Common polarization channel unbalance calibration method based on circular polarization reflection symmetry |
CN115718284B (en) * | 2023-01-10 | 2023-04-04 | 中国科学院空天信息创新研究院 | Common polarization channel unbalance calibration method based on circular polarization reflection symmetry |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107422324A (en) | Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information | |
CN106405547B (en) | A kind of multipolarization SAR condenses polarization mode calibrating method | |
CN105445708B (en) | A kind of calibrating method of polarimetric synthetic aperture radar | |
CN109633576A (en) | A kind of polarization SAR high accuracy Scaling algorithm under channel noise | |
CN108375770B (en) | Polarization calibration method of full-polarization SAR data based on rotational symmetry | |
DE102009046499A1 (en) | radar device | |
CN107607948B (en) | Motion compensation and imaging method of medium orbit satellite SAR | |
CN109975755B (en) | Short-wave multi-station direct positioning method under condition of existence of correction source | |
CN107689064A (en) | Take the strict geometry imaging model construction method of satellite optical of aberration correction into account | |
CN105158759B (en) | HRWS SAR channel phases deviation correction methods based on clutter cancellation | |
CN104931958A (en) | Double-difference beam angle measurement method and system | |
CN109143237A (en) | The PFA inclination of wave front bearing calibration of biradical Spotlight SAR Imaging suitable for any platform track | |
CN103645466B (en) | Polarization calibration method based on platform attitude time variation compensation | |
CN115061156A (en) | Array antenna satellite navigation deception resisting method and system based on integrated navigation | |
CN108427104A (en) | A kind of complete polarization multiple-input and multiple-output synthetic aperture radar(FP-MIMO-SAR)Polarimetric Calibration method | |
CN103728601A (en) | Radar signal motion disturbance spatial-polarizational domain combined stable filtering method | |
CN107576951A (en) | Wave arrival direction estimating method based on nested type Electromagnetic Vector Sensor Array | |
CN105652255A (en) | Spatial aligning method for radar networking system | |
CN102135613B (en) | Antenna cross polarization correction method for satellite-borne polarized related completely-polarized microwave radiometer | |
CN110161496B (en) | Target surface normal vector obtaining method and system based on polarized millimeter wave radiation | |
CN101900805B (en) | Spherical wave imaging mathematical model and compensation method of near-field effect | |
CN108107427A (en) | Airborne/missile-borne array radar forword-looking imaging method based on super resolution technology | |
Deng et al. | A processing framework for airborne microwave photonic SAR with resolution up to 0.03 m: motion estimation and compensation | |
CN108776347B (en) | High-dynamic GNSS interference suppression method of dual-polarized antenna array based on null broadening technology | |
CN109696651A (en) | It is a kind of based on M estimation low number of snapshots under Wave arrival direction estimating method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171201 |
|
RJ01 | Rejection of invention patent application after publication |