CN102520405B - 基于图像域墙体影响补偿的穿墙成像方法 - Google Patents
基于图像域墙体影响补偿的穿墙成像方法 Download PDFInfo
- Publication number
- CN102520405B CN102520405B CN 201110425730 CN201110425730A CN102520405B CN 102520405 B CN102520405 B CN 102520405B CN 201110425730 CN201110425730 CN 201110425730 CN 201110425730 A CN201110425730 A CN 201110425730A CN 102520405 B CN102520405 B CN 102520405B
- Authority
- CN
- China
- Prior art keywords
- wall
- compensation
- imaging
- image
- image area
- 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.)
- Active
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种基于图像域墙体影响补偿的穿墙成像方法。基本思路是:首先,通过传统成像算法对回波数据进行成像处理。然后在无墙体先验知识的情况下,利用获取的图像及其图像域滤波结果确定补偿因子中的最佳参数。最后利用该补偿因子对成像结果进行墙体影响的补偿。具体包括下述步骤:第一步,确定成像坐标系;第二步,对回波数据进行初始成像处理;第三步,计算图像域补偿因子;第四步,进行图像域墙体补偿处理。本发明仅需对墙体补偿进行一次成像处理便可同时完成多目标的补偿处理,从而满足实际应用高计算效率要求。
Description
技术领域
本发明属于超宽带合成孔径雷达成像技术领域,特别是一种基于图像域墙体影响补偿的成像方法。
背景技术
超宽带合成孔径雷达将超宽带技术与合成孔径技术相结合,发射的电磁波具有低频穿透以及优异的方位向分辨能力,实现对墙后目标甚至建筑物结构的高分辨成像,极大增强了目标信息获取能力。由于墙体遮挡的影响,电磁波在穿透墙体时后产生能量衰减、反射、折射和速度变慢等现象,利用天线接收回波进行成像就会出现定位偏差、图像散焦甚至产生虚假目标等问题,因而在穿墙成像中,需要进行墙体影响的补偿处理。
墙体影响的补偿是超宽带合成孔径雷达穿墙成像中的关键技术。为了补偿墙体造成的影响,常用的处理方法是利用已知的墙体电磁参数对天线接收回波进行补偿。而在实际穿墙中,墙体电磁参数一般是未知的,在这种情况下,穿墙成像补偿的关键一般依靠墙体参数的估计进行处理。目前关于墙体影响的补偿,通常在回波域进行,该方法将墙体对天线接收回波的影响在回波域进行补偿,补偿后再利用成像算法进行成像。该方法存在两个主要问题。一是不能对墙后多目标同时补偿。因为当墙后有多个目标时,不同目标的入射角是不同的,即不同目标对应的回波域补偿因子是不同的。二是该方法计算量很大。因为在墙体厚度确定中,需要进行多次成像,计算量随之大幅度增加。
发明内容
本发明提出一种基于图像域墙体影响补偿的穿墙成像方法,可实现对多个目标同时补偿后成像,同时具有运算量小、成像精度高的特点。
本发明的基本思路是:首先,通过传统成像算法对回波数据进行成像处理。然后在无墙体先验知识的情况下,利用获取的图像及其图像域滤波结果确定补偿因子中的最佳参数。最后利用该补偿因子对成像结果进行墙体影响的补偿。
本发明的技术方案包括以下处理步骤:
第一步,确定成像坐标系
确定成像坐标系:以超宽带合成孔径雷达形成的合成孔径中心点为坐标原点O,以合成孔径所在直线为X轴,在合成孔径中心点处做垂直于Y轴并垂直指向墙体的直线,记为Y轴。
第二步,对回波数据进行初始成像处理
回波数据的初始成像处理是后续补偿的基础。记s(t,u)为天线接收回波,其中t表示接收时间,u表示天线接收数据时在成像坐标系对应的横坐标。采用已知成像算法对接收到的回波数据进行初始成像处理,得到初始成像结果g(x,y),g(x,y)表示在成像坐标系中坐标为(x,y)处的像素点的值。
第三步,计算图像域补偿因子
首先利用合成孔径中心点的接收回波数据,计算墙体外表面回波与内表面回波的延时差估计值Δτ′W。
Δτ′W=t″w-t′w (1)
其中,t′w、t″w分别为墙体外表面与墙体内表面的响应时间,即对应墙体回波数据的第一个和第二个包络峰值。此时得到墙体等效估计厚度为
其中c为自由空间下光速。
步骤②:构造含参变量D的图像域补偿因子
对初始成像结果g(x,y)进二维傅里叶变换,则得到G(kx,ky),其中kx和ky分别为x和y的傅里叶变换对偶变量,称为方位向波数和距离向波数。构造图像域补偿矩阵F(kx,ky)|D为
其中θi=arctan(|kx/ky|),D为墙体厚度。
步骤③:估计墙体厚度
g′(x,y)|D=FT-1[G(kx,ky)F(kx,ky)|D] (4)
式中,FT-1表示逆傅里叶变换。
依次计算各补偿后的图像中目标主瓣在方位向的3dB宽度作为目标方位向分辨率。取目标方位分辨率最优时参变量D的取值d为墙体厚度的估计值。
第四步,进行图像域墙体补偿处理
由于不同的目标具有相同的图像域补偿因子,因而本发明可以同时对多目标进行补偿处理。
I(x,y)=g′(x,y)|D=d (5)
I(x,y)即为补偿后的最终成像结果。
本发明的有益效果:本发明通过提出基于图像域墙体影响补偿的穿墙成像方法,可在无墙体参数先验知识的情况下,获得对墙体影响的高精度补偿后图像。同时该发明仅需对墙体补偿进行一次成像处理并可同时完成多目标的补偿处理,从而满足实际应用高计算效率要求。
附图说明
图1为本发明处理流程示意图;
图2为本发明验证仿真的三维场景图;
图3为初始成像结果;
图4为一维回波信号包络图像;
图5为目标方位向分辨率与墙体厚度的对应关系图像;
图6为本发明补偿后的成像结果。
具体实施方式
基于图像域墙体影响补偿的穿墙成像方法分为四个步骤,如图1所示。下面结合XFDTD电磁仿真软件获取的数据对本发明做进一步解释。
在仿真中,墙体厚度0.2m,相对介电常数为3.24。天线距离墙体垂直距离8m,目标放置与墙体后方2m处。输入激励源选择一阶微分高斯信号,整个计算区域利用完全匹配层截断,以满足吸收边界条件。
第一步:确定成像坐标系。图2左侧从上至下依次为场景的XY、YZ和ZX的二维侧视图,右侧为场景的三维图。以天线孔径点所在直线为X轴,以垂直于X轴并垂直指向墙体的直线为Y轴,坐标原点为合成孔径中心点,确定成像坐标写如图2所示。
第二步:对回波数据进行初始成像处理。如采用BP(Back Project,后向投影)成像算法对接收回波数据进行成像处理,成像结果见图3。
在该步处理中,同样可以采用其他成像算法进行初始成像处理,如ωK算法等。
从初始成像结果可以看到,目标发生了定位误差和散焦,目标最大值出现在(0,10.15)处,距离向定位误差为0.15m(理论定位误差为0.16m)。
第三步,计算图像域补偿因子。
步骤①,计算墙体等效估计厚度在合成孔径中心点的接收回波数据中,墙体外表面响应的包络峰值出现在66.83ns,墙体内表面的包络峰值出现在69.24ns,此时墙体外表面和内表面的延时差为2.41ns。利用公式2得到墙体等效厚度为0.36m。
步骤②,构造含参变量D的图像域补偿因子。
步骤③,取D从0至0.36m进行遍历,用对应的图像补偿因子对初始成像结果进行补偿,依次计算各补偿结果中目标主瓣在方位向的3dB宽度作为目标方位向分辨率。得到不同D值与目标方位向分辨率对应结果如图5所示。从图中可以得到当D=0.2m时,目标方位向分辨率取得最小值,此时目标方位向分辨率最优。
第四步,进行图像域墙体补偿处理。
利用第三步中得到的图像域补偿因子进行补偿处理,得到最终成像结果如图6所示。从图中可以发现,补偿后的目标最大值位于(0,10)处,与仿真中目标实际所处位置完全吻合,有效的实现了目标定位误差和散焦的校正。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于图像域墙体影响补偿的穿墙成像方法,其特征在于,包括下述步骤:
第一步,确定成像坐标系;
以超宽带合成孔径雷达形成的合成孔径中心点为坐标原点O,以合成孔径所在直线为X轴,在合成孔径中心点处做垂直于X轴并垂直指向墙体的直线,记为Y轴;
第二步,对回波数据进行初始成像处理;
记s(t,u)为天线接收回波,其中t表示接收时间,u表示天线接收数据时在成像坐标系对应的横坐标;对接收到的回波数据进行初始成像处理,得到初始成像结果g(x,y);
第三步,计算图像域补偿因子;
首先利用合成孔径中心点的接收回波数据,计算墙体外表面回波与内表面回波的延时差估计值△τ′W;
△τ′W=t′′w-t′w
其中,t′w、t′′w分别为墙体外表面与墙体内表面的响应时间,即对应墙体回波数据的第一个和第二个包络峰值;此时得到墙体等效估计厚度为
其中c为自由空间下光速;
步骤②:构造含参变量D的图像域补偿因子;
对初始成像结果g(x,y)进行二维傅里叶变换,则得到G(kx,ky),其中kx和ky分别为x和y的傅里叶变换对偶变量,x和y分别为方位向坐标和距离向坐标;构造图像域补偿矩阵F(kx,ky)|D为
其中θi=arctan(|kx/ky|),D为墙体厚度;
步骤③:估计墙体厚度;
g′(x,y)|D=FT-1[G(kx,ky)F(kx,ky)|D]
式中,FT-1表示逆傅里叶变换;
依次计算各补偿后的图像中目标主瓣在方位向的3dB宽度作为目标方位向分辨率;取目标方位分辨率最优时参变量D的取值d为墙体厚度的估计值;
第四步,进行图像域墙体补偿处理;
得到进行图像域墙体补偿处理后的成像结果为:
I(x,y)=g′(x,y)|D=d。
2.根据权利要求1所述的基于图像域墙体影响补偿的穿墙成像方法,其特征在于,采用后向投影成像算法对接收回波数据进行成像处理。
3.根据权利要求1所述的基于图像域墙体影响补偿的穿墙成像方法,其特征在于,采用ωK算法对接收回波数据进行成像处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110425730 CN102520405B (zh) | 2011-12-16 | 2011-12-16 | 基于图像域墙体影响补偿的穿墙成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110425730 CN102520405B (zh) | 2011-12-16 | 2011-12-16 | 基于图像域墙体影响补偿的穿墙成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102520405A CN102520405A (zh) | 2012-06-27 |
CN102520405B true CN102520405B (zh) | 2013-08-28 |
Family
ID=46291388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110425730 Active CN102520405B (zh) | 2011-12-16 | 2011-12-16 | 基于图像域墙体影响补偿的穿墙成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102520405B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278818B (zh) * | 2013-06-13 | 2014-10-08 | 中国人民解放军国防科学技术大学 | 基于相位统计参数加权的穿墙雷达成像方法 |
CN103399302B (zh) * | 2013-06-13 | 2014-12-10 | 中国人民解放军国防科学技术大学 | 基于穿墙阵列空间对称性的墙体杂波抑制方法 |
CN103675810B (zh) * | 2013-11-13 | 2016-03-02 | 中国科学院电子学研究所 | 穿墙雷达成像的方法 |
CN103969647A (zh) * | 2014-05-29 | 2014-08-06 | 电子科技大学 | 基于图像域滤波的稀疏阵列建筑布局成像方法 |
CN104391285A (zh) * | 2014-12-08 | 2015-03-04 | 中国人民解放军国防科学技术大学 | 基于图像域的低频超宽带sar自聚焦方法 |
CN104502911B (zh) * | 2014-12-25 | 2017-04-26 | 湖南华诺星空电子技术有限公司 | 一种穿墙成像雷达的墙壁参数估计方法 |
CN106680809B (zh) * | 2016-12-27 | 2019-02-12 | 中国人民解放军国防科学技术大学 | 一种穿墙雷达自聚焦压缩感知成像方法 |
CN107300694B (zh) * | 2017-06-14 | 2020-04-17 | 南京信息工程大学 | 一种基于电磁波透射系数的未知墙体参数估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7339516B2 (en) * | 2005-08-11 | 2008-03-04 | Realtronics Corporation | Method to provide graphical representation of Sense Through The Wall (STTW) targets |
CN101320085A (zh) * | 2008-07-21 | 2008-12-10 | 哈尔滨工业大学 | 基于后向投影算法的超宽带穿墙点目标定位成像方法 |
CN101441268A (zh) * | 2007-11-23 | 2009-05-27 | 电子科技大学 | 墙壁参数测试方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7142153B2 (en) * | 2004-06-08 | 2006-11-28 | Raytheon Company | Short pulse/stepped frequency radar system |
-
2011
- 2011-12-16 CN CN 201110425730 patent/CN102520405B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 电子科技大学 | 墙壁参数测试方法 |
CN101320085A (zh) * | 2008-07-21 | 2008-12-10 | 哈尔滨工业大学 | 基于后向投影算法的超宽带穿墙点目标定位成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102520405A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102520405B (zh) | 基于图像域墙体影响补偿的穿墙成像方法 | |
CN110426707B (zh) | 基于轨道角动量的涡旋sar成像方法及成像系统 | |
CN102393518B (zh) | 一种适用于大斜视角的机载sar成像方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN109541548B (zh) | 一种基于匹配场的空气声呐定位方法 | |
CN102749620B (zh) | 一种弹载/机载雷达单脉冲前视成像处理方法 | |
CN109581388B (zh) | 一种实时三维成像声纳的近场宽视角波束形成方法 | |
CN113281727B (zh) | 一种基于水平线列阵的输出增强的波束形成方法及其系统 | |
CN108845325A (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
CN104833972A (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN103018739A (zh) | 一种校正多通道幅相误差的微波三维成像方法 | |
CN104698431A (zh) | 模糊分量空间角的估计及多通道sar方位解模糊方法 | |
CN108226925A (zh) | 一种适用于弹载大前斜视时变参数sar的重叠子孔径算法 | |
CN103487796B (zh) | 一种利用水声信道统计不变特征实现被动测距的方法 | |
CN106950531A (zh) | 基于频域时间反转的雷电多源定向方法 | |
CN106125075B (zh) | 一种双基地前视合成孔径雷达的运动误差估计方法 | |
CN104182942B (zh) | Sar图像方位模糊抑制方法 | |
CN111551942A (zh) | 一种基于解卷积算法的水下自主航行器回坞方法 | |
CN103513238A (zh) | 一种规整化最小二乘子空间相交的目标方位测向方法 | |
CN103809164B (zh) | 一种重构参考通道的后多普勒最优杂波抑制方法 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
CN103513249B (zh) | 一种宽带相干模基信号处理方法及系统 | |
CN113219458B (zh) | 双基地合成孔径雷达盲定位方法 | |
CN113126029B (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 |
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 |