CN110837127A - 一种基于柱面雷达成像装置的稀疏天线布局方法 - Google Patents

一种基于柱面雷达成像装置的稀疏天线布局方法 Download PDF

Info

Publication number
CN110837127A
CN110837127A CN201911174015.9A CN201911174015A CN110837127A CN 110837127 A CN110837127 A CN 110837127A CN 201911174015 A CN201911174015 A CN 201911174015A CN 110837127 A CN110837127 A CN 110837127A
Authority
CN
China
Prior art keywords
array element
array
radar
transmitting
region
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
Application number
CN201911174015.9A
Other languages
English (en)
Other versions
CN110837127B (zh
Inventor
谭维贤
胡礼珍
黄平平
徐伟
洪文
李亚超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inner Mongolia University of Technology
Original Assignee
Inner Mongolia University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Inner Mongolia University of Technology filed Critical Inner Mongolia University of Technology
Priority to CN201911174015.9A priority Critical patent/CN110837127B/zh
Publication of CN110837127A publication Critical patent/CN110837127A/zh
Application granted granted Critical
Publication of CN110837127B publication Critical patent/CN110837127B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于柱面雷达成像装置的稀疏天线阵元布局方法。包括：对柱面雷达进行区域划分，获得若干阵元集合区域；对各所述阵元集合区域进行区域划分获得若干阵元区域；确定各所述阵元集合区域中位于首行和尾行的阵元区域为发射阵元区域；确定各所述阵元集合区域中位于首列和尾列的、除了所述发射阵元区域之外的阵元区域为接收阵元区域；在所述发射阵元区域处设置发射阵元；在所述接收阵元区域处设置与所述发射阵元对应的接收阵元；本发明的基于柱面雷达成像装置的稀疏天线阵元布局方法，能够减少柱面阵列雷达中的阵元的数量，提高数据采集时间，进而提高成像效率。

Description

一种基于柱面雷达成像装置的稀疏天线布局方法
技术领域
本发明涉及雷达成像技术领域，特别涉及一种基于柱面雷达成像装置的稀疏天线布局方法。
背景技术
为了保护公众安全，防范各种威胁，对公共场所特别是火车站、机场、客运站等来往人员安检必不可少。众多安检器中，微波成像是利用最为广泛。但是传统的安检器无法做到全方位探测，因而国内外研究机构把目光转向了柱面阵列雷达成像。
现有的柱面阵列雷达有两种工作方式，一种是利用收发天线在曲面内进行机械扫描；这种方法虽然阵元数量少，系统简单，但数据采集时间长，难以满足实际需求。另一种是利用实孔径或实孔径和合成孔径相结合的方式获取数据，但这种方法需耗费大量天线阵元且系统复杂度高。因此现有的柱面阵列雷达成像装置在满足奈奎斯特采样定理条件下，柱面阵列雷达阵元数量庞大，难以应用于实际系统。
发明内容
本发明实施例的目的在于提供一种基于柱面雷达成像装置的稀疏天线阵元的布局方法，用于解决现有技术中柱面阵列雷达天线阵元数量庞大，难以应用于实际系统的问题。
为了解决上述技术问题，本发明的实施例采用了如下技术方案：一种基于柱面雷达成像装置的稀疏天线阵元的布局方法，包括如下步骤：
为了解决上述技术问题，本发明实施例提供一种对柱面雷达进行区域划分，获得若干阵元集合区域；
对各所述阵元集合区域进行区域划分获得若干阵元区域；
确定各所述阵元集合区域中位于首行和尾行的阵元区域为发射阵元区域；
确定各所述阵元集合区域中位于首列和尾列的、除了所述发射阵元区域之外的阵元区域为接收阵元区域；
在所述发射阵元区域处设置发射阵元；
在所述接收阵元区域处设置与所述发射阵元对应的接收阵元。
可选的，所述对柱面阵列雷达进行区域划分，获得若干阵元集合区域，具体包括：
确定柱面阵列雷达的高度；
按照预定的第一高度间距沿着所述柱面阵列雷达的高度方向对柱面阵列雷达的柱面进行区域划分，获得若干个第一区域；
按照预定的第一弧长沿第一区域的周长方向对所述第一区域进行区域划分，获得若干个阵元集合区域。
可选的，所述对各所述阵元集合区域进行区域划分获得若干阵元区域，具体包括：
确定各所述阵元集合区域的高度和弧长；
按照预定的第一份数确定第二高度间距；
按照预定的第二份数确定第二弧长间距；
按照所述第二高度间距沿着所述柱面阵列雷达的高度方向对阵元集合区域进行区域划分，获得若干第二区域；
按照所述第二弧长间距沿所述第二区域的弧长方向对各所述第二区域进行区域划分，获得若干阵元区域。
可选的，所述第二高度间距Δh0满足：
Figure BDA0002289498320000021
Figure BDA0002289498320000022
其中，BWmax表示z轴方向最小波数支持的带宽；
fmax表示柱面阵列雷达装置的最大工作频率；
φ表示PO2与xy平面的夹角；
θ表示PO2与xy平面上的投影与x轴夹角；
P点表示目标点；
O2表示柱面阵列最上层阵元区域的圆平面圆心，其空间位置坐标为
Figure BDA0002289498320000031
H0表示柱面阵列雷达的高度。
可选的，所述第二弧长所述对应的柱面雷达的圆心角Δθ0满足：
Figure BDA0002289498320000032
其中，
Figure BDA0002289498320000033
b＝Rr；
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)，
Figure BDA0002289498320000066
表示第kA＝(i-1)NA+MA个阵元集合区域；
步骤S12：阵元布局；将步骤S11中的获取的每一个阵元集合区域Cell，按照步骤S1中的划分方式将各阵元集合区域划分成NS*MS个阵元区域，每个阵元集合区域中的阵元区域之间的间距满足采样定理，相应的阵元区域中放置阵元，具体步骤如下：
步骤S121：阵元集合区域高度划分；根据阵元集合区域的实际高度h0，利用公式(1)，将步骤S12中的每个阵元集合区域在z轴上，分成NS(NS＞1，NS为正整数)份，则在z轴上的高度间距为Δh0,
Figure BDA0002289498320000061
高度间距Δh0，需满足采样定理：
Figure BDA0002289498320000062
Figure BDA0002289498320000063
其中
Figure BDA0002289498320000064
表示Z方向最小波数支持带宽，
Figure BDA0002289498320000065
满足公式(3)，fmax表示柱面阵列雷达装置的最大工作频率，如图3和图4所示，O(0,0,0)表示柱面整列中间层阵元圆平面圆心，
Figure BDA0002289498320000071
表示柱面整列最上层阵元圆平面圆心，
Figure BDA0002289498320000072
表示柱面整列最上层阵元圆平面圆心，O'(0,0,zCell)表示所测阵元所处层的平面圆心，φ表示PO2与xy平面的夹角，θ表示PO2在xy平面上的投影与x轴夹角；
步骤S122：子阵角度划分；以z轴为中心，根据公式(4)将xy平面等分成MS(MS＞1，MS为正整数)份，每份的角度间距为Δθ0，角度间距Δθ0所对应的弧长即为第二弧长间距。
角度间距Δθ0需满足公式(5)，
Figure BDA0002289498320000074
其中b＝Rr，c表示光速，如图4所示，R表示等效相位中心到目标点P的距离，是一个系统参数，P点可以表示成P(rp,θp,zp)或P(xP,yP,zP)，rp表示P点到z轴的距离，zp表示P点到Oxy平面的距离；r表示目标点到原点的距离。
本实施例中，较佳的NS和MS最佳选取为|MS-NS|＝0，例如MS＝100个，NS最优选择为NS＝100。
步骤S123：获取阵元集合方阵；如图5所示，步骤S121和步骤S122将每个阵元集合区域分成NS×MS份等大的阵元区域，每个阵元区域的高度为Δh0，水平角度间距为Δθ0，将每个阵元集合区域构成的阵元集方阵用Cell(即
Figure BDA0002289498320000076
)表示，Cell是一个NS行MS列的阵元集合方阵，阵元方格
Figure BDA0002289498320000077
(阵元区域)表示阵元集Cell中第ic(ic≤NS ic为正整数)行第jc(jc≤MS,jc为正整数)列元素(即阵元)；
步骤S124：阵元布局；如图6所示，在Cell(即
Figure BDA0002289498320000081
)中的第一行、最后一行、第一列、最后一列阵元方格中布置阵元，Cell第一行阵元
Figure BDA0002289498320000082
(jc≤MS,jc为正整数)和最后一行阵元
Figure BDA0002289498320000083
用于发射信号，Cell第一列阵元(2≤ic≤NS-1ic为正整数)和最后一列阵元
Figure BDA0002289498320000085
用于接收信号，其它阵元方格不放置阵元，图6中标有“T”，表示该阵元用于发射阵元，标有“R”，表示该阵元用于接收阵元，其余阵元方格不用放置阵元，每个子阵(阵元集方阵)布置的阵元总数为2MS+2NS-4。
步骤S125：获取发射阵元矩阵；将步骤S124中的发射阵元构建1行2MS列的行矩阵
Figure BDA0002289498320000086
具体表达式如公式(6)，
Figure BDA0002289498320000087
表示Trans矩阵第kT(kT≤2Ms,k为正整数)个元素(阵元)，发射阵元的标号如图6所示；
Figure BDA0002289498320000088
步骤S126：获取接收阵元矩阵；将步骤S124中的接收阵元构建1行2(NS-2)列的行矩阵
Figure BDA0002289498320000089
具体表达式如公式(7)，
Figure BDA00022894983200000810
表示Rece矩阵第kR(kR≤2(Ns-2),kR为正整数)个元素(阵元)；
Figure BDA00022894983200000811
采用本发明中的布局方法，布置的阵元个数为(2Ms+2Ns-4)*NA*MA。而采用传统方式，柱面阵列雷达需要布置的阵元的个数为NS*MS*NA*MA，由此可见，本发明实施例中的天线阵元的布局方法能够大大减少了阵元的布置数量，节约了成本，同时能够提高数据采集的效率。
本发明实施例中在完成天线阵元的布置之后，就可以利用柱面阵列雷达来进行信号的发射和接收，以获得若干回波信号，具体获取回波信号的步骤包括：
步骤一、根据柱面阵列雷达的最小工作频率和最大工作频率获取目标频点矩阵；
本步骤中，在具体实施过程中具体包括：根据最小工作频率、最大工作频率，利用第一计算公式计算获得频点差值；
根据所述最小工作频率以及所述频点差值利用第二计算公式计算获得所述目标频点矩阵中的频点元素，以获得目标频点矩阵。
步骤二、同时循环选取各阵元集合方阵中的一个发射阵元作为第一发射阵元；
本实施中，各所述阵元集合方阵中分别设置有若干发射阵元和若干接收阵元。
步骤三、为各所述第一发射阵元分别依次分配所述目标频点矩阵中的各频点，以作为第一发射阵元的发射频点；其中，同一时刻不同阵元集合方阵中的第一发射阵元的发射频点不同；
步骤四、利用第一发射阵元根据分配的发射频点发射信号；
步骤五、利用接收阵元接收所述发射信号，以获得回波数据。
本步骤中，接收所述发射信号的接收阵元与发射所述发射信号的第一发射阵元位于相同的阵元集合方阵。
具体的，根据柱面阵列雷达的最小工作频率和最大工作频率获取目标频点矩阵，具体包括：
根据最小工作频率、最大工作频率，利用第一计算公式计算获得频点差值；根据所述最小工作频率以及所述频点差值利用第二计算公式计算获得所述目标频点矩阵中的频点元素，以获得目标频点矩阵。具体的所述第一计算公式为：其中，Δf为频点差值；NA×MA为阵元集合方阵的个数；Q为正整数；fmin为最小工作频率；fmax为最大工作频率。
本实施例中，所述第二计算公式为：其中
Figure BDA0002289498320000093
为第if个频点元素；Q为正整数，Q大于或等于阵元集合方阵的个数；fmin为最小工作频率；if为正整数；其中，所述最小工作频率fmin的范围为：fmin≥1GHz；所述最大工作频率fmax的范围为：fmax≤10THz。
具体的，在获得所述频点元素后，还包括：利用所述频点元素构建获得第一频点矩阵；利用若干所述第一频点矩阵构建获得目标频点矩阵。
本实施例在上述实施例的基础上，作为一种具体的实施方式，获取回波信号具体过程为：
步骤S2：频点设置，获取频点矩阵；将工作频段设置成多个频点，获取频点矩阵，具体如下：
步骤S21：频点设置；根据柱面阵列雷达最小工作频率fmin(fmin≥1GHz)和最大工作频率fmax(fmax≤10THz)设置Q(Q≥N*M，Q为正整数)个等差分布的频点，利用公式(8)(即第一计算公式)，计算Δf，将计算结果带入公式(9)，求出
Figure BDA0002289498320000102
Figure BDA0002289498320000103
用中所有元素构建一个Q×1的列矩阵
Figure BDA0002289498320000105
(即第一频点矩阵)；
步骤S22：获取目标频点矩阵Freq；利用步骤S21中F矩阵构建目标频点矩阵Freq，Freq矩阵大小为2Q×1，
Figure BDA0002289498320000106
(ifreq≤2Q，ifreq为正整数)表示频点矩阵Freq中第ifreq个元素(频点)，Freq表达式如下所示，
Figure BDA0002289498320000111
步骤S3：收发方式和策略；选取不同阵元集合方阵中的发射阵元，匹配阵元发射频点，发射、接收信号，流程图如图7所示；具体如下：
步骤S31：选取工作阵元；选取阵元集合方阵Aij中Trans的第kT个阵元
Figure BDA0002289498320000112
作为本次发射的发射阵元(即第一发射阵元)，其它发射阵元不工作，kT从1开始循环，循环选择发射阵元，直至kT＝2Ms才结束循环，每选一个阵元执行一次步骤S33，阵元发射不同频点的信号；
步骤S32：发射信号的频点匹配；根据阵元以及循环次数选取频点；工作的阵元所发射频点从目标频点矩阵Freq中选择；
根据步骤S31中选取Aij中
Figure BDA0002289498320000113
发射阵元，选取目标频点矩阵Freq中的元素
Figure BDA0002289498320000114
对应的频点满足公式(11),
ifreq＝(i-1)*MS+j+circulation-1 (11)
其中i和j表示阵元集合方阵A对应的阵元集合方阵Aij，circulation表示循环次数，circulation初始值为1，每次循环circulation＝circulation+1，直至循环到circulation＞Q，也就是每个阵元循环发射每个频点，同一时刻不同阵元发射的频点不同，不同时刻同一阵元发射的频点不同，保证每个阵元都发射所有不同大小的频点，而且仅发射一次；
步骤S33：阵元发射、接收信号；具体如下：
步骤S331：阵元发射信号；根据步骤S31选取Aij中
Figure BDA0002289498320000115
发射阵元、步骤S32选取的频点
Figure BDA0002289498320000116
发射信号，信号表达式如下：
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)，
Figure BDA0002289498320000121
图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的距离，如下所示；
Figure BDA0002289498320000122
Figure BDA0002289498320000123
接收阵元接收的信号可以写成公式(17)：
Figure BDA0002289498320000124
Figure BDA0002289498320000131
其中，A表示阵元集合方阵；Trans表示阵元集合方阵A中的发射阵元；Rece表示阵元集合方阵A中的接收阵元；θT表示发射阵元的空间位置坐标中的方位角；θR表示接收阵元的空间位置坐标中的方位角；R0表示柱面阵列的底面半径；x表示目标三维像素点的空间位置坐标中的x轴的坐标；y表示目标三维像素点的空间位置坐标中的y轴的坐标；z表示目标三维像素点的空间位置坐标中的z轴的坐标；Kω表示波数；f表示发射频点；c表示光速。
本实施例中通过采用上述的信号收方式来获得回波信号，能够减少数据的采集时间，为后续快速成像提供了基础。
本发明实施例在上述实施例的基础上，在获得了回波信号之后就可以进行三维成像，具体成像的过程为：
步骤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：计算匹配滤波函数(即每个接收阵元接收的回波信号对应一个滤波函数)；循环选取的阵元集合方阵
Figure BDA0002289498320000132
发射阵元接收阵元以及像素点Im，根据选取的阵元集合方阵
Figure BDA0002289498320000135
发射阵元
Figure BDA0002289498320000136
接收阵元
Figure BDA0002289498320000141
及像素点Im，计算匹配滤波函数
Figure BDA0002289498320000142
(A,Trans,Rece,Im)，具体如下：
步骤S421：系统初始化；对像素点、阵元集合方阵
Figure BDA0002289498320000143
发射阵元Trans以及接收阵元Rece进行初始化操作，具体如下：
步骤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：计算匹配滤波函数；根据选取的阵元集合方阵、发射阵元、接收阵元以及像素点计算匹配滤波函数，所述滤波函数的表达式为：
Figure BDA0002289498320000151
其中：
Figure BDA0002289498320000153
Figure BDA0002289498320000154
其中，A表示阵元集合方阵；Trans表示阵元集合方阵A中的发射阵元；Rece表示阵元集合方阵A中的接收阵元；Im表示第m个像素点；RTI表示发射阵元与目标三维像素点之间的距离；RRI表示接收阵元与目标三维像素点之间的距离；Kω表示波数；f表示发射频点；c表示光速。
步骤S43：将匹配滤波函数
Figure BDA0002289498320000155
(A,Trans,Rece,Im)与回波数据S(A,Trans,Rece,Kω)相乘，得到匹配信号
Figure BDA0002289498320000156
(A,Trans,Rece,Kω,Im)，具体如下：
步骤S44：匹配信号逆傅里叶变换；将匹配信号
Figure BDA0002289498320000158
(A,Trans,Rece,Kω,Im)进行逆傅里叶变换，得到距离向压缩信号
Figure BDA0002289498320000159
(A,Trans,Rece,Im)，具体如下：
Figure BDA00022894983200001510
步骤S45：取距离向压缩信号的峰值；距离向压缩信号
Figure BDA00022894983200001511
(A,Trans,Rece,Im)是一个sinc函数，且距离向压缩信号的峰值
Figure BDA0002289498320000161
(A,Trans,Rece,Im)位于距离向压缩信号第一个采样点处，具体如下；
Figure BDA0002289498320000162
|first表示距离向压缩信号(A,Trans,Rece,Im)第一个采样点；
步骤S46：距离向压缩信号峰值处理；根据方位向合成孔径和高程向合成孔径范围，对距离向压缩信号的峰值
Figure BDA0002289498320000164
(A,Trans,Rece,Im)进行处理，具体如下：
步骤S461：计算方位向筛选条件；方位向筛选条件如下所示，
Figure BDA0002289498320000165
angle(um+vmi)函数用来求解复数um+vmi在不同象限中的相角；θA表示合成孔径方位向范围，是系统参数；其中，θT表示发射阵元的空间位置坐标中的方位角；θR表示接收阵元的空间位置坐标中的方位角；um表示第m个三维像素点的空间位置坐标中的x轴的坐标；vm表示第m个三维像素点的空间位置坐标中的y轴的坐标。
步骤S462：计算高程向筛选条件；高程向筛选条件如下所示，
Figure BDA0002289498320000166
θH表示合成孔径高程向范围，是系统参数；um表示第m个三维像素点的空间位置坐标中的x轴的坐标；vm表示第m个三维像素点的空间位置坐标中的y轴的坐标；wm表示第m个三维像素点的空间位置坐标中的z轴的坐标；R0表示柱面阵列的底面半径；zT表示发射阵元的空间位置坐标中的z轴的坐标；zR表示接收阵元的空间位置坐标中的z轴的坐标；
步骤S463：距离向压缩信号的峰值处理；根据方位向筛选条件和高程向筛选条件，令不符合条件的距离向压缩信号的峰值
Figure BDA0002289498320000171
(A,Trans,Rece,Im)等于0，具体如下：
Figure BDA0002289498320000172
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＝[um vm wm σsum(Im)]矩阵，其中第一列表示像素点Im的x横坐标，第二列表示像素点Im的y坐标，第三列表示像素点Im的z坐标，第四列表示像素点Im的散射强度，利用Tar_Sca矩阵绘出柱面阵列雷达三维复图像。
本实施例中的柱面阵列雷达成像方法，能够实现360°全方位探测，并且所采用的阵元相对较少，降低成本，还能减少数据量，提高了成像效率。
以上实施例仅为本发明的示例性实施例，不用于限制本发明，本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内，对本发明做出各种修改或等同替换，这种修改或等同替换也应视为落在本发明的保护范围内。

Claims (6)

1.一种基于柱面雷达成像装置的稀疏天线阵元布局方法，其特征在于，包括如下步骤：
对柱面雷达进行区域划分，获得若干阵元集合区域；
对各所述阵元集合区域进行区域划分获得若干阵元区域；
确定各所述阵元集合区域中位于首行和尾行的阵元区域为发射阵元区域；
确定各所述阵元集合区域中位于首列和尾列的、除了所述发射阵元区域之外的阵元区域为接收阵元区域；
在所述发射阵元区域处设置发射阵元；
在所述接收阵元区域处设置与所述发射阵元对应的接收阵元。
2.如权利要求1所述的方法，其特征在于，所述对柱面阵列雷达进行区域划分，获得若干阵元集合区域，具体包括：
确定柱面阵列雷达的高度；
按照预定的第一高度间距沿着所述柱面阵列雷达的高度方向对柱面阵列雷达的柱面进行区域划分，获得若干个第一区域；
按照预定的第一弧长沿第一区域的周长方向对所述第一区域进行区域划分，获得若干个阵元集合区域。
3.如权利要求1所述的方法，其特征在于，所述对各所述阵元集合区域进行区域划分获得若干阵元区域，具体包括：
确定各所述阵元集合区域的高度和弧长；
按照预定的第一份数确定第二高度间距；
按照预定的第二份数确定第二弧长间距；
按照所述第二高度间距沿着所述柱面阵列雷达的高度方向对阵元集合区域进行区域划分，获得若干第二区域；
按照所述第二弧长间距沿所述第二区域的弧长方向对各所述第二区域进行区域划分，获得若干阵元区域。
4.如权利要求3所述的方法，其特征在于，所述第二高度间距Δh0满足：
Figure FDA0002289498310000011
Figure FDA0002289498310000012
其中，
Figure FDA0002289498310000024
表示z轴方向最小波数支持的带宽；
fmax表示柱面阵列雷达装置的最大工作频率；
φ表示PO2与xy平面的夹角；
θ表示PO2与xy平面上的投影与x轴夹角；
P点表示目标点；
O2表示柱面阵列最上层阵元区域的圆平面圆心，其空间位置坐标为
Figure FDA0002289498310000021
H0表示柱面阵列雷达的高度。
5.如权利要求1所述的方法，其特征在于，所述第二弧长所述对应的柱面雷达的圆心角Δθ0满足：
Figure FDA0002289498310000022
其中，b＝Rr；
c表示光速表示等效相位中心到目标点的距离，是一个系统参数；
R表示等效相位中心到目标点P的距离，是一个系统参数；
P点的空间位置坐标表示为P(rp,θp,zp)，其中rp表示P点到z轴的距离；
zp表示P点到Oxy平面的距离；
r表示目标点到原点的距离。
6.如权利要求3所述的方法，其特征在于，所述第一份数等于所述第二份数。
CN201911174015.9A 2019-11-26 2019-11-26 一种基于柱面雷达成像装置的稀疏天线布局方法 Active CN110837127B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911174015.9A CN110837127B (zh) 2019-11-26 2019-11-26 一种基于柱面雷达成像装置的稀疏天线布局方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911174015.9A CN110837127B (zh) 2019-11-26 2019-11-26 一种基于柱面雷达成像装置的稀疏天线布局方法

Publications (2)

Publication Number Publication Date
CN110837127A true CN110837127A (zh) 2020-02-25
CN110837127B CN110837127B (zh) 2021-09-10

Family

ID=69577554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911174015.9A Active CN110837127B (zh) 2019-11-26 2019-11-26 一种基于柱面雷达成像装置的稀疏天线布局方法

Country Status (1)

Country Link
CN (1) CN110837127B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111812642A (zh) * 2020-05-25 2020-10-23 北京理工大学 一种圆柱孔径mimo阵列天线、成像方法及补偿方法
CN112799064A (zh) * 2020-12-30 2021-05-14 内蒙古工业大学 柱面孔径非线性渐进式相位迭代成像的方法及装置
CN112835040A (zh) * 2020-12-30 2021-05-25 内蒙古工业大学 球面孔径分区域渐进式相位迭代成像的方法及装置
CN112835039A (zh) * 2020-12-30 2021-05-25 内蒙古工业大学 平面孔径分区域非线性渐进式相位迭代成像的方法及装置

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040090359A1 (en) * 2001-03-16 2004-05-13 Mcmakin Douglas L. Detecting concealed objects at a checkpoint
WO2008054350A2 (en) * 2005-07-12 2008-05-08 Trex Enterprises Corporation Multi-channel millimeter wave image system
CN102866401A (zh) * 2012-08-06 2013-01-09 西北工业大学 一种基于mimo技术的三维成像方法
CN103454631A (zh) * 2013-08-31 2013-12-18 西北工业大学 一种基于mimo技术和稀疏阵技术的三维成像方法
CN103645469A (zh) * 2013-12-18 2014-03-19 中国科学院电子学研究所 确定微波三维成像中相位中心时空分布的方法
CN104269612A (zh) * 2014-10-15 2015-01-07 内蒙古工业大学 近距离平面阵列多输入多输出成像天线布局方法
CN104269658A (zh) * 2014-10-21 2015-01-07 内蒙古工业大学 用于mimo-sar成像的弧形阵列天线
CN105891827A (zh) * 2015-10-30 2016-08-24 中国人民解放军空军工程大学 一种机载mimo-sar下视三维稀疏成像方法
CN107219527A (zh) * 2017-05-27 2017-09-29 吉林大学 一种周期型双阵列通道式成像系统的单快拍快速成像方法
CN107422310A (zh) * 2017-09-05 2017-12-01 芜湖华创光电科技有限公司 一种用于方位与俯仰二维测向的稀疏阵列设计方法
CN206848472U (zh) * 2017-05-27 2018-01-05 吉林大学 周期型双阵列通道式毫米波主动三维运动成像系统
CN109281651A (zh) * 2017-07-19 2019-01-29 中国科学院声学研究所 一种应用于柱面超声阵列的超声井壁成像方法
US20190057957A1 (en) * 2016-12-21 2019-02-21 PhantaField, Inc. Augmented reality display system
CN109782367A (zh) * 2019-01-14 2019-05-21 河北华讯方舟太赫兹技术有限公司 一种安检系统以及安检方法
CN109884622A (zh) * 2019-03-11 2019-06-14 内蒙古工业大学 柱面阵列天线三维成像的方法
CN110045367A (zh) * 2019-03-11 2019-07-23 内蒙古工业大学 柱面阵列天线目标体三维成像的装置

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040090359A1 (en) * 2001-03-16 2004-05-13 Mcmakin Douglas L. Detecting concealed objects at a checkpoint
WO2008054350A2 (en) * 2005-07-12 2008-05-08 Trex Enterprises Corporation Multi-channel millimeter wave image system
CN102866401A (zh) * 2012-08-06 2013-01-09 西北工业大学 一种基于mimo技术的三维成像方法
CN103454631A (zh) * 2013-08-31 2013-12-18 西北工业大学 一种基于mimo技术和稀疏阵技术的三维成像方法
CN103645469A (zh) * 2013-12-18 2014-03-19 中国科学院电子学研究所 确定微波三维成像中相位中心时空分布的方法
CN104269612A (zh) * 2014-10-15 2015-01-07 内蒙古工业大学 近距离平面阵列多输入多输出成像天线布局方法
CN104269658A (zh) * 2014-10-21 2015-01-07 内蒙古工业大学 用于mimo-sar成像的弧形阵列天线
CN105891827A (zh) * 2015-10-30 2016-08-24 中国人民解放军空军工程大学 一种机载mimo-sar下视三维稀疏成像方法
US20190057957A1 (en) * 2016-12-21 2019-02-21 PhantaField, Inc. Augmented reality display system
CN206848472U (zh) * 2017-05-27 2018-01-05 吉林大学 周期型双阵列通道式毫米波主动三维运动成像系统
CN107219527A (zh) * 2017-05-27 2017-09-29 吉林大学 一种周期型双阵列通道式成像系统的单快拍快速成像方法
CN109281651A (zh) * 2017-07-19 2019-01-29 中国科学院声学研究所 一种应用于柱面超声阵列的超声井壁成像方法
CN107422310A (zh) * 2017-09-05 2017-12-01 芜湖华创光电科技有限公司 一种用于方位与俯仰二维测向的稀疏阵列设计方法
CN109782367A (zh) * 2019-01-14 2019-05-21 河北华讯方舟太赫兹技术有限公司 一种安检系统以及安检方法
CN109884622A (zh) * 2019-03-11 2019-06-14 内蒙古工业大学 柱面阵列天线三维成像的方法
CN110045367A (zh) * 2019-03-11 2019-07-23 内蒙古工业大学 柱面阵列天线目标体三维成像的装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FORTUNY-GUASCH等: ""Extension of the 3-D range migration algorithm to cylindrical and spherical scanning geometries"", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 *
叶恺 等: ""基于柱形抛物面天线的MIMO SAR研究"", 《电子与信息学报》 *
温馨 等: ""主动式毫米波近距离圆柱扫描三维成像系统"", 《系统工程与电子技术》 *
田鹤 等: ""频域稀疏毫米波人体安检成像处理和快速成像稀疏阵列设计"", 《雷达学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111812642A (zh) * 2020-05-25 2020-10-23 北京理工大学 一种圆柱孔径mimo阵列天线、成像方法及补偿方法
CN112799064A (zh) * 2020-12-30 2021-05-14 内蒙古工业大学 柱面孔径非线性渐进式相位迭代成像的方法及装置
CN112835040A (zh) * 2020-12-30 2021-05-25 内蒙古工业大学 球面孔径分区域渐进式相位迭代成像的方法及装置
CN112835039A (zh) * 2020-12-30 2021-05-25 内蒙古工业大学 平面孔径分区域非线性渐进式相位迭代成像的方法及装置
CN112835039B (zh) * 2020-12-30 2023-09-08 内蒙古工业大学 平面孔径分区域非线性渐进式相位迭代成像的方法及装置

Also Published As

Publication number Publication date
CN110837127B (zh) 2021-09-10

Similar Documents

Publication Publication Date Title
CN110837127B (zh) 一种基于柱面雷达成像装置的稀疏天线布局方法
CN110837128B (zh) 一种柱面阵列雷达的成像方法
CN103682677B (zh) 一种艇载雷达共形稀疏阵列天线及其信号处理方法
CN105759262B (zh) 一种基于mimo体制的近程三维成像装置及其成像方法
CN112099011B (zh) 基于focuss算法的全息sar子孔径三维重建方法
US9638795B2 (en) Device for the measurement of microwave signals and a method for the configuration of the same
CN109765562A (zh) 一种三维前视声像声纳系统和方法
CN110794471B (zh) 一种毫米波稀疏阵列远程监视成像方法及系统
CN103731189B (zh) 共形天线阵动态子阵划分方法及波达方向估计方法
CN104133213A (zh) 一种结合rm算法与bp算法的柱面近场三维rcs成像方法
CN209821372U (zh) 用于主动式微波毫米波安检设备的电磁成像装置和安检设备
CN102520408B (zh) 一种圆柱阵面三维成像系统的三维成像方法
CN107153179B (zh) 一种雷达目标rcs与散射中心同步测试方法
CN112859072A (zh) 一种基于平面mimo阵列的雷达系统及三维成像方法
CN103675817A (zh) 一种基于变换域稀疏的合成孔径雷达侧视三维成像方法
CN111142164B (zh) 一种柱面雷达成像系统
CN110837126B (zh) 一种用于柱面阵列雷达成像的信号收发方法及装置
CN112649806B (zh) 一种mimo雷达近场三维成像方法
WO2020147437A1 (zh) 一种天线面阵、安检装置以及安检方法
CN108717188B (zh) 适用于mimo穿墙雷达成像的快速墙体补偿方法
CN113359196B (zh) 基于子空间法和dbf的多目标生命体征探测方法
CN108761457A (zh) 基于mimo阵列合成孔径的高精度三维快速成像方法及装置
CN107678009A (zh) 干涉处理艇载雷达阵列形变误差补偿和目标探测方法
CN109613533B (zh) 双站微波凝视关联成像方法及装置、存储介质及电子设备
Craddock et al. Development of a hemi-spherical wideband antenna array for breast cancer imaging

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