主权项 |
一种基于非完全数据解卷积的实孔径前视成像方法,包括如下步骤:(1)通过雷达波束的等间距发射,对地面监视区域进行顺序重叠扫描,获取雷达回波数据g<sub>T</sub>(x,y):g<sub>T</sub>(x,y)=f(x,y)*h(x,y)+n(x,y),其中,f(x,y)表示地面场景,h(x,y)表示天线方向图,n(x,y)表示噪声,*表示卷积,x表示矩阵元素的横坐标,y表示矩阵元素的纵坐标;(2)利用线性插值扩展雷达回波数据g<sub>T</sub>(x,y),得到近似完全卷积数据<img file="FDA0000452328160000011.GIF" wi="302" he="69" />i和j表示[0,1)区间内相互独立的浮点数;(3)对近似完全卷积数据<img file="FDA0000452328160000012.GIF" wi="280" he="69" />进行傅里叶变换,得到近似数据的频域响应<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000452328160000013.GIF" wi="181" he="82" /></maths>(4)定义频域高斯低通滤波器<img file="FDA0000452328160000016.GIF" wi="426" he="74" />其中,D<sub>0</sub>为截止参数,即截止频率距离频率中心点的距离,u和v分别为矩阵元素的横坐标和纵坐标,D(u,v)是(u,v)点与频率矩形中心的距离,当D(u,v)=D<sub>0</sub>时,高斯低通滤波器由最大值1降为0.607;(5)将近似数据的频域响应<img file="FDA0000452328160000014.GIF" wi="150" he="82" />与频域高斯低通滤波器S(u,v)相乘,得到处理后回波数据的频域响应G(u,v);(6)对天线方向图h(x,y)进行傅里叶变换,得到频域响应H(u,v);(7)采用空间不变的滤波算法,利用信号功率谱sf(u,v)和噪声功率谱sn(u,v)的比值作为先验知识,计算传输函数H<sub>w</sub>(u,v):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>H</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>H</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow></mrow><mrow><msup><mrow><mo>|</mo><mi>H</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>sn</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>/</mo><mi>sf</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000452328160000015.GIF" wi="823" he="156" /></maths>其中,(·)<sup>*</sup>表示共轭,|·|<sup>2</sup>表示模的平方,sf(u,v)表示信号功率谱,sn(u,v)表示噪声功率谱;(8)用传输函数H<sub>w</sub>(u,v)乘以处理后回波数据的频域响应G(u,v),得到解卷积后的复原图像的频域响应<img file="FDA0000452328160000021.GIF" wi="179" he="81" />(9)将复原图像的频域响应<img file="FDA0000452328160000022.GIF" wi="148" he="81" />进行逆傅里叶变换,得到复原图像<img file="FDA0000452328160000023.GIF" wi="180" he="83" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>IFFT</mi><mo>[</mo><mover><mi>F</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000452328160000024.GIF" wi="527" he="103" /></maths>其中,IFFT[·]表示逆傅里叶变换。 |