基于STFT和Hough变换高频雷达机动目标检测.docVIP

基于STFT和Hough变换高频雷达机动目标检测 　　高频雷达回波中存在强大的地海杂波及各种高频干扰，一般通过空、时、频三维处理实现目标的检测。对于机动目标，该方法在两方面存在缺陷。其一，机动目标回波在相干积累时间内相位呈非线性变化，导致多普勒频谱展宽、检测性能下降。其二，上述处理过程不能获得目标的加速度。本文针对上述缺陷提出了基于STFT和Hough变换的高频雷达机动目标检测方法。其中，STFT用于获得目标回波的时频谱，而Hough变换通过对目标的时频“谱脊”进行非相干积累实现目标的检测并估计其速度和加速度。理论分析和数据处理结果均表明该方法在高频雷达机动目标检测方面性能良好。 　　【关键词】高频雷达 机动目标 STFT Hough变换 　　1 引言 　　高频雷达回波中存在强大的地海杂波及各种高频干扰，一般通过空、时、频三维处理实现目标的分辨。其中，距离上的分辨通过脉冲压缩来实现；方位上的分辨通过DBF来实现；速度上的分辨通过基于FFT的相干积累实现。上述处理方法在非机动目标检测方面性能良好。对于机动目标，该方法在两方面存在缺陷。其一，机动目标在相干积累时间内速度发生变化，导致其回波信号相位呈非线性变化，相干积累性能下降，表现为多普勒频谱展宽、SNR降低、系统的检测性能下降。设高频雷达以f0为载频发射连续线性调频信号，N为发射脉冲个数，T为脉冲重复周期，B为信号带宽，K=B/T为调频斜率。 　　图1为含机动目标的某高频雷达数据频谱及检测结果。其二，上述处理过程不能获得机动目标的加速度，而在某些情况下该参数对一些典型目标的识别具有重要意义。针对这两方面的问题，本文提出了基于STFT和Hough变换的高频雷达机动目标检测方法。其中，STFT用于获得目标回波的时频谱，而Hough变换通过对目标的时频“谱脊”进行非相干积累实现目标的检测并估计其速度和加速度。 　　2 基于STFT的回波时频谱 　　用傅利叶变换的方法提取信号频谱时，需要利用信号的全部时域信息。然而信号的局部变化会影响到信号的整个频谱。为了获得频谱关于时间定位的信息，可以用一个具有适当宽度的窗函数从信号中取出一段作傅利叶分析，于是得到信号在这段时间内的局部频谱。如果再让窗函数沿时间轴不断移动，便能对信号逐段进行频谱分析，从而获得目标回波的时频谱，这就是短时傅利叶变换的基本思想。 　　内的频谱信息。因此，用时域窗在以t=b为中心，宽为Dt的局部时间范围内考察信号得到的信息，也可以用频域窗在以为中心，宽为Dw的局部频率范围内考察该信号的频谱来得到。需要指出的是，由不定性原理可知，时域窗和频域窗不能同时任意窄。对于给定窗函数的STFT，其时间分辨率与频率分辨率之积是恒定的，它们不可能同时提高，只能以一种分辨率的降低换取另一种分辨率的提高。 　　本文所用的某雷达数据的发射脉冲个数为512，每次处理时窗口包含的脉冲周期数为128，窗口每次向前移动一个周期，所得到的处理结果如图2所示。 　　从图2可以看出，杂波的能量集中在零频附近，目标回波的时频谱为一条倾斜的直线，其斜率与目标的加速度对应，而直线与频率轴的交点对应着目标的初速度，即目标的时频 “谱脊”包含了其全部信息。在后续处理中，本文引入Hough变换来检测目标并估计其速度和加速度。 　　3 Hough变换平面上的目标检测 　　可见图像上的每一点都对应Hough参数空间中的一条正弦曲线，而处于同一条直线上的点所对应的一簇曲线必然交于一点。将参数空间划分为许多小区间，建立各区间对应的累加器，并将各累加器的初值置为0。对图像空间中的每个点均进行式（6）所示的变换，若所得变换结果落入Hough平面某区间中，则将该区间对应的累加器值加1，由此完成图像空间到参数空间的变换。 　　为了降低计算量和虚警率，在按本文所提方法进行Hough变换前，必须首先进行杂波抑制，去除零频附近的地海杂波，即使得杂波点不参与变换。抑制杂波后的时频谱见图3所示。 　　同时，为了进一步降低计算复杂度，在变换之前设定某一阈值，Hough变换只针对时频谱中高于该阈值的图像点进行。根据Hough变换原理，时频谱上机动目标对应 “谱脊”上的像素均落在同一累加器中，使得该累加器在Hough变换平面上表现为峰值，从而将时频谱上机动目标谱脊的检测转换为Hough平面上峰值点的检测。图4为数据的Hough变换结果。 　　在对图4所示的变换结果进行检测时，由于采样频率和相干积累时间对时频谱中机动目标 “谱脊”的倾斜角度有一定限制，故在Hough平面进行检测时极角维的检测范围不用覆盖。 　　需要注意的是，机动目标的 “谱脊”在频率维上存在一定展宽，因此，进行Hough变换时可能同一条 “谱脊”上某些边缘点会映射到参数空间中真实映射点附近的几个单元，即