[matplotlib 3D] 28. Pythonで3D voxelグラフ

はじめに

matplotlibのmplot3dでボクセルグラフを表示する方法について解説する。

3Dボクセルグラフは、データを立方体で表示するグラフのことであり、２次元のピクセルに対して、3次元のボクセルということになっている。[1]

np.indicesで各要素のindexの値をもつ配列ができる。例えば、np.indices((2, 2, 2))とすれば、

x= [[[0 0]  [0 0]] [[1 1]  [1 1]]]
y= [[[0 0]  [1 1]] [[0 0]  [1 1]]]
z= [[[0 1]  [0 1]] [[0 1]  [0 1]]]

となる。[2]

8×8×8の配列でcube1の条件に当てはまるとこだけ、Trueとなる。cube2, linkも同様。

で各配列を比較して、Trueの場合はTrueになる。[3]

np.emptyで適当な同じ形の配列をつくって、linkの場合、linkがTrueの部分を’red’とする。他も同様に色を設定する。

ax.voxels(voxels)でボクセルグラフを表示できる。
facecolorで色を設定して、edgecolorを黒としてプロットした。

defaultでshade=Trueなので、影無しにしたい場合、shade=Falseとする必要がある。