CN109270536B - 一种多普勒天气雷达的信号处理器 - Google Patents
一种多普勒天气雷达的信号处理器 Download PDFInfo
- Publication number
- CN109270536B CN109270536B CN201811287767.1A CN201811287767A CN109270536B CN 109270536 B CN109270536 B CN 109270536B CN 201811287767 A CN201811287767 A CN 201811287767A CN 109270536 B CN109270536 B CN 109270536B
- Authority
- CN
- China
- Prior art keywords
- signal
- gain
- channel
- horizontal
- vertical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001914 filtration Methods 0.000 claims abstract description 73
- 238000012545 processing Methods 0.000 claims abstract description 71
- 238000003908 quality control method Methods 0.000 claims abstract description 60
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 238000011161 development Methods 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 37
- 238000012937 correction Methods 0.000 claims description 28
- 230000006978 adaptation Effects 0.000 claims description 22
- 238000012360 testing method Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000002310 reflectometry Methods 0.000 claims description 9
- 239000013256 coordination polymer Substances 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 claims description 8
- 230000021615 conjugation Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 239000013307 optical fiber Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 101800000863 Galanin message-associated peptide Proteins 0.000 claims 1
- 102100028501 Galanin peptides Human genes 0.000 claims 1
- 230000006870 function Effects 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 3
- NUHSROFQTUXZQQ-UHFFFAOYSA-N isopentenyl diphosphate Chemical compound CC(=C)CCO[P@](O)(=O)OP(O)(O)=O NUHSROFQTUXZQQ-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多普勒天气雷达的信号处理器,包括信号处理单元;所述信号处理单元包括:信号接收模块、信号存储模块、信号分析模块、滤波模块、信号处理模块、质量控制模块、参数配置模块。本发明解决了信号处理算法的研制开发和仿真调试的难度大,以及信号处理算法的改进和功能升级困难的问题,采用通用服务器作为信号处理单元的硬件平台,运算速度快、扩展性强,采用软件化设计,通过软件的重构和不断改进,以提升或扩展雷达系统功能,实现快速多元化探测。
Description
技术领域
本发明涉及雷达信号处理器的开发领域,尤其是一种多普勒天气雷达的信号处理器。
背景技术
现有技术中,雷达的信号处理部分一般采用专用DSP芯片结合FPGA实现,由于FPGA、DSP芯片编程的复杂度较高,使得信号处理算法的研制开发、仿真调试的难度增大,特别是受限于芯片内部资源,很难在后期对信号处理算法进行改进和功能升级;且雷达信号处理部分采用背板的专用差分总线或并行数据总线进行接收IQ信号,数据传输的速率受限于背板总线速率。
发明内容
为了克服上述现有技术的缺陷,本发明提供一种多普勒天气雷达的信号处理器,解决了信号处理算法的研制开发和仿真调试的难度大,以及信号处理算法的改进和功能升级困难的问题,采用通用服务器作为信号处理单元的硬件平台,运算速度快、扩展性强,采用软件化设计,通过软件的重构和不断改进,以提升或扩展雷达系统功能,实现快速多元化探测。
为实现上述目的,本发明采用以下技术方案,包括:
一种多普勒天气雷达的信号处理器,包括信号处理单元,所述信号处理单元包括以下组成部分:
信号接收模块、信号存储模块、信号分析模块、滤波模块、信号处理模块、质量控制模块、参数配置模块;其中,
所述信号接收模块接收来源于雷达接收机的IQ信号,并将IQ信号分别发送至信号存储模块和信号分析模块;所述IQ信号为一种正交的基带信号;
所述信号分析模块对IQ信号进行杂波识别,判断IQ信号是否为地物信号或气象信号,并对IQ信号进行标记,地物信号标记1,气象信号标记0;所述信号分析模块将标记后的IQ信号发送至滤波模块;
所述滤波模块对标记后的IQ信号进行滤波以滤除地物信号,并将滤波前的IQ信号和滤波后的IQ信号均发送至信号处理模块;所述滤波前的IQ信号为原始信号即标记后的IQ信号即包括地物信号和气象信号;所述滤波后的IQ信号为气象信号;
所述信号处理模块分别对滤波前和滤波后的IQ信号均进行多普勒参数计算,分别得到IQ信号在滤波前和滤波后的强度、速度、谱宽数据,并分别将IQ信号在滤波前和滤波后的强度、速度、谱宽数据发送至质量控制模块;
所述质量控制模块接收IQ信号在滤波前和滤波后的强度、速度、谱宽数据,并根据所接收的数据计算该IQ信号的质量控制因子,所述质量控制模块根据质量控制因子门限对该IQ信号进行判断,若该IQ信号的质量控制因子大于质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;所述质量控制模块将属于气象信息的IQ信号在滤波后的强度、速度、谱宽数据发送至终端显控系统;
所述参数配置模块用于向质量控制模块提供质量控制因子门限,所述质量控制因子门限为用户所设置。
所述信号接收模块接收的IQ信号包括四个通道的IQ信号,分别为:水平低增益通道的IQ信号Vh_low、水平高增益通道的IQ信号Vh_high、垂直低增益通道的IQ信号Vv_low、垂直高增益通道的IQ信号Vv_high;
所述信号分析模块先对水平低增益通道的IQ信号Vh_low和垂直低增益通道的IQ信号Vv_low均进行修正即进行增益补偿;
所述修正具体包括以下步骤:
S1,分别计算水平通道增益差GΔH和垂直通道增益差GΔV,且水平通道增益差GΔH和垂直通道增益差GΔV均为一个复数,表达式分别为:GΔH=IΔH+j*QΔH;GΔV=IΔV+j*QΔV;
根据水平通道增益差GΔH的表达式分别计算水平通道的功率值即幅值AΔH和相位θΔH,AΔH=IΔH 2+QΔH 2,θΔH=arctan(QΔH/IΔH);
根据垂直通道增益差GΔV的表达式分别计算直通道对应的功率值即幅值AΔV和相位θΔV,AΔV=IΔV 2+QΔV 2,θΔV=arctan(QΔV/IΔV);其中,
S11,计算水平通道增益差GΔH,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的水平低增益通道的IQ信号计算每个距离库上的水平低增益通道的功率谱:
其中,M表示脉冲点数;Vh_low(m)表示M个脉冲点中的第m个脉冲点的水平低增益通道的IQ信号;Ph_low表示水平低增益通道的功率谱
对测试库中的水平低增益通道的功率谱求平均值Rh0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Ph_low(i)表示Nd个距离库中的第i个距离库上的水平低增益通道的功率谱;
计算水平通道增益差GΔH:
其中,conj(·)表示共轭;Vh_low(i)表示M个脉冲点下的Nd个距离库的水平低增益通道的IQ信号中的第i个水平低增益通道的IQ信号;Vh_high(i)表示M个脉冲点下的Nd个距离库的水平高增益通道的IQ信号中的第i个水平高增益通道的IQ信号;
S12,计算垂直通道增益差GΔV,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的垂直低增益通道的IQ信号计算每个距离库上的垂直低增益通道的功率谱:
其中,M表示脉冲点数;Vv_low(m)表示M个脉冲点中的第m个脉冲点的垂直低增益通道的IQ信号;Pv_low表示垂直低增益通道的功率谱;
对测试库中的垂直低增益通道的功率谱求平均值Rv0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Pv_low(i)表示Nd个距离库中的第i个距离库上的垂直低增益通道的功率谱;
计算垂直通道增益差GΔV:
其中,conj(·)表示共轭;Vv_low(i)表示M个脉冲点下的Nd个距离库的垂直低增益通道的IQ信号中的第i个垂直低增益通道的IQ信号;Vv_high(i)表示M个脉冲点下的Nd个距离库的垂直高增益通道的IQ信号中的第i个垂直高增益通道的IQ信号;
S2,将水平通道的增益补偿功率值A1与水平通道的功率值AΔH相乘,得到修正后的水平通道的功率值A′ΔH,即A′ΔH=A1*AΔH;将水平通道的增益补偿相位θ1与水平通道的相位θΔH相加,得到修正后的水平通道的相位θ′ΔH,即θ′ΔH=θ1+θΔH;根据修正后的水平通道的功率值A′ΔH和修正后的水平通道的相位θ′ΔH得到水平通道增益值GH,表达式为GH=IH+j*QH,其中,IH=A′ΔH*cosθ′ΔH,QH=A′ΔH*sinθ′ΔH;
将垂直通道的增益补偿功率值A2与垂直通道的功率值AΔV相乘,得到修正后的垂直通道的功率值A′ΔV,即A′ΔV=A2*AΔV;将垂直通道的增益补偿相位θ2与垂直通道的相位θΔV相加,得到修正后的垂直通道的相位θ′ΔV,即θ′ΔV=θ2+θΔV;根据修正后的垂直通道的功率值A′ΔV和修正后的垂直通道的相位θ′ΔV得到垂直通道增益值GV,表达式为GV=IV+j*QV,其中,IV=A′ΔV*cosθ′ΔV,QV=A′ΔV*sinθ′ΔV;
其中,水平通道的增益补偿功率值A1和增益补偿相位θ1,以及垂直通道的增益补偿功率值A2和增益补偿相位θ2均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S3,计算水平通道增益值GH与水平通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于水平通道增益故障门限,则使用水平通道的增益适配值作为水平通道的修正值;否则,使用水平通道增益值作为水平通道的修正值;
计算垂直通道增益值GV与垂直通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于垂直通道增益故障门限,则使用垂直通道的增益适配值作为垂直通道的修正值;否则,使用垂直通道增益值作为垂直通道的修正值;
其中,水平通道的增益适配值和垂直通道的增益适配值,以及水平通道增益故障门限和垂直通道增益故障门限也均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S4,水平低增益通道的IQ信号Vh_low的修正:对于水平低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以水平通道的修正值;
垂直低增益通道的IQ信号Vv_low的修正:对于垂直低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以垂直通道的修正值。
所述信号分析模块对每个距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号进行水平高低增益拼接,在同一距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的水平通道信号;
其中,第j个距离库的水平高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择水平高增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的水平高增益通道的IQ信号的强度值是否低于水平高增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的水平低增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的修正后的水平低增益通道的IQ信号的强度值是否低于水平低增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的水平通道信号;
所述信号分析模块对每个距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号进行垂直高低增益拼接,在同一距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的垂直通道信号;
且所述垂直高低增益拼接的方式与所述水平高低增益拼接的方式相同,
其中,第j个距离库的垂直高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择垂直高增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的垂直高增益通道的IQ信号的强度值是否低于垂直高增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的垂直低增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的修正后的垂直低增益通道的IQ信号的强度值是否低于垂直低增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的垂直通道信号;
所述信号分析模块对高低增益拼接后的IQ信号利用CMD算法进行杂波识别;所述CMD算法为:通过计算特征量反射率纹理TDBZ、杂波相位阵列校准值CPA、反映反射率因子在径向梯度变化的SPIN,以计算地物概率CP,并将地物概率CP超过设定的门限的IQ信号判断为地物信号。
所述质量控制因子包括:LOG1、LOG2、SQI、CCOR、SIG;其中,LOG1用于估计滤波前的IQ信号的强度,LOG2用于估计滤波后的IQ信号的强度,SQI用于测量对IQ信号相关性,CCOR用于杂波校准,SIG用于估计IQ信号的信噪比;所述质量控制因子的具体计算方式如下所示:
其中,T0表示IQ信号在滤波前的强度;N表示噪声信号的幅度值;R0表示IQ信号在滤波后的零阶自相关;R1表示IQ信号在滤波后的一阶自相关;
所述质量控制因子门限包括:LOG1的门限、LOG2的门限、SQI的门限、CCOR的门限、SIG的门限;
若该IQ信号的四个质量控制因子均大于其相对应的质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;
所述质量控制模块在将属于气象信息的IQ信号发送至终端显控系统前,还将属于气象信息的IQ信号通过孤噪处理器进行孤噪处理,去除处理属于气象信息的IQ信号中的离散噪声点。
LOG1的门限和LOG2的门限取值均为0.75db;SQI的门限的取值范围为0.4~0.5;CCOR的门限取值为-18db;SIG的门限取值为5db。
用户还通过所述参数配置模块为所述滤波模块选择滤波方式;所述滤波方式包括IIR滤波、固定宽度和自适应谱处理滤波、CMD结合GMAP的滤波。
所述信号处理单元还包括信号回放模块,若终端显控系统所显示的IQ信号的强度、速度、谱宽数据未达到设定的要求时,则用户在信号回放模块中对信号处理参数进行重新设置,所述信号回放模块利用新的信号处理参数对信号存储模块中的IQ信号进行重新处理,并将重新处理后的属于气象信息的IQ信号发送至终端显控系统;所述信号处理参数为各模块进行信号处理时所需的参数,包括:质量控制因子门限、滤波方式、增益差补偿、增益适配值、增益故障门限。
所述信号处理器采用通用服务器作为信号处理单元的硬件平台,所述信号处理器采用Qt工具作为信号处理单元的开发工具;所述通用服务器的操作系统采用中标麒麟操作系统。
所述信号处理服务器具有千兆以太网卡,所述信号接收模块通过光纤并按巨型帧格式接收雷达接收机所发送的IQ信号的数据。
所述信号处理器的CPU具有多个内核,所述信号处理单元为每个模块均创建线程,且每个线程均运行在不同的内核上。
本发明的优点在于:
(1)采用通用服务器作为信号处理单元的硬件平台,运算速度快、扩展性强。
(2)通过参数设置模块并根据雷达实际回波选择合适的质量因子门限,更为有效便捷的加强了数据质量控制。
(3)通过孤噪处理器使得在质量控制因子门限降低的情况下,保证系统有更高的灵敏度和更少的虚警点。
(4)本发明采用软件化设计,开放大量的雷达信号处理的参数供用户设置,针对各模块进行处理时所需的参数,用户通过终端显控系统的人机界面即参数配置模块进行设置,通过软件的重构和不断改进,以提升或扩展雷达系统功能,实现快速多元化探测。
(5)信号处理器具有千兆以太网卡,按巨型帧格式接收IQ信号的数据,减少了网络中数据包的个数,减轻了网络设备处理包头的额外的资源占用,提高了信号处理效率;通过光纤来接收IQ信号的数据,提高了数据传输速率,缩减了整机外部的电缆数量,提高了数据传输的抗干扰能力。
(6)本发明为多个模块创建不同的线程运行在不同的内核上,方便了二次雷达的信号处理算法的研制开发,降低了仿真调试的难度。
附图说明
图1为本发明的整体架构图。
图2为通道转换的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的一种多普勒天气雷达的信号处理器,所述信号处理器采用通用服务器作为信号处理单元的硬件平台,所述通用服务器的操作系统采用中标麒麟操作系统。本实施例中,所述信号处理器采用麒麟系统服务器,型号为HPDL380Gen9;CPU为2颗/十核/2.5GHz以上;内存为64G以上ECC DDR3;硬盘为3T以上;2块千兆以太网卡。
所述信号处理器采用Qt工具作为信号处理单元的开发工具,Qt中的C++编程语言灵活,运算符的数据结构丰富、具有结构化控制语句、程序执行效率高,而且同时具有高级语言与汇编语言的优点,与其它语言相比,可以直接访问物理地址,与汇编语言相比又具有良好的可读性和可移植性。
所述信号处理器调用IPP库,IPP库为一套跨平台的软件函数库,信号处理单元的各个模块的相关计算调用IPP库中的算法函数,并将各算法函数封装成独立的插件,具有标准的输入输出接口,用户在进行算法优化时:第一步按照插件定义的输入输出进行算法设计和编码,第二步进行软件模块的编译,第三部进行相应算法函数的替换;极大改进了信号处理单元中各模块的性能。
所述信号处理器具有千兆以太网卡,所述信号处理器通过光纤并按巨型帧格式接收雷达接收机所发送的IQ信号的数据,减少了网络中数据包的个数。
所述信号处理器预留双偏振接口,当多普勒雷达升级成双偏振雷达时,雷达发射水平和垂直偏振状态电磁波时,信号处理器对接收的IQ信号仍可直接进行信号处理。
由图1所示,一种多普勒天气雷达的信号处理器,包括信号处理单元,
所述信号处理单元包括如下组成部分:信号接收模块1、信号分析模块2、信号存储模块3、滤波模块4、信号处理模块5、质量控制模块6、晴空图模块7、参数配置模块8、信号回放模块9。
信号处理单元为每个模块均创建线程或线程池,且每个线程或线程池均运行在不同的内核上。
所述信号接收模块1接收来源于雷达接收机的IQ信号,并将IQ信号分别发送至信号存储模块3和信号分析模块2;所述IQ信号为一种正交的基带信号。
所述信号接收模块1接收的IQ信号包括四个通道的IQ信号,分别为:水平低增益通道的IQ信号Vh_low、水平高增益通道的IQ信号Vh_high、垂直低增益通道的IQ信号Vv_low、垂直高增益通道的IQ信号Vv_high。
所述信号存储模块3用于存储所述信号接收模块1所发送的IQ信号。
由图2所示,所述信号分析模块2对IQ信号进行通道转换,将发送来的按距离库排列的四通道的IQ信号,转成按通道排列的每个距离库的IQ信号。
所述信号分析模块2对水平低增益通道的IQ信号Vh_low和垂直低增益通道的IQ信号Vv_low均进行修正即进行增益补偿;且所述水平高增益通道的IQ信号Vh_high和所述垂直高增益通道的IQ信号Vv_high不进行修正;
所述修正具体包括以下步骤:
S1,分别计算水平通道增益差GΔH和垂直通道增益差GΔV,且水平通道增益差GΔH和垂直通道增益差GΔV均为一个复数,表达式分别为:GΔH=IΔH+j*QΔH;GΔV=IΔV+j*QΔV;
根据水平通道增益差GΔH的表达式分别计算水平通道的功率值即幅值AΔH和相位θΔH,AΔH=IΔH 2+QΔH 2,θΔH=arctan(QΔH/IΔH);
根据垂直通道增益差GΔV的表达式分别计算直通道对应的功率值即幅值AΔV和相位θΔV,AΔV=IΔV 2+QΔV 2,θΔV=arctan(QΔV/IΔV);其中,
S11,计算水平通道增益差GΔH,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的水平低增益通道的IQ信号计算每个距离库上的水平低增益通道的功率谱:
其中,M表示脉冲点数;Vh_low(m)表示M个脉冲点中的第m个脉冲点的水平低增益通道的IQ信号;Ph_low表示水平低增益通道的功率谱
对测试库中的水平低增益通道的功率谱求平均值Rh0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Ph_low(i)表示Nd个距离库中的第i个距离库上的水平低增益通道的功率谱;
计算水平通道增益差GΔH:
其中,conj(·)表示共轭;Vh_low(i)表示M个脉冲点下的Nd个距离库的水平低增益通道的IQ信号中的第i个水平低增益通道的IQ信号;Vh_high(i)表示M个脉冲点下的Nd个距离库的水平高增益通道的IQ信号中的第i个水平高增益通道的IQ信号;
S12,计算垂直通道增益差GΔV,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的垂直低增益通道的IQ信号计算每个距离库上的垂直低增益通道的功率谱:
其中,M表示脉冲点数;Vv_low(m)表示M个脉冲点中的第m个脉冲点的垂直低增益通道的IQ信号;Pv_low表示垂直低增益通道的功率谱;
对测试库中的垂直低增益通道的功率谱求平均值Rv0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Pv_low(i)表示Nd个距离库中的第i个距离库上的垂直低增益通道的功率谱;
计算垂直通道增益差GΔV:
其中,conj(·)表示共轭;Vv_low(i)表示M个脉冲点下的Nd个距离库的垂直低增益通道的IQ信号中的第i个垂直低增益通道的IQ信号;Vv_high(i)表示M个脉冲点下的Nd个距离库的垂直高增益通道的IQ信号中的第i个垂直高增益通道的IQ信号;
S2,将水平通道的增益补偿功率值A1与水平通道的功率值AΔH相乘,得到修正后的水平通道的功率值A′ΔH,即A′ΔH=A1*AΔH;将水平通道的增益补偿相位θ1与水平通道的相位θΔH相加,得到修正后的水平通道的相位θ′ΔH,即θ′ΔH=θ1+θΔH;根据修正后的水平通道的功率值A′ΔH和修正后的水平通道的相位θ′ΔH得到水平通道增益值GH,表达式为GH=IH+j*QH,其中,IH=A′ΔH*cosθ′ΔH,QH=A′ΔH*sinθ′ΔH;
将垂直通道的增益补偿功率值A2与垂直通道的功率值AΔV相乘,得到修正后的垂直通道的功率值A′ΔV,即A′ΔV=A2*AΔV;将垂直通道的增益补偿相位θ2与垂直通道的相位θΔV相加,得到修正后的垂直通道的相位θ′ΔV,即θ′ΔV=θ2+θΔV;根据修正后的垂直通道的功率值A′ΔV和修正后的垂直通道的相位θ′ΔV得到垂直通道增益值GV,表达式为GV=IV+j*QV,其中,IV=A′ΔV*cosθ′ΔV,QV=A′ΔV*sinθ′ΔV;
其中,水平通道的增益补偿功率值A1和增益补偿相位θ1,以及垂直通道的增益补偿功率值A2和增益补偿相位θ2均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S3,计算水平通道增益值GH与水平通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于水平通道增益故障门限,则使用水平通道的增益适配值作为水平通道的修正值;否则,使用水平通道增益值作为水平通道的修正值;
计算垂直通道增益值GV与垂直通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于垂直通道增益故障门限,则使用垂直通道的增益适配值作为垂直通道的修正值;否则,使用垂直通道增益值作为垂直通道的修正值;
其中,水平通道的增益适配值和垂直通道的增益适配值,以及水平通道增益故障门限和垂直通道增益故障门限也均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S4,水平低增益通道的IQ信号Vh_low的修正:对于水平低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以水平通道的修正值;
垂直低增益通道的IQ信号Vv_low的修正:对于垂直低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以垂直通道的修正值。
所述信号分析模块2对每个距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号进行水平高低增益拼接,在同一距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的水平通道信号;
其中,第j个距离库的水平高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择水平高增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的水平高增益通道的IQ信号的强度值是否低于水平高增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的水平低增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的修正后的水平低增益通道的IQ信号的强度值是否低于水平低增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的水平通道信号;
所述信号分析模块2对每个距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号进行垂直高低增益拼接,在同一距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的垂直通道信号;
且所述垂直高低增益拼接的方式与所述水平高低增益拼接的方式相同,
其中,第j个距离库的垂直高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择垂直高增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的垂直高增益通道的IQ信号的强度值是否低于垂直高增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的垂直低增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的修正后的垂直低增益通道的IQ信号的强度值是否低于垂直低增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的垂直通道信号。
所述信号分析模块2利用CMD算法对修正后的IQ信号进行杂波识别,判断IQ信号为地物信号或气象信号,并对IQ信号进行标记,地物信号标记1,气象信号标记0;所述信号分析模块2将标记后的IQ信号发送至滤波模块4;所述CMD算法具体方式如下所示:通过计算特征量反射率纹理TDBZ、杂波相位阵列校准值CPA、反映反射率因子在径向梯度变化的SPIN,以计算地物概率CP,并将地物概率CP超过0.5的IQ信号标识为地物信号。
所述滤波模块4对标记后的IQ信号进行滤波以滤除地物信号,并将滤波前的IQ信号和滤波后的IQ信号均发送至信号处理模块5,还将滤波后的IQ信号发送至晴空图模块7;所述滤波前的IQ信号为原始信号即标记后的IQ信号即包括地物信号和气象信号;所述滤波后的IQ信号为气象信号;所述信号处理单元为所述滤波模块4创建了滤波线程池,且滤波线程池中包括八个滤波线程,每个滤波线程处理八分之一的距离库上的IQ信号。
所述晴空图模块7根据滤波后的IQ信号建立晴空图和估算气象等级,并再终端显控系统进行显示。
所述信号处理模块5分别对滤波前的IQ信号和滤波后的IQ信号进行多普勒参数计算,分别得到IQ信号在滤波前的强度、速度、谱宽数据和IQ信号在滤波后的强度、速度、谱宽数据,并分别IQ信号在滤波前和滤波后的强度、速度、谱宽数据发送至质量控制模块6。
所述质量控制模块6接收IQ信号在滤波前和滤波后的强度、速度、谱宽数据,并根据所接收的数据计算该IQ信号的质量控制因子,所述质量控制模块根据质量控制因子门限对该IQ信号进行判断,若该IQ信号的质量控制因子大于质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;所述质量控制模块还将属于气象信息的IQ信号通过孤噪处理器进行孤噪处理,去除处理属于气象信息的IQ信号中的离散噪声点;所述孤噪处理为采用一维或二维的处理方式;所述质量控制模块6将属于孤噪处理后的气象信息的IQ信号发送至终端显控系统;
所述质量控制因子包括:LOG1、LOG2、SQI、CCOR、SIG;所述质量控制因子的具体计算方式如下所示:
其中,T0表示IQ信号在滤波前的强度;N表示噪声信号的幅度值;R0表示IQ信号在滤波后的零阶自相关;R1表示IQ信号在滤波后的一阶自相关;
所述质量控制因子门限包括:LOG1的门限、LOG2的门限、SQI的门限、CCOR的门限、SIG的门限;其中,LOG1的门限和LOG2的门限取值均为0.75db;SQI的门限的取值范围为0.4~0.5;CCOR的门限取值为-18db;SIG的门限取值为5db;
若该IQ信号的四个质量控制因子均大于其相对应质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;
其中,LOG1用于估计滤波前的IQ信号的强度,LOG2用于估计滤波后的IQ信号的强度,LOG1的门限、LOG2的门限的设置用于反射率的质量控制;SQI用于测量对IQ信号相关性,SQI的门限一般用于速度和谱宽的检测;CCOR用于杂波校准,CCOR的门限的设置用于删除杂波非常强的回波信号,当杂波非常强,计算出的CCOR是一个非常小的负值;SIG用于估计IQ信号的信噪比,SIG的门限的设置主要用于谱宽的质量控制,以确保气象功率足够大,谱宽测量的精度可以满足要求。
用户从终端显控系统并通过所述参数配置模块8设置所述质量控制模块所需的质量控制因子门限;用户还通过所述参数配置模块8为所述滤波模块4选择滤波方式;所述滤波方式包括IIR滤波、固定宽度和自适应谱处理滤波、CMD算法结合GMAP算法的滤波;所述CMD算法通过计算特征量反射率纹理TDBZ、杂波相位阵列校准值CPA、反映反射率因子在径向梯度变化的SPIN,以计算地物概率CP,并将地物概率CP超过0.5的IQ信号标识为地物信号;所述GMAP(Gaussian model adaptive processing)算法是在2004年,由SIGMET公司的两位工程师A.D.Siggia和R.E.Passarelli提出的,GMAP算法旨在更彻底的滤除地物杂波,同时保留零频附近的气象回波;所述CMD结合GMAP的滤波是在滤除地物杂波的基础上,保留零频附近的气象回波。
当所述显控终端显示的IQ信号的强度、速度、谱宽数据未达到设定的要求时,用户通过所述信号回放模块9中对信号处理参数进行重新设置,所述信号回放模块利用新的信号处理参数对信号存储模块3中的IQ信号进行重新处理,并将重新处理后的属于气象信息的IQ信号发送至终端显控系统;所述信号处理参数为各模块进行信号处理时所需的参数,包括:质量控制因子门限、滤波方式、增益差补偿、增益适配值、增益故障门限、水平高增益门限值、水平低增益门限值、垂直高增益门限值、垂直低增益门限值。
以上仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明创造的保护范围之内。
Claims (9)
1.一种多普勒天气雷达的信号处理器,其特征在于,包括信号处理单元,所述信号处理单元包括以下组成部分:
信号接收模块、信号存储模块、信号分析模块、滤波模块、信号处理模块、质量控制模块、参数配置模块;其中,
所述信号接收模块接收来源于雷达接收机的IQ信号,并将IQ信号分别发送至信号存储模块和信号分析模块;所述IQ信号为一种正交的基带信号;
所述信号分析模块对IQ信号进行杂波识别,判断IQ信号是否为地物信号或气象信号,并对IQ信号进行标记,地物信号标记1,气象信号标记0;所述信号分析模块将标记后的IQ信号发送至滤波模块;
所述滤波模块对标记后的IQ信号进行滤波以滤除地物信号,并将滤波前的IQ信号和滤波后的IQ信号均发送至信号处理模块;所述滤波前的IQ信号为原始信号即标记后的IQ信号即包括地物信号和气象信号;所述滤波后的IQ信号为气象信号;
所述信号处理模块分别对滤波前和滤波后的IQ信号均进行多普勒参数计算,分别得到IQ信号在滤波前和滤波后的强度、速度、谱宽数据,并分别将IQ信号在滤波前和滤波后的强度、速度、谱宽数据发送至质量控制模块;
所述质量控制模块接收IQ信号在滤波前和滤波后的强度、速度、谱宽数据,并根据所接收的数据计算该IQ信号的质量控制因子,所述质量控制模块根据质量控制因子门限对该IQ信号进行判断,若该IQ信号的质量控制因子大于质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;所述质量控制模块将属于气象信息的IQ信号在滤波后的强度、速度、谱宽数据发送至终端显控系统;
所述参数配置模块用于向质量控制模块提供质量控制因子门限,所述质量控制因子门限为用户所设置;
所述信号接收模块接收的IQ信号包括四个通道的IQ信号,分别为:水平低增益通道的IQ信号Vh_low、水平高增益通道的IQ信号Vh_high、垂直低增益通道的IQ信号Vv_low、垂直高增益通道的IQ信号Vv_high;
所述信号分析模块先对水平低增益通道的IQ信号Vh_low和垂直低增益通道的IQ信号Vv_low均进行修正即进行增益补偿;
所述修正具体包括以下步骤:
S1,分别计算水平通道增益差GΔH和垂直通道增益差GΔV,且水平通道增益差GΔH和垂直通道增益差GΔV均为一个复数,表达式分别为:GΔH=IΔH+j*QΔH;GΔV=IΔV+j*QΔV;
根据水平通道增益差GΔH的表达式分别计算水平通道的功率值即幅值AΔH和相位θΔH,AΔH=IΔH 2+QΔH 2,θΔH=arctan(QΔH/IΔH);
根据垂直通道增益差GΔV的表达式分别计算直通道对应的功率值即幅值AΔV和相位θΔV,AΔV=IΔV 2+QΔV 2,θΔV=arctan(QΔV/IΔV);其中,
S11,计算水平通道增益差GΔH,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的水平低增益通道的IQ信号计算每个距离库上的水平低增益通道的功率谱:
其中,M表示脉冲点数;Vh_low(m)表示M个脉冲点中的第m个脉冲点的水平低增益通道的IQ信号;Ph_low表示水平低增益通道的功率谱
对测试库中的水平低增益通道的功率谱求平均值Rh0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Ph_low(i)表示Nd个距离库中的第i个距离库上的水平低增益通道的功率谱;
计算水平通道增益差GΔH:
其中,conj(·)表示共轭;Vh_low(i)表示M个脉冲点下的Nd个距离库的水平低增益通道的IQ信号中的第i个水平低增益通道的IQ信号;Vh_high(i)表示M个脉冲点下的Nd个距离库的水平高增益通道的IQ信号中的第i个水平高增益通道的IQ信号;
S12,计算垂直通道增益差GΔV,具体计算方式如下所示:
根据每个距离库上的M个脉冲点的垂直低增益通道的IQ信号计算每个距离库上的垂直低增益通道的功率谱:
其中,M表示脉冲点数;Vv_low(m)表示M个脉冲点中的第m个脉冲点的垂直低增益通道的IQ信号;Pv_low表示垂直低增益通道的功率谱;
对测试库中的垂直低增益通道的功率谱求平均值Rv0:
其中,Nd表示测试库的库数,所述测试库为从所有距离库中选择参数配置模块指定的Nd个距离库作为测试库;Pv_low(i)表示Nd个距离库中的第i个距离库上的垂直低增益通道的功率谱;
计算垂直通道增益差GΔV:
其中,conj(·)表示共轭;Vv_low(i)表示M个脉冲点下的Nd个距离库的垂直低增益通道的IQ信号中的第i个垂直低增益通道的IQ信号;Vv_high(i)表示M个脉冲点下的Nd个距离库的垂直高增益通道的IQ信号中的第i个垂直高增益通道的IQ信号;
S2,将水平通道的增益补偿功率值A1与水平通道的功率值AΔH相乘,得到修正后的水平通道的功率值A′ΔH,即A′ΔH=A1*AΔH;将水平通道的增益补偿相位θ1与水平通道的相位θΔH相加,得到修正后的水平通道的相位θ′ΔH,即θ′ΔH=θ1+θΔH;根据修正后的水平通道的功率值A′ΔH和修正后的水平通道的相位θ′ΔH得到水平通道增益值GH,表达式为GH=IH+j*QH,其中,IH=A′ΔH*cosθ′ΔH,QH=A′ΔH*sinθ′ΔH;
将垂直通道的增益补偿功率值A2与垂直通道的功率值AΔV相乘,得到修正后的垂直通道的功率值A′ΔV,即A′ΔV=A2*AΔV;将垂直通道的增益补偿相位θ2与垂直通道的相位θΔV相加,得到修正后的垂直通道的相位θ′ΔV,即θ′ΔV=θ2+θΔV;根据修正后的垂直通道的功率值A′ΔV和修正后的垂直通道的相位θ′ΔV得到垂直通道增益值GV,表达式为GV=IV+j*QV,其中,IV=A′ΔV*cosθ′ΔV,QV=A′ΔV*sinθ′ΔV;
其中,水平通道的增益补偿功率值A1和增益补偿相位θ1,以及垂直通道的增益补偿功率值A2和增益补偿相位θ2均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S3,计算水平通道增益值GH与水平通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于水平通道增益故障门限,则使用水平通道的增益适配值作为水平通道的修正值;否则,使用水平通道增益值作为水平通道的修正值;
计算垂直通道增益值GV与垂直通道的增益适配值之间的差值的绝对值,若该差值的绝对值大于垂直通道增益故障门限,则使用垂直通道的增益适配值作为垂直通道的修正值;否则,使用垂直通道增益值作为垂直通道的修正值;
其中,水平通道的增益适配值和垂直通道的增益适配值,以及水平通道增益故障门限和垂直通道增益故障门限也均为用户通过参数配置模块进行设置的,并存储在参数配置模块中;
S4,水平低增益通道的IQ信号Vh_low的修正:对于水平低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以水平通道的修正值;
垂直低增益通道的IQ信号Vv_low的修正:对于垂直低增益通道上的所有距离库的所有脉冲点数的IQ信号均乘以垂直通道的修正值。
2.根据权利要求1所述的一种多普勒天气雷达的信号处理器,其特征在于,
所述信号分析模块对每个距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号进行水平高低增益拼接,在同一距离库上的水平高增益通道的IQ信号和修正后的水平低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的水平通道信号;
其中,第j个距离库的水平高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择水平高增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的水平高增益通道的IQ信号的强度值是否低于水平高增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的水平低增益通道的IQ信号作为该距离库上的水平通道信号,则判断第j个距离库即当前距离库上的修正后的水平低增益通道的IQ信号的强度值是否低于水平低增益门限值,若低于,则选择水平高增益通道的IQ信号作为第j个距离库的水平通道信号;否则,选择修正后的水平低增益通道的IQ信号作为第j个距离库的水平通道信号;其中,用户通过参数配置模块设置所述水平低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的水平通道信号;
所述信号分析模块对每个距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号进行垂直高低增益拼接,在同一距离库上的垂直高增益通道的IQ信号和修正后的垂直低增益通道的IQ信号中仅选出一个IQ信号作为此距离库上的垂直通道信号;
且所述垂直高低增益拼接的方式与所述水平高低增益拼接的方式相同,
其中,第j个距离库的垂直高低增益拼接的具体方式为:
若第j-1个距离库即上一个距离库选择垂直高增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的垂直高增益通道的IQ信号的强度值是否低于垂直高增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直高增益门限值;
若第j-1个距离库即上一个距离库选择修正后的垂直低增益通道的IQ信号作为该距离库上的垂直通道信号,则判断第j个距离库即当前距离库上的修正后的垂直低增益通道的IQ信号的强度值是否低于垂直低增益门限值,若低于,则选择垂直高增益通道的IQ信号作为第j个距离库的垂直通道信号;否则,选择修正后的垂直低增益通道的IQ信号作为第j个距离库的垂直通道信号;其中,用户通过参数配置模块设置所述垂直低增益门限值;
且当j=1时,即第一个距离库时,选择高增益通道的IQ信号作为第一个距离库的垂直通道信号;
所述信号分析模块对高低增益拼接后的IQ信号利用CMD算法进行杂波识别;所述CMD算法为:通过计算特征量反射率纹理TDBZ、杂波相位阵列校准值CPA、反映反射率因子在径向梯度变化的SPIN,以计算地物概率CP,并将地物概率CP超过设定的门限的IQ信号判断为地物信号。
3.根据权利要求1所述的一种多普勒天气雷达的信号处理器,其特征在于,所述质量控制因子包括:LOG1、LOG2、SQI、CCOR、SIG;其中,LOG1用于估计滤波前的IQ信号的强度,LOG2用于估计滤波后的IQ信号的强度,SQI用于测量对IQ信号相关性,CCOR用于杂波校准,SIG用于估计IQ信号的信噪比;所述质量控制因子的具体计算方式如下所示:
其中,T0表示IQ信号在滤波前的强度;N表示噪声信号的幅度值;R0表示IQ信号在滤波后的零阶自相关;R1表示IQ信号在滤波后的一阶自相关;
所述质量控制因子门限包括:LOG1的门限、LOG2的门限、SQI的门限、CCOR的门限、SIG的门限;
若该IQ信号的四个质量控制因子均大于其相对应的质量控制因子门限,则该IQ信号属于气象信息,否则属于噪声信息;
所述质量控制模块在将属于气象信息的IQ信号发送至终端显控系统前,还将属于气象信息的IQ信号通过孤噪处理器进行孤噪处理,去除处理属于气象信息的IQ信号中的离散噪声点。
4.根据权利要求3所述的一种多普勒天气雷达的信号处理器,其特征在于,LOG1的门限和LOG2的门限取值均为0.75db;SQI的门限的取值范围为0.4~0.5;CCOR的门限取值为-18db;SIG的门限取值为5db。
5.根据权利要求1所述的一种多普勒天气雷达的信号处理器,其特征在于,用户还通过所述参数配置模块为所述滤波模块选择滤波方式;所述滤波方式包括IIR滤波、固定宽度和自适应谱处理滤波、CMD结合GMAP的滤波。
6.根据权利要求1所述的一种多普勒天气雷达的信号处理器,其特征在于,所述信号处理单元还包括信号回放模块,若终端显控系统所显示的IQ信号的强度、速度、谱宽数据未达到设定的要求时,则用户在信号回放模块中对信号处理参数进行重新设置,所述信号回放模块利用新的信号处理参数对信号存储模块中的IQ信号进行重新处理,并将重新处理后的属于气象信息的IQ信号发送至终端显控系统;所述信号处理参数为各模块进行信号处理时所需的参数,包括:质量控制因子门限、滤波方式、增益差补偿、增益适配值、增益故障门限。
7.根据权利要求1所述的一种多普勒天气雷达的信号处理器,其特征在于,所述信号处理器采用通用服务器作为信号处理单元的硬件平台,所述信号处理器采用Qt工具作为信号处理单元的开发工具;所述通用服务器的操作系统采用中标麒麟操作系统。
8.根据权利要求7所述的一种多普勒天气雷达的信号处理器,其特征在于,所述信号处理服务器具有千兆以太网卡,所述信号接收模块通过光纤并按巨型帧格式接收雷达接收机所发送的IQ信号的数据。
9.根据权利要求7所述的一种多普勒天气雷达的信号处理器,其特征在于,所述信号处理器的CPU具有多个内核,所述信号处理单元为每个模块均创建线程,且每个线程均运行在不同的内核上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811287767.1A CN109270536B (zh) | 2018-10-31 | 2018-10-31 | 一种多普勒天气雷达的信号处理器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811287767.1A CN109270536B (zh) | 2018-10-31 | 2018-10-31 | 一种多普勒天气雷达的信号处理器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109270536A CN109270536A (zh) | 2019-01-25 |
CN109270536B true CN109270536B (zh) | 2020-09-01 |
Family
ID=65191301
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811287767.1A Active CN109270536B (zh) | 2018-10-31 | 2018-10-31 | 一种多普勒天气雷达的信号处理器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109270536B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112444814B (zh) * | 2020-11-11 | 2022-07-05 | 安徽四创电子股份有限公司 | 基于pcie光纤采集卡的数字阵列天气雷达信号处理器 |
CN112666559B (zh) * | 2021-01-27 | 2022-09-02 | 苏州市气象局 | 一种多普勒天气雷达数据质量控制方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0613022A2 (en) * | 1993-02-26 | 1994-08-31 | AlliedSignal Inc. | Airborne weather radar system with aircraft altitude and antenna tilt angle dependent sensitivity time control |
US6307500B1 (en) * | 1999-08-13 | 2001-10-23 | University Corporation For Atmospheric Research | Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time |
CN201285802Y (zh) * | 2008-07-03 | 2009-08-05 | 冯超 | X波段双偏振雷达天线 |
CN103323850A (zh) * | 2013-05-28 | 2013-09-25 | 芜湖航飞科技股份有限公司 | 一种双线偏振多普勒天气雷达系统 |
CN103969635A (zh) * | 2014-04-30 | 2014-08-06 | 上海航天电子通讯设备研究所 | 低空监视雷达的气象信号处理ip核及其实时数据排序方法 |
CN104678369A (zh) * | 2015-01-20 | 2015-06-03 | 南京大学 | 一种基于非固定金属球的双偏振天气雷达标校方法 |
CN105548984A (zh) * | 2015-12-10 | 2016-05-04 | 中国人民解放军理工大学 | 一种基于模糊逻辑的双偏振多普勒天气雷达地物杂波抑制方法 |
CN105891833A (zh) * | 2016-04-15 | 2016-08-24 | 兰州大学 | 基于多普勒雷达信息识别暖云降水率的方法 |
CN107843895A (zh) * | 2017-10-20 | 2018-03-27 | 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷检测技术中心) | 一种双多普勒雷达三维风场反演方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5523759A (en) * | 1992-04-14 | 1996-06-04 | Honeywell Inc. | In flight doppler weather radar wind shear detection system |
CN100520444C (zh) * | 2006-09-26 | 2009-07-29 | 南京大桥机器有限公司 | 便携式x波段多普勒天气雷达信号处理方法及装置 |
-
2018
- 2018-10-31 CN CN201811287767.1A patent/CN109270536B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0613022A2 (en) * | 1993-02-26 | 1994-08-31 | AlliedSignal Inc. | Airborne weather radar system with aircraft altitude and antenna tilt angle dependent sensitivity time control |
US6307500B1 (en) * | 1999-08-13 | 2001-10-23 | University Corporation For Atmospheric Research | Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time |
CN201285802Y (zh) * | 2008-07-03 | 2009-08-05 | 冯超 | X波段双偏振雷达天线 |
CN103323850A (zh) * | 2013-05-28 | 2013-09-25 | 芜湖航飞科技股份有限公司 | 一种双线偏振多普勒天气雷达系统 |
CN103969635A (zh) * | 2014-04-30 | 2014-08-06 | 上海航天电子通讯设备研究所 | 低空监视雷达的气象信号处理ip核及其实时数据排序方法 |
CN104678369A (zh) * | 2015-01-20 | 2015-06-03 | 南京大学 | 一种基于非固定金属球的双偏振天气雷达标校方法 |
CN105548984A (zh) * | 2015-12-10 | 2016-05-04 | 中国人民解放军理工大学 | 一种基于模糊逻辑的双偏振多普勒天气雷达地物杂波抑制方法 |
CN105891833A (zh) * | 2016-04-15 | 2016-08-24 | 兰州大学 | 基于多普勒雷达信息识别暖云降水率的方法 |
CN107843895A (zh) * | 2017-10-20 | 2018-03-27 | 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷检测技术中心) | 一种双多普勒雷达三维风场反演方法 |
Non-Patent Citations (1)
Title |
---|
Radar rainfall estimation for the identification of debris flow occurrence threshold;Francesco Marra;《journal of hydrology 》;20141130;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109270536A (zh) | 2019-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102624468B (zh) | 基于双fft的宽带自动检测方法 | |
CN109270536B (zh) | 一种多普勒天气雷达的信号处理器 | |
CN110297218B (zh) | 基于生成对抗网络的雷达信号未知调制方式检测方法 | |
CN115664568A (zh) | 测试天线性能的方法、装置、电子设备及存储介质 | |
CN109613503A (zh) | 雷达回波信号的标校方法和装置 | |
CN107045122A (zh) | 一种目标检测系统及其检测方法 | |
CN116224274A (zh) | 基于动态门限因子的低慢小目标检测方法、装置及设备 | |
CN109061632A (zh) | 一种无人机识别方法 | |
CN112444814B (zh) | 基于pcie光纤采集卡的数字阵列天气雷达信号处理器 | |
CN111539372A (zh) | 病虫害分布的监测方法、设备、存储介质及装置 | |
CN112230212B (zh) | 雷达测距信号处理方法和装置 | |
CN114254265A (zh) | 基于统计流形距离的卫星通信干扰几何分析方法 | |
CN111371436B (zh) | 雷达天线扫描周期测量方法、装置、设备及存储介质 | |
JP2024009334A (ja) | レーダーシステム | |
CN116961799A (zh) | 一种基于时频域分布特征的信号干扰检测方法 | |
CN110991378A (zh) | 一种功率放大器个体识别方法及装置 | |
CN106375039A (zh) | 一种提高矢量网络分析仪接收机动态范围的方法 | |
CN111830479B (zh) | 雷达回波双分量混合幅度分布模型参数估计方法及装置 | |
CN116125497A (zh) | 卫星导航天线阵列抗干扰方法及装置 | |
CN116400302B (zh) | 一种雷达信号接收处理方法 | |
CN113532588A (zh) | 水位获取方法、装置、设备及存储介质 | |
CN111736146A (zh) | 一种基于速度滤波的双基地检测前跟踪方法及装置 | |
CN109782233A (zh) | 一种基于傅里叶变换的雷达工作方法和系统 | |
CN116996138B (zh) | 多分辨率载波信号检测方法、系统、电子设备及存储介质 | |
CN115833941B (zh) | Dwdm光模块apd调试方法、电子设备及计算机可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |