37 lines
1.0 KiB
Python
37 lines
1.0 KiB
Python
|
#!/usr/bin/env python
|
||
|
|
from __future__ import print_function
|
||
|
|
import numpy as np
|
||
|
|
from astropy.io import ascii
|
||
|
|
import matplotlib.pyplot as plt
|
||
|
|
import astropy.units as u
|
||
|
|
from astropy.coordinates import SkyCoord,FK5,AltAz,EarthLocation
|
||
|
|
from astropy.time import Time
|
||
|
|
|
||
|
|
table=ascii.read("imgstat.csv",format="csv")
|
||
|
|
|
||
|
|
t=Time(table['mjd'],format="mjd",scale="utc")
|
||
|
|
pos=SkyCoord(ra=table['ra'],dec=table['de'],frame="icrs",unit="deg")
|
||
|
|
|
||
|
|
loc=EarthLocation(lat=52.8344*u.deg,lon=6.3785*u.deg,height=10*u.m)
|
||
|
|
|
||
|
|
pa=pos.transform_to(AltAz(obstime=t,location=loc))
|
||
|
|
|
||
|
|
mjd0=np.floor(np.min(table['mjd']))
|
||
|
|
|
||
|
|
plt.subplot(411)
|
||
|
|
|
||
|
|
plt.plot(table['mjd']-mjd0,table['mean'],label='Brightness')
|
||
|
|
plt.plot(table['mjd']-mjd0,table['std'],label='Variation')
|
||
|
|
plt.legend()
|
||
|
|
plt.subplot(412)
|
||
|
|
plt.plot(table['mjd']-mjd0,pa.az.degree,label='Azimuth')
|
||
|
|
plt.legend()
|
||
|
|
plt.subplot(413)
|
||
|
|
plt.plot(table['mjd']-mjd0,pa.alt.degree,label='Altitude')
|
||
|
|
plt.legend()
|
||
|
|
plt.subplot(414)
|
||
|
|
plt.plot(table['mjd']-mjd0,table['rmsx'],label='RA rms')
|
||
|
|
plt.plot(table['mjd']-mjd0,table['rmsy'],label='Dec rms')
|
||
|
|
plt.legend()
|
||
|
|
plt.show()
|