CN109213759B - 一种大气风场急流线的检测方法 - Google Patents
一种大气风场急流线的检测方法 Download PDFInfo
- Publication number
- CN109213759B CN109213759B CN201810868066.0A CN201810868066A CN109213759B CN 109213759 B CN109213759 B CN 109213759B CN 201810868066 A CN201810868066 A CN 201810868066A CN 109213759 B CN109213759 B CN 109213759B
- Authority
- CN
- China
- Prior art keywords
- point
- grid
- wind
- points
- skeleton
- 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.)
- Expired - Fee Related
Links
Images
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、根据公式
进一步地,步骤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、根据公式
步骤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、根据公式
3.根据权利要求2所述的大气风场急流线的检测方法,其特征在于:所述步骤S3中的相干性阈值为π2/18;其中π为常数。
4.根据权利要求3所述的大气风场急流线的检测方法,其特征在于:所述步骤S4中风速阈值为16m/s。
5.根据权利要求4所述的大气风场急流线的检测方法,其特征在于:所述步骤S6中九宫格的外围格子的赋值方法为:
将九宫格的外围格子从左上角开始按逆时针方向依次赋值为1、2、4、16、128、64、32、8。
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 CN109213759A (zh) | 2019-01-15 |
CN109213759B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110263412B (zh) * | 2019-06-04 | 2023-04-07 | 南京大学 | 副热带急流或极锋急流强度和径向位置协同变化的表征方法 |
CN110221359B (zh) * | 2019-06-11 | 2020-06-05 | 成都信息工程大学 | 一种自动识别切变线的方法 |
CN110413656B (zh) * | 2019-06-20 | 2021-05-18 | 南京信大气象科技有限公司 | 一种基于高空间分辨率数值模式的急流轴自动识别方法 |
CN111239852B (zh) * | 2020-01-13 | 2021-09-07 | 成都信息工程大学 | 一种利用风场数据自动识别高空横槽的方法 |
CN111399084B (zh) * | 2020-03-02 | 2021-12-24 | 中国电子科技集团公司第二十八研究所 | 一种基于三维风场数据的高空急流提取方法 |
Citations (5)
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短期风速预测风电场的扇区管理方法 |
CN106844983A (zh) * | 2017-01-26 | 2017-06-13 | 厦门理工学院 | 一种提高建筑物防台风能力的方法 |
CN106919792A (zh) * | 2017-02-24 | 2017-07-04 | 天津大学 | 基于高精度数值风场资料的涡旋中心自动识别方法 |
CN107957321A (zh) * | 2016-10-17 | 2018-04-24 | 中国人民解放军装甲兵工程学院 | 一种高压风场模拟系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2577532A4 (en) * | 2010-05-28 | 2016-04-06 | Geostellar Inc | SYSTEM AND METHOD FOR BASIC GEOMATIC MODELING OF VARIOUS RESOURCES ON LARGE LANDSCAPES |
-
2018
- 2018-08-02 CN CN201810868066.0A patent/CN109213759B/zh not_active Expired - Fee Related
Patent Citations (5)
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短期风速预测风电场的扇区管理方法 |
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 |
---|
ECMWF细网格数值预报产品在乌鲁木齐东南大风预报中的释用;万瑜 等;《第32届中国气象学会年会S3 军用数值天气预报技术及应用》;20151014;第274-286页 * |
Lagrangian coherent structures near a subtropical jet stream;Wenbo Tang 等;《Journal of the Atmospheric Sciences》;20100701;第2307–2319页 * |
低涡、急流、切变线的微机自动识别;林振敏 等;《广西气象》;19971230;第47-48页 * |
地基遥感大气温湿风垂直廓线观测方法综述;王志诚等;《气象水文海洋仪器》;20180615(第02期);第109-116页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109213759A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109213759B (zh) | 一种大气风场急流线的检测方法 | |
CN110781967A (zh) | 一种基于可微分二值化的实时文本检测方法 | |
CN101315663B (zh) | 一种基于区域潜在语义特征的自然场景图像分类方法 | |
CN110222769B (zh) | 一种基于YOLOV3-tiny的改进目标检测方法 | |
CN111179217A (zh) | 一种基于注意力机制的遥感图像多尺度目标检测方法 | |
CN104463870A (zh) | 一种图像显著区域检测方法 | |
CN102915541B (zh) | 多尺度的图像分割方法 | |
CN102831604A (zh) | 一种灰度图像的二维Renyi熵阈值分割方法 | |
CN107341470A (zh) | 一种基于航拍图像的输电线检测方法 | |
CN110097598A (zh) | 一种基于pvfh特征的三维物体位姿估计方法 | |
CN111915558A (zh) | 一种高压输电线销钉状态检测方法 | |
CN116452850A (zh) | 一种基于数据挖掘和深度学习的道路积水区域识别方法 | |
CN106204719A (zh) | 基于二维邻域检索的三维场景中海量模型实时调度方法 | |
CN111368865B (zh) | 遥感影像储油罐检测方法、装置、可读存储介质及设备 | |
Huang et al. | Survey of target detection algorithms in SAR images | |
CN112101113B (zh) | 一种轻量化的无人机图像小目标检测方法 | |
CN113936041A (zh) | 一种面向光伏发电的云团跟踪及运动趋势预测方法 | |
CN109063834A (zh) | 一种基于卷积特征响应图的神经网络剪枝方法 | |
CN108615240B (zh) | 一种结合邻域信息与距离权重的非参贝叶斯过分割方法 | |
CN110910417A (zh) | 一种基于超像素邻帧特征对比的弱小运动目标检测方法 | |
Jin et al. | Filtering processing of LIDAR point cloud data | |
CN113920273B (zh) | 图像处理方法、装置、电子设备和存储介质 | |
Huang et al. | Multi-scale feature pair based R-CNN method for defect detection | |
CN103745493B (zh) | 一种基于图像处理的二维内阴影生成方法 | |
CN113256642B (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 |
Granted publication date: 20210907 |
|
CF01 | Termination of patent right due to non-payment of annual fee |