を始めたばかりですPyEphem
。私の当面の仕事は、日付と時刻が与えられたときに、地球上の亜太陽点をlatitude-longitude
値で計算することです。これを解決するために掘り下げますPyEphem
が、誰かがすでにこれを行っている場合は、サンプルコードをいただければ幸いです。
2 に答える
今朝の場所から実際のコードをテストすることはできませんが、赤緯ϕのオブジェクトは、緯度ϕを持つ地球上の一連の場所の真上に常に乗る必要があるため、緯度番号は身体の属性によって直接与えられ.dec
ます (または.a_dec
または.g_dec
アプリケーションによって異なります)。
では、経度はどうでしょうか。
状況を想像してみてください。ほぼ 1 日に 1 回、経度 0°のグリニッジが見上げると、真上に赤経 0°の線が見えます。その瞬間、PyEphem の場合と同様に、赤経θ の空にある天体は経度 θを見下ろすことになります。
では、グリニッジがゼロ以外の赤経線を見上げているとしたらどうでしょうか? それから、経度を作るために天体の赤経からそれを差し引く必要があるように私には思えます。なぜなら、日が進み、地球が回転し、赤経線がグリニッジ上を通過するにつれて、赤経線はますます大きくなり、それらに割り当てられます。任意の天体は地球を横切って西に通過し、その経度は減少し、西半球を通過するときにマイナスになります。
任意の時点でのグリニッジの頭上の赤経線は、経度 0° でオブザーバーを作成し.sidereal_time()
、クイック リファレンスを正しく思い出せるかどうかを尋ねることで決定できます。したがって、体の下の経度は次のようになると思います。
lon = body.ra - greenwich.sidereral_time()
合理的な数値が出てくるかどうかを確認するために、今日の後半にこれを使って簡単なテストを行います.