非均匀的数据变换到频域，从而 在频域对回波数据进行处理。首先对雷达接收的回波信号进行去调频处理，再进行傅立叶 变换实现距离向压缩。通过特显点法估计出雷达在方位向采样的位置，然后利用基于最小 二乘法的非均匀傅立叶变换将方位向数据变换到频域。最后在频域通过stolt插值以及二 维傅立叶反变换得到最终的成像结果。与传统的补偿方法相比，本方法在成像过程中有效 减小了由于雷达运动不均匀造成的方位向散焦现象。
[0042] 本发明采用基于非均匀傅立叶变换的成像方法对方位向数据进行非均匀傅立叶 变换，将时域非均匀分布的数据变换到频域均匀分布的数据。这种方法直接在成像过程中 进行运动误差的补偿，避免了多普勒调频率的估计和成像后的相位校正。
[0043] 综上所述，本发明方法与传统方法相比，本方法直接在成像过程中进行运动误差 的补偿，避免了多普勒调频率的估计和成像后的相位校正，即减小了由于雷达运动不理想 造成的方位向散焦现象。
【附图说明】
[0044] 图1为本发明基于非均匀傅立叶变换的太赫兹SAR成像方法流程图；
[0045] 图2为方位向直接傅立叶变换得到的点目标成像结果图；
[0046] 图3为方位向采用非均匀傅立叶变换得到的点目标成像结果图。
【具体实施方式】
[0047] 下面结合附图对本发明做进一步的说明。
[0048]采用前述方法，流程具体如图1所示，并设定仿真实验中心频率为320GHz，带宽为 28.8GHz，雷达平台在轨道上做直线运动，运动速度均值为lm/s，同时叠加有正弦形式的速 度误差，目标设置为一个点目标。雷达平台运动的误差导致方位向的数据采集不均匀。对采 集得到的数据进行成像，距离向做距离压缩，方位向基于波数域变换进行聚焦处理。图2为 方位向直接傅立叶变换得到的成像结果图，从图中可以看到，成像结果是目标在方位向存 在严重的散焦现象，不能识别为一个点目标。图3为本方法的成像结果图。从图中可以看到， 目标点在方位向聚焦效果良好，成像结果表现为一个点。
[0049]对比图2和图3,可见本发明方法直接在成像过程中进行运动误差的补偿，避免了 多普勒调频率的估计和成像后的相位校正。
【主权项】
1. 一种非匀速太赫兹雷达平台SAR成像方法，包括以下步骤： 步骤1:雷达发射太赫兹频段0.1 THz~ΙΟΤΗζ线性调频信号其中，τ为快时间，t为慢时间，TP为脉冲宽度，f。为中心频率，选取在太赫兹频段，γ为调点目标坐标为(xP，yP，0)，在慢时间时刻t，雷达的位置为(xa(t)，0，h)，其中方位向坐标 ' (〇 = W + .外)，v为理想状态下的雷达速度，玲；)为雷达速度不均匀导致的位置偏离；则t 时刻，点目标与雷达的瞬时距离R(t)为雷达接收到的回波信号为以雷达运动轨迹与场景中心的垂直距离作为参考距离Rref，做去调频处理，输出信号为其中AR=R( t) -Rref，式中相位项是去调频处理产生的剩余视频相位，在距离 压缩前将其去除；对去除剩余视频相位后的信号做距离向傅立叶变换，实现距离向的压缩，得到在距离 频域上各目标的sine状脉冲信号； 步骤2:通过回波信号计算雷达在方位向各个采样时刻的位置 对步骤1中距离向傅立叶变换后的数据进行截取，获得特显点目标对应的数据矩阵，对 所得数据矩阵在一个方位向上搜索，得到该方位向上的最大值以及最大值对应的位置，对 每个方位向重复以上搜索得到所有方位向上的最大值和其对应的位置； 由上述位置信息得到点目标与雷达的瞬时距离取其最小值为点目标与雷达的最 近距离i?min =min{郎；)}，由公式计算出每个方位向采样时刻(to，11，…，tN )对应的方位向位置，进而得到对应的采样序 列号k步骤3:基于步骤2得到的采样序列号，对步骤1中得到的信号数据x(c)做方位向，即慢 时间域的非均匀傅立叶变换，得到方位向频域均匀分布的信号，〇 = -Ν/2，···，Ν/2-1表示方 位向序列号； 步骤3-1:构造矩阵其中，Ν为方位向数据长度，ω =e#/mN，m为信号过采样率，q为插值点数；，.其中{mk} =mk-[mk]，[mk]是mk的整数部分，n =由方程组Fzn(k) =an(k)，求得在不同采样位置k处对应的q+Ι个插值系数zn(k); 〇, ··· ,q； s。为定标因子，取余弦定标因子 步骤3-2:对方位向数据x(c)计算傅立叶变换系数約， 再利用快速傅立叶变换方法计算通过定标因子进行数据变换，得到方位向频域数据通过以上变换将回波信号的方位向变换到频域，即得到回波信号在方位向的频域形 式； 步骤4:信号做STOLT插值处理并进行二维傅立叶反变换得到SAR图像； 将步骤3得到的数据在对应的波数域分布为(Kx，KR)域，其中，Kx表示方位向波数，K R表示 目标与雷达的视线波数，二者并非垂直关系，stolt变换公式为通过距离向sine插值即可得到变换结果S(Kx，Ky); 步骤5:将步骤4中的结果S(Kx，Ky)做二维逆傅立叶变换，即得到最终的成像结果。
【专利摘要】本发明属于雷达成像技术，涉及一种非匀速太赫兹雷达平台SAR成像方法。首先对接收到的线性调频信号做去调频处理，然后进行距离向傅立叶变换，得到距离向压缩的数据。通过特显点的回波信号估计出雷达由于非匀速运动造成的误差，从而得到方位向数据的采样位置，然后根据估计得到的方位向采样位置，利用非均匀傅立叶变换，将回波数据在方位向变换到波数域，再经过STOLT变换以及傅立叶反变换，得到最终的成像结果。该成像方法将方位向采样不均匀的数据准确的变换到波数域，在波数域进行插值操作，完成方位向聚焦。由于没有进行距离近似，所以得到的成像结果是最优的，适用于高分辨率的太赫兹SAR成像。
【IPC分类】G01S7/41, G01S13/90
【公开号】CN105676223
【申请号】CN201610049975
【发明人】闵锐, 郝居博, 李晋, 曹宗杰, 皮亦鸣
【申请人】电子科技大学
【公开日】2016年6月15日
【申请日】2016年1月25日