[matplotlib 3D] 35. 円柱(円筒座標)

はじめに

matplotlib mplot3dで円柱を表示する方法を解説する。ここでは円筒座標系でデータを生成する。

figを作成して、fig.gca(projection=’3d’)で3Dグラフにする。

半径2, 高さ10の円柱にする。
theta,zでmeshgridを作成して[1]、X, Yは円筒座標なので、

$$x=rcos\theta\\y=rsin\theta$$

とする[2]。

color=’g’で緑色、alpha=0.5で透明度0.5のsurfaceグラフを作成する。

Circleで半径2の中心(0,0)の円を作成する。
ax.path_patchでグラフに追加する。
pathpatch_2d_to_3dで3dグラフ上に2dの円が表示可能となる。[3]
表示する場所は、蓋なのでz=10, zdir=”z”とする。

蓋と同様に底も作成する。
蓋の位置は、z=0, zdir=”z”となる。

ax.set_box_aspect((1,1,1))でアスペクト比を揃えている。

作製した図を回転させると以下のようになる。

回転アニメーションは以下の記事を参考にして作成した。