主权项 |
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>Δ</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>ρ</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>ρ</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>ρ</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>ρ</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>ρ</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>Δ</mi><msub><mover><mi>r</mi><mo>^</mo></mover><mi>k</mi></msub><mo>·</mo><mfrac><mi>πc</mi><mrow><msub><mi>r</mi><mi>e</mi></msub><mo>·</mo><mrow><mo>(</mo><msubsup><mi>λ</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>λ</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的取值根据实际需要确定。 |