CN109856634A - 一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 - Google Patents
一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 Download PDFInfo
- Publication number
- CN109856634A CN109856634A CN201910166018.1A CN201910166018A CN109856634A CN 109856634 A CN109856634 A CN 109856634A CN 201910166018 A CN201910166018 A CN 201910166018A CN 109856634 A CN109856634 A CN 109856634A
- Authority
- CN
- China
- Prior art keywords
- polar coordinates
- image
- imaging
- unit
- graing lobe
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法,包括以下步骤:S1、设定极坐标成像距离单元和方位角角度单元,并根据角度和距离的划分构建极坐标成像网格矩阵,S2、通过目标散射回波进行极坐标后向投影成像;S3、对极坐标成像图像进行单元均值恒虚警检测,得到二值图像;S4、用二值图像对成像图像进行掩膜,得到掩膜图像;S5、提取掩膜图像像素连通区域的质心并按照质心的距离值对质心进行聚类;S6、比较各个聚类中连通区域质心的像素值大小。本发明利用极坐标后投影成像中主瓣和栅瓣位于同一距离的特点,并通过对极坐标成像图像进行单元平均恒虚警检测和方位向投影聚类,实现了主瓣和栅瓣的有效分离,具有计算简单灵活和效率高的特点。
Description
技术领域
本发明涉及到阵列成像技术领域,尤其涉及一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法。
背景技术
超宽带成像雷达具有良好的穿透能力、优良的分辨性能,在灾害生命探测救援、反恐维稳穿墙探测以及可疑物品检测等领域有着广泛的应用。然而,为了满足这些应用领域对雷达便携性的要求,超宽带成像雷达往往采用孔径大小有限且阵元数量有限的稀疏阵列,成像质量受限于稀疏阵列引入的方位栅瓣的干扰。
常见的成像方法均是在笛卡尔坐标系下,栅瓣和主瓣的位置位于椭圆弧线上(如图2所示),难以利用栅瓣和主瓣的位置关系来区分栅瓣和主瓣;基于CF加权的栅瓣抑制方法,仅仅是将栅瓣抑制到一定水平,并不能真正去除栅瓣;若后续的目标检测处理不当,依旧无法排除栅瓣的干扰。
鉴于此,研究一种实现主瓣和栅瓣分离的超宽带成像的方法是本技术领域人员亟待解决的技术问题。
发明内容
本发明的目的是提供一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法,该方法利用极坐标成像中的栅瓣能量低于主瓣能量的特点,对成像图像进行单元平均恒虚警检测和掩膜聚类,实现了栅瓣的去除,计算灵活,有效避免了栅瓣对成像质量的干扰,从而提高成像图像的质量,有利于超宽带成像雷达对目标的探测,抑制虚假目标。
为解决上述技术问题,本发明提供一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法,所述方法包括以下步骤:
S1、在极坐标中,设定成像距离单元的距离范围[Rmin,Rmax]和距离间隔dR,以及成像方位角角度单元的角度范围[θmin,θmax]和角度间隔dθ,并根据极坐标成像距离单元和成像方位角角度单元的划分构建成像网络矩阵G;
S2、基于K个收发通道的目标散射回波,对所述步骤S1中的网络矩阵G上的极坐标进行后向投影成像,得到极坐标成像图像P;
S3、对所述步骤S2中极坐标成像图像P进行单元平均恒虚警检测,得到检测结果的二值图像W;
S4、通过所述步骤S3中的二值图像W对所述步骤S2中极坐标成像图像P进行掩膜操作,得到极坐标成像图像P的掩膜图像Q;
S5、提取所述步骤S4中掩膜图像Q像素连通区域的质心,并通过K-means聚类方法对所述掩膜图像Q连通区域的质心坐标中的距离值进行聚类,得到多个连通区域质心距离值的聚类;
S6、分别对所述步骤S5中的每一个聚类中的连通区域像素值大小进行比较,选择像素值最大的连通区域作为主瓣区域,余下的连通区域作为栅瓣区域予以删除。
优选地,所述步骤S1中成像方位角角度单元的角度范围满足条件-90°≤θmin<θmax≤90°。
优选地,所述步骤S2中的散射回波为雷达回波。
优选地,所述步骤S2中极坐标成像图像P的计算公式可表示为:
式(1)中,m表示极坐标成像图像P中方位角角度单元的像素点坐标索引,n表示极坐标成像图像P中距离单元的像素点坐标索引,s(k,t)表示雷达回波,t表示雷达回波的传播时间,(Txk,Tyk)表示第k个通道中发射天线的坐标,(Rxk,Ryk)表示第k个通道中接收天线的坐标,θm表示网格矩阵G中第m个角度单元的角度值,Rn表示网格矩阵G中第n个距离单元的距离值,j表示虚数单位。
优选地,所述步骤S4中对极坐标成像图像P进行掩膜操作,掩膜操作公式可表示为:
Q(m,n)=W(m,n)×P(m,n) (2)
式(2)中,W(m,n)表示点(m,n)的极坐标二值图像像素,P(m,n)表示点(m,n)的极坐标成像图像像素。
与现有技术比较,本发明的超宽带稀疏阵列极坐标成像下的栅瓣去除方法,利用极坐标后投影成像中主瓣和栅瓣位于同一距离,且栅瓣能量低于主瓣能量的特点,并通过对极坐标成像图像进行单元平均恒虚警检测和方位向投影聚类,实现了主瓣和栅瓣的有效分离,具有计算简单灵活和效率高的特点。
附图说明
图1是本发明一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法流程图,
图2是本发明仿真实验中笛卡尔坐标后向投影成像结果示意图,
图3是本发明仿真实验中极坐标后向投影成像结果,
图4为本发明仿真实验中极坐标后向投影成像检测结果示意图,
图5为本发明仿真实验中笛卡尔坐标后向投影成像掩膜结果示意图,
图6为本发明仿真实验中的栅瓣去除结果示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明的技术方案,下面结合附图对本发明作进一步的详细说明。
如图1所示,一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法,所述方法包括以下步骤:
S1、在极坐标中,设定成像距离单元的距离范围[Rmin,Rmax]和距离间隔dR,以及成像方位角角度单元的角度范围[θmin,θmax]和角度间隔dθ,并根据极坐标成像距离单元和成像方位角角度单元的划分构建成像网络矩阵G;
S2、基于K个收发通道的目标散射回波,对所述步骤S1中的网络矩阵G上的极坐标进行后向投影成像,得到极坐标成像图像P;
S3、对所述步骤S2中极坐标成像图像P进行单元平均恒虚警检测,得到检测结果的二值图像W;
S4、通过所述步骤S3中的二值图像W对所述步骤S2中极坐标成像图像P进行掩膜操作,得到极坐标成像图像P的掩膜图像Q;
S5、提取所述步骤S4中掩膜图像Q像素连通区域的质心,并通过K-means聚类方法对所述掩膜图像Q连通区域的质心坐标中的距离值进行聚类,得到多个连通区域质心距离值的聚类;
S6、分别对所述步骤S5中的每一个聚类中的连通区域像素值大小进行比较,选择像素值最大的连通区域作为主瓣区域,余下的连通区域作为栅瓣区域予以删除。
本实施例中,利用极坐标后投影成像中主瓣和栅瓣位于同一距离,且栅瓣能量低于主瓣能量的特点,并通过对极坐标成像图像P进行单元平均恒虚警检测和方位向投影聚类,实现了主瓣与栅瓣的有效分离,排除了栅瓣对成像质量的影响,具有计算简单灵活和效率高的特点。
如图1所示,所述步骤S1中成像方位角角度单元的角度范围满足条件-90°≤θmin<θmax≤90°。
如图1所示,所述步骤S2中的散射回波为雷达回波。
如图1所示,所述步骤S2中极坐标成像图像P的计算公式可表示为:
式(1)中,m表示极坐标成像图像P中方位角角度单元的像素点坐标索引,n表示极坐标成像图像P中距离单元的像素点坐标索引,s(k,t)表示雷达回波,t表示雷达回波的传播时间,(Txk,Tyk)表示第k个通道中发射天线的坐标,(Rxk,Ryk)表示第k个通道中接收天线的坐标,θm表示网格矩阵G中第m个角度单元的角度值,Rn表示网格矩阵G中第n个距离单元的距离值,j表示虚数单位。
如图1所示,所述步骤S4中对极坐标成像图像P进行掩膜操作,掩膜操作公式可表示为:
Q(m,n)=W(m,n)×P(m,n) (2)
式(2)中,W(m,n)表示点(m,n)的极坐标二值图像像素,P(m,n)表示点(m,n)的极坐标成像图像像素。
为了更好的理解本发明的工作原理和技术效果,下面结合一个仿真实验予以说明。
首先,设定超宽带稀疏阵的参数,具体如表1所示:
表1超宽带稀疏阵列的参数
超宽带信号起始频率 | 超宽带信号带宽 | 稀疏阵列长度 | 稀疏阵列阵元间距 |
1GHz | 400MHz | 2m | 0.5m |
S1、在极坐标中,设定成像距离单元的距离范围为[0,10],其距离间隔dR=0.01,以及成像方位角角度单元的角度范围为[-60°,60°],其角度间隔dθ=0.02°,并根据极坐标成像距离单元和成像方位角角度单元的划分构建成像网络矩阵G;
S2、基于K(K=8)个收发通道的目标散射回波,对所述网络矩阵G上的极坐标进行后向投影成像,得到极坐标成像图像P,从图3中可以看出,极坐标后向投影成像中同样存在栅瓣,且栅瓣的能量比目标主瓣的能量低。但与图2中笛卡尔坐标后向投影成像结果不同的是,极坐标后向投影成像中目标和栅瓣均位于同一距离上;
S3、对所述极坐标成像图像P进行单元平均恒虚警检测,得到检测结果的二值图像W,从图4中可以看出,栅瓣极易被认定为目标,引起虚警;
S4、通过所述二值图像W对所述极坐标成像图像P进行掩膜操作,得到极坐标成像图像P的掩膜图像Q;
S5、提取所述掩膜图像Q像素连通区域的质心,并通过K-means聚类方法对所述掩膜图像Q连通区域的质心坐标中的距离值进行聚类,得到两个连通区域质心距离值的聚类,如图5所示,图中的星号标记分别表示掩膜图像中六个连通区域的质心,即可将六个连通区域聚为两类;
S6、分别对所述两个聚类中的连通区域像素值大小进行比较,选择像素值最大的连通区域作为主瓣区域,余下的连通区域作为栅瓣区域予以删除,从图6中可以看出,通过对极坐标后投影成像进行恒虚警检测和掩膜聚类,能够有效去除栅瓣,实现主瓣与栅瓣的有效分离,提升成像的质量。
以上对本发明所提供的一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (5)
1.一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法,其特征在于,所述方法包括以下步骤:
S1、在极坐标中,设定成像距离单元的距离范围[Rmin,Rmax]和距离间隔dR,以及成像方位角角度单元的角度范围[θmin,θmax]和角度间隔dθ,并根据成像距离单元和成像方位角角度单元的划分构建成像网络矩阵G;
S2、基于K个收发通道的目标散射回波,对所述步骤S1中的网络矩阵G上的极坐标进行后向投影成像,得到极坐标成像图像P;
S3、对所述步骤S2中极坐标成像图像P进行单元平均恒虚警检测,得到检测结果的二值图像W;
S4、通过所述步骤S3中的二值图像W对所述步骤S2中极坐标成像图像P进行掩膜操作,得到极坐标成像图像P的掩膜图像Q;
S5、提取所述步骤S4中掩膜图像Q像素连通区域的质心,并通过K-means聚类方法对所述掩膜图像Q连通区域的质心坐标中的距离值进行聚类,得到多个连通区域质心距离值的聚类;
S6、分别对所述步骤S5中的每一个聚类中的连通区域像素值大小进行比较,选择像素值最大的连通区域作为主瓣区域,余下的连通区域作为栅瓣区域予以删除。
2.如权利要求1所述的超宽带稀疏阵列极坐标成像下的栅瓣去除方法,其特征在于,所述步骤S1中成像方位角角度单元的角度范围满足条件-90°≤θmin<θmax≤90°。
3.如权利要求2所述的超宽带稀疏阵列极坐标成像下的栅瓣去除方法,其特征在于,所述步骤S2中的散射回波为雷达回波。
4.如权利要求3所述的超宽带稀疏阵列极坐标成像下的栅瓣去除方法,其特征在于,所述步骤S2中极坐标成像图像P的计算公式可表示为:
式(1)中,m表示极坐标成像图像P中方位角角度单元的像素点坐标索引,n表示极坐标成像图像P中距离单元的像素点坐标索引,s(k,t)表示雷达回波,t表示雷达回波的传播时间,(Txk,Tyk)表示第k个通道中发射天线的坐标,(Rxk,Ryk)表示第k个通道中接收天线的坐标,θm表示网格矩阵G中第m个角度单元的角度值,Rn表示网格矩阵G中第n个距离单元的距离值,j表示虚数单位。
5.如权利要求4所述的超宽带稀疏阵列极坐标成像下的栅瓣去除方法,其特征在于,所述步骤S4中对极坐标成像图像P进行掩膜操作,掩膜操作公式可表示为:
Q(m,n)=W(m,n)×P(m,n) (2)
式(2)中,W(m,n)表示点(m,n)的极坐标二值图像像素,P(m,n)表示点(m,n)的极坐标成像图像像素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910166018.1A CN109856634B (zh) | 2019-03-06 | 2019-03-06 | 一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910166018.1A CN109856634B (zh) | 2019-03-06 | 2019-03-06 | 一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109856634A true CN109856634A (zh) | 2019-06-07 |
CN109856634B CN109856634B (zh) | 2022-12-13 |
Family
ID=66899865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910166018.1A Active CN109856634B (zh) | 2019-03-06 | 2019-03-06 | 一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109856634B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112014835A (zh) * | 2020-09-01 | 2020-12-01 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4649390A (en) * | 1983-08-05 | 1987-03-10 | Hughes Aircraft Company | Two dimension radar system with selectable three dimension target data extraction |
US20050200541A1 (en) * | 2004-03-09 | 2005-09-15 | The Boeing Company | System and method for preferentially controlling grating lobes of direct radiating arrays |
CN101566689A (zh) * | 2009-06-10 | 2009-10-28 | 电子科技大学 | 一种消除三维合成孔径雷达图像栅瓣的方法 |
CN101813764A (zh) * | 2010-03-25 | 2010-08-25 | 电子科技大学 | 一种针对均匀超稀疏阵列天线波束指向模糊的抑制方法 |
CN104391276A (zh) * | 2014-10-08 | 2015-03-04 | 西安电子工程研究所 | 一种收发分置平面阵相控阵雷达天线布阵及波束形成器 |
CN104808200A (zh) * | 2015-01-27 | 2015-07-29 | 湖南华诺星空电子技术有限公司 | 一种超宽带雷达成像方位向栅瓣的抑制方法 |
CN105334508A (zh) * | 2015-09-24 | 2016-02-17 | 哈尔滨工程大学 | 一种稀疏阵列宽带波束形成的栅瓣抑制方法 |
US20170164835A1 (en) * | 2014-06-10 | 2017-06-15 | Ithera Medical Gmbh | Device and method for hybrid optoacoustic tomography and ultrasonography |
US20180151961A1 (en) * | 2016-11-28 | 2018-05-31 | Mando Corporation | Radar apparatus and antenna apparatus therefor |
CN108594233A (zh) * | 2018-04-24 | 2018-09-28 | 森思泰克河北科技有限公司 | 一种基于mimo汽车雷达的速度解模糊方法 |
CN108872970A (zh) * | 2018-06-14 | 2018-11-23 | 苏州桑泰海洋仪器研发有限责任公司 | 适用于一般等间距稀疏阵列单频信号波束形成的栅瓣判别方法 |
JP2019030623A (ja) * | 2016-11-15 | 2019-02-28 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
-
2019
- 2019-03-06 CN CN201910166018.1A patent/CN109856634B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4649390A (en) * | 1983-08-05 | 1987-03-10 | Hughes Aircraft Company | Two dimension radar system with selectable three dimension target data extraction |
US20050200541A1 (en) * | 2004-03-09 | 2005-09-15 | The Boeing Company | System and method for preferentially controlling grating lobes of direct radiating arrays |
CN101566689A (zh) * | 2009-06-10 | 2009-10-28 | 电子科技大学 | 一种消除三维合成孔径雷达图像栅瓣的方法 |
CN101813764A (zh) * | 2010-03-25 | 2010-08-25 | 电子科技大学 | 一种针对均匀超稀疏阵列天线波束指向模糊的抑制方法 |
US20170164835A1 (en) * | 2014-06-10 | 2017-06-15 | Ithera Medical Gmbh | Device and method for hybrid optoacoustic tomography and ultrasonography |
CN104391276A (zh) * | 2014-10-08 | 2015-03-04 | 西安电子工程研究所 | 一种收发分置平面阵相控阵雷达天线布阵及波束形成器 |
CN104808200A (zh) * | 2015-01-27 | 2015-07-29 | 湖南华诺星空电子技术有限公司 | 一种超宽带雷达成像方位向栅瓣的抑制方法 |
CN105334508A (zh) * | 2015-09-24 | 2016-02-17 | 哈尔滨工程大学 | 一种稀疏阵列宽带波束形成的栅瓣抑制方法 |
JP2019030623A (ja) * | 2016-11-15 | 2019-02-28 | 炭 親良 | ビームフォーミング方法、計測イメージング装置、及び、通信装置 |
US20180151961A1 (en) * | 2016-11-28 | 2018-05-31 | Mando Corporation | Radar apparatus and antenna apparatus therefor |
CN108594233A (zh) * | 2018-04-24 | 2018-09-28 | 森思泰克河北科技有限公司 | 一种基于mimo汽车雷达的速度解模糊方法 |
CN108872970A (zh) * | 2018-06-14 | 2018-11-23 | 苏州桑泰海洋仪器研发有限责任公司 | 适用于一般等间距稀疏阵列单频信号波束形成的栅瓣判别方法 |
Non-Patent Citations (3)
Title |
---|
JUN HU等: "Grating_lobe_mitigation_based_on_extended_coherence_factor_in_sparse_MIMO_UWB_array", 《IEEE》 * |
LONG ZHUANG等: "Coherent_synthesis_sparse_aperture_radar_with_grating_lobes_suppressed_using_frequency_MIMO_technique", 《IEEE》 * |
涂新等: "穿墙成像雷达旁瓣/栅瓣抑制方法", 《雷达科学与技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112014835A (zh) * | 2020-09-01 | 2020-12-01 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
CN112014835B (zh) * | 2020-09-01 | 2023-05-26 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109856634B (zh) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104020451B (zh) | 基于聚类的外辐射源雷达目标航迹处理方法 | |
US9002286B2 (en) | Method and system for identification and mitigation of errors in non-line-of-sight distance estimation | |
DE112013002115B4 (de) | Orientierung eines Ultraschallsignals | |
CN103197302B (zh) | 一种适用于穿墙雷达成像的目标位置提取方法 | |
CN106597368B (zh) | 室内干扰源的定位方法及系统 | |
CN105044684B (zh) | 基于射频隐身的mimo跟踪雷达发射波束的形成方法 | |
CN103630910A (zh) | 一种gnss接收机设备的抗干扰方法 | |
DE60109990T2 (de) | Monopuls-Radar-Prozessor zur Bestimmung des Azimuth- und Elevationswinkels zweier Ziele über eine Matrix von Amplitudenverhältnissen | |
CN105891815A (zh) | 一种基于广播信号无源定位的联合估计算法 | |
CN111812630B (zh) | 干扰剩余时外辐射源雷达目标检测与doa估计系统及方法 | |
DE112017006884T5 (de) | Radarvorrichtung | |
CN107479043A (zh) | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 | |
WO2006122316A2 (en) | Method for detecting navigation beacon signals using two antennas or equivalent thereof | |
EP3398262A1 (en) | Micro-route characterization and selection | |
De Wilde et al. | Authentication by polarization: A powerful anti-spoofing method | |
CN111580099A (zh) | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 | |
CN109856634A (zh) | 一种超宽带稀疏阵列极坐标成像下的栅瓣去除方法 | |
CN108872977A (zh) | 基于单通道超宽带雷达的生命体双站协同探测方法 | |
CN109143198A (zh) | 副瓣抑制方法及系统 | |
CN112327292B (zh) | 一种二维稀疏阵列doa估计方法 | |
Barott et al. | Cochannel interference in ATSC passive radar | |
Jagadesh et al. | Modeling Target Detection and Performance Analysis of Electronic Countermeasures for Phased Radar. | |
CN115575921B (zh) | 一种基于俯仰向多通道多干扰基压制干扰抑制方法 | |
WO2023109755A1 (zh) | 感知方法、装置及通信设备 | |
CN114111456B (zh) | 一种电力变电站及重要设施无人机防御系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |