[go: up one dir, main page]

CN106556874B - 一种近距离微波成像方法及系统 - Google Patents

一种近距离微波成像方法及系统 Download PDF

Info

Publication number
CN106556874B
CN106556874B CN201610931503.XA CN201610931503A CN106556874B CN 106556874 B CN106556874 B CN 106556874B CN 201610931503 A CN201610931503 A CN 201610931503A CN 106556874 B CN106556874 B CN 106556874B
Authority
CN
China
Prior art keywords
echo data
grid
echo
imaging
data
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
Application number
CN201610931503.XA
Other languages
English (en)
Other versions
CN106556874A (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.)
Shenzhen Zhongtou Huaxun Terahertz Technology Co ltd
Original Assignee
Shenzhen Institute of Terahertz Technology and Innovation
Shenzhen Huaxun Ark Technology Co Ltd
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 Shenzhen Institute of Terahertz Technology and Innovation, Shenzhen Huaxun Ark Technology Co Ltd filed Critical Shenzhen Institute of Terahertz Technology and Innovation
Priority to CN201610931503.XA priority Critical patent/CN106556874B/zh
Publication of CN106556874A publication Critical patent/CN106556874A/zh
Priority to US16/345,971 priority patent/US11215730B2/en
Priority to PCT/CN2017/096099 priority patent/WO2018076884A1/zh
Application granted granted Critical
Publication of CN106556874B publication Critical patent/CN106556874B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/38—Processing data, e.g. for analysis, for interpretation, for correction
    • 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
    • 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/2806—Employing storage or delay devices which preserve the pulse form of the echo signal, e.g. for comparing and combining echoes received during different periods
    • 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
    • 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

Landscapes

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

Abstract

本发明提供一种近距离微波成像方法及系统，通过控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描目标区域，在圆弧轨道上的每个方位角位置点，控制线性阵列天线采集预设个数的回波数据，并将预设个数的回波数据构成的回波数据集发送给信号处理设备，直到线性阵列天线在圆弧轨道上的预设方位角位置点完成回波数据的采集时止；控制信号处理设备每接收到一次回波数据集，则对回波数据集进行成像处理，可以有效提升近距离微波成像的处理速度、缩短处理时间，实现实时成像。

Description

一种近距离微波成像方法及系统
技术领域
本发明实施例属于近距离微波成像技术领域，尤其涉及一种近距离微波成像方法及系统。
背景技术
近距离微波成像技术采用主动发射微波波段的电磁波的方式对物体进行扫描成像，因其能够穿透物体表面而检测藏匿于物体内部的金属或非金属违禁品，且由于其具有辐射剂量小、属于非电离辐射等特点，广泛应用于安检系统，用于执行人体安全检查任务。
然而，目前较为流行的近距离微波成像技术，均要求首先得到所有等效天线采集位置处的回波数据组成的回波数据集，再传送将回波数据集发送至安检系统的信号处理设备进行统一成像处理，才能得到目标点物体的图像，无法在数据采集的过程中同时进行成像处理，数据处理效率低下。
发明内容
本发明实施例提供一种近距离微波成像方法及系统，能够在数据采集的过程中同时进行成像处理，可以有效提升近距离微波成像的处理速度、缩短处理时间，实现实时成像。
本发明实施例一方面提供一种近距离微波成像方法，其包括：
将目标区域划分为预设分辨率大小的栅格集；
建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果；
控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止；
控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
本发明实施例另一方面还提供一种近距离微波成像系统，其包括：
栅格化模块，用于将目标区域划分为预设分辨率大小的栅格集；
分区存储模块，用于建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果；
数据采集模块，用于控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止；
数据处理模块，用于控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
本发明实施例通过控制信号处理设备每接收到一次回波数据集，则对所述回波数据集进行成像处理，可以有效提升近距离微波成像的处理速度、缩短处理时间，实现实时成像。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明的一个实施例提供的笛卡尔坐标系下的柱面扫描系统的示意图。
图2是本发明的一个实施例提供的近距离微波成像方法的流程框图；
图3是本发明的一个实施例提供的步骤S40的流程框图；
图4是本发明的一个实施例提供的近距离微波成像系统的结构框图；
图5是本发明的一个实施例提供的数据处理模块的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“包括”以及它们任何变形，意图在于覆盖不排他的包含。例如包含一系列步骤或单元的过程、方法或系统、产品或设备没有限定于已列出的步骤或模块，而是可选地还包括没有列出的步骤或模块，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或模块。
本发明所有实施例基于近距离微波成像系统实现，该系统由线性排列的预设个数的天线组成的线性阵列天线、信号收发设备和机械旋转设备和信号处理设备。
在具体应用中，信号处理设备可以采用图形处理器(Graphics Processing Unit，GPU)。
机械旋转设备控制线性阵列天线旋转，信号收发设备通过线性阵列天线在旋转过程中发射信号并接收待成像物体反射的回波数据，信号收发设备将回波数据发送给信号处理设备处理为成像结果实现对待成像物体上各采样点的扫描成像。根据线性阵列天线旋转过程中的等效相位位置所形成的扫描面的形状，可以将近距离微波成像系统分为平面扫描系统或柱面扫描系统。
以下着重介绍柱面扫描系统的工作原理：
如图1所示，设定待成像物体为人体，人体自然站立的平面为水平面，设定与X轴正交且处于同一水平面的轴为Y轴，设定X轴、Y轴和Z轴相交于所述集合中轴线上的一点O，以O为原点建立包括X轴、Y轴和Z轴的笛卡尔坐标系X-Y-Z-O。
N个天线在垂直于水平面的Z轴方向上组成线性阵列天线，线性阵列天线距柱面扫描系统的几何中轴线的距离为ρ，该线性阵列天线绕所述几何中轴线旋转一周后，等效形成在柱面上分布的等方位角间隔、等高度间隔的多个等效天线位置，该等效天线位置称为天线的等效相位中心。
以柱面坐标表示某个天线等效相位中心的位置，其中表示以等方位角间隔采样的条件下，第m个方位角位置处的采样点的方位角坐标值，其中，m的范围为0≤m≤M-1，M表示同一天线绕所述几何中轴线旋转一周的采样点的总数，方位角是指某一天线的X轴坐标的变化角度；zn表示在Z轴方向上等间隔采样的条件下，Z轴方向上第n个采样点的Z轴坐标值，其中，n的范围为0≤n≤N-1，N表示线性阵列天线在Z轴方向上的采样点的总数。在某个方位角坐标值处，天线的微波信号的发射和接收是通过控制微波开关的通断电来实现的，线性阵列天线所包括的第0个天线的等效相位中心依次向第N-1个天线的等效相位中心扫描时，微波开关的切换速度在纳秒级别，信号的发射和接收的时间在几十纳秒级别，因此，线性阵列天线在扫描过程中不同天线之间的扫描时间差异可以忽略不计，可以认为线性天线阵列的所有天线均是同时发射和接收信号。
在具体应用中，机械旋转设备控制线性阵列天线沿着与线性阵列天线正交的圆弧轨道旋转运动，以对待测物体进行扫描，扫描一周即完成0～2π(包括端点值)角度的旋转，所需时间为5～6秒钟。
设定，某一天线的等效相位中心处发射的微波信号st(t)为线性调频信号，其表达式为：
其中，t表示时间，A表示发射的信号幅度，Ts表示发射的信号宽度，fc表示发射的信号的中心频率，B表示发射的信号的带宽。
假设人体目标区域中某个位置处的目标点P在笛卡尔坐标系中的坐标为(xP,yP,zP)，其散射系数表示为σ(xP,yP,zP)，则对于某个柱面坐标处的天线的等效相位中心，微波信号st(t)照射在目标点P上，经目标点P散射后被坐标处的天线接收的信号表示为：
其中，0≤m≤M-1，0≤n≤N-1，A′表示接收的信号的幅度，C表示电磁波在真空中传播速度，F时表示目标点P是否在当前天线照射范围内的因子，当目标点P处于当前天线的照射范围时，F值取1，否则F值取0，V表示目标点P属于人体目标区域，R表示目标点P和当前天线的等效相位中心的瞬时斜距，R的表达式为：
接收到的信号sr(t)在信号收发设备中经过模拟解线性频调(analogue de-chirp)操作后，表示为如下形式：
其中，f表示信号的频率，其范围是f∈[fc-B/2,fc+B/2]。
上述信号再经过模数采样后即为当前柱面坐标处的天线的等效相位中心所采集到的一条回波数据，设模数采样后共在信号频率轴f上采集了L个采样点，则在柱面坐标处的天线的等效相位中心位置采集到的回波数据的表达式为：
其中，fl表示第l个频点的频率。通常，三维成像处理的目的是通过多条回波数据去反演或重建人体目标区域V上所有目标点的散射系数的空间分布σ(xP,yP,zP)，其中，(xP,yP,zP)∈V。
在固定的方位角坐标值处，可以获得线性阵列天线上N个天线的N个等效相位中心位置采集到的由N条回波数据组成的回波数据集
其中l＝0,1,...,L-1。
现有技术中，通常通过获取所有天线的效相位中心位置采集到的回波数据组成的数据集然后对数据集中所有回波数据所对应的采样点的时间频率信息进行成像处理，从而生成人体目标区域的三维图像，耗费的数据处理时间较长。
本发明实施例在于提供一种近距离微波成像方法和系统，通过将线性阵列天线在每个方位角位置处采集得到的回波数据集
及时发送给信号处理设备进行成像处理，不需要等到所有数据采集完成再集中进行成像处理，增加了数据处理的实时性，提高了工作效率，节约了数据处理时间。
如图2所示，本发明的一个实施例提供的一种近距离微波成像方法，其包括：
步骤S10：将目标区域划分为预设分辨率大小的栅格集。
在本实施例中，目标区域是指待成像物体所在区域，该区域的大小具体由近距离微波成像系统的线性阵列天线旋转一周所能覆盖的扫描区域的大小决定，所述栅格集中的栅格即为待成像的目标点。
具体应用中，待成像物体可以是人体，也可以是其他有机体或无机体，待成像物体可以是一维、二维或者三维物体，对应的所述栅格集的维度也可以是一维、二维或者三维。
步骤S20：建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果。
在具体应用中，预设分辨率可以根据实际的成像精度要求进行设置。
本实施例中，通过预先建立存储区域，方便对后续过程中信号处理设备每次处理得到的成像结果进行保存。
在具体应用中，矩阵中各元素的初始值设置为0。
步骤S30：控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止。
在本实施例中，所述圆弧轨道为与所述线性阵列天线正交的圆弧轨道，所述预设个数的回波数据为由预设个数的栅格反射的回波信号经模数转换得到的数据。
步骤S40：控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
在具体应用中，线性阵列天线采集每个栅格的回波数据，信号处理设备对每个栅格的回波数据进行计算处理，得到以与每个栅格对应的矩阵元素的值的形式表示的成像结果。
本发明实施例通过控制信号处理设备每接收到一次回波数据集，则对所述回波数据集进行成像处理，可以有效提升近距离微波成像的处理速度、缩短处理时间，实现实时成像。
如图2所示，在本发明的一个实施例中，步骤S40具体包括：
步骤S401：分别对所述回波数据集中的每个回波数据进行快速傅里叶变换。
在具体应用中，所述快速傅里叶变换可以为并行快速傅里叶变换，即在每一次获取到回波数据集时，同时对该回波数据集中的所有回波数据进行快速傅里叶变换。
步骤S402：根据所述快速傅里叶变换的结果，通过后向投影算法计算得到与所述每个回波数据对应的元素的值。
在具体应用中，所述后向投影算法可以为并行后向投影算法，即同时对所述回波数据集所述包括的所有回波数据的快速傅里叶变换结果进行后向投影算法计算，得到与所述每个回波数据对应的元素的值，该值即为与相应元素对应的栅格的成像结果。
在一个实施例中，步骤S402具体包括：
分别对所述回波数据集中的每个回波数据沿频率轴方向进行快速傅里叶变换，得到与所述回波数据集对应的数据集，所述数据集包括与所述回波数据集中的每个回波数据对应的数据；
步骤S403：根据每个元素的值更新所述每个元素的初始值。
在具体应用中，进行近距离微波成像之前，可将用于表示成像结果的矩阵中的所有元素的值初始化为0，在后续计算得到每个元素的值时，再将每个元素的初始值更新为计算得到的值。
在一个实施例中，步骤S403具体包括：
根据所述数据集，通过基于插值法的后向投影算法计算得到与所述每个回波数据对应的元素的值。
本发明实施例通过对回波数据进行并行处理，可以进一步提高数据处理效率，减少处理时间。
在本发明的一个实施例中，将图1所对应的实施例所提供的近距离微波成像方法用于获取人体的三维成像结果时，该方法具体可以包括如下步骤：
步骤S1：将人体目标区域划分为分辨率大小为NX×NY×NZ的三维栅格，其中，NX、NY和NZ分别表示所述三维栅格在笛卡尔坐标系下平行于水平面的横轴(本实施例中表示为X轴)方向上的栅格数量、纵轴(本实施例中表示为Y轴)方向上的栅格数量和垂直于所述水平面方向的竖轴(本实施例中表示为Z轴)方向上的栅格数量。
在具体应用中，所述水平面是指人体自然站立时、垂直于人体高度方向上的中轴线的平面。
在一个实施例中，步骤S1中三维栅格划分的方法可以包括以下步骤：
步骤S11：通过用户预先设置的划分规则或者用户当前输入的划分规则，确定笛卡尔坐标系下人体目标区域的尺寸。
本实施例中，人体目标区域在X轴、Y轴和Z轴三个方向上的待成像的尺寸分别表示为ΔX、ΔY和ΔZ。通常情况下，自然站立的人体位于柱面扫描系统的圆弧轨道的圆心所在的中轴线上，即柱面扫描系统的几何中轴线。
本实施例中，定义X轴是水平面上由人体的背部指向胸部的方向，在具体应用中，ΔX可以设置为大于人体厚度的数值，例如0.5米；定义Y轴是水平面上由人体一侧肩部指向另一侧肩部的方向，在具体应用中，ΔY可以设置为大于人体宽度的数值，例如1.0米；Z轴是垂直于水平面且穿过圆弧轨道圆心的方向，在具体应用中，ΔY可以设置为大于人体身高的数值，例如2.2米。
步骤S12：设定三维栅格集中的栅格在X轴、Y轴和Z轴方向上的尺寸δX、δY和δZ。
在具体应用中，δX＝δY＝δZ＝0.5厘米；
步骤S13：计算得到三维栅格集在X轴、Y轴和Z轴方向上的栅格个数NX、NY和NZ。
在具体应用中，三维栅格集在X轴、Y轴和Z轴方向上的栅格个数分别为
和其中符号表示向上取整算子。
步骤S13：根据NX、NY和NZ将人体目标区域划分为分辨率大小为NX×NY×NZ的三维栅格集，三维栅格集中索引为(i,j,k)的栅格在笛卡尔坐标系X-Y-Z-O下的坐标为：
步骤S2：根据步骤S1中的三维栅格划分结果，控制信号处理设备划分一块存储区域用于保存人体目标区域的三维成像结果，三维成像结果表示为分辨率大小为NX×NY×NZ的三维矩阵I，NX、NY和NZ分别表示所述三维矩阵I在横向、纵向和竖向上的元素个数，将三维矩阵I所有元素的值均初始化为0，即I(i,j,k)＝0，其中1≤i≤NX，1≤j≤NY，1≤k≤NZ。
在具体应用中，三维栅格集与三维矩阵I中脚标相同的栅格和元素一一对应。
步骤S3：械旋转设备控制线性阵列天线沿着与线性阵列天线正交的圆弧轨道旋转运动，在方位角位置处，控制所述线性阵列天线采集得到由该线性阵列天线包含的N个天线的N个等效相位中心位置处的由共N条回波数据组成的回波数据集将回波数据集Dm发送给信号处理设备，控制线性阵列天线继续沿所述圆弧轨道进行机械旋转运动，至下一个方位角位置处采集回波数据，直至圆弧轨道上预设的第M-1个方位角位置处的回波数据被采集时止；、传输、处理完毕，从而最终的人体目标点成像结果保存在三维成像结果矩阵I中；
步骤S4：控制所述信号处理设备接收回波数据集
并对回波数据集中的每条回波数据进行快速傅里叶逆变换(IFFT)，再通过后向投影方式计算三维矩阵I中每个元素的值I(i,j,k)，其中1≤i≤NX，1≤j≤NY，1≤k≤NZ，且i、j和k均为整数，根据计算的到每个元素的值I(i,j,k)更新该元素的初始值。
在一个实施例中，步骤S4具体包括：
步骤S41：控制所述信号处理设备接收回波数据集
并对回波数据集中的每条回波数据沿其信号频率轴方向进行快速傅里叶逆变换(IFFT)，例如，对第n条回波数据进行IFFT后得到新的数据
其中，n＝0,1,...,N，表示回波数据在其信号频率轴方向上的快速傅里叶逆变换算子，tl表示回波数据在时间域上的采样轴，l＝0,1,...,L，其中L表示快速傅里叶变换之前，每条回波数据在其信号频率轴方向上的采样点的总数；
步骤S42：对回波数据集Dm中的共N条回波数据通过快速傅里叶变换后，将得到的所有新的数据以集合形式表示为数据集D′m，
步骤S43：根据新的数据集D′m中的数据，通过后向投影公式
计算并更新三维矩阵I中每个元素的值I(i,j,k)；
其中，符号“+＝”表示将该符号右端的值与该符号左端变量的值相加作为该符号左端变量的更新值，R′表示柱面坐标为的等效相位总线与索引为(i,j,k)的栅格的几何中心的连线长度，C表示电磁波在真空中传播速度，fc表示天线发射的信号的中心频率。
在一个实施例中，步骤S43中的后向投影公式中的的值是通过对快速傅里叶变换后得到的新的数据集D′m中第n条数据进行插值得到。优选地，具体可以通过三次样条插值方式实现。
如图4所示，本发明的一个实施例提供一种近距离微波成像系统100，该系统用于执行图2所对应的实施例中的方法步骤，该系统包括：
栅格化模块10，用于将目标区域划分为预设分辨率大小的栅格集；
分区存储模块20，用于建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果；
数据采集模块30，用于控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止；
数据处理模块40，用于控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
本发明实施例通过控制信号处理设备每接收到一次回波数据集，则对所述回波数据集进行成像处理，可以有效提升近距离微波成像的处理速度、缩短处理时间，实现实时成像。
如图5所示，在本发明的一个实施例中，数据处理模块40用于执行图3所对应的实施例中的方法步骤，该数据处理模块40包括：
傅里叶变换单元401，用于分别对所述回波数据集中的每个回波数据进行快速傅里叶变换；
后向投影计算单元402，用于根据所述快速傅里叶变换的结果，通过后向投影算法计算得到与所述每个回波数据对应的元素的值；
更新单元403，用于根据每个元素的值更新所述每个元素的初始值。
在一个实施例中，傅里叶变换单元401具体用于：
分别对所述回波数据集中的每个回波数据沿频率轴方向进行快速傅里叶变换，得到与所述回波数据集对应的数据集；所述数据集包括与所述回波数据集中的每个回波数据对应的数据；
后向投影计算单元402具体用于：
根据所述数据集，通过基于插值法的后向投影算法计算得到与所述每个回波数据对应的元素的值。
在具体应用中，所述插值法可以为样条插值法，所述快速傅里叶变换可以为为并行快速傅里叶变换，所述后向投影算法可以为并行后向投影算法。
本发明所有实施例中的模块或单元，可以通过通用集成电路，例如CPU(CentralProcessing Unit，中央处理器)，或通过ASIC(Application Specific IntegratedCircuit，专用集成电路)来实现。
本发明实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减。
本发明实施例装置中的模块或单元可以根据实际需要进行合并、划分和删减。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种近距离微波成像方法，其特征在于，所述近距离微波成像方法包括：
将目标区域划分为预设分辨率大小的栅格集，所述栅格集中的栅格即为待成像的目标点，所述目标区域和所述栅格集的维度是一维、二维或者三维；
建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果；
控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止，所述预设个数的回波数据为由预设个数的栅格反射的回波信号经模数转换得到的数据；
控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
2.如权利要求1所述的近距离微波成像方法，其特征在于，所述对所述回波数据集进行成像处理，包括：
分别对所述回波数据集中的每个回波数据进行快速傅里叶变换；
根据所述快速傅里叶变换的结果，通过后向投影算法计算得到与所述每个回波数据对应的元素的值；
根据每个元素的值更新所述每个元素的初始值。
3.如权利要求2所述的近距离微波成像方法，其特征在于，所述分别对所述回波数据集中的每个回波数据进行快速傅里叶变换，包括：
分别对所述回波数据集中的每个回波数据沿频率轴方向进行快速傅里叶变换，得到与所述回波数据集对应的数据集，所述数据集包括与所述回波数据集中的每个回波数据对应的数据；
根据所述快速傅里叶变换的结果，通过后向投影算法计算得到与所述每个回波数据对应的元素的值，包括：
根据所述数据集，通过基于插值法的后向投影算法计算得到与所述每个回波数据对应的元素的值。
4.如权利要求3所述的近距离微波成像方法，其特征在于，插值法为样条插值法。
5.如权利要求2或3所述的近距离微波成像方法，其特征在于，所述快速傅里叶变换为并行快速傅里叶变换，所述后向投影算法为并行后向投影算法。
6.一种近距离微波成像系统，其特征在于，所述近距离微波成像系统包括：
栅格化模块，用于将目标区域划分为预设分辨率大小的栅格集，所述栅格集中的栅格即为待成像的目标点，所述目标区域和所述栅格集的维度是一维、二维或者三维；
分区存储模块，用于建立存储所述目标区域的成像结果的存储区域，所述成像结果以预设分辨率大小的矩阵形式表示，所述栅格集中的每个栅格均对应所述矩阵中的一个元素，元素的值即为其对应的栅格的成像结果；
数据采集模块，用于控制由预设个数的天线组成的线性阵列天线绕预设的圆弧轨道旋转扫描所述目标区域，在所述圆弧轨道上的每个方位角位置点，控制所述线性阵列天线采集预设个数的回波数据，并将所述预设个数的回波数据构成的回波数据集发送给信号处理设备，直到所述线性阵列天线在所述圆弧轨道上的预设方位角位置点完成回波数据的采集时止，所述预设个数的回波数据为由预设个数的栅格反射的回波信号经模数转换得到的数据；
数据处理模块，用于控制所述信号处理设备每接收到一次所述回波数据集，则对所述回波数据集进行成像处理，并将成像结果保存在所述矩阵中。
7.如权利要求6所述的近距离微波成像系统，其特征在于，所述数据处理模块包括：
傅里叶变换单元，用于分别对所述回波数据集中的每个回波数据进行快速傅里叶变换；
后向投影计算单元，用于根据所述快速傅里叶变换的结果，通过后向投影算法计算得到与所述每个回波数据对应的元素的值；
更新单元，用于根据每个元素的值更新所述每个元素的初始值。
8.如权利要求7所述的近距离微波成像系统，其特征在于，所述傅里叶变换单元具体用于：
分别对所述回波数据集中的每个回波数据沿频率轴方向进行快速傅里叶变换，得到与所述回波数据集对应的数据集；所述数据集包括与所述回波数据集中的每个回波数据对应的数据；
所述后向投影计算单元具体用于：
根据所述数据集，通过基于插值法的后向投影算法计算得到与所述每个回波数据对应的元素的值。
9.如权利要求8所述的近距离微波成像系统，其特征在于，插值法为样条插值法。
10.如权利要求7或8所述的近距离微波成像系统，其特征在于，所述快速傅里叶变换为并行快速傅里叶变换，所述后向投影算法为并行后向投影算法。
CN201610931503.XA 2016-10-31 2016-10-31 一种近距离微波成像方法及系统 Active CN106556874B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201610931503.XA CN106556874B (zh) 2016-10-31 2016-10-31 一种近距离微波成像方法及系统
US16/345,971 US11215730B2 (en) 2016-10-31 2017-08-04 Close range microwave imaging method and system
PCT/CN2017/096099 WO2018076884A1 (zh) 2016-10-31 2017-08-04 一种近距离微波成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610931503.XA CN106556874B (zh) 2016-10-31 2016-10-31 一种近距离微波成像方法及系统

Publications (2)

Publication Number Publication Date
CN106556874A CN106556874A (zh) 2017-04-05
CN106556874B true CN106556874B (zh) 2018-10-23

Family

ID=58443210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610931503.XA Active CN106556874B (zh) 2016-10-31 2016-10-31 一种近距离微波成像方法及系统

Country Status (3)

Country Link
US (1) US11215730B2 (zh)
CN (1) CN106556874B (zh)
WO (1) WO2018076884A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11561300B2 (en) 2017-04-28 2023-01-24 China Communication Technology Co., Ltd. Close-range microwave imaging method and system

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106556874B (zh) * 2016-10-31 2018-10-23 华讯方舟科技有限公司 一种近距离微波成像方法及系统
CN106680814B (zh) * 2016-12-09 2018-07-27 华讯方舟科技有限公司 旋转扫描三维成像系统
CN107167771B (zh) * 2017-04-28 2018-10-26 深圳市无牙太赫兹科技有限公司 一种微波成像系统的直达波抑制方法及系统
CN109597070B (zh) * 2018-12-29 2022-10-04 内蒙古工业大学 一种螺旋扫描式弧形阵列微波成像的方法及装置
CN109786977A (zh) * 2019-01-14 2019-05-21 河北华讯方舟太赫兹技术有限公司 一种天线面阵、安检装置以及安检方法
CN111812634B (zh) * 2020-06-05 2023-12-05 森思泰克河北科技有限公司 警戒线目标监测方法、装置和系统
CN112986986A (zh) * 2020-12-25 2021-06-18 中国航空工业集团公司沈阳飞机设计研究所 一种多发多收天线阵及其sar成像方法
CN112965061B (zh) * 2020-12-29 2022-06-21 北京理工大学 一种基于柱面mimo面阵的成像系统及成像方法
CN112835039B (zh) * 2020-12-30 2023-09-08 内蒙古工业大学 平面孔径分区域非线性渐进式相位迭代成像的方法及装置
CN112835040B (zh) * 2020-12-30 2023-05-23 内蒙古工业大学 球面孔径分区域渐进式相位迭代成像的方法及装置
CN112924965B (zh) * 2021-02-01 2022-04-12 桂林电子科技大学 基于聚类相参叠加的频率分集阵列雷达目标成像方法
CN115267766B (zh) * 2022-07-29 2024-07-26 中铁第六勘察设计院集团有限公司 利用地质雷达监测裂缝深部裂隙发育的装置
CN116977297B (zh) * 2023-07-14 2025-07-22 广西中烟工业有限责任公司 一种基于抑制因子的ffbp算法的烟箱缺条检测方法
CN117930231B (zh) * 2023-08-23 2024-08-06 苏州威陌电子信息科技有限公司 一种超近距离目标成像方法及装置
CN118859199A (zh) * 2024-06-26 2024-10-29 重庆邮电大学 一种沿轴旋转的毫米波mimo阵列扫描系统及成像方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5557283A (en) * 1991-08-30 1996-09-17 Sheen; David M. Real-time wideband holographic surveillance system
CN102520408A (zh) * 2011-12-30 2012-06-27 北京华航无线电测量研究所 一种圆柱阵面三维成像系统的三维成像方法
CN102540183A (zh) * 2010-12-14 2012-07-04 中国科学院电子学研究所 基于柱面几何的三维微波成像方法
CN102565795A (zh) * 2011-12-30 2012-07-11 北京华航无线电测量研究所 一种圆柱阵面三维成像系统的多角度成像显示方法
CN103018738A (zh) * 2011-09-20 2013-04-03 中国科学院电子学研究所 基于旋转天线阵列的微波三维成像方法
CN104133213A (zh) * 2014-07-23 2014-11-05 中国电子科技集团公司第四十一研究所 一种结合rm算法与bp算法的柱面近场三维rcs成像方法

Family Cites Families (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3713156A (en) * 1970-10-12 1973-01-23 R Pothier Surface and subsurface detection device
US5859609A (en) * 1991-08-30 1999-01-12 Battelle Memorial Institute Real-time wideband cylindrical holographic surveillance system
US7365672B2 (en) * 2001-03-16 2008-04-29 Battelle Memorial Institute Detection of a concealed object
US7405692B2 (en) * 2001-03-16 2008-07-29 Battelle Memorial Institute Detecting concealed objects at a checkpoint
US7248204B2 (en) * 2001-09-28 2007-07-24 Trex Enterprises Corp Security system with metal detection and mm-wave imaging
US6870162B1 (en) * 2003-01-31 2005-03-22 Millivision, Inc. Weighted noise compensating method and camera used in millimeter wave imaging
US7385549B2 (en) * 2003-08-12 2008-06-10 Trex Enterprises Corp Millimeter wave portal imaging system
US7212153B2 (en) * 2003-12-05 2007-05-01 Safeview, Inc. Millimeter-wave active imaging system with fixed array
US7205926B2 (en) * 2004-04-14 2007-04-17 Safeview, Inc. Multi-source surveillance system
US7253766B2 (en) * 2004-09-24 2007-08-07 Battelle Memorial Institute Three-dimensional surface/contour processing based on electromagnetic radiation interrogation
US6965340B1 (en) * 2004-11-24 2005-11-15 Agilent Technologies, Inc. System and method for security inspection using microwave imaging
US7844081B2 (en) * 2006-05-15 2010-11-30 Battelle Memorial Institute Imaging systems and methods for obtaining and using biometric information
US7804442B2 (en) * 2007-01-24 2010-09-28 Reveal Imaging, Llc Millimeter wave (MMW) screening portal systems, devices and methods
IL186884A (en) * 2007-10-24 2014-04-30 Elta Systems Ltd Object simulation system and method
DE102007055217B4 (de) * 2007-11-19 2013-01-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Verfahren und Vorrichtung zum ortsaufgelösten Erfassen und Rekonstruieren von Objekten mittels Mikrowellen
WO2009148637A1 (en) * 2008-02-14 2009-12-10 Brijot Imaging Systems, Inc. Millimeter wave concealed object detection system using portal deployment
DE102008011350A1 (de) * 2008-02-27 2009-09-03 Loeffler Technology Gmbh Vorrichtung und Verfahren zur Echtzeiterfassung von elektromagnetischer THz-Strahlung
US20100013920A1 (en) * 2008-07-21 2010-01-21 Safeview, Inc. Surveillance imaging with upsampling
CN104375141A (zh) * 2013-08-15 2015-02-25 同方威视技术股份有限公司 毫米波三维全息扫描成像设备及检查方法
CN104375143A (zh) * 2013-08-15 2015-02-25 同方威视技术股份有限公司 毫米波三维全息扫描成像设备及人体或物品检查方法
US10732276B2 (en) * 2013-10-21 2020-08-04 Sony Corporation Security system, method and device
DE102014210227A1 (de) * 2014-01-23 2015-07-23 Rohde & Schwarz Gmbh & Co. Kg System und Verfahren zur effizienten Abtastung von Gegenständen
AU2015227069B2 (en) * 2014-03-07 2020-05-14 Rapiscan Systems, Inc. Ultra wide band detectors
US9791553B2 (en) * 2014-04-03 2017-10-17 Evolv Technologies, Inc. Partitioning for radar systems
DK3129804T3 (da) * 2014-04-07 2019-08-12 Levitection Ltd Elektromagnetisk søgning og identifikation, i nærfeltsområder
AU2016244018B2 (en) * 2015-04-03 2020-10-08 Evolv Technologies, Inc. Modular imaging system
CN105699494B (zh) * 2015-12-28 2018-08-17 深圳市太赫兹科技创新研究院 毫米波全息三维成像检测系统和方法
CN105759315B (zh) * 2016-04-26 2018-10-23 华讯方舟科技有限公司 扫描机构及具有该扫描机构的安检仪
CN106556874B (zh) * 2016-10-31 2018-10-23 华讯方舟科技有限公司 一种近距离微波成像方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5557283A (en) * 1991-08-30 1996-09-17 Sheen; David M. Real-time wideband holographic surveillance system
CN102540183A (zh) * 2010-12-14 2012-07-04 中国科学院电子学研究所 基于柱面几何的三维微波成像方法
CN103018738A (zh) * 2011-09-20 2013-04-03 中国科学院电子学研究所 基于旋转天线阵列的微波三维成像方法
CN102520408A (zh) * 2011-12-30 2012-06-27 北京华航无线电测量研究所 一种圆柱阵面三维成像系统的三维成像方法
CN102565795A (zh) * 2011-12-30 2012-07-11 北京华航无线电测量研究所 一种圆柱阵面三维成像系统的多角度成像显示方法
CN104133213A (zh) * 2014-07-23 2014-11-05 中国电子科技集团公司第四十一研究所 一种结合rm算法与bp算法的柱面近场三维rcs成像方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11561300B2 (en) 2017-04-28 2023-01-24 China Communication Technology Co., Ltd. Close-range microwave imaging method and system

Also Published As

Publication number Publication date
WO2018076884A1 (zh) 2018-05-03
US11215730B2 (en) 2022-01-04
CN106556874A (zh) 2017-04-05
US20190243022A1 (en) 2019-08-08

Similar Documents

Publication Publication Date Title
CN106556874B (zh) 一种近距离微波成像方法及系统
CN107132510B (zh) 一种微波成像系统的幅相校正方法及系统
CN109143218B (zh) 一种基于vhf雷达的闪电定位系统及闪电定位方法
CN107102324B (zh) 一种近景微波成像方法及系统
JP6911861B2 (ja) 物体検知装置および物体検知方法
CN107167771B (zh) 一种微波成像系统的直达波抑制方法及系统
CN109813967A (zh) 一种阵列天线方向图的测量方法、设备、系统以及计算机可读存储介质
CN105068074B (zh) 一种窄带被动雷达三维成像方法
CN103616667A (zh) 一种用于散射成像的二维天线阵列布置方法
CN107390216A (zh) 基于波数域相干因子的高速超分辨率驻点扫描成像方法
CN106291476B (zh) 机载三维异构阵的雷达地面杂波回波获取方法
RU2507531C1 (ru) Способ пассивного определения координат источников излучения
CN110837127A (zh) 一种基于柱面雷达成像装置的稀疏天线布局方法
CN106646664B (zh) 基于gpu的人体微波回波模拟方法及系统
CN104698442A (zh) 机载下视阵列三维合成孔径雷达分布式三维场景仿真方法
CN107942326A (zh) 一种具有高通用性的二维主动毫米波成像方法
Nechaev et al. Evaluation and minimization of Cramer-Rao bound for conformal antenna arrays with directional emitters for DOA-estimation
CN115128607A (zh) 一种十字mimo阵列雷达系统及其三维成像方法
CN103969642B (zh) 一种用于多探头阵列成像的相位补偿方法
CN108427111A (zh) 一种雷达测距方法及装置
US11855338B2 (en) Radar device
CN108303685A (zh) 一种无源雷达超分辨三维成像方法及系统
CN109613533B (zh) 双站微波凝视关联成像方法及装置、存储介质及电子设备
CN110837126B (zh) 一种用于柱面阵列雷达成像的信号收发方法及装置
CN113176566B (zh) 一种用于mimo-sar的三维快速成像方法

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
TR01 Transfer of patent right

Effective date of registration: 20200915

Address after: 518000 Guangdong Province, Baoan District Xixiang street Shenzhen City Tian Yi Lu Chen Tian Bao Industrial District thirty-seventh building 430 room

Co-patentee after: SHENZHEN THZ SCIENCE AND TECHNOLOGY INNOVATION INSTITUTE

Patentee after: Shenzhen Zhongtou Huaxun Terahertz Technology Co.,Ltd.

Address before: 518000 404, building 37, chentian Industrial Zone, chentian community, Xixiang street, Bao'an District, Shenzhen City, Guangdong Province

Co-patentee before: SHENZHEN THZ SCIENCE AND TECHNOLOGY INNOVATION INSTITUTE

Patentee before: Shenzhen Huaxun ark Photoelectric Technology Co.,Ltd.

Effective date of registration: 20200915

Address after: 518000 404, building 37, chentian Industrial Zone, chentian community, Xixiang street, Bao'an District, Shenzhen City, Guangdong Province

Co-patentee after: SHENZHEN THZ SCIENCE AND TECHNOLOGY INNOVATION INSTITUTE

Patentee after: Shenzhen Huaxun ark Photoelectric Technology Co.,Ltd.

Address before: 518000 Guangdong city of Shenzhen province Baoan District Xixiang Bao Tian Yi Lu Chen Tian Industrial District thirty-seventh Building 1 floor and 2 floor on the West

Co-patentee before: SHENZHEN THZ SCIENCE AND TECHNOLOGY INNOVATION INSTITUTE

Patentee before: CHINA COMMUNICATION TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210927

Address after: 518000 Shenzhen, Baoan District, Xixiang, Guangdong Xixiang street, thirty-seventh fields, 430 estate.

Patentee after: Shenzhen Zhongtou Huaxun Terahertz Technology Co.,Ltd.

Address before: 518000 Shenzhen, Baoan District, Xixiang, Guangdong Xixiang street, thirty-seventh fields, 430 estate.

Patentee before: Shenzhen Zhongtou Huaxun Terahertz Technology Co.,Ltd.

Patentee before: SHENZHEN THZ SCIENCE AND TECHNOLOGY INNOVATION INSTITUTE

TR01 Transfer of patent right