はじめに
plot_trisurfは、非構造化データのsurface plotを表示するための関数です。通常のsurface plotは格子状に整列したデータを必要としますが、plot_trisurfは不規則に分布するデータポイントからなめらかな表面を生成できます。
コード

解説
モジュールのインポート
データの生成
rが0だとデータがダブっているのがたくさんできるので0.125からの始まる配列とする。endpoint=Falseで2πを含まない数列、[…, np.newaxis]で次元を追加している。[…, np.newaxis]がないとraddi*np.cos(angles)の計算ができない。
x = np.append(0, (radii*np.cos(angles)).flatten())
np.appendで配列を追加。raddiとnp.cos(angles)をかけることで配列をブロードキャストし、flattenで一次元化しているので、8*36 = 288個の配列となる。0をいれると289個となる。
trisurf plotの作成
plot_trisurfは内部でDelaunay三角形分割を使用して、点群から三角形メッシュを生成します。これにより、規則的なグリッドを持たないデータでも滑らかな表面を作成できます。plot_trisurfを使用すると三角形のsurfaceグラフを作成できます。この関数ではx, y, zはmeshgridのような構造化データではなく、単純な1次元配列として扱われます。
z=0とした場合

回転アニメーション
作製した図を回転させると以下のようになります。
回転アニメーションは以下の記事を参考にして作成した。

まとめ
plot_trisurfは、非構造化データや不規則に分布するデータポイントから3D表面プロットを生成するための強力なツールです。特に地形データや科学計算の結果など、規則的なグリッドを持たないデータの可視化に適しています。
参考


コメント