发明名称 一种受Faraday旋转变化影响的雷达目标后向散射模拟方法
摘要 本发明提供了一种受Faraday旋转变化影响的雷达目标后向散射模拟方法,包括:获取同一场景的四幅校准各种误差影响后的全极化SAR复图像;获取目标的极化散射矩阵S;计算每个脉冲发射时刻由电离层TEC值引起的Faraday旋转角Ω;计算受Faraday旋转变化影响的极化散射矩阵O;获取受Faraday旋转影响的后向散射系数;将受Faraday旋转变化影响的后向散射系数带入到SAR回波表达式,做成像处理及相应的分析。采用本方法,可模拟任意Faraday旋转变化对后向散射系数及SAR成像的影响。而且,基于已有的SAR图像模拟受Faraday旋转变化影响的雷达目标后向散射系数,简单有效。
申请公布号 CN103792521B 申请公布日期 2016.05.04
申请号 CN201410022836.1 申请日期 2014.01.17
申请人 西安空间无线电技术研究所 发明人 王旭艳;王伟伟;赵泓懿;朱雅琳;李财品;李光廷;黎薇萍;刘波
分类号 G01S7/41(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 一种受Faraday旋转变化影响的SAR目标后向散射模拟方法,其特征在于,包括以下步骤:第一步、获取同一场景的四幅全极化SAR复图像,且该四幅图像是校准各种误差影响后的图像;第二步、获取目标的极化散射矩阵S:以分辨单元为单位,在四幅复图像中分别取得该分辨单元对应的复像素值,作为极化散射矩阵S中相应的元素,该极化散射矩阵S表示如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>S</mi><mrow><mi>h</mi><mi>h</mi></mrow></msub></mtd><mtd><msub><mi>S</mi><mrow><mi>h</mi><mi>v</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mrow><mi>v</mi><mi>h</mi></mrow></msub></mtd><mtd><msub><mi>S</mi><mrow><mi>v</mi><mi>v</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000814152870000011.GIF" wi="372" he="158" /></maths>上式中,S为一个2×2阶的矩阵,S<sub>hh</sub>、S<sub>hv</sub>、S<sub>vh</sub>、S<sub>vv</sub>为分别在四幅复图像中取得的复像素值,h表示水平极化,v表示垂直极化;第三步、根据下面的公式计算每个脉冲发射时刻由电离层TEC引起的Faraday旋转角Ω:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Omega;</mi><mo>=</mo><mn>2.36</mn><mo>&times;</mo><msup><mn>10</mn><mn>4</mn></msup><mfrac><msub><mi>B</mi><mrow><mi>a</mi><mi>v</mi></mrow></msub><mrow><msup><msub><mi>f</mi><mi>c</mi></msub><mn>2</mn></msup></mrow></mfrac><mi>T</mi><mi>E</mi><mi>C</mi><mo>,</mo></mrow>]]></math><img file="FDA0000814152870000012.GIF" wi="496" he="135" /></maths>式中,B<sub>av</sub>为地磁场沿传播路径上的平均值,f<sub>c</sub>为发射信号的载频,TEC为SAR信号传播路径上的电离层总电子含量;第四步、根据下面的公式计算受Faraday旋转变化影响的极化散射矩阵O:O=F·S·F,式中,S为第二步中获取的极化散射矩阵,F表示Faraday旋转对目标散射特性的影响,其表达式为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&Omega;</mi></mrow></mtd><mtd><mrow><mi>s</mi><mi>i</mi><mi>n</mi><mi>&Omega;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&Omega;</mi></mrow></mtd><mtd><mrow><mi>cos</mi><mi>&Omega;</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000814152870000013.GIF" wi="461" he="159" /></maths>O的表达式为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>O</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>O</mi><mrow><mi>h</mi><mi>h</mi></mrow></msub></mtd><mtd><msub><mi>O</mi><mrow><mi>h</mi><mi>v</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>O</mi><mrow><mi>v</mi><mi>h</mi></mrow></msub></mtd><mtd><msub><mi>O</mi><mrow><mi>v</mi><mi>v</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000814152870000014.GIF" wi="333" he="157" /></maths>其中,h表示水平极化,v表示垂直极化;第五步、根据下面的公式计算受Faraday旋转影响的观测到的后向散射系数,即,两个同极化通道观测到的后向散射系数P<sub>hh</sub>、P<sub>vv</sub>,两个交叉极化通道观测到的后向散射系数P<sub>hv</sub>、P<sub>vh</sub>:P<sub>hh</sub>=O<sub>hh</sub>,P<sub>vv</sub>=O<sub>vv</sub>,P<sub>vh</sub>=O<sub>vh</sub>,P<sub>hv</sub>=O<sub>hv</sub>;第六步、将第五步中获取的观测到的后向散射系数带入下面的SAR回波信号表达式中,对该回波信号表达式表述的回波信号进行成像处理,并对成像结果做相应的分析,下式中,X分别代表P<sub>hh</sub>、P<sub>vh</sub>、P<sub>hv</sub>、P<sub>vv</sub>,以分别分析Faraday旋转变化对水平极化发射水平极化接收信号成像、对水平极化发射垂直极化接收信号成像、对垂直极化发射水平极化接收信号成像、对垂直极化发射垂直极化接收信号成像的影响:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mrow><mo>(</mo><mrow><mi>t</mi><mo>,</mo><msub><mi>t</mi><mi>m</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mi>X</mi><mo>&CenterDot;</mo><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mrow><mo>&lsqb;</mo><mrow><mi>t</mi><mo>-</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><msub><mi>R</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow></mrow><mi>c</mi></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>&CenterDot;</mo><msub><mi>w</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>w</mi><mi>r</mi></msub><mrow><mo>&lsqb;</mo><mrow><mi>t</mi><mo>-</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><msub><mi>R</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow></mrow><mi>c</mi></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>&CenterDot;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>exp</mi><mrow><mo>{</mo><mrow><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mrow><mo>&lsqb;</mo><mrow><mi>t</mi><mo>-</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><msub><mi>R</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow></mrow><mi>c</mi></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>+</mo><mi>j</mi><mi>&pi;</mi><mi>&gamma;</mi><msup><mrow><mo>&lsqb;</mo><mrow><mi>t</mi><mo>-</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><msub><mi>R</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow></mrow><mi>c</mi></mfrac></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow><mo>}</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000814152870000021.GIF" wi="1283" he="365" /></maths>上式中,j为虚数单位,j=sqrt(‑1);t、t<sub>m</sub>分别为距离向时间、方位向时间;c为光速;w<sub>a</sub>(t<sub>m</sub>)、<img file="FDA0000814152870000022.GIF" wi="404" he="159" />分别为方位向天线方向图、距离向天线方向图;f<sub>c</sub>为载频;γ为调频率;R<sub>B</sub>为目标与SAR之间的最近距离;R(t<sub>m</sub>,R<sub>B</sub>)为时刻t<sub>m</sub>时SAR与目标之间的瞬时斜距。
地址 710100 陕西省西安市长安区西街150号
您可能感兴趣的专利