CN103675810B - 穿墙雷达成像的方法 - Google Patents
穿墙雷达成像的方法 Download PDFInfo
- Publication number
- CN103675810B CN103675810B CN201310572048.5A CN201310572048A CN103675810B CN 103675810 B CN103675810 B CN 103675810B CN 201310572048 A CN201310572048 A CN 201310572048A CN 103675810 B CN103675810 B CN 103675810B
- Authority
- CN
- China
- Prior art keywords
- echo data
- wall
- domain echo
- space
- wave number
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005259 measurement Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 2
- 230000010363 phase shift Effects 0.000 claims description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 12
- 230000008030 elimination Effects 0.000 abstract description 4
- 238000003379 elimination reaction Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 abstract description 4
- 238000001228 spectrum Methods 0.000 abstract description 4
- 230000003993 interaction Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 7
- 238000005070 sampling Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2928—Random or non-synchronous interference pulse cancellers
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种穿墙雷达成像的方法。该方法在雷达成像处理过程中,对频谱数据乘以墙体补偿因子,补偿因子使用了精确的物理模型,考虑了电磁波与墙体之间的相互作用,从而消除了墙体对成像位置的影响,保证了成像准确性。此外,由于墙体影响的消除在频域进行,只涉及乘法运算,并且补偿因子计算无需求解非线性方程,因此计算量较小。
Description
技术领域
本发明涉及雷达技术领域,尤其涉及一种穿墙雷达成像的方法。
背景技术
穿墙成像(TWI)主要用来实现对隐藏在墙体或障碍物后面的目标进行检测或定位。近年来,穿墙成像已在生命检测,反恐维稳和灾害援救等方面有着广泛的应用,这些应用对穿墙成像算法的准确性和高效性提出了客观要求。因此,高效、准确的穿墙成像算法的研究具有重要的应用价值,已成为该领域研究的热点。
穿墙成像是通过发射电磁波,对接收到的回波数据进行处理而成像的。电磁波在穿透墙体时会发生折射和多次散射。电磁波与墙体的这种相互作用会造成成像结果产生偏差,即目标位置发生偏移、几何形状发生畸变等。因此,必须去除墙体对成像结果的影响。去除墙体影响的方法不同,成像的时效性和准确性也不同。
在实现本发明的过程中,申请人发现现有技术穿墙雷达成像方法不能有效补偿墙体对成像位置和精度的影响。此外,由于算法较为繁琐,造成成像速度较低,时效性较差。
发明内容
(一)要解决的技术问题
鉴于上述技术问题,本发明提供了一种穿墙雷达成像的方法,以补偿墙体对成像位置和精度的影响,并且尽可能不增加算法的复杂度。
(二)技术方案
根据本发明的一个方面,提供了一种穿墙雷达成像的方法。该方法包括:步骤A,在墙体外发射雷达波,并接收该雷达波在z=0处的时域回波数据;步骤B,将时域回波数据变换到频域,得到频域回波数据;步骤C,将频域回波数据进行相位偏移,得到整个成像区域的频域回波数据;步骤D,对整个成像区域的频域回波数据乘以一与墙体影响相关的复数因子;步骤E,对乘以复数因子后的频域回波数据,将频域回波数据中的变量ω转换到变量kz,得到波数域回波数据,其中,ω为时间t对应的角频率,kz为空域变量z在波数域对应的波数变量;步骤F,对波数域回波数据关于kx项求逆傅里叶变换得到空间波数域回波数据,其中,kx为空域变量x在波数域对应的波数变量;以及步骤G,对空间-波数域回波数据对其关于kz项变换到空时域,从而得到消除墙体影响后的雷达成像结果。
(三)有益效果
从上述技术方案可以看出,本发明穿墙雷达成像的方法具有以下有益效果:
(1)成像处理过程中,对频谱数据乘以墙体补偿因子,补偿因子使用了精确的物理模型,考虑了电磁波与墙体之间的相互作用,从而消除了墙体对成像位置的影响,保证了成像准确性;
(2)由于墙体影响的消除在频域进行,只涉及乘法运算,并且补偿因子计算无需求解非线性方程,因此计算量较小;
(3)使用NUFFT方法实现kz方向非均匀采样波数域据变换到空域,计算复杂度为O(Nlog2N)。使用插值再求FFT方法,计算复杂度为O(Nlog2N),但插值会影响精度。使用求和方法,计算复杂度为O(N2)。本发明使用NUFFT方法既能保证成像精度,又能兼顾计算效率。
附图说明
图1A为根据本发明实施例一相应场景二维空间三层媒质的几何结构示意图;
图1B为根据本发明实施例一针对二维空间三层媒质的穿墙雷达成像方法的流程图;
图2为二维空间,三层媒质点目标成像,未消除墙体影响,NUFFT处理,成像结果;
图3为二维空间,三层媒质点目标成像,采用本实施例方法消除墙体影响后,NUFFT处理,成像结果;
图4为二维空间,三层媒质点目标成像,消除墙体影响后,插值加FFT处理,成像结果;
图5A为根据本发明实施例二相应场景三维空间三层媒质的几何结构示意图;
图5B为根据本发明实施例二针对三维空间三层媒质的穿墙雷达成像方法的流程图;
图6为三维空间,三层媒质点目标成像,未消除墙体影响,NUFFT处理,成像结果;
图7为三维空间,三层媒质点目标成像,消除墙体影响后,NUFFT处理,成像结果。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。
本发明穿墙雷达成像方法中,对频谱数据乘以一个复数因子,补偿墙体对成像位置和几何形状造成的误差,并用FFT变换均匀间隔的波数域数据(kx,ky)到空间域(x,y),NUFFT变换不均匀间隔kz空间波数域回波数据到空域(z方向),之后再进行成像。
实施例一
在本发明的第一个示例性实施例中,提供了一种针对二维空间三层媒质的穿墙雷达成像方法。请参照示意图1,本实施例穿墙雷达成像方法包括:
步骤S102,在墙体外发射雷达波,并接收该雷达波在测量线z=0处的时域回波数据d′1(x,z=0,t);
如图1所示,T为发射电磁波雷达装置,R为接收时域回波数据的雷达装置,两者均位于沿墙体外表面。至少接收时域回波数据的雷达装置沿与该墙体外表面平行的测线移动。其中,x和z为二维空间的两个坐标轴。x坐标轴的正方向为沿测量线向外的方向,z坐标轴的正方向为垂直于测量线向测量目标的方向。
步骤S104,将时域回波数据d′1(x,z=0,t)变换到频域中,变换后的频域回波数据为D1(kx,z=0,ω);
本步骤中,将时域数据转换为频域数据采用的为本领域通用的FFT的方法,此处不再赘述。其中,kx为空域变量x在波数域对应的波数变量,ω为时间t对应的角频率。
步骤S106,按照以下公式对频域回波数据D1(kx,z=0,ω)进行相位偏移,得到整个成像区域的频域回波数据D1(kx,z,ω):
其中,kx为空域变量x在波数域对应的波数变量,kz为空域变量z在波数域对应的波数变量。由于本实施例的场景是二维空间数据,所以ky=0,kx,kz满足频散关系c为光速。
步骤S108,对整个成像区域的频域回波数据D1(kx,z,ω)乘以一个复数因子以补偿墙体影响,得到频域回波数据:
其中为分层媒质的传输系数,||·||表示取模。(·)*表示取共轭。d1为测线到第一层介质与第二层介质交界处的厚度,d2为测线到第二层介质与第三层介质交界面之间的厚度,为菲涅尔反射系数, k1为第一层介质对应的波数,k2为第二层介质对应的波数。本实施例中,第一层介质为空气;所述第二层介质为墙体;所述第三层介质为空气。
步骤S110,对进行墙体影响补偿后的频域回波数据D′1(kx,z,ω),由Stolt变换由频域回波数据D′1(kx,z,ω)得到波数域数据
本步骤中,ky=0,因此
步骤S112,按照以下公式对其关于kx项求逆傅里叶变换得到空间波数域回波数据
本步骤中,由于kx为均匀采样间隔,因此直接对关于kx项求逆傅里叶变换即可得到空间波数域回波数据D1(x,z,kz)。
步骤S114,对空间-波数域数据D1(x,z,kz)按照以下公式对其关于kz项采用NUFFT求逆变换到空时域d′1(x,z,t=0),从而得到消除墙体影响后的雷达成像结果:
本步骤中,由于kz为非均匀采样间隔,因此需要对数据D1(x,z,kz)关于kz项采用NUFFT求逆变换,才可以得到空时域数据d′1(x,z,t=0)并进行雷达成像。
为验证本实施例方法的效果,申请人进行了对比试验。图2为二维空间,三层媒质点目标成像,采用现有技术方法未消除墙体影响,NUFFT处理,成像结果。图3为二维空间,三层媒质点目标成像,采用本实施例方法消除墙体影响后,NUFFT处理,成像结果。对比图2与图3可知,消除墙体影响后,两维空间三层媒质点目标成像结果比未消除墙体影响的成像结果更清晰,聚焦作用更明显。
此外,需要说明的是,本实施例中,步骤S114利用NUFFT对kz方向数据求逆变换得到空时域数据,除此之外,本发明该步骤也可采用插值加FFT的传统方法对kz方向数据求逆变换得到空时域数据。利用NUFFT对kz方向数据求逆变换得到穿墙雷达成像结果见图3,利用插值加FFT对kz方向数据求逆变换得到穿墙雷达成像结果见图4。图3与图4相比,运用NUFFT求kz方向数据逆变换得到的成像结果比用插值加FFT求kz方向数据逆变换成像结果精度更高,并且计算复杂度均为O(Nlog2N)。
实施例2
在本发明的第二个示例性实施例中,提供了一种针对三维空间三层媒质的穿墙雷达成像方法。请参照示意图5A,本实施例穿墙雷达成像方法包括:
步骤S202,在墙体外侧用雷达发射电磁波,获取该电磁波在测量面z=0处时域回波数据d′2(x,y,z=0,t);
如示意图5A所示,T为发射电磁波雷达装置,R为接收时域回波数据的雷达装置,两者均位于沿墙体外表面。至少接收时域回波数据的雷达装置沿与该墙体外表面平行的测量面内移动。其中,x、y和z为三维空间的三个坐标轴。x坐标轴的正方向为沿测量线向外的方向;y坐标轴的正方向为垂直于测量线向上的方向;z坐标轴的正方向为垂直于测量线向测量目标的方向。
步骤S204,对时域回波数据d′2(x,y,z=0,t)取三维傅里叶变换到频域得到频域回波数据:D2(kx,ky,z=0,ω);
其中,kx、ky和kz分别为空域x、y、z的波数。
步骤S206,对测量面z=0处的频域回波数据D2(kx,ky,z=0,ω)进行相位偏移,得到整个成像目标区域的频域回波数据D2(kx,ky,z,ω),其中,相位偏移的公式如下:
其中,kx,ky,kz,ω满足频散关系:(ω/c)2=k=kx 2+ky 2+kz 2,c为光速。
步骤S208,对整个成像目标区域的频域回波数据频谱数据D2(kx,ky,z,ω)乘以复数因子以补偿墙体对成像结果的影响,得到不同位置z处频域回波数据:
其中,
TTE,TTM分别是TE波和TM波的传输系数;
RTE,RTM别为TE波和TM波的反射系数;εr=ε1/ε0为分层介质两侧介电常数之比;
步骤S210,对进行墙体影响补偿后的频域回波数据D′2(kx,ky,z,ω),由Stolt变换 得到波数域回波数据
本步骤中,
步骤S212,按照以下公式对其关于kx,ky项求逆傅里叶变换得到空间波数域回波数据D2(x,y,z,kz):
本步骤中,由于kx,ky为均匀采样间隔,因此直接对关于kx,ky项求逆傅里叶变换即可得到数据D2(x,y,z,kz)。
步骤S214,对空间-波数域回波数据D2(x,y,z,kz),按照以下公式对其关于kz项采用NUFFT求逆变换到空时域d′2(x,y,z,t=0),从而得到消除墙体影响后的雷达成像结果:
本步骤中,由于kz为非均匀采样间隔,因此需要对数据D2(x,y,z,kz)关于kz项采用NUFFT求逆变换,才可以得到空时域数据d′2(x,y,z,t=0)并进行雷达成像。
为验证本实施例方法的效果,申请人进行了三维空间点目标对比试验。图6为三维空间,三层媒质点目标成像,采用现有技术方法未消除墙体影响,NUFFT处理,成像结果。图7为三维空间,三层媒质点目标成像,采用本实施例方法消除墙体影响后,NUFFT处理,成像结果。对比图6与图7可知,消除墙体影响后,三维空间三层媒质点目标成像结果比未消除墙体影响的成像结果更清晰,聚焦作用更明显。
本发明可用于二维空间两层至多层媒质穿墙雷达成像,三维空间两层至多层媒质穿墙雷达成像。
至此,已经结合附图对本发明两个实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明穿墙雷达成像的方法有了清楚的认识。
此外,上述对各元件和方法的定义并不仅限于实施方式中提到的各种具体结构、形状或方式,本领域的普通技术人员可对其进行简单地熟知地替换,例如:
(1)实施例1中,二维空间三层媒质还可以扩展到多层媒质形式,如两层至多层墙体。
(2)实施例2中三维空间三层媒质可以用多层媒质来代替。
综上所述,本发明提供一种新型穿墙雷达成像方法,以补偿墙体对成像位置和精度的影响,提高成像速度和时效,使穿墙成像结果准确、有效、实时。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种穿墙雷达成像的方法,其特征在于,包括:
步骤A,在墙体外发射雷达波,并接收该雷达波在z=0处的时域回波数据d′2(x,y,z=0,t);
步骤B,将所述时域回波数据取三维傅里叶变换,得到频域回波数据D2(kx,ky,z=0,ω);
步骤C,采用如下相位偏移公式,所述频域回波数据D2(kx,ky,z=0,ω)进行相位偏移,得到整个成像区域的频域回波数据D2(kx,ky,z,ω):
步骤D,对整个成像区域的频域回波数据D2(kx,ky,z,ω)乘以一个与墙体影响相关的复数因子以补偿墙体对成像结果的影响,得到频率频域回波数据:
其中,c为电磁波在自由空间的传播速度,||·||2表示取模,(·)*表示取共轭;
TTE,TTM分别是TE波和TM波的传输系数;
RTE,RTM分别为TE波和TM波的反射系数;εr=ε1/ε0为分层介质两侧介电常数之比;
步骤E,将乘以复数因子后的频域回波数据中的变量ω转换到变量kz,得到波数域回波数据,其中,ω为时间t对应的角频率,kz为空域变量z在波数域对应的波数变量;
步骤F,将所述波数域回波数据对其关于kx项求逆傅里叶变换得到空间-波数域回波数据,其中,kx为空域变量x在波数域对应的波数变量;以及
步骤G,对所述空间-波数域回波数据对其关于kz项变换到空时域,从而得到消除墙体影响后的雷达成像结果。
2.根据权利要求1所述的方法,其特征在于,应用于三维场景下,其中:
所述步骤E包括:对进行墙体影响补偿后的频域回波数据D′2(kx,ky,z,ω),由Stolt变换 得到波数域回波数据
所述步骤F包括:按照以下公式对波数域回波数据关于kx,ky项求逆傅里叶变换得到空间-波数域回波数据D2(x,y,z,kz):
所述步骤G包括:将空间-波数域回波数据D2(x,y,z,kz)变换到空时域d′2(x,y,z,t=0),从而得到消除墙体影响后的雷达成像结果。
3.根据权利要求2所述的方法,其特征在于,所述步骤G具体包括:
对数据D2(x,y,z,kz),按照以下公式对其关于kz项采用NUFFT求逆变换到空时域d′2(x,y,z,t=0),从而得到消除墙体影响后的雷达成像结果:
4.根据权利要求2所述的方法,其特征在于,所述步骤G具体包括:
对数据D2(x,y,z,kz),利用插值加FFT对其关于kz项求逆变换将其变换到空时域d′2(x,y,z,t=0),得到消除墙体影响后的雷达成像结果。
5.一种穿墙雷达成像的方法,其特征在于,包括:
步骤A,在墙体外发射雷达波,并接收该雷达波在z=0处的时域回波数据d′1(x,z=0,t);
步骤B,将所述时域回波数据变换到频域,得到频域回波数据D1(kx,z=0,ω);
步骤C,按照以下公式,将所述频域回波数据D1(kx,z=0,ω)进行相位偏移,得到整个成像区域的频域回波数据D1(kx,z,ω):
步骤D,对整个成像区域的频域回波数据D1(kx,z,ω)乘以一个与墙体影响相关的复数因子以补偿墙体影响,得到频域回波数据:
其中为分层媒质的传输系数,||·||2表示取模,d1为测线到第一层介质与第二层介质交界处的厚度,d2为测线到第二层介质与第三层介质交界处的厚度,为菲涅尔反射系数,k1为第一层介质对应的波数,k2为第二层介质对应的波数;
其中,所述第一层介质为空气;所述第二层介质为墙体;所述第三层介质为空气;
步骤E,将乘以复数因子后的频域回波数据中的变量ω转换到变量kz,得到波数域回波数据,其中,ω为时间t对应的角频率,kz为空域变量z在波数域对应的波数变量;
步骤F,将所述波数域回波数据对其关于kx项求逆傅里叶变换得到空间-波数域回波数据,其中,kx为空域变量x在波数域对应的波数变量;以及
步骤G,对所述空间-波数域回波数据对其关于kz项变换到空时域,从而得到消除墙体影响后的雷达成像结果。
6.根据权利要求5所述的方法,其特征在于,应用于二维场景下,其中:
所述步骤E中,对进行墙体影响补偿后的频域回波数据D′1(kx,z,ω),由Stolt变换由频域回波数据D′1(kx,z,ω)得到波数域数据
所述步骤F中,按照以下公式对其关于kx项求逆傅里叶变换得到空间-波数域回波数据D1(x,z,kz):
所述步骤G中,对空间-波数域数据D1(x,z,kz)对其求逆变换到空时域d′1(x,z,t=0),从而得到消除墙体影响后的雷达成像结果。
7.根据权利要求6所述的方法,其特征在于,所述步骤G具体包括:
对空间-波数域数据D1(x,z,kz)按照以下公式对其关于kz项采用NUFFT求逆变换到空时域d′1(x,z,t=0),从而得到消除墙体影响后的雷达成像结果:
8.根据权利要求6所述的方法,其特征在于,所述步骤G具体包括:利用插值加FFT对kz方向数据求逆变换得到穿墙雷达成像结果。
9.根据权利要求1至8中任一项所述的方法,其特征在于,所述步骤A中,发射电磁波的雷达装置和接收时域回波数据的雷达装置均位于沿墙体外表面,至少接收时域回波数据的雷达装置沿与该墙体外表面平行的z=0的测量线/测量面移动。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310572048.5A CN103675810B (zh) | 2013-11-13 | 2013-11-13 | 穿墙雷达成像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310572048.5A CN103675810B (zh) | 2013-11-13 | 2013-11-13 | 穿墙雷达成像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675810A CN103675810A (zh) | 2014-03-26 |
CN103675810B true CN103675810B (zh) | 2016-03-02 |
Family
ID=50313947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310572048.5A Active CN103675810B (zh) | 2013-11-13 | 2013-11-13 | 穿墙雷达成像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675810B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995256B (zh) * | 2014-05-29 | 2016-03-30 | 电子科技大学 | 一种适用于穿墙雷达成像的多面墙体快速补偿方法 |
CN108061890B (zh) * | 2016-11-08 | 2021-07-13 | 北京机电工程研究所 | 一种sar成像方法 |
CN108181625A (zh) * | 2018-01-30 | 2018-06-19 | 电子科技大学 | 一种基于非均匀频谱估计的钻孔雷达成像方法 |
CN108761446A (zh) * | 2018-04-09 | 2018-11-06 | 中国科学院电子学研究所 | 频率步进探地雷达的建模方法 |
CN108957450B (zh) * | 2018-07-10 | 2022-12-09 | 西安恒帆电子科技有限公司 | 一种毫米波雷达gpu实时三维成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858975A (zh) * | 2009-08-14 | 2010-10-13 | 电子科技大学 | 一种基于穿墙雷达成像的目标定位方法 |
CN101995573A (zh) * | 2009-08-12 | 2011-03-30 | 中国科学院电子学研究所 | 一种用于超宽带穿墙雷达的目标边界成像方法 |
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
CN102540184A (zh) * | 2011-12-30 | 2012-07-04 | 陕西延长石油(集团)有限责任公司研究院 | 探地雷达频域成像方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5594946B2 (ja) * | 2008-08-21 | 2014-09-24 | 三菱電機株式会社 | レーダ画像処理装置 |
US8207886B2 (en) * | 2008-09-11 | 2012-06-26 | Lawrence Livermore National Security, Llc | Model-based tomographic reconstruction |
KR100991348B1 (ko) * | 2008-11-28 | 2010-11-01 | 영남대학교 산학협력단 | Pso를 이용한 스텝주파수 레이더 영상의 요동 보상 방법 |
-
2013
- 2013-11-13 CN CN201310572048.5A patent/CN103675810B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101995573A (zh) * | 2009-08-12 | 2011-03-30 | 中国科学院电子学研究所 | 一种用于超宽带穿墙雷达的目标边界成像方法 |
CN101858975A (zh) * | 2009-08-14 | 2010-10-13 | 电子科技大学 | 一种基于穿墙雷达成像的目标定位方法 |
CN102520405A (zh) * | 2011-12-16 | 2012-06-27 | 中国人民解放军国防科学技术大学 | 基于图像域墙体影响补偿的穿墙成像方法 |
CN102540184A (zh) * | 2011-12-30 | 2012-07-04 | 陕西延长石油(集团)有限责任公司研究院 | 探地雷达频域成像方法 |
Non-Patent Citations (1)
Title |
---|
SAR穿墙成像中墙体影响图像域补偿方法;金添等;《电子学报》;20120731;第40卷(第7期);1423-1428 * |
Also Published As
Publication number | Publication date |
---|---|
CN103675810A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675810B (zh) | 穿墙雷达成像的方法 | |
CN107132520B (zh) | 一种基于水声超短基线定位系统的声线修正方法及系统 | |
CN104459633B (zh) | 结合局部频率估计的小波域InSAR干涉相位滤波方法 | |
CN103048642B (zh) | 基于频域最小二乘法的水声脉冲信号匹配场定位方法 | |
CN104360251B (zh) | 一种变压器局部放电的超声波信号时延估计方法 | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN102520405B (zh) | 基于图像域墙体影响补偿的穿墙成像方法 | |
CN102621549A (zh) | 多基线/多频段干涉相位解缠频域快速算法 | |
CN105022050A (zh) | 一种多传感器阵列的水声信道离散噪声源抑制方法 | |
CN103728591A (zh) | 一种mimo雷达近场目标高效实波束方位聚焦方法 | |
CN103018739A (zh) | 一种校正多通道幅相误差的微波三维成像方法 | |
CN104502911A (zh) | 一种穿墙成像雷达的墙壁参数估计方法 | |
CN104977565A (zh) | 一种基于分布式多阵列实现的水下目标三维定位方法 | |
CN103630903B (zh) | 基于顺轨干涉sar测量海面流场径向速度的方法 | |
CN103809180B (zh) | 用于InSAR地形测量的方位向预滤波处理方法 | |
CN103487796B (zh) | 一种利用水声信道统计不变特征实现被动测距的方法 | |
CN103487791B (zh) | 一种基于散射中心矩阵的rcs转换方法 | |
CN108681622B (zh) | 用于探地雷达波形优化的方法 | |
Flückiger et al. | Optimization of receiver arrangements for passive emitter localization methods | |
CN102680956A (zh) | 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统 | |
CN102932034B (zh) | 一种快速的宽带相干源方位估计方法 | |
CN104636551B (zh) | 一种网状反射面天线金属丝网等效电磁参数推演方法 | |
CN102183755A (zh) | 一种基于柯西-高斯模型的新型高分辨方位估计方法 | |
JP6211329B2 (ja) | 超音波厚み測定方法および超音波厚み測定システム | |
CN104090275A (zh) | 一种用于多探头阵列成像的幅度补偿方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |