Skip to content

Instantly share code, notes, and snippets.

View pije76's full-sized avatar

panjianom pije76

View GitHub Profile
Traceback (most recent call last):
File "C:/Users/taiko/PycharmProjects/untitled/solar_system_toy.py", line 125, in <module>
orbit_list[idx] = orbit_list[idx].propagate(1*u.day,method=cowell,ad=ad)
File "C:\Users\taiko\AppData\Local\Programs\Python\Python37-32\lib\site-packages\poliastro\twobody\orbit.py", line 999, in propagate
cartesian = propagate(self, time_of_flight, method=method, rtol=rtol, **kwargs)
File "C:\Users\taiko\AppData\Local\Programs\Python\Python37-32\lib\site-packages\poliastro\twobody\propagation.py", line 468, in propagate
orbit.attractor.k, orbit.r, orbit.v, time_of_flight.to(u.s), rtol=rtol, **kwargs
File "C:\Users\taiko\AppData\Local\Programs\Python\Python37-32\lib\site-packages\poliastro\twobody\propagation.py", line 105, in cowell
dense_output=True,
File "C:\Users\taiko\AppData\Local\Programs\Python\Python37-32\lib\site-packages\scipy\integrate\_ivp\ivp.py", line 477, in solve_ivp
@pije76
pije76 / observe_skyfield.py
Last active July 14, 2021 20:14 — forked from jberlanga/observe_skyfield.py
Observing excercise with Skyfield
#!/usr/bin/env python3.4
"""
observe_skyfield.py
Observing exercise with Skyfield.
"""
import numpy as np
@pije76
pije76 / iss_transit_skyfield.py
Created July 14, 2021 20:13 — forked from cbassa/iss_transit_skyfield.py
Reproducing the ISS->Venus->Solar transit of June 8, 2004 (Skyfield)
#!/usr/bin/env python3
import numpy as np
from skyfield.api import Topos, load
from skyfield.api import EarthSatellite
import matplotlib.pyplot as plt
from matplotlib.patches import Circle
import astropy.units as u
from astropy.wcs import wcs
if __name__ == "__main__":
@pije76
pije76 / pdiddy_island_mars_selfie.py
Created July 14, 2021 20:12 — forked from nathanielatom/pdiddy_island_mars_selfie.py
Find out where Mars is, in the sky at P'diddy Island, PEI, ON, Canada.
#!/usr/bin/env python
"""
Find out where Mars is, in the sky at P'diddy Island, PEI, ON, Canada.
"""
import datetime
# ! pip install skyfield
from skyfield.api import load, Topos
from skyfield.api import Topos, load
import socket
from time import sleep
import re
import array as arr
#connecting to server
s = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
s.connect(('where.satellitesabove.me', 5021))
print("connecting to where.satellitesabove.me:5021")
from tkinter import *
from skyfield.api import Topos, load
import motor
ts = load.timescale()
#Motor al laser
m1 = motor.Motor([18, 22, 24, 26], mode = 3)
m1.rpm = 10
#Motor a la base
m2 = motor.Motor([17,25,12,16], mode = 3)
@pije76
pije76 / iss.py
Created July 14, 2021 20:11 — forked from rbs-tim/iss.py
iss.py
# python 3.7 suggested
import datetime
import time
from skyfield.api import load, Topos, EarthSatellite
TLE_FILE = "https://celestrak.com/NORAD/elements/active.txt"
ISS_NAME = "ISS (ZARYA)"
LONGITUDE = -73.9940702
LATITUDE = 40.7546283
from skyfield.api import load
from jdutil import *
planets = load('de432s.bsp')
earth, jupiter = planets['EARTH'], planets['JUPITER BARYCENTER']
ts = load.timescale()
t = ts.now()
timenow = jd_to_datetime(float(t.tt))
astrometric = earth.at(t).observe(jupiter)
ra, dec, distance = astrometric.radec()
#!/usr/bin/env python3
from datetime import datetime, timezone
from math import atan, atan2, degrees, floor, pi, radians, sin, sqrt
from skyfield.api import earth, JulianDate, now, sun
def earth_latlon(x, y, z, time):
"""
@pije76
pije76 / moon.py
Created July 14, 2021 20:11 — forked from 910JQK/moon.py
A program to calculate exact time of the four special moon phases.
#!/usr/bin/env python3
import sys
from skyfield.api import load, utc
from datetime import datetime
TS_ZERO = datetime(1970, 1, 1)
NAMES = ['朔', '上弦', '望', '下弦']