发明名称 一种星载SAR数据的电离层色散效应校正方法
摘要 本发明提出一种电离层色散效应校正方法,应用于星载SAR数据的电离层色散效应校正。技术方案是以星载SAR图像及其载频、带宽、采样率为输入,利用下面5个步骤进行迭代,第一步:对输入星载SAR图像分别进行两次距离向带通滤波得到两个子带图像;第二步,通过两个子带图像间的距离向实相关计算距离向实相关函数;第三步,利用距离向实相关函数估计两个子带图像的相对偏移量;第四步,利用两个子带图像的偏移量估计值估计数据中的TEC;第五步,对输入星载SAR图像中的电离层色散效应进行校正。本发明不依赖于场景中的强散射体,也不依赖于其他测量设备,可以使受电离层色散效应影响而散焦的图像重新聚焦,正确定位。
申请公布号 CN103760534A 申请公布日期 2014.04.30
申请号 CN201410029577.5 申请日期 2014.01.19
申请人 中国人民解放军国防科学技术大学 发明人 董臻;张永胜;李力;余安喜;何峰;赵腾飞;黄海风;孙造宇;金光虎;何志华;陈祺;马喜乐
分类号 G01S7/40(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 1.一种星载SAR数据的电离层色散效应校正方法,已知星载SAR系统载频为f<sub>c</sub>,带宽为B,采样频率为F<sub>s</sub>;星载SAR图像I为M×N维复数矩阵,其中M是方位向点数,N是距离向点数,其特征在于包括下述步骤:设k=1,星载SAR图像I<sub>k</sub>=I,然后进行下述步骤:第一步、带通滤波,获取子带图像对星载SAR图像I<sub>k</sub>进行两次带通滤波,分别获得两个带宽为B2的子带图像<img file="FDA0000458701190000011.GIF" wi="50" he="79" />和<img file="FDA0000458701190000012.GIF" wi="80" he="79" />其中子带图像<img file="FDA0000458701190000013.GIF" wi="52" he="79" />的载频为f<sub>c</sub>-B4,子带图像<img file="FDA0000458701190000014.GIF" wi="51" he="79" />的载频为f<sub>c</sub>+B4;第二步、求距离向实相关函数计算子带图像<img file="FDA0000458701190000015.GIF" wi="58" he="79" />和子带图像<img file="FDA0000458701190000016.GIF" wi="57" he="79" />的距离向相关函数ρ<sub>k</sub>(l),其中l=-(N-1)…-2,-1,0,1,2…N-1;第三步、估计实相关函数最大值位置利用下式计算子带图像<img file="FDA0000458701190000017.GIF" wi="50" he="78" />和子带图像<img file="FDA0000458701190000018.GIF" wi="52" he="79" />间的距离向偏移量估计值<img file="FDA0000458701190000019.GIF" wi="102" he="73" /><maths num="0001"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mover><mi>r</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mi>x</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>&rho;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中x∈[-(N-2),(N-2)],并且取值为整数,同时满足下式:ρ<sub>k</sub>(x)=max(ρ<sub>k</sub>(l)<sub>-(N-2)…-2,-1,0,1,2…N-2</sub>)第四步、估计TEC利用下式计算TEC的估计值<img file="FDA00004587011900000111.GIF" wi="152" he="79" /><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>TEC</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mi>&Delta;</mi><msub><mover><mi>r</mi><mo>^</mo></mover><mi>k</mi></msub><mo>&CenterDot;</mo><mfrac><mi>&pi;c</mi><mrow><msub><mi>r</mi><mi>e</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>&lambda;</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&lambda;</mi><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中π是圆周率;c是真空中的光速;r<sub>e</sub>是经典电子半径;λ<sub>1</sub>=c(f<sub>c</sub>-B4),λ<sub>2</sub>=c(f<sub>c</sub>+B4);第五步、校正电离层色散效应利用下述计算第k次校正电离层色散效应后的星载SAR图像<img file="FDA00004587011900000113.GIF" wi="61" he="55" /><img file="FDA0000458701190000021.GIF" wi="426" he="79" />上式中,E是M×N维矩阵,E(m,n)是矩阵E第m行、第n列的元素,且<img file="FDA0000458701190000022.GIF" wi="727" he="97" />m∈[1,M],n∈[1,N];⊙是Hadamard积运算;FT<sub>r</sub>{·}表示距离向的傅立叶变换,IFT<sub>r</sub>{·}表示距离向的逆傅立叶变换;<img file="FDA0000458701190000023.GIF" wi="183" he="83" />表示虚数单位;然后进行下述判断:如果<img file="FDA0000458701190000024.GIF" wi="300" he="79" />或k=Q,则输出星载SAR图像<img file="FDA0000458701190000025.GIF" wi="76" he="71" />本技术方案结束;否则令<img file="FDA0000458701190000026.GIF" wi="200" he="72" />k=k+1,返回第一步;迭代精度门限ΔTEC和迭代终止次数Q的取值根据实际需要确定。
地址 410073 湖南省长沙市开福区德雅路109号