Plotten Ellipsoid mit Matplotlib

Hat jemand Beispiel-code für das Plotten Ellipsoide? Es ist eine für die sphere auf matplotlib Ort, aber nichts für Ellipsoide. Ich bin versucht zu zeichnen

x**2 + 2*y**2 + 2*z**2 = c

wo c ist eine Konstante (wie 10), der definiert, ein ellipsoid. Ich habe versucht, die meshgrid(x,y) route, überarbeitet die Gleichung so z ist auf der einen Seite, aber die sqrt ist ein problem. Die matplotlib Bereich z.B. arbeiten mit Winkeln, u,v, aber ich bin nicht sicher, wie Sie Sie zu arbeiten, die für das ellipsoid.

Die Referenz, die Sie geben, beschreibt die Darstellung eines ellipse, während die Frage über eine Ellipse, welches die drei-dimensionale Entsprechung einer ellipse.
EOL: Du hast Recht, werde ich entfernen meinen Kommentar.

InformationsquelleAutor user423805 | 2011-10-19

  1. 18

    Hier ist, wie Sie es tun können, über sphärische Koordinaten:

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

fig = plt.figure(figsize=plt.figaspect(1))  # Square figure
ax = fig.add_subplot(111, projection='3d')

coefs = (1, 2, 2)  # Coefficients in a0/c x**2 + a1/c y**2 + a2/c z**2 = 1 
# Radii corresponding to the coefficients:
rx, ry, rz = 1/np.sqrt(coefs)

# Set of all spherical angles:
u = np.linspace(0, 2 * np.pi, 100)
v = np.linspace(0, np.pi, 100)

# Cartesian coordinates that correspond to the spherical angles:
# (this is the equation of an ellipsoid):
x = rx * np.outer(np.cos(u), np.sin(v))
y = ry * np.outer(np.sin(u), np.sin(v))
z = rz * np.outer(np.ones_like(u), np.cos(v))

# Plot:
ax.plot_surface(x, y, z,  rstride=4, cstride=4, color='b')

# Adjustment of the axes, so that they all have the same span:
max_radius = max(rx, ry, rz)
for axis in 'xyz':
    getattr(ax, 'set_{}lim'.format(axis))((-max_radius, max_radius))

plt.show()

    Das resultierende Diagramm ist ähnlich

    Plotten Ellipsoid mit Matplotlib

    Dem oben genannten Programm tatsächlich produziert ein schöner Blick "square" Grafik.

    Diese Lösung ist stark inspiriert von der Beispiel in Matplotlib-Galerie.

    vielen Dank für diese Antwort. Die Achse Anpassung Zeilen gab mir eine Fehlermeldung, jedoch "AttributeError: 'Axes3DSubplot' - Objekt hat kein Attribut 'set_zlim'". Ist es, weil ich bin auf eine ältere version von matplotlib (ich bin auf 1.0.1)
    yep, das war das problem. Ich habe ein Upgrade auf 1.1.0 und es funktioniert alles jetzt.
    rx, ry, rz = 1/np.sqrt(coef) werden Sollte, rx, ry, rz = 1/np.sqrt(coefs)
    Wenn ich semi-major axis a und c, semi-Moll-Achse b, wird die coeff (1/a2,1/b2,1/c**2) ?
    Ja. Natürlich in diesem Fall können Sie direkt zu tun rx, ry, rz = a, b, c (keine Notwendigkeit für die Koeffizienten coefs).

    InformationsquelleAutor Eric Lebigot

  2. 9

    Gebäude auf EOL Antwort. Manchmal muss man ein ellipsoid in matrix format:

    A und c, Wobei A die ellipsoid-matrix und c ein Vektor, repräsentiert das Zentrum des Ellipsoids.

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

# your ellispsoid and center in matrix form
A = np.array([[1,0,0],[0,2,0],[0,0,2]])
center = [0,0,0]

# find the rotation matrix and radii of the axes
U, s, rotation = linalg.svd(A)
radii = 1.0/np.sqrt(s)

# now carry on with EOL's answer
u = np.linspace(0.0, 2.0 * np.pi, 100)
v = np.linspace(0.0, np.pi, 100)
x = radii[0] * np.outer(np.cos(u), np.sin(v))
y = radii[1] * np.outer(np.sin(u), np.sin(v))
z = radii[2] * np.outer(np.ones_like(u), np.cos(v))
for i in range(len(x)):
    for j in range(len(x)):
        [x[i,j],y[i,j],z[i,j]] = np.dot([x[i,j],y[i,j],z[i,j]], rotation) + center

# plot
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_wireframe(x, y, z,  rstride=4, cstride=4, color='b', alpha=0.2)
plt.show()
plt.close(fig)
del fig

    So, nicht allzu viel neues dabei, aber hilfreich, wenn Sie haben ein ellipsoid in matrix-form, die gedreht wird und vielleicht nicht zentriert auf die Koordinaten 0,0,0 und zeichnen möchten.

    Hi, warum Sie benötigen, um zu berechnen, die inverse Quadratwurzel der singulären Werte, statt nur die Quadratwurzel, um die Radien? Ich meine, warum Radien = 1.0/np.sqrt(s) anstelle von Radien = np.sqrt(s). Ich denke, der später gibt die richtigen Radien.

    InformationsquelleAutor minillinim

Schreibe einen Kommentar