CN107656256A - 利用雷达数据识别跟踪暴雨的方法 - Google Patents

利用雷达数据识别跟踪暴雨的方法 Download PDF

Info

Publication number
CN107656256A
CN107656256A CN201711061564.6A CN201711061564A CN107656256A CN 107656256 A CN107656256 A CN 107656256A CN 201711061564 A CN201711061564 A CN 201711061564A CN 107656256 A CN107656256 A CN 107656256A
Authority
CN
China
Prior art keywords
radar
radar map
map
heavy rain
lattice point
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
Application number
CN201711061564.6A
Other languages
English (en)
Other versions
CN107656256B (zh
Inventor
孙海元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou fengfenglei Technology Co., Ltd.
Original Assignee
Guangzhou Science And Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangzhou Science And Technology Co Ltd filed Critical Guangzhou Science And Technology Co Ltd
Priority to CN201711061564.6A priority Critical patent/CN107656256B/zh
Publication of CN107656256A publication Critical patent/CN107656256A/zh
Application granted granted Critical
Publication of CN107656256B publication Critical patent/CN107656256B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明涉及一种利用雷达数据识别跟踪暴雨的方法，依据第一预定时间的雷达格点数据获取第一雷达图；以及从第一雷达图中获取雷达回波强度大于预设阈值的若干第二雷达图；依据当前的雷达格点数据获取第三雷达图；以及从第三雷达图中获取雷达回波强度大于预设阈值的若干第四雷达图；获取各个第二雷达图和第四雷达图的质心点以及平均回波强度值；对各个第二雷达图和各个第四雷达图进行图形匹配，得到满足预定条件的若干组雷达图，依据各组雷达图中的第二雷达图和第四雷达图获取暴雨的移动速度和暴雨的移动方向；进而获取暴雨在第二预定时间覆盖的区域；对所述覆盖区域进行暴雨预警。

Description

利用雷达数据识别跟踪暴雨的方法
技术领域
本发明涉及气象预警技术领域，特别是涉及一种利用雷达数据识别跟踪暴雨的方法。
背景技术
突发的暴雨是主要气象灾害之一，气象雷达是监测暴雨最有效的手段，目前的天气预报都是大范围大时段的预警，不能直接预报(预警)未来30分钟内是否有雷暴雨。而企业和公众非常需要获取未来30分钟左右的精细降雨信息，以便采取适当的防护措施。比如，预报下午广州有暴雨，企业不可能整个下午停产。只有针对某家企业精确的预报或预警，企业据此停止危险左右，这样既保证了安全，又大幅减少了雷雨对正常生产的影响。
发明内容
本发明的目的是提供一种利用雷达数据识别跟踪暴雨的方法，可以提高雷雨预警的实时性和准确性。
一种利用雷达数据识别跟踪暴雨的方法，所述方法包括：
依据第一预定时间的雷达格点数据获取第一雷达图；以及从第一雷达图中获取雷达回波强度大于预设阈值的若干第二雷达图；
依据当前的雷达格点数据获取第三雷达图；以及从第三雷达图中获取雷达回波强度大于预设阈值的若干第四雷达图；
获取各个第二雷达图的质心点以及平均回波强度值；
获取各个第四雷达图的质心点以及平均回波强度值；
对各个第二雷达图和各个第四雷达图进行图形匹配，得到重叠度以及图形面积满足预定条件的若干组雷达图，每组雷达图包括一个第二雷达图和一个第四雷达图；
依据各组雷达图中的第二雷达图和第四雷达图的质心点和所述第一预定时间，获取暴雨的移动速度；
依据各组雷达图中的第二雷达图和第四雷达图的质心点获取暴雨的移动方向；
依据暴雨的移动速度和移动方向，获取暴雨在第二预定时间覆盖的区域；
对所述覆盖区域进行暴雨预警。
本发明中利用不同时间(如前后6分钟)获取的雷达图，依据其质心点和平均回波强度值进行图形匹配，获得重叠度以及图形面积满足预定条件的若干组雷达图，其中包括的第四雷达图即是第二雷达图经过6分钟的移动后所在的位置，然后依据第二雷达图和第四雷达图获取的暴雨的移动方向和暴雨移动速度，得到未来一段时间暴雨的覆盖区域，对覆盖区域进行暴雨预警，可以提高雷雨预警的实时性和准确性。
附图说明
附图是用来提供对本发明的进一步理解，并构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但不应构成对本发明的限制。在附图中，
图1为一个利用雷达数据识别跟踪暴雨的方法的流程示意图；
图2为本发明实施例中获取暴雨移动方向的示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步的详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不限定本发明的保护范围。
参考图1，一种利用雷达数据识别跟踪暴雨的方法，包括步骤：
S1、依据第一预定时间的雷达格点数据获取第一雷达图；以及从第一雷达图中获取雷达回波强度大于预设阈值的若干第二雷达图；
具体的可以是获取6分钟前的雷达格点数据。读取3km层的雷达格点数据，雷达6分钟扫描一次，扫描由雷达专用软件转为格点数据，1km*1km一个格点，每个格点数据包含经纬度、回波强度两个数据，将回波强度值用1km*1km的正方形颜色块表示按照其经纬度坐标画到电子地图对应的位置上，生成第一雷达图。
在一个实施例中，从第一雷达图中，获取雷达回波强度大于50DB的若干第二雷达图，如图2所示。对若干第二雷达图标号1、2、3……
S2、依据当前的雷达格点数据获取第三雷达图；以及从第三雷达图中获取雷达回波强度大于预设阈值的若干第四雷达图；
同理，从第三雷达图中，获取雷达回波强度大于50DB的若干第四雷达图，如图2所示。对若干第四雷达图标号A、B、C……
S3、获取各个第二雷达图的质心点以及平均回波强度值；
获得各个第二雷达图的质心点；
获取各个第二雷达图的平均回波强度值时：获取各个第二雷达图中包含的雷达格点总数，以及其中每个雷达格点的回波强度值；依据每个雷达格点的回波强度值和雷达格点总数，获得各个第二雷达图的平均回波强度值。具体的，可以依据如下方式得到：先从雷达格点数据(雷达数据每平方公里一个格点)里统计几何图形里(第二雷达图)包含的雷达格点总数n，然后读出每个雷达格点的回波强度值：db1…dbn，那么第二雷达图的回波平均值为：D1＝(db1+db2+…+dbn)/n。
S4、获取各个第四雷达图的质心点以及平均回波强度值；
获取各个第四雷达图的平均回波强度值的步骤包括：获取各个第四雷达图中包含的雷达格点总数，以及其中每个雷达格点的回波强度值；依据每个雷达格点的回波强度值和雷达格点总数，获得各个第四雷达图的平均回波强度值。具体的，可以依据如下方式得到：先从雷达格点数据(雷达数据每平方公里一个格点)里统计几何图形里(第四雷达图)包含的雷达格点总数n，然后读出每个雷达格点的回波强度值：db1…dbn，那么第四雷达图的回波平均值为：D2＝(db1+db2+…+dbn)/n。
S5、对各个第二雷达图和各个第四雷达图进行图形匹配，得到重叠度以及图形面积满足预定条件的若干组雷达图，每组雷达图包括一个第二雷达图和一个第四雷达图；
具体的，将各个第二雷达图(图2中编号为1、2、3…的雷达图)分别与各个第四雷达图(图2中编号为A、B、C…)进行匹配，相互之间重叠最大的而且面积相差最小的为匹配的雷达图，比如，如果第二雷达图1和第四雷达图A匹配，则说明第四雷达图A是从第二雷达图1移过来的。
S6、依据各组雷达图中的第二雷达图和第四雷达图的质心点和所述第一预定时间，获取暴雨的移动速度；
具体的，获取匹配的第二雷达图和第四雷达图的质心点之间的距离eg，两个质心点之间的移动时间是6分钟，距离eg除以6分钟即为该雷达图移动的速度，也是暴雨的移动速度。
S7、依据各组雷达图中的第二雷达图和第四雷达图的质心点获取暴雨的移动方向；
具体的，以正北为零度，见图2，efg组成一个直角三角形，用三角形方法计算暴雨移动的方向角度a。
S8、依据暴雨的移动速度和移动方向，获取暴雨在第二预定时间覆盖的区域；
具体的，第二预定时间可以是30分钟。
S9、对所述覆盖区域进行暴雨预警。
本发明中利用不同时间(如前后6分钟)获取的雷达图，依据其质心点和平均回波强度值进行图形匹配，获得重叠度以及图形面积满足预定条件的若干组雷达图，其中包括的第四雷达图即是第二雷达图经过6分钟的移动后所在的位置，然后依据第二雷达图和第四雷达图获取的暴雨的移动方向和暴雨移动速度，得到未来一段时间暴雨的覆盖区域，对覆盖区域进行暴雨预警，可以提高雷雨预警的实时性和准确性。
进一步地，由于雷达扫描时间为6分钟，6分钟内也有几何面积变化很大的雷达图，这些图形无法用上述方法完成匹配，则在对各个第二雷达图和各个第四雷达图进行图形匹配的步骤之后，还包括步骤：对于未完成图形匹配的第二雷达图和第四雷达图，依据第二雷达图以及第四雷达图的质心，进行两两循环计算，获得质心距离最近的若干组第二雷达图和第四雷达图，据此得到所述若干组雷达图。
进一步地，由于暴雨移动过程中可能增强也可能减弱，如果只要在暴雨路径上就发出预报预警信息，则可能存在误报。因此，在获得所述若干组雷达图之后还包括步骤：依据各组雷达图中第二雷达图的平均回波强度值、雷达格点总数，和第四雷达图的平均回波强度值、雷达格点总数，以及在所述第一预定时间监测到的风速和当前监测到的风速，获取暴雨增减系数；其中，第一预设时间可以是6分钟。具体的，可以依据如下方式获得暴雨增减系数T：第四雷达图的回波平均值为D2，与该第四雷达图匹配的第二雷达图的回波平均值为D1，则，T＝(D2*(V1/b)-D1*(V2/c))/6，其中，V1为当前风速，V2为6分钟前风速，b为第四雷达图中的雷达格点总数，c为上一张(6分钟前)第二雷达图中的雷达格点总数，T为暴雨增减系数。当T大于1时表明暴雨在增强。
在所述对覆盖区域进行暴雨预警的步骤中，对暴雨增减系数大于1的覆盖区域进行暴雨预警。而且当T大于1，向这个区域内的用户发出预警信息。表明暴雨在增强，30分钟内会漂移到用户上方。用T系数判断暴雨移动过程中随着距离和时间增强或减弱程度，据此，决定给暴雨移动路径上的企业和个人发布预报预警信息，减少误报率。
只要不违背本发明创造的思想，对本发明的各种不同实施例进行任意组合，均应当视为本发明公开的内容；在本发明的技术构思范围内，对技术方案进行多种简单的变型及不同实施例进行的不违背本发明创造的思想的任意组合，均应在本发明的保护范围之内。

Claims (5)

1.一种利用雷达数据识别跟踪暴雨的方法，其特征在于，包括步骤：
依据第一预定时间的雷达格点数据获取第一雷达图；以及从第一雷达图中获取雷达回波强度大于预设阈值的若干第二雷达图；
依据当前的雷达格点数据获取第三雷达图；以及从第三雷达图中获取雷达回波强度大于预设阈值的若干第四雷达图；
获取各个第二雷达图的质心点以及平均回波强度值；
获取各个第四雷达图的质心点以及平均回波强度值；
对各个第二雷达图和各个第四雷达图进行图形匹配，得到重叠度以及图形面积满足预定条件的若干组雷达图，每组雷达图包括一个第二雷达图和一个第四雷达图；
依据各组雷达图中的第二雷达图和第四雷达图的质心点和所述第一预定时间，获取暴雨的移动速度；
依据各组雷达图中的第二雷达图和第四雷达图的质心点获取暴雨的移动方向；
依据暴雨的移动速度和移动方向，获取暴雨在第二预定时间覆盖的区域；
对所述覆盖区域进行暴雨预警。
2.根据权利要求1所述的方法，其特征在于，获取各个第二雷达图的平均回波强度值的步骤包括：获取各个第二雷达图中包含的雷达格点总数，以及其中每个雷达格点的回波强度值；依据每个雷达格点的回波强度值和雷达格点总数，获得各个第二雷达图的平均回波强度值。
3.根据权利要求1所述的方法，其特征在于，获取各个第四雷达图的平均回波强度值的步骤包括：获取各个第四雷达图中包含的雷达格点总数，以及其中每个雷达格点的回波强度值；依据每个雷达格点的回波强度值和雷达格点总数，获得各个第四雷达图的平均回波强度值。
4.根据权利要求1所述的方法，其特征在于，在对各个第二雷达图和各个第四雷达图进行图形匹配的步骤之后，还包括步骤：
对于未完成图形匹配的第二雷达图和第四雷达图，依据第二雷达图以及第四雷达图的质心，进行两两循环计算，获得质心距离最近的若干组第二雷达图和第四雷达图，据此得到所述若干组雷达图。
5.根据权利要求1至4任一项所述的方法，其特征在于，获得所述若干组雷达图之后还包括步骤：依据各组雷达图中第二雷达图的平均回波强度值、雷达格点总数，和第四雷达图的平均回波强度值、雷达格点总数，以及在所述第一预定时间监测到的风速和当前监测到的风速，获取暴雨增减系数；
在所述对覆盖区域进行暴雨预警的步骤中，对暴雨增减系数大于1的覆盖区域进行暴雨预警。
CN201711061564.6A 2017-11-02 2017-11-02 利用雷达数据识别跟踪暴雨的方法 Active CN107656256B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711061564.6A CN107656256B (zh) 2017-11-02 2017-11-02 利用雷达数据识别跟踪暴雨的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711061564.6A CN107656256B (zh) 2017-11-02 2017-11-02 利用雷达数据识别跟踪暴雨的方法

Publications (2)

Publication Number Publication Date
CN107656256A true CN107656256A (zh) 2018-02-02
CN107656256B CN107656256B (zh) 2019-06-28

Family

ID=61096196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711061564.6A Active CN107656256B (zh) 2017-11-02 2017-11-02 利用雷达数据识别跟踪暴雨的方法

Country Status (1)

Country Link
CN (1) CN107656256B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108896995A (zh) * 2018-08-03 2018-11-27 中国航空工业集团公司雷华电子技术研究所 一种机载气象雷达雷暴识别方法
CN109359617A (zh) * 2018-10-30 2019-02-19 西北大学 一种基于格网降雨数据的识别暴雨的方法
CN110850414A (zh) * 2019-11-26 2020-02-28 北京天元创新科技有限公司 降雨预警方法及装置
CN111095033A (zh) * 2018-03-05 2020-05-01 株式会社东芝 气象雷达装置及暴雨预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237890A (zh) * 2014-09-03 2014-12-24 天津大学 一种由“列车效应”引起的暴雨识别及预报方法
CN105842756A (zh) * 2016-03-30 2016-08-10 广州市积雨云信息科技有限公司 一种雷雨信息预报方法
CN106872981A (zh) * 2017-02-17 2017-06-20 水利部南京水利水文自动化研究所 雨量雷达的降水强中心跟踪与预报方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237890A (zh) * 2014-09-03 2014-12-24 天津大学 一种由“列车效应”引起的暴雨识别及预报方法
CN105842756A (zh) * 2016-03-30 2016-08-10 广州市积雨云信息科技有限公司 一种雷雨信息预报方法
CN106872981A (zh) * 2017-02-17 2017-06-20 水利部南京水利水文自动化研究所 雨量雷达的降水强中心跟踪与预报方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
RAFFAELE BOLLA ET AL.: ""The Tracking and Prediction of High Intensity Rainstorms"", 《REMOTE SENSING REVIEWS》 *
THOMAS EINFALT ET AL.: ""A RADAR RAINFALL FORECASTING METHOD DESIGNED FOR HYDROLOGICAL PURPOSES"", 《JOURNAL OF HYDROLOGY》 *
杨洪平 等: ""新一代天气雷达灾害性天气警报和短时预报系统"", 《2006 CAMS ANNUAL REPORT》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111095033A (zh) * 2018-03-05 2020-05-01 株式会社东芝 气象雷达装置及暴雨预测方法
CN111095033B (zh) * 2018-03-05 2022-08-30 株式会社东芝 气象雷达装置及暴雨预测方法
CN108896995A (zh) * 2018-08-03 2018-11-27 中国航空工业集团公司雷华电子技术研究所 一种机载气象雷达雷暴识别方法
CN109359617A (zh) * 2018-10-30 2019-02-19 西北大学 一种基于格网降雨数据的识别暴雨的方法
CN109359617B (zh) * 2018-10-30 2021-09-14 西北大学 一种基于格网降雨数据的识别暴雨的方法
CN110850414A (zh) * 2019-11-26 2020-02-28 北京天元创新科技有限公司 降雨预警方法及装置

Also Published As

Publication number Publication date
CN107656256B (zh) 2019-06-28

Similar Documents

Publication Publication Date Title
CN107656256B (zh) 利用雷达数据识别跟踪暴雨的方法
CN108961790B (zh) 基于四维实景交通仿真的恶劣天气预警管理系统及方法
US20220036738A1 (en) Method and system for assessing and early warning ship collision risk
CN103730031B (zh) 内河桥区船载航行主动避碰系统及避碰方法
KR102006847B1 (ko) 다중규모 광학 흐름을 이용한 레이더 초단기 강수 예측을 위한 장치 및 방법
US7391358B2 (en) Weather radar echo tops forecast generation
CN108761484A (zh) 一种基于多源卫星遥感数据的海雾监测方法
CN105548984A (zh) 一种基于模糊逻辑的双偏振多普勒天气雷达地物杂波抑制方法
KR100963532B1 (ko) 기상레이더의 강수량 추정 방법
CN108873013A (zh) 一种采用多线激光雷达的道路可通行区域获取方法
Li et al. Numerical simulation study of the effect of buildings and complex terrain on the low-level winds at an airport in typhoon situation
CN102221389A (zh) 结合统计模型与动力模型的乘潮水位预报方法
EP2246715A2 (en) Method for compiling and displaying atmospheric uncertainty information
CN110203204A (zh) 一种车辆周边环境感知方法
CN104102905A (zh) 一种车道线自适应检测的方法
CN108828691A (zh) 短临降水预报方法与装置
CN104570161A (zh) 基于ec/jma全球格点预报资料的台风自动化预报方法
CN115494521B (zh) 一种基于激光雷达的机场跑道低空风切变预警方法
CN109946648B (zh) 一种车路协同下基于超宽带的车辆高精度定位方法
CN104849777A (zh) 基于交叉相关追踪法的电网设备区域对流天气预报方法
CN107193060A (zh) 一种多路径台风风暴潮快速预测方法及系统
CN116776592A (zh) 一种输电线路的风偏分析方法及系统
CN116224339A (zh) 一种s、c和x波段雷达降水组网融合方法
CN113176546B (zh) 基于nurbs建模的路侧雷达的实时滤除背景的方法
CN113176548B (zh) 基于多边形建模的路侧雷达的实时滤除背景的方法

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
CP01 Change in the name or title of a patent holder

Address after: Room 1001, room 23, No. 23, Tianan headquarters, No. 555, Panyu Avenue, east ring street, Panyu District, Panyu District, Guangdong

Patentee after: Guangzhou fengfenglei Technology Co., Ltd.

Address before: Room 1001, room 23, No. 23, Tianan headquarters, No. 555, Panyu Avenue, east ring street, Panyu District, Panyu District, Guangdong

Patentee before: Guangzhou science and Technology Co., Ltd..

CP01 Change in the name or title of a patent holder