2次元リストの使い方(python)

Question

###Hを作成するプログラム

orient = np.random.randint(0, 360, 100000)
speed = np.random.randint(0, 900, 100000)

x_edges = np.arange(0, 901, 100)
y_edges = np.linspace(0, 360, 17)
H = np.histogram2d(np.mod(orient + 11.25, 360), speed,
                   bins=(y_edges, x_edges))[0].astype(int)
print(H.shape)  # (16, 9)

sns.set()
# 目盛りのラベルを設定する。
xlabels = ['{} ~ {}'.format(start, end)
           for start, end in zip(x_edges, x_edges[1:])]
ylabels = ['N', 'NNW', 'NW', 'WNW',
           'W', 'WSW', 'SW', 'SSW',
           'S', 'SSE', 'SE', 'ESE',
           'E', 'ENE', 'NE', 'NNE']
# ヒートマップを作成する。
sns.heatmap(H, annot=True, cmap='Reds', fmt='d',
            xticklabels=xlabels, yticklabels=ylabels, cbar=False)
plt.show()

前提・実現したいこと

上図のような１６×９の配列Hを利用して、１マスずつ以下のような計算を行い、計算結果を同じように１６×９の配列に収納したいと考えております。
下のプログラムを実行したところvel.append(ve[m,n])の行にinvalid index to scalar variable.のエラーが帰ってきました。これは１次元配列に２次元配列を代入しようとしたためと考えましたが、この捉え方であっているのでしょうか。その場合二次元の空配列のリストを作る方法などはあるのでしょうか。
もしくは、リストを作成する以外にいい方法があれば教えていただきたいです。（その後のプログラムで計算結果を行ごとにグラフにしたいため、配列のように列や行で扱いやすいものが良いです。）

発生している問題・エラーメッセージ

#流向の割合を求める

dir = list()
dir = []

for i in range(16):
    di=sum(H[i,:])/len(orient)*100
    dir.append(di)
#print(dir)

#流速の割合を求める
vel = list()
vel =[]
for m in range (16):
    for n in range (9):
        ve=H[n,m]/sum(H[m,:])*dir[m]
        vel.append(ve[m,n])
#print(vel)

Comments

[1T2R3M4]

H[n,m]

[quickquip]

想像して答えることはできるのでしょうけど、回答者から見ると「H の正体(定義とか、型とか)が書いてない状態でどうすればいいのか?」という質問になっちゃってますよ。

Answer 1

配列Hに入るデータが「すべて数値、必要な精度が一緒、可変長配列を使わない」ならnumpyを使うことを強く勧めます。pythonのnumpyモジュールはとても強力で計算が早くなる工夫がいっぱいあります。

numpyはいろいろなところで解説があるので見てください。

Answer 2

        ve=H[n,m]/sum(H[m,:])*dir[m]
        vel.append(ve[m,n])

Hが2D arrayとして、veはスカラーの値になりますね（型はたぶんnp.float64かな）。

スカラーの[m,n]は取れないでしょう・・・。

Answer 3

veはスカラー値です。
それに対して配列操作ve[m,n]を試みているので提示エラーが発生しています。