一种大气风场急流线的检测方法技术

技术编号:20117505 阅读:32 留言:0更新日期:2019-01-16 12:01
本发明专利技术公开了一种大气风场急流线的检测方法,其包括以下步骤:S1、获取大气风场格点矢量场数据,并建立格点位置二维图;S2、获取各个格点处的风向相干性系数;S3、选取急流区域候选格点;S4、选取急流区域的格点;S5、扩充二维图边界;S6、建立九宫格赋值表;S7、调整九宫格赋值;S8、根据赋值计算结果并在索引表中查询索引结果,提取骨架点;S9、补充骨架点;S10、获取每个骨架点的风力重心;S11、将每个骨架点的风力重心连成线,得到大气风场的急流线。本发明专利技术根据大气风场格点矢量数据可以自动识别急流区域并绘制急流线,加速了气象业务现代化和提高了灾害性天气预报准确率的进程。

A Method for Detecting Jet Line in Atmospheric Wind Field

The invention discloses a detection method of atmospheric wind field jet line, which includes the following steps: S1, acquiring grid point vector field data of atmospheric wind field and establishing two-dimensional map of grid point position; S2, acquiring wind direction coherence coefficient at each grid point; S3, selecting candidate grid points in jet area; S4, selecting grid points in jet area; S5, expanding two-dimensional map boundary; S6, establishing Nine-palace grid assignment. Tables; S7, adjust the nine-grid assignment; S8, according to the result of assignment calculation and query the index results in the index table, extract skeleton points; S9, supplement skeleton points; S10, get the wind gravity center of each skeleton point; S11, connect the wind gravity center of each skeleton point into a line, get the jet line of the atmospheric wind field. According to the lattice vector data of atmospheric wind field, the invention can automatically identify the jet area and draw the jet line, thus speeding up the modernization of meteorological operation and improving the accuracy of disastrous weather prediction.

【技术实现步骤摘要】
一种大气风场急流线的检测方法
本专利技术涉及风场检测
,具体涉及一种大气风场急流线的检测方法。
技术介绍
随着气象知识的逐步累积,知识表示的不断丰富,计算机知识获取能力的日益加强,特别是随着计算机智能硬件环境的日益加强,气象人工智能在加速气象业务现代化和提高灾害性天气预报准确率的进程中,必将发挥愈加明显的作用。但目前急流线的提取还处于人用手工绘制阶段,并不能画出风场又能智能分析出急流线。
技术实现思路
针对现有技术中的上述不足,本专利技术提供的一种大气风场急流线的检测方法可以根据风场数据自动画出急流线。为了达到上述专利技术目的,本专利技术采用的技术方案为:提供一种大气风场急流线的检测方法,其包括以下步骤:S1、获取大气风场格点矢量场数据,并建立格点位置二维图;S2、根据大气风场格点矢量场数据获取各个格点处的风向相干性系数;S3、将相干性系数小于相干性阈值的格点作为急流区域内的候选格点;S4、根据大气风场格点矢量场数据获取各个候选格点处的风速,将风速大于风速阈值的候选格点作为急流区域内的格点;将风速小于等于风速阈值的候选格点作为不在急流区域内的格点;S5、将格点位置二维图的边缘扩充一圈格点,将扩充格点作为不在急流区域内的格点;S6、建立九宫格赋值表,将九宫格的中心格子对应为任一风场格点,将九宫格的外围格子对应为该风场格点周围的风场格点,并对九宫格的外围格子进行赋值;S7、将不在急流区域内的格点所对应的赋值变为0,将在急流区域内的格点取其对应格子中的赋值;S8、将对应为中心格子的风场格点周围的风场格点的取值进行累加,根据累加数值在索引表中查询索引结果,并根据索引结果提取骨架点;S9、通过对风向进行分区,并根据格点风向所在分区的规则补充判断当前格点是否为骨架点;S10、获取每个骨架点的风力重心;S11、将每个骨架点的风力重心连成线,得到大气风场的急流线。进一步地,步骤S2的具体方法包括以下步骤:S2-1、获取任一点(x0,y0)处的风向与(x0+1,y0)点处的风向构成的夹角fx1;获取点(x0,y0)处的风向与(x0-1,y0)点处的风向构成的夹角fx2;获取点(x0-1,y0)处的风向与(x0+1,y0)点处的风向构成的夹角fx3;获取点(x0,y0)处的风向与(x0,y0+1)点处的风向构成的夹角fy1;获取点(x0,y0)处的风向与(x0,y0-1)点处的风向构成的夹角fy2;获取点(x0,y0-1)处的风向与(x0,y0+1)点处的风向构成的夹角fy3;S2-2、比较fx1、fx2和fx3的大小,并将最大值记为fx;比较fy1、fy2和fy3的大小,并将最大值记为fy;S2-3、根据公式获取格点(x0,y0)的风向相干性系数;其中λ1和λ2是结构矩阵的特征值。进一步地,步骤S3中的相干性阈值为π2/18;其中π为常数。进一步地,步骤S4中风速阈值为16m/s。进一步地,步骤S6中九宫格的外围格子的赋值方法为:将九宫格的外围格子从左上角开始按逆时针方向依次赋值为1、2、4、16、128、64、32、8。进一步地,步骤S8的具体方法包括以下步骤:S8-1、将累加数值在如下索引表中按从左至右的方式逐行索引;S8-2、将索引结果为0的风场格点作为骨架点,将索引结果为1的风场格点略去。进一步地,步骤S9的具体方法包括以下步骤:S9-1、将位于第一象限中且与X轴正方向夹角为67.5°的射线作为起始边、以45°为间隔范围顺时针将格点位置二维图依次分成8个区域;其中射线的起点为坐标轴原点;S9-2、对于第一个区域中的任一点(x1,y1),若满足点(x1,y1+1)、点(x1+1,y1+1)和点(x1+1,y1)不是急流区域内的格点,且点(x1-1,y1+1)和点(x1+1,y1-1)不是骨架点;或满足点(x1-1,y1)、点(x1-1,y1-1)和点(x1,y1-1)不是急流区域内的格点,且点(x1-1,y1+1)和点(x1+1,y1-1)不是骨架点,则将点(x1,y1)补充为骨架点;对于第二个区域中的任一点(x2,y2),若满足点(x2+1,y2)、点(x2+1,y2+1)和点(x2+1,y2-1)不是急流区域内的格点,且点(x2,y2+1)和点(x2,y2-1)不是骨架点;或满足点(x2-1,y2)、点(x2-1,y2-1)和点(x2-1,y2+1)不是急流区域内的格点,且点(x1-1,y1+1)和点(x1+1,y1-1)不是骨架点,则将点(x2,y2)补充为骨架点;对于第三个区域中的任一点(x3,y3),若满足点(x3+1,y3)、点(x3+1,y3-1)和点(x3,y3-1)不是急流区域内的格点,且点(x3+1,y3+1)和点(x3-1,y3-1)不是骨架点;或满足点(x3-1,y3)、点(x3-1,y3+1)和点(x3,y3+1)不是急流区域内的格点,且点(x3+1,y3+1)和点(x3-1,y3-1)不是骨架点,则将点(x3,y3)补充为骨架点;对于第四个区域中的任一点(x4,y4),若满足点(x4-1,y4+1)、点(x4,y4+1)和点(x4+1,y4+1)不是急流区域内的格点,且点(x4-1,y4)和点(x4+1,y4)不是骨架点;或满足点(x4-1,y4-1)、点(x4,y4-1)和点(x4+1,y4-1)不是急流区域内的格点,且点(x4-1,y4)和点(x4+1,y4)不是骨架点,则将点(x4,y4)补充为骨架点;对于第五个区域中的任一点(x5,y5),若满足点(x5,y5+1)、点(x5+1,y5+1)和点(x5+1,y5)不是急流区域内的格点,且点(x5+1,y5-1)和点(x5-1,y5+1)不是骨架点;或满足点(x5-1,y5)、点(x5-1,y5-1)和点(x5,y5-1)不是急流区域内的格点,且点(x5+1,y5-1)和点(x5-1,y5+1)不是骨架点,则将点(x5,y5)补充为骨架点;对于第六个区域内的任一点(x6,y6),若满足点(x6+1,y6)、点(x6+1,y6+1)和点(x6+1,y6-1)不是急流区域内的格点,且点(x6,y6+1)和点(x6,y6-1)不是骨架点;或满足点(x6-1,y6+1)、点(x6-1,y6)和点(x6-1,y6-1)不是急流区域内的格点,且点(x6,y6+1)和点(x6,y6-1)不是骨架点,则将点(x6,y6)补充为骨架点;对于第七区域内的任一点(x7,y7),若满足点(x7+1,y7)、点(x7+1,y7-1)和点(x7,y7-1)不是急流区域内的格点,且点(x7+1,y7+1)和点(x7-1,y7-1)不是骨架点;或满足点(x7-1,y7)、点(x7-1,y7+1)和点(x7,y7+1)不是急流区域内的格点,且点(x7+1,y7+1)和点(x7-1,y7-1)不是骨架点,则将点(x7,y7)补充为骨架点;对于第八区域内的任一点(x8,y8),若满足点(x8-1,y8+1)、点(x8,y8+1)和点(x8+1,y8+1)不是急流区域内的格点,且点(x8+1,y8)和点(x8-1,y8)不是骨架点;或满足点(x8-1,y8-1)、点(x8,y8-1)和点(x8+1,y8-1)不是急流区域内的格点,且点(x8+1,y8)和点(x8-1,y本文档来自技高网...

【技术保护点】
1.一种大气风场急流线的检测方法,其特征在于:包括以下步骤:S1、获取大气风场格点矢量场数据,并建立格点位置二维图;S2、根据大气风场格点矢量场数据获取各个格点处的风向相干性系数;S3、将相干性系数小于相干性阈值的格点作为急流区域内的候选格点;S4、根据大气风场格点矢量场数据获取各个候选格点处的风速,将风速大于风速阈值的候选格点作为急流区域内的格点;将风速小于等于风速阈值的候选格点作为不在急流区域内的格点;S5、将格点位置二维图的边缘扩充一圈格点,将扩充格点作为不在急流区域内的格点;S6、建立九宫格赋值表,将九宫格的中心格子对应为任一风场格点,将九宫格的外围格子对应为该风场格点周围的风场格点,并对九宫格的外围格子进行赋值;S7、将不在急流区域内的格点所对应的赋值变为0,将在急流区域内的格点取其对应格子中的赋值;S8、将对应为中心格子的风场格点周围的风场格点的取值进行累加,根据累加数值在索引表中查询索引结果,并根据索引结果提取骨架点;S9、通过对风向进行分区,并根据格点风向所在分区的规则补充判断当前格点是否为骨架点;S10、获取每个骨架点的风力重心;S11、将每个骨架点的风力重心连成线,得到大气风场的急流线。...

【技术特征摘要】
1.一种大气风场急流线的检测方法,其特征在于:包括以下步骤:S1、获取大气风场格点矢量场数据,并建立格点位置二维图;S2、根据大气风场格点矢量场数据获取各个格点处的风向相干性系数;S3、将相干性系数小于相干性阈值的格点作为急流区域内的候选格点;S4、根据大气风场格点矢量场数据获取各个候选格点处的风速,将风速大于风速阈值的候选格点作为急流区域内的格点;将风速小于等于风速阈值的候选格点作为不在急流区域内的格点;S5、将格点位置二维图的边缘扩充一圈格点,将扩充格点作为不在急流区域内的格点;S6、建立九宫格赋值表,将九宫格的中心格子对应为任一风场格点,将九宫格的外围格子对应为该风场格点周围的风场格点,并对九宫格的外围格子进行赋值;S7、将不在急流区域内的格点所对应的赋值变为0,将在急流区域内的格点取其对应格子中的赋值;S8、将对应为中心格子的风场格点周围的风场格点的取值进行累加,根据累加数值在索引表中查询索引结果,并根据索引结果提取骨架点;S9、通过对风向进行分区,并根据格点风向所在分区的规则补充判断当前格点是否为骨架点;S10、获取每个骨架点的风力重心;S11、将每个骨架点的风力重心连成线,得到大气风场的急流线。2.根据权利要求1所述的大气风场急流线的检测方法,其特征在于:所述步骤S2的具体方法包括以下步骤:S2-1、获取任一点(x0,y0)处的风向与(x0+1,y0)点处的风向构成的夹角fx1;获取点(x0,y0)处的风向与(x0-1,y0)点处的风向构成的夹角fx2;获取点(x0-1,y0)处的风向与(x0+1,y0)点处的风向构成的夹角fx3;获取点(x0,y0)处的风向与(x0,y0+1)点处的风向构成的夹角fy1;获取点(x0,y0)处的风向与(x0,y0-1)点处的风向构成的夹角fy2;获取点(x0,y0-1)处的风向与(x0,y0+1)点处的风向构成的夹角fy3;S2-2、比较fx1、fx2和fx3的大小,并将最大值记为fx;比较fy1、fy2和fy3的大小,并将最大值记为fy;S2-3、根据公式获取格点(x0,y0)的风向相干性系数;其中λ1和λ2是结构矩阵的特征值。3.根据权利要求2所述的大气风场急流线的检测方法,其特征在于:所述步骤S3中的相干性阈值为π2/18;其中π为常数。4.根据权利要求3所述的大气风场急流线的检测方法,其特征在于:所述步骤S4中风速阈值为16m/s。5.根据权利要求4所述的大气风场急流线的检测方法,其特征在于:所述步骤S6中九宫格的外围格子的赋值方法为:将九宫格的外围格子从左上角开始按逆时针方向依次赋值为1、2、4、16、128、64、32、8。6.根据权利要求5所述的大气风场急流线的检测方法,其特征在于:所述步骤S8的具体方法包括以下步骤:S8-1、将累加数值在如下索引表中按从左至右的方式逐行索引;S8-2、将索引结果为0的风场格点作为骨架点,将索引结果为1的风场格点略去。7.根据权利要求6所述的大气风场急流线的检测方法,其特征在于:所述步骤S9的具体方法包括以下步骤:S9-1、将位于第一象限中且与X轴正方向夹角为67.5°的射线作为起始边、以45°为间隔范围顺时针将格点位置二维图依次分成8个区域;其中射线的起点为坐标轴原点;S9-2、对于第一个区域中的任一点(x1,y1),若满足点(x1,y1+1)、点(x1+1,y1+1)和点(x1+1,y1)不是急流区域内的格点,且点(x1-1,y1+1)和点(x1+1,y1-1)不是骨架点;或满足点(x1-1,y1)、点(x1-1,y1-1)和点(x1,y1-1)不是急流区域内的格点,且点(x1-1,y1+1)和点(x1+1,y1-1)不是骨架点,则将点(x1,y1)补充为骨架点;对于第二个区域中的任一点(x2,y2),若满足点(x2+1,y2)、点(x2+1,y2+1)和点(x2+1,y2-1)不是急流区域内的格点...

【专利技术属性】
技术研发人员:甘建红舒红平胡文东罗飞赵卓宁刘敦龙何童丽王仲乾
申请(专利权)人:成都信息工程大学成都淞幸科技有限责任公司
类型:发明
国别省市:四川,51

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1