一种暴雨相似天气自动检索方法

本发明涉及气象监测技术领域，具体的，涉及一种暴雨相似天气自动检索方法。

背景技术：

天气预报是气象工作为国民经济和国防建设服务的重要手段，今天的天气预报主要是使用气象局收集的大量格点数据(气温、湿度、风向和风速、气压等等)，当进行某区域的气象预报时，首先根据该区域的地理位置，选择关键气象要素，气压、和/或温度、和/或湿度、和/或风场等，然后分别将各气象要素的格点数据与对应的历史数据进行比对，综合考虑各气象要素的对比结果，得到最终的预报结果。在实际的气压要素相似比对过程中，数据比对花费时间长，而且预测结果容易受到局部数据的影响，导致预测结果不准。

技术实现要素：

本发明提出一种暴雨相似天气自动检索方法，解决了现有技术中暴雨相似检索比对结果不准确的问题。

本发明的技术方案如下：包括待预报区域气压比对的步骤，具体为：

获得待预报区域气压数据、以及与所述待预报区域气压数据对应的历史数据，所述历史数据包括待预报区域内m个降雨天气的历史气压数据；

将待预报区域气压数据转换为待预报图像；

将m个降雨天气的历史气压数据分别转换为m个历史图像；

通过图像比对，计算待预报图像与m个历史图像的相似度，按照相似度从大到小的顺序，将m个历史图像进行排序，得到待预报区域气压比对结果。

进一步，将待预报区域气压数据转换为待预报图像的步骤，与将m个降雨天气的历史气压数据分别转换为m个历史图像的步骤相同，其中，将待预报区域气压数据转换为待预报图像的步骤具体为：

s210：将待预报区域气压数据转换为等值线图，并对等值线图进行填充，得到色斑图；

s220：通过特征提取，将色斑图缩小至n×n个像素点，得到n×n图片，其中n为整数，取值范围为16～64，

s230：将n×n图片转换为灰度图gp，作为待预报图像。

进一步，通过图像比对，计算待预报图像与m个历史图像的相似度，按照相似度从大到小的顺序，将m个历史图像进行排序，具体为：

计算待预报图像的哈希值hp，具体包括：

s411：将待预报图像进行dct变化得到系数矩阵，并取系数矩阵左上角的k×k矩阵，其中，k为整数，k：n的范围为1/5～3/10，

s412：计算k×k矩阵中所有元素的平均值eavg，

s413：按照元素下标从小到大的顺序，将k×k矩阵中每个元素ej与平均值eavg进行比较，如果ej≥eavg，则二进制位hpj＝1，否则hpj＝0，k²个二进制位hpj组成的字符串作为待预报图像的哈希值hp，其中，j＝1，2…k²；

采用与步骤s411～s413相同的方法，分别计算m个历史图像的哈希值hc1，hc2，…hcm；

计算待预报图像的哈希值hp与每个历史图像的哈希值之间的海明距离，得到m个海明距离ha1，ha2，…ham，按照海明距离从小到大的顺序，将m个历史图像排序。

进一步，步骤s210具体为：采用python自带填充函数contourf()对待预报区域气压数据进行处理，得到色斑图。

进一步，还包括全局气压比对的步骤，具体为：

根据地理位置，获得待预报区域的上级区域气压数据、以及与所述上级区域气压数据对应的全局历史数据，所述全局历史数据包括上级区域内m个降雨天气的历史气压数据；

采用与将待预报区域气压数据转换为待预报图像相同的步骤，将上级区域气压数据转换为上级区域图像；

采用与计算待预报图像的哈希值hp相同的步骤，计算上级区域图像的哈希值hg；

采用与将待预报区域气压数据转换为待预报图像相同的步骤，将全局历史数据转换为m个全局历史图像；

采用与计算待预报图像的哈希值hp相同的步骤，计算m个全局历史图像的哈希值hcg1，hcg2…hcgm；

计算上级区域图像的哈希值hg与每个全局历史图像的哈希值之间的全局海明距离，得到m个全局海明距离hag1，hag2，…hagm；

根据海明距离ha1，ha2，…ham和全局海明距离hag1，hag2，…hagm，得到m个关联海明距离haz1，haz2，…hazm，具体为：

haz1＝corr1×ha1+corr2×hag1，

haz2＝corr1×ha2+corr2×hag2，

：

hazm＝corr1×ham+corr2×hagm，

其中，corr1，corr2均为系数，corr1的取值范围为0.5-1，corr2的取值范围为0-0.5；

按照关联海明距离从小到大的顺序，将m个历史图像和m个全局历史图像分别排序，得到全局气压比对结果。

进一步，获得待预报区域气压数据、以及与所述待预报区域气压数据对应的历史数据，具体为：

根据经度和纬度信息，从实时监测数据表中提取待预报区域气压数据，实时监测数据表由用户输入；

根据经度和纬度信息，从历史数据库中提取历史数据，历史数据库预先保存在存储器中。

进一步，还包括：选择温度、湿度、风场中的一种或几种，分别进行相似温度数据比对、和/或相似湿度数据比对和/或相似风场数据比对，并将待预报区域气压比对结果和全局气压比对结果，与相似温度数据比对结果、和/或相似湿度数据比对结果、和/或相似风场数据比对结果相结合，得到暴雨相似天气自动检索结果。

本发明的工作原理及有益效果为：

本发明通过对历史数据进行处理，建立暴雨历史数据库；然后将暴雨历史数据库的数据进行图像化处理，得到历史图像；然后采用同样的方法将待预报区域气压数据图像化，得到待预报图像；之后将待预报图像和历史图像进行比对，将历史图像按照相似度进行排序，用于后续的暴雨预报。通过图像比对，可以大大提高相似天气的比对速度，而且图像数据在直观上保存了气象的环流特征，大大提高了相似天气比对结果的准确性。

附图说明

下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1为本发明中待预报区域气压比对方法流程图；

图2为本发明中待预报数据转换为待预报图像方法流程图；

图3为本发明中图像对比方法流程图；

图4为本发明中相似全局气压图像检索方法流程图；

图5为本发明中待预报区域原始数据示意图；

图6为待预报图像示意图；

图7为全局气压对比结果示意图；

具体实施方式

下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都涉及本发明保护的范围。

如图1所示，本实施例一种暴雨相似天气自动检索方法包括待预报区域气压比对的步骤，具体为：

s100：获得待预报区域气压数据、以及与所述待预报区域气压数据对应的历史数据，所述历史数据包括待预报区域内m个降雨天气的历史气压数据；

其中，从气象局获得的格点数据，包括经度、纬度、气压、时间等信息，根据待预报区域的经度和纬度信息，可以从格点数据中提取待预报区域的气压数据，与待预报区域气压数据对应的历史数据也可以从气象局得到，并预先保存在历史数据库中，历史数据包括待预报区域内m个降雨天气的历史气压数据。

s200：将待预报区域气压数据转换为待预报图像；

s300：将m个降雨天气的历史气压数据分别转换为m个历史图像；

s400：通过图像比对，计算待预报图像与m个历史图像的相似度，按照相似度从大到小的顺序，将m个历史图像进行排序，得到待预报区域气压比对结果。

本实施例中，通过对历史数据进行处理，建立暴雨历史数据库；然后将暴雨历史数据库的数据进行图像化处理，得到历史图像；然后采用同样的方法将待预报区域气压数据图像化，得到待预报图像；之后将待预报图像和历史图像进行比对，将历史图像按照相似度进行排序，用于后续的暴雨预报。通过图像比对，可以大大提高相似天气的比对速度，而且图像数据在直观上保存了气象的环流特征，大大提高了相似天气比对结果的准确性。

进一步，步骤s200与s300采用相同的方法，如图2所示，步骤s200具体为：

s210：将待预报区域气压数据转换为等值线图，并对等值线图进行填充，得到色斑图；

发明人经过大量试验研究得到，采用色斑图的形式与单纯的采用等值线图的形式相比，不仅能保留原始数据的数值特征，更能直观的显示天气的环流特征，因此，在得到等值线图之后，还对等值线进行了填充，得到色斑图，用于后续的图像比对。

s220：通过特征提取，将色斑图缩小至n×n个像素点，得到n×n图片，其中n为整数，取值范围为16～64，其中，n＝32时，能够取得较好的预测效果。

s230：将n×n图片转换为灰度图gp，作为待预报图像。

进一步，如图3所示，步骤s400中通过图像比对，计算待预报图像与m个历史图像的相似度，按照相似度从大到小的顺序，将m个历史图像进行排序，得到待预报区域气压比对结果，具体为：

s410：计算待预报图像的哈希值hp，具体包括：

s411：将待预报图像进行dct变化得到系数矩阵，并取系数矩阵左上角的k×k矩阵，其中，k为整数，k：n的范围为1/5～3/10，最优值为k＝8，n＝32。

利用dct变化(离散余弦变化)将待预报图像从像素域变换为频率域，将待预报图像表示为由不同频率分量叠加而成的形式。待预报图像转换到频率域后，大部分频率分量对应的系数很小，只有少数频率分量对应的系数较大，且系数较大的值一般集中在系数矩阵的左上角，沿对角线方向越往右下角则系数矩阵元素值越小，左上角部分保留了图片的绝大多数信息。因此，取系数矩阵左上角的k×k矩阵来进行图像哈希值的计算，能增加更多的图像细节。dct变换如下面公式所示：

其中x,y是像素域中元素的坐标，f(x,y)是对应元素的值,是像素的值，n是像素矩阵的阶。u,v是频率域中元素的坐标，f(x,y)是转换后频域的系数矩阵的元素，将系数矩阵记为鼠mn×n。

s412：计算k×k矩阵中所有元素的平均值eavg，

s413：按照元素下标从小到大的顺序，将k×k矩阵中每个元素ej与平均值eavg进行比较，如果ej≥eavg，则二进制位hpj＝1，否则hpj＝0，k²个二进制位hpj组成的字符串作为待预报图像的哈希值hp，其中，j＝1，2…k²；

s420：采用与步骤s411～s413相同的方法，分别计算m个历史图像的哈希值hc1，hc，…hcm；

s430：计算待预报图像的哈希值hp与每个历史图像的哈希值之间的海明距离，得到m个海明距离ha1，ha2，…ham，按照海明距离从小到大的顺序，将m个历史图像排序。

两张图片哈希值的海明距离越小图片越相似，距离越大图片差异性越大。

进一步，步骤s210具体为：采用python自带填充函数contourf()对待预报区域气压数据进行处理，得到色斑图。

采用python自带填充函数contourf()对待预报区域气压数据进行处理，生成的等值线填充图(即色斑图)保留了原始数据的数值特征，以可视化的方式表达了原始数据中的天气环流系统。

进一步，发明人在试验过程中发现，待预报区域的上一级区域的气压数据(即上级区域气压数据)也会对待预报区域的降雨情况产生一定的影响，因此，本实施例将上级区域气压数据考虑在内，进行全局气压比对，如图4所示，具体为：

s510：根据地理位置，获得待预报区域的上级区域气压数据、以及与所述上级区域气压数据对应的全局历史数据，所述全局历史数据包括上级区域内m个降雨天气的历史气压数据；

根据上级区域的经度和纬度信息，可以从格点数据中提取上级区域的气压数据，与上级区域气压数据对应的全局历史数据也可以从气象局得到，并预先保存在历史数据库中，全局历史数据包括上级区域内m个降雨天气的历史气压数据。

s520：采用与将待预报区域气压数据转换为待预报图像相同的步骤(步骤s210～s230)，将上级区域气压数据转换为上级区域图像；

s530：采用与计算待预报图像的哈希值hp相同的步骤(步骤s411～s413)，计算上级区域图像的哈希值hg；

s540：采用与将待预报区域气压数据转换为待预报图像相同的步骤(步骤s210～s230)，将m组全局数据分别转换为m个全局历史图像；

s550：采用与计算待预报图像的哈希值hp相同的步骤(步骤s411～s413)，计算m个全局历史图像的哈希值hcg1，hcg2…hcgm；

s560：计算上级区域图像的哈希值hg与每个全局历史图像的哈希值之间的全局海明距离，得到m个全局海明距离hag1，hag2，…hagm；

s570：根据海明距离ha1，ha2，…ham和全局海明距离hag1，hag2，…hagm，得到m个关联海明距离haz1，haz2，…hazm，具体为：

haz1＝corr1×ha1+corr2×hag1，

haz2＝corr1×ha2+corr2×hag2，

hazm＝corr1×ham+corr2×hagm，

其中，corr1，corr2均为系数，corr1的取值范围为0.5-1，corr2的取值范围为0-0.5；

s570：按照关联海明距离从小到大的顺序，将m个历史图像和m个全局历史图像分别排序，得到全局气压比对结果。

本实施例中，海明距离ha1表明待预报图片与第一个历史图片的相似度，以此类推，海明距离ham表明待预报图片与第m个历史图片的相似度；同理，全局海明距离hag1表明上级区域图片与第一个全局历史图片的相似度，以此类推，全局海明距离hagm表明上级区域图片与第m个全局历史图片的相似度。

通过为ha1、hag1分配不同的权重，进行关联海明距离haz1的计算，以此类推，通过为ham、hagm分配不同的权重，进行关联海明距离hazm的计算，得到m个关联海明距离haz1，haz2，…hazm；然后通过关联海明距离haz1，haz2，…hazm的比较，分别对历史图像、全局历史图像进行排序，这样将待预报区域气压数据和上级区域气压数据均考虑进气压图像比对中，有利于得到更加准确的预测结果。经发明人的大量实验研究表明，当海明距离ha1的权重系数corr1为0.8，全局海明距离ha2的权重系数corr2为0.2时，有利于得到准确的预报结果。

进一步，本实施方法还包括：

选择温度、湿度、风场中的一种或几种，分别进行相似温度数据比对、和/或相似湿度数据比对和/或相似风场数据比对，并将待预报区域气压比对结果和全局气压比对结果，与相似温度数据比对结果、和/或相似湿度数据比对结果、和/或相似风场数据比对结果相结合，得到暴雨相似天气自动检索结果。

除了气压因素外，还需要考虑与降雨有关的其他气象要素，例如温度、湿度、风场等，根据专家意见以及各种预报资料，选择温度、湿度、风场等气象要素，分别进行比对，综合考虑各气象要素的比对结果，得到最终的预报结果。

本实施例中，发明人主要研究的是石家庄地区的气象预报，以北纬10-80，东经30-160作为上级区域，北纬30-60，东经90-130作为待预报区域，待预报区域的原始气压数据如图4所示。根据待预报区域的原始气压数据得到色斑图，如图5所示。通过全局气压比对、计算关联海明距离之后，按照关联海明距离从大到小的顺序，将各历史图像排序，如图7a-7h所示。其中，以图7a为例，左边两幅图像分别为上级区域图像和全局历史图像，图像文件名后两位均为aa；右边两幅图像分别为待预报图像和历史图像，图像文件名后两位均为la。分别计算上级区域图像和各全局历史图像的全局海明距离、待预报图像和各历史图像的海明距离之后，计算关联海明距离，即为各图像的相似度。

各图像相似度计算结果具体为：图7a相似度0.809，图7b相似度0.750，图7c相似度0.723，图7d相似度0.714，图7e相似度0.703，图7f相似度0.699，图7g相似度0.691，图7h相似度0.681。

以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。