本发明提出了一种基于多视角成像的三维成像方法，所述方法包括如下步骤：步骤1：在同一笛卡尔坐标系下生成N个不同角度范围的三维图像；步骤2：生成与N个不同角度范围的三维图像对应的N个窗函数，并将该N个窗函数由圆柱坐标系转换到笛卡尔坐标系；步骤3：将所述N个不同角度范围的三维图像与所述N个窗函数分别对应相乘并求和，从而获得人体三维成像数据。本发明提出了一种实现360度全分辨三维成像的方法，在一幅三维图像中包含了360度角度的目标信息，有效的解决了边缘图像分辨率下降的问题，为后续的图像处理提供更丰富、准确和全面的信息。
1.一种基于多视角成像的三维成像方法，其特征在于，所述方法包括如下步骤：步骤1：在同一笛卡尔坐标系下生成N个不同角度范围的三维图像，具体包括：
1.1首先获取去调频后的三维回波数据，在圆柱坐标系下对三维回波数据做去除残余视频相位项的处理；
1.2将经过去除残余视频相位项处理后的数据分别在方位维和高度维做傅立叶变换，得到三维波数域数据；
1.3将所述三维波数域数据与匹配滤波函数相乘，获得匹配滤波后的三维波数域数据；
1.4对所述匹配滤波后的三维波数域数据在方位维做一维逆傅立叶变换后获得圆柱坐标系下的三维波数域数据；
1.5将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下；
1.6将在360度扫描范围的所述插值到笛卡尔坐标系下的三维波数域数据等分成N个成像区域；设n表示所述N个成像区域中的第n个成像区域，将选定的第n个成像区域以外置零；
1.7最后通过对N个所述第n个成像区域内的三维波数域数据做三维逆傅立叶变换，从而获得对应的N个所述第n个区域内的三维图像；
步骤2：生成与N个不同角度范围的三维图像对应的N个窗函数，并将该N个窗函数由圆柱坐标系转换到笛卡尔坐标系；
步骤3：将所述N个不同角度范围的三维图像与所述N个窗函数分别对应相乘并求和，从而获得人体三维成像结果。
2.如权利要求1所述的方法，其特征在于，所述步骤1.1中，去除残余视频相位项的方法为：将所述去调频后的三维回波数据在距离维做傅立叶变换后，乘以补偿项，然后做逆傅立叶变换。
3.如权利要求2所述的方法，其特征在于，所述步骤1.1中，所述补偿项如式(1)所示：其中，fi为频率，且fi∈[-fs/2,fs/2-fs/Nf]，fs为采样率，Nf为采样点数，γ为调频信号的调频斜率。
4.如权利要求1所述的方法，其特征在于，所述步骤1.3中，所述匹配滤波函数的具体表达式如下：其中， 为ε阶汉克尔函数，kr为r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，R0为天线中心到目标中心的参考距离，ε为扫描角度θ所对应的频率域。
5.如权利要求1所述的方法，其特征在于，所述步骤1.4中，对所述匹配滤波后的波数域数据在方位维做一维逆傅立叶变换后获得圆柱坐标系下的三维波数域数据，具体包括：所述圆柱坐标系下的三维波数域数据中任意一点的坐标可表示为(kr,θ,kz)，其中kr为r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
6.如权利要求1-5中任一项所述的方法，其特征在于，所述步骤1.5中，将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下，具体包括：选择正对目标中心的笛卡尔坐标系。
7.如权利要求1-5中任一项所述的方法，其特征在于，所述步骤1.5中，将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下，具体包括：在笛卡尔坐标系下的坐标为(kx,ky,kz)，其中kx＝kr*cos(θ)
其中kx为在笛卡尔坐标系下方位维x对应的波数域坐标，ky为在笛卡尔坐标系下距离维y对应的波数域坐标，kz为高度维z对应的波数域坐标，kr为在柱坐标系下r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
一种基于多视角成像的三维成像方法
技术领域
[0001] 本发明涉及雷达技术领域，特别涉及一种基于多视角成像的三维成像方法。
背景技术
[0002] 随着雷达成像技术的飞速发展，雷达探测模式从之前的一维测距逐渐发展到二维的测距测角、合成孔径雷达成像(SAR)、干涉SAR，以及三维都具备高分辨能力的平面孔径三维成像和圆柱孔径三维成像等模式。为了满足日益迫切的民用和军用需求，雷达的探测精度越来越高，成像水平已接近于光学水平。
[0003] 目前的圆柱扫描成像，逐渐成为人体安全检查的主流技术，多采取在空间域或频域分区域进行成像。该成像方法存在如下缺点：图像中间部分像素的分辨率较高，而边缘部分分辨率会下降，不能实现360度全分辨的三维成像，真实的反映出目标的信息，不利于后续图像处理和目标特征提取。
发明内容
[0004] 针对现有技术中的上述问题，本发明提供了一种微波圆柱扫描成像系统的人体三维成像方法。该方法通过对多角度三维回波数据的处理，可获得同时包含多角度信息的三维图像，从而为后续的图像处理提供更丰富更全面的三维信息。
[0005] 本发明的技术方案如下：
[0006] 一种基于多视角成像的三维成像方法，所述方法包括如下步骤：
[0007] 步骤1：在同一笛卡尔坐标系下生成N个不同角度范围的三维图像；
[0008] 步骤2：生成与N个不同角度范围的三维图像对应的N个窗函数，并将该N个窗函数由圆柱坐标系转换到笛卡尔坐标系；
[0009] 步骤3：将所述N个不同角度范围的三维图像与所述N个窗函数分别对应相乘并求和，从而获得人体三维成像结果。
[0010] 进一步的，所述步骤1：在同一笛卡尔坐标系下生成N个不同角度范围的三维图像；具体包括：
[0011] 1.1首先获取去调频后的三维回波数据，在圆柱坐标系下对三维回波数据做去除残余视频相位项(RVP)的处理；
[0012] 1.2将经过去除残余视频相位项处理后的数据分别在方位维和高度维做傅立叶变换，得到三维波数域数据；
[0013] 1.3将所述三维波数域数据与匹配滤波函数相乘，获得匹配滤波后的三维波数域数据；
[0014] 1.4对所述匹配滤波后的波数域数据在方位维做一维逆傅立叶变换后获得圆柱坐标系下的三维波数域数据；
[0015] 1.5将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下；
[0016] 1.6将在360度扫描范围的所述插值到笛卡尔坐标系下的三维波数域数据等分成N个成像区域；设n表示所述N个成像区域中的第n个成像区域，将选定的第n个成像区域以外置零；其中n的值为N个成像区域的顺序编号；
[0017] 1.7最后通过对N个所述第n个成像区域内的三维波数域数据做三维逆傅立叶变换，从而获得对应的N个所述第n个区域内的三维图像。
[0018] 进一步的，所述步骤1.1中，去除残余视频相位项的方法为：将所述去调频后的三维回波数据在距离维做傅立叶变换后，乘以补偿项，然后做逆傅立叶变换。
[0019] 进一步的，所述步骤1.1中，所述补偿项如式(1)所示：
[0020]
[0021] 其中，fi为频率，且fi∈[-fs/2,fs/2-fs/Nf]，fs为采样率，Nf为采样点数，γ为调频信号的调频斜率。
[0022] 进一步的，所述步骤1.3中，所述匹配滤波函数的具体表达式如下：
[0023]
[0024] 其中， 为ε阶汉克尔函数，kr为r方向所对应的波数域表示，且k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，R0为天线中心到目标中心的参考距离，ε为扫描角度θ所对应的频率域。
[0025] 进一步的，所述步骤1.4中，对所述匹配滤波后的波数域数据在方位维做一维逆傅立叶变换后获得圆柱坐标系下的三维波数域数据，具体包括：所述圆柱坐标系下的三维波数域数据中任意一点的坐标可表示为(kr,θ,kz)，其中kr为r方向所对应的波数域表示，且k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
[0026] 进一步的，所述步骤1.5中，将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下，具体包括：选择正对目标中心的笛卡尔坐标系。
[0027] 进一步的，所述步骤1.5中，将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下，具体包括：在笛卡尔坐标系下的坐标为(kx,ky,kz)，其中
[0028] kx＝kr*cos(θ)
[0029] ky＝kr*sin(θ)
[0030] kz＝kz (3)
[0031] 其中kx为在笛卡尔坐标系下方位维x对应的波数域坐标，ky为在笛卡尔坐标系下距离维y对应的波数域坐标，kz为高度维z对应的波数域坐标，kr为r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
[0032] 有益效果：
[0033] 本发明提出了一种实现360度全分辨三维成像的方法，在一幅三维图像中包含了360度角度的目标信息，有效的解决了边缘图像分辨率下降的问题，为后续的图像处理提供更丰富、准确和全面的信息。
附图说明
[0034] 图1微波人体三维成像方法流程图
[0035] 图2去除RVP项方法流程图
[0036] 图3波数域坐标变换关系图
[0037] 图4大波束宽度成像区域示意图
[0038] 图5小波束宽度成像区域示意图
具体实施方式
[0039] 下面结合附图和具体实施方式对本发明的技术方案做详细描述。
[0040] 本发明通过将去调频后的三维回波数据变换到频域，在频域将数据插值到同一笛卡尔坐标系下，并分区域完成补零、插值和加窗运算后，求和，完成多角度三维数据的融合处理，该三维融合数据中包含了多个视角的三维信息。
[0041] 图1为本发明的方法流程图。
[0042] 一种基于多视角成像的三维成像方法，所述方法包括如下步骤：
[0043] 步骤1：在同一笛卡尔坐标系下生成N个不同角度范围的三维图像；
[0044] 1.1首先获取去调频后的三维回波数据，在圆柱坐标系下对三维回波数据做去[0045] 除残余视频相位项(RVP)的处理。
[0046] 去除残余视频相位项处理的方法流程图如图2所示。去除残余视频相位项的方[0047] 法为：将所述去调频后的三维回波数据在距离维做傅立叶变换后，乘以补偿项，[0048] 然后做逆傅立叶变换。补偿项如式(1)所示：
[0049]
[0050] 其中，fi为频率，且fi∈[-fs/2,fs/2-fs/Nf]，fs为采样率，Nf为采样点[0051] 数，γ为调频信号的调频斜率。
[0052] 1.2将经过去除残余视频相位项处理后的数据在方位维和高度维分别做傅立叶[0053] 变换，得到三维波数域数据；
[0054] 1.3将所述三维波数域数据与匹配滤波函数相乘，获得匹配滤波后的三维波数域数据。所述匹配滤波函数可通过驻定相位原理推导获得，具体表达式如下：
[0055]
[0056] 其中， 为ε阶汉克尔函数，kr为r方向所对应的波数域表示，且k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，R0为天线中心到目标中心的参考距离，ε为扫描角度θ所对应的频率域。
[0057] 1.4对所述匹配滤波后的波数域数据在方位维做一维逆傅立叶变换后获得圆柱坐标系下的三维波数域数据，所述圆柱坐标系下的三维波数域数据中任意一点的坐标可表示为(kr,θ,kz)，其中kr为r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
[0058] 1.5将所述圆柱坐标系下的三维波数域数据插值到笛卡尔坐标系下(通常选择正对目标中心的笛卡尔坐标系)，三维波数域数据坐标变换关系图如图3所示。坐标变换关系如式(3)所示。在笛卡尔坐标系下的坐标为(kx,ky,kz)，其中
[0059]
[0060] 其中kx为在笛卡尔坐标系下方位维x对应的波数域坐标，ky为在笛卡尔坐标系下距离维y对应的波数域坐标，kz为高度维z对应的波数域坐标，kr为r方向所对应的波数域表示，且 k为波数，且k＝2*π*f/c，kz为z方向所对应的波数域表示，f为系统的工作频率，c为光速，θ为波数域数据的扫描角度。
[0061] 1.6将在360度扫描范围的所述插值到笛卡尔坐标系下的三维波数域数据等分成N个成像区域；设n表示所述N个成像区域中的第n个成像区域，将选定的第n个成像区域以外置零。其中n的值为N个成像区域的顺序编号，即，当对N个成像区域中的第1个成像区域进行处理时，n＝1；当对N个成像区域中的第2个成像区域进行处理时，n＝2；依此类推，当对N个成像区域中的第N个成像区域进行处理时，n＝N。
[0062] 所述N个成像区域的划分可根据天线单元的波束宽度、目标成像区域的范围、天线到目标中心的距离R0、圆柱成像系统本身的结构特点以及微波不能穿透人体的特性来设定。具体如下：
[0063] 1.6.1首先确定成像区域中心半径r上的点获得最优分辨率的时候所对应的圆心角φ。
[0064] 下面分两种情况来分析。
[0065] 第一种情况：如果天线的波束宽度θantenna如图4所示满足式(4)：
[0066]
[0067] 其中半径r为目标上任意一点到目标中心的距离，距离R0为天线到目标中心的距离，且孔径的长度大于B点到C点的长度BC，即成像区域的圆周角φ满足式(5)：
[0068]
[0069] 则A点获得最优的分辨率。由式(4)和式(5)可知，r越小，获得最优分辨率所需要的天线波束宽度越小，所需成像区域的圆周角越大。因此需要根据r的最小值rmin确定成像区域的范围，即成像区域的圆周角 由于微波不能穿透人体，rmin一般选为0.1～0.2m。
[0070] 第二种情况：如果天线的波束宽度θantenna如图5所示，满足式(6)：
[0071]
[0072] 如图5所示，根据余弦定理得到式(7)所示的关系：
[0073]
[0074] ，其中A′C′为A′点到C′点的距离。
[0075] 由式(7)推导得式(8)
[0076]
[0077] 如图5所示，根据正弦定理得到式(7)所示的关系：
[0078]
[0079] 由式(8)和式(9)可知，半径r越小，A′点得到最优分辨率，所需的成像区域越大，所以选择成像区域时，以r较小的值计算最优分辨率所需的成像区域，由于微波不能穿透人体，rmin一般选为0.1～0.2m。
[0080] 综合第一种情况和第二种情况，选择成像区域时，以r最小处的点目标获得最优分辨率所需的成像区域为标准。其中rmin一般选为0.1～0.2m。
[0081] 1.6.2设需要获得最优分辨率的区域为φ0，则保留的数据区域范围为φ0+φ，即区域的起始点前移φ/2，区域的结束点后移φ/2。这样才能确保φ0内所有的点满足最优分辨率。
[0082] 1.6.3确定需要分解的成像区域个数N如式(9)所示
[0083]
[0084] 则N个成像区域的范围分别为：(nφ0-φ/2):(nφ0+φ/2)。
[0085] 1.6.4将选定的第n个成像区域之外置零。
[0086] 1.6.5对n从1到N进行遍历，从而对应得到N个所述第n个成像区域内的三维波数域数据。
[0087] 1.7最后通过对所述N个第n个成像区域内的三维波数域数据做三维逆傅立叶变换，从而获得对应的N个第n个区域内的三维图像。
[0088] 步骤2:生成N个窗函数，并将该N个窗函数由圆柱坐标系转换到笛卡尔坐标系；
[0089] 根据天线的波束宽度以及人体成像区域的范围，确定第n个成像区域能够全分辨成像的区域。将能够全分辨成像的区域，窗函数的幅值设为1，两边区域的窗函数幅值则逐渐衰减。将窗函数从圆柱坐标系(r,θ,z)变换到笛卡尔坐标系(x,y,z)下。坐标变换关系如式(10)所示：
[0090]
[0091] 步骤3：将N个角度的三维图像与N个窗函数分别对应相乘并求和，从而获得人体三维成像数据。
[0092] 根据上述具体实施方式介绍可见，本发明提出了一种融合多视角信息的三维成像方法，通过N个角度三维成像、加窗、求和等方法，实现三维数据的融合处理，为后续图像处理提供更丰富更全面的三维信息。
[0093] 上述具体实施方式仅用于解释和说明本发明的技术方案，并不能构成对权利要求保护范围的限定。本领域技术人员应当清楚，在上述技术方案的基础上进行任何简单的修改和替换而得到的新的技术方案，均将落入本发明的保护范围之内。
