CN104316929A - 一种功率谱识别方法及装置 - Google Patents
一种功率谱识别方法及装置 Download PDFInfo
- Publication number
- CN104316929A CN104316929A CN201410466622.3A CN201410466622A CN104316929A CN 104316929 A CN104316929 A CN 104316929A CN 201410466622 A CN201410466622 A CN 201410466622A CN 104316929 A CN104316929 A CN 104316929A
- Authority
- CN
- China
- Prior art keywords
- point
- power spectrum
- fft
- end points
- maximum
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 740
- 238000000034 method Methods 0.000 title claims abstract description 61
- 230000011664 signaling Effects 0.000 claims description 133
- 230000002902 bimodal effect Effects 0.000 claims description 61
- 238000001514 detection method Methods 0.000 abstract 1
- 238000001556 precipitation Methods 0.000 description 49
- 239000002245 particle Substances 0.000 description 4
- 230000000996 additive effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000006854 communication Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种功率谱识别方法及装置,涉及气象雷达探测领域,用于通过在双峰谱中识别出降水谱及湍流谱,提高获取大气湍流信息及降水信息的准确性。所述方法,包括:获取第一功率谱;确定所述第一功率谱是否为双峰谱;其中,所述双峰谱中包含有两个极大值点;若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点;将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱;所述第一功率谱包括所述第二功率谱及第三功率谱;在所述第一功率谱中获取所述第三功率谱。
Description
技术领域
本发明涉及气象雷达探测领域,尤其涉及一种功率谱识别方法及装置。
背景技术
风廓线雷达是一种气象雷达,可以向空中发射不同方向的电磁波束,电磁波束在传播过程中,由于大气中折射率的分布不均匀或者遇到颗粒物等而引起散射,从而产生了散射回波。风廓线雷达可以接收到散射回来的散射回波,并根据接收到的散射回波生成功率谱。
电磁波在大气传播过程中,因为大气湍流造成的折射率分布不均匀而产生散射称为湍流散射,而因为降水粒子造成的散射称为降水散射。在晴空条件下,由于大气中不存在降水粒子,所以电磁波仅因大气湍流造成的折射率分布不均匀而产生散射,此时,风廓线雷达仅接收到湍流散射的散射回波,此时风廓线雷达生成的功率谱仅与大气中的湍流散射有关,风廓线雷达生成的功率谱为单峰谱。在阴雨条件下,风廓线雷达不仅能够接收到湍流散射的散射回波,还可以接收到电磁波因降水粒子而产生降水散射的散射回波。这样,风廓线雷达根据湍流散射的散射回波及降水散射的散射回波得到功率谱,此时风廓线雷达得到的功率谱是湍流谱和降水谱叠加的双峰谱。
在阴雨条件下,风廓线雷达获取的双峰谱是降水谱与湍流谱的叠加功率谱,若通过此双峰谱获取仅与湍流谱有关的大气湍流信息时会存在较大误差。若通过此双峰普获取仅与降水谱有关的降水信息也会存在较大误差。所以在现有技术中,由于无法在双峰谱中识别出降水谱和湍流谱,从而无法通过此双峰谱准确获取大气湍流信息及降水信息。
发明内容
本发明提供了一种功率谱识别方法及装置,用于通过在双峰谱中识别出降水谱及湍流谱,提高获取大气湍流信息及降水信息的准确性。
为达到上述目的,本发明的实施例采用如下技术方案:
第一方面,本发明实施例提供了一种功率谱识别方法,包括:获取第一功率谱,所述第一功率谱是指去噪声后的雷达功率谱,确定第一功率谱是否为双峰谱,其中,所述双峰谱中包含有两个极大值点。若确定出第一功率谱为双峰谱,则确定出第一双峰谱的第一对称极值点。将第一功率谱以第一对称极值点为对称中心进行对称处理,获取第二功率谱,第一功率谱包括了第二功率谱及第三功率谱,在第一功率谱中获取第三功率谱。
第二方面,本发明实施例提供了一种功率谱识别装置,包括:获取单元,用于获取第一功率谱,所述第一功率谱是指去噪声后的雷达功率谱。确定单元,用于确定所述第一功率谱是否为双峰谱,其中,所述双峰谱中包含有两个极大值点。所述确定单元,还用于若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点。处理单元,用于将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱。所述第一功率谱包括所述第二功率谱及第三功率谱。所述获取单元,还用于在所述第一功率谱中获取所述第三功率谱。
本发明实施例提供了一种功率谱识别方法及装置,获取第一功率谱,第一功率谱是指去噪声后的雷达功率谱。确定第一功率谱是否为双峰谱,其中,所述双峰谱中包含有两个极大值点。若确定出第一功率谱为双峰谱,则确定出第一双峰谱的第一对称极值点。将第一功率谱以第一对称极值点为对称中心进行对称处理,获取第二功率谱。第一功率谱包括了所述第二功率谱及第三功率谱,进而可以在第一功率谱中获取第三功率谱。这样,在阴雨条件下,当风廓线雷达获取的功率谱为降水谱与湍流谱叠加的双峰谱时,通过上述方法可以准确的识别出此双峰谱中的降水谱及湍流谱,从而可以仅根据降水谱获取到相应的降水信息,仅根据湍流谱获取大气湍流信息,提高了获取大气湍流信息及降水信息的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种功率谱识别方法的流程示意图;
图2为本发明实施例提供的一种确定第一功率谱是否为双峰谱的流程示意图;
图3为本发明实施例提供的一种确定第一功率谱信号区间的流程示意图;
图4为本发明实施例提供的一种第一功率谱及其信号区间的示例图;
图5为本发明实施例提供的另一种第一功率谱及其信号区间的示例图;
图6为本发明实施例提供的一种功率谱识别装置的功能结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
另外,本文中术语“系统”和“网络”在本文中常被可互换使用。本文中术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。
本发明实施例提供了一种功率谱识别方法,如图1所示,包括:
101、获取第一功率谱。
其中,所述第一功率谱是指由雷达功率谱去掉噪声功率谱后得到的功率谱。
具体的,功率谱识别装置获取第一功率谱包括:获取雷达功率谱。对所述雷达功率谱进行去噪声处理,得到第一功率谱。
也就是说,功率谱识别装置通过散射回来的电磁波可以获取到雷达功率谱,在获取了雷达功率谱后,对此雷达功率谱进行去噪声处理,即为将雷达功率谱与预设的噪声门限初始值进行比较,由于低于噪声门限初始值的为噪声,所以对雷达功率谱中高于此噪声门限初始值的部分进行裁剪,将雷达功率谱中的噪声部分保留,从而可以获取到相应的噪声功率谱。
需要说明的是,从雷达功率谱中,获取相应的噪声功率谱的方法与现有技术中相同,在此不再赘述。
在获取到噪声功率谱后,可以将雷达功率谱中的噪声功率谱消除掉,进而获取第一功率谱。
需要说明的是,第一功率谱和噪声功率谱相互独立,且满足可加性,即为本发明实施例中的所述的雷达功率谱是由第一功率谱及噪声功率谱叠加而成。
需要说明的是,在本发明所有实施例中,噪声功率谱均为白噪声功率谱。
102、确定所述第一功率谱是否为双峰谱。
其中,所述双峰谱中包含有两个极大值点。
具体的,功率谱识别装置在获取到第一功率谱后,需要确定此功率谱是双峰谱还是单峰谱。
其中,确定第一功率谱是否为双峰谱的方法,如图2所示,包括以下几个步骤:
1021、确定所述第一功率谱对应的信号区间。
具体的,功率谱识别的装置在获取到第一功率谱后,可以根据第一功率谱中的信号点确定出其对应的信号区间,确定第一功率谱对应的信号区间有以下方法,如图3所示,包括:
10211、确定信号区间的起始端点。
具体的,确定信号区间的起始端点的方法如下:
确定所述第一FFT(Fast Fourier Transformation,快速傅氏变换)点;若确定出第一FFT点至其后连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第一FFT点确定为信号区间的起始端点。若确定出第一FFT点至其后连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第i个FFT点对应的第一功率谱的信号点确定为第一功率谱的起始信号点,并将第i个FFT点确定为第一FFT点,直至确定出第一FFT点至其后连续的n-1个对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零。
所述第一FFT点是第一功率谱的起始信号点对应的FFT点。所述第i个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第一FFT点之后的,距离所述第一FFT点最近的点;所述n为大于0的整数,q为预先设置的,大于0的整数,所述i=1,2,......,n;其中,p为大于n与q之差的整数。
也就是说,功率谱识别装置在第一功率谱的信号点对应的FFT点中,先确定出第一FFT点。此时,功率谱识别装置将第一功率谱的起始信号点即为第一个信号点确定为第一FFT点。在确定出第一FFT点后,确定第一FFT点及其后的连续n-1个FFT点对应的第一功率谱的信号点的功率中,若有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第一FFT点确定为第一功率谱对应的信号区间的起始端点。若第一FFT点及其后的连续n-1个FFT点对应的第一功率谱的信号点的功率中,有p个FFT点对应的第一功率谱的信号点的功率均不为零,则需要重新设置第一功率谱的起始端点,即为将第i个FFT点对应的信号点确定为第一功率谱的起始端点,并将第i个FFT点确定为第一FFT点。例如:确定出第一个信号点也对应的FFT点及其后的连续n-1个FFT点对应的第一功率谱的信号点的功率中,只有p个FFT点对应的第一功率谱的信号点的功率均不为零,那么重新确定第一FFT点,将第一个信号点接下来的最近的不为零的信号点也就是第二个信号点对应的FFT点作为第一FFT点,依次类推,直到确定出满足要求的第一FFT点及信号起始端点。
10212、确定信号区间的终止端点。
具体的,确定信号区间的终止端点的方法如下:
确定所述第二FFT点,若确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第二FFT点确定为信号区间的终止端点,所述第二FFT点是第一功率谱的终止信号点对应的FFT点。若确定出第二FFT点至其前连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第m个FFT点对应的第一功率谱的信号点确定为第一功率谱的终止信号点,并将第m个FFT点确定为第二FFT点,直至确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零。
所述第二FFT点是第一功率谱的终止信号点对应的FFT点。第m个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第二FFT点之前的,距离所述第二FFT点最近的点。
也就是说,功率谱识别装置先在第一功率谱的信号点对应的FFT点中,先确定出第一FFT点。此时,功率谱识别装置将第一功率谱的终止信号点即为最后一个信号点确定为第二FFT点。在确定出第二FFT点后,确定第二FFT点及其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,若有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第二FFT点确定为第一功率谱对应的信号区间的终止端点。若第二FFT点及其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有p个FFT点对应的第一功率谱的信号点的功率均不为零,则需要重新设置第一功率谱的终止端点,即为将第i个FFT点对应的信号点确定为第一功率谱的起始端点,并将第i个FFT点确定为第二FFT点。例如:确定出最后一个信号点对应的FFT点及其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,只有p个FFT点对应的第一功率谱的信号点的功率均不为零,那么重新确定第二FFT点,将最后一个信号点前边接下来的最近的一个不为零的信号点对应的FFT点作为第二FFT点,依次类推,直到确定出满足要求的第二FFT点及信号终止端点。
10213、根据信号区间的起始端点和终止端点确定所述信号区间。
具体的,功率谱识别装置在确定出信号区间的起始端点和终止端点之后,将信号区间的起始端点和终止端点之间的区间确定为所述第一功率谱的信号区间。
示例性的,假设功率谱识别装置确定出第一功率谱的信号点对应的FFT点的范围为1-256,如图4所示。其中,在此FFT点中,FFT点1-90对应的第一功率谱的信号点的功率值为0,FFT点162-256对应的第一功率谱的信号点的功率值也为0。此时,在图4中并没有表示出FFT点1-72及FFT点172-256对应的第一功率谱的信号点,仅表示出FFT点90-162对应的第一功率谱的信号点。
功率谱识别装置先在第一功率谱的信号点对应的FFT点中,先确定出第一FFT点和第二FFT点。
假设n为10,p等于4,q等于5。那么将图4所示FFT点是95的FFT点记为第一FFT点,表示信号区间的起始端点。FFT点是162的FFT点记为第二FFT点表示信号区间的终止端点。确定出95的FFT点和162的FFT点之间的横坐标区间为信号区间。
需要说明的是,通过公式FFT点对应的速度值=C*FFT点,其中,C表示最小速度间隔系数,可以将上述横坐标FFT点转换成如图5所示的FFT点对应的速度值的横坐标形式,在这种横坐标形式下,确定信号区间的方法和横坐标是FFT点时的方法一致。
1022、在所述信号区间内,对所述第一功率谱进行插值及滑动平均处理,确定出第一功率谱的极大值点的个数。
需要说明的是,插值法是指在离散数据之间补充一些数据,使这组离散数据能够符合某个连续函数。所述滑动平均处理,可以使用为三点滑动平均、五点滑动平均等,本发明对此不做限制。
具体的,功率谱识别装置在确定出第一功率谱的信号区间后,在此信号区间内,确定出每相邻的两个功率值不为零的信号点间是否有功率值为零的信号点。若确定出相邻的两个功率值不为零的信号点间有功率值为零的信号点,则确定出此相邻的两个功率值不为零的信号点间功率值为零的信号点的个数。并根据确定出的此相邻的两个功率值不为零的信号点间功率值为零的信号点的个数,利用插值法,计算出此相邻的两个功率值不为零的信号点间需要插入的信号点;用计算出的需要插入的信号点替换此相邻的两个功率值不为零的信号点间的功率值为零的信号点,以使得相邻的两个功率值不为零的信号点间的信号点的功率均不为零。
示例性的,以插值法为线性插值法为例进行说明。假设第一功率谱中有功率值不为零的两个相邻信号点A点和B点,且A点和B点之间有两个功率值为零的点,这时可以根据A点和B点的功率值,通过线性插值法计算,获取需插入的两个信号点。此需插入的两个信号点间的功率值与A点、B点的功率值成等差关系,将计算出的需要插入的信号点替换A点和B点之间的功值为零的信号点,以使得A、B两个功率值不为零的信号点间的信号点的功率均不为零。
需要说明的是,插值法还可以是其他插值法,例如拉格朗日插值,最小二乘插值。根据插值法不同,则获取需要插入的信号点的方法不同,即为对第一功率谱进行插值处理的方法不同,本发明对此不做限制。
进一步的,功率谱识别装置可以检测出第一功率谱中极大值点的个数。
1023、根据所述第一功率谱极大值点的个数确定所述第一功率谱是否为双峰谱。
具体的,功率谱识别装置根据检测出的第一功率谱中极大值点的个数,确定第一功率谱是否为双峰谱。即当检测到所述第一功率谱中极大值点个数是一个时,确定所述第一功率谱是单峰谱。若极大值点个数是两个时,确定所述第一功率谱为双峰谱。
若确定第一功率谱为单峰谱,则可确定获取的第一功率谱即为所需的功率谱,从而可以根据此第一功率谱获取相应的信息。
103、若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点。
其中,所述第一对称极值点是用于获取第二功率谱的极大值点。
需要说明的是,第二功率谱及第三功率谱叠加构成第一功率谱。
具体的,功率谱识别装置确定出第一功率谱为双峰谱时,说明第一功率谱是由第二功率谱和第三功率谱叠加而成,此时需要将构成此第一功率谱的第二功率谱和第三功率谱识别出来。此时,功率谱识别装置可以先获取第一功率谱的第一对称极值点,进而可以根据此第一对称极值点确定出第二功率谱,及第三功率谱,从而实现将第二功率谱及第三功率谱识别出来。
进一步的,所述确定出第一对称极值点的方法如下:
确定所述第一功率谱的第一极大值点与所述信号区间的起始端点间的距离。确定所述第一功率谱的第二极大值点与所述信号区间的起始端点间的距离。
若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第一极大值点或第二极大值点确定为所述第一对称极值点。
其中,所述第一功率谱包括第一极大值点和第二极大值点。所述第一极大值点是所述第一功率谱中与所述信号区间的起始端点间的距离最小的极大值点。所述第二极大值点是所述第一功率谱中与所述信号区间的终止端点间的距离最小的极大值点。
需要说明的是,当利用第一功率谱的第一极大值点作对称获取第二功率谱时,第一极大值点是第一对称极值点;当利用第一功率谱的第二极大值点作对称获取第二功率谱时,第二极大值点是第一对称极值点。
也就是说,功率谱识别装置确定出第一功率谱的第一极大值点与信号区间的起始端点间的距离,及第二极大值点与信号区间的起始端点间的距离。此时,根据第一极大值点与信号区间的起始端点间的距离与第一预设值间关系,及第二极大值点与信号区间的起始端点间的距离与第二预设值间的关系不同,则确定出第一对称极值点的方法不同。若第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则第一对称极值点可以是第一极大值点,也可以是第二极大值点。
需要说明的是,第一预设值及第二预设值是预先设置的。
示例性的,假设第一预设值是信号区间起始端点与终止端点间的距离的5%,第二预设值是信号区间起始端点与终止端点间的距离的95%。若第一极大值点与所述信号区间的起始端点间的距离大于等于信号区间起始端点与终止端点间的距离的5%,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于信号区间起始端点与终止端点间的距离的95%,则将所述第一极大值点或第二极大值点确定为所述第一对称极值点。
进一步的,若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将所述第一极大值点确定为所述第一对称极值点。
也就是说,功率谱识别装置确定出第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,说明第二极大值点接近于信号区间的终止端点,此时若将第二极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第二极大值点接近于信号区间的终止端点,使得对称获取的功率谱中的信号点过少,无法有效利用。因此,第二极大值点并不能作为第一对称极值点,将第一极大值点作为第一对称极值点。
示例性的,假设第一预设值是信号区间起始端点与终止端点间的距离的5%,第二预设值是信号区间起始端点与终止端点间的距离的95%。若第一极大值点与所述信号区间的起始端点间的距离大于等于信号区间起始端点与终止端点间的距离的5%,且所述第二极大值点与所述信号区间的起始端点间的距离大于信号区间起始端点与终止端点间的距离的95%,则将所述第一极大值点确定为所述第一对称极值点。
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第二极大值点确定为所述第一对称极值点;
也就是说,当功率谱识别装置确定出第一极大值点与信号区间的起始端点间的距离小于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,说明第一极大值点接近于信号区间的起始端点,此时若将第一极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第一极大值点接近于信号区间的起始端点,使得对称获取的功率谱中的信号点过少,无法有效利用。因此,第一极大值点并不能作为第一对称极值点,将第二极大值点作为第一对称极值点。
示例性的,假设第一预设值是信号区间起始端点与终止端点间的距离的5%,第二预设值是信号区间起始端点与终止端点间的距离的95%。若第一极大值点与所述信号区间的起始端点间的距离小于信号区间起始端点与终止端点间的距离的5%,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于信号区间起始端点与终止端点间的距离的95%,则将所述第二极大值点确定为所述第一对称极值点。
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将功率值最大对应的极值点确定为所述第一对称极值点。
也就是说,功率谱识别装置确定出第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,说明第一极大值点接近于信号区间的起始端点,且第二极大值点接近于信号区间的终止端点,此时确定信号点对应功率值最大对应的极值点为所述第一极大值即当第一极大值点是区间范围内的最大值点时,确定第一极大值点为第一对称极值点,当第二极大值点是区间范围内的最大值点时,确定第二极大值点为第一对称极值点。需要说明的是,在实际功率谱中,当功率谱为双峰谱时,第一极大值点接近于信号区间的起始端点的,同时第二极大值点接近于信号区间的终止端点的情况几乎不存在。
示例性的,假设第一预设值是信号区间起始端点与终止端点间的距离的5%,第二预设值是信号区间起始端点与终止端点间的距离的95%。若第一极大值点与所述信号区间的起始端点间的距离信号区间起始端点与终止端点间的距离的5%,且所述第二极大值点与所述信号区间的起始端点间的距离大于信号区间起始端点与终止端点间的距离的95%,若第一极值点为信号区间内的最大功率值,则确定第一极大值点为第一对称极值点,若第二极值点为信号区间内的最大功率值,则确定第二极大值点为第一对称极值点。
需要说明的是,第一功率谱可以是由湍流功率谱及降水功率谱叠加而成,即为第二功率谱可以是湍流功率谱,此时,第三功率谱可以是降水功率谱,第二功率谱也可是降水谱,此时第三功率谱是湍流功率谱。当然,第一功率谱可以是由其他功率谱叠加而成,这是第二功率谱及第三功率谱是其他叠加成第一功率谱的功率谱,本发明对此不做限制。
104、将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱。
其中,所述第一功率谱包括所述第二功率谱及第三功率谱。
具体的,若将所述功率谱的第一极大值点作为第一对称极值点,则以第一极大值点为对称中心,以信号起始端点为起点,从区间起始端点到终止端点方向作对称,对称后获取的功率谱为第二功率谱。
此时,若第一功率谱由湍流功率谱及降水功率谱叠加而成时,第二功率谱为湍流功率谱。
若将所述功率谱的第二极大值点作为第一对称极值点,则以第二极大值点为对称中心,以信号终止端点为起点,从区间终止端点到起始端点方向作对称,对称后获取的功率谱为第二功率谱。
此时,若第一功率谱由湍流功率谱及降水功率谱叠加而成时,第二功率谱为降水功率谱。
105、在所述第一功率谱中获取所述第三功率谱。
具体的,根据第一极大值点与信号区间的起始端点间的距离与第一预设值间的关系,及第二极大值点与所述信号区间的起始端点间的距离与第二预设值关系不同,在所述第一功率谱中获取所述第三功率谱的方法不同。具体如下:
在第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,将所述第一功率谱以所述第一功率谱的第二对称极值点为对称中心进行对称处理,获取所述第三功率谱。
其中,所述第二对称极值点是所述第一功率谱中除所述第一对称极值点之外的极值点。
也就是说,当第一极大值点与区间的区间起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,第一对称极值点可以是第一极大值点或者第二极大值点。功率谱识别装置在将第一极大值点作为第一对称极值点时,第二极大值点是第二对称极值点,功率谱识别装置以第二对称极值点为对称中心进行对称处理,获取的功率谱为第三功率谱。功率谱识别装置在将第二极大值点作为第一对称极值点时,第一极大值点是第二对称极值点,功率谱识别装置以第二对称极值点为对称中心进行对称处理,获取的功率谱为第三功率谱。
或者,所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,或所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
即为,当第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
也就是说,当第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的区间端点间的距离大于第二预设值时,说明第二极大值点接近于信号区间的终止端点,此时若将第二极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第二极大值点接近于信号区间的终止端点,使得对称获取的功率谱中的信号点过少,无法有效利用,所以将第一极大值点确定为所述第一对称极值点。功率谱识别装置根据第一对称极值点可以获取到第二功率谱。在获取了第二功率谱后,可以根据第一功率谱和第二功率谱,利用差值法,获取第三功率谱,即为在所述第一功率谱中将第二功率谱消除掉,从而获取第三功率谱。
当第一极大值点与信号区间的起始端点间的距离小于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
也就是说,当第一极大值点与信号区间的起始端点间的距离小于第一预设值,且第二极大值点与所述信号区间的区间端点间的距离小于等于第二预设值时,说明第一极大值点接近于信号区间的起始端点,此时若将第一极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第一极大值点接近于信号区间的起始端点,使得对称获取的功率谱中的信号点过少,无法有效利用,所以将第二极大值点确定为所述第一对称极值点。功率谱识别装置根据第一对称极值点可以获取到第二功率谱。在获取了第二功率谱后,可以根据第一功率谱和第二功率谱,利用差值法,获取第三功率谱,即为在所述第一功率谱中将第二功率谱消除掉,从而获取第三功率谱。
当第一极大值点与信号区间的起始端点间的距离小于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,将信号功率值最大值点作为第一对称极值点,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
也就是说,当第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的区间端点间的距离大于第二预设值时,将信号功率值最大值点作为第一对称极值点。若功率值最大值点为第二极大值点,则将第二极大值点作为第一对称极值点,功率谱识别装置根据第一对称极值点可以获取到第二功率谱。在获取了第二功率谱后,可以根据第一功率谱和第二功率谱,利用差值法,获取第三功率谱,即为在所述第一功率谱中将第二功率谱消除掉,从而获取第三功率谱。若功率值最大值点为第一极大值点,则将第一极大值点作为第一对称极值点,功率谱识别装置根据第一对称极值点可以获取到第二功率谱。在获取了第二功率谱后,可以根据第一功率谱和第二功率谱,利用差值法,获取第三功率谱,即为在所述第一功率谱中将第二功率谱消除掉,从而获取第三功率谱。
需要说明的是,若第一功率谱由湍流功率谱及降水功率谱叠加而成,且第二功率谱为湍流功率谱时,则第三功率谱为降水功率谱。若第一功率谱由湍流功率谱及降水功率谱叠加而成,且第二功率谱为降水功率谱时,则第三功率谱为湍流功率谱。
本发明实施例提供了一种功率谱的分离方法,获取第一功率谱,第一功率谱是指无噪声功率谱。确定第一功率谱是否为双峰谱,其中,所述双峰谱中包含有两个极大值点。若确定出第一功率谱为双峰谱,则确定出第一双峰谱的第一对称极值点。将第一功率谱以第一对称极值点为对称中心进行对称处理,获取第二功率谱。第一功率谱包括了所述第二功率谱及第三功率谱,进而可以在第一功率谱中获取第三功率谱。这样,在阴雨条件下,在风廓线雷达获取的功率谱为降水谱与湍流谱叠加的双峰谱时,通过上述方法可以准确的识别出此双峰谱中的降水谱及湍流谱,从而可以仅根据降水谱获取到相应的降水信息,仅根据湍流谱获取大气湍流信息,提高了获取大气湍流信息及降水信息的准确性。
本发明实施例提供了一种功率谱识别装置,如图6所示,包括:
获取单元601,用于获取第一功率谱。
其中,所述第一功率谱是指指由雷达功率谱去掉噪声功率谱后得到的功率谱。
具体的,所述获取单元601,具体用于获取雷达功率谱,并对所述雷达功率谱进行去噪声处理,得到第一功率谱。
也就是说,获取单元601通过散射回的电磁波信号可以获取到雷达功率谱,在获取了雷达功率谱后,对此雷达功率谱进行去噪声处理,即为将雷达功率谱与预设的噪声门限初始值进行比较,由于低于噪声门限初始值的信号为噪声信号,所以对雷达功率谱中高于此噪声门限初始值的部分进行裁剪,将雷达功率谱中的噪声信号部分保留,从而可以获取到相应的噪声功率谱。
需要说明的是,第一功率谱和噪声功率谱相互独立,且满足可加性,即为本发明实施例中的所述的雷达功率谱是由第一功率谱及噪声功率谱叠加而成。
确定单元602,用于确定所述第一功率谱是否为双峰谱。
其中,所述双峰谱中包含有两个极大值点。
具体的,确定单元602需要在获取单元601获取到第一功率谱后,确定此功率谱是双峰谱还是单峰谱信号。
进一步的,确定单元602,具体用于确定所述第一功率谱对应的信号区间;在所述信号区间内,对所述第一功率谱进行插值及滑动平均处理;确定出插值及滑动平均处理后的第一功率谱对应的极大值点的个数;根据所述极大值点的个数确定所述插值及滑动平均处理后的第一功率谱是否为双峰谱。
其中,确定单元602用于确定所述第一功率谱对应的信号区间的具体方式如下:
所述确定单元602,具体用于确定所述第一FFT点,若确定出第一FFT点至其后连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第一FFT点确定为信号区间的起始端点。若确定出第一FFT点至其后连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第i个FFT点对应的第一功率谱的信号点确定为第一功率谱的起始信号点,并将第i个FFT点确定为第一FFT点,直至确定出第一FFT点至其后连续的n-1个对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零。
确定所述第二FFT点,若确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第二FFT点确定为信号区间的终止端点所述第二FFT点是第一功率谱的终止信号点对应的FFT点。若确定出第二FFT点至其前连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第m个FFT点对应的第一功率谱的信号点确定为第一功率谱的终止信号点,并将第m个FFT点确定为第二FFT点,直至确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零。
其中,所述第一FFT点是第一功率谱的起始信号点对应的FFT点。所述第i个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第一FFT点之后的,距离所述第一FFT点最近的点。所述n为大于0的整数,q为预先设置的,大于0的整数,所述i=1,2,......,n;其中,p为大于n与q之差的整数。所述第二FFT点是第一功率谱的终止信号点对应的FFT点。第m个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第二FFT点之前的,距离所述第二FFT点最近的点。
所述确定单元602用于在所述信号区间内,对所述第一功率谱进行插值及滑动平均处理的具体方式如下:在确定出第一功率谱的信号区间后,在此信号区间内,确定出每相邻的两个功率值不为零的信号点间是否有功率值为零的信号点。若确定出相邻的两个功率值不为零的信号点间有功率值为零的信号点,则确定出此相邻的两个功率值不为零的信号点间功率值为零的信号点的个数。并根据确定出的此相邻的两个功率值不为零的信号点间功率值为零的信号点的个数,利用插值法,计算出此相邻的两个功率值不为零的信号点间需要插入的信号点;用计算出的需要插入的信号点替换此相邻的两个功率值不为零的信号点间的功率值为零的信号点,以使得相邻的两个功率值不为零的信号点间的信号点的功率均不为零。确定单元602再对插值后的信号区间范围内的第一功率谱点进行滑动平均处理。
需要说明的是,插值法是指在离散数据之间补充一些数据,使这组离散数据能够符合某个连续函数。所述滑动平均处理,可以使用为三点滑动平均、五点滑动平均等,本发明对此不做限制。
所述确定单元602用于根据所述第一功率谱极大值点的个数确定所述第一功率谱是否为双峰谱的具体方式如下:
具体的,确定单元602可以检测出第一功率谱中极大值点的个数,当检测到所述第一功率谱中极大值点个数是一个时,确定所述第一功率谱是单峰谱信号。若极大值点个数是两个时,确定所述第一功率谱为双峰谱。
进一步的,若确定单元602确定第一功率谱为单峰谱信号,则可确定获取的第一功率谱即为所需的功率谱,从而可以根据此第一功率谱获取相应的信息。
所述确定单元602,还用于若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点。
其中,所述第一对称极值点是用于获取第二功率谱的极大值点。
需要说明的是,第二功率谱及第三功率谱叠加构成第一功率谱。
具体的,确定单元602确定出第一功率谱为双峰谱时,说明第一功率谱是由第二功率谱及第三功率谱叠加而成,此时需要将构成此第一功率谱的第二功率谱及第三功率谱识别出来。此时,确定单元602可以先确定出插值及滑动平均处理后的第一功率谱第一对称极值点,进而可以根据此第一对称极值点确定出第二功率谱,及第三功率谱,从而实现将第二功率谱及第三功率谱识别出来。
其中,所述确定单元602用于确定出第一对称极值点的具体方式如下:
所述确定单元602,具体用于确定所述第一功率谱的第一极大值点与所述信号区间的起始端点间的距离。并确定所述第一功率谱的第二极大值点与所述信号区间的起始端点间的距离。
其中,所述第一功率谱包括第一极大值点和第二极大值点。所述第一极大值点是所述第一功率谱中与所述信号区间的起始端点间的距离最小的极大值点。所述第二极大值点是所述第一功率谱中与所述信号区间的终止端点间的距离最小的极大值点。
若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第一极大值点或第二极大值点确定为所述第一对称极值点。
也就是说,当确定单元602确定出第一功率谱的极大值点与信号区间的起始端点间的最小距离,及极大值点与信号区间的终止端点间的最小距离,即为确定出第一功率谱的第一极大值点与信号区间的起始端点间的距离,及第二极大值点与信号区间的起始端点间的距离。此时,根据第一极大值点与信号区间的起始端点间的距离与第一预设值间关系,及第二极大值点与信号区间的起始端点间的距离与第二预设值间的关系不同,则确定出第一对称极值点的方法不同。若第一极大值点与所述信号区间的区间端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的终止端点间的距离小于等于第二预设值,则第一对称极值点可以是第一极大值点,也可以是第二极大值点。
需要说明的是,第一预设值及第二预设值是预先设置的。
或者,若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将所述第一极大值点确定为所述第一对称极值点。
也就是说,当确定单元602确定出第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,说明第二极大值点接近于信号区间的终止端点,此时若将第二极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第二极大值点接近于信号区间的终止端点,使得对称获取的功率谱中的信号点过少,无法有效利用。因此,第二极大值点并不能作为第一对称极值点,将第一极大值点作为第一对称极值点。
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第二极大值点确定为所述第一对称极值点。
也就是说,当确定单元602确定出第一极大值点与信号区间的起始端点间的距离小于第一预设值,且第二极大值点与所述信号区间的起始端点间的小于等于第二预设值时,说明第一极大值点接近于信号区间的起始端点,此时若将第一极大值点作为第一对称极值点时,在以第一对称极值点为中心进行对称处理时,由于第一极大值点接近于信号区间的起始端点,使得对称获取的功率谱中的信号点过少,无法有效利用。因此,第一极大值点并不能作为第一对称极值点,将第二极大值点作为第一对称极值点。
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将功率值最大对应的极值点确定为所述第一对称极值点。
也就是说,当确定单元602确定出第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,说明第一极大值点接近于信号区间的起始端点,且第二极大值点接近于信号区间的终止端点,此时确定信号点对应功率值最大对应的极值点为所述第一极大值即当第一极大值点是区间范围内的最大值点时,确定第一极大值点为第一对称极值点,当第二极大值点是区间范围内的最大值点时,确定第二极大值点为第一对称极值点。需要说明的是,在实际功率谱中,当功率谱为双峰谱时,第一极大值点接近于信号区间的起始端点的,同时第二极大值点接近于信号区间的终止端点的情况几乎不存在。
需要说明的是,第一功率谱可以是由湍流功率谱及降水功率谱叠加而成,即为第二功率谱可以是湍流功率谱,此时,第三功率谱可以是降水功率谱,第二功率谱也可是降水谱,此时第三功率谱是湍流功率谱。当然,第一功率谱可以是由其他功率谱叠加而成,这是第二功率谱及第三功率谱是其他叠加成第一功率谱的功率谱,本发明对此不做限制。
处理单元603,用于将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱。
其中,所述第一功率谱包括所述第二功率谱及第三功率谱。
具体的,若确定单元602将所述第一功率谱的第一极大值点作为第一对称极值点,则处理单元603以第一极大值点为对称中心,以信号起始端点为起点,从信号区间起始端点到终止端点方向作对称,对称后获取的功率谱为第二功率谱。
此时,若第一功率谱由湍流功率谱及降水功率谱叠加而成时,第二功率谱为湍流功率谱。
若确定单元602将所述第一功率谱的第二极大值点作为第一对称极值点,则处理单元603以第二极大值点为对称中心,以信号终止端点为起点,从信号区间终止端点到起始端点方向作对称,对称后获取的功率谱为第二功率谱。
此时,若第一功率谱由湍流功率谱及降水功率谱叠加而成时,第二功率谱为降水功率谱。
所述获取单元601,还用于在所述第一功率谱中获取所述第三功率谱。
具体的,获取单元601根据第一极大值点与信号区间的起始端点间的距离与第一预设值间的关系,及第二极大值点与所述信号区间的起始端点间的距离与第二预设值关系不同,在所述第一功率谱中获取所述第三功率谱的方法不同。具体如下:
所述获取单元,具体用于在第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,将所述第一功率谱以所述第一功率谱的第二对称极值点为对称中心进行对称处理,获取所述第三功率谱。
其中,所述第二对称极值点是所述第一功率谱中除所述第一对称极值点之外的极值点。
也就是说,当第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,第一对称极值点可以是第一极大值点或者第二极大值点。确定单元602在将第一极大值点确定为第一对称极值点时,那么第二极大值点确定为第二对称极值点,获取单元601以第二对称极值点为对称中心进行对称处理,获取的功率谱为第三功率谱。确定单元602在将第二极大值点确定为第一对称极值点时,那么第一极大值点是第二对称极值点,获取单元601以第二对称极值点为对称中心进行对称处理,获取的功率谱为第三功率谱。
或者,所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,或所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
即是,当所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
也就是说,获取单元601在第一极大值点与信号区间的起始端点间的距离大于等于第一预设值,且第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
当第一极大值点与信号区间的起始端点间的距离小第一预设值,且第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
当第一极大值点与信号区间的起始端点间的距离小第一预设值,且第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,将功率值最大值点作为第一对称极值点,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
需要说明的是,若第一功率谱由湍流功率谱及降水功率谱叠加而成,且第二功率谱为湍流功率谱时,则第三功率谱为降水功率谱。若第一功率谱由湍流功率谱及降水功率谱叠加而成,且第二功率谱为降水功率谱时,则第三功率谱为湍流功率谱。
本发明实施例提供了一种功率谱识别方法及装置,获取第一功率谱,第一功率谱是指去噪声后的雷达功率谱。确定第一功率谱是否为双峰谱,其中,所述双峰谱中包含有两个极大值点。若确定出第一功率谱为双峰谱,则确定出第一双峰谱的第一对称极值点。将第一功率谱以第一对称极值点为对称中心进行对称处理,获取第二功率谱。第一功率谱包括了所述第二功率谱及第三功率谱,进而可以在第一功率谱中获取第三功率谱。这样,在阴雨条件下,在风廓线雷达获取的功率谱为降水谱与湍流谱叠加的双峰谱时,通过上述方法可以准确的识别出此双峰谱中的降水谱及湍流谱,从而可以仅根据降水谱获取到相应的降水信息,仅根据湍流谱获取大气湍流信息,提高了获取大气湍流信息及降水信息的准确性。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种功率谱识别方法,其特征在于,包括:
获取第一功率谱;所述第一功率谱是指由雷达功率谱去掉噪声功率谱后得到的功率谱;
确定所述第一功率谱是否为双峰谱;其中,所述双峰谱中包含有两个极大值点;
若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点;
将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱;所述第一功率谱包括所述第二功率谱及第三功率谱;
在所述第一功率谱中获取所述第三功率谱。
2.根据权利要求1所述的方法,其特征在于,所述确定所述第一功率谱是否为双峰谱包括:
确定所述第一功率谱对应的信号区间;
在所述信号区间内,对所述第一功率谱进行插值及滑动平均处理,
确定出第一功率谱的极大值点的个数;
根据所述第一功率谱极大值点的个数确定所述第一功率谱是否为双峰谱。
3.根据权利要求2所述的方法,其特征在于,所述确定出第一对称极值点包括:
确定所述第一功率谱的第一极大值点与所述信号区间的起始端点间的距离;所述第一功率谱包括第一极大值点和第二极大值点;所述第一极大值点是所述第一功率谱中与所述信号区间的起始端点间的距离最小的极大值点;
确定所述第一功率谱的第二极大值点与所述信号区间的起始端点间的距离;所述第二极大值点是所述第一功率谱中与所述信号区间的终止端点间的距离最小的极大值点;
若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第一极大值点或第二极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将所述第一极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第二极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将功率值最大对应的极值点确定为所述第一对称极值点;
所述第一功率谱中获取所述第三功率谱包括:
在所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,将所述第一功率谱以所述第一功率谱的第二对称极值点为对称中心进行对称处理,获取所述第三功率谱;其中,所述第二对称极值点是所述第一功率谱中除所述第一对称极值点之外的极值点;
或者,所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,或所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
4.根据权利要求2或3所述的方法,其特征在于,所述确定所述第一功率谱对应的信号区间包括:
确定所述第一快速傅氏变换FFT点;所述第一FFT点是第一功率谱的起始信号点对应的FFT点;
若确定出第一FFT点至其后连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第一FFT点确定为信号区间的起始端点;所述n为大于0的整数;q为预先设置的,大于0的整数;
若确定出第一FFT点至其后连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第i个FFT点对应的第一功率谱的信号点确定为第一功率谱的起始端点,并将第i个FFT点确定为第一FFT点,直至确定出第一FFT点至其后连续的n-1个对应的第一功率谱的点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零;所述i=1,2,……,n;其中,p为大于n与q之差的整数;第i个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第一FFT点之后的,距离所述第一FFT点最近的点;
确定所述第二FFT点;所述第二FFT点是第一功率谱的终止信号点对应的FFT点;
若确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第二FFT点确定为信号区间的终止端点;
若确定出第二FFT点至其前连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第m个FFT点对应的第一功率谱的信号点确定为第一功率谱的终止信号点,并将第m个FFT点确定为第二FFT点,直至确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零;所述m=1,2,……,n;第m个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第二FFT点之前的,距离所述第二FFT点最近的点。
5.根据权利要求1-4任一项所述的方法,其特征在于,
所述获取第一功率谱包括:
获取雷达功率谱;
对所述雷达功率谱进行去噪声处理,得到第一功率谱。
6.一种功率谱识别装置,其特征在于,包括:
获取单元,用于获取第一功率谱;所述第一功率谱是指由雷达功率谱去掉噪声功率谱后得到的功率谱;
确定单元,用于确定所述第一功率谱是否为双峰谱;其中,所述双峰谱中包含有两个极大值点;
所述确定单元,还用于若确定出所述第一功率谱为双峰谱,则确定出所述第一功率谱的第一对称极值点;
处理单元,用于将所述第一功率谱以所述第一功率谱的第一对称极值点为对称中心进行对称处理,获取第二功率谱;所述第一功率谱包括所述第二功率谱及第三功率谱;
所述获取单元,还用于在所述第一功率谱中获取所述第三功率谱。
7.根据权利要求6所述的装置,其特征在于,
所述确定单元,具体用于确定所述第一功率谱对应的信号区间;
在所述信号区间内,对所述第一功率谱进行插值及滑动平均处理,
确定第一功率谱的极大值点的个数;
根据所述第一功率谱极大值点的个数确定所述第一功率谱是否为双峰谱。
8.根据权利要求7所述的装置,其特征在于,
所述确定单元,具体用于确定所述第一功率谱的第一极大值点与所述信号区间的起始端点间的距离;所述第一功率谱包括第一极大值点和第二极大值点;所述第一极大值点是所述第一功率谱中与所述信号区间的起始端点间的距离最小的极大值点;
确定所述第一功率谱的第二极大值点与所述信号区间的起始端点间的距离;所述第二极大值点是所述第一功率谱中与所述信号区间的终止端点间的距离最小的极大值点;
若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第一极大值点或第二极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将所述第一极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值,则将所述第二极大值点确定为所述第一对称极值点;
或者,若所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值,则将功率值最大对应的极值点确定为所述第一对称极值点;
所述获取单元,具体用于在所述第一极大值点与所述信号区间的起始端点间的距离大于等于第一预设值,且所述第二极大值点与所述信号区间的起始端点间的距离小于等于第二预设值时,将所述第一功率谱以所述第一功率谱的第二对称极值点为对称中心进行对称处理,获取所述第三功率谱;其中,所述第二对称极值点是所述第一功率谱中除所述第一对称极值点之外的极值点;
或者,所述第一极大值点与所述信号区间的起始端点间的距离小于第一预设值,或所述第二极大值点与所述信号区间的起始端点间的距离大于第二预设值时,根据所述第二功率谱,在所述第一功率谱中获取所述第三功率谱。
9.根据权利要求7或8所述的装置,其特征在于,
所述确定单元,具体用于确定所述第一FFT点;所述第一FFT点是第一功率谱的起始信号点对应的FFT点;
若确定出第一快速傅氏变换FFT点至其后连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第一FFT点确定为信号区间的起始端点;所述n为大于0的整数;q为预先设置的,大于0的整数;
若确定出第一FFT点至其后连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第i个FFT点对应的第一功率谱的信号点确定为第一功率谱的起始信号点,并将第i个FFT点确定为第一FFT点,直至确定出第一FFT点至其后连续的n-1个对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零;所述i=1,2,……,n;其中,p为大于n与q之差的整数;第i个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第一FFT点之后的,距离所述第一FFT点最近的点;
确定所述第二FFT点;所述第二FFT点是第一功率谱的终止信号点对应的FFT点;
若确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零,则将第二FFT点确定为信号区间的终止端点;
若确定出第二FFT点至其前连续的n-1个FFT点中,有p个FFT点对应的第一功率谱的信号点的功率为零,则将第m个FFT点对应的第一功率谱的信号点确定为第一功率谱的终止信号点,并将第m个FFT点确定为第二FFT点,直至确定出第二FFT点至其前连续的n-1个FFT点对应的第一功率谱的信号点的功率中,有q个FFT点对应的第一功率谱的信号点的功率均不为零;所述m=1,2,……,n;第m个FFT点是对应的第一功率谱的信号点的功率不为零的,所述第二FFT点之前的,距离所述第二FFT点最近的点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466622.3A CN104316929B (zh) | 2014-09-12 | 2014-09-12 | 一种功率谱识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466622.3A CN104316929B (zh) | 2014-09-12 | 2014-09-12 | 一种功率谱识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104316929A true CN104316929A (zh) | 2015-01-28 |
CN104316929B CN104316929B (zh) | 2017-02-15 |
Family
ID=52372182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410466622.3A Expired - Fee Related CN104316929B (zh) | 2014-09-12 | 2014-09-12 | 一种功率谱识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316929B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614579A (zh) * | 2018-12-24 | 2019-04-12 | 雷象科技(北京)有限公司 | 相控阵天气雷达最大值提取信号方法 |
CN111505569A (zh) * | 2020-05-20 | 2020-08-07 | 浙江大华技术股份有限公司 | 声源定位方法以及相关设备、装置 |
CN112526548A (zh) * | 2020-11-18 | 2021-03-19 | 董晶晶 | 一种基于测风激光雷达的降雨识别方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4005424A (en) * | 1975-10-24 | 1977-01-25 | The United States Of America As Represented By The Secretary Of The Navy | Simulator for doppler radar returns from precipitation |
CN102426352A (zh) * | 2011-10-17 | 2012-04-25 | 中国人民解放军理工大学气象学院 | 风廓线雷达风计算方法 |
CN102508219A (zh) * | 2011-10-17 | 2012-06-20 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测方法 |
CN202433521U (zh) * | 2011-10-17 | 2012-09-12 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测处理板 |
-
2014
- 2014-09-12 CN CN201410466622.3A patent/CN104316929B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4005424A (en) * | 1975-10-24 | 1977-01-25 | The United States Of America As Represented By The Secretary Of The Navy | Simulator for doppler radar returns from precipitation |
CN102426352A (zh) * | 2011-10-17 | 2012-04-25 | 中国人民解放军理工大学气象学院 | 风廓线雷达风计算方法 |
CN102508219A (zh) * | 2011-10-17 | 2012-06-20 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测方法 |
CN202433521U (zh) * | 2011-10-17 | 2012-09-12 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测处理板 |
Non-Patent Citations (3)
Title |
---|
何平等: "风廓线雷达探测降水过程的初步研究", 《应用气象学报》 * |
何越等: "双高斯拟合的风廓线雷达降水时的风场反演", 《科学技术与工程》 * |
何越等: "基于双高斯拟合的风廓线雷达反演雨滴谱", 《创新驱动发展 提高气象灾害防御能力——S1第五届气象综合探测技术研讨会》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614579A (zh) * | 2018-12-24 | 2019-04-12 | 雷象科技(北京)有限公司 | 相控阵天气雷达最大值提取信号方法 |
CN109614579B (zh) * | 2018-12-24 | 2022-05-17 | 雷象科技(北京)有限公司 | 相控阵天气雷达最大值提取信号方法 |
CN111505569A (zh) * | 2020-05-20 | 2020-08-07 | 浙江大华技术股份有限公司 | 声源定位方法以及相关设备、装置 |
CN111505569B (zh) * | 2020-05-20 | 2022-04-19 | 浙江大华技术股份有限公司 | 声源定位方法以及相关设备、装置 |
CN112526548A (zh) * | 2020-11-18 | 2021-03-19 | 董晶晶 | 一种基于测风激光雷达的降雨识别方法及装置 |
CN112526548B (zh) * | 2020-11-18 | 2024-02-23 | 南京泰爱信科技有限公司 | 一种基于测风激光雷达的降雨识别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104316929B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931949B (zh) | 雷达扫描模式下风轮机杂波背景下飞机目标检测方法 | |
CN103714541A (zh) | 一种利用山体轮廓区域约束识别定位建筑物的方法 | |
CN112881993B (zh) | 自动识别雷达分布杂波引起虚假航迹的方法 | |
CN103592620B (zh) | 长基线时差定位系统中高重频信号解定位模糊方法 | |
CN102636775B (zh) | 基于模糊逻辑识别的风廓线雷达回波谱重构方法 | |
CN108562885B (zh) | 一种高压输电线路机载LiDAR点云提取方法 | |
CN113376648B (zh) | 基于激光雷达探测的高速非合作目标轨迹提取方法 | |
CN104035090A (zh) | 一种基于雷达目标检测的自由滑窗式cfar控制方法 | |
CN106443662B (zh) | 一种低重频体制下速度模糊时的目标稳健关联方法 | |
CN104316929A (zh) | 一种功率谱识别方法及装置 | |
CN105445731A (zh) | 一种基于高精度测速的雷达全自动跟踪方法 | |
CN112986964A (zh) | 基于噪声邻域密度的光子计数激光点云自适应去噪方法 | |
CN110208789A (zh) | 一种多帧联合目标航迹识别方法 | |
CN102879774B (zh) | 高频地波超视距雷达短航迹合成方法和装置 | |
CN110147716A (zh) | 一种基于频域与空域相结合的sar图像内波区域检测方法 | |
CN108107410B (zh) | 一种异型雷达联合判决目标检测方法 | |
CN109828275A (zh) | 一种零度层亮带的识别方法和系统 | |
CN106950527B (zh) | 一种多基线干涉仪测向体制下的脉冲信号分选方法 | |
KR101652193B1 (ko) | 레이더를 이용한 타겟 검출 방법 및 그 장치 | |
CN109283507B (zh) | 一种基于时频域特征的雷达目标识别方法及系统 | |
JP3127945B2 (ja) | 移動音源識別方法及び装置 | |
CN109239677B (zh) | 一种环境自适应恒虚警检测门限确定方法 | |
CN106371082B (zh) | 线性速度拖引干扰识别方法 | |
CN112083502B (zh) | 基于并联单稳随机共振的磁异信号检测方法 | |
CN115932831A (zh) | 基于雷达的目标分段跟踪方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 |