CN109697317B - 计算风向的装置及方法 - Google Patents
计算风向的装置及方法 Download PDFInfo
- Publication number
- CN109697317B CN109697317B CN201811577329.9A CN201811577329A CN109697317B CN 109697317 B CN109697317 B CN 109697317B CN 201811577329 A CN201811577329 A CN 201811577329A CN 109697317 B CN109697317 B CN 109697317B
- Authority
- CN
- China
- Prior art keywords
- wind direction
- value
- data
- wind
- angle
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
一种计算风向的装置及方法,该方法包括:采集或输入风向数据并存储;利用采集或输入的风向数据创建若干数据表,用于分别存储不同时间级别的风向角度数据;从所述若干数据表中调取一选择时间段内的若干风向角度,结合风向数计算并输出风向的方位值,所述方位值为将风向集合中的风向按照连续的角度范围依次排列时进行连续标记的序号。本发明在一个计算过程内即可得到风向的方位值,不用来回调用其他方法函数,可以计算出风向的方位值。
Description
技术领域
本发明涉及风向测量技术领域,尤其涉及一种计算风向的装置及方法。
背景技术
在大气的近地面层,由于风场结构的湍流性质,使风向风速成为在时间和空间上都不断变化的量。风向可以用16个方位来进行描述,也可用8个方位来进行描述。当我们做某项定性定量的研究时就需要对风向风速进行标准、准确的计算,只有在准确的风向下计算研究出的结果才具有可信性,准确性,客观性。下面列举两种平均风向的计算方法。
(1)算术平均法:
算术平均法是一项最简单的计算方法,算法公式如下:
式中Aa是简单的算术平均风向,Ai为第i个风向样本的风向方位角,N为样本个数。
(2)标量平均法
标量平均法是由美国环保局推荐使用的一种算法,算法公式如下:
算式中式中Aa是简单的算术平均风向,Ai为第i个风向样本的风向方位角,N为样本个数,当i=1时:Di=Ai;当i>1,δi<-180时:Di=Di-1+360+δi;当i>1,δi<|180|时Di=Di-1+δi;当i>1,δi>180时:Di=Di-1-360+δi;当i>1时:δi=Ai-Di-1。
目前针对大量的样本的平均风向,会利用上述的算法进行大量的繁琐的计算,且需要对每个计算结果进行依次比对才能去除不合理的结果,得到最后的方位值。
发明内容
有鉴于此,本发明的主要目的在于提供一种计算风向的装置及方法,以期至少部分地解决上述提及的技术问题中的至少之一。
为实现上述目的,本发明的技术方案如下:
作为本发明的一个方面,提供了一种计算风向的方法,包括以下步骤:
步骤A:采集或输入风向数据并存储;
步骤B:利用采集或输入的风向数据创建若干数据表,用于分别存储不同时间级别的风向角度数据;
步骤C:从所述若干数据表中调取一选择时间段内的若干风向角度,结合风向数计算并输出方位值,其中,将一风向集合中的风向按照连续的角度范围依次排列并对所述风向连续地标记序号,所述方位值能对应于由该风向角度确定的风向的序号,所述风向数为风向集合的长度。
作为本发明的另一个方面,提供了一种计算风向的装置,包括采集单元与输入/输出单元中的至少之一,以及处理单元,其中:
采集单元,用于采集风向数据并存储;
输入/输出单元,用于输入风向数据,输出计算得到的风向的方位值;
处理单元,用于对采集单元采集的或者输入/输出单元输入的风向数据进行处理,计算得到方位值之后再通过输入/输出单元输出;其中所述处理单元进一步包括存储表模块和计算模块,其中:
存储表模块,包括若干数据表,用于分别存储不同时间级别的风向角数据;
计算模块,用于从所述存储表模块中调取在一选择时间段内的若干风向角度,结合风向数计算得到方位值,其中,将一风向集合中的风向按照连续的角度范围依次排列并对所述风向连续地标记序号,所述方位值能对应于由该风向角度确定的风向的序号,所述风向数为风向集合的长度。
基于上述技术方案,本发明的有益效果在于:
(1)本发明将不同时间级别的风向角度数据存储在数据表中,根据选择的时间段调用相应的风向角数据,结合风向数可以直接计算风向的方位值,加快了计算进程。
(2)进一步通过精简计算过程,无需在计算过程中进行循环判断,较少了数据丢失,提高精确性。
(3)本发明在一个计算过程内即可计算得到风向的方位值,不用来回调用其他方法函数,减少来回往返的次数,提高了执行效率。
附图说明
图1是本发明实施例计算风向的方法的流程示意图;
图2是本发明实施例采用第一种方式计算风向的方位值的流程示意图;
图3是本发明实施例采用第二种方式计算风向的方位值的流程示意图;
图4是本发明实施例计算风向的装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。
本发明公开了一种计算风向的装置及方法,将不同时间级别的风向角度数据存储在数据表中,根据选择的时间段调用相应的风向角数据,结合风向数可以直接计算风向的方位值,加快了计算进程。
具体地,作为一示例性实施例,本发明提供了一种计算风向的方法,包括以下步骤:
步骤A:采集或输入风向数据并存储;
一般城市中分布有几百个空气超级检测站点,来监测环境中的风向风速情况,每个站点每分钟会上传一个含实时风向角度的数据包到服务器进行存储。
步骤B:利用采集或输入的风向数据创建若干数据表,用于分别存储不同时间级别的风向角度数据;该若干数据表包括分钟数据表、小时数据表以及日数据表,其中,分钟数据表用于存储监测站点采集的每分钟实时风向角度,小时数据表用于存储一小时内的每分钟平均风向角度,日数据表用于存储一日内每小时的平均风向角度。
具体地,首先创建分钟数据表,将上传的每分钟实时风向角度数据逐条地插入到分钟数据表中,然后创建小时数据表,小时数据作业例如可在每小时的第2分钟执行,根据分钟数据表的数据计算一小时内的每分钟平均风向角度插入到小时数据表中,然后创建日数据表,日数据作业例如可在每天的00:20分执行,根据小时数据表的数据计算一天内的每小时平均风向角度。可以理解,此处的平均风向角度可以采用现有的平均风向算法如算数平均法、标量平均法等进行计算。
步骤C:从若干数据表中调取一选择时间段内的若干风向角度,结合风向数计算并输出方位值,其中将一风向集合中的风向按照连续的角度范围依次排列并对所述风向连续地标记序号,方位值能对应于由该风向角度确定的风向的序号。
风向例如可以用4、8、16个方位来进行描述,如表1和表2中所示分别为16方位和8方位的风向符号和角度范围:
表1
表2
序号 | 缩写 | 风向 | 角度范围 |
0 | N | 北风 | 0-22.5;337.6-360 |
1 | EN | 东北风 | 22.6-67.5 |
2 | E | 东风 | 67.6-112.5 |
3 | ES | 东南风 | 112.6-157.5 |
4 | S | 南风 | 157.6-202.5 |
5 | WS | 西南风 | 202.6-247.5 |
6 | W | 西风 | 247.6-292.5 |
7 | WN | 西北风 | 292.6-337.5 |
以8方位的风向集合为例,该风向集合表示为{“东风”、“东南风”、“南风”、“西南风”、...“北风”},集合长度为8,也即风向数C为8。本实施例中以中心角度0°的风向作为起始风向,将各个风向按照连续的角度范围由小至大依次排列,并依次标记为序号0、1、…、7,如表2中所示。
本发明具体提供了两种方式来进行风向的方位值计算,其中第一种计算过程如图2所示,包括如下步骤:
子步骤C1:根据以下公式计算序号i,有i=(A+180°/C)/(360°/C),其中A为风向角度。子步骤C2:对i进行向下取整得到value值;子步骤C3:判断value值与序号i之间、以及风向角度与临界值之间的关系,对value值进行修正得到方位值。
更具体地,子步骤C3包括:若风向角度为临界值即t=1,value值与序号i相等且非零,则将value值减一,作为方位值,反之则判断value值是否等于风向数C,若value值等于风向数C,则将value值赋值为0,作为方位值,反之则直接将value值作为方位值。
举例来说,以风向角度A=180度为例,可以得到i=4.5,进行向下取整后得到value=4;由于风向角度A为非临界值,value值不等于8,因此方位值即为4。
第二种计算过程如图3所示,根据以下算法计算风向的方位值:
方位值=[(A+180°/C)/(360°/C)]%C,其中A为风向角度,%表示取余。可以理解取余的除数和被除数均为整数,此公式中的第二个“/”表示取整。
同样以风向角度180度为例,可以得到方位值=4%8=4,可以看出相较于第一种计算过程,无需循环计算判断,且无需进行向下取整的操作,精简了计算过程,减少了数据丢失,使得计算结果更为精确。
作为本发明的另一个方面,本发明还提供了一种计算风向的装置,包括采集单元与输入/输出单元中的至少之一,以及处理单元,其中:
采集单元,用于采集风向数据并存储;
输入/输出单元,用于输入风向数据,输出计算得到的风向的方位值;
处理单元,用于对采集单元采集的或者输入/输出单元输入的风向数据进行处理,计算得到方位值之后再通过输入/输出单元输出;其中处理单元进一步包括存储表模块和计算模块,如图4所示,其中:
存储表模块,包括若干数据表,用于分别存储不同时间级别的风向角数据;
计算模块,用于从存储表模块中调取在一选择时间段内的若干风向角度,结合风向数计算并输出方位值,其中,将一风向集合中的风向按照连续的角度范围依次排列并对所述风向连续地标记序号,方位值能对应于由该风向角度确定的风向的序号。
其中,采集单元例如可以包括各气象站、空气超级检测站点、对卫星气象云图等中的风向数据进行提取的装置等。
其中,输入/输出单元例如可以包括各种通讯模块或者数据接口,用于接收从外部输入的风向数据。
其中,处理单元例如是基于单片机、可编程门阵列(FPGA)、可编程逻辑控制器(PLA)、台式机、笔记本电脑、智能手机、服务器等可以运行程序的设备,具体的算法语言不限,只要能执行本发明所述的算法即可。
综上所述,本发明的计算风向的装置及方法在一个计算过程内即可得到风向的方位值,不用来回调用其他方法函数,就可以计算出风向的方位值。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种计算风向的方法,包括以下步骤:
步骤A:采集或输入风向数据并存储;
步骤B:利用采集或输入的风向数据创建若干数据表,用于分别存储不同时间级别的风向角度数据;
步骤C:从所述若干数据表中调取一选择时间段内的若干风向角度,结合风向数计算并输出风向的方位值,其中,以中心角度0°的风向作为起始风向,将一风向集合中的风向按照连续的角度范围由小至大依次排列并对所述风向连续地标记序号0、1、…、C-1,所述方位值能对应于由该风向角度确定的风向的序号,所述风向数为风向集合的长度;
其中,所述方位值通过以下算法计算得到:根据以下公式计算序号i,有i=(A+180°/C)/(360°/C),其中A为风向角度;对i进行向下取整得到value值;判断value值与序号i之间、以及风向角度与临界值之间的关系,对value值进行修正得到方位值;
或者所述方位值通过以下算法计算得到:[(A+180°/C)/(360°/C)]%C,其中A为风向角度,%表示取余。
2.根据权利要求1所述的方法,其特征在于,步骤B中,所述若干数据表包括分钟数据表、小时数据表以及日数据表,其中,所述分钟数据表用于存储监测站点采集的每分钟实时风向角度,所述小时数据表用于存储一小时内的每分钟平均风向角度,所述日数据表用于存储一日内每小时的平均风向角度。
3.根据权利要求1所述的方法,其特征在于,所述判断value值与序号i之间、以及风向角度与临界值之间的关系,对value值进行修正得到方位值具体包括:
若风向角度为临界值,value值与序号i相等且非零,则将value值减一,作为方位值,反之则判断value值是否等于风向数C,若value值等于风向数C,则将value值赋值为0,作为方位值,反之则直接将value值作为方位值。
4.一种计算风向的装置,包括采集单元与输入/输出单元中的至少之一,以及处理单元,其特征在于:
采集单元,用于采集风向数据并存储;
输入/输出单元,用于输入风向数据,输出计算得到的风向的方位值;
处理单元,用于对所述采集单元采集的或者所述输入/输出单元输入的风向数据进行处理,计算得到方位值之后再通过输入/输出单元输出;其中所述处理单元进一步包括:
存储表模块,包括若干数据表,用于分别存储不同时间级别的风向角数据;
计算模块,用于从所述存储表模块中调取在一选择时间段内的若干风向角度,结合风向数计算得到方位值,其中,以中心角度0°的风向作为起始风向,将一风向集合中的风向按照连续的角度范围由小至大依次排列并对所述风向连续地标记序号0、1、…、C-1,所述方位值能对应于由该风向角度确定的风向的序号,所述风向数为风向集合的长度;
其中,所述方位值通过以下算法计算得到:根据以下公式计算序号i,有i=(A+180°/C)/(360°/C),其中A为风向角度;对i进行向下取整得到value值;判断value值与序号i之间、以及风向角度与临界值之间的关系,对value值进行修正得到方位值;
或者所述方位值通过以下算法计算得到:[(A+180°/C)/(360°/C)]%C,其中A为风向角度,%表示取余。
5.根据权利要求4所述的装置,其特征在于,所述若干数据表包括分钟数据表、小时数据表以及日数据表,其中,所述分钟数据表用于存储监测站点采集的每分钟实时风向角度,所述小时数据表用于存储一小时内的每分钟平均风向角度,所述日数据表用于存储一日内每小时的平均风向角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811577329.9A CN109697317B (zh) | 2018-12-21 | 2018-12-21 | 计算风向的装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811577329.9A CN109697317B (zh) | 2018-12-21 | 2018-12-21 | 计算风向的装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109697317A CN109697317A (zh) | 2019-04-30 |
CN109697317B true CN109697317B (zh) | 2020-12-01 |
Family
ID=66232085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811577329.9A Active CN109697317B (zh) | 2018-12-21 | 2018-12-21 | 计算风向的装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109697317B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780104A (zh) * | 2016-11-18 | 2017-05-31 | 中国电建集团成都勘测设计研究院有限公司 | 一种基于概率统计的平均风向计算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6426718B2 (ja) * | 2013-05-20 | 2018-11-21 | プリンシプル・パワー・インコーポレーテツド | オフショア浮体式風力タービン・プラットフォームを制御するシステムおよび方法 |
CN103500370B (zh) * | 2013-10-21 | 2016-08-17 | 华北电力大学 | 动态风电场风向坐标预计算系统建立方法 |
CN104021309B (zh) * | 2014-06-20 | 2016-10-12 | 中国能源建设集团江苏省电力设计院有限公司 | 基于方位偏转法的测风数据风向订正方法 |
CN104763585B (zh) * | 2015-02-06 | 2017-12-22 | 湘潭大学 | 基于分布式数据收集的风电机组动态重构方法 |
CN107609280B (zh) * | 2017-09-15 | 2021-03-02 | 深圳市卡普瑞环境科技有限公司 | 一种时间区间的确定方法以及终端 |
-
2018
- 2018-12-21 CN CN201811577329.9A patent/CN109697317B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780104A (zh) * | 2016-11-18 | 2017-05-31 | 中国电建集团成都勘测设计研究院有限公司 | 一种基于概率统计的平均风向计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109697317A (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109580003B (zh) | 一种静止气象卫星热红外数据估算近地面大气温度方法 | |
CN111538957B (zh) | 大气污染物来源贡献度的获取方法、装置、设备及介质 | |
CN111126662B (zh) | 基于大数据的灌溉决策制定方法、装置、服务器和介质 | |
CN111551964B (zh) | 一种基于导航卫星系统的水汽观测方法 | |
CN111950396B (zh) | 一种仪表读数神经网络识别方法 | |
CN110390343B (zh) | 一种空间气象数据的订正方法及系统 | |
CN112417734B (zh) | 一种基于风电场地理信息的风速订正方法及装置 | |
CN116245230B (zh) | 一种电站设备的运行巡检与趋势分析方法及系统 | |
CN115239105A (zh) | 一种在役风电场的风资源的评估方法以及装置 | |
CN106845785A (zh) | 一种基于多目标函数的多星测控资源调度方法 | |
CN115759365A (zh) | 光伏发电功率的预测方法及相关设备 | |
CN109697317B (zh) | 计算风向的装置及方法 | |
CN112985867B (zh) | 舵机测试方法、装置、设备及存储介质 | |
LU503778B1 (en) | The method for wind speed downscaling of wind farms based on super-resolution reconstruction technology | |
CN108021774B (zh) | 数据处理方法以及装置 | |
AU2023202496A1 (en) | An on-line monitoring method and system for fugitive carbon emission in cement industrial plant area | |
CN107607913B (zh) | 基于对数累积量的海杂波Pareto分布参数估计方法 | |
CN112163695A (zh) | 基于copula函数的风力光伏发电预测方法、系统和介质 | |
CN108196087B (zh) | 数据处理装置 | |
CN113159426A (zh) | 天气型相似性判断方法、装置、电子设备及可读存储介质 | |
CN113094638A (zh) | 一种气象网格化数据的多源场融合方法和系统 | |
CN111027019B (zh) | 统计分析高精度定位结果置信度的方法及装置 | |
CN111540418A (zh) | 一种植物中砷过量的概率值的预测方法及系统 | |
CN110705924A (zh) | 一种基于风向扇区的测风塔测风数据处理方法及装置 | |
CN116468205B (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 |