3

与えられた日時の地球の位置（太陽に対する）と軸の回転を計算したいと思います。太陽が0,0,0座標で静止していると仮定しても問題ありません。たとえば、月の引力によるごくわずかなたわみも無視できます。ある程度の範囲内で正確なものであれば十分です。

これを達成するのに役立つライブラリ/ソース/データはありますか？

4

3 に答える 3

4

ここからダウンロードできるaa-56コードには、おそらくあなたのニーズを満たす太陽暦が含まれています。高精度の作業には、JPLのDE421のようなより正確なものが必要ですが、不便なほど大きな係数のテーブルが含まれているため、1度の精度に満足している場合はおそらく極端にやり過ぎです。

与えられた時間での地球の自転は、グリニッジ恒星時によって与えられます。
Jean Meeusの「AstronomicalAlgorithms」（この種の計算に役立つ参考資料です！）はtheta0、ユリウス日に関する（累積回転角（度単位））の式を示していますJD。

T = (JD - 2451545.0 ) / 36525

theta0 = 280.46061837 + 360.98564736629*(JD-2451545.0) + 
           0.000387933*T*T - T*T*T/38710000.0

theta0 =0度mod360は、グリニッジ子午線が天球座標で赤経0:00に位置合わせされた瞬間を表します。

于 2010-12-02T00:08:19.047 に答える
2

はい。必要なのはエフェメリスパッケージです。

JPLには、計算を行うオンラインのエフェメリスサービスがあるため、Webに対応している場合は、それを利用できます。

無料のEphermerisパッケージも見つけましたが、JPLサイトにヒットしたようです。ただし、JPLのデータベースをダウンロードして、それを機能させるためのリンクがあります。したがって、オフラインで作業したい場合で、JPLからデータベースを手動で更新することを気にしない場合は、それがオプションになる可能性があります。

于 2010-12-01T23:55:35.653 に答える
2

Pythonでは、ephemライブラリを使用します。

>>> import ephem
>>> sun = ephem.Sun()
>>> sun.compute(ephem.now())
>>> sun.hlong, sun.hlat, sun.earth_distance
(69:41:32.6, 0:00:00.2, 0.98602390289306641)

ephem地球を体として便利に表現することはできませんが、地球sun.hlongのsun.hlat地動説の経度と緯度を示します。これはより適切に文書化される可能性があります。

地球の自転については、ここで探している値を言うことができます。（通常は時刻を使用しますが、それを把握する方法は一般的に明らかだと思います！）

于 2010-12-02T00:11:58.653 に答える