发明名称 一种适用于穿墙雷达成像的多面墙体快速补偿方法
摘要 本发明提供一种适用于穿墙雷达成像的多面墙体快速补偿方法,基于已知或估计的墙体参数,考虑天线贴墙、多面平行墙体建筑结构内隐蔽目标成像问题,假设这多面墙体具有相同相对介电常数和不同厚度,首先,根据折射定律,将多面墙体平行移到第一面墙体边缘,多面墙体结合形成一个等效的单面墙体,同时保证成像聚焦延迟不变。通过上述等效处理,将多面墙体电磁传播涉及的多折射问题,等效为传统单面墙体单折射问题,简化了问题难度。然后,基于最短时间法搜索等效墙体上折射点位置,计算传播路径延迟。最后,利用后向投影成像算法得到多面墙体补偿后目标聚焦图像,准确补偿目标位置误差的优势,具有很强的实用性。
申请公布号 CN103995256B 申请公布日期 2016.03.30
申请号 CN201410232091.1 申请日期 2014.05.29
申请人 电子科技大学 发明人 崔国龙;张鹏;刘剑刚;姚雪;孔令讲;胡飞;杨晓波;殷光强;张征宇
分类号 G01S7/40(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 邹裕蓉
主权项 一种适用于穿墙雷达成像的多面墙体快速补偿方法,其特征在于,包括以下步骤:1)探测步骤:穿墙雷达利用一个单基天线对多面平行墙后的目标进行合成孔径成像探测,天线移动过程中始终处于同一高度且紧靠墙体,设天线移动合成N元的等效均匀线阵,N元的等效均匀线阵接收的回波信号为[s<sub>1</sub>(t),s<sub>2</sub>(t),……,s<sub>N</sub>(t)]<sup>T</sup>;2)虚拟折射点位置估计步骤:取成像区域内一像素点z<sub>h</sub>,在探测范围内遍历像素点z<sub>h</sub>与第k个位置天线之间的等效传播路径中的虚拟折射点E<sub>k</sub>坐标为:<img file="FDA0000886299460000011.GIF" wi="247" he="135" />其中,L为多面墙的总数,d<sub>i</sub>为第i面墙的厚度;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>k</mi><mo>&prime;</mo></msubsup><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>x</mi></munder><mrow><mo>&lsqb;</mo><mrow><mfrac><mrow><mn>2</mn><msqrt><mrow><msup><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>1</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mi>L</mi></munderover><msub><mi>d</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow><mrow><mi>c</mi><mo>/</mo><msqrt><mi>&epsiv;</mi></msqrt></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><msqrt><mrow><msup><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mi>L</mi></munderover><msub><mi>d</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow><mi>c</mi></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000886299460000012.GIF" wi="1444" he="221" /></maths>其中,ε为墙体的相对介电常数,c为光速,arg min表示使目标函数取最小值时的虚拟折射点横坐标变量x的值,(x<sub>h</sub>,y<sub>h</sub>)为像素点z<sub>h</sub>的坐标,k=1,2,…,N;3)聚焦延迟计算步骤:第k个位置天线A<sub>k</sub>(x<sub>k</sub>,‑d<sub>1</sub>)与像素z<sub>h</sub>(x<sub>h</sub>,y<sub>h</sub>)之间的聚焦延时τ<sub>hk</sub>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mrow><mi>h</mi><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msqrt><mrow><msup><mrow><mo>(</mo><mrow><msubsup><mi>x</mi><mi>k</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>x</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>1</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mi>L</mi></munderover><msub><mi>d</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow><mrow><mi>c</mi><mo>/</mo><msqrt><mi>&epsiv;</mi></msqrt></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><msqrt><mrow><msup><mrow><mo>(</mo><mrow><msubsup><mi>x</mi><mi>k</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>x</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mi>L</mi></munderover><msub><mi>d</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>h</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow><mi>c</mi></mfrac></mrow>]]></math><img file="FDA0000886299460000013.GIF" wi="1276" he="222" /></maths>4)投影成像补偿步骤:利用接收到的回波数据以及聚焦延时进行合成孔径成像:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>I</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>h</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>s</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><msub><mi>&tau;</mi><mrow><mi>h</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow></msub></mrow>]]></math><img file="FDA0000886299460000014.GIF" wi="527" he="94" /></maths>其中,︱表示参数赋值,I(z<sub>h</sub>)表示经聚焦延迟补偿后像素点z<sub>h</sub>的像素值,判读成像区域内所有像素点是否完成投影成像补偿,如否,则返回步骤2),如是,成像区域内所有像素点遍历完毕,得到多面墙体补偿后的整幅图像。
地址 611731 四川省成都市高新区(西区)西源大道2006号