Моделирование астероида с помощью Matplotlib с использованием поверхности и каркаса

Question

Я пытаюсь смоделировать астероид, используя plot_surface а также plot_wireframe, У меня есть значения x y и z для точек на поверхности астероида. Каркас точно соответствует форме астероида, но график поверхности не соответствует каркасу. Как я могу заставить график поверхности соответствовать каркасу или как я могу использовать каркас, чтобы получить трехмерную твердотельную модель? Вот мой код для модели:

from mpl_toolkits.mplot3d import axes3d
import matplotlib.pyplot as plt
import numpy as np
from matplotlib import cm
from matplotlib.mlab import griddata

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
data = np.genfromtxt('data.txt')
x = data[:,0]
y = data[:,1]
z = data[:,2]

ax.plot_wireframe(x, y, z, rstride=1, cstride=1, alpha=1)

xi = np.linspace(min(x), max(x))
yi = np.linspace(min(y), max(y))

X, Y = np.meshgrid(xi, yi)
Z = griddata(x, y, z, xi, yi)


surf = ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=cm.coolwarm,
        linewidth=0, antialiased=False)
ax.set_zlim(-1.01, 1.01)

plt.show()

Данные в этом формате, хотя в исходном файле намного больше строк:

-1.7738946051191869E-002   4.3461451610545973E-002   1.3393057231408241     
-0.29733561550902488       0.32305812106837900        1.3393057231408241     
-0.29733561550902488       0.16510132228266330        1.3548631099230350     
-0.21872587865015569        2.4170900455101410E-002   1.3610011616437809     
1.4452975249810950E-002 -0.20900795344486520        1.3610011616437809    
1.5732454381265970E-002 -0.20900795344486520        1.3608751439485580    
-0.34501536374240321       0.51320241386595655        1.3158820995876130     
-0.40193014435941982       0.45628763324893978        1.3158820995876130     
-0.42505849480150409       0.28183419537116011        1.3307863198123011     
-0.18994178462386799      -0.19294290416565860        1.3424523041534830     
1.4452975249810939E-002 -0.39733766403933751        1.3424523041534830     
5.8021940902131752E-002 -0.57108837516584876        1.3210481842104100     
9.3746267961881152E-002 -0.61017602710257668        1.3136798474111200     
0.26609469681891229      -0.43782759824554562        1.3136798474111200     
  0.17938460413447810       0.39179924148155021        1.2357401964919650     
   8.9613011902522258E-002  0.42818009222325598        1.2584008460875080     
  0.33671539027096409      -0.47165177581327772        1.2965073126705291     
  0.53703772594296528      -0.47165177581327777        1.2357401964919561     
 -0.19242375014122229       0.71021685426700043        1.2584008460875080     
 -0.34501536374240321       0.66763766324752027        1.2904902860951690 

Надеюсь, вы можете помочь

Answer 1

Не могли бы вы предоставить представление о том, что происходит? Я предполагаю, что вы получаете изображение, где поверхность астероида, кажется, прыгает повсюду. Это верно? Если это так, это может быть вызвано тем, что плоттер не знает порядок точек.

Если у нас есть набор точек, который содержит все точки единичного круга, то мы ожидаем, что рисование этих точек создаст единичный круг. Однако, если вы решили соединить каждую точку с двумя другими точками, это не обязательно будет выглядеть как круг. Если (по какой-то причине) вы подключили одну точку к другой точке на другой стороне круга и продолжали делать это до тех пор, пока каждая точка не была соединена с двумя другими точками, она может или не может выглядеть как окружность, потому что каждая точка не обязательно подключается к соседним точкам.

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

Следующий пример круга должен проиллюстрировать мою точку зрения:

import math
import matplotlib.pylab as plt
import random

thetaList = range(360)
random.shuffle(thetaList)
degToRad = lambda x: float(x) * math.pi / float(180)
x = [math.cos(degToRad(theta)) for theta in thetaList]
y = [math.sin(degToRad(theta)) for theta in thetaList]       

#plot the cirlce
plt.plot(x,y)
plt.show()