发明名称 基于联合多波位的高分辨DBS成像方法
摘要 本发明公开了一种基于联合多波位的高分辨DBS成像方法。主要解决现有DBS成像分辨率较低的问题。其实现过程为:(1)利用雷达照射地面点目标的时间计算雷达相干积累脉冲数和需要联合的波位数;(2)根据需要联合的波位数对雷达多个波位的数据进行合并;(3)对合并后的数据进行距离向的脉冲压缩和方位向的傅里叶变换,得到多普勒通道数据;(4)取出多普勒通道数据中雷达主波束范围对应的部分,并计算该对应部分的数据在地面上的坐标;(5)利用求得的地面坐标将雷达主波束范围数据标记到地面上得到图像。本发明能在飞机飞行过程通过雷达进行机械扫描,对地面进行成像,极大地提高了图像的分辨率,可用于对地面运动目标的检测和定位。
申请公布号 CN103605131B 申请公布日期 2015.09.30
申请号 CN201310647776.8 申请日期 2013.12.04
申请人 西安电子科技大学 发明人 王彤;吴晓明;吴建新;吴亿锋
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于联合多波位的高分辨DBS成像方法,包括如下步骤:(1)将雷达安装在飞机上,雷达在飞机飞行过程中进行机械扫描,且扫描波位随着时间变化,雷达共扫描N个波位,N≥2;将雷达开机工作扫描的起始波位作为第1个波位,每扫描一个波位雷达录取一次地面回波数据;(2)建立第1个波位对应的地面坐标系X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>,读取飞机在正东方向速度v<sub>E1</sub>,正北方向速度v<sub>N1</sub>,得到飞机的航向角θ<sub>1</sub>=atan(v<sub>E1</sub>/v<sub>N1</sub>);(3)建立雷达在第i个波位照射地面点目标p的几何模型;(4)利用步骤(3)中建立的模型求雷达在第i个波位照射地面点目标的时间t<sub>i</sub>,进而求得雷达相干积累脉冲数K<sub>i</sub>=floor(t<sub>i</sub>·f<sub>r</sub>),确定要联合的波位数为:B<sub>i</sub>=floor(K<sub>i</sub>/K),其中,f<sub>r</sub>为雷达发射脉冲重复频率,K为每个波位发射的脉冲数,floor代表向下取整;i为扫描波位数,其起始为1,终止为N;(5)取出第i到i+B<sub>i</sub>‑1共B<sub>i</sub>个波位雷达录取的回波数据,将取出的B<sub>i</sub>个波位的数据合并,得到波位数据矩阵A;(6)利用相关函数法估计出波位数据矩阵A的多普勒中心频率f<sub>dc</sub>,得到多普勒中心通道;(7)对上述波位数据矩阵A分别进行距离向的脉冲压缩和方位向的加权傅里叶变换,得到波位数据矩阵A的距离‑多普勒图;(8)根据步骤(6)中得到的多普勒中心通道,分别取多普勒中心通道两侧各L个通道,得到覆盖主波束范围的2L+1个多普勒通道的数据矩阵B;(9)建立第i个波位对应的地面坐标系X<sub>i</sub>O<sub>i</sub>Y<sub>i</sub>,求得多普勒通道数据矩阵B中任意一个点b在该地面坐标系X<sub>i</sub>O<sub>i</sub>Y<sub>i</sub>下的坐标(x<sub>bi</sub>,y<sub>bi</sub>)为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>bi</mi></msub><mo>=</mo><mfrac><mrow><mi>&lambda;</mi><msub><mi>R</mi><mi>b</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>i</mi></msub></mrow></mfrac><msub><mi>f</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>bi</mi></msub><mo>=</mo><msqrt><msubsup><mi>R</mi><mi>b</mi><mn>2</mn></msubsup><mo>-</mo><msup><mi>H</mi><mn>2</mn></msup><mo>-</mo><mfrac><mrow><msup><mi>&lambda;</mi><mn>2</mn></msup><msubsup><mi>R</mi><mi>b</mi><mn>2</mn></msubsup></mrow><mrow><mn>4</mn><msup><msub><mi>V</mi><mi>i</mi></msub><mn>2</mn></msup></mrow></mfrac><msubsup><mi>f</mi><mi>d</mi><mn>2</mn></msubsup></msqrt></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000698246910000021.GIF" wi="596" he="325" /></maths>其中,λ是雷达发射信号波长,R<sub>b</sub>表示雷达与点b之间的距离,f<sub>d</sub>是点b的多普勒频率,H是飞机高度,V<sub>i</sub>是飞机在第i个波位的速度;(10)读取飞机在第i个波位的正东方向速度v<sub>Ei</sub>,正北方向速度v<sub>Ni</sub>,得到飞机在该波位的航向角θ<sub>i</sub>=atan(v<sub>Ei</sub>/v<sub>Ni</sub>),求出飞机航向角在第i个波位与第1个波位之间的夹角α<sub>i</sub>:α<sub>i</sub>=θ<sub>i</sub>‑θ<sub>1</sub>=atan(v<sub>Ei</sub>/v<sub>Ni</sub>)‑atan(v<sub>E1</sub>/v<sub>N1</sub>);(11)求得点b在第1个波位坐标系X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>下的坐标(x<sub>b1</sub>,y<sub>b1</sub>)为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mo>=</mo><mi>cos</mi><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>x</mi><mi>bi</mi></msub><mo>-</mo><mi>sin</mi><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>y</mi><mi>bi</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mo>=</mo><mi>sin</mi><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>x</mi><mi>bi</mi></msub><mo>+</mo><mi>cos</mi><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>y</mi><mi>bi</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000698246910000022.GIF" wi="553" he="182" /></maths>(12)利用步骤(11)的公式,求出多普勒通道数据矩阵B中所有点在第1个波位对应的地面坐标系X<sub>1</sub>O<sub>1</sub>Y<sub>1</sub>下的坐标,根据求出的地面坐标,将数据矩阵B中的所有点标记到其坐标对应的地面位置上;(13)令i=i+1,重复步骤(3)~(12),直到处理完所有雷达扫描波位的数据,即i=N时,显示出全部波位数据在地面的DBS图像。
地址 710071 陕西省西安市太白南路2号