問題タブ [pyephem]

これは python/pyephem の質問ではなく、軌道力学の質問ですが、次のようになります。

Pyephem は、春分が Thu 2014-03-20 12:57 EST に発生すると報告しています。これは、日の出が真東 (90°) に最も近くなり、日の入りが真東 (270°) に最も近くなる日であり、ほぼ世界中で行われています。

その日の経度 -78:30、緯度 35:30 標高 124m の時刻と方位角の日の出/設定を計算します。水平線は「-0:34」、気圧は 0 です (米国海軍天文台にできるだけ近い数値を取得するため)。私は得ます（現地時間、ESTに合わせて調整されます）：

太陽正午が分点の少し後に発生すると、日の出が 90:00:00.0 に近づくよりも日没が 270:00:00.0 に近づくと予想されますが、その逆です。日の出は真東から約 4.5 分、日の入りは真西から約 19 分です。なんで？

私は PyEphem を使用して、星座設計プロジェクトのある地点での衛星の通過を予測しています。

私はそれが軌道エポックであると私の質問に答えますが、私は確信したいと思います. 最初に疑問を投げかけた例を次に示します。

機能するもの

PyEphem (アクセスごとに 1 つの日付のみ) で取得した結果を、STK と比較して以下に示します (取得方法については、以下を参照してください)。

PyEphem のセットアップ

STK のセットアップは次のとおりです (実験の日付は同じです): STK 衛星のセットアップ

STK と PyEphem エポックの間になぜこの 1 時間の違いがあるのか​​ 理解できなくても、UTCG の何かを疑っています。

しないもの

ここで、傾きを 49 に上げると、次の結果が得られます (最初に PyEphem、次に Stk)。

エポックを に設定すると2014/02/14 11:35:00、期待どおりの結果が得られます。

最初の例は満足のいく結果を示しているので、私のモデルが違いを引き起こしたとは言えませんが、違いはありました。だから私の質問は、このエポック時間で何が間違っているのですか?

この質問は、おそらく次の質問に関連しています: Pyephem を使用した範囲レートの誤り

無線リンクの衛星通過予測とドップラー シフト計算を行う目的で、pyephem の実験を始めたところです。私はこれに非常に慣れていないので、まだコードを書いていません.Pythonインタラクティブインタープリターに手を出しただけです。

私の質問は、観測者の高度を 0mb に設定することで無効にできる大気屈折の補正についてです。pyephem で next_pass() を使用したパス予測を gnu predict からのパス予測と比較すると、AOS 時間と LOS 時間の間に 30 秒のオーダーの違いがあることに気付きました。大気屈折の補正が無効になりました。これは、predict と pyephem で使用される同一の TLE によるものです。

Q1) 最初に 30 秒程度の大きな差が生じたのは、通過予測で大気の屈折が補正されていないことが予測されたためではないでしょうか? これはもっともらしい説明のように思えます。

Q2) ほとんどの地上電波伝搬モデルは、地球半径の 4/3 のファッジ係数を使用して大気の屈折を補正し、曲線経路を効果的にまっすぐにして、計算された範囲に影響を与え、したがって範囲レートに影響を与えることを、私の (おそらく弱い) 理解です。pyephem の大気屈折補正は、このアプローチと同等ですか?

Q3) オブザーバーの圧力を 1010.0mb から 0.0mb に変更しても (したがって、大気屈折補正を無効にします)、range と range_rate の値は影響を受けていないように見えますが、next_pass() を呼び出した結果は影響を受けます。これは事実ですか？range と range_rate の値が変わると思っていたのですが、実際には、「見かけの距離」と「見かけの距離率」を、大気の屈折補正が適用されたときに計算できる追加の値と見なすべきでしょうか?

これは、スタック オーバーフローに関する私の最初の質問です。クエリが適切に形成されていない場合は、お詫び申し上げます。ありがとう。

更新 - 要求に応じた一部のコードと一部の結果

ここにいくつかのコードがあります...

ここにいくつかの結果があります：

Python で気象観測所用の小さなプログラムを作成しましたが、すべて問題なく動作していますが、小さな問題が 1 つあります。

出力は次のように表示されます。

のように表示させたい

これを行う簡単な方法はありますか?文字を削除することを考えていましたが、変換するよりもボッジのように聞こえます。

助けてくれてありがとう

私は pyEphem を使用して衛星の RA/Decs を計算していますが、計算されて記述されたさまざまな値に混乱しています

http://rhodesmill.org/pyephem/radec.html

このコードのビット

与える

a_ra= 0:52:40.75 a_dec= -3:15:23.7 g_ra= 1:14:10.55 g_dec= 0:06:09.8 ra= 0:53:23.57 dec= -3:10:50.5

オブザーバーの場所だけを変更して、gatech.lon、gatech.lat = '-5.47322222'、'19.82561111' と言うと

私は得る

a_ra= 1:15:36.95 a_dec= -2:32:29.9 g_ra= 1:14:10.55 g_dec= 0:06:09.8 ra= 1:16:19.75 dec= -2:28:04.6

オブザーバーの位置は sat.ra と sat.dec の計算に入るだけだと思っていたので、a_ra と a_dec が変更されているのを見て驚きました。

私は何が欠けていますか？

ありがとう広告

同じオブジェクト、同じ観測者の長さ、緯度でコードを計算すると、エフェムが 2014 年 4 月 6 日より前の日付の立ち上がり時間よりも前の設定時間を返すことに気付きました。

私のコードでは、sys.argv[1] の日付はローカルであり、emphem を通過する前に最初に utc に変換され、ep.localtime() を使用した後に出力を返します。

ここの出力が示すように: 1) 2014-04-06 の設定時刻は 2014-04-05 の 17:06:10 であり、上昇時刻は同じ日付ですが 23:16:29 です。2) 2014 年 4 月 7 日のコードを実行すると、設定時間が 4 分ではなく 8 分進みました。しかし、今では設定時間は立ち上がり時間と一致しています。

[13]: CygArecurazalt.py 2014-04-06 を実行

上昇:2014-04-05 23:16:29.000005 セット:2014-04-05 17:06:10.000005

[14]: CygArecurazalt.py 2014-04-07 を実行

上昇:2014-04-06 23:12:33.000006 セット:2014-04-07 16:58:22

誰も同じ問題に遭遇しましたか? それに対する解決策は何ですか？立ち上がり時間 (現地時間に変換) が sys.argv[1] 日付 (現地時間) よりも後になる方法はありますか?

私はおそらくユニークなユース ケースに取り組んでいます。Skyfield を使用して、架空の星系でいくつかの計算を行いたいと考えています。これを行うには、独自のエフェメリスを作成し、実際のエフェメリスの代わりにそれを使用します。私が見つけている問題は、エフェメリスを自分のものに置き換えるための API に関するドキュメントが見つからないことです。

ドキュメントはありますか？スカイフィールドは、私が試みていることを実行するのに十分な柔軟性を備えていますか?

編集：私が求めていることを明確にするために、重力モデリングを行う必要があることを理解しています（そして、この家のすべてのコンピューター、タブレット、ケーブルボックス、トースターを構成して、これらの数値を数日間処理することを完全に望んでいます： ) ですが、実際に詳しく説明する前に、データがどのように見えるかを知りたいと思いました。それが名前付きの numpy 2d 配列を含む単なるモジュールである場合...それはかなり簡単になりますが、これがどこにも文書化されているのを見たことがありません。