CN109213759A - 一种大气风场急流线的检测方法 - Google Patents
一种大气风场急流线的检测方法 Download PDFInfo
- Publication number
- CN109213759A CN109213759A CN201810868066.0A CN201810868066A CN109213759A CN 109213759 A CN109213759 A CN 109213759A CN 201810868066 A CN201810868066 A CN 201810868066A CN 109213759 A CN109213759 A CN 109213759A
- Authority
- CN
- China
- Prior art keywords
- point
- lattice
- wind
- skeletal
- torrent
- 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
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种大气风场急流线的检测方法,其包括以下步骤:S1、获取大气风场格点矢量场数据,并建立格点位置二维图;S2、获取各个格点处的风向相干性系数;S3、选取急流区域候选格点;S4、选取急流区域的格点;S5、扩充二维图边界;S6、建立九宫格赋值表;S7、调整九宫格赋值;S8、根据赋值计算结果并在索引表中查询索引结果,提取骨架点;S9、补充骨架点;S10、获取每个骨架点的风力重心;S11、将每个骨架点的风力重心连成线,得到大气风场的急流线。本发明根据大气风场格点矢量数据可以自动识别急流区域并绘制急流线,加速了气象业务现代化和提高了灾害性天气预报准确率的进程。
Description
技术领域
本发明涉及风场检测技术领域,具体涉及一种大气风场急流线的检测方法。
背景技术
随着气象知识的逐步累积,知识表示的不断丰富,计算机知识获取能力的日益加强,特别是随着计算机智能硬件环境的日益加强,气象人工智能在加速气象业务现代化和提高灾害性天气预报准确率的进程中,必将发挥愈加明显的作用。但目前急流线的提取还处于人用手工绘制阶段,并不能画出风场又能智能分析出急流线。
发明内容
针对现有技术中的上述不足,本发明提供的一种大气风场急流线的检测方法可以根据风场数据自动画出急流线。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种大气风场急流线的检测方法,其包括以下步骤:
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,y8)不是骨架点,则将点(x8,y8)补充为骨架点。
进一步地,步骤S10的具体方法包括以下步骤:
S10-1、根据任一骨架点左、右、上、下、左上、左下、右上和右下共八个风场格点的大气风场格点矢量场数据分别得到所述八个风场格点的风力;
S10-2、将所述八个风场格点中不在急流区域内的格点的风力取值变为0,并保持在急流区域内的格点的风力值作为其风力取值;
S10-3、根据公式
和
得到任一骨架点i的风力重心坐标(Xi,Yi),即骨架点的风力重心;其中j表示左、右、上、下、左上、左下、右上和右下八个方位中的任一个;xi,j表示位于骨架点i的j方位的风场格点的横坐标;yi,j表示位于骨架点i的j方位的风场格点的纵坐标;Fi,j表示位于骨架点i的j方位的风场格点的风力取值。
本发明的有益效果为:本发明根据大气风场格点矢量数据可以自动识别急流区域并绘制急流线,相比人工处理效率更高、准确度更高,加速了气象业务现代化和提高了灾害性天气预报准确率的进程。
附图说明
图1为本发明的流程示意图;
图2为九宫格赋值示意图;
图3为格点位置二维图分区示意图;
图4为风场格点构成急流区域的示意图;
图5为骨架点及骨架点的风力重心连成急流线的示意图;
图6为图5中框中部分的局部放大示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1和图4所示,该大气风场急流线的检测方法包括以下步骤:
S1、通过读取micaps第11类数据等方法获取大气风场格点矢量场数据,并建立格点位置二维图;
S2、根据大气风场格点矢量场数据获取各个格点处的风向相干性系数;
S3、将相干性系数小于相干性阈值的格点作为急流区域内的候选格点;
S4、根据大气风场格点矢量场数据获取各个候选格点处的风速,将风速大于风速阈值的候选格点作为急流区域内的格点;将风速小于等于风速阈值的候选格点作为不在急流区域内的格点;其中风速阈值为16m/s;
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;其中π为常数。
如图2所示,步骤S6中九宫格的外围格子的赋值方法为:将九宫格的外围格子从左上角开始按逆时针方向依次赋值为1、2、4、16、128、64、32、8。
步骤S8的具体方法包括以下步骤:
S8-1、将累加数值在如下索引表中按从左至右的方式逐行索引;
S8-2、将索引结果为0的风场格点作为骨架点,将索引结果为1的风场格点略去。
如图3所示,步骤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,y8)不是骨架点,则将点(x8,y8)补充为骨架点。
步骤S10的具体方法包括以下步骤:
S10-1、根据任一骨架点左、右、前、后、左上、左下、右上和右下共八个风场格点的大气风场格点矢量场数据分别得到所述八个风场格点的风力;
S10-2、将所述八个风场格点中不在急流区域内的格点的风力取值变为0,并保持在急流区域内的格点的风力值作为其风力取值;
S10-3、根据公式
和
得到任一骨架点i的风力重心坐标(Xi,Yi),即骨架点的风力重心;其中j表示左、右、前、后、左上、左下、右上和右下八个方位中的任一个;xi,j表示位于骨架点i的j方位的风场格点的横坐标;yi,j表示位于骨架点i的j方位的风场格点的纵坐标;Fi,j表示位于骨架点i的j方位的风场格点的风力取值。
在本发明的一个实施例中,如图4所示,图中圆形点为急流区域内的格点,圆形格点连成的线即为急流区域;如图5和图6所示,图中三角形的点为骨架点,五角星形为骨架点的风力重心,粗线即为急流线。
综上所述,本发明根据大气风场格点矢量数据可以自动识别急流区域并绘制急流线,相比人工处理效率更高、准确度更高,加速了气象业务现代化和提高了灾害性天气预报准确率的进程。
Claims (8)
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)不是急流区域内的格点,且点(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,y8)不是骨架点,则将点(x8,y8)补充为骨架点。
8.根据权利要求7所述的大气风场急流线的检测方法,其特征在于:所述步骤S10的具体方法包括以下步骤:
S10-1、根据任一骨架点左、右、上、下、左上、左下、右上和右下共八个风场格点的大气风场格点矢量场数据分别得到所述八个风场格点的风力;
S10-2、将所述八个风场格点中不在急流区域内的格点的风力取值变为0,并保持在急流区域内的格点的风力值作为其风力取值;
S10-3、根据公式
和
得到任一骨架点i的风力重心坐标(Xi,Yi),即骨架点的风力重心;其中j表示左、右、上、下、左上、左下、右上和右下八个方位中的任一个;xi,j表示位于骨架点i的j方位的风场格点的横坐标;yi,j表示位于骨架点i的j方位的风场格点的纵坐标;Fi,j表示位于骨架点i的j方位的风场格点的风力取值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810868066.0A CN109213759B (zh) | 2018-08-02 | 2018-08-02 | 一种大气风场急流线的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810868066.0A CN109213759B (zh) | 2018-08-02 | 2018-08-02 | 一种大气风场急流线的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109213759A true CN109213759A (zh) | 2019-01-15 |
CN109213759B CN109213759B (zh) | 2021-09-07 |
Family
ID=64987904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810868066.0A Expired - Fee Related CN109213759B (zh) | 2018-08-02 | 2018-08-02 | 一种大气风场急流线的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109213759B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110221359A (zh) * | 2019-06-11 | 2019-09-10 | 成都信息工程大学 | 一种自动识别切变线的方法 |
CN110263412A (zh) * | 2019-06-04 | 2019-09-20 | 南京大学 | 副热带急流或极锋急流强度和径向位置协同变化的表征方法 |
CN110413656A (zh) * | 2019-06-20 | 2019-11-05 | 南京信大气象科技有限公司 | 一种基于高空间分辨率数值模式的急流轴自动识别方法 |
CN111239852A (zh) * | 2020-01-13 | 2020-06-05 | 成都信息工程大学 | 一种利用风场数据自动识别高空横槽的方法 |
CN111399084A (zh) * | 2020-03-02 | 2020-07-10 | 中国电子科技集团公司第二十八研究所 | 一种基于三维风场数据的高空急流提取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102216811A (zh) * | 2008-11-12 | 2011-10-12 | 三菱重工业株式会社 | 用于创建气流场数据库的系统和用于预测扩散物质的扩散状况的系统 |
CN104699936A (zh) * | 2014-08-18 | 2015-06-10 | 沈阳工业大学 | 基于cfd短期风速预测风电场的扇区管理方法 |
US20160179994A1 (en) * | 2010-05-28 | 2016-06-23 | Geostellar, Inc. | System and Method for Geomatic Modeling of a Diverse Resource Base Across Broad Landscapes |
CN106844983A (zh) * | 2017-01-26 | 2017-06-13 | 厦门理工学院 | 一种提高建筑物防台风能力的方法 |
CN106919792A (zh) * | 2017-02-24 | 2017-07-04 | 天津大学 | 基于高精度数值风场资料的涡旋中心自动识别方法 |
CN107957321A (zh) * | 2016-10-17 | 2018-04-24 | 中国人民解放军装甲兵工程学院 | 一种高压风场模拟系统 |
-
2018
- 2018-08-02 CN CN201810868066.0A patent/CN109213759B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102216811A (zh) * | 2008-11-12 | 2011-10-12 | 三菱重工业株式会社 | 用于创建气流场数据库的系统和用于预测扩散物质的扩散状况的系统 |
US20160179994A1 (en) * | 2010-05-28 | 2016-06-23 | Geostellar, Inc. | System and Method for Geomatic Modeling of a Diverse Resource Base Across Broad Landscapes |
CN104699936A (zh) * | 2014-08-18 | 2015-06-10 | 沈阳工业大学 | 基于cfd短期风速预测风电场的扇区管理方法 |
CN107957321A (zh) * | 2016-10-17 | 2018-04-24 | 中国人民解放军装甲兵工程学院 | 一种高压风场模拟系统 |
CN106844983A (zh) * | 2017-01-26 | 2017-06-13 | 厦门理工学院 | 一种提高建筑物防台风能力的方法 |
CN106919792A (zh) * | 2017-02-24 | 2017-07-04 | 天津大学 | 基于高精度数值风场资料的涡旋中心自动识别方法 |
Non-Patent Citations (4)
Title |
---|
WENBO TANG 等: "Lagrangian coherent structures near a subtropical jet stream", 《JOURNAL OF THE ATMOSPHERIC SCIENCES》 * |
万瑜 等: "ECMWF细网格数值预报产品在乌鲁木齐东南大风预报中的释用", 《第32届中国气象学会年会S3 军用数值天气预报技术及应用》 * |
林振敏 等: "低涡、急流、切变线的微机自动识别", 《广西气象》 * |
王志诚等: "地基遥感大气温湿风垂直廓线观测方法综述 ", 《气象水文海洋仪器》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110263412A (zh) * | 2019-06-04 | 2019-09-20 | 南京大学 | 副热带急流或极锋急流强度和径向位置协同变化的表征方法 |
CN110263412B (zh) * | 2019-06-04 | 2023-04-07 | 南京大学 | 副热带急流或极锋急流强度和径向位置协同变化的表征方法 |
CN110221359A (zh) * | 2019-06-11 | 2019-09-10 | 成都信息工程大学 | 一种自动识别切变线的方法 |
CN110413656A (zh) * | 2019-06-20 | 2019-11-05 | 南京信大气象科技有限公司 | 一种基于高空间分辨率数值模式的急流轴自动识别方法 |
CN110413656B (zh) * | 2019-06-20 | 2021-05-18 | 南京信大气象科技有限公司 | 一种基于高空间分辨率数值模式的急流轴自动识别方法 |
CN111239852A (zh) * | 2020-01-13 | 2020-06-05 | 成都信息工程大学 | 一种利用风场数据自动识别高空横槽的方法 |
CN111399084A (zh) * | 2020-03-02 | 2020-07-10 | 中国电子科技集团公司第二十八研究所 | 一种基于三维风场数据的高空急流提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109213759B (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109213759A (zh) | 一种大气风场急流线的检测方法 | |
CN107330903B (zh) | 一种人体点云模型的骨架提取方法 | |
CN103810739B (zh) | 一种图像文字变形动画的生成方法 | |
CN109344874A (zh) | 一种基于深度学习的染色体自动分析方法及系统 | |
CN104463870A (zh) | 一种图像显著区域检测方法 | |
CN107516316A (zh) | 一种在fcn中引入聚焦机制对静态人体图像进行分割的方法 | |
CN104166999B (zh) | 一种基于地基云图强度分层的云团提取方法 | |
CN104899883A (zh) | 一种深度图像场景的室内对象立方体检测方法 | |
CN103218600B (zh) | 一种实时人脸检测算法 | |
CN104821005A (zh) | 可模拟笔势风格的汉字书法作品书写过程动态复现的方法 | |
CN102136133B (zh) | 一种图像处理方法及图像处理装置 | |
CN108875127B (zh) | 一种计算机气象软件中基于风场数据的槽线修正方法 | |
CN108961283A (zh) | 基于特征对应的图像变形方法及装置 | |
CN102129576A (zh) | 全天空极光图像占空比参数的提取方法 | |
CN106469138B (zh) | 词云的生成方法及装置 | |
CN105550682A (zh) | 钟鼎碑刻拓印方法 | |
CN103413306B (zh) | 一种自适应阈值的Harris角点检测方法 | |
CN106202852B (zh) | 一种植被生态系统气候敏感带类型的空间定量识别方法 | |
CN104346798B (zh) | 目标轮廓多边形逼近方法及其装置 | |
CN103093241B (zh) | 基于同质化处理的光学遥感图像非匀质云层判别方法 | |
CN109543700B (zh) | 一种抗遮挡的服饰关键点检测方法 | |
CN107862749A (zh) | 一种点云几何细节特征去除方法 | |
CN108830896A (zh) | 一种基于全卷积网络和建筑阴影几何模型的遥感影像水体斑块提取方法 | |
CN109117886B (zh) | 一种图像帧中目标尺度与区域估计的方法 | |
CN106780577A (zh) | 一种基于组特征的直线匹配方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210907 |