一种基于毫米波测云雷达的自适应频域检测方法技术

本申请公开一种基于毫米波测云雷达的自适应频域检测方法，首先进行数据预处理和噪声功率自适应估计；其次依据检测噪声门限对功率谱进行逐频点检测，并对保护窗外的噪声点进行抑制；再次完成功率谱谱矩计算，并进行SNR估计和检测，得到检测后的逻辑矩阵；最后，对中值滤波后的检测结果进行评估，完成检测门限的自适应修正迭代。本申请所记载的技术方案减小了噪声对计算的影响，计算简便，实现简便，利于在实际工程中移植，可以应用于毫米波测云雷达等气象探测设备中。

【技术实现步骤摘要】
一种基于毫米波测云雷达的自适应频域检测方法
本专利技术涉及雷达
，特别涉及一种基于毫米波测云雷达的自适应频域检测方法。
技术介绍
检测是指确认雷达测量值是否存在目标，即测量值仅为干扰，或者测量值为干扰与目标回波之和。检测的本质是对样本空间进行划分，完成气象目标检测是毫米波测云雷达的一项重要数据处理功能。检测时针对的样本数据信噪比越大，检测效果越明显。毫米波测云雷达整个信号处理过程中，能够提升信噪比的步骤包括相干积累、周期图法功率谱估计和非相干积累。因此非相干积累后的距离-多普勒数据是信号处理过程中信噪比最大的样本数据，本文阐述的检测方法主要针对该阶段数据进行处理。检测中的关键步骤包括噪声电平功率估计和检测门限确定。其主要估计方法包括固定噪声功率法、远区距离库法、客观噪声电平法、最大速度区法以及分段噪声电平计算方法。固定噪声功率法在早期的天气雷达中通过设定一个固定值或谱峰下的某固定值(Battan,1964；Donaldson,1967；Sekhonetal,1971；Donaldsonetal,1972),但由于雷达接收到的功率谱并不是固定不变的，在不同天气背景下差异较大，该方法效果并不理想，而且不同设备间噪声功率各不相同，因此该方法局限性较大。风廓线雷达和天气雷达中常采用远距离库法，在垂直探测时认为雷达的远距离库没有气象信号，功率谱中仅仅包含噪声信号，因此，把远距离库的噪声平均值作为一根径向上所有的噪声电平(王莎等，2012)。毫米波雷达功率谱噪声在不同高度的其分布亦，即使在相同距离单元，不同时间或者温度湿度环境条件稍有变化噪声功率也有所差别。1974年，Hildebrand和Sekhon基于高斯白噪声特点，提出了客观噪声电平计算法，该方法主要建立在两个假设上一是噪声幅度在频带范围内满足均匀分布，二是噪声幅度的概率分布满足高斯分布，因此该方法具有一定的局限性，而且计算复杂，计算量较大实时性差。在风廓线雷达功率谱研究中，Monique提出一种分段的噪声电平计算方法，该方法多用于风廓线雷达的功率谱处理中，在毫米波雷达中应用需要配合后续质量控制，否则检测后杂波点过多。另外，刘黎平等在国产的行波管毫米波雷达中，将功率谱大速度大值区的功率作为噪声电平，提出了最大速度法，但该方法估计的噪声电平波动较大，影响了探测性能，在速度模糊的情况下容易产生误判。有鉴于此，提出一种基于毫米波测云雷达的自适应频域检测方法。
技术实现思路
本专利技术的目的在于提供一种基于毫米波测云雷达的自适应频域检测方法。为达到上述目的，本专利技术采用下述技术方案，包括：S1.对数据进行预处理和噪声功率自适应估计；S2.依据检测噪声门限对功率谱进行逐频点检测，并对保护窗外的噪声点进行抑制；S3.再次完成功率谱谱矩计算进行SNR估计和检测，对得到检测后的逻辑矩阵中值滤波；S4.对中值滤波后的检测结果进行评估，完成检测门限的自适应修正迭代。进一步地，S1所述对数据进行预处理包括:对功率谱数据进行转化，所述功率谱数据为对数值PdB(fk)，单位为dB，得到所测功率的线性值P(fk)；优选的，转化公式为：P(fk)＝10^[PdB(fk)/10]，fk＝(k-FFTNum/2)/FFTNum×Prf/CNum，假设FFT点数为FFTNum，则第k条谱线fk采用上述公式计算，其中k的范围为0～FFTNum-1，其中Prf表示脉冲重复频率，CNum表示相干积累点数。进一步地，S1所述噪声功率自适应估计包括：将距离单元内的频谱等分，计算每一份的平均功率得到其集合，找出其中最小的一个作为噪声功率估计NoiseBase。进一步地，S2所述逐频点检测包括：S201a.计算确定检测噪声门限NoiseThreshold，NoiseThreshold＝NoiseDeep*NoiseBase，NoiseBase为S1中计算得到的噪声功率，NoiseDeep为上一帧计算得到的噪声深度；S201b.对功率谱进行逐频点检测，将频谱信号强度小于噪声门限值的频点检测掉，保留高于噪声门限的频点。进一步地，S2所述对保护窗外的噪声点进行抑制包括：S202a.找到功率谱最大的谱线，以其为中点设置保护窗，保护窗外的孤立噪点不受保护，即使该谱线超过门限也可将其去除；S202b.将叠加在目标信号谱线上的噪声去除；S202c.估计得到该距离单元内的噪声功率；S202d.得到检测后的功率谱信号。进一步地，S3所述再次完成功率谱谱矩计算包括：求解零阶矩，以得到平均功率估计公式为：其中Y(fk)表示第k根谱线功率谱，FFTNum表示FFT点数，对其取对数可以得到对数值。进一步地，S3所述进行SNR估计和检测包括：S301.获取噪声功率和信号平均功率对其做差求得信噪比SNR，其表达式为：S302.设置输出信噪比门限SNRThreshold，然后对输出信噪比进行检测，如果目标SNR大于SNRThreshold，则认为目标存在，输出所估计距离库的功率谱谱矩参数，对检测目标逻辑矩阵Detection赋值为1；如果目标SNR小于SNRThreshold，则认为目标不存在，所计算的功率谱谱矩参数为无效值，对检测目标逻辑矩阵Detection赋值为0。进一步地，S3所述中值滤波包括：S301a.假设雷达回波的距离单元个数为M，射线数为N，在经过前述步骤后得到二维Detection[M,N]逻辑矩阵；S301b.以二维窗对Detection[M,N]进行滑窗中值滤波处理；优选的，其处理方法为：以3×3的二维窗进行滑窗中值滤波处理，二维窗内共9个点，如果其中有5个点及以上为逻辑1，则表示中心位置目标存在，反之则目标不存在，经过处理后可以得到矩阵Detection_Med[M,N]。进一步地，S4所述检测结果进行评估包括：对Detection[M,N]中的0和1的个数进行统计，其结果假设分别为A0个和B0个，其中A0+B0＝M*N。对Detection_Med[M,N]中的0和1的个数进行统计，其结果假设分别为A1个和B1个，其中A1+B1＝M*N。其中虚警概率可以表示为：Pf＝(B0-B1)/A1。进一步地，S4所述检测门限的自适应修正迭代包括：预设期望虚警概率为Pf0，若满足以下条件：(Pf-Pf0)>|ε|，则对噪声深度NoiseDeep进行如下调整：NoiseDeep＝NoiseDeep+NoiseDeep*0.01；若满足以下条件：(Pf-Pf0)<-|ε|，则对噪声深度NoiseDeep进行如下调整：NoiseDeep＝NoiseDeep-NoiseDeep*0.01。附图说明下面结合附图对本专利技术的具体实施方式作进一步详细的说明：图1为本专利技术的处理流程图；图2为噪声基底估计示意图；图3为逐频点门限检测示意图；图4为保护窗外噪声抑制示意图；图5为距离-多普勒域数据示例；图6为最终检测结果；图7为各阶矩气象参数；图8为反射率因子估计结果。具体实施方式为了更清楚地说明本专利技术，下面结合优选实施例和附图对本专利技术做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本专利技术的保护范围。北京XX研究所生产的垂直本文档来自技高网...

【技术保护点】
一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，包括：S1.对数据进行预处理和噪声功率自适应估计；S2.依据检测噪声门限对功率谱进行逐频点检测，并对保护窗外的噪声点进行抑制；S3.再次完成功率谱谱矩计算进行SNR估计和检测，对得到检测后的逻辑矩阵中值滤波；S4.对中值滤波后的检测结果进行评估，完成检测门限的自适应修正迭代。

【技术特征摘要】
1.一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，包括：S1.对数据进行预处理和噪声功率自适应估计；S2.依据检测噪声门限对功率谱进行逐频点检测，并对保护窗外的噪声点进行抑制；S3.再次完成功率谱谱矩计算进行SNR估计和检测，对得到检测后的逻辑矩阵中值滤波；S4.对中值滤波后的检测结果进行评估，完成检测门限的自适应修正迭代。2.根据权利要求1所述的一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，S1所述对数据进行预处理包括:对功率谱数据进行转化，所述功率谱数据为对数值PdB(fk)，单位为dB，得到所测功率的线性值P(fk)；优选的，转化公式为：P(fk)＝10^[PdB(fk)/10]，fk＝(k-FFTNum/2)/FFTNum×Prf/CNum，假设FFT点数为FFTNum，则第k条谱线fk采用上述公式计算，其中k的范围为0～FFTNum-1，其中Prf表示脉冲重复频率，CNum表示相干积累点数。3.根据权利要求1所述的一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，S1所述噪声功率自适应估计包括：将距离单元内的频谱等分，计算每一份的平均功率得到其集合，找出其中最小的一个作为噪声功率估计NoiseBase。4.根据权利要求1所述的一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，S2所述逐频点检测包括：S201a.计算确定检测噪声门限NoiseThreshold，NoiseThreshold＝NoiseDeep*NoiseBase，NoiseBase为S1中计算得到的噪声功率，NoiseDeep为上一帧计算得到的噪声深度；S201b.对功率谱进行逐频点检测，将频谱信号强度小于噪声门限值的频点检测掉，保留高于噪声门限的频点。5.根据权利要求1所述的一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，S2所述对保护窗外的噪声点进行抑制包括：S202a.找到功率谱最大的谱线，以其为中点设置保护窗，保护窗外的孤立噪点不受保护，即使该谱线超过门限也可将其去除；S202b.将叠加在目标信号谱线上的噪声去除；S202c.估计得到该距离单元内的噪声功率；S202d.得到检测后的功率谱信号。6.根据权利要求1所述的一种基于毫米波测云雷达的自适应频域检测方法，其特征在于，S3所述再次完成功率谱谱矩计算包括：求解零阶矩，以得到平均功率估计公式为：其...

【专利技术属性】
技术研发人员：张鹏，王志锐，周亭亭，
申请(专利权)人：北京无线电测量研究所，
类型：发明
国别省市：北京,11