Skip to content

Instantly share code, notes, and snippets.

View pije76's full-sized avatar

panjianom pije76

View GitHub Profile
@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 = ['朔', '上弦', '望', '下弦']
@pije76
pije76 / orbit_compare.py
Created July 14, 2021 13:40
Compare PyEphem and Skyfield modelled position to actual position
# Sentinel 1a position data taken from https://qc.sentinel1.eo.esa.int
# POD Precise Orbit Ephemerides AUX_POEORB (only last 366 days)
#
# S1A_OPER_AUX_POEORB_OPOD_20151028T122426_V20151007T225943_20151009T005943
#
# !! Coordinate system unknow, indicates fixed but not which one assume ITRF
#
# <OSV>
# <TAI>TAI=2015-10-08T21:51:09.000000</TAI>
# <UTC>UTC=2015-10-08T21:50:33.000000</UTC>