本发明涉及雷达成像技术领域,特别涉及一种基于柱面雷达成像装置的稀疏天线布局方法。
背景技术:
为了保护公众安全,防范各种威胁,对公共场所特别是火车站、机场、客运站等来往人员安检必不可少。众多安检器中,微波成像是利用最为广泛。但是传统的安检器无法做到全方位探测,因而国内外研究机构把目光转向了柱面阵列雷达成像。
现有的柱面阵列雷达有两种工作方式,一种是利用收发天线在曲面内进行机械扫描;这种方法虽然阵元数量少,系统简单,但数据采集时间长,难以满足实际需求。另一种是利用实孔径或实孔径和合成孔径相结合的方式获取数据,但这种方法需耗费大量天线阵元且系统复杂度高。因此现有的柱面阵列雷达成像装置在满足奈奎斯特采样定理条件下,柱面阵列雷达阵元数量庞大,难以应用于实际系统。
技术实现要素:
本发明实施例的目的在于提供一种基于柱面雷达成像装置的稀疏天线阵元的布局方法,用于解决现有技术中柱面阵列雷达天线阵元数量庞大,难以应用于实际系统的问题。
为了解决上述技术问题,本发明的实施例采用了如下技术方案:一种基于柱面雷达成像装置的稀疏天线阵元的布局方法,包括如下步骤:
为了解决上述技术问题,本发明实施例提供一种对柱面雷达进行区域划分,获得若干阵元集合区域;
对各所述阵元集合区域进行区域划分获得若干阵元区域;
确定各所述阵元集合区域中位于首行和尾行的阵元区域为发射阵元区域;
确定各所述阵元集合区域中位于首列和尾列的、除了所述发射阵元区域之外的阵元区域为接收阵元区域;
在所述发射阵元区域处设置发射阵元;
在所述接收阵元区域处设置与所述发射阵元对应的接收阵元。
可选的,所述对柱面阵列雷达进行区域划分,获得若干阵元集合区域,具体包括:
确定柱面阵列雷达的高度;
按照预定的第一高度间距沿着所述柱面阵列雷达的高度方向对柱面阵列雷达的柱面进行区域划分,获得若干个第一区域;
按照预定的第一弧长沿第一区域的周长方向对所述第一区域进行区域划分,获得若干个阵元集合区域。
可选的,所述对各所述阵元集合区域进行区域划分获得若干阵元区域,具体包括:
确定各所述阵元集合区域的高度和弧长;
按照预定的第一份数确定第二高度间距;
按照预定的第二份数确定第二弧长间距;
按照所述第二高度间距沿着所述柱面阵列雷达的高度方向对阵元集合区域进行区域划分,获得若干第二区域;
按照所述第二弧长间距沿所述第二区域的弧长方向对各所述第二区域进行区域划分,获得若干阵元区域。
可选的,所述第二高度间距δh0满足:
其中,bwmax表示z轴方向最小波数支持的带宽;
fmax表示柱面阵列雷达装置的最大工作频率;
φ表示po2与xy平面的夹角;
θ表示po2与xy平面上的投影与x轴夹角;
p点表示目标点;
o2表示柱面阵列最上层阵元区域的圆平面圆心,其空间位置坐标为
可选的,所述第二弧长所述对应的柱面雷达的圆心角δθ0满足:
其中,
c表示光速表示等效相位中心到目标点的距离,是一个系统参数;
r表示等效相位中心到目标点p的距离,是一个系统参数;
p点的空间位置坐标表示为p(rp,θp,zp),其中rp表示p点到z轴的距离;zp表示p点到oxy平面的距离;
r表示目标点到原点的距离。
可选的,所述第一份数等于所述第二份数。
本发明实施例的有益效果在于:本发明中的基于柱面雷达成像装置的稀疏天线布局方法能够减少柱面阵列雷达装置阵元数量,减少数据采集时间,节约成本,提高成像效率。
附图说明
图1为本发明实施例基于柱面雷达成像装置的稀疏天线阵元的布局方法的流程图;
图2为本发明实施例中柱面阵列雷达中的阵元集合方阵示意图;
图3为本发明实施例中的柱面阵列雷达的示意图;
图4为本发明实施例中柱面阵列雷达的几何示意图;
图5为本发明实施例中部分阵元集合方阵中的阵元布局图;
图6为本发明实施例中阵元集合方阵的具体结构示意图;
图7为本发明实施例中多频点正交发射多通道信号接收的流程图;
图8为本发明实施例中三维成像流程图;
图9为本发明实施例中监测区域网的进行划分的立体示意图;
图10为图9的俯视图。
具体实施方式
此处参考附图描述本申请的各种方案以及特征。
应理解的是,可以对此处申请的实施例做出各种修改。因此,上述说明书不应该视为限制,而仅是作为实施例的范例。本领域的技术人员将想到在本申请的范围和精神内的其他修改。
包含在说明书中并构成说明书的一部分的附图示出了本申请的实施例,并且与上面给出的对本申请的大致描述以及下面给出的对实施例的详细描述一起用于解释本申请的原理。
通过下面参照附图对给定为非限制性实例的实施例的优选形式的描述,本申请的这些和其它特性将会变得显而易见。
还应当理解,尽管已经参照一些具体实例对本申请进行了描述,但本领域技术人员能够确定地实现本申请的很多其它等效形式,它们具有如权利要求所述的特征并因此都位于借此所限定的保护范围内。
当结合附图时,鉴于以下详细说明,本申请的上述和其他方面、特征和优势将变得更为显而易见。
此后参照附图描述本申请的具体实施例;然而,应当理解,所申请的实施例仅仅是本申请的实例,其可采用多种方式实施。熟知和/或重复的功能和结构并未详细描述以避免不必要或多余的细节使得本申请模糊不清。因此,本文所申请的具体的结构性和功能性细节并非意在限定,而是仅仅作为权利要求的基础和代表性基础用于教导本领域技术人员以实质上任意合适的详细结构多样地使用本申请。
本说明书可使用词组“在一种实施例中”、“在另一个实施例中”、“在又一实施例中”或“在其他实施例中”,其均可指代根据本申请的相同或不同实施例中的一个或多个。
本发明实施例提供一种基于柱面雷达成像装置的稀疏天线阵元布局方法,如图1所示,包括如下步骤:
步骤一,对柱面雷达进行区域划分,获得若干阵元集合区域;
步骤二,对各所述阵元集合区域进行区域划分获得若干阵元区域;
步骤三,确定各所述阵元集合区域中位于首行和尾行的阵元区域为发射阵元区域;
步骤四,确定各所述阵元集合区域中位于首列和尾列的、除了所述发射阵元区域之外的阵元区域为接收阵元区域;
步骤五,在所述发射阵元区域处设置发射阵元;
步骤六,在所述接收阵元区域处设置与所述发射阵元对应的接收阵元。
本发明实施例中,如图2所示通过对柱面雷达进行区域划分,获得若干阵元集合区域cell(即aij),然后再进一步对阵元集合区域cell进行划分获得若干阵元区域(阵元方格)。对于每个阵元集合区域cell,在首行和尾行的位置设置发射阵元,在首列和尾列的位置设置接收阵元,其余阵元区域不设置阵元,由此能够减少阵元数量,提高数据采集的时间。
本实施例提供一种基于柱面雷达成像装置的稀疏天线阵元布局方法,包括如下步骤:
步骤s1:阵元布局;根据柱面阵列雷达装置的大小,划分子阵,获取所有子阵所构成的子阵集方阵a,每个阵元集合方阵中布置阵元,获取每个阵元集合方阵的反射阵元矩阵。
步骤s1具体包括如下:
步骤s11:划分子阵;根据图2所示,将柱面阵列雷达划分成若干阵元集合区域
步骤s111:柱面阵列雷达高度划分;如图3所示,根据柱面阵列雷达的实际大小h0,将柱面阵列雷达在z轴上,分成na(na>1,na为正整数)份,则在y轴上的高度间距为h0,
步骤s112:柱面阵列雷达角度划分;以z轴为中心,xy平面等分成ma(ma>1,ma为正整数)份,每份的角度间距为θ0。θ0所对应的弧长即为第一弧长。
步骤s113:在步骤s111和步骤s112将柱面分成sub=na×ma份等大的阵元集合区域(子阵)后,每份阵元集合区域的高度为h0,水平角度间距为θ0,用所有阵元集合方阵构建一个na×ma的方阵,如图2所示,每一个小方格表示一个阵元集合区域用字母a表示,aij表示方阵中第i行第j列所对应的阵元集合区域(也可称为cell),
步骤s12:阵元布局;将步骤s11中的获取的每一个阵元集合区域cell,按照步骤s1中的划分方式将各阵元集合区域划分成ns*ms个阵元区域,每个阵元集合区域中的阵元区域之间的间距满足采样定理,相应的阵元区域中放置阵元,具体步骤如下:
步骤s121:阵元集合区域高度划分;根据阵元集合区域的实际高度h0,利用公式(1),将步骤s12中的每个阵元集合区域在z轴上,分成ns(ns>1,ns为正整数)份,则在z轴上的高度间距为δh0,
高度间距δh0,需满足采样定理:
其中
步骤s122:子阵角度划分;以z轴为中心,根据公式(4)将xy平面等分成ms(ms>1,ms为正整数)份,每份的角度间距为δθ0,角度间距δθ0所对应的弧长即为第二弧长间距。
角度间距δθ0需满足公式(5),
其中
本实施例中,较佳的ns和ms最佳选取为|ms-ns|=0,例如ms=100个,ns最优选择为ns=100。
步骤s123:获取阵元集合方阵;如图5所示,步骤s121和步骤s122将每个阵元集合区域分成ns×ms份等大的阵元区域,每个阵元区域的高度为δh0,水平角度间距为δθ0,将每个阵元集合区域构成的阵元集方阵用cell(即
步骤s124:阵元布局;如图6所示,在cell(即
步骤s125:获取发射阵元矩阵;将步骤s124中的发射阵元构建1行2ms列的行矩阵
步骤s126:获取接收阵元矩阵;将步骤s124中的接收阵元构建1行2(ns-2)列的行矩阵
采用本发明中的布局方法,布置的阵元个数为(2ms+2ns-4)*na*ma。而采用传统方式,柱面阵列雷达需要布置的阵元的个数为ns*ms*na*ma,由此可见,本发明实施例中的天线阵元的布局方法能够大大减少了阵元的布置数量,节约了成本,同时能够提高数据采集的效率。
本发明实施例中在完成天线阵元的布置之后,就可以利用柱面阵列雷达来进行信号的发射和接收,以获得若干回波信号,具体获取回波信号的步骤包括:
步骤一、根据柱面阵列雷达的最小工作频率和最大工作频率获取目标频点矩阵;
本步骤中,在具体实施过程中具体包括:根据最小工作频率、最大工作频率,利用第一计算公式计算获得频点差值;
根据所述最小工作频率以及所述频点差值利用第二计算公式计算获得所述目标频点矩阵中的频点元素,以获得目标频点矩阵。
步骤二、同时循环选取各阵元集合方阵中的一个发射阵元作为第一发射阵元;
本实施中,各所述阵元集合方阵中分别设置有若干发射阵元和若干接收阵元。
步骤三、为各所述第一发射阵元分别依次分配所述目标频点矩阵中的各频点,以作为第一发射阵元的发射频点;其中,同一时刻不同阵元集合方阵中的第一发射阵元的发射频点不同;
步骤四、利用第一发射阵元根据分配的发射频点发射信号;
步骤五、利用接收阵元接收所述发射信号,以获得回波数据。
本步骤中,接收所述发射信号的接收阵元与发射所述发射信号的第一发射阵元位于相同的阵元集合方阵。
具体的,根据柱面阵列雷达的最小工作频率和最大工作频率获取目标频点矩阵,具体包括:
根据最小工作频率、最大工作频率,利用第一计算公式计算获得频点差值;根据所述最小工作频率以及所述频点差值利用第二计算公式计算获得所述目标频点矩阵中的频点元素,以获得目标频点矩阵。具体的所述第一计算公式为:
本实施例中,所述第二计算公式为:
具体的,在获得所述频点元素后,还包括:利用所述频点元素构建获得第一频点矩阵;利用若干所述第一频点矩阵构建获得目标频点矩阵。
本实施例在上述实施例的基础上,作为一种具体的实施方式,获取回波信号具体过程为:
步骤s2:频点设置,获取频点矩阵;将工作频段设置成多个频点,获取频点矩阵,具体如下:
步骤s21:频点设置;根据柱面阵列雷达最小工作频率fmin(fmin≥1ghz)和最大工作频率fmax(fmax≤10thz)设置q(q≥n*m,q为正整数)个等差分布的频点,利用公式(8)(即第一计算公式),计算δf,将计算结果带入公式(9),求出
用
步骤s22:获取目标频点矩阵freq;利用步骤s21中f矩阵构建目标频点矩阵freq,freq矩阵大小为2q×1,
步骤s3:收发方式和策略;选取不同阵元集合方阵中的发射阵元,匹配阵元发射频点,发射、接收信号,流程图如图7所示;具体如下:
步骤s31:选取工作阵元;选取阵元集合方阵aij中trans的第kt个阵元
步骤s32:发射信号的频点匹配;根据阵元以及循环次数选取频点;工作的阵元所发射频点从目标频点矩阵freq中选择;
根据步骤s31中选取aij中
ifreq=(i-1)*ms+j+circulation-1(11)
其中i和j表示阵元集合方阵a对应的阵元集合方阵aij,circulation表示循环次数,circulation初始值为1,每次循环circulation=circulation+1,直至循环到circulation>q,也就是每个阵元循环发射每个频点,同一时刻不同阵元发射的频点不同,不同时刻同一阵元发射的频点不同,保证每个阵元都发射所有不同大小的频点,而且仅发射一次;
步骤s33:阵元发射、接收信号;具体如下:
步骤s331:阵元发射信号;根据步骤s31选取aij中
st(a,trans,freq)=σ(x,y,z)e-j2πft(12)
其中a表示阵元集合方阵,trans表示发射阵元,freq表示目标频点矩阵,σ(x,y,z)表示p目标像素点的后向散射特性,f表示选取的阵元集合方阵a中阵元trank发射的频点,选取方式,t表示信号传播的时间;
步骤s332:接收阵元接收信号以获得回波信号。接收阵元只接收本阵元集合方阵中的发射阵元发射的发射信号,阵元接收的信号表达式为公式(13),
s(a,trans,rece,freq)=∫∫∫vσ(x,y,z)e-j2πfτdxdydz(13)
v表示柱面阵列雷达观测区域,τ表示时间延迟,其表达式为公式(14),
图4中左边小方格阴影区表示一个阵元集合方阵,阵元集合方阵中含有发射阵元和接收阵元,发射阵元的空间位置坐标为trans(r0,θt,zt)或表示成trans(xt,yt,zt),接收阵元的空间位置坐标为rece(r0,θr,zr)或表示成rece(xr,yr,zr),p点(即目标像素点)的空间位置坐标为p(r,θ,z)或表示成p(x,y,z),rt表示发射阵元与p点的距离,rr表示接收阵元与p的距离,如下所示;
接收阵元接收的信号可以写成公式(17):
本实施例中通过采用上述的信号收方式来获得回波信号,能够减少数据的采集时间,为后续快速成像提供了基础。
本发明实施例在上述实施例的基础上,在获得了回波信号之后就可以进行三维成像,具体成像的过程为:
步骤s4:三维成像;根据接收阵元接收的回波信号,应用共焦投影、bp等三维成像算法实现三维成像,步骤流程图如图8所示,具体如下:
步骤s41:监测区域网格化分,具体如下:
步骤s411:监测区域网格化分;计算监测区域的空间坐标和点数;把监测区域等角度、等距离、等高度划分网格。如图9和图10所示,图9中表示划分网格三维示意图,图10表示二维划分网格二维示意图(即图9的俯视图),黑点表示三维像素点位置,网格坐标用i(u,v,w)表示,u表示x轴坐标,v表示y轴坐标,w表示z轴坐标,每一个网格称为一个三维像素点,距离方向向等分成mu,角度方向上等分成mv,高度方向上等分成份mw,共有m=mumvmw个像素点,第m(0≤m≤m)个像素点im(即目标像素点)的坐标表示为(um,vm,wm);
步骤s412:三维像素点散射强度初始化;令所有三维像素点的散射强度σsum(im)(σsum(im)表示像素点im的散射强度)为0,即σsum(im)=0;
步骤s42:计算匹配滤波函数(即每个接收阵元接收的回波信号对应一个滤波函数);循环选取的阵元集合方阵
步骤s421:系统初始化;对像素点、阵元集合方阵
步骤s4211:像素点初始化;选取第一个像素点,即令m=1,得到im=i1;
步骤s4212:阵元集合方阵初始化;选取第一个阵元集合方阵,即令ka=1,得到
步骤s4213:发射阵元初始化;选取第一个发射阵元,即令kt=1,得到
步骤s4214:接收阵元初始化;选取第一个接收阵元,即令kr=1,得到trans,系统初始化结束,执行步骤s423;
步骤s422:循环选取;主要有像素点i循环、阵元集合方阵a循环、发射阵元trans循环以及接收阵元rece循环,具体如下:
步骤s4221:接收阵元rece循环;若kr≤2(ns-2),需将接收阵元循环到下一个(kr=kr+1),执行步骤s423;反之则执行步骤s4222;
步骤s4222:发射阵元trans循环;若kt≤2ms,需将发射阵元循环到下一个(kt=kt+1),且对接收阵元进行初始化操作,最后执行步骤s423;反之则执行步骤s4223;
步骤s4223:阵元集合方阵a循环;若ka≤sub,需将阵元集合方阵循环到下一个(ka=ka+1),且对发射阵元trans和接收阵元rece进行初始化操作,再执行步骤s423;反之则执行步骤s4224;
步骤s4224:保存像素点im的散射强度;
步骤s4225:像素点i循环;若m≤m,首先需将像素点循环到下一个(im=im+1),然后选取σsum中第im元素,对阵元集合方阵a、发射阵元trans和接收阵元rece进行初始化操作,最后执行步骤s423;反之则表示已经计算完成所有像素点im的散射强度σsum(im),执行步骤s48实现三维成像;
步骤s423:计算匹配滤波函数;根据选取的阵元集合方阵、发射阵元、接收阵元以及像素点计算匹配滤波函数,所述滤波函数的表达式为:
其中:
步骤s43:将匹配滤波函数
步骤s44:匹配信号逆傅里叶变换;将匹配信号
步骤s45:取距离向压缩信号的峰值;距离向压缩信号
|first表示距离向压缩信号
步骤s46:距离向压缩信号峰值处理;根据方位向合成孔径和高程向合成孔径范围,对距离向压缩信号的峰值
步骤s461:计算方位向筛选条件;方位向筛选条件如下所示,
angle(um+vmi)函数用来求解复数um+vmi在不同象限中的相角;θa表示合成孔径方位向范围,是系统参数;其中,θt表示发射阵元的空间位置坐标中的方位角;θr表示接收阵元的空间位置坐标中的方位角;um表示第m个三维像素点的空间位置坐标中的x轴的坐标;vm表示第m个三维像素点的空间位置坐标中的y轴的坐标。
步骤s462:计算高程向筛选条件;高程向筛选条件如下所示,
θh表示合成孔径高程向范围,是系统参数;um表示第m个三维像素点的空间位置坐标中的x轴的坐标;vm表示第m个三维像素点的空间位置坐标中的y轴的坐标;wm表示第m个三维像素点的空间位置坐标中的z轴的坐标;r0表示柱面阵列的底面半径;zt表示发射阵元的空间位置坐标中的z轴的坐标;zr表示接收阵元的空间位置坐标中的z轴的坐标;
步骤s463:距离向压缩信号的峰值处理;根据方位向筛选条件和高程向筛选条件,令不符合条件的距离向压缩信号的峰值
true表示满足方位向筛选条件、高程向筛选条件的距离向压缩信号的峰值;false表示只满足其中一个筛选条件,或都不满足筛选条件;
步骤s47:相干叠加;根据像素点im、方位向筛选条件、距离向筛选条件,将距离向压缩信号的峰值进行相干叠加σsum(im),计算像素点im的散射强度:
sub表示阵元集合方阵的个数,ka表示第ka个阵元集合方阵;2ms表示一个阵元集合方阵中发射阵元的个数,kt表示ka阵元集合方阵中第kt个发射阵元;2(ns-2)表示一个阵元集合方阵中接收阵元的个数,kr表示ka阵元集合方阵中第kr个接收阵元;
步骤s48:三维复图像;根据像素点im的坐标和计算像素点im的对应的散射强度σsum(im)(步骤s4224),构建tar_sca=[umvmwmσsum(im)]矩阵,其中第一列表示像素点im的x横坐标,第二列表示像素点im的y坐标,第三列表示像素点im的z坐标,第四列表示像素点im的散射强度,利用tar_sca矩阵绘出柱面阵列雷达三维复图像。
本实施例中的柱面阵列雷达成像方法,能够实现360°全方位探测,并且所采用的阵元相对较少,降低成本,还能减少数据量,提高了成像效率。
以上实施例仅为本发明的示例性实施例,不用于限制本发明,本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内,对本发明做出各种修改或等同替换,这种修改或等同替换也应视为落在本发明的保护范围内。