построение точечного графика в python с использованием matplotlib

Question

У меня есть файл с 3000 точек в нем.Я должен построить трехмерный график рассеяния с этими точками. Точки таковы:

6.655597594660206395e-01,-5.205175889492101859e-01,4.583497554501108073e-01
3.357418723194116605e-02,-2.482341476533565849e+00,2.009030294705723030e+00
-1.398411818716352728e-01,-1.990250936356241063e+00,2.325394845551588929e+00

есть 3000 таких точек.

У меня есть код, и он показывает ошибку

TypeError: Cannot cast array data from dtype('float64') to dtype('<U32') according to the rule 'safe'

код, который я написал

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

x=[]
y=[]
z=[]

with open('3dpd.out') as f:
    lines = f.readlines()
    for line in lines :
        x.append(line.split(',')[0])
        y.append(line.split(',')[1])
        z.append(line.split(',')[2])


fig = plt.figure()
ax = fig.add_subplot(111,projection='3d')

ax.set_title("Scatter Plot")
ax.set_xlabel('X Label')
ax.set_ylabel('Y Label')
ax.set_zlabel('Z Label')


ax.scatter(x, y, z, c='r')

plt.show()

Answer 1

Когда строки читаются readlines, они считываются как строки, поэтому вам необходимо преобразовать эти строки в числа.

Также обратите внимание, что я разделил каждый шаг в отдельной строке кода (например, split, float и т. Д.), Что обычно помогает при отладке.(И вам нужно только вызвать split один раз для каждой строки.)

Вот пример, который работает:

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

x=[]
y=[]
z=[]

with open('3dpd.out') as f:
    lines = f.readlines()
    for line in lines :
        xs, ys, zs =  line.split(',')
        xv, yv, zv = [float(v) for v in (xs, ys, zs)]
        x.append(xv)
        y.append(yv)
        z.append(zv)

print type(xs)
print type(xv)

fig = plt.figure()
ax = fig.add_subplot(111,projection='3d')

ax.set_title("Scatter Plot")
ax.set_xlabel('X Label')
ax.set_ylabel('Y Label')
ax.set_zlabel('Z Label')


ax.scatter(x, y, z, c='r')

plt.show()