はじめに
skimage.morphologyのoctahedronで、大きさの異なる八面体の構造化要素の図をadd_subplotで作成し、それぞれの図の回転アニメーションを表示する。
コード
解説
モジュールのインポートなど
バージョン
構造化要素の生成
辞書形式で大きさの異なる八面体を作成する。
octahedron(5)で半径が5の八面体を作成できる。
構造化要素の表示
最初に空のリスト(axes)を作成する。
plt.figure(figsize=(12,8))でfigを作成し、fig.add_subplot(2, 3, idx)で図を順次生成していく。projection=Axes3D.nameとすることで3Dグラフを表示できる。
ax.voxels(struc,ec=’w’)によりボクセルグラフで立方体を表示する。ecで各セルのエッジカラーを設定でき、ここではwで白とした。fcで各セルの面の色を設定できる。
なお、ボクセルグラフについては下記で解説した。
[matplotlib 3D] 28. Pythonで3D voxelグラフ
3D voxelグラフの作成について解説する。
ax.set(xlim=(-1,struc.shape[0]+1)などで各構造化要素の大きさに応じた軸範囲設定を行う。
作成した図はaxes.append(ax)でリストに加える。
idx += 1によりidxに1を足して次の図の作成を行う。
アニメーションの関数
view_initをリスト内の各要素に適用して図の見える角度を変化させて回転アニメーションとする。
アニメーションの表示
FuncAnimationでアニメーションを表示する。frame数はデフォルトで100なのでintervalが100msの場合、10秒のアニメーションとなる。HTML(ani.to_html5_video())
により、jupyter notebook またはjupyter lab上にアニメーションを表示できる。
参考
Generate footprints (structuring elements) — skimage 0.25.0rc2.dev0 documentation
skimage.morphology — skimage 0.24.1rc0.dev0 documentation
[matplotlib 3D] 14. 3Dグラフの回転アニメーション
matplotlib mplot3d のグラフの回転アニメーション
コメント