用于太赫兹高斯波束全息成像的三维大景深图像重建方法[发明专利]

(12)发明专利申请

(10)申请公布号(10)申请公布号 CN 104732584 A (43)申请公布日(43)申请公布日 2015.06.24

(21)申请号 201510119167.4(22)申请日 2015.03.18

(71)申请人中国科学院电子学研究所

地址100080 北京市海淀区北四环西路19

号(72)发明人李超 刘玮 张群英 方广有(74)专利代理机构北京理工大学专利中心

11120

代理人张瑜 仇蕾安(51)Int.Cl.

G06T 17/00(2006.01)G01S 13/(2006.01)

权利要求书2页 说明书7页 附图2页

()发明名称

用于太赫兹高斯波束全息成像的三维大景深图像重建方法(57)摘要

本发明提供一种用于太赫兹高斯波束全息成像的三维大景深图像重建方法，属于图像处理技术领域。该方法通过对太赫兹雷达原始回波信号的三维频谱进行参考距离匹配滤波，距离方位向的解耦合，方位向的匹配滤波得到的太赫兹三维全息成像结果。本发明方法能精确恢复大景深分布的目标图像，处理过程简单快速，适用于大带宽实时应用成像中。

C N 1 0 4 7 3 2 5 8 4 A CN 104732584 A

权 利 要 求 书

1/2页

1.一种用于太赫兹高斯波束全息成像的三维大景深图像重建方法，其特征在于，包括：

S1：如果雷达发射信号属于线性调频连续波信号则先补偿掉该发射方式导致的剩余视频相位，消除剩余视频相位后得到和步进频连续波信号相同的的回波信号后，再进行以下处理；如果雷达发射信号属于步进频连续波信号，则直接雷达发射信号进行以下处理：

对接收到的回波信号作方位向二维傅里叶变换得到三维频谱S1(kx,ky,k)，并建立坐标系XYZ，定义X、Y为扫描方向，Z为雷达照射方向；kx、ky为方位向波数，k为自由空间的波数，三维频谱S1(kx,ky,k)表示为：

式中σ(x′,y′,z′)为坐标系XYZ中位于(x′,y′,z′)处目标的散射系数，w0为雷达天线发射波束的束腰大小，

S2：对S1所得的三维频谱S1(kx,ky,k)乘以参考距离匹配滤波函数H(kx,ky,zo,k)，获得三维频谱S2(kx,ky,k)，其中参考距离匹配滤波函数H(kx,ky,zo,k)表示为：

z0为目标所处位置的参考距离；

S3：对S2所得三维频谱S2(kx,ky,k)进行距离方位向解耦合，得到解耦合频谱结果即三维频谱S7(kx,ky,z)，其中，z为目标在坐标系XYZ中Z方向的取值；

具体的：S31：根据太赫兹雷达系统参数和目标最大分布距离比率γ选取参数α，且参数α选择需满足：

式中B为发射信号带宽，Δf为回波数据的频率采样间隔，目标最大分布距离比率γ的取值范围为γ∈(0,1)，kc为发射信号中心频率波数；c为电磁波在自由空间的传播速度；

S32：对步骤S2所得的三维频谱S2(kx,ky,k)乘以线性调频信号H1(kx,ky,kb,α)并对乘积中的kb进行逆傅里叶变换，获得三维频谱S3(kx,ky,z)，线性调频信号H1(kx,ky,kb,α)表达式为：

2

CN 104732584 A

权 利 要 求 书

2/2页

其中kb为基频波数；

S33：对三维频谱S3(kx,ky,z)乘以线性调频信号H2(kx,ky,z,α)并对乘积中的z进行傅里叶变换，获得三维频谱S4(kx,ky,kb)，线性调频信号H2(kx,ky,z,α)表达式为：

S34：根据太赫兹雷达系统参数和选取的参数γ、α得到补零个数N，并对三维频谱S4(kx,ky,kb)的基频波数域进行N点补零得到S5(kx,ky,kb)，且N满足：

式中c为电磁波在自由空间的传播速度，ceil(a)表示为大于a的最小整数；

S35：对三维频谱S5(kx,ky,kb)乘以线性调频信号H3(kx,ky,kb,α)并对乘积中的kb进行逆傅里叶变换，获得三维频谱S6(kx,ky,z)，线性调频信号H3(kx,ky,kb,α)表达式为：

S37：对三维频谱S6(kx,ky,z)乘以线性调频信号H4(kx,ky,z,α)，获得三维频谱

线性调频信号H4(kx,ky,z,α)表达式为：S7(kx,ky,z)，

S4：对S3所得三维频谱S7(kx,ky,z)进行方位向匹配滤波，乘以方位向匹配滤波函数

H5(kx,ky,z′-zo)，得到三维频谱S8(kx,ky,z)，方位向匹配滤波函数H5(kx,ky,z′-zo)表达式为：

步骤S5：对步骤S4所得的三维频谱S8(kx,ky,z)作方位向二维逆傅里叶变换，得到三维频谱S9(x,y,z)，根据三维频谱S9(x,y,z)进行图像处理得到三维聚焦图像；S9(x,y,z)表达式为：

式中

3

CN 104732584 A

说 明 书

1/7页

用于太赫兹高斯波束全息成像的三维大景深图像重建方法

技术领域

本发明属于图像处理技术领域，尤其涉及一种用于太赫兹高斯波束全息成像的三维大景深图像重建方法。

[0001]

背景技术

近年来，随着太赫兹辐射源、检测器及其他相关元器件的逐渐成熟，太赫兹全息成

像技术越来越受到人们的关注。太赫兹全息成像技术利用太赫兹波可以穿透很多非极性材料(如皮革、塑料、泡沫等)的特性对人体或其他物体进行高分辨率穿透成像。太赫兹波与X射线相比不会因为光致电离而破坏被检测的物质，对人体来说太赫兹全息成像是基本无害，具有较好的安全性，因此比较适用于公共区域的安全检测领域。[0003] 随着人们对公共安全检测的要求越来越高，对大景深目标区域进行高分辨率快速成像的应用越来越广泛。为了实现同时对一定距离范围内的所有目标的高分辨率成像，需要能够处理太赫兹雷达系统接收到的包含有目标信息的大带宽，高频点数的回波数据的太赫兹三维全息成像方法。

[0004] 目前有两种方式，一种是太赫兹三维相位偏移算法(Phase-Shift Migration,EPSM)，能够对大带宽、大景深数据进行太赫兹三维全息图像重建，但是耗时长，不能满足实时性要求。另一种是太赫兹三维快速增强相位偏移算法(Enhanced Phase-Shift Migration,EPSM)，该方法通过对相位偏移因子进行合理的近似来实现利用快速傅里叶变换和近似的相位偏移因子替代非常耗时的相位偏移累加操作，提高了运算效率，能够实时显示精确聚焦的太赫兹全息图像。但是由于在三维全息重建过程中解决了非常耗时的操作，而且其方法流程中只利用了快速傅里叶变换和复数乘法操作，太赫兹三维快速增强相位偏移算法处理大带宽，高频点数的回波数据时，由于相位偏移因子近似的合理性不再满足，造成重建的太赫兹三维全息图像出现散焦，并且随着带宽和频点数的增加和景深的扩大，散焦现象越来越严重，会对三维全息图像的图像重建造成很大的影响，经常无法重建得到清晰的目标图像，不再适用于大景深目标区域实时应用成像中。因此太赫兹三维快速增强相位偏移算法能够在一定带宽和景深范围内快速而精确地进行太赫兹三维全息图像重建。

[0002]

发明内容

为解决上述问题，本发明提供一种用于太赫兹高斯波束全息成像的三维大景深图

像重建方法，不仅能够对大带宽、大景深数据进行太赫兹三维全息图像重建，而且满足实时性要求。

[0006] 本发明的用于太赫兹高斯波束全息成像的三维大景深图像重建方法，其括：

[0005]

S1：如果雷达发射信号属于线性调频连续波信号则先补偿掉该发射方式导致的剩余视频相位，消除剩余视频相位后得到和步进频连续波信号相同的的回波信号后，再进行以下处理；如果雷达发射信号属于步进频连续波信号，则直接雷达发射信号进行以下处

[0007]

4

CN 104732584 A

说 明 书

2/7页

理：

对接收到的回波信号作方位向二维傅里叶变换得到三维频谱S1(kx,ky,k)，并建立坐标系XYZ，定义X、Y为扫描方向，Z为雷达照射方向；kx、ky为方位向波数，k为自由空间的波数，三维频谱S1(kx,ky,k)表示为：

[0008]

[0009]

[0010] 式中σ(x′,y′,z′)为坐标系XYZ中位于(x′,y′,z′)处目标的散射系

数，w0为雷达天线发射波束的束腰大小，

[0011]

S2：对S1所得的三维频谱S1(kx,ky,k)乘以参考距离匹配滤波函数

获得三维频谱S2(kx,ky,k)，其中参考距离匹配滤波函数H(kx,ky,zo,k)表示

为：

[0012]

z0为目标所处位置的参考距离；[0014] S3：对S2所得三维频谱S2(kx,ky,k)进行距离方位向解耦合，得到解耦合频谱结果即三维频谱S7(kx,ky,z)，其中，z为目标在坐标系XYZ中Z方向的取值；[0015] 具体的：[0016] S31：根据太赫兹雷达系统参数和目标最大分布距离比率γ选取参数α，且参数α选择需满足：

[0013]

[0017]

式中B为发射信号带宽，Δf为回波数据的频率采样间隔，目标最大分布距离比率

γ的取值范围为γ∈(0,1)，kc为发射信号中心频率波数；c为电磁波在自由空间的传播速度；

[0019] S32：对步骤S2所得的三维频谱S2(kx,ky,k)乘以线性调频信号H1(kx,ky,kb,α)并对乘积中的kb进行逆傅里叶变换，获得三维频谱S3(kx,ky,z)，线性调频信号H1(kx,ky,kb,α)表达式为：

[0018] [0020] [0021]

其中kb为基频波数；

5

CN 104732584 A[0022]

说 明 书

3/7页

S33：对三维频谱S3(kx,ky,z)乘以线性调频信号H2(kx,ky,z,α)并对乘积中的z进行傅里叶变换，获得三维频谱S4(kx,ky,kb)，线性调频信号H2(kx,ky,z,α)表达式为：

[0023]

S34：根据太赫兹雷达系统参数和选取的参数γ、α得到补零个数N，并对三维频谱S4(kx,ky,kb)的基频波数域进行N点补零得到S5(kx,ky,kb)，且N满足：

[0024] [0025]

式中c为电磁波在自由空间的传播速度，ceil(a)表示为大于a的最小整数；[0027] S35：对三维频谱S5(kx,ky,kb)乘以线性调频信号H3(kx,ky,kb,α)并对乘积中的kb进行逆傅里叶变换，获得三维频谱S6(kx,ky,z)，线性调频信号H3(kx,ky,kb,α)表达式为：

[0026] [0028]

S37：对三维频谱S6(kx,ky,z)乘以线性调频信号H4(kx,ky,z,α)，获得三维频谱S7(kx,ky,z)，线性调频信号H4(kx,ky,z,α)表达式为：

[0029] [0030]

S4：对S3所得三维频谱S7(kx,ky,z)进行方位向匹配滤波，乘以方位向匹配滤波函数H5(kx,ky,z′-zo)，得到三维频谱S8(kx,ky,z)，方位向匹配滤波函数H5(kx,ky,z′-zo)表达式为：

[0031] [0032] [0033]

步骤S5：对步骤S4所得的三维频谱S8(kx,ky,z)作方位向二维逆傅里叶变换，得到

三维频谱S9(x,y,z)，根据三维频谱S9(x,y,z)进行图像处理得到三维聚焦图像；S9(x,y,z)表达式为：

[0034]

[0035] [0036]

式中

有益效果：

[0037] 本发明通过变标操作和快速傅里叶变换来代替耗时的操作累加，不涉及满足或不满足相位偏移因子近似的合理性问题，从而能够满足实时性问题也不会造成图像散焦的问题。具体的：利用快速傅里叶变换和实数相乘来实现参考距离匹配滤波，然后通过距离方位向的解耦合和方位向的匹配滤波得到的太赫兹三维全息成像结果；从而在不降低运算效率的情况下，可以有效的处理大带宽，高频点数的太赫兹雷达三维数据，对大景深目标区域进

6

CN 104732584 A

说 明 书

4/7页

行高分辨率快速成像。附图说明

[0038] [0039] [0040] [0041] [0042]

图1为本发明的使用的雷达系统的扫描示意图；

图2为本发明的待成像目标的光学图像；

图3为现有技术的太赫兹三维快速增强相位偏移算法处理结果示意图；图4为本发明的本方法处理结果示意图；图5为两种方法处理结果对比。

具体实施方式

[0043] 本发明的目的是提出一种用于太赫兹高斯波束全息成像的三维快速大景深图像重建方法，该方法通过对太赫兹雷达原始回波信号的三维频谱进行参考距离匹配滤波，距离方位向的解耦合，方位向的匹配滤波得到的太赫兹三维全息成像结果。本发明方法能精确恢复大景深分布的目标图像，处理过程简单快速，适用于大带宽实时应用成像中。[0044] 为了达到上述目的，本发明的技术解决方案是：[0045] 步骤S1：如果雷达发射信号为线性调频连续波信号，则需要先补偿掉该发射方式导致的剩余视频相位，消除剩余视频相位后得到和步进频连续波信号相同的的回波信号再进行如下处理；如果雷达发射信号为步进频连续波信号，直接进行如下处理：对接收到的回波信号作方位向二维傅里叶变换，得到三维频谱S1(kx,ky,k)，如图1所示，定义x、y为扫描方向，z为雷达照射方向，kx、ky为方位向波数，k为自由空间的波数，三维频谱S1(kx,ky,k)可以表示为：

[0046]

[0047]

式中σ(x′,y′,z′)为位于(x′,y′,z′)处目标的散射系数，w0为雷达天

线发射波束的束腰大小，

步骤S2：消除回波信号三维频谱S1(kx,ky,k)的参考距离处的距离徙动。对步骤S1

所得的结果S1(kx,ky,k)乘以参考距离匹配滤波函数H(kx,ky,zo,k)，消除了三维频谱中的距离不变项的相位，获得频谱S2(kx,ky,k)，其中参考距离匹配滤波函数H(kx,ky,zo,k)可以表示为：

[0048] [0049]

其中z0为参考距离。

[0051] 步骤S3：对步骤S2所得频谱S2(kx,ky,k)进行距离方位向解耦合，得到解耦合频谱

[0050]

7

CN 104732584 A

说 明 书

5/7页

结果S7(kx,ky,z)。具体包括子步骤：[0052] 子步骤S3(a)：为了避免子步骤S3(b)对频谱S2(kx,ky,k)造成过大的频谱扩展，利用奈奎斯特采样定理，可根据太赫兹雷达系统参数和目标最大分布距离比率γ选取参数α。参数α选择需满足：

[0053]

式中B为发射信号带宽，Δf为回波数据的频率采样间隔，目标最大分布距离比率

γ的取值范围为γ∈(0,1)。[0055] 子步骤S3(b)：对步骤S2所得的结果S2(kx,ky,k)乘以线性调频信号H1(kx,ky,kb,α)对不同距离处回波信号进行置斜处理，并对置斜结果进行关于kb的逆傅里叶变换，获得频谱S3(kx,ky,z)，其中kb为基频波数。线性调频信号H1(kx,ky,kb,α)表达式为：

[00] [0056]

式中kc为发射信号中心频率波数

[0058] 子步骤S3(c)：补偿三维频谱的差分距离单元徙动。对子步骤S3(b)所得的结果S3(kx,ky,z)乘以线性调频信号H2(kx,ky,z,α)并对z进行傅里叶变换，消除三维频谱的距离方位向解耦合，获得频谱S4(kx,ky,kb)。线性调频信号H2(kx,ky,z,α)表达式为：

[0057] [0059]

子步骤S3(d)：消除子步骤S3(b)造成的频谱扩展造成的混叠，根据太赫兹雷达系统参数和选取的参数γ、α得到补零个数N，并对子步骤S3(c)所得的结果S4(kx,ky,kb)的基频波数域进行N点补零得到S5(kx,ky,kb)。根据奈奎斯特采样定理可得，其中补零个数N满足：

[0060] [0061]

式中c为电磁波在自由空间的传播速度，ceil(a)表示为大于a的最小整数。

[0063] 子步骤S3(e)：为了补偿子步骤S3(b)引入的置斜操作，对子步骤S3(d)所得的结果S5(kx,ky,kb)乘以线性调频信号H3(kx,ky,kb,α)并对kb进行逆傅里叶变换，获得频谱S6(kx,ky,z)。线性调频信号H3(kx,ky,kb,α)表达式为：

[0062] [00]

子步骤S3(f)：补偿子步骤S3(c)引入的剩余视频相位，对子步骤S3(e)所得的

结果S6(kx,ky,z)乘以线性调频信号H4(kx,ky,z,α)，获得频谱S7(kx,ky,z)。线性调频信号H4(kx,ky,z,α)表达式为：

[0065] [0066]

8

CN 104732584 A[0067]

说 明 书

6/7页

步骤S4：补偿方位向剩余相位。对步骤S3所得结果S7(kx,ky,z)进行方位向匹配滤波，乘以方位向匹配滤波函数H5(kx,ky,z′-zo)，得到S8(kx,ky,z)。方位向匹配滤波函数H5(kx,ky,z′-zo)表达式为：

[0068] [0069]

步骤S5：对步骤S4所得的结果S8(kx,ky,z)作方位向二维逆傅里叶变换，得到三维

聚焦图像S9(x,y,z)，其表达式为：

[0070]

[0071] [0072]

式中

此外，上述对各元件和方法的定义并不仅限于实施例中提到的各种具体结构、形状或方式，本领域普通技术人员可对其进行简单地更改或替换，例如：[0073] 对于目标与雷达天线距离比较近的情况下可以设置z0＝0，因此可以省略步骤S2。

[0074] 下面通过实测数据对该算法进行了验证。高斯波束照射下太赫兹逐点扫描三维全息成像雷达系统参数如表1所示：[0075] 表1

[0076]

9

CN 104732584 A[0077]

说 明 书

7/7页

图2为待成像目标的光学图像。图3为太赫兹三维快速增强相位偏移算法处理结果，图4为本方法通过参考距离匹配滤波，距离方位向的解耦合，方位向的匹配滤波得到处理结果。可以看出本方法能对目标进行很好的聚焦成像，而太赫兹三维快速增强相位偏移算法处理结果有一定程度的散焦。为了进一步比较两种方法的成像效果，如图5所示，为对重建结果中目标左侧七个金属条的Y方向切片进行显示的结果，可以看出本方法可以清楚地分辨七个金属条，而太赫兹三维快速增强相位偏移算法的重建结果的分辨率变差，对距离比较近的第五六七根金属条不能很好的分辨，而前四根金属条随能够分辨，但分辨率已经出现明显恶化。[0078] 当然，本发明还可有其他多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

10

CN 104732584 A

说 明 书 附 图

1/2页

图1

图2

11

图3

CN 104732584 A

说 明 书 附 图

2/2页

图4

图5

12