CN101995573A - 一种用于超宽带穿墙雷达的目标边界成像方法 - Google Patents
一种用于超宽带穿墙雷达的目标边界成像方法 Download PDFInfo
- Publication number
- CN101995573A CN101995573A CN2009100904715A CN200910090471A CN101995573A CN 101995573 A CN101995573 A CN 101995573A CN 2009100904715 A CN2009100904715 A CN 2009100904715A CN 200910090471 A CN200910090471 A CN 200910090471A CN 101995573 A CN101995573 A CN 101995573A
- Authority
- CN
- China
- Prior art keywords
- wall
- object boundary
- ultra
- epsiv
- imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
一种用于超宽带穿墙雷达的目标边界成像方法,涉及在有限距离下超宽带穿墙雷达系统对墙后隐藏目标边界成像技术。该方法包括:建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数;输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号;利用三次样条函数对上步输出信号进行拟合插值运算,提取准波前;根据准波前集合计算曲线簇,计算曲线簇的包络线;对上步中包络线进行拟合,从而得到目标边界图像。此方法将墙体的存在对电磁波传播的影响考虑在内,可以实现对墙后隐藏目标边界准确成像,且具有较高的运算效率。
Description
技术领域
本发明涉及一种超宽带穿墙成像系统中的成像方法,具体地说涉及一种在有限距离下超宽带穿墙雷达系统对墙后隐藏目标边界成像技术。
背景技术
国内外普遍采用非相干方法和相干方法来实现超宽带穿墙成像。非相干方法采用三圆定位(Trilateration)原理来确定目标的位置,虽然硬件简单,但是需要求解一个复杂的方程组,并且存在多个目标时难以求解。相干方法一般是利用时域后向投影(Back-Projection,BP)算法成像,这种算法硬件简单,精度高,但是计算量较大,很难做到实时处理,成像效果也不尽人意。
近年来,有学者提出了基于逆边界散射变换的成像算法,该算法主要是利用了目标边界形状和接收脉冲延迟之间存在逆边界散射变换(IBST,Inverse Boundary Scattering Transform)的特性,其最大的优越性在于直接完成由准波前到目标边界的转换计算,计算时间非常短,且可以对目标边界形状清晰成像。如图1所示全向天线沿x轴方向移动扫描,天线扫描到位置(X,0)时,回波沿与目标边界垂直方向返回,根据每一个天线扫描位置处的回波可以得到发射脉冲的相应延迟t,接收回波准波前Y=ct/2,提取出所有天线回波准波前即构成了图2中数据空间,由图中几何关系可以得到BST变换,它描述了目标边界点x,y)与准波前(X,Y)之间的关系:
那么,其逆变换就解决了从接收准波前到目标边界的重建问题。
BST的逆变换(IBST)为:
因此,从接收回波中提取出目标边界对应的准波前信息进行IBST,便可以得到目标边界点(x,y),从而实现对目标边界成像。
但是在穿墙雷达应用中,由于墙体的存在,电磁波在墙体与空气交界处发生折射,使得这种成像算法中的几何关系不再成立,并且每个收/发复用天线入射角及折射角均不同,电磁波在墙体中速度也发生改变,此时使用IBST变换计算出来的目标边界存在较大的误差,边界形状也会发生一定变形。
发明内容
本发明的目的是提供一种用于超宽带穿墙雷达的目标边界成像方法,以克服公知技术中穿墙成像方法中计算效率低的缺点,
为实现上述目的,本发明考虑到穿墙成像模型中墙体对电磁波传播的影响,在提取目标回波准波前后,根据建立的墙体模型计算出准波前所对应的曲线簇的包络线,即为目标边界,由此提供了一种用于超宽带穿墙雷达的目标边界成像算法。
具体地说,本发明提供的用于超宽带穿墙雷达的目标边界成像方法,其主要步骤如下:
1)建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数;
2):输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号;
3)利用三次样条函数对步骤2中输出的高信噪比信号进行拟合插值,提取准波前;
4)根据步骤3中拟合插值后的准波前集合计算曲线簇,并计算曲线簇的包络线;
5)对步骤4中得到的包络线进行拟合,从而得到目标边界图像。
本发明克服了传统穿墙成像方法中计算效率低的缺点,实现了对隐藏在墙后目标边界的准确有效成像,满足了在城市巷战、反恐斗争、公安防暴、灾害救援等特殊场合下穿墙成像的需要。
附图说明
图1是基于逆边界散射变换成像方法的实际空间;
图2是基于逆边界散射变换成像方法的数据空间;
图3是本发明穿墙成像墙体折射几何图;
图4是本发明穿墙成像几何模型;
图5是本发明圆柱目标边界仿真成像结果;
图6是本发明圆柱目标回波提取的准波前;
图7是本发明圆柱目标回波准波前对应曲线簇与包络线示意图。
图8是本发明的计算步骤流程图。
具体实施方式
下面将结合附图和具体实施方式对本发明做进一步详细的解释。
本发明中收发复用天线阵列沿方位向紧贴墙壁放置,则可以应用本发明所阐述的目标边界成像方法对天线阵列回波信号进行成像处理,形成墙后隐藏目标边界图像(请参阅图8)。
第一步,建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数。
第二步,输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号。
第三步,利用三次样条函数对上一步输出的高信噪比信号进行拟合插值,提取准波前。
第四步,根据拟合插值后的准波前集合计算曲线簇,并计算曲线簇的包络线。
设墙壁厚度为d,相对介电常数为εr,如图3所示,天线位置坐标为(X,-d),电磁波在墙壁与空气分界面上的折射关系满足Snell折射定理:
其中,θi和θr分别为入射角和折射角,设折射点坐标为(X1,0),则由上式可得:
化简得到关于X1的四次多项式方程:
(x-X1)2[(X1-X)2+d2]-εr(X1-X)2[(x-X1)2+y2]=0
由上述方程可解出折射点位置,为了节省计算量,根据工程应用常使用的一种两层介质折射点近似确定方法有以下关系式:
设天线接收回波中发射脉冲的延迟为t,那么有:
其中,c为光速。
令M=(x-X)2,则化简上式得到:
令
则方程化简为关于M的一元二次方程:
(A+BM)2=CM+D
其解为:
另一根
超出成像区域范围,舍去。
在实际计算中,曲线簇是由离散点构成的,因此不可能所有的(xe,ye)满足等于零,只需满足即可,Δ取值范围在[0,0.1]之间,就可算出满足条件的(xe,ye)构成的包络线,即真实目标边界。
第五步,对步骤四中得到的包络线进行拟合,从而得到目标边界图像。
图5给出了圆柱目标边界仿真成像结果。图6为从天线阵列接收的圆柱目标回波信号中提取的准波前,图7为根据准波前集合计算处的曲线簇及其包络线。
理论分析和仿真实验结果都表明,本发明所述成像方法可以实现对墙后隐藏目标边界准确成像,算法效率高。
本发明提出的超宽带穿墙雷达的目标边界成像方法通过仿真实验进行了理论验证,理论分析和仿真实验结果证明了本发明的有效性。
Claims (7)
1.一种用于超宽带穿墙雷达的目标边界成像方法,其主要步骤如下:
1)建立穿墙雷达目标边界成像场景中墙体模型,设置成像所需系统参数;
2)输入阵列天线回波信号,并对输入的回波信号进行累积平均和匹配滤波处理,得到高信噪比输出信号;
3)对高信噪比信号进行拟合插值,提取准波前;
4)根据拟合插值后的准波前集合计算曲线簇,并计算曲线簇的包络线;
5)对包络线进行拟合,从而得到目标边界图像。
2.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,天线阵列沿方位向紧贴墙壁放置,并且为收发复用天线。
3.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,步骤1中的设置成像所需系统参数包括:
墙壁厚度及相对介电常数;
天线位置坐标;
电磁波在墙壁与空气分界面上的入射角和折射角,及
折射点坐标。
4.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,步骤2中,在座标(X,0)处天线接收回波中为s′(X,Y),首先去除墙壁反射波,为了提高信噪比,天线阵原始回波数据经过累积平均和匹配滤波处理后输出为s(X,Y),其中,Y为(X,-d)位置处天线的接收回波准波前Y=ct/2,t为(X,-d)位置处天线的接收回波相对发射脉冲的时间延迟,c为光速。
5.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,步骤3中是用三次样条函数对高信噪比信号进行拟合插值。
6.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,步骤4中根据拟合插值后的准波前集合计算曲线簇的计算步骤为:
化简得到关于X1的四次多项式方程
(x-X1)2[(X1-X)2+d2]-εr(X1-X)2[(x-X1)2+y2]=0,
根据两层介质折射点近似确定方法计算处折射点位置坐标
其中,(x,y)为目标边界点坐标;
6-2)根据步骤3中得到的准波前集合计算曲线簇,计算步骤为:
由天线方位向坐标X和步骤6-1中计算出的折射点位置坐标得到准波前
对于每一个成像区域范围内的y计算出其对应x,对应于所有(X,Y)的(x,y)的集合构成一簇曲线,
6-3)计算步骤6-2中得到的曲线簇的包络线,即满足条件
7.根据权利要求1所述的用于超宽带穿墙雷达的目标边界成像方法,其中,步骤6-3中计算曲线簇的包络线时,需满足Δ取值范围在[0,0.1]之间,算出满足条件的(xe,ye)构成的包络线,即目标边界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910090471A CN101995573B (zh) | 2009-08-12 | 2009-08-12 | 一种用于超宽带穿墙雷达的目标边界成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910090471A CN101995573B (zh) | 2009-08-12 | 2009-08-12 | 一种用于超宽带穿墙雷达的目标边界成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101995573A true CN101995573A (zh) | 2011-03-30 |
CN101995573B CN101995573B (zh) | 2012-09-26 |
Family
ID=43785979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910090471A Expired - Fee Related CN101995573B (zh) | 2009-08-12 | 2009-08-12 | 一种用于超宽带穿墙雷达的目标边界成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101995573B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605129A (zh) * | 2013-11-26 | 2014-02-26 | 桂林电子科技大学 | 基于椭圆包络线的超宽带穿墙雷达成像方法 |
CN103675810A (zh) * | 2013-11-13 | 2014-03-26 | 中国科学院电子学研究所 | 穿墙雷达成像的方法 |
CN103954967A (zh) * | 2014-05-19 | 2014-07-30 | 么彬 | 一种图像声纳用近场声纳图像快速成像方法 |
CN104224188A (zh) * | 2014-09-23 | 2014-12-24 | 中国科学院电子学研究所 | 用于生命探测雷达的微弱呼吸信号检测方法 |
CN105182322A (zh) * | 2015-08-31 | 2015-12-23 | 西北大学 | 基于反射信号相位差的被动式定位方法 |
CN106226765A (zh) * | 2016-09-12 | 2016-12-14 | 桂林电子科技大学 | 一种建筑物布局成像方法及系统 |
CN107576969A (zh) * | 2017-08-08 | 2018-01-12 | 中国科学院西安光学精密机械研究所 | 基于gpu并行计算的大场景隐藏目标成像系统与方法 |
CN108226861A (zh) * | 2017-11-24 | 2018-06-29 | 百年金海科技有限公司 | 超宽带定位穿墙自补偿方法及自补偿系统 |
CN110456344A (zh) * | 2019-08-13 | 2019-11-15 | 成都电科慧安科技有限公司 | 穿墙雷达成像中对墙体参数的估计方法 |
CN115372919A (zh) * | 2022-08-30 | 2022-11-22 | 中国船舶集团有限公司第七二三研究所 | 一种基于t检验的雷达目标回波模拟性能评估方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107656272B (zh) * | 2017-09-25 | 2019-08-20 | 厦门大学 | 一种在高阶时域算法下的电磁波三维逆时偏移成像方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1711844B1 (en) * | 2004-01-20 | 2009-10-14 | BAE SYSTEMS Information and Electronic Systems Integration Inc. | Multiple frequency through-the-wall motion detection and ranging using difference-based estimation technique |
US7339516B2 (en) * | 2005-08-11 | 2008-03-04 | Realtronics Corporation | Method to provide graphical representation of Sense Through The Wall (STTW) targets |
CN101441268A (zh) * | 2007-11-23 | 2009-05-27 | 电子科技大学 | 墙壁参数测试方法 |
CN101320085B (zh) * | 2008-07-21 | 2012-07-25 | 哈尔滨工业大学 | 基于后向投影算法的超宽带穿墙点目标定位成像方法 |
-
2009
- 2009-08-12 CN CN200910090471A patent/CN101995573B/zh not_active Expired - Fee Related
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675810A (zh) * | 2013-11-13 | 2014-03-26 | 中国科学院电子学研究所 | 穿墙雷达成像的方法 |
CN103675810B (zh) * | 2013-11-13 | 2016-03-02 | 中国科学院电子学研究所 | 穿墙雷达成像的方法 |
CN103605129A (zh) * | 2013-11-26 | 2014-02-26 | 桂林电子科技大学 | 基于椭圆包络线的超宽带穿墙雷达成像方法 |
CN103954967B (zh) * | 2014-05-19 | 2018-02-02 | 北京海卓同创科技有限公司 | 一种图像声纳用近场声纳图像快速成像方法 |
CN103954967A (zh) * | 2014-05-19 | 2014-07-30 | 么彬 | 一种图像声纳用近场声纳图像快速成像方法 |
CN104224188A (zh) * | 2014-09-23 | 2014-12-24 | 中国科学院电子学研究所 | 用于生命探测雷达的微弱呼吸信号检测方法 |
CN104224188B (zh) * | 2014-09-23 | 2016-06-01 | 中国科学院电子学研究所 | 用于生命探测雷达的微弱呼吸信号检测方法 |
CN105182322A (zh) * | 2015-08-31 | 2015-12-23 | 西北大学 | 基于反射信号相位差的被动式定位方法 |
CN106226765A (zh) * | 2016-09-12 | 2016-12-14 | 桂林电子科技大学 | 一种建筑物布局成像方法及系统 |
CN107576969A (zh) * | 2017-08-08 | 2018-01-12 | 中国科学院西安光学精密机械研究所 | 基于gpu并行计算的大场景隐藏目标成像系统与方法 |
CN107576969B (zh) * | 2017-08-08 | 2019-06-28 | 中国科学院西安光学精密机械研究所 | 基于gpu并行计算的大场景隐藏目标成像系统与方法 |
CN108226861A (zh) * | 2017-11-24 | 2018-06-29 | 百年金海科技有限公司 | 超宽带定位穿墙自补偿方法及自补偿系统 |
CN108226861B (zh) * | 2017-11-24 | 2021-10-29 | 郑州国测智能科技有限公司 | 超宽带定位穿墙自补偿方法及自补偿系统 |
CN110456344A (zh) * | 2019-08-13 | 2019-11-15 | 成都电科慧安科技有限公司 | 穿墙雷达成像中对墙体参数的估计方法 |
CN115372919A (zh) * | 2022-08-30 | 2022-11-22 | 中国船舶集团有限公司第七二三研究所 | 一种基于t检验的雷达目标回波模拟性能评估方法 |
CN115372919B (zh) * | 2022-08-30 | 2023-06-09 | 中国船舶集团有限公司第七二三研究所 | 一种基于t检验的雷达目标回波模拟性能评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101995573B (zh) | 2012-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101995573B (zh) | 一种用于超宽带穿墙雷达的目标边界成像方法 | |
CN102778681B (zh) | 固定发射站双基地前视合成孔径雷达成像方法 | |
CN102967858A (zh) | 雷达前视超分辨成像方法 | |
CN102313887B (zh) | 一种星机联合双基地合成孔径雷达成像方法 | |
CN108710111B (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN103777178B (zh) | 一种同步误差补偿方法、设备及系统 | |
CN103983974B (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN101581780A (zh) | 一种用于侧视层析合成孔径雷达的三维聚焦成像方法 | |
CN104020456B (zh) | 一种基于多特显点线阵成像雷达系统幅相误差校正方法 | |
CN102393518A (zh) | 一种适用于大斜视角的机载sar成像方法 | |
CN109358328B (zh) | 机动平台双基地前视sar的极坐标格式成像方法 | |
CN105137424A (zh) | 一种杂波背景下实波束扫描雷达角超分辨方法 | |
CN102520405B (zh) | 基于图像域墙体影响补偿的穿墙成像方法 | |
CN105759262A (zh) | 一种基于mimo体制的近程三维成像装置及其成像方法 | |
CN105137430A (zh) | 一种前视阵列sar的回波稀疏获取及其三维成像方法 | |
CN103293528B (zh) | 扫描雷达超分辨成像方法 | |
CN103605118B (zh) | 一种利用极地探冰雷达提取极地冰层位方法 | |
CN102129072A (zh) | 基于逆边界散射变换的超宽带雷达人体运动目标成像方法 | |
CN106772361A (zh) | 一种基于fpga的超宽带穿墙雷达成像算法的实现方法 | |
CN104007437A (zh) | 基于fpga和多dsp的sar实时成像处理方法 | |
CN103472450A (zh) | 基于压缩感知的非均匀空间构形分布式sar动目标三维成像方法 | |
CN109471097B (zh) | 一种穿墙雷达信号优化处理方法及装置 | |
CN102540184B (zh) | 探地雷达频域成像方法 | |
CN103412306B (zh) | 基于多站多波段带宽融合的超宽带相参处理方法 | |
CN103630905A (zh) | 阵列天线sar极坐标交叠子孔径成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120926 Termination date: 20130812 |