30 lines
878 B
Python
30 lines
878 B
Python
import numpy as np
|
|
import plotly.graph_objects as go
|
|
from ellipsoide import EllipsoidTriaxial
|
|
import winkelumrechnungen as wu
|
|
from dashboard import ellipsoid_figure
|
|
|
|
u = np.linspace(0, 2*np.pi, 51)
|
|
v = np.linspace(0, np.pi, 51)
|
|
ell = EllipsoidTriaxial.init_name("BursaSima1980round")
|
|
points = []
|
|
lines = []
|
|
for u_i, u_value in enumerate(u):
|
|
for v_i, v_value in enumerate(v):
|
|
cart = ell.ell2cart(u_value, v_value)
|
|
if u_i != 0 and v_i != 0:
|
|
lines.append((points[-1], cart, "red"))
|
|
points.append(cart)
|
|
points = []
|
|
for v_i, v_value in enumerate(v):
|
|
for u_i, u_value in enumerate(u):
|
|
cart = ell.ell2cart(u_value, v_value)
|
|
if u_i != 0 and v_i != 0:
|
|
lines.append((points[-1], cart, "blue"))
|
|
points.append(cart)
|
|
ax = ell.ax
|
|
ay = ell.ay
|
|
b = ell.b
|
|
|
|
figu = ellipsoid_figure(ax, ay, b, lines=lines)
|
|
figu.show() |