5

多数の PyEphem Web ページを検索し、コードをダウンロードし、オブジェクトを調査しましたが、公称慣性座標系で地球衛星の地球中心の直交座標を取得する方法が見つかりませんでした。明らかな何かを見逃しましたか?「従来の」x、y、z および x-dot、y-dot、z-dot を探しています。ご提案いただきありがとうございます。

4

1 に答える 1

5

良い質問!

内部の衛星位置ルーチンは内部で座標をlibastro使用しますが、Python がそれらを見る機会を得る前にそれらを破棄します。パッチを適用してから PyEphem を試行してデータを利用可能にすることもできましたが、質問を調査しているうちに、SGP4 衛星追跡アルゴリズムを更新し、テスト スイートを提供するという素晴らしい最近の作業が行われたことがわかりました。x y zlibastrolibastro

http://www.celestrak.com/publications/AIAA/2006-6753/

この発見に刺激されて、私ははるかに良いことをしました:sgp4このリファレンス実装に基づいて新しい純粋な Python パッケージを作成するために週末を費やしました!

http://pypi.python.org/pypi/sgp4

すぐに、この衛星予測エンジンを中心に Skyfield という名前の天文学ライブラリ全体を構築しますが、当面は Python Package Index (上記のリンク) から直接ダウンロードして、従来の座標を取得するために生で呼び出すことができるはずです。

from sgp4.earth_gravity import wgs72
from sgp4.io import twoline2rv

line1 = ('1 00005U 58002B   00179.78495062  '
         '.00000023  00000-0  28098-4 0  4753')
line2 = ('2 00005  34.2682 348.7242 1859667 '
         '331.7664  19.3264 10.82419157413667')

satellite = twoline2rv(line1, line2, wgs72)
position, velocity = satellite.propagate(
    2000, 6, 29, 12, 50, 19)

試してみて、あなたのマシンで動作するかどうか教えてください!

于 2012-08-26T19:28:27.250 に答える