python绘制三维图像球_【python】读取三维点云球坐标数据并动态生成三维图像与着色...

import numpy as np

import math

import matplotlib.pyplot as plt

from mpl_toolkits.mplot3d import Axes3D

from matplotlib.animation import FuncAnimation

f=open('data.txt','r')

point=f.read()

f.close()

l1=point.replace('\n',',')

l2=l1.replace(' ',',')

l3=l2.split(',')

#print(l3)

m1=np.array(l3)

m2=m1.reshape(2664,6)

m3=[]

for each in m2:

each_line=list(map(lambda x:float(x),each))

m3.append(each_line)

m4=np.array(m3)

print(m4)

yaw=[i[0] for i in m4]

pitch=[i[1] for i in m4]

radius=[i[2] for i in m4]

c1=[i[3:6] for i in m4]

c2=np.array(c1)

c3=c2.reshape(2664,3)

print(c3)

x=[]

y=[]

z=[]

for i in range(len(m4)):

x.append(radius[i]*math.sin(0.0174532924*pitch[i])*math.sin(0.0174532924*yaw[i]))

y.append(radius[i]*math.cos(0.0174532924*pitch[i]))

z.append(radius[i]*math.sin(0.0174532924*pitch[i])*math.cos(0.0174532924*yaw[i]))

# print(x)

# print(y)

# print(z)

def animate():

return point

def init():

return point

fig=plt.figure(figsize=(16,9),dpi=120)

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

plt.title('point')

ax.set_xlabel('X Label')

ax.set_ylabel('Z Label')

ax.set_zlabel('Y Label')

anim = FuncAnimation(fig,animate, frames=np.arange(100), init_func=init,

interval=100, blit=True)

for i in range(2664):

point=ax.scatter(z[i],x[i],y[i],c=(c3[i]/255),marker='.',s=10,linewidth=1,alpha=1,cmap='spectral')

plt.ion()

plt.pause(0.01)

plt.close

plt.show()

视频地址：点击打开链接

数据点未进行着色时的情况显示如下：