Любой простой способ построить 3D-разброс в Python, который я могу вращать? - PullRequest
Новая
       33

Любой простой способ построить 3D-разброс в Python, который я могу вращать?

18 голосов
/ 19 января 2011

В настоящее время я использую matplotlib для построения трехмерного рассеяния, и пока он выполняет свою работу, я не могу найти способ повернуть его, чтобы лучше видеть мои данные.

Вот пример: 

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 также легче.

Ответы [ 2 ]

27 голосов
/ 19 января 2011

Ну, во-первых, вам нужно определить, что вы подразумеваете под «лучше видеть мои данные» ...

Вы можете вращать и увеличивать график, используя мышь, если вы хотите работать в интерактивном режиме.

Если вы просто хотите программно поворачивать оси, используйте ax.view_init(elev, azim), где elev и azim - это углы возвышения и азимута (в градусах), с которых вы хотите просмотреть свой график.

Кроме того, вы можете использовать свойства ax.elev, ax.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

Надеюсь, это немного поможет!

11 голосов
/ 19 января 2011

Используя 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()

Графический интерфейс пользователя позволяет вращать с помощью щелчка и перетаскивания и увеличения / уменьшения с помощью щелчка и перетаскивания правой кнопкой мыши.

alt text

...