CN106125057B - 一种雷达回波移动矢量场处理方法 - Google Patents
一种雷达回波移动矢量场处理方法 Download PDFInfo
- Publication number
- CN106125057B CN106125057B CN201610452661.7A CN201610452661A CN106125057B CN 106125057 B CN106125057 B CN 106125057B CN 201610452661 A CN201610452661 A CN 201610452661A CN 106125057 B CN106125057 B CN 106125057B
- Authority
- CN
- China
- Prior art keywords
- mobile vector
- vector
- echo
- mobile
- radar return
- 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
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
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)
Abstract
本发明提供了一种雷达回波移动矢量场处理方法。其包括如下步骤:A、计算有移动矢量场的回波的主导移动矢量;B、从回波的移动矢量场中分解变形矢量;C、对所述变形矢量进行一致性检验;D、一致性检验符合预定的要求后,将所述变形矢量及主导移动矢量合成移动矢量。经过处理后的移动矢量场令所有雷达回波均能具有一定的移动矢量,而且移动矢量的分布均匀性较好,避免了凌乱移动矢量给回波外推预测带来的困难。
Description
技术领域
本发明涉及雷达回波技术领域,尤其涉及一种雷达回波移动矢量场处理方法。
背景技术
高质量的雷达回波移动矢量场是利用雷达回波开展准确的短时临近预报工作的基础。目前,有多种方法来进行移动矢量场的计算,主要包括:交叉相关法、变分法、TITAN技术等。但是,在进行雷达回波的移动矢量场计算过程中,会存在部分区域无有效移动矢量,部分移动矢量零散、一致性差等问题,从而导致雷达回波的短时临近预报准确性不高。
因此,现有技术还有待发展。
发明内容
鉴于上述现有技术的不足之处,本发明的目的在于提供一种雷达回波移动矢量场处理方法,旨在解决现有技术中移动矢量场中部分区域缺失、矢量场零散、一致性差,雷达回波的短时临近预报准确性不高的问题。
为了达到上述目的,本发明采取了以下技术方案:
一种雷达回波移动矢量场处理方法,其中,包括:
A、计算有移动矢量场的回波的主导移动矢量;
B、从回波的移动矢量场中分解变形矢量;
C、对所述变形矢量进行一致性检验;
D、一致性检验符合预定的要求后,将所述变形矢量及主导移动矢量合成移动矢量。
所述的雷达回波移动矢量场处理方法,其中,所述方法还包括:
采用主导移动矢量填补没有移动矢量场的回波;
根据步骤D计算获得的合成矢量结果进行回波外推预测。
所述的雷达回波移动矢量场处理方法,其中,所述分解变形矢量方法具体包括:
将回波块内各个点的移动矢量减去主导移动矢量,获得所述变形矢量。
所述的雷达回波移动矢量场处理方法,其中,所述步骤C具体包括:
C1、预设一致性系数的阈值;
C2、计算平均一致性系数,并且当所述平均一致性系数大于等于所述预设的阈值时,判断为符合一致性检验要求。
所述的雷达回波移动矢量场处理方法,其中,所述相邻两点之间的一致性系数如下所示:
ρim=Dim×Lim
其中,Dim为移动方向的一致性系数,Lim为移动速度的一致性系数。
所述的雷达回波移动矢量场处理方法,其中,所述Dim=cosθim;其中,θim为两移动矢量之间的方向差。
所述的雷达回波移动矢量场处理方法,其中,所述其中,和分别为相邻两点上的移动矢量。
所述的雷达回波移动矢量场处理方法,其中,所述步骤A具体包括:
以32个方位为标准来进行计算,选取回波块内各个点的主导移动矢量方向;
将回波块内出现频率最高的主导移动矢量方向作为回波块的主导移动矢量方向;
计算在选取的主导移动矢量方向上的各个点的移动距离平均值作为所述回波块的主导移动矢量的移动距离。
有益效果:本发明提供的一种雷达回波移动矢量场处理方法,经过处理后的移动矢量场令所有雷达回波均能具有一定的移动矢量,而且移动矢量的分布均匀性较好,避免了凌乱移动矢量给回波外推预测带来的困难。
附图说明
图1为本发明具体实施例的一种雷达回波移动矢量场处理方法的方法流程图。
图2a为采用现有方法的原始预报结果示意图。
图2b为使用本发明所述方法进行再处理的预报结果示意图。
图2c为实况结果示意图。
具体实施方式
本发明提供一种雷达回波移动矢量场处理方法。为使本发明的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,为本发明具体实施例的雷达回波移动矢量场处理方法。所述方法包括如下步骤:
S100、计算有移动矢量场的回波的主导移动矢量。具体的,在所述步骤S100之前,包括将回波分为有移动矢量场的回波和无移动矢量的回波。然后,采用主导移动矢量进行填补没有移动矢量场的区域。
S200、从回波的移动矢量场中分解变形矢量。
S300、对所述变形矢量进行一致性检验。
S400、一致性检验符合预定的要求后,将所述变形矢量及主导移动矢量合成移动矢量。
在本发明的具体实施例中,所述方法还包括:
S500、计算获得的合成矢量结果与所述没有移动矢量场的回波,从而获得最终用于进行回波外推预测的移动矢量场。
在实际应用过程中,可以将雷达回波的运动分为整体移动与内部变形两部分,其移动矢量可以分解为回波块移动矢量与变形矢量两部分,即移动矢量
其中,是指整个回波块的移动矢量,其方向是采用的回波块中出现频次最大的主导移动矢量方向。其移动距离值则选取的主导移动矢量方向上各个移动距离值的平均值。
具体的,回波块主导移动矢量方向可以以32个方位为标准来进行计算,选取出主导的移动矢量方位。
通过上述方法计算获得回波块移动矢量方向和移动距离后,可以据此计算变形矢量。所述变形矢量的计算方法具体为:将回波块内各个点的移动矢量减去主导移动矢量,获得所述变形矢量,用于表征回波块内部运动不均匀性。
在本实施例中,上述步骤S300的具体包括:预设一致性系数的阈值。当计算获得的结果大于阈值时,判断为符合一致性检验要求,可以认为移动矢量是合理的,不需要进行进一步处理。而当计算结果小于阈值时,需要进行进一步处理,改用回波块的主导移动矢量作为改点修正以后的移动矢量。
其计算过程具体如下所示:
首先,所述相邻两点i和m之间的一致性系数如下算式所示:
ρim=Dim×Lim (1)
其中,Dim为移动方向的一致性系数,Lim为移动速度的一致性系数。
算式(1)中Dim=cosθim;其中,θim为两移动矢量之间的方向差;算式(1)中其中,和分别为相邻两点i和m上的移动矢量。
为计算某一点移动矢量与周围移动矢量场的总一致性,以点i为中心,以R为影响半径,分别计算点i处的移动矢量与周围各个矢量的一致性系数,然后进行加权平均,从而得到总的平均一致性系数。由于两点越近,移动矢量的相关性越好,因此权重函数可以设为距离的函数:
rim是点i与点m间的距离,因而平均的一致性系数为
对于点i处的移动矢量,可以计算以它为中心以R为半径的水平范围内所有的移动矢量与它的一致性。
如上所述,判断计算所得的平均一致性系数与阈值ρ0(临界值)之间的大小。在本实施例中,所述阈值ρ0=0.5。
上述处理方法,通过计算主导移动矢量,变形矢量一致性检验以及最后的合成步骤,最终实现回波移动矢量场的再处理,能够使得所有的雷达回波内均有一定的移动矢量,并且移动矢量分布均匀性较好,避免了零乱移动矢量给回波外推预测带来的困难。
如图2a-2c所示,分别为使用现有方法的原始预报结果,使用上述方法进行再处理的预报结果以及实况结果。
可以理解的是,对本领域普通技术人员来说,可以根据本发明的技术方案及本发明构思加以等同替换或改变,而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.一种雷达回波移动矢量场处理方法,其特征在于,包括:
A、计算有移动矢量场的回波的主导移动矢量;
B、从回波的移动矢量场中分解变形矢量;
C、对所述变形矢量进行一致性检验;具体包括:C1、预设一致性系数的阈值;C2、计算平均一致性系数,并且当所述平均一致性系数大于等于所述预设一致性系数的阈值时,判断为符合一致性检验要求;
D、一致性检验要求符合后,将所述变形矢量及主导移动矢量合成移动矢量。
2.根据权利要求1所述的雷达回波移动矢量场处理方法,其特征在于,所述方法还包括:
采用主导移动矢量填补没有移动矢量场的回波;
根据步骤D计算获得的合成移动矢量结果与所述没有移动矢量场的回波,计算获得用于进行回波外推预测的移动矢量场。
3.根据权利要求1所述的雷达回波移动矢量场处理方法,其特征在于,所述分解变形矢量方法具体包括:
将回波块内各个点的移动矢量减去主导移动矢量,获得所述变形矢量。
4.根据权利要求3所述的雷达回波移动矢量场处理方法,其特征在于,相邻两点之间的一致性系数如下所示:
其中,为移动方向的一致性系数,为移动速度的一致性系数。
5.根据权利要求4所述的雷达回波移动矢量场处理方法,其特征在于,所述;其中,为两移动矢量之间的方向差。
6.根据权利要求5所述的雷达回波移动矢量场处理方法,其特征在于,所述;其中,和分别为相邻两点上的移动矢量。
7.根据权利要求1所述的雷达回波移动矢量场处理方法,其特征在于,所述步骤A具体包括:
以32个方位为标准来进行计算,选取回波块内各个点的主导移动矢量方向;
将回波块内出现频率最高的主导移动矢量方向作为回波块的主导移动矢量方向;
计算在选取的主导移动矢量方向上的各个点的移动距离平均值作为所述回波块的主导移动矢量的移动距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610452661.7A CN106125057B (zh) | 2016-06-20 | 2016-06-20 | 一种雷达回波移动矢量场处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610452661.7A CN106125057B (zh) | 2016-06-20 | 2016-06-20 | 一种雷达回波移动矢量场处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125057A CN106125057A (zh) | 2016-11-16 |
CN106125057B true CN106125057B (zh) | 2018-11-02 |
Family
ID=57470429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610452661.7A Active CN106125057B (zh) | 2016-06-20 | 2016-06-20 | 一种雷达回波移动矢量场处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125057B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538013A (zh) * | 2020-05-12 | 2020-08-14 | 上海眼控科技股份有限公司 | 雷达回波外推方法、装置、计算机设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002107462A (ja) * | 2000-09-28 | 2002-04-10 | Foundation Of River & Basin Integrated Communications Japan | 降雨洪水予測システム |
CN102662172A (zh) * | 2012-03-29 | 2012-09-12 | 天津大学 | 一种基于多普勒雷达反射率图像的风暴云团的外推方法 |
CN104849777A (zh) * | 2014-12-08 | 2015-08-19 | 国家电网公司 | 基于交叉相关追踪法的电网设备区域对流天气预报方法 |
CN104977584A (zh) * | 2015-06-29 | 2015-10-14 | 深圳市气象台 | 一种对流天气的临近预报方法及系统 |
-
2016
- 2016-06-20 CN CN201610452661.7A patent/CN106125057B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002107462A (ja) * | 2000-09-28 | 2002-04-10 | Foundation Of River & Basin Integrated Communications Japan | 降雨洪水予測システム |
CN102662172A (zh) * | 2012-03-29 | 2012-09-12 | 天津大学 | 一种基于多普勒雷达反射率图像的风暴云团的外推方法 |
CN104849777A (zh) * | 2014-12-08 | 2015-08-19 | 国家电网公司 | 基于交叉相关追踪法的电网设备区域对流天气预报方法 |
CN104977584A (zh) * | 2015-06-29 | 2015-10-14 | 深圳市气象台 | 一种对流天气的临近预报方法及系统 |
Non-Patent Citations (2)
Title |
---|
TREC算法在人工影响天气作业指挥系统中的应用;马建立;《气象科技》;20110430;第236-239页 * |
基于雷达回波外推和中尺度模式预报的短时降水对比分析;王丹;《高原气象》;20140630;第811-822页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106125057A (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108122234B (zh) | 卷积神经网络训练及视频处理方法、装置和电子设备 | |
US20210042929A1 (en) | Three-dimensional object detection method and system based on weighted channel features of a point cloud | |
CN107578430B (zh) | 一种基于自适应权值和局部熵的立体匹配方法 | |
CN106991650A (zh) | 一种图像去模糊的方法和装置 | |
WO2020125637A1 (zh) | 一种立体匹配方法、装置和电子设备 | |
CN111784588A (zh) | 图像数据增强方法、装置、计算机设备和存储介质 | |
CN111105452B (zh) | 基于双目视觉的高低分辨率融合立体匹配方法 | |
CN103218809A (zh) | 一种珍珠长度参数的图像测量方法 | |
CN106125057B (zh) | 一种雷达回波移动矢量场处理方法 | |
CN115601385A (zh) | 气泡形态处理方法、设备、介质及程序产品 | |
CN109903347B (zh) | 一种颜色混合的方法、系统、计算机设备及存储介质 | |
CN106323190B (zh) | 可定制深度测量范围的深度测量方法及深度图像的系统 | |
JP2020118536A (ja) | 風速分布推定装置及び風速分布推定方法 | |
CN114445473B (zh) | 基于深度学习算子的立体匹配方法和系统 | |
Lu et al. | Real-time stereo matching: A cross-based local approach | |
CN107424583B (zh) | 异形图像的显示数据处理方法和系统 | |
US20120121208A1 (en) | Image processing device and method | |
CN108717699B (zh) | 一种基于连续最小割的超声波图像分割方法 | |
CN111079523A (zh) | 物体检测方法、装置、计算机设备和存储介质 | |
CN106470342A (zh) | 全局运动估计方法及装置 | |
CN113096231B (zh) | 一种图像处理方法、装置、电子设备及存储介质 | |
CN107493471A (zh) | 一种视频传输质量的计算方法和装置 | |
CN110910436B (zh) | 基于图像信息增强技术的测距方法、装置、设备和介质 | |
CN104243947B (zh) | 视差估计方法和装置 | |
Höller et al. | Automatic object annotation in streamed and remotely explored large 3D reconstructions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |