3

试图计算的是卫星上升到观察者地平线以上的时间 (TO)、卫星达到地平线以上最大仰角的时间 (TCA) 以及卫星在一天内经过观察者地平线以下的时间 (LOS) (2013/06/01)，我一次路过发现了一个错误

import ephem
import math

line1 = "ISS (ZARYA)"
line2 = "1 25544U 98067A   13150.58028350  .00009895  00000-0  16910-3 0  9060"
line3 = "2 25544 051.6474 195.6448 0009976 014.6542 345.5254 15.50593990236071"
iss = ephem.readtle(line1, line2, line3)

campos = ephem.Observer()
campos.lon,campos.lat,campos.elevation = '-41:30.00','-21:75.00', 0 

campos.date = '2013/6/1'
fim = campos.date+1

print """Date/Time (UTC)       Alt/Azim   Lat/Long  Elev"""
print """====================================================="""
while  campos.date < fim:

   info = campos.next_pass(iss) 
   print info
   for p in range(0,6,2):

      campos.date = info[p]   
      iss.compute(campos)
      print "%s | %4.1f %5.1f | %4.1f %+6.1f | %5.1f" % \
    (info[p], math.degrees(iss.alt), math.degrees(iss.az), 
     math.degrees(iss.sublat), math.degrees(iss.sublong),iss.elevation/1000.)
   print " "

输出：

日期/时间 (UTC) Alt/Azim Lat/Long Elev ===================================== ================ (2013/6/1 10:01:51, 31:57:07.9, 2013/6/1 10:04:21, 1:57:44.6 , 2013/6/1 10:06:52, 86:00:15.8) 2013/6/1 10:01:51 | 0.0 32.0 | -4.2 -30.7 | 422.7 2013/6/1 10:04:21 | 2.0 59.1 | -11.8 -25.1 | 423.7 2013/6/1 10:06:52 | -0.0 86.0 | -19.2 -19.3 | 424.5

(2013/6/1 11:34:59, 323:43:50.2, 2013/6/1 11:40:35, 76:09:23.5, 2013/6/1 11:46:11, 138:14:56.1)
 2013/6/1 11:34:59 |  0.2 323.7 | -5.2  -53.5 | 422.8
 2013/6/1 11:40:35 | 76.2  68.9 | -21.8  -40.6 | 424.7
 2013/6/1 11:46:11 |  0.0 138.2 | -36.9  -24.1 | 424.9

 (2013/6/1 13:13:25, 266:18:49.6, 2013/6/1 11:43:23, 13:48:06.0, 2013/6/1 11:46:14, 138:15:55.2)
 2013/6/1 13:13:25 |  0.0 266.3 | -21.9  -64.1 | 424.7
 2013/6/1 11:43:23 | 13.8 136.3 | -29.7  -33.0 | 425.0
 2013/6/1 11:46:14 | -0.2 138.3 | -37.0  -23.9 | 424.9

 (2013/6/1 13:13:25, 266:18:50.8, 2013/6/1 13:17:32, 6:57:09.7, 2013/6/1 13:21:40, 171:25:48.7)
 2013/6/1 13:13:25 |  0.0 266.3 | -21.9  -64.1 | 424.7
 2013/6/1 13:17:32 |  7.0 218.6 | -33.2  -52.4 | 425.0
 2013/6/1 13:21:40 |  0.0 171.4 | -42.9  -37.3 | 424.5

 (2013/6/1 19:48:36, 195:26:26.5, 2013/6/1 12:42:57, -62:57:07.1, 2013/6/1 13:21:39, 171:33:53.8)
 2013/6/1 19:48:36 |  0.0 195.4 | -42.2  -48.9 | 418.2
 2013/6/1 12:42:57 | -63.0 318.8 | 51.6 -164.5 | 405.3
 2013/6/1 13:21:39 |  0.0 171.6 | -42.9  -37.4 | 424.5

(2013/6/1 19:48:36, 195:27:10.8, 2013/6/1 12:42:57, -62:57:11.0, 2013/6/1 13:21:39, 171:33:31.1)
2013/6/1 19:48:36 |  0.0 195.5 | -42.2  -48.9 | 418.2
2013/6/1 12:42:57 | -63.0 318.8 | 51.6 -164.5 | 405.3
2013/6/1 13:21:39 |  0.0 171.6 | -42.9  -37.4 | 424.5

(2013/6/1 19:48:36, 195:27:08.8, 2013/6/1 12:42:57, -62:57:10.9, 2013/6/1 13:21:39, 171:33:32.1)
2013/6/1 19:48:36 |  0.0 195.5 | -42.2  -48.9 | 418.2
2013/6/1 12:42:57 | -63.0 318.8 | 51.6 -164.5 | 405.3
2013/6/1 13:21:39 |  0.0 171.6 | -42.9  -37.4 | 424.5

(2013/6/1 19:48:36, 195:27:08.9, 2013/6/1 12:42:57, -62:57:10.9, 2013/6/1 13:21:39,  171:33:32.0)
2013/6/1 19:48:36 |  0.0 195.5 | -42.2  -48.9 | 418.2
2013/6/1 12:42:57 | -63.0 318.8 | 51.6 -164.5 | 405.3
2013/6/1 13:21:39 |  0.0 171.6 | -42.9  -37.4 | 424.5

(2013/6/1 19:48:36, 195:27:08.9, 2013/6/1 12:42:57, -62:57:10.9, 2013/6/1 13:21:39, 171:33:32.0)
2013/6/1 19:48:36 |  0.0 195.5 | -42.2  -48.9 | 418.2
2013/6/1 12:42:57 | -63.0 318.8 | 51.6 -164.5 | 405.3
2013/6/1 13:21:39 |  0.0 171.6 | -42.9  -37.4 | 424.5

......

我哪里错了？

4

0 回答 0