私は以前に pyephem を使用したことがなく、衛星測位の専門家でもありません。pyephem を利用して、TLE を使用して衛星の位置を計算したいと考えています。次のように、非常に簡単なことをしなければなりません。
tle=["ISS (ZARYA)","1 25544U 98067A 03097.78853147 .00021906 00000-0 28403-3 0 8652","2 25544 51.6361 13.7980 0004256 35.6671 59.2566 15.58778559250029"]
iss = ephem.readtle(*tle)
observer = ephem.Observer()
observer.lon, observer.lat = ('-84.39733', '33.775867')
observer.date = ephem.Date('2002/4/23 10:10:00.000')
iss.compute(observer)
print iss.alt, iss.az, iss.range
-40:06:46.3 199:08:24.3 8834968.0
これらの 3 つの変数は、地平線基準系における衛星の位置を提供します。pyephem がこの値をどのように計算するかは明確ではありません。リファレンス ガイドを読みました: http://rhodesmill.org/pyephem/radec
ドキュメントを読むと、pyephem は歳差運動と章動を適用しているように見えますが、ドキュメントの最後の 2 行には次のように書かれています。
"最後の 2 つの座標セットのいずれにも歳差運動が適用されていないことに注意してください。最初の座標にのみ適用されます。これは、「アストロメトリック」位置のみが星図の線に対応することを意味します。他の位置は「 epoch-of-date」座標であり、観測自体のまさにその日の天の極と天の赤道の方向から測定されます。
az と alt に地球の歳差運動は適用されますか?
さらに、歳差運動と章動運動に pyephem が使用するモデルの種類を知りたいです (本当に参照が必要です)。Xephem と libastro へのリンクがありますが、アルゴリズムについては何も見つかりません。何か提案はありますか?
どうもありがとうございました!