web-dev-qa-db-ja.com

Pythonで回転できる3D散布図をプロットする簡単な方法はありますか?

現在、私はmatplotlibを使用して3D散布図をプロットしていますが、3D散布図を作成している間、それを回転させてデータを見やすくする方法を見つけることができません。

次に例を示します。

import pylab as p
import mpl_toolkits.mplot3d.axes3d as p3

#data is an ndarray with the necessary data and colors is an ndarray with
#'b', 'g' and 'r' to Paint each point according to its class

...

fig=p.figure()
ax = p3.Axes3D(fig)
ax.scatter(data[:,0], data[:,2], data[:,3], c=colors)
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')
fig.add_axes(ax)
p.show()

実行時に実行できる解決策が欲しいのですが、回転させることができ、それが短い/迅速である限り、問題はありません。

PCAを虹彩データセットに適用した後に生成されたプロットの比較を次に示します。
1。マヤビ
alt text
2。 matplotlib
alt text

Mayaviを使用するとデータの視覚化が容易になりますが、MatPlotLibはよりプロフェッショナルに見えます。 Matplotlibも軽量です。

22
pnodbnda

mayavi を使用すると、このようなプロットを作成できます。

import enthought.mayavi.mlab as mylab
import numpy as np
x, y, z, value = np.random.random((4, 40))
mylab.points3d(x, y, z, value)
mylab.show()

GUIでは、クリックアンドドラッグによる回転、および右クリックアンドドラッグによるズームイン/ズームアウトが可能です。

alt text

14
unutbu

さて、最初に「私のデータをよりよく見る」とはどういう意味かを定義する必要があります...

インタラクティブに作業したい場合は、マウスを使用してプロットを回転および拡大できます。

プログラムで軸を回転させたいだけの場合は、ax.view_init(elev, azim)を使用します。ここで、elevazimは、表示する仰角と方位角(度単位)です。あなたの陰謀。

または、ax.elevax.azim、およびax.dist現在の視点の標高、方位角、距離を取得/設定するためのプロパティ。

ソース から この例 を借りる:

import numpy as np
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt

def randrange(n, vmin, vmax):
    return (vmax-vmin)*np.random.Rand(n) + vmin

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
n = 100
for c, m, zl, zh in [('r', 'o', -50, -25), ('b', '^', -30, -5)]:
    xs = randrange(n, 23, 32)
    ys = randrange(n, 0, 100)
    zs = randrange(n, zl, zh)
    ax.scatter(xs, ys, zs, c=c, marker=m)

ax.set_xlabel('X Label')
ax.set_ylabel('Y Label')
ax.set_zlabel('Z Label')

plt.show()

いい散布図が得られます。 alt text

以下のように、プログラムで軸を回転できます。

import numpy as np
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt

def randrange(n, vmin, vmax):
    return (vmax-vmin)*np.random.Rand(n) + vmin

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
n = 100
for c, m, zl, zh in [('r', 'o', -50, -25), ('b', '^', -30, -5)]:
    xs = randrange(n, 23, 32)
    ys = randrange(n, 0, 100)
    zs = randrange(n, zl, zh)
    ax.scatter(xs, ys, zs, c=c, marker=m)

ax.set_xlabel('X Label')
ax.set_ylabel('Y Label')
ax.set_zlabel('Z Label')

ax.azim = 200
ax.elev = -45

plt.show()

alt text

少しお役に立てば幸いです。

29
Joe Kington