[发明专利]一种基于稀疏频点的三维全息成像的重建方法

摘要

本发明提供一种基于稀疏频点的微波、毫米波、太赫兹三维全息成像的重建方法，属于图像处理技术领域。该方法通过对微波、毫米波、太赫兹雷达原始稀疏频点及其个数的选择，稀疏频点数据重建满频点数据，满频点数据处理得到的微波、毫米波、太赫兹三维全息成像结果。本发明方法能有效消除距离模糊现象，降低对雷达系统采样率和储存深度的要求，适合于实时应用成像中。

主权项

一种基于稀疏频点的微波、毫米波、太赫兹三维全息成像的重建方法，适用于微波、毫米波、太赫兹波的雷达发射信号，其特征在于，包括：S1：确定雷达发射信号进行采样的稀疏频点个数N，所述稀疏频点个数N满足：N≥max{(B+1+2B)r,2BfΔRc}---(1)]]>其中，B为设定中间量，且B＝lnM‑ln(1‑β)，M为待恢复的频点个数，即满频点的个数；β为最大可忍受的峰值旁瓣比r的置信系数；Bf为雷达发射信号的带宽，c为电磁波在自由空间的传播速度，ΔR为在雷达天线有效波束角的照射范围内的重建目标表面在距离向最大分布范围长度；S2：如果雷达发射信号属于线性调频连续波信号则先补偿掉该发射方式导致的剩余视频相位，消除剩余视频相位后得到和步进频连续波信号相同的回波信号后，再进行以下处理；如果雷达发射信号属于步进频连续波信号，则直接雷达发射信号进行以下处理：在带宽内选取均匀随机分布的N个频点f1,f2,…,fN，定义fsparse＝(f1,f2,…,fN)；采样得到N个稀疏频点的回波数据S1(x,y,ksparse)，建立坐标系XYZ；其中定义X、Y为扫描方向，Z为雷达照射方向，(x,y,0)为在坐标系XYZ中的天线相位中心位置，ksparse为fsparse对应的自由空间的波数；利用式(2)获得回波信号S1(x,y,ksparse)的转置S1T(x,y,ksparse)=[s1(x,y,k1),s1(x,y,k2),...,s1(x,y,kN)]---(2)]]>其中，m＝1,2…N，σ(x′,y′,z′)为在坐标系XYZ中位于(x′,y′,z′)处目标的散射系数，aθ(·)为雷达波束角的窗函数，S3：对S2所得与设定的参考信号利用式(4)求相关：Corr(Rj′)=S1T(x,y,ksparse)S2*(x,y,ksparse,Rj′)---(4)]]>其中，为的共轭矩阵；将Corr(R′j)结果中的最大值max[Corr(R′j)]对应的距离R′j作为在雷达天线有效波束角的照射范围内的目标表面在距离向估计值，记为距离R′max；S4：对S2所得回波数据进行频谱搬移，利用式(5)结合S3获得的距离R′max将其频谱搬移至基带，得到式中表示为矩阵的Hadamard积；S5：对S4所得结果利用式(6)对波数进行线性插值得到S5T(x,y,k)=interpksparse[S4T(x,y,ksparse)]---(6)]]>式中k＝(k′1,k′2,…,k′M)为满频点的各频率对应的自由空间的波数，其中k′i＝2π[fc‑Bf/2+(i‑1)Bf/(M‑1)]/c,i＝1,2,…,M，fc为雷达发射信号的中心频率，M为满频点的个数；interpX[Y(X)]为插值函数；S6：对S5所得结果进行频谱搬移，结合S3获得的距离R′max利用式(7)将其频谱的中心频率搬移至其实际所处的频率，得到式中S7：利用相位偏移算法对S6所得结果进行图像处理，得到恢复图像。