发明名称 星载SAR方位成像中的电离层空变效应影响判决方法
摘要 本发明提供一种星载SAR方位成像中的电离层空变效应影响判决方法。技术方案是:以星载SAR系统参数为输入,第一步:针对观测场景中任意目标点P,计算得到方位向慢时间t<sub>n</sub>时刻的电离层穿刺点坐标<img file="DDA0000984906900000011.GIF" wi="323" he="101" />第二步:针对所有的电离层穿刺点坐标<img file="DDA0000984906900000012.GIF" wi="323" he="95" />利用在线IRI模型获取t<sub>0</sub>时刻的STEC<img file="DDA0000984906900000013.GIF" wi="371" he="96" />第三步:计算电离层STEC沿方位向的平均空间变化率ΔSTEC;第四步:判断ΔSTEC与门限值的大小,得到判决结果。本发明针对星载SAR方位成像处理,提出了电离层空变效应的影响判决方法,给出了判决门限和判决方法,在星载SAR电离层影响分析与校正处理中有广泛应用。
申请公布号 CN106019279A 申请公布日期 2016.10.12
申请号 CN201610300802.3 申请日期 2016.05.09
申请人 中国人民解放军国防科学技术大学 发明人 董臻;张启雷;计一飞;张永胜;余安喜;何志华;黄海风;何峰;孙造宇;金光虎
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 一种星载SAR方位成像中的电离层空变效应影响判决方法,SAR是指合成孔径雷达,已知星载SAR的系统参数包括:载频f<sub>c</sub>,合成孔径时间T<sub>s</sub>,方位向分辨率ρ<sub>a</sub>;观测场景中任意目标点P在东北天坐标系下的坐标为(x<sub>0</sub>,y<sub>0</sub>,z<sub>0</sub>),该点对应的合成孔径中心时刻为t<sub>0</sub>;星载SAR在方位向慢时间t<sub>n</sub>时刻的东北天坐标系轨道坐标数据为<img file="FDA0000984906870000011.GIF" wi="299" he="102" />n的取值为任意整数并且|t<sub>n</sub>‑t<sub>0</sub>|≤T<sub>s</sub>/2;电离层的高度为z<sub>i</sub>;其特征在于,包括下述步骤:第一步,针对观测场景中任意目标点P,利用空间直线方程,计算得到方位向慢时间t<sub>n</sub>时刻的星载SAR波束在电离层上的东北天坐标系穿刺点坐标<img file="FDA0000984906870000012.GIF" wi="318" he="95" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></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><mi>n</mi></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><mi>n</mi></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="FDA0000984906870000013.GIF" wi="651" he="151" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></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><mi>n</mi></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><mi>n</mi></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="FDA0000984906870000014.GIF" wi="667" he="145" /></maths>第二步,针对所有的穿刺点坐标<img file="FDA0000984906870000015.GIF" wi="323" he="95" />利用在线国际参考电离层模型获取获取t0时刻的电离层垂直向TEC电子总量值<img file="FDA0000984906870000016.GIF" wi="483" he="79" />再利用下式得到空变的<img file="FDA0000984906870000017.GIF" wi="498" he="103" />的值:<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><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mn>0</mn></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><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mn>0</mn></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="FDA0000984906870000018.GIF" wi="1165" he="94" /></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="FDA0000984906870000019.GIF" wi="1054" he="182" /></maths>第三步,计算电离层STEC沿方位向的平均空间变化率ΔSTEC:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>S</mi><mi>T</mi><mi>E</mi><mi>C</mi><mo>=</mo><mi>E</mi><mo>{</mo><mfrac><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><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><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><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub></mrow></msub><mo>,</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub></mrow></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><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></msub><mo>-</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mfrac><mo>}</mo></mrow>]]></math><img file="FDA00009849068700000110.GIF" wi="1420" he="222" /></maths>其中,E{·}代表求平均;第四步,如果下述不等式成立,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>|</mo><mi>&Delta;</mi><mi>S</mi><mi>T</mi><mi>E</mi><mi>C</mi><mo>|</mo><mo>&le;</mo><mfrac><msubsup><mi>f</mi><mi>c</mi><mn>2</mn></msubsup><mrow><mn>2</mn><msub><mi>Kz</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&rho;</mi><mi>a</mi></msub></mrow>]]></math><img file="FDA0000984906870000021.GIF" wi="462" he="159" /></maths>则对观测场景中任意目标点P的方位成像时忽略电离层空变效应的影响,其中K=40.28m<sup>3</sup>/s<sup>2</sup>;否则,必须考虑其影响。
地址 410073 湖南省长沙市开福区德雅路109号