CN104569971A - 一种近距离三维全息成像方法及系统 - Google Patents
一种近距离三维全息成像方法及系统 Download PDFInfo
- Publication number
- CN104569971A CN104569971A CN201410827394.8A CN201410827394A CN104569971A CN 104569971 A CN104569971 A CN 104569971A CN 201410827394 A CN201410827394 A CN 201410827394A CN 104569971 A CN104569971 A CN 104569971A
- Authority
- CN
- China
- Prior art keywords
- msub
- signals
- mrow
- electromagnetic wave
- signal
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000033001 locomotion Effects 0.000 claims abstract description 19
- 238000005070 sampling Methods 0.000 claims description 29
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 11
- 230000007704 transition Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 abstract 2
- 230000000694 effects Effects 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- 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/35—Details of non-pulse systems
- G01S7/352—Receivers
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
本发明涉及一种近距离三维全息成像方法及系统,包括如下步骤:发射连续电磁波雷达信号;在以时间、圆周角及Z轴方向形成三维域中测得回波信号;对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像。本发明一种近距离三维全息成像方法及系统,通过在以时间、圆周角及Z轴方向形成三维域中测得回波信号,根据重构的目标散射强度信号进行三维全息成像,在连续波信号成像过程中未进行运动补偿,实现了对目标物体较好的三维全息成像。
Description
技术领域
本发明涉及一种成像方法及系统,尤其涉及一种近距离三维全息成像方法及系统。
背景技术
调频连续波信号与不同波长信号成像技术的融合,促进了一个宽带、有效、低耗、高质量的成像系统形成,特别是在安全探测系统的应用中,当天线阵列连续的发射和接收调频连续波信号,其不间断的运动影响将不再会被忽略,因此,传统的在合成孔径成像算法中的间断式方法需要在调频连续波成像处理中被优化改进,传统算法,例如波数域算法、频率缩放算法以及范围多普勒算法等,均是集中在调频连续波孔径成像数据的优化上。现有的图像成像处理方法注重成像数据的优化,并没有考虑到信号传送过程中运动的影响,因此,严重影响电磁波信号探测的成像效果。
发明内容
本发明解决的技术问题是:构建一种近距离三维全息成像方法及系统,克服现有技术连续波信号成像过程中未考虑运动的影响,成像效果不佳的技术问题。
本发明的技术方案是:提供一种近距离三维全息成像方法,包括如下步骤:
发射雷达信号:沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号;
获取采样信号:在以时间、圆周角及Z轴方向形成三维域中测得回波信号;
信号转换:利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;
重构回波信号:对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像。
本发明的进一步技术方案是:沿待成像物表面环绕发射连续电磁波雷达信号。
本发明的进一步技术方案是:在相对方分别环绕发射连续电磁波雷达信号。
本发明的进一步技术方案是:所述双线性插值运算包括对于三维空间波数域内的回波信号进行非均匀采样和均匀采样的插值运算。
本发明的进一步技术方案是:对于空间波数域内的非均匀采样,还包括在空间波数域内进行向均匀采样过度的差值运算。
本发明的技术方案是:构建一种近距离三维全息成像系统,包括电磁波发射源、获取采样信号的信号采样模块、进行信号转换的信号转换模块、重构回波信号的重构模块、成像模块,所述电磁波发射源沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号,所述信号采样模块在以时间、圆周角及Z轴方向形成三维域中测得回波信号;所述信号转换模块利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;所述重构模块对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,所述成像模块根据重构的目标散射强度信号进行三维全息成像。
本发明的进一步技术方案是:所述电磁波发射源为多个,多个所述电磁波发射源排列成阵列。
本发明的进一步技术方案是:所述电磁波发射源分别沿待成像物表面环绕发射连续电磁波雷达信号。
本发明的进一步技术方案是:所述电磁波发射源为至少两个,所述电磁波发射源在相对方分别环绕发射连续电磁波雷达信号。
本发明的进一步技术方案是:还包括所述重构模块对于三维空间波数域内的回波信号进行非均匀采样和均匀采样的插值运算。
本发明的技术效果是:构建一种近距离三维全息成像方法及系统,包括如下步骤:沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号;在以时间、圆周角及Z轴方向形成三维域中测得回波信号;利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像。本发明一种近距离三维全息成像方法及系统,通过在以时间、圆周角及Z轴方向形成三维域中测得回波信号,对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像,在连续波信号成像过程中未进行运动补偿,实现了对目标物体较好的三维全息成像。
附图说明
图1为本发明的成像系统模型。
图2为本发明的成像系统结构模块图。
具体实施方式
下面结合具体实施例,对本发明技术方案进一步说明。
如图1、图2所示,本发明的具体实施方式是:提供一种近距离三维全息成像方法,包括如下步骤:定义被成像物体区域为(X0,Y0,Z0)=(R0cosθ,R0sinθ,Z)的圆柱体,其中R0为需要成像区域半径,θ是圆柱坐标系统中的角度,θ∈[0,2π],天线阵列长度即沿着Z轴方向的合成孔径长度为LZ,孔径中心位置z=ZC的平面。成像过程中,天线阵列绕着被成像物体或者部分的绕着被成像物体转动,形成圆周θ方向的合成孔径。采样位置为(R,θ,Z),物体任意成像位置Pn的坐标为(xn,yn,zn),其对应的散射强度为σ(xn,yn,zn)。
电磁波发射源1沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号。沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号。具体实施例中,将待成像物体视为柱状,沿待成像物表面发射连续电磁波雷达信号时,可以将雷达信号发射头沿待成像物环绕转动一周,环绕过程中,同时发射连续电磁波雷达信号。也可以将雷达信号发射头环绕转动一定弧度即可,只要其发射的雷达信号覆盖待成像物体即可。
具体实施过程如下:发射雷达信号为p(t),其中f0是基本频率,t是单个信号发射周期内的时间变量,K是对发射信号的频率进行扫描的速率,我们假设雷达信号的发射时间是τ,接收时间为τ+τd,其中τd是双向延迟时间,天线元件与靶目标之间的瞬时距离范围是R(τ)~R(τ+τd)之间。双向延迟时间可以表示为
其中,c是光速,
其中,
τ=nTθ+mTzv+t=τn+τm+t, (3)
n是阵列元素的个数,m是沿着正面图方向的采样个数,Tθ是沿着弧形阵列元件在方位域中的信号发射周期,
获取采样信号:信号采样模块2在以时间、圆周角及Z轴方向形成三维域中测得回波信号。若天线阵列与目标之间的距离非常短,(1)式可以近似表示为:
雷达的辐射图在聚束目标区上是不变的,在(t,θ,z)域中测得单点靶目标Pn点的回波信号为
SR(t,θ,z)=σ(xn,yn,zn)·p(t-τd) (5)
信号转换模块3利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换。
具体实施过程如下:在去调频接收的合成孔径系统中,为了减少采样要求以及数据传输速率,对接收到的信号可用一个参考信号进行最大化,我们假设参考信号是一个延时时间为τi的发射信号,则回波信号可表示为:
SF(t,θ,z)=σ(xn,yn,zn)·exp[-j2πf0(τd-τi)]exp[-j2πK(τd-τi)(t-τi)] (6)
f=K(t-τi),将其代入(6)式中,可以得到
SF(f,θ,z)=σ(xn,yn,zn)exp[-j2πK(f+f0)(τd-τc)] (7)
由(3)式代入(7)式可得
其中
带着空间变量zm(zm=vτm=vmTy),对(8)式进行一维傅里叶变换即可得
用固定相位法对(9)式进行估计,偏导后取极值,同时结合式(8),取 最终得
当被检测体体积较大时,我们定义t=f/K+2Ri/c,
信号模型可表示为
其中,
重构模块4对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,成像模块5根据重构的目标散射强度信号进行三维全息成像。
基于Parseval定理以及利用圆形对称函数的傅里叶性质推导后可得出:
推导可知:
其中:
对于式(14),kzz0表示阵列的原始位置,kzvτn表示由于第n个天线元件的运动引起的仰角范围变换所对应的相位变化,表示在一个扫描时间内由于阵列的运动引起的空间非变量值,和2krRi表示方位角和距离的常数变化变量。
对于空间波数域内的非均匀采样,需要在空间波数域内进行向均匀采样过度的差值运算,圆柱采样数Jkxy(kxy,θ,kz)通过一个插值算法转换成J(kx,ky,kz),其中kx=kxycosθ,ky=kxysinθ,σ(x,y,z)和J(kx,ky,kz)组成了一个傅里叶变换对,进行三维反傅里叶变换最终得到直角坐标系下重构的目标散射强度为
σ(x,y,z)=v∫∫∫J(kx,ky,kz)exp[j(kxx-kyy+kzz)]dkxdkydkz (15)
得到直角坐标系下重构的目标散射强度信号,所述成像模块根据重构的目标散射强度信号进行三维全息成像。
本发明的具体实施方式是:构建一种近距离三维全息成像系统,包括电磁波发射源1、获取采样信号的信号采样模块2、进行信号转换的信号转换模块3、重构回波信号的重构模块4、成像模块5,所述电磁波发射源1沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号,所述信号采样模块2在以时间、圆周角及Z轴方向形成三维域中测得回波信号;所述信号转换模块3利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;所述重构模块4对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,所述成像模块5根据重构的目标散射强度信号进行三维全息成像。
本发明的技术效果是:构建一种近距离三维全息成像方法及系统,包括如下步骤:沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号;在以时间、圆周角及Z轴方向形成三维域中测得回波信号;利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像。本发明一种近距离三维全息成像方法及系统,通过在以时间、圆周角及Z轴方向形成三维域中测得回波信号,对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像,在连续波信号成像过程中未进行运动补偿,实现了对目标物体较好的三维全息成像。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1. 一种近距离三维全息成像方法,包括如下步骤:
发射雷达信号:沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号;
获取采样信号: 在以时间、圆周角及Z轴方向形成三维域中测得回波信号;
信号转换:利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;
重构回波信号并成像:对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,根据重构的目标散射强度信号进行三维全息成像。
2.根据权利要求1所述近距离三维全息成像方法,其特征在于,沿待成像物表面环绕发射连续电磁波雷达信号。
3.根据权利要求1所述近距离三维全息成像方法,其特征在于,在相对方分别环绕发射连续电磁波雷达信号。
4.根据权利要求1所述近距离三维全息成像方法,其特征在于,所述双线性插值运算包括对于三维空间波数域内的回波信号进行非均匀采样和均匀采样的插值运算。
5.根据权利要求4所述近距离三维全息成像方法,其特征在于,对于空间波数域内的非均匀采样,还包括在空间波数域内进行向均匀采样过度的差值运算。
6.一种近距离三维全息成像系统,其特征在于,包括电磁波发射源、获取采样信号的信号采样模块、进行信号转换的信号转换模块、重构回波信号的重构模块、成像模块,所述电磁波发射源沿待成像物表面发射连续电磁波雷达信号,所述连续电磁波雷达信号包括连续波雷达探测信号和连续波雷达参考信号,所述信号采样模块在以时间、圆周角及Z轴方向形成三维域中测得回波信号;所述信号转换模块利用参考信号对接收到的回波信号进行最大化,对最大化的回波信号进行傅里叶变换,再利用相位固定法,实现回波信号的时域向频域的转换;所述重构模块对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿,得到直角坐标系下重构的目标散射强度信号,所述成像模块根据重构的目标散射强度信号进行三维全息成像。
7.根据权利要求6所述近距离三维全息成像系统,其特征在于,所述电磁波发射源为多个,多个所述电磁波发射源排列成阵列。
8.根据权利要求6所述近距离三维全息成像系统,其特征在于,所述电磁波发射源沿待成像物表面环绕发射连续电磁波雷达信号。
9.根据权利要求6所述近距离三维全息成像系统,其特征在于,所述电磁波发射源为至少两个,所述电磁波发射源在相对方分别环绕发射连续电磁波雷达信号。
10.根据权利要求6所述近距离三维全息成像系统,其特征在于,还包括所述重构模块对于三维空间波数域内的回波信号进行非均匀采样和均匀采样的插值运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827394.8A CN104569971A (zh) | 2014-12-25 | 2014-12-25 | 一种近距离三维全息成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827394.8A CN104569971A (zh) | 2014-12-25 | 2014-12-25 | 一种近距离三维全息成像方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104569971A true CN104569971A (zh) | 2015-04-29 |
Family
ID=53086489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410827394.8A Pending CN104569971A (zh) | 2014-12-25 | 2014-12-25 | 一种近距离三维全息成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104569971A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563958A (zh) * | 2017-08-23 | 2018-01-09 | 广州视源电子科技股份有限公司 | 全息图像转换方法和系统 |
CN112114310A (zh) * | 2020-08-27 | 2020-12-22 | 博微太赫兹信息科技有限公司 | 一种基于三维分解的微波毫米波全息图像重建方法 |
-
2014
- 2014-12-25 CN CN201410827394.8A patent/CN104569971A/zh active Pending
Non-Patent Citations (1)
Title |
---|
B.L.REN ET AL.: "Modified Cylindrical Holographic Algorithm For Three-Dimensional Millimeter-Wave Imaging", 《PROGRESS IN ELECTROMAGNETICS RESEARCH》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563958A (zh) * | 2017-08-23 | 2018-01-09 | 广州视源电子科技股份有限公司 | 全息图像转换方法和系统 |
CN112114310A (zh) * | 2020-08-27 | 2020-12-22 | 博微太赫兹信息科技有限公司 | 一种基于三维分解的微波毫米波全息图像重建方法 |
CN112114310B (zh) * | 2020-08-27 | 2022-08-02 | 博微太赫兹信息科技有限公司 | 一种基于三维分解的微波毫米波全息图像重建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105044719B (zh) | 一种基于圆周sar的太赫兹高精度垂直曲面成像方法 | |
Wang et al. | Nonuniform frequency diverse array for range-angle imaging of targets | |
CN108107431B (zh) | 一种圆柱扫描sar三维成像快速实现方法 | |
CN104515989A (zh) | 一种近距离毫米波三维全息成像方法及系统 | |
Wei et al. | Sparse array microwave 3-D imaging: Compressed sensing recovery and experimental study | |
Wang et al. | Bistatic synthetic aperture radar imaging using ultranarrowband continuous waveforms | |
CN102540183B (zh) | 基于柱面几何的三维微波成像方法 | |
CN112162326B (zh) | 全息成像安检系统以及安检方法 | |
Li et al. | Efficient near-field imaging using cylindrical MIMO arrays | |
Yanik et al. | Millimeter-wave near-field imaging with two-dimensional SAR data | |
Zheng et al. | ISAR imaging for fluctuating ships based on a fast bilinear parameter estimation algorithm | |
Vouras et al. | An overview of advances in signal processing techniques for classical and quantum wideband synthetic apertures | |
CN106855619B (zh) | 一种获取mimo成像雷达系统各方向的分辨率的方法 | |
Setsu et al. | Super-Resolution Doppler Velocity Estimation by Kernel-Based Range–$\tau $ Point Conversions for UWB Short-Range Radars | |
KR101103416B1 (ko) | 역합성 개구면 레이더 영상의 거리횡단 방향을 스케일링하는 방법 및 그 장치 | |
CN104569971A (zh) | 一种近距离三维全息成像方法及系统 | |
Akiyama et al. | Low complexity algorithm for range-point migration-based human body imaging for multistatic UWB radars | |
Li et al. | A high-resolution handheld millimeter-wave imaging system with phase error estimation and compensation | |
CN104569999A (zh) | 一种近距离红外三维全息成像方法及系统 | |
CN104765034A (zh) | 一种近距离太赫兹三维全息成像方法及系统 | |
Hu et al. | A Fast Wavenumber Domain 3D Near-Field Imaging Algorithm for Cross MIMO Array | |
Zhao et al. | 0.14 THz imaging system for security and surveillance | |
JPWO2019180767A1 (ja) | 物体検知装置、物体検知方法、及びプログラム | |
Yarman et al. | Passive synthetic aperture radar imaging with single frequency sources of opportunity | |
Yamaryo et al. | Range-point migration-based image expansion method exploiting fully polarimetric data for UWB short-range radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150429 |
|
WD01 | Invention patent application deemed withdrawn after publication |