spacecruft
/
sattools
1
0
Fork 0
Code Releases Activity

Fix astropy issue

pull/13/merge
Cees Bassa 2018-03-04 10:16:38 +01:00
parent cbebb43231
commit 311db65f59
2 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
python/extract_tracks.py
View File

@ -158,9 +158,10 @@ if __name__ == '__main__':
# extract_tracks("2018-02-26T05:26:15.801.fits",trkrmin,drdtmin,trksig,ntrkmin)
files=sorted(glob.glob("2*.fits"))
files=sorted(glob.glob("2018-03-04T02*.fits"))
for file in files:
print(file)
extract_tracks(file,trkrmin,drdtmin,trksig,ntrkmin)

8
python/stio.py
View File

@ -33,8 +33,12 @@ class observation:
self.ymax=self.y0+self.dydt*(self.tmax-self.tmid)
# Correct for rotation
hobs=Time(ff.mjd+0.5*ff.texp/86400.0,format='mjd',scale='utc').sidereal_time("mean",longitude=0.0).degree
hmid=Time(self.mjd,format='mjd',scale='utc').sidereal_time("mean",longitude=0.0).degree
tobs=Time(ff.mjd+0.5*ff.texp/86400.0,format='mjd',scale='utc')
tobs.delta_ut1_utc=0
hobs=tobs.sidereal_time("mean",longitude=0.0).degree
tmid=Time(self.mjd,format='mjd',scale='utc')
tmid.delta_ut1_utc=0
hmid=tmid.sidereal_time("mean",longitude=0.0).degree
# Compute ra/dec
world=ff.w.wcs_pix2world(np.array([[self.x0,self.y0]]),1)