前提・実現したいこと
以下のようなcsvファイルを読み取ってデータをまとめる作業を行おうと考えています。
ダウンロードした時刻:2019/11/11 16:04:33
寺泊 寺泊
年 月 日 時 風速(m/s) 風向
2016 12 23 5 9.6 南南東
2016 12 23 6 9.8 南東
2016 12 23 7 10.6 東
2016 12 23 8 10.4 北東
2016 12 23 9 10.5 西
2016 12 23 10 9.1 西
2016 12 23 11 8.1 西北西
2016 12 23 12 7.8 西北西
2016 12 23 13 6.8 西北西
2016 12 23 14 6.3 北西
2016 12 23 15 6.2 北西
2016 12 23 16 6.5 北北西
2016 12 23 17 6.3 北北西
2016 12 23 18 6.3 北
2016 12 23 19 5.4 南西
2016 12 23 20 3.9 西南西
2016 12 23 21 4 南南西
2016 12 23 22 4.6 南
2016 12 23 23 4 南
処理の方法は
➀風速の値を2乗する。
➁次に角度補正を行いたいため、風向によってcosθを風速の2乗の値にかける
ex)上記のデータの場合
・風向は16方位に区切られているので、隣り合う方位の角度の間隔は22.5°
西を正面とし
西北西と西南西は22.5°、北西と南西は45°、北北西と南南西は67.5°
・この西を基準とした角度を基に、風速の2乗値にそれぞれcosθをかける
(10.5)^2×cos(0)+(9.1)^2×cos(0)+
(8.1)^2×cos(22.5)+(7.8)^2×cos(22.5)+(6.8)^2×cos(22.5)+
(6.3)^2×cos(45)+(6.2)^2×cos(45)+
(6.5)^2×cos(67.5)+(6.3)^2×cos(67.5)+
(5.4)^2×cos(45)+
(3.9)^2×cos(22.5)+
(4)^2×cos(67.5)
➂ ➁の値で、特定の風向(今回の場合'南南西','南西','西南西','西','西北西','北西','北北西')を満たすもの全部足し合わせる。
このような処理を行いたいと考えています。
発生している問題・エラーメッセージ
現在特定の風向を満たす、風速の値を2乗して、それを足し合わせる
というところまでは動きましたが
➁のような、ある風向にはcos(22.5)をかけ、ある風向にはcos(45)をかけるといった処理ができずに悩んでおります。
該当のソースコード
今現在のプログラムを以下に示します。
長い文章になって申し訳ございません。
そこで今現在このような流れで動かしていこうかなと思っています。
以下に記します。
import csv
import math
cos0 = math.cos(math.radians(0))
cos225 = math.cos(math.radians(22.5))
cos45 = math.cos(math.radians(45))
cos675 = math.cos(math.radians(67.5))
f = open("test.csv", "r")
reader = csv.reader(f)
next(reader); next(reader); next(reader); next(reader); next(reader)
Direction = ['南南西','南西','西南西','西','西北西','北西','北北西']
sum = 0
for row in reader:
if not row[4].strip(): # 空行スキップ
continue
colE = float(row[4])
#この先からつまずいています。
この後に
➀ colEを全て二乗
➁ Directionの7方位で
方位が西のcolE(2乗したもの)にcos0を掛ける
方位が西南西・西北西のcolE(2乗したもの)にcos225を掛ける
方位が南西・北西のcolE(2乗したもの)にcos45を掛ける
方位が南南西・北北西のcolE(2乗したもの)にcos675を掛ける
➂2乗値を補正した値を足し合わせる
といった流れでできないかと考えています。