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 PDF

Info

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
Application number
CN201710672671.6A
Other languages
Chinese (zh)
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.)
SURVEYING MAPPING TECHNOLOGY SERVICE CENTER SICHUAN BUREAU OF SURVEYING MAPPING AND GEOINFORMATION
Original Assignee
SURVEYING MAPPING TECHNOLOGY SERVICE CENTER SICHUAN BUREAU OF SURVEYING MAPPING AND GEOINFORMATION
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 SURVEYING MAPPING TECHNOLOGY SERVICE CENTER SICHUAN BUREAU OF SURVEYING MAPPING AND GEOINFORMATION filed Critical SURVEYING MAPPING TECHNOLOGY SERVICE CENTER SICHUAN BUREAU OF SURVEYING MAPPING AND GEOINFORMATION
Priority to CN201710672671.6A priority Critical patent/CN107422324A/en
Publication of CN107422324A publication Critical patent/CN107422324A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9076Polarimetric features in SAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means 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

Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information
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
CN201710672671.6A 2017-08-08 2017-08-08 Reflect the synthetic aperture radar polarimetric calibration method of asymmetrical information Pending CN107422324A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
廖露: "机载SAR系统极化定标方法的优化研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
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