本发明涉及天气预测技术领域，具体的说是一种基于多尺度卷积神经网络的强降雨及雷暴预报方法和系统。
背景技术：
目前，利用机器学习模型进行强降雨或者雷暴预报的模型主要包括神经网络和支持向量机。这两种模型在进行预报的时候都需要用到一些气象指标，如地面温度、气压、高空温度、风场、散度、垂直速度等，以及由气象指标计算得到的一些大气不稳定因子，如k指数、ct指数、vt指数等。预测的结果在时间和空间上比较粗糙，在准确性上也精度不高。在时间上通常是预报24小时内是否有雷暴，在空间上通常预报某个城市范围内的雷暴情况。
随着人们对生活水平要求的提高，更希望能知道下一小时或者下一分钟的天气变化，不仅知道自己城市的天气变化，还想知道自己活动区域范围的天气变化。在受天气影响的某些特殊行业，如航空业，对时间和空间上的精度要求较高，通常希望能预报1小时内的雷暴天气，空间范围能缩小到1公里范围。这些靠传统的气象指标和模型已经很难做到。
技术实现要素：
针对现有技术中存在的上述不足之处，本发明要解决的技术问题是提供一种基于多尺度卷积神经网络的强降雨及雷暴预报方法和系统，以解决传统强降雨及雷暴预报系统在时间和空间上粒度较粗以及准确性不高的问题。
本发明为实现上述目的所采用的技术方案是：一种基于多尺度卷积神经网络的强降雨及雷暴预报方法，包括以下步骤：
获取原始雷达回波数据；
对所述原始雷达回波数据进行预处理，得到雷达回波时间序列图像；
取任意三帧连续发生的雷达回波时间序列图像，通过构建的雷达回波卷积神经网络推演下一帧雷达回波图；
通过光流场法和线性插值法来推断出一分钟内的雷达回波变化。
所述二维雷达回波数据在每一点的取值为该点在第三维度上的最大值。
所述对原始雷达回波数据进行预处理，得到雷达回波时间序列图像，包括以下步骤：
将三维空间中的雷达回波映射到二维地理空间，得到二维雷达回波数据；
按时间序列对所述二维雷达回波数据进行排序，得到雷达回波时间序列图像。
所述雷达回波卷积神经网络的构建，包括以下步骤：
在雷达回波时间序列图像中取连续的四帧历史雷达回波序列；
取前面三帧作为输入，第四帧作为输出；
构造卷积网络，通过前面三帧图像的变化，来推演第三帧图像到第四帧图像的变化。
还包括：将推演得到的下一帧雷达回波图加入到所述雷达回波时间序列图像中，并删除所述雷达回波时间序列图像中最早时间点的雷达回拨图。
一种基于多尺度卷积神经网络的强降雨及雷暴预报系统，包括：
获取模块，用于获取原始雷达回波数据；
预处理模块，用于对所述原始雷达回波数据进行预处理，得到雷达回波时间序列图像；
推演模块，用于取任意三帧连续发生的雷达回波时间序列图像，通过构建的雷达回波卷积神经网络推演下一帧雷达回波图；
差值计算模块，用于通过光流场法和线性插值法来推断出一分钟内的雷达回波变化。
所述预处理模块包括：
维度转换模块，用于将三维空间中的雷达回波映射到二维地理空间，得到二维雷达回波数据；
排序模块，用于按时间序列对所述二维雷达回波数据进行排序，得到雷达回波时间序列图像。
还包括神经网络构建模块，用于雷达回波卷积神经网络的构建；所述神经网络构建模块包括：
提取模块，用于在雷达回波时间序列图像中取连续的四帧历史雷达回波序列；取前面三帧作为输入，第四帧作为输出；
构造模块，用于构造卷积网络，通过前面三帧图像的变化，来推演第三帧图像到第四帧图像的变化。
还包括更新模块，用于将推演得到的下一帧雷达回波图加入到所述雷达回波时间序列图像中，并删除所述雷达回波时间序列图像中最早时间点的雷达回拨图。
本发明具有以下优点及有益效果：
1、本发明能够预报1小时内的雷暴天气，空间范围能缩小到1公里范围内。
2、本发明突破以前的经验模型，利用了大数据和神经网络，可以针对特定区域建立预测模型，使得模型更有针对性，更加准确。
附图说明
图1为本发明一个实施例的二维雷达回波图；
图2为本发明一实施例的方法流程图；
图3为本发明的卷积神经网络中任意一层的结构图；
图4为本发明的系统框架图；
图5为本发明一个实施例中的预测值和真实值的对比图，其中(a)为真实回波图，(b)为预测回波图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
图1是2016年6月23日8点06分江苏地区的一副雷达回波图，其中反射率的大小用颜色进行了标注，颜色较深的地方，如红色和紫色部分，属于反射率较强，通常该区域具有强降雨或者雷暴，一个像素在图像所占大小表示地表一公里范围的大小，雷达回波图的产生周期通常是5分钟得到一幅回波图。用雷达回波图得到的数据和传统的气象指标数据相比，在短时天气预报中，具有实时性较好、准确性高和成本较低等特点。
虽然雷达回波图是反应当前天气状况，但是我们通过多个前后相关的雷达回波图就可以构成一个空间和时间尺度上相互关联的数据，通过该数据就能预测未来某个时刻天气的状况。传统方法是通过外推法，来跟踪回波图上所表示的可能雷暴产生区域。由于外推法模型简单，在线性变化情况下预测效果较好，但是在雷暴预报等复杂情况下，精度较低。
如图2所示，本发明一个实施例的基本流程是：
对原始雷达回波图进行处理，把三维空间中的雷达回波映射到二维地理空间。二维空间中的每一点值取为在第三维度上的最大值。把二维雷达回波图按时间排序，得到雷达回波时间序列图像。
由于雷达回波图反映了某一地区和某一时刻的天气情况，我们把天气的预测转换为雷达回波图的预测问题，也就是由已有的时间序列滑动窗口，预测未来下一个时间点的序列值。这也可以看成是雷达图的重构问题，即重新构造未来某个时间点的雷达回波图。雷达回波卷积神经网络构建具体为：构造在时间上连续的四帧历史雷达回波序列；取前面三帧作为输入，第四帧作为输出；构造卷积网络，通过前面三帧图像的变化，来推演第三帧图像到第四帧图像的变化。卷积网络如图3所示，前三帧图像作为卷积网络输入，卷积网络的输出就是推演得到的第四帧图像。
迭代预测未来任意时间的雷达回波变化：取任意三帧连续发生的历史雷达回波序列；由之前得到的雷达回波卷积神经网络，推演下一帧雷达回波图。再由推演得到的雷达回波图和最近两帧的雷达回波图构成三帧连续的雷达回波序列，再重复上面过程。
这样就能利用构造的雷达回波图加入已有的时间序列预测下一个时间点的雷达回波图，如此迭代使用，可以预测未来任意时间点的雷达回波图，从而得到未来时间点的天气。
雷达回波图在地理空间中的精度通常是1公里范围，所以能得到一公里精度范围内的天气情况。时间序列的精度最多大概是5分钟，所以通过深度学习而产生的回波图也是下一个5分钟的预测结果。但我们可以通过光流场法和线性插值法(可参考论文《改进的光流法在回波外推预报中的应用》)来推断出一分钟内的雷达回波变化。最终能够得到在任意地点1公里范围内，时间精度在1分钟的强降雨或者雷暴状况。最高能得到空间精度为0.5公里，时间精度为1分钟的未来天气状况。
本发明创新点是通过深度学习的卷积神经网络，将依次变化的4个时间点雷达回波数据结合起来作为一个滑动观察窗口，来生成下一时刻的雷达回波数据，再通过该雷达回波数据和前三个时间点数据组合成新的滑动观察窗口，得到下一时刻的天气状况，并判断在这一时刻下，该点是否有强降雨或者雷暴，以此类推。
本发明是通过多尺度的2d卷积神经网络对雷达回波图进行滤波处理，并得到未来时刻的雷达回波图，并作为强降雨或者雷暴的预报依据。由于未来的雷达回波图和现有的雷达回波图具有复杂的非线性关系，所以我们基于较为复杂的resnext的卷积神经网络进行设计，其任意一层的网络模型如图3所示。该模型在输入端通过多个分离的通道对输入进行卷积处理，在输出端把处理结果进行合并，从而实现了多尺度的卷积效果。
如图4所示，本发明的系统包括：获取模块，用于获取原始雷达回波数据；预处理模块，用于对所述原始雷达回波数据进行预处理，得到雷达回波时间序列图像；推演模块，用于取任意三帧连续发生的雷达回波时间序列图像，通过构建的雷达回波卷积神经网络推演下一帧雷达回波图；差值计算模块，用于通过光流场法和线性插值法来推断出一分钟内的雷达回波变化。
所述预处理模块包括：维度转换模块，用于将三维空间中的雷达回波映射到二维地理空间，得到二维雷达回波数据；排序模块，用于按时间序列对所述二维雷达回波数据进行排序，得到雷达回波时间序列图像。
还包括神经网络构建模块，用于雷达回波卷积神经网络的构建；所述神经网络构建模块包括：提取模块，用于在雷达回波时间序列图像中取连续的四帧历史雷达回波序列；取前面三帧作为输入，第四帧作为输出；构造模块，用于构造卷积网络，通过前面三帧图像的变化，来推演第三帧图像到第四帧图像的变化。还包括更新模块，用于将推演得到的下一帧雷达回波图加入到所述雷达回波时间序列图像中，并删除所述雷达回波时间序列图像中最早时间点的雷达回拨图。
在本发明的一个具体实施例中，通过迭代重构雷达回波图获取某一地区1公里范围内，1分钟的天气变化。重构出来的雷达回波图如图5所示，对比实际测试的雷达回波图和模型重构的雷达回波图，其准确性能达到98％。传统的天气预测和雷暴预报方式通常是通过光流场外推得到。在论文《改进的光流法在回波外推预报中的应用》(张蕾，魏鸣，李南等.科学技术与工程，2014，14(32):133-137)中，其提出的改进的光流法和传统的光流法在预测的准确性上都没有超过70％。此外还有一些指标是对于某个地区是否出现雷暴来判断系统的准确性，其空间和阈值的精度要求上都要低于本发明提出的指标。