CN104459655A - 混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 - Google Patents
混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 Download PDFInfo
- Publication number
- CN104459655A CN104459655A CN201410743634.6A CN201410743634A CN104459655A CN 104459655 A CN104459655 A CN 104459655A CN 201410743634 A CN201410743634 A CN 201410743634A CN 104459655 A CN104459655 A CN 104459655A
- Authority
- CN
- China
- Prior art keywords
- error
- transmitting terminal
- max
- phi
- omega
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
- G01S7/4008—Means for monitoring or calibrating of parts of a radar system of transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Communication System (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种混合极化SAR系统π/4简缩极化模式的发射端误差求解的方法。该方法首先推导出最大归一化误差表达式,而后将误差源代入该最大归一化误差表达式中,求取该误差源带来的误差。本发明可获得发射端各误差:发射波Jones矢量误差、通道失衡、通道串扰和/或法拉第旋转角对混合极化SAR系统的影响。
Description
技术领域
本发明涉及简缩极化雷达遥感领域,尤其涉及一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法。
背景技术
近年来随着雷达技术的快速发展,雷达已广泛应用于航天控制,军事探测,目标检测,林业遥感,环境监督等众多领域。为了更好地实现雷达遥感应用与目标探测,全极化雷达应运而生且在过去的60年来蓬勃发展。
但目前全极化雷达在应用上仍面临若干缺点:(1)同极化与交叉极化最优动态范围差异较大,需要复杂的通道增益控制和内定标回路控制;(2)较强的同极化距离模糊回波及较弱的交叉极化距离模糊回波需要控制入射角范围及观测带宽;(3)PRF是单极化系统的2倍;(4)全极化4个通道加重数据存储及传输。
为此,混合极化SAR系统应运而生,混合极化SAR系统具有三个大优势,1)其可实现简缩极化模式还可实现全极化模式;2)据有较好的抑制距离模糊的特点;3)通道间具有较高均衡性降低系统的复杂性,且使定标变得相对简单。混合极化SAR系统的三个简缩极化模式有:π/4,CL,CC简缩极化模式。其各模式特点如下。π/4简缩极化模式特点为发射1路由线性水平极化(H)与垂直极化(V)合成的线性极化波,同时接收H、V两路极化波。CL简缩极化模式为发射一路圆极化波,同时接收H、V两路极化波。此外还有发射圆极化波接收2路正交圆极化波的CC简缩极化模式。
混合极化技术是近年来的一个热点研究。2014年,加拿大雷达星座计划(Radar Constellation Mission,RCM)开始采用混合极化模式。混合极化系统设计,混合极化数据处理及其应用成为近年来研究的热点问题。
混合极化SAR系统存在一个棘手的问题,就是它的简缩极化模式的发射端误差不能进行定标,即发射端误差不能消除,这些误差包括,发射 波Jones矢量误差(φ偏差),通道失衡,通道串扰,以及法拉第旋转角。这些误差都将会影响到获取的混合极化数据质量,进而影响到后续应用。
发明内容
(一)要解决的技术问题
鉴于上述技术问题,本发明提供了一种混合极化SAR系统π/4简缩极化模式发射端误差求解的方法。
(二)技术方案
本发明混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法包括:步骤A:接收发射端误差源的误差变量,其中,所述误差变量为:Jones矢量中误差参数φ、通道失衡参数kt、通道串扰参数∈h,∈v,和/或法拉第旋转角Ω;以及步骤B:按照下式由误差变量求解发射端误差源引起的最大归一化误差:其中,λmax为(PD-P0)H(PD-P0)的最大特征值;P为发射波Jones矢量误差的矩阵形式, D为包含通道失衡误差、通道串扰误差及法拉第旋转角误差的矩阵形式, P0为无误差的矩阵形式,
(三)有益效果
从上述技术方案可以看出,本发明混合极化SAR系统π/4简缩极化模式的发射端误差求解的方法具有以下有益效果:
(1)建立CL简缩极化模式下误差评估模型及详细的误差求解方法。可获得发射端各误差:发射波Jones矢量误差(φ偏差)、通道失衡、通道串扰以及法拉第旋转角对混合极化SAR系统的影响。
(2)在限定系统误差的条件下,可计算出各误差参数的上限值,为混合极化SAR系统的设计提供一定的指导。
附图说明
图1为根据本发明实施例混合极化SAR系统π/4简缩极化模式的发射端误差求解方法的示意流程图;
图2为发射波Jones矢量偏差(椭圆方向角误差)所引起的系统最大归一化误差图;
图3为通道失衡参数所引起的系统最大归一化误差图;
图4为通道串扰所引起的系统最大归一化误差;
图5为法拉第旋转角所引起的系统最大归一化误差。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是在可接受的误差容限或设计约束内近似于相应的值。
本发明提供了一种混合极化SAR系统π/4简缩极化模式的发射端误差求解的方法。该方法将首先推导出最大归一化误差表达式e,而后将误差源代入该最大归一化误差表达式e中,求取该误差源带来的误差。
一、实例一
在本发明的第一个示例性实施例中,提供了混合极化SAR系统π/4简缩极化模式的发射端误差求解的方法。请参照图1,本实施包括:
步骤A:接收发射端误差源的误差变量;
其中,当考虑发射波Jones矢量误差时,误差变量为Jones矢量误差参数φ;当考虑考虑通道失衡参数的影响时,误差变量为通道失衡参数kt,当考虑通道串扰参数的影响时,误差变量为通道串扰参数∈h,∈v;当考虑法拉第旋转角的影响时,误差变量为法拉第旋转角Ω。
步骤B:按照下式求解发射端误差源引起的最大归一化误差:
其中,λmax为(PD-P0)H(PD-P0)的最大特征值。P为发射波Jones矢量 误差的矩阵形式, D为包含通道失衡误差、通道串扰误差及法拉第旋转角误差的矩阵形式, P0为无误差的矩阵形式, 代表矩阵内积。
以下给出式1的推导过程:
步骤S102:根据全极化SAR系统观测模型Sm推导混合极化SAR系统π/4简缩极化模式的误差观测模型M;
该步骤S102进一步包括:
子步骤S102a:依据下式,根据全极化SAR散射矩阵S推导带有误差源的全极化SAR系统π/4简缩极化模式的观测模型Sm;
Sm=XrSXt (1)
其中,Xr表示接收端失真矩阵,Xt表示发射端失真矩阵。在H/V极化基下,Xt与Xr可以分别表示为:
其中,Ω为法拉第旋转角,kt,kr为通道失衡参数,εh,εv,gh,gv为通道串扰参数。
子步骤S102b:依据下式,根据全极化SAR系统观测模型Sm推导混合极化SAR系统π/4模式简缩极化模式的误差观测模型M:
M=PrSmpt=prXrSXtpt (3)
其中,pt=[cosφsinφ]T为混合极化π/4简缩极化模式发射端Jones矢量,φ为表征发射波极化状态的椭圆方向角,上标T表示矩阵转置。pr为两个接收通道Jones矢量组成的2×2矩阵。对于理想π/4简缩极化模式,发射波为 45°线极化波,即φ=45°,因此pr为单位阵。
由于发射端误差可以通过定标消除,因此我们在这里只评估发射端误差,因此忽略pr,Xr,式(3)可简化为:
M=SXtpt (4)
步骤S104:按照辞书式排列,将混合极化观测矩阵M转换成矢量形式M;
M=vec(M)=PDS (5)
其中,S为全极化散射矩阵矢量形式,定义为S=[SHH SHV SVH SVV]T,SHH,SHV,SVH,SVV表示全极化散射矩阵S的HH、HV、VH、VV四通道数据。并且,
矩阵P为2×4矩阵,包含发射波Jones矢量;矩阵D为4×4矩阵,它包含了通道失衡,通道串扰及法拉第旋转角在内的三种误差源。
步骤S106:建立无误差的混合极化SAR观测矢量V;
首先,建立无误差的混合极化SAR观测矩阵V:
然后将式7转化为矢量形式,得到混合极化SAR观测矢量V:
V=P0 S (8)
其中:
步骤S108:根据有误差观测矢量M及无误差矢量V建立误差项e表达式;
e=M-V=(PD-P0)S (10)
步骤S110:根据误差矢量e建立最大归一化误差表达式e:
可见,最大归一化误差转化为求解矩阵的2范数,按照通常的2范数求解方法进行求解,即求解(PD-P0)H(PD-P0)矩阵的最大特征值λmax;根据矩阵论知识可知,得到最大归一化误差:
需要说明的是,误差源可以有一个,也可以为多个。对于存在多个误差源的情况,直接按照式1进行求解即可。以下给出应用上述理论求取仅包含一种误差源场景时的四种典型实施例。
二、实施例二
在本发明的第二个示例性实施例中,提供了求解发射波Jones矢量误差(τ偏差)引起π/4简缩极化模式下的最大归一化误差的方法。
本实施例中,仅考虑发射波Jones矢量误差时,其它误差设置为零,因此矩阵D变成单位阵列,则最大归一化误差表达式e为:
其中,λmax[(P-P0)H(P-P0)]表示(P-P0)H(P-P0)的最大特征值。
将发射波Jones矢量误差的误差参数φ代入式13中,即可以得到其引起的最大归一化误差,具体过程如下:
子步骤S202,求解(P-P0)H(P-P0)的最大特征值λmax。
其中, φ的标准无误差值为45°。任何偏离45°将会带来误差。
子步骤S204,根据下式,求解出最大化归一化误差e:
P中误差变量τ所引起的最大归一化误差如图2所示。图2中,横坐标代表椭圆方向角,单位(°),纵坐标代表最大归一化误差值,单位为dB。由图2可知,若使系统误差不大于-20dB,则椭圆率τ误差就控制在±5°之间。
三、实施例三
在本发明的第三个示例性实施例中,提供了求解通道失衡参数引起的最大归一化误差的方法。
本实施例中,仅考虑通道失衡的影响,则发射波Jones矢量无误差,即P=P0,矩阵D中,Ω,∈h,∈v均为零。误差变量仅有通道失衡参数kt。误差表达式为:
将通道失衡参数代入式14中,即可以得到其引起的最大归一化误差,具体过程如下:
子步骤S302,求解(P0D-P0)H(P0D-P0)的最大特征值λmax。
其中, kt为通道失衡参数。代表矩阵直积。
子步骤S304,根据下式,求解出最大化归一化误差值e:
通道失衡参数引起的最大归一化误差如图3所示。图3中,横坐标代表通道失衡参数的幅度,单位为dB,纵坐标代表通道失衡参数的相位偏移,单位为(°),图中为最大归一化误差的等值线,单位为dB。由图3可知,若使系统误差不大于-20dB,则通道失衡参数幅度应控制在±1.2dB,相位偏差控制在±8°以内。
四、实施例四
在本发明的第四个示例性实施例中,提供了求解通道串扰引起的最大归一化误差的方法。
本实施例中,仅考虑通道串扰的影响,则发射波Jones矢量无误差,即P=P0,矩阵D中,Ω,kt均为零。误差变量仅有通道串扰参数∈h,∈v,且 一般假定∈h=∈v。误差表达式同式15。
将通道串扰参数代入式15中,即可以得到其引起的最大归一化误差,具体过程如下:
子步骤S402,求解(P0D-P0)H(P0D-P0)的最大特征值λmax。
其中, εh,εv代表通道串扰参数。 代表矩阵直积。
子步骤S404,根据下式,求解出最大化归一化误差e:
通道串扰引起的最大归一化误差如图4所示。在图4中,横坐标代表通道串扰,单位为dB,纵坐标代表最大归一化误差,单位为dB。由图4可知,若使系统误差不大于-20dB,则串扰幅度应小于-20dB。
五、实施例五
在本发明的第五个示例性实施例中,提供了求解法拉第旋转角引起的最大归一化误差的方法。
本实施例中,仅考虑法拉第旋转角,则发射波Jones矢量无误差,即P=P0,矩阵D中,∈h,∈v,kt均为零。误差变量仅有法拉第旋转角Ω,误差表达式同式(15)。
将法拉第旋转角误差代入式15中,即可以得到其引起的最大归一化误差,具体过程如下:
子步骤S502,求解(P0D-P0)H(P0D-P0)的最大特征值λmax。
其中, Ω代表法拉第旋转角。代表矩阵直积。
子步骤S504,根据下式,求解出最大化归一化误差e。
法拉第旋转角引起的最大归一化误差如图5所示。图5中,横坐标代表法拉第旋转角,单位为(°),纵坐标代表最大归一化误差,单位为dB。由图5可知,若使系统误差不大于-20dB,则法拉第旋转角不得超过4°。
至此,已经结合附图对本发明实施例进行了详细描述。依据以上描述,本领域技术人员应对本发明混合极化SAR系统π/4简缩极化模式的发射端误差求解的方法有了清楚的认识。
综上所述,本发明提供了混合极化SAR系统π/4简缩极化模式发射端误差求解的方法。对混合极化SAR发射端误差源:包括发射波Jones矢量偏差,通道失衡,通道串扰,法拉第旋转角误差所引起的系统误差进行清晰明了的计算与分析,并且给出了在系统误差一定的前提下各误差源的上限值。该发明为混合极化SAR系统误差分析及实际系统设计提供了新的技术手段及支持。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法,其特征在于,包括:
步骤A:接收发射端误差源的误差变量,其中,所述误差变量为:Jones矢量中误差参数φ、通道失衡参数kt、通道串扰参数∈h,∈v,和/或法拉第旋转角Ω;以及
步骤B:按照下式由误差变量求解发射端误差源引起的最大归一化误差:
其中,λmax为(PD-P0)H(PD-P0)的最大特征值;P为发射波Jones矢量误差的矩阵形式, D为包含通道失衡误差、通道串扰误差及法拉第旋转角误差的矩阵形式, P0为无误差的矩阵形式,
2.一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法,其特征在于,误差源仅为发射波Jones矢量误差,该方法包括:
步骤A:接收发射端Jones矢量误差参数φ;
步骤B:按照下式求解由发射波Jones矢量误差引起的最大归一化误差:
其中,λmax为(P-P0)H(P-P0)的最大特征值;P为发射波Jones矢量误差的矩阵形式, P0为无误差的矩阵形式,
3.一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法,其特征在于,误差源仅为通道失衡,该方法包括:
步骤A:接收通道失衡参数kt;
步骤B:按照下式求解通道失衡引起的最大归一化误差:
其中,λmax为(P0D-P0)H(P0D-P0)的最大特征值,P0为无误差的矩阵形式, D为包含通道失衡误差的矩阵形式,
4.一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法,其特征在于,误差源仅为通道串扰,该方法包括:
步骤A:接收通道串扰参数∈h,∈v;
步骤B:按照下式求解通道串扰引起的最大归一化误差:
其中,λmax为(P0D-P0)H(P0D-P0)的最大特征值,P0为无误差的矩阵形式, D为包含通道串扰误差的矩阵形式,
5.一种混合极化SAR系统π/4简缩极化模式下发射端误差求解的方法,其特征在于,误差源仅为法拉第旋转角,该方法包括:
步骤A:接收法拉第旋转角Ω;
步骤B:按照下式求解法拉第旋转角引起的最大归一化误差:
其中,λmax为(P0D-P0)H(P0D-P0)的最大特征值,P0为无误差的矩阵形式, D为包含法拉第旋转角引起的误差的矩阵形式,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410743634.6A CN104459655B (zh) | 2014-12-08 | 2014-12-08 | 混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410743634.6A CN104459655B (zh) | 2014-12-08 | 2014-12-08 | 混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459655A true CN104459655A (zh) | 2015-03-25 |
CN104459655B CN104459655B (zh) | 2017-04-26 |
Family
ID=52906008
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410743634.6A Active CN104459655B (zh) | 2014-12-08 | 2014-12-08 | 混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459655B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137433A (zh) * | 2015-08-27 | 2015-12-09 | 北京航空航天大学 | 一种星载π/4模式简缩极化SAR相位重构方法 |
CN108051790A (zh) * | 2017-12-13 | 2018-05-18 | 厦门大学 | 一种用于混合紧缩极化探地雷达系统的极化校准方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110175771A1 (en) * | 2007-05-08 | 2011-07-21 | Raney Russell K | Synthetic Aperture Radar Hybrid-Quadrature-Polarity Method and Architecture for Obtaining the Stokes Parameters of Radar Backscatter |
CN102253374A (zh) * | 2011-04-20 | 2011-11-23 | 北京航空航天大学 | 一种长波长星载π/4模式简缩极化SAR的定标方法 |
CN103941245A (zh) * | 2014-05-09 | 2014-07-23 | 中国科学院电子学研究所 | 简缩极化干涉数据的Freeman-Durden目标分解方法 |
-
2014
- 2014-12-08 CN CN201410743634.6A patent/CN104459655B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110175771A1 (en) * | 2007-05-08 | 2011-07-21 | Raney Russell K | Synthetic Aperture Radar Hybrid-Quadrature-Polarity Method and Architecture for Obtaining the Stokes Parameters of Radar Backscatter |
CN102253374A (zh) * | 2011-04-20 | 2011-11-23 | 北京航空航天大学 | 一种长波长星载π/4模式简缩极化SAR的定标方法 |
CN103941245A (zh) * | 2014-05-09 | 2014-07-23 | 中国科学院电子学研究所 | 简缩极化干涉数据的Freeman-Durden目标分解方法 |
Non-Patent Citations (3)
Title |
---|
YANTING WANG,等: "Assessment of System Polarization Quality for Polarimetric SAR Imagery and Target Decomposition", 《IEEE TRANSACTIONS ON GEOSCIENCR AND REMOTE SENSTNG》 * |
陈琳,等: "简缩极化干涉合成孔径雷达三分量目标分解方法", 《中国科学院大学学报》 * |
陈琳: "单发双收SAR系统通用极化定标算法", 《雷达学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137433A (zh) * | 2015-08-27 | 2015-12-09 | 北京航空航天大学 | 一种星载π/4模式简缩极化SAR相位重构方法 |
CN108051790A (zh) * | 2017-12-13 | 2018-05-18 | 厦门大学 | 一种用于混合紧缩极化探地雷达系统的极化校准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104459655B (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405547B (zh) | 一种多极化sar简缩极化模式定标方法 | |
CN103067096B (zh) | 准确测试天线两个正交线极化分量的校准方法 | |
CN102080988B (zh) | 一种实时检测单光子偏振量子态的装置与方法 | |
CN111948615B (zh) | 一种星载全极化sar数据的极化定标方法及装置 | |
CN104111448A (zh) | Mimo雷达收发阵列误差的联合校正方法 | |
CN106970361A (zh) | 干扰条件下双极化单脉冲雷达目标角度的估计方法 | |
US20110248881A1 (en) | Radar with High Angular Accuracy, Notably for the Obstacle Sensing and Avoidance Function | |
CN105353229A (zh) | 一种基于一维旋转的相控阵幅相误差近场校准方法 | |
CN102183743B (zh) | 一种长波长星载ctlr模式简缩极化sar的定标方法 | |
CN101251597A (zh) | 一种多输入多输出雷达系统阵列误差自校正的方法 | |
CN104267383A (zh) | 一种针对雷达电磁信号的极化参数自适应测量装置 | |
CN103323832B (zh) | 一种相控阵三维摄像声纳系统换能器阵列的幅相误差校正方法 | |
CN106227701A (zh) | 一种阵列信号的接收通道的幅相误差的自校正方法 | |
CN108110419A (zh) | 一种修正天线罩瞄准误差的方法 | |
CN104459655A (zh) | 混合极化SAR系统π/4简缩极化模式发射端误差求解的方法 | |
CN103941245A (zh) | 简缩极化干涉数据的Freeman-Durden目标分解方法 | |
CN106980104A (zh) | 用于传感器阵列的信号波达方向自校正方法 | |
CN103926555A (zh) | 一种利用非圆信号测定天线阵列接收机幅相响应的方法 | |
CN103576151B (zh) | 基于压缩感知的方位向多通道sar成像方法及系统 | |
CN105046072B (zh) | 基于压缩感知理论的二维到达角估计方法 | |
CN103760518B (zh) | 多个信号方向未知的天线阵列方向向量的测定方法 | |
Touzi et al. | Assessment of PALSAR-2 compact non-circularity using Amazonian rainforests | |
CN102520407A (zh) | 基于相干系数线性度最优的极化干涉sar图像配准方法 | |
Weishaupt et al. | Calibration and signal processing of polarimetric radar data in automotive applications | |
CN106771690A (zh) | 一种固定式准光法拉第旋转器性能测定方法与装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |