发明名称 一种用于超宽带穿墙雷达的目标边界成像方法
摘要 一种用于超宽带穿墙雷达的目标边界成像方法,涉及在有限距离下超宽带穿墙雷达系统对墙后隐藏目标边界成像技术。该方法包括:建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数;输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号;利用三次样条函数对上步输出信号进行拟合插值运算,提取准波前;根据准波前集合计算曲线簇,计算曲线簇的包络线;对上步中包络线进行拟合,从而得到目标边界图像。此方法将墙体的存在对电磁波传播的影响考虑在内,可以实现对墙后隐藏目标边界准确成像,且具有较高的运算效率。
申请公布号 CN101995573B 申请公布日期 2012.09.26
申请号 CN200910090471.5 申请日期 2009.08.12
申请人 中国科学院电子学研究所 发明人 黄琼;陈洁;孟升卫;方广有
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 周长兴
主权项 1.一种用于超宽带穿墙雷达的目标边界成像方法,其主要步骤如下:步骤1:建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数;步骤2:输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号;步骤3:对高信噪比输出信号进行拟合插值,提取准波前;步骤4:根据拟合插值后的准波前集合计算曲线簇,首先根据电磁波在墙壁与空气分界面上的折射关系满足Snell折射定理,计算(X,-d)位置处天线的发射/接收的电磁波在墙壁与空气分界面的折射点坐标,根据Snell定理<maths num="0001"><![CDATA[<math><mrow><mi>sin</mi><msub><mi>&theta;</mi><mi>r</mi></msub><mo>/</mo><mi>sin</mi><msub><mi>&theta;</mi><mi>i</mi></msub><mo>=</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub></msqrt></mrow></math>]]></maths>得到<maths num="0002"><![CDATA[<math><mrow><mo>[</mo><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>/</mo><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></msqrt><mo>]</mo><mo>/</mo><mo>[</mo><mrow><mo>(</mo><msub><mi>X</mi><mn>1</mn></msub><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mo>/</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>1</mn></msub><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>d</mi><mn>2</mn></msup></msqrt><mo>]</mo><mo>=</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub></msqrt><mo>,</mo></mrow></math>]]></maths>其中,θ<sub>i</sub>和θ<sub>r</sub>分别为入射角和折射角,墙壁厚度为d,墙壁相对介电常数为ε<sub>r</sub>,化简得到关于X<sub>1</sub>的四次多项式方程(x-X<sub>1</sub>)<sup>2</sup>[(X<sub>1</sub>-X)<sup>2</sup>+d<sup>2</sup>]-ε<sub>r</sub>(X<sub>1</sub>-X)<sup>2</sup>[(x-X<sub>1</sub>)<sup>2</sup>+y<sup>2</sup>]=0,根据两层介质折射点近似确定方法计算折射点位置坐标<maths num="0003"><![CDATA[<math><mrow><mrow><msub><mi>X</mi><mn>1</mn></msub><mo>=</mo><mi>X</mi><mo>+</mo><mi>d</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mo>/</mo><mo>[</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub></msqrt><mrow><mo>(</mo><mi>y</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,(x,y)为目标边界点坐标;由天线方位向坐标X和折射点位置坐标得到准波前<maths num="0004"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>Y</mi><mo>=</mo><mi>ct</mi><mo>/</mo><mn>2</mn><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>X</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></msqrt><mo>+</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub><mo>[</mo><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>1</mn></msub><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>d</mi><mn>2</mn></msup><mo>]</mo></msqrt></mtd></mtr><mtr><mtd><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>{</mo><mn>1</mn><mo>-</mo><mi>d</mi><mo>/</mo><mo>[</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub></msqrt><mrow><mo>(</mo><mi>y</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup><mo>+</mo><mi>d</mi><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub><mo>+</mo><msup><mrow><mo>[</mo><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mi>y</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></msqrt></msqrt></mtd></mtr></mtable></mfenced></math>]]></maths>其中,t为天线接收回波中发射脉冲的延迟,c为光速,对于每一个成像区域范围内的y计算出其对应x,对应于所有(X,Y)的(x,y)的集合构成一簇曲线,<maths num="0005"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>{</mo><mn>1</mn><mo>-</mo><mi>d</mi><mo>/</mo><mo>[</mo><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub></msqrt><mrow><mo>(</mo><mi>y</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></msqrt><mo>+</mo><mi>d</mi><msqrt><msub><mi>&epsiv;</mi><mi>r</mi></msub><mo>+</mo><msup><mrow><mo>[</mo><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>X</mi><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mi>y</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><mi>Y</mi><mo>=</mo><mn>0</mn><mo>;</mo></mrow></math>]]></maths>步骤5:计算曲线簇的包络线,即满足条件{(x<sub>e</sub>,y<sub>e</sub>)|F(x<sub>e</sub>,y<sub>e</sub>,X,Y)=0且<img file="FDA00001726610700016.GIF" wi="447" he="58" />的点为目标边界点的集合,其中Δ取值范围在[0,0.1]之间,满足条件的点(x<sub>e</sub>,y<sub>e</sub>)构成的包络线,即目标边界;步骤6:对包络线进行拟合,从而得到目标边界图像。
地址 100080 北京市海淀区北四环西路19号
您可能感兴趣的专利