[发明专利]一种长波长星载π/4模式简缩极化SAR的定标方法无效
申请号: | 201110099474.2 | 申请日: | 2011-04-20 |
公开(公告)号: | CN102253374A | 公开(公告)日: | 2011-11-23 |
发明(设计)人: | 陈杰;尹巽军;李卓;李春升 | 申请(专利权)人: | 北京航空航天大学 |
主分类号: | G01S7/40 | 分类号: | G01S7/40;G01S13/90 |
代理公司: | 北京永创新实专利事务所 11121 | 代理人: | 赵文利 |
地址: | 100191*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种长波长星载π/4模式简缩极化SAR的定标方法,包括以下几个步骤:步骤一:获取发射通道不平衡度的估计值;步骤二:获取垂直发射水平接收分量的通道串扰系数的估计值;步骤三:获取法拉第旋转角的估计值;步骤四:利用全球导航卫星系统提供的TEC数据获取无模糊的法拉第旋转角估计值;步骤五:获取接收通道不平衡度的估计值;步骤六:获取水平发射垂直接收分量的通道串扰系数的估计值。本发明在没有定标器误差及系统噪声时能进行精确的定标,在存在定标器误差及系统噪声时仍能进行高精度的定标。本发明提出的方法对以P波段为代表的长波长星载π/4模式简缩极化SAR数据处理具有十分重要的应用价值。 | ||
搜索关键词: | 一种 波长 模式 简缩 极化 sar 定标 方法 | ||
【主权项】:
1.一种长波长星载π/4模式简缩极化SAR的定标方法,基于π/4模式简缩极化SAR的误差模型,模型如式(1)所示:M EH M EV = 1 2 1 δ 2 δ 1 f 1 · cos Ω sin Ω - sin Ω cos Ω · S HH S HV S VH S VV · cos Ω + f 2 sin Ω - sin Ω + f 2 cos Ω - - - ( 1 ) ]]> 式中,MEV和MEH分别表示测量散射矢量的接收垂直极化分量和接收水平极化分量,SHH、SHV、SVH、SVV构成的矩阵表示目标的散射矩阵,分别表示矩阵中发射水平极化接收水平极化分量、发射垂直极化接收水平极化分量、发射水平极化接收垂直极化分量、发射垂直极化接收垂直极化分量;Ω表示法拉第旋转角,f1表示接收通道不平衡度,f2表示发射通道不平衡度,δ1表示水平发射垂直接收分量的通道串扰系数,δ2表示垂直发射水平接收分量的通道串扰系数;四个主动极化定标器表示为X、Y、A和B,其散射矩阵分别为:S X = S HH S HV S VH S VV = 0 0 1 0 ]]>S Y = S HH S HV S VH S VV = 0 1 0 0 ]]>S A = S HH S HV S VH S VV = 1 0 0 0 ]]>S B = S HH S HV S VH S VV = 0 0 0 1 - - - ( 2 ) ]]> 其特征在于,具体包括以下几个步骤;步骤一:获取发射通道不平衡度的估计值根据式(3)及式(4)计算接收通道不平衡度的估计值式中,|x|为取复数x的模运算,且:x 1 = 1 2 [ ( M EH X + M EH Y ) + ( M EH X + M EH Y ) 2 - 4 ( 1 - M EH A + M EH B ) ] x 2 = 1 2 [ ( M EH X + M EH Y ) - ( M EH X + M EH Y ) 2 - 4 ( 1 - M EH A + M EH B ) ] - - - ( 4 ) ]]> 式中:表示测量得到的定标器X散射矢量的水平接收分量,表示测量得到的定标器Y散射矢量的水平接收分量,表示测量得到的定标器A散射矢量的水平接收分量,表示测量得到的定标器B散射矢量的水平接收分量;步骤二:获取垂直发射水平接收分量的通道串扰系数的估计值根据式(5)计算垂直发射水平接收分量的通道串扰系数的估计值δ ^ 2 = M EH X + M EH Y - f ^ 2 - - - ( 5 ) ]]> 步骤三:获取法拉第旋转角的估计值根据式(6)及式(7)计算法拉第旋转角的估计值Ω ^ F = 1 2 arg ( cos 2 Ω ^ F + j sin 2 Ω ^ F ) - - - ( 6 ) ]]> 式中,arg(x)为取复数x的幅角运算,且cos 2 Ω ^ F = real [ ( M EH A + M EH B ) ( 2 - M EH A + M EH B ) - ( M EH X - M EH Y ) ( 2 f ^ 2 - M EH X - M EH Y ) ( M EH A + M EH B ) 2 + ( M EH X - M EH Y ) 2 ] sin 2 Ω ^ F = real [ ( M EH X - M EH Y ) ( 2 - M EH A + M EH B ) + ( M EH A + M EH B ) ( 2 f ^ 2 - M EH X - M EH Y ) ( M EH A + M EH B ) 2 + ( M EH X - M EH Y ) 2 ] - - - ( 7 ) ]]> 式中,real(x)为取复数x的实部运算;步骤四:利用全球卫星导航系统提供的TEC数据获取无模糊的法拉第旋转角估计值利用全球导航卫星系统提供的电离层TEC观测数据,并结合第10代国际参考地磁场模型地磁计算模型,根据式(8)求得法拉第旋转角的粗略估计值Ω ^ GNSS ≈ K f 0 2 · [ B cos ψ · sec θ ] 400 · TEC - - - ( 8 ) ]]> 式中,f0表示SAR系统的工作频率,单位为Hz,K为常量2.365×104A·m2/kg,B表示地球磁场强度,单位为Wb/m2,θ表示星载SAR天线的视角,ψ表示地球磁场方向与雷达电磁波传播方向的夹角,TEC为在垂直于地面方向上的电离层电子总含量,单位为TECU,1TECU=1016m-2;地球磁场因子[Bcosψ·secθ]400在400公里高度上计算得到;采用式(9)对法拉第旋转角估计值的进行解模糊,获取无模糊的法拉第旋转角估计值Ω ^ = Ω ^ F + round ( Ω ^ GNSS - Ω ^ F π / 2 ) · π 2 - - - ( 9 ) ]]> 式中,round(x)为取与x最接近的整数值;步骤五:获取接收通道不平衡度的估计值根据式(10)计算接收通道不平衡度的估计值f ^ 1 = 1 2 [ ( M EV X - M EV Y ) cos 2 Ω ^ - ( M EV A + M EV B ) sin 2 Ω ^ + ( M EV A + M EV B ) ] - - - ( 10 ) ]]> 式中:表示测量得到的定标器X散射矢量的垂直接收分量,表示测量得到的定标器Y散射矢量的垂直接收分量,表示测量得到的定标器A散射矢量的垂直接收分量,表示测量得到的定标器B散射矢量的垂直接收分量;步骤六:获取水平发射垂直接收分量的通道串扰系数的估计值根据式(11)计算水平发射垂直接收分量的通道串扰系数的估计值δ ^ 1 = 1 2 [ ( M EV X - M EV Y ) sin 2 Ω ^ + ( M EV A + M EV B ) cos 2 Ω ^ + ( M EV A - M EV B ) ] - - - ( 11 ) ]]> 经过以上六个步骤,完成了对法拉第旋转角的估计、发射及接收通道不平衡度误差的估计,及通道串扰系数的估计,完成π/4模式简缩极化SAR定标处理,实现定标。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201110099474.2/,转载请声明来源钻瓜专利网。