发明名称 GEOSAR方位成像中的电离层时变效应影响判决方法
摘要 本发明提供一种GEOSAR方位成像中的电离层时变效应影响判决方法。技术方案是:以GEOSAR系统参数为输入,第一步:针对观测场景中任意目标点P,计算得到其合成孔径中心时刻t<sub>0</sub>的电离层穿刺点坐标<img file="DDA0000984905920000011.GIF" wi="323" he="101" />第二步:利用在线IRI模型获取方位向慢时间t<sub>n</sub>时刻穿刺点坐标处的电离层<img file="DDA0000984905920000012.GIF" wi="518" he="96" />第三步:对得到的<img file="DDA0000984905920000013.GIF" wi="497" he="95" />进行多项式拟合,得到一阶分量系数k<sub>1</sub>的值和二阶分量系数k<sub>2</sub>的值;第四步:判断k<sub>1</sub>和k<sub>2</sub>与门限值的大小,得到判决结果。本发明针对星载SAR方位成像处理,提出了电离层时变效应的影响判决方法,给出了判决门限和判决方法,在星载SAR电离层影响分析与校正处理中有广泛应用。
申请公布号 CN105785336A 申请公布日期 2016.07.20
申请号 CN201610301841.5 申请日期 2016.05.09
申请人 中国人民解放军国防科学技术大学 发明人 张启雷;计一飞;张永胜;董臻;余安喜;何志华;黄海风;何峰;孙造宇;金光虎
分类号 G01S7/40(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 一种GEOSAR方位成像中的电离层时变效应的影响判决方法,GEOSAR是指地球同步轨道合成孔径雷达,已知GEOSAR的系统参数包括:载频f<sub>c</sub>,合成孔径时间T<sub>s</sub>;观测场景中任意目标点P在东北天坐标系下的坐标为(x<sub>0</sub>,y<sub>0</sub>,z<sub>0</sub>),该点对应的合成孔径中心时刻为t<sub>0</sub>;GEOSAR在方位向慢时间t<sub>n</sub>时刻的东北天坐标系轨道坐标数据为<img file="FDA0000984905890000011.GIF" wi="299" he="95" />n的取值为任意整数并且|t<sub>n</sub>‑t<sub>0</sub>|≤T<sub>s</sub>/2;电离层的高度为z<sub>i</sub>;其特征在于,包括下述步骤:第一步,针对观测场景中任意目标点P,利用空间直线方程,计算得到合成孔径中心时刻t<sub>0</sub>的GEOSAR波束在电离层上的穿刺点坐标<img file="FDA0000984905890000012.GIF" wi="323" he="95" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow><mrow><msub><mi>z</mi><msub><mi>t</mi><mn>0</mn></msub></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>x</mi><msub><mi>t</mi><mn>0</mn></msub></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>x</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000984905890000013.GIF" wi="645" he="142" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow><mrow><msub><mi>z</mi><msub><mi>t</mi><mn>0</mn></msub></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>y</mi><msub><mi>t</mi><mn>0</mn></msub></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>y</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000984905890000014.GIF" wi="662" he="143" /></maths>第二步,针对穿刺点坐标<img file="FDA0000984905890000015.GIF" wi="323" he="95" />利用在线国际参考电离层模型获取方位向慢时间t<sub>n</sub>时刻的电离层垂直向电子总量的值<img file="FDA0000984905890000016.GIF" wi="482" he="79" />再利用下式得到时变的<img file="FDA0000984905890000017.GIF" wi="499" he="95" />值:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>S</mi><mi>T</mi><mi>E</mi><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>T</mi><mi>E</mi><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&gamma;</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000984905890000018.GIF" wi="1165" he="102" /></maths>其中,几何变换因子γ(t<sub>0</sub>)为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&gamma;</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mrow><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000984905890000019.GIF" wi="1052" he="189" /></maths>第三步,对得到的<img file="FDA00009849058900000110.GIF" wi="499" he="102" />进行多项式拟合,得到一阶分量系数k<sub>1</sub>的值和二阶分量系数k<sub>2</sub>的值;第四步,如果下述不等式同时成立,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>|</mo><msub><mi>k</mi><mn>1</mn></msub><mo>|</mo><mo>&le;</mo><mfrac><mrow><mn>0.886</mn><mo>&times;</mo><msub><mi>cf</mi><mi>c</mi></msub></mrow><mrow><mn>4</mn><mi>K</mi></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mn>1</mn><msub><mi>T</mi><mi>s</mi></msub></mfrac></mrow>]]></math><img file="FDA0000984905890000021.GIF" wi="462" he="149" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>|</mo><msub><mi>k</mi><mn>2</mn></msub><mo>|</mo><mo>&le;</mo><mfrac><mrow><msub><mi>cf</mi><mi>c</mi></msub></mrow><mrow><mn>4</mn><mi>K</mi></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mn>1</mn><msubsup><mi>T</mi><mi>s</mi><mn>2</mn></msubsup></mfrac></mrow>]]></math><img file="FDA0000984905890000022.GIF" wi="318" he="150" /></maths>则对观测场景中任意目标点P的方位向成像时忽略电离层时变效应的影响,其中,c=3×10<sup>8</sup>m/s为真空中光速,K=40.28m<sup>3</sup>/s<sup>2</sup>;否则,必须考虑其影响。
地址 410073 湖南省长沙市开福区德雅路109号