一种基于多项式拟合的多普勒天气雷达速度退模糊方法技术

本发明专利技术公开了一种基于多项式拟合的多普勒天气雷达速度退模糊方法，具体包括如下步骤：对雷达探测到的数据进行处理，获得雷达探测到的不同仰角的极坐标格式的雷达探测速度V

【技术实现步骤摘要】
一种基于多项式拟合的多普勒天气雷达速度退模糊方法
本专利技术涉及气象雷达
，具体为一种基于多项式拟合的多普勒天气雷达速度退模糊方法。
技术介绍
多普勒天气雷达可利用多普勒效应来探测云雨的强度和移动速度等信息，有“超级千里眼”之称。多普勒天气雷达的探测资料具有很高的时空分辨率和精准度，在中尺度气象业务和研究中具有其他观测资料无可比拟的优势。然而，受到发射波波长和脉冲重复频率的限制，多普勒天气雷达能探测到的最大速度VN是有限的，当云团的真实移动速度Vr大于这个最大可探测速度后，雷达观测到的速度Vo将出现模糊，即Vo＝Vr+2nVN，其中n为-3、-2、-1、0、1、2、3这七个整数中的其中之一，且必须满足条件|Vo|≤VN，其中正速度代表远离雷达，负速度代表朝向雷达。若雷达的最大可正确探测速度VN＝32m/s，而云团的真实移动速度Vr＝35m/s，则雷达的观测速度Vo＝-29m/s，这里n必须为-1，因为其他的n值都会使得|Vo|＞VN而不满足条件。可见，雷达速度模糊会使观测速度从真实的正的极大速度突变至负值，或从真实的负的极大速度突变为正值(如附图1)。雷达探测速度的模糊扭曲了云团移动的真实信息，给多普勒雷达速度的应用带来了极大的麻烦。因此，如何对雷达探测速度中的模糊速度进行处理，使其还原成云团移动的真实移动速度(即退模糊)，已在大气探测和遥感领域被广泛关注。多普勒天气雷达速度退模糊处理是先选定好未受速度模糊影响的参考速度Vref，然后再通过不断调整n值，使得参考速度附近的速度值Vr＝Vo+2nVN与参考速度的差的绝对值小于0.8倍最大可探测速度，即|Vr-Vref|＜0.8VN。这样，将退模糊处理后的Vr作为新的参考速度Vref，继续进行点靠点退模糊处理，直到整个雷达探测平面都退模糊完成。在多普勒天气雷达速度退模糊过程中，最开始的参考速度(初始参考速度)的选择是最为关键的过程，不同的选择方法也衍生出了不同的退模糊方法。目前，现有的退模糊方法多是通过搜索弱风参考径向(即某一方位角上，不同距离的探测速度都较弱)来作为初始参考速度。然而，当遇到由地物杂波、速度模糊(如：Vr＝64m/s，VN＝32m/s，根据Vo＝Vr+2nVN，则n＝-1，Vo＝0m/s)等原因造成的假零速度线时，这些退模糊方法就可能选择出错误的参考速度，导致退模糊失败甚至错误。此外，当雷达速度场中因建筑物或地形阻挡而出现2个或以上较大方位断裂时(如附图1)，因为方位上的不连续性会导致现有方法退模糊提前结束，导致退模糊不完整。
技术实现思路
本专利技术提供了一种基于多项式拟合的多普勒天气雷达速度退模糊方法，通过初始参考速度选择方式的改进，解决了现有方法因雷达速度场存在假零速度线而造成退模糊失败的问题和因雷达速度场存在较大方位断裂而引起退模糊不完整的问题。为实现上述目的，本专利技术提供如下技术方案：一种基于多项式拟合的多普勒天气雷达速度退模糊方法，具体包括如下步骤：步骤S1：对雷达探测到的数据进行处理，获得雷达探测到的不同仰角的极坐标格式的雷达探测速度Vo；步骤S2：选取初始参考距离圈；根据不同时刻的雷达速度图，找出一个没有速度模糊且有效探测较多(雷达探测速度点密集)的雷达速度图，采用不同幂指数M的多项式对该雷达速度图中某一距离圈上的雷达探测速度Vo进行拟合，选出拟合效果较好且最小的幂指数M作为多项式的幂指数(M确定后是固定不变的，不随雷达速度数据的改变而改变)；然后利用确定了的幂指数M的多项式对所有时刻不同的距离圈的雷达探测速度Vo随方位角变化进行拟合，并将拟合值作为参考速度Vref，对雷达探测速度Vo进行退模糊处理，筛选出不同时刻对应的初始参考距离圈；步骤S3：径向退模糊处理；找出初始参考距离圈后，先对该距离圈上的雷达速度探测速度Vo进行退模糊处理，再将该距离圈上的多项式拟合值作为临近距离圈的参考速度，进行退模糊处理；依次将退模糊处理完的距离圈的多项式拟合值作为临近距离圈新的参考速度，继续执行退模糊处理，直至不同距离的距离圈全部退模糊完毕；步骤S4：全局退模糊处理；在方位距离平面内寻找没有退模糊处理的雷达探测速度Vo，对没有退模糊处理的进行退模糊处理。进一步的方案是，所述步骤S2选取初始参考距离圈具体包括如下步骤：具体地，首先根据不同时刻的雷达探测速度图，找出一个没有速度模糊且有效探测较多的雷达速度图，选取其某一距离圈(即某一距离上不同方位组成的圈)上的雷达探测速度Vo作为拟合对象，采用不同幂指数M的多项式进行拟合，其中为y为拟合得到的雷达速度，ai为多项式系数，β为方位角，找出拟合曲线与雷达探测速度Vo拟合效果较好且最小的幂指数M，作为多项式的幂指数；其次，利用确定了幂指数M的多项式对不同时刻雷达探测速度图上的不同距离圈上的雷达探测速度Vo进行拟合，获得不同距离圈上不同方位角的拟合值；然后，将某一距离圈上的拟合值作为参考速度Vref，对该距离圈上相同方位角的雷达探测速度Vo进行退模糊处理，即运用Vr＝Vo+2nVN，不断调整n使得|Vr-Vref|＜0.8VN；其中n为-3、-2、-1、0、1、2、3这七个整数中的其中之一；最后，当该距离圈退模糊处理完成后，若该距离圈上的雷达速度不存在相连两个方位的速度之差的绝对值≥0.8VN的情况，则该距离圈被视为好距离圈，在众多好距离圈中，雷达有效探测数最多的距离圈(雷达探测速度点最密集的距离圈)被作为初始参考距离圈。进一步的方案是，所述步骤S3径向退模糊处理具体包括如下步骤：步骤S31：对初始参考距离圈上的雷达探测速度Vo进行多项式拟合，并将拟合值作为该距离圈上对应方位角处雷达探测速度Vo的参考速度Vref，进行退模糊处理，将模糊速度还原，即运用Vr＝Vo+2nVN，不断调整n使得|Vr-Vref|＜0.8VN，其中n为-3、-2、-1、0、1、2、3这七个整数中的其中之一，并将退模糊处理完成的点标记为1；步骤S32：将拟合值也作为相连距离圈的参考速度Vref，沿朝向雷达和远离雷达中心两个方向继续进行退模糊处理，即运用Vr＝Vo+2nVN，不断调整n使得|Vr-Vref|＜0.8VN，依次将退模糊处理完的距离圈的多项式拟合值作为临近距离圈新的参考速度Vref，继续执行退模糊处理，直至不同距离的距离圈全部退模糊完毕，并将处理完成的点标记为1。进一步的方案是，所述步骤步骤S4全局退模糊处理具体步骤为：在方位距离平面内寻找没有退模糊处理的雷达探测速度Vo，没有经过退模糊处理的点没有被标记为1，找到未退模糊速度点后，以该点为中心搜索前后5个径向、左右2个方位范围内的所有已经退模糊处理的云团真实移动速度Vr，并将这些已经退模糊的云团真实移动Vr的平均值作为参考速度，对该点进行退模糊处理。进一步的方案是，所述步骤步骤S1中雷达探测到的数据为雷达基数据，雷达基数据的型号(CINRADSA/SB/SC/CB/CC/CD等)，对雷达基数据进行处理为对基数据解码并提取雷达探测速度Vo。进一步本文档来自技高网...

【技术保护点】
1.一种基于多项式拟合的多普勒天气雷达速度退模糊方法，其特征在于，具体包括如下步骤：/n步骤S1：对雷达探测到的数据进行处理，获得雷达探测到的不同仰角的极坐标格式的雷达探测速度V

【技术特征摘要】
1.一种基于多项式拟合的多普勒天气雷达速度退模糊方法，其特征在于，具体包括如下步骤：
步骤S1：对雷达探测到的数据进行处理，获得雷达探测到的不同仰角的极坐标格式的雷达探测速度Vo；
步骤S2：选取初始参考距离圈；
根据不同时刻的雷达探测速度图，找出一个没有速度模糊且有效探测较多的雷达速度图，采用不同幂指数M的多项式对该雷达速度图中某一距离圈上的雷达探测速度Vo进行拟合，选出拟合效果较好且最小的幂指数M作为多项式的幂指数；然后利用确定了的幂指数M的多项式对所有时刻不同的距离圈的雷达探测速度Vo随方位角变化进行拟合，并将拟合值作为参考速度Vref，对雷达探测速度Vo进行退模糊处理，筛选出不同时刻对应的初始参考距离圈；
步骤S3：径向退模糊处理；
找出初始参考距离圈后，先对该距离圈上的雷达速度探测速度Vo进行退模糊处理，再将该距离圈上的多项式拟合值作为临近距离圈的参考速度，进行退模糊处理；依次将退模糊处理完的距离圈的多项式拟合值作为临近距离圈新的参考速度，继续执行退模糊处理，直至不同距离的距离圈全部退模糊完毕；
步骤S4：全局退模糊处理；
在方位距离平面内寻找没有退模糊处理的雷达探测速度Vo，对没有退模糊处理的进行退模糊处理。


2.根据权利要求1所述的一种基于多项式拟合的多普勒天气雷达速度退模糊方法，其特征在于，所述步骤S2选取初始参考距离圈具体包括如下步骤：
具体地，首先根据不同时刻的雷达探测速度图，找出一个没有速度模糊且有效探测较多的雷达速度图，选取其某一距离圈(即某一距离上不同方位组成的圈)上的雷达探测速度Vo作为拟合对象，采用不同幂指数M的多项式进行拟合，其中为y为拟合得到的雷达速度，ai为多项式系数，β为方位角，找出拟合曲线与雷达探测速度Vo拟合效果较好且最小的幂指数M，作为多项式的幂指数；
其次，利用确定了幂指数M的多项式对不同时刻雷达探测速度图上的不同距离圈上的雷达探测速度Vo进行拟合，获得不同距离圈上不同方位角的拟合值；
然后，将某一距离圈上的拟合值作为参考速度Vref，对该距离圈上相同方位角的雷达探测速度Vo进行退模糊处理，即运用Vr＝Vo+2nVN，不断调整n使得|Vr-Vref|＜0.8VN；其中n为-3、-2、-1、0、1、2、3这七个整数中的其中之一；
最...

【专利技术属性】
技术研发人员：邹海波，吴珊珊，单九生，
申请(专利权)人：江西师范大学，
类型：发明
国别省市：江西;36