Untitled
Never
def makecubelimits(axis, centers=None, hw=None): lims = ax.get_xlim(), ax.get_ylim(), ax.get_zlim() if centers == None: centers = [0.5*sum(pair) for pair in lims] if hw == None: widths = [pair[1] - pair[0] for pair in lims] hw = 0.5*max(widths) ax.set_xlim(centers[0]-hw, centers[0]+hw) ax.set_ylim(centers[1]-hw, centers[1]+hw) ax.set_zlim(centers[2]-hw, centers[2]+hw) print("hw was None so set to:", hw) else: try: hwx, hwy, hwz = hw print("ok hw requested: ", hwx, hwy, hwz) ax.set_xlim(centers[0]-hwx, centers[0]+hwx) ax.set_ylim(centers[1]-hwy, centers[1]+hwy) ax.set_zlim(centers[2]-hwz, centers[2]+hwz) except: print("nope hw requested: ", hw) ax.set_xlim(centers[0]-hw, centers[0]+hw) ax.set_ylim(centers[1]-hw, centers[1]+hw) ax.set_zlim(centers[2]-hw, centers[2]+hw) return centers, hw TLE = """1 25544U 98067A 19161.11058730 .00000940 00000-0 23612-4 0 9994 2 25544 51.6452 36.4196 0007986 32.0937 125.3020 15.51186091174156""" L1, L2 = TLE.splitlines() from skyfield.api import Loader, EarthSatellite from skyfield.timelib import Time import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D halfpi, pi, twopi = [f*np.pi for f in [0.5, 1, 2]] degs, rads = 180/pi, pi/180 load = Loader('~/Documents/fishing/SkyData') data = load('de421.bsp') ts = load.timescale() planets = load('de421.bsp') earth = planets['earth'] Roadster = EarthSatellite(L1, L2) print(Roadster.epoch.tt) hours = np.arange(0, 3, 0.01) time = ts.utc(2018, 2, 7, hours) Rpos = Roadster.at(time).position.km Rposecl = Roadster.at(time).ecliptic_position().km print(Rpos.shape) re = 6378. theta = np.linspace(0, twopi, 201) print("Theta: ",theta) cth, sth, zth = [f(theta) for f in [np.cos, np.sin, np.zeros_like]] lon0 = re*np.vstack((cth, zth, sth)) lons = [] for phi in rads*np.arange(0, 180, 15): cph, sph = [f(phi) for f in [np.cos, np.sin]] lon = np.vstack((lon0[0]*cph - lon0[1]*sph, lon0[1]*cph + lon0[0]*sph, lon0[2]) ) lons.append(lon) lat0 = re*np.vstack((cth, sth, zth)) lats = [] for phi in rads*np.arange(-75, 90, 15): cph, sph = [f(phi) for f in [np.cos, np.sin]] lat = re*np.vstack((cth*cph, sth*cph, zth+sph)) lats.append(lat) if True: fig = plt.figure(figsize=[10, 8]) # [12, 10] ax = fig.add_subplot(1, 1, 1, projection='3d') x, y, z = Rpos ax.plot(x, y, z) for x, y, z in lons: ax.plot(x, y, z, '-k') for x, y, z in lats: ax.plot(x, y, z, '-k') centers, hw = makecubelimits(ax) print("centers are: ", centers) print("hw is: ", hw) plt.show() r_Roadster = np.sqrt((Rpos**2).sum(axis=0)) alt_roadster = r_Roadster - re if True: plt.figure() plt.plot(hours, r_Roadster) plt.plot(hours, alt_roadster) plt.xlabel('hours', fontsize=14) plt.ylabel('Geocenter radius or altitude (km)', fontsize=14) plt.show()
Raw Text
-
Lily and Redhead RyAnne play Strip Mogadishu
22 min ago
-
dsadsa
31 min ago
-
fg
34 min ago
-
ghjklñ
35 min ago
-
dsadsda
38 min ago
-
fdsfdsfds
40 min ago
-
Belal Hamideh Law - Abogados de Accidentes
47 min ago
-
fghjklñ
50 min ago
-
Stepsister Lost Her Pussy At Cards. Unclasp Your Pussy, We'll Fuck!
52 min ago
-
jkljhgf
1 hour ago