質問編集履歴
1
説明の追加
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -1,4 +1,28 @@
|
|
|
1
1
|
|
|
2
|
+
###Hを作成するプログラム
|
|
3
|
+
```
|
|
4
|
+
orient = np.random.randint(0, 360, 100000)
|
|
5
|
+
speed = np.random.randint(0, 900, 100000)
|
|
6
|
+
|
|
7
|
+
x_edges = np.arange(0, 901, 100)
|
|
8
|
+
y_edges = np.linspace(0, 360, 17)
|
|
9
|
+
H = np.histogram2d(np.mod(orient + 11.25, 360), speed,
|
|
10
|
+
bins=(y_edges, x_edges))[0].astype(int)
|
|
11
|
+
print(H.shape) # (16, 9)
|
|
12
|
+
|
|
13
|
+
sns.set()
|
|
14
|
+
# 目盛りのラベルを設定する。
|
|
15
|
+
xlabels = ['{} ~ {}'.format(start, end)
|
|
16
|
+
for start, end in zip(x_edges, x_edges[1:])]
|
|
17
|
+
ylabels = ['N', 'NNW', 'NW', 'WNW',
|
|
18
|
+
'W', 'WSW', 'SW', 'SSW',
|
|
19
|
+
'S', 'SSE', 'SE', 'ESE',
|
|
20
|
+
'E', 'ENE', 'NE', 'NNE']
|
|
21
|
+
# ヒートマップを作成する。
|
|
22
|
+
sns.heatmap(H, annot=True, cmap='Reds', fmt='d',
|
|
23
|
+
xticklabels=xlabels, yticklabels=ylabels, cbar=False)
|
|
24
|
+
plt.show()
|
|
25
|
+
```
|
|
2
26
|
### 前提・実現したいこと
|
|
3
27
|
|
|
4
28
|
上図のような１６×９の配列Hを利用して、１マスずつ以下のような計算を行い、計算結果を同じように１６×９の配列に収納したいと考えております。
|