PyEphem中的地球衛星的矩形座標

Question

我已經通過大量的PyEphem網頁，下載的代碼和調查對象進行了搜索，但是我還沒有找到一種方法來以某種名義獲得地球衛星的地球中心直角座標慣性座標系。我錯過了明顯的東西嗎？我正在尋找'傳統'x，y，z和x點，y點，z點。非常感謝您的任何建議。

Answer 1

好問題！

的libastro內的衛星位置程序不使用xyz內部協調，但libastro拋出他們離開的Python得到一個機會，看看在他們面前。我可以嘗試修補libastro然後PyEphem使數據可用，但在研究你的問題，我發現最近宏偉的工作已經完成更新SGP4衛星跟蹤算法，併爲其提供一個測試套件：

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

受此發現的啓發，我做了一些更好的事情：我花了整個週末的時間創建了一個基於此參考實現的新的純Python sgp4包！

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

不久，我將放棄這一新的封裝成PyEphem本身作爲其新的衛星預測引擎，但目前你應該能夠直接從Python包指數（以上鍊接）和呼叫下載它原始爲了得到你的傳統座標：

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) 

請試試看，讓我知道它是否可以在你的機器上運行！