python ISS跟踪

Posted 2021-05-09

#!/usr/bin/python

import sys
import math
import ephem

# https://www.celestrak.com/NORAD/elements/stations.txt
iss = ephem.readtle("ISS (ZARYA)",
	"1 25544U 98067A   17219.16976944  .00001980  00000-0  36996-4 0  9991",
	"2 25544  51.6417 144.1699 0006441 110.1117 332.4530 15.54266961 69681")

obs = ephem.Observer()
obs.lat = '50.06143'
obs.long = '19.93658'

for p in range(3):
	tr, azr, tt, altt, ts, azs = obs.next_pass(iss)
	print("Date/Time (UTC)       Alt/Azim	  Lat/Long	Elev")
	print("=====================================================")
	while tr < ts:
		obs.date = tr
		iss.compute(obs)
		print("%s | %4.1f %5.1f | %4.1f %+6.1f | %5.1f" % \
			(tr, 
			 math.degrees(iss.alt), 
			 math.degrees(iss.az), 
			 math.degrees(iss.sublat), 
			 math.degrees(iss.sublong), 
			 iss.elevation/1000.))
		tr = ephem.Date(tr + 20.0 * ephem.second)
	print
	obs.date = tr + ephem.minute