[matplotlib 3D] 18. 3D surfaceプロットで球を表示

はじめに

matplotlibのmplot3dをつかって、3Dグラフ上にsurface plotで球を表示する方法について解説する。

uは0から2πの数列、vは0からπの数列。いずれも100個ある。
x,y,zは球面の方程式を極座標で表した場合で、

$$x=rcos\theta sin\phi, y=rsin\theta sin\phi, z=rcos\phi$$

となる。
なお直交座標系で半径rの球面は、

$$x^2+y^2+z^2=r^2$$

で表される。

データとしては真球であるのに対して、見た目が楕円球となっているので図のアスペクト比を揃えることを試みた。最初に下記のコードを実行すると、

ax.set_aspect('equal')

でエラー(NotImplementedError)となり、うまく行かない。そこでfigsizeでアスペクト比を揃えてグラフを表示することで真球を表示することができた。

matplotlib3.3以降では、set_box_aspectにより3Dグラフのアスペクト比を変えることができる。

[…] [matplotlib 3D] 18. 3D surfaceプロットで球を表示matplotlibのmplot3dをつかって、3Dグラフ上にsurface plotで球を表示する方法について解説する。sabopy.com2018.10.17 […]

返信

返信