回答編集履歴
1
d
test
CHANGED
|
@@ -97,3 +97,67 @@
|
|
|
97
97
|
|
|
98
98
|
|
|
99
99
|
|
|
100
|
+
|
|
101
|
+
|
|
102
|
+
|
|
103
|
+
## 追記
|
|
104
|
+
|
|
105
|
+
|
|
106
|
+
|
|
107
|
+
速度ごとの各方角の割合を出したい場合は以下のようにしてください。
|
|
108
|
+
|
|
109
|
+
|
|
110
|
+
|
|
111
|
+
```python
|
|
112
|
+
|
|
113
|
+
# 速度の割合
|
|
114
|
+
|
|
115
|
+
sum_along_speed = H.sum(axis=0) # 速度ごとの総和
|
|
116
|
+
|
|
117
|
+
orient_ratio = np.divide(
|
|
118
|
+
|
|
119
|
+
H, sum_along_speed,
|
|
120
|
+
|
|
121
|
+
out=np.zeros(H.shape, dtype=float),
|
|
122
|
+
|
|
123
|
+
where=sum_along_speed != 0)
|
|
124
|
+
|
|
125
|
+
|
|
126
|
+
|
|
127
|
+
# H.sum(axis=0) は列 (速度) ごとの総和なので、形状は (速度の数,)
|
|
128
|
+
|
|
129
|
+
# H の形状は (方角の数, 速度の数) なので、numpy のブロードキャストにより、各速度ごとの方向の割合が出る。
|
|
130
|
+
|
|
131
|
+
sns.heatmap(orient_ratio, annot=True, cmap='Reds', fmt='.1%',
|
|
132
|
+
|
|
133
|
+
xticklabels=xlabels, yticklabels=ylabels, cbar=False)
|
|
134
|
+
|
|
135
|
+
plt.show()
|
|
136
|
+
|
|
137
|
+
```
|
|
138
|
+
|
|
139
|
+
|
|
140
|
+
|
|
141
|
+
|
|
142
|
+
|
|
143
|
+
|
|
144
|
+
|
|
145
|
+
```
|
|
146
|
+
|
|
147
|
+
np.divide(
|
|
148
|
+
|
|
149
|
+
H, sum_along_speed,
|
|
150
|
+
|
|
151
|
+
out=np.zeros(H.shape, dtype=float),
|
|
152
|
+
|
|
153
|
+
where=sum_along_speed != 0)
|
|
154
|
+
|
|
155
|
+
```
|
|
156
|
+
|
|
157
|
+
|
|
158
|
+
|
|
159
|
+
0割り (division by zero) が発生すると、数学ではこのような演算は定義されないので、結果は NaN になります。
|
|
160
|
+
|
|
161
|
+
そのため、`where=sum_along_speed != 0` で割る数`sum_along_speed`の値が0でない要素だけ除算するようにします。
|
|
162
|
+
|
|
163
|
+
割る数が0の場合は、結果も0とするので、`out=np.zeros(H.shape, dtype=float)`としています。
|