この質問は、おそらく次の質問に関連しています: Pyephem を使用した範囲レートの誤り
無線リンクの衛星通過予測とドップラー シフト計算を行う目的で、pyephem の実験を始めたところです。私はこれに非常に慣れていないので、まだコードを書いていません.Pythonインタラクティブインタープリターに手を出しただけです。
私の質問は、観測者の高度を 0mb に設定することで無効にできる大気屈折の補正についてです。pyephem で next_pass() を使用したパス予測を gnu predict からのパス予測と比較すると、AOS 時間と LOS 時間の間に 30 秒のオーダーの違いがあることに気付きました。大気屈折の補正が無効になりました。これは、predict と pyephem で使用される同一の TLE によるものです。
Q1) 最初に 30 秒程度の大きな差が生じたのは、通過予測で大気の屈折が補正されていないことが予測されたためではないでしょうか? これはもっともらしい説明のように思えます。
Q2) ほとんどの地上電波伝搬モデルは、地球半径の 4/3 のファッジ係数を使用して大気の屈折を補正し、曲線経路を効果的にまっすぐにして、計算された範囲に影響を与え、したがって範囲レートに影響を与えることを、私の (おそらく弱い) 理解です。pyephem の大気屈折補正は、このアプローチと同等ですか?
Q3) オブザーバーの圧力を 1010.0mb から 0.0mb に変更しても (したがって、大気屈折補正を無効にします)、range と range_rate の値は影響を受けていないように見えますが、next_pass() を呼び出した結果は影響を受けます。これは事実ですか？range と range_rate の値が変わると思っていたのですが、実際には、「見かけの距離」と「見かけの距離率」を、大気の屈折補正が適用されたときに計算できる追加の値と見なすべきでしょうか?
これは、スタック オーバーフローに関する私の最初の質問です。クエリが適切に形成されていない場合は、お詫び申し上げます。ありがとう。
更新 - 要求に応じた一部のコードと一部の結果
ここにいくつかのコードがあります...
import ephem
g0hww = ephem.Observer()
g0hww.lat='52:14:15.70'
g0hww.lon='0:43:24.49'
g0hww.elevation=40.0
g0hww.pressure = 0.0
g0hww.date = "2014/5/9 08:48:53"
stations = open("/home/darren/cronjobs/stations.txt","r")
line1=stations.readline()
line2=stations.readline()
line3=stations.readline()
print "Using keps for:" + line1.strip()
print line2.strip()
print line3.strip()
print "\n"
iss = ephem.readtle(line1,line2,line3)
iss.compute(g0hww)
print g0hww
print "iss azimuth: " + str(iss.az)
print "iss elevation: " + str(iss.alt)
print "iss range: " + str(iss.range)
print "iss range rate: " + str(iss.range_velocity)
print g0hww.next_pass(iss)
print "\n"
g0hww.pressure = 1010.0
iss.compute(g0hww)
print g0hww
print "iss azimuth: " + str(iss.az)
print "iss elevation: " + str(iss.alt)
print "iss range: " + str(iss.range)
print "iss range rate: " + str(iss.range_velocity)
print g0hww.next_pass(iss)
ここにいくつかの結果があります：
Using keps for:ISS (ZARYA)
1 25544U 98067A 14126.92299264 .00007994 00000-0 14853-3 0 2319
2 25544 51.6500 304.1061 0002752 337.8472 139.1065 15.49877967884923
<ephem.Observer date='2014/5/9 08:48:52' epoch='2000/1/1 12:00:00' lon=0:43:24.5 lat=52:14:15.7 elevation=40.0m horizon=0:00:00.0 temp=15.0C pressure=0.0mBar>
iss azimuth: 90:47:19.1
iss elevation: 21:02:43.0
iss range: 999640.75
iss range rate: 5735.83300781
(2014/5/9 10:18:09, 269:24:46.2, 2014/5/9 08:48:44, 22:40:03.7, 2014/5/9 08:52:14, 79:47:08.3)
<ephem.Observer date='2014/5/9 08:48:52' epoch='2000/1/1 12:00:00' lon=0:43:24.5 lat=52:14:15.7 elevation=40.0m horizon=0:00:00.0 temp=15.0C pressure=1010.0mBar>
iss azimuth: 90:47:19.1
iss elevation: 21:05:11.0
iss range: 999640.75
iss range rate: 5735.83300781
(2014/5/9 10:18:04, 269:25:01.5, 2014/5/9 08:48:49, 21:49:15.9, 2014/5/9 08:52:29, 79:31:12.8)