发明名称 基于非完全数据解卷积的实孔径前视成像方法
摘要 本发明公开了一种基于非完全数据解卷积的实孔径前视成像方法,主要解决飞行轨迹正前方场景的成像问题。其处理过程为:1)对地面监视区域进行顺序重叠扫描,获取雷达回波数据;2)利用线性插值扩展雷达回波数据,得到近似完全卷积数据;3)计算近似完全卷积数据的频域响应;4)定义频域高斯低通滤波器;5)将近似数据的频域响应与频域高斯低通滤波器相乘,得到处理后回波数据的频域响应;6)计算天线方向图的频域响应;7)采用空间不变的滤波算法,计算传输函数;8)计算解卷积后复原图像的频域响应;9)求取复原图像。本发明具有提高成像质量和检测速度的优点,可用于机载雷达监视系统实现航迹线方向地面场景成像的数据处理领域。
申请公布号 CN102608597B 申请公布日期 2014.07.23
申请号 CN201210072515.3 申请日期 2012.03.19
申请人 西安电子科技大学 发明人 杨志伟;廖桂生;杨凯新;刘笑菲;李军;刘志凌;杜文韬
分类号 G01S13/89(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于非完全数据解卷积的实孔径前视成像方法,包括如下步骤:(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[·]表示逆傅里叶变换。
地址 710071 陕西省西安市太白南路2号