Where are the star trackers pointing?¶
It is possible to use PandoraSpacecraft to see where the star trackers have pointed during prior Pandora observations.
In [1]:
Copied!
import numpy as np
import matplotlib.pyplot as plt
import pandoraspacecraft as psc
from astropy.time import Time
import astropy.units as u
import numpy as np
import matplotlib.pyplot as plt
import pandoraspacecraft as psc
from astropy.time import Time
import astropy.units as u
In [2]:
Copied!
ps = psc.PandoraSpacecraft()
ps = psc.PandoraSpacecraft()
In [3]:
Copied!
t = Time(np.linspace(2461135, 2461136, 10000), format='jd')
t = Time(np.linspace(2461135, 2461136, 10000), format='jd')
In [4]:
Copied!
st1_to_earth = ps.get_angle_to_body(t, direction='st1', body='earth')
st2_to_earth = ps.get_angle_to_body(t, direction='st2', body='earth')
st1_to_earth = ps.get_angle_to_body(t, direction='st1', body='earth')
st2_to_earth = ps.get_angle_to_body(t, direction='st2', body='earth')
In [5]:
Copied!
fig, ax = plt.subplots(facecolor='white', dpi=150)
ax.plot(t.jd, st1_to_earth.value, label='ST1', c='r')
ax.plot(t.jd, st2_to_earth.value, label='ST2', c='b')
ax.set(xlabel='Time [JD]', ylabel='Angle to Earth [deg]')
ax.legend();
fig, ax = plt.subplots(facecolor='white', dpi=150)
ax.plot(t.jd, st1_to_earth.value, label='ST1', c='r')
ax.plot(t.jd, st2_to_earth.value, label='ST2', c='b')
ax.set(xlabel='Time [JD]', ylabel='Angle to Earth [deg]')
ax.legend();
In [6]:
Copied!
st1_to_sun = ps.get_angle_to_body(t, direction='st1', body='sun')
st2_to_sun = ps.get_angle_to_body(t, direction='st2', body='sun')
st1_to_sun = ps.get_angle_to_body(t, direction='st1', body='sun')
st2_to_sun = ps.get_angle_to_body(t, direction='st2', body='sun')
In [7]:
Copied!
fig, ax = plt.subplots(facecolor='white', dpi=150)
ax.plot(t.jd, st1_to_sun.value, label='ST1', c='r')
ax.plot(t.jd, st2_to_sun.value, label='ST2', c='b')
ax.set(xlabel='Time [JD]', ylabel='Angle to Sun [deg]')
ax.legend();
fig, ax = plt.subplots(facecolor='white', dpi=150)
ax.plot(t.jd, st1_to_sun.value, label='ST1', c='r')
ax.plot(t.jd, st2_to_sun.value, label='ST2', c='b')
ax.set(xlabel='Time [JD]', ylabel='Angle to Sun [deg]')
ax.legend();