CN103149552A - 一种多普勒天气雷达径向速度场的退模糊方法 - Google Patents
一种多普勒天气雷达径向速度场的退模糊方法 Download PDFInfo
- Publication number
- CN103149552A CN103149552A CN2013100418408A CN201310041840A CN103149552A CN 103149552 A CN103149552 A CN 103149552A CN 2013100418408 A CN2013100418408 A CN 2013100418408A CN 201310041840 A CN201310041840 A CN 201310041840A CN 103149552 A CN103149552 A CN 103149552A
- Authority
- CN
- China
- Prior art keywords
- noise
- velocity field
- radial velocity
- weather radar
- fuzzy
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多普勒天气雷达径向速度场的退模糊方法,属于气象雷达数据质量控制领域。该方法包括如下步骤:第一步,用噪声分离方法将噪声数据从径向速度场中分离出来;第二步,用多曲线拟合的方法退速度模糊;第三步,将第一步中误删除的非噪声数据恢复到原位置。本发明能直接改善天气雷达速度场的质量,间接提升短时临近预报、气象灾害预警、人工影响天气等业务的水平,有较高的应用价值和较好的应用前景。
Description
技术领域
本发明涉及一种多普勒天气雷达径向速度场的退模糊方法,属于气象雷达数据质量控制领域。
背景技术
速度场是多普勒天气雷达获取的观测数据之一,它广泛应用于短时临近预报、气象灾害预警、人工影响天气、数值预报模式的资料同化等领域。雷达通过测量相邻脉冲间的相位变化获取径向速度,测速范围为( ,),是最大不模糊速度。当真实的速度值超出测速范围时,会折叠到(,)内,产生错误的速度值,这就是速度模糊。速度模糊是天气雷达广泛存在的一个问题,严重限制了速度场的应用。
把测量值恢复成真实值,这被称作速度退模糊。速度退模糊可以通过硬件和软件两种方法来实现。硬件的方法通过双脉冲重复频率或多脉冲重复频率来实现。例如:何平等(2012)提出用多脉冲重复频率TPRF法,扩展的方法。但是,硬件方法存在非均一采样所引起的数据质量下降和现有天气雷达网硬件升级的问题。另外,即使扩展了,当遇到大风时,如台风、龙卷,仍然会出现模糊。所以,软件方法成了低成本解决这一问题的途径。
近三十年来,国内外学者提出了多个软件速度退模糊方法。有沿径向一维的《依据径向数据的连续性判断模糊》、沿切向一维的《依据切向数据的连续性判断模糊》、沿径向和切向二维的、沿径向、切向、垂直、时间四维的方法。这些方法都是在假设速度场是连续的,通过判断速度值的突然变化,来消除模糊。当速度场不连续时,方法会产生错误。退模糊错误分为两种:第一种是模糊数据没有被校正,第二种是不模糊数据变成模糊数据(或模糊数据变成更模糊数据)。后一种错误是退模糊方法的负作用,称之为“污染”。“污染”使退模糊后的速度场更难以理解和应用,也是方法在应用中遇到的最大问题。产生“污染”的主要原因是噪声的干扰。
从数据分布来看,噪声分为两种:孤立噪声和连续噪声。当一个噪声点的周围,好点(非噪声点包括缺测点)数大于噪声点数时,该点为孤立噪声;好点数小于等于噪声点数时,该点为连续噪声。两种噪声对退模糊算法的影响和消除的难度是不同的。孤立噪声可以参照周围点的分布特征来消除,所以它对退模糊算法的影响较小。然而,连续噪声的消除是非常困难的,它对算法的影响是严重的。连续噪声通常出现在地物区、低信噪比区和高谱宽区。由于连续噪声会使某一区域内的噪声点数大于好点数,因此,退模糊方法中常用的中值滤波、局部平滑的消除方法对连续噪声不仅无效而且会使其变得更连续。另外,退模糊方法中噪声的消除方法,在消除连续噪声的同时也会误删大量的非噪声点。所以,这类方法在退模糊算法中的应用时是比较谨慎的,它们对连续噪声的抑制能力也比较有限。
因此,如何在退模糊过程中有效地抑制噪声,同时又不损失风场信息是解决天气雷达速度模糊问题的关键。
发明内容:
本发明提出了一种多普勒天气雷达径向速度场的退模糊方法,采用“分离—恢复”无损的噪声抑制技术,可以在不损失任何有效数据的前提下,正确校正速度场中的模糊。
本发明为解决其技术问题采用如下技术方案:
一种多普勒天气雷达径向速度场的退模糊方法,包括如下步骤:第一步,用噪声分离方法将噪声数据从径向速度场中分离出来;第二步,用多曲线拟合的方法退速度模糊;第三步,将第一步中误删除的非噪声数据恢复到原位置。
所述噪声分离方法包括处理地物噪声、低信噪比噪声和高谱宽噪声方法。
本发明的有益效果如下:
本发明可在无需升级现有雷达硬件的前提下,解决速度模糊的问题。经过3年4站的连续观测数据验证(包括台风、晴空、层云、弱对流、强对流、龙卷等类型模糊),退模糊正确率接近90%(以体扫文件为单位),对台风、强对流的模糊正确率>94%。由于退模糊过程中融合了连续噪声抑制技术,方法性能明显优于传统方法,正确率提升约30%。因此,本发明能直接改善天气雷达速度场的质量,间接提升短时临近预报、气象灾害预警、人工影响天气等业务的水平,有较高的应用价值和较好的应用前景。
附图说明
图1是本发明的方法步骤流程图。
图2是曲线拟合的所有参与点的示意图。
具体实施方式
下面结合附图对本发明创造做进一步详细说明。
本发明是一个三维的速度退模糊方法。采用“分离—恢复”噪声的无损抑制方案,由噪声分离、曲线退模糊、噪声恢复三个步骤组成。图1是方法流程图,其中VAD是Velocity Azimuth Display的缩写,是一种沿切向拟合的简谐曲线;MVAD是 Modified VAD 的缩写,是一种改进的VAD拟合曲线。方法的第一步是噪声分离。设计了三个噪声分离方法,使用严格的阈值把噪声点分离出来。第二步是曲线退模糊。使用三条拟合的曲线的方法,在退模糊过程中进一步抑制残留的噪声,同时用多曲线协调技术适应不同尺度风场。第三步是噪声恢复。第一步分离出来的噪声中包含非噪声点,这一步将噪声逐点恢复到原位置,然后校正误删除的非噪声点的模糊。
第一步、噪声分离
噪声的产生原因是多方面的,可由地物、低信噪比、气象目标的高脉动(高谱宽)、生物、电磁干扰、距离折叠、超折射、测量误差等引起。有些噪声并非真正意义的噪声,如地物区噪声,而是真实的测量值,但由于其与云雨的径向速度差异较大,因此认为其是噪声。噪声分离由三个分离方法实现,分别用于分离地物区、低信噪比区和高谱宽区的连续噪声。
1) 地物区噪声分离
地物区噪声较多,尤其是在低仰角的速度场中。定义三个条件:一是回波点所在高度小于阈值;二是回波点的强度大于阈值;三是回波点的径向速度绝对值小于阈值。将满足这三个条件的所有回波点分离出来,相应位置用缺测值代替。
2) 低信噪比区噪声分离
由于回波功率低,低信噪比区的速度值常常是不可靠的,如远距离回波边缘的速度点。分离低信噪比区的噪声,首先将反射率转换成SNR(信号噪声比),然后将SNR小于信号噪声比阈值的回波点分离出来,相应位置用缺测值代替。
3) 高谱宽区噪声分离
第二步、曲线退模糊
1) 基本流程
退模糊顺序在垂直向上从最高仰角到最低仰角,在方位向上是从初始方位起顺时针执行,在径向上是从雷达中心到最远处。退模糊一个点时,首先计算当前位置的参考值(详细描述在后面),然后依据当前点与参考值的差值,判断当前点是否模糊。如果当前点是模糊的,则恢复它的真实速度值。最后,对当前点进行错误检查。如果当前点未通过错误检查,则认为它是残留噪声,将它分离出来,当前位置用缺测值代替。
2) 参考值的计算
用已正确校正的不模糊的数据拟合三条曲线:沿径向的中直线(中尺度:20~200km)、中直线(中尺度:2~20km)和VAD(Velocity Azimuth Display,是一种沿切向拟合的简谐曲线)曲线。三条曲线拟合时使用的数据点见图2,“EL0”表示当前层;“EL1”表示上一层;箭头指向的“A”点是表示当前点,即当前将要处理的点,“P”点表示当前径向(“A”点所在的径向)上已退模糊的点;“U”点表示在上一层中与“A”点同方位角的径向上的点;“R”点表示在上一层中与“A”点同一高度的点。用“P” 点和“U”点在不同距离(尺度)上拟合中直线和中直线;用“R”点拟合曲线。然后,用三条曲线分别计算当前点位置的估计值,、。以三个曲线的标准方差,,为权重通过公式3,得到当前点的退模糊参考值。
3) 初始层的处理
退模糊算法初始时,由于没有有效的参考,所以算法容易出现错误。最高仰角层的是方法的初始层,由于没有上层数据作参考,所以处理方法与其它层不同。退模糊时,沿方位向拟合一条VAD曲线和一条MVAD(Modified VAD)曲线,选择标准方差小的一条作为参考曲线。利用这条参考曲线计算当前点位置的参考值,依据恢复当前点的真实速度值并进行更严格的错误检查。经过噪声分离后,最高仰角层的噪声已经大量减少,且最高仰角层的方位向数据率覆盖率比较高。噪声少,方位向数据覆盖率高,这两个条件可以满足VAD和MVAD的要求。
4) 错误检查
为了防止退模糊的错误扩散,AND算法在退模糊一个点后,立即进行错误检查。方法是计算退模糊后的速度值与参考值的差值,如果差值的绝对值大于阈值,则认为当前点是残留噪声,将其分离出来,当前位置用缺测值填补。
第三步、噪声恢复
Claims (2)
1.一种多普勒天气雷达径向速度场的退模糊方法,其特征在于,包括如下步骤:第一步,用噪声分离方法将噪声数据从径向速度场中分离出来;第二步,用多曲线拟合的方法退速度模糊;第三步,将第一步中误删除的非噪声数据恢复到原位置。
2.根据权利要求1所述一种多普勒天气雷达径向速度场的退模糊方法,其特征在于,所述噪声分离方法包括处理地物噪声、低信噪比噪声和高谱宽噪声方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310041840.8A CN103149552B (zh) | 2013-02-04 | 2013-02-04 | 一种多普勒天气雷达径向速度场的退模糊方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310041840.8A CN103149552B (zh) | 2013-02-04 | 2013-02-04 | 一种多普勒天气雷达径向速度场的退模糊方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103149552A true CN103149552A (zh) | 2013-06-12 |
CN103149552B CN103149552B (zh) | 2015-11-18 |
Family
ID=48547738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310041840.8A Active CN103149552B (zh) | 2013-02-04 | 2013-02-04 | 一种多普勒天气雷达径向速度场的退模糊方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103149552B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105988114A (zh) * | 2015-02-02 | 2016-10-05 | 福建省气象台 | 多普勒雷达风速资料的填补方法 |
CN107526067A (zh) * | 2017-03-07 | 2017-12-29 | 中国气象局武汉暴雨研究所 | 全自动多普勒天气雷达径向速度退模糊算法 |
CN108459323A (zh) * | 2018-01-31 | 2018-08-28 | 天津大学 | 基于多普勒天气雷达数据的强辐合场识别及速度订正方法 |
CN108535730A (zh) * | 2018-04-16 | 2018-09-14 | 青海大学 | 一种多普勒气象雷达解速度模糊方法和系统 |
CN108562903A (zh) * | 2017-12-25 | 2018-09-21 | 天津大学 | 基于多普勒天气雷达的强对流系统动力场结构识别方法 |
CN109283535A (zh) * | 2018-11-16 | 2019-01-29 | 北京无线电测量研究所 | 非均匀风场的探测方法 |
CN111983570A (zh) * | 2020-07-16 | 2020-11-24 | 北京宏锐星通科技有限公司 | 对sar-gmti模式的干扰方法 |
CN113835091A (zh) * | 2021-07-05 | 2021-12-24 | 浙江省气象信息网络中心 | 适应多种天气事件的多普勒雷达径向速度自动退模糊方法 |
CN114488033A (zh) * | 2022-04-14 | 2022-05-13 | 南京信息工程大学 | 一种瞬时宽带多频复合雷达信号多通道时频积累处理方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017203A (zh) * | 2006-09-26 | 2007-08-15 | 南京大桥机器有限公司 | 便携式x波段多普勒天气雷达信号处理方法及装置 |
CN101706571A (zh) * | 2009-11-05 | 2010-05-12 | 宁波市气象信息中心 | 一种雷达回波中非气象杂波的处理方法 |
EP2290396A2 (en) * | 2009-08-31 | 2011-03-02 | Kabushiki Kaisha Toshiba | Doppler radar apparatus and method of calculating Doppler velocity |
CN102385054A (zh) * | 2010-09-02 | 2012-03-21 | 何平 | 多普勒天气雷达退速度模糊的tprf方法 |
CN102508219A (zh) * | 2011-10-17 | 2012-06-20 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测方法 |
-
2013
- 2013-02-04 CN CN201310041840.8A patent/CN103149552B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017203A (zh) * | 2006-09-26 | 2007-08-15 | 南京大桥机器有限公司 | 便携式x波段多普勒天气雷达信号处理方法及装置 |
EP2290396A2 (en) * | 2009-08-31 | 2011-03-02 | Kabushiki Kaisha Toshiba | Doppler radar apparatus and method of calculating Doppler velocity |
CN101706571A (zh) * | 2009-11-05 | 2010-05-12 | 宁波市气象信息中心 | 一种雷达回波中非气象杂波的处理方法 |
CN102385054A (zh) * | 2010-09-02 | 2012-03-21 | 何平 | 多普勒天气雷达退速度模糊的tprf方法 |
CN102508219A (zh) * | 2011-10-17 | 2012-06-20 | 中国人民解放军理工大学气象学院 | 风廓线雷达湍流目标检测方法 |
Non-Patent Citations (1)
Title |
---|
刘淑媛等: "一种简易的多普勒雷达速度模糊纠正技术", 《应用气象学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105988114A (zh) * | 2015-02-02 | 2016-10-05 | 福建省气象台 | 多普勒雷达风速资料的填补方法 |
CN107526067A (zh) * | 2017-03-07 | 2017-12-29 | 中国气象局武汉暴雨研究所 | 全自动多普勒天气雷达径向速度退模糊算法 |
CN108562903A (zh) * | 2017-12-25 | 2018-09-21 | 天津大学 | 基于多普勒天气雷达的强对流系统动力场结构识别方法 |
CN108562903B (zh) * | 2017-12-25 | 2021-10-01 | 天津大学 | 基于多普勒天气雷达的强对流系统动力场结构识别方法 |
CN108459323A (zh) * | 2018-01-31 | 2018-08-28 | 天津大学 | 基于多普勒天气雷达数据的强辐合场识别及速度订正方法 |
CN108459323B (zh) * | 2018-01-31 | 2022-04-15 | 天津大学 | 基于多普勒天气雷达数据的强辐合场识别及速度订正方法 |
CN108535730B (zh) * | 2018-04-16 | 2020-09-29 | 青海大学 | 一种多普勒气象雷达解速度模糊方法和系统 |
CN108535730A (zh) * | 2018-04-16 | 2018-09-14 | 青海大学 | 一种多普勒气象雷达解速度模糊方法和系统 |
CN109283535A (zh) * | 2018-11-16 | 2019-01-29 | 北京无线电测量研究所 | 非均匀风场的探测方法 |
CN109283535B (zh) * | 2018-11-16 | 2022-12-20 | 北京无线电测量研究所 | 非均匀风场的探测方法 |
CN111983570A (zh) * | 2020-07-16 | 2020-11-24 | 北京宏锐星通科技有限公司 | 对sar-gmti模式的干扰方法 |
CN113835091A (zh) * | 2021-07-05 | 2021-12-24 | 浙江省气象信息网络中心 | 适应多种天气事件的多普勒雷达径向速度自动退模糊方法 |
CN113835091B (zh) * | 2021-07-05 | 2024-05-28 | 浙江省气象信息网络中心 | 适应多种天气事件的多普勒雷达径向速度自动退模糊方法 |
CN114488033A (zh) * | 2022-04-14 | 2022-05-13 | 南京信息工程大学 | 一种瞬时宽带多频复合雷达信号多通道时频积累处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103149552B (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103149552B (zh) | 一种多普勒天气雷达径向速度场的退模糊方法 | |
CN110458055B (zh) | 一种障碍物检测方法及系统 | |
CN104036146B (zh) | 一种用于雷达目标点迹凝聚的点迹聚类方法 | |
CN106093890A (zh) | 建立剩余杂波图抑制雷达剩余杂波的方法 | |
CN103439697B (zh) | 基于动态规划的目标检测方法 | |
CN103777178B (zh) | 一种同步误差补偿方法、设备及系统 | |
CN107884822A (zh) | 一种提高矿用微震震源定位精度的方法 | |
CN106054169A (zh) | 基于跟踪信息的多站雷达信号融合检测方法 | |
CN103487788B (zh) | 一种序列脉冲信号快速自动提取方法 | |
CN105022057A (zh) | 基于改进Radon变换与多帧联合处理的目标检测方法 | |
CN105223614A (zh) | 一种基于dwt_sta/lta的含噪信号p波初至峰度拾取方法 | |
CN101813725A (zh) | 一种低频率信号的相位差测量方法 | |
CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及系统 | |
CN105676205A (zh) | 一种机载LiDAR波形数据高斯分解方法 | |
CN108957403B (zh) | 一种基于广义互相关的高斯拟合包络时延估计方法及系统 | |
CN103197297B (zh) | 基于认知框架的雷达动目标检测方法 | |
CN104570108A (zh) | 估算等效品质因子方法及用其估算地层品质因子的方法 | |
CN104765064A (zh) | 一种微地震干涉成像的方法 | |
CN104197933A (zh) | 一种望远镜视场中高星等滑动恒星的增强和提取方法 | |
CN107219554A (zh) | 陆地地震资料的剩余静校正量的自动获取方法 | |
CN104881561A (zh) | 一种基于Hough变换的多维参数的检测前跟踪方法 | |
CN104656131A (zh) | 一种基于波形恢复技术的提高地震层位自动追踪精度的方法 | |
CN104503432A (zh) | 一种基于小波能量的自主式水下机器人故障辨识方法 | |
CN105116387A (zh) | 基于位置和多普勒速度信息的pd雷达抗速度拖引方法 | |
CN106990403A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200720 Address after: Room 701-40, building 2, Zhicheng Park, no.6, Zhida Road, Jiangbei new district, Nanjing, Jiangsu Province, 210032 Patentee after: NANJING XINDA METEOROLOGICAL SCIENCE AND TECHNOLOGY Co.,Ltd. Address before: 210044 Nanjing Ning Road, Jiangsu, No. six, No. 219 Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY |