問題タブ [pyephem]

入力軌道要素値のセットと関連するエポックが与えられた場合、pyephem が特定の日付/時刻の「摂動」軌道要素値を計算する関数を提供するかどうかは誰にもわかりませんか? 現在、pyslalib パッケージの sla_pertel 関数を使用して軌道要素を摂動させていますが、それは正常に動作しますが、pyephem に同じことを行うために使用できるものが既にある場合は、そのパッケージを使用する必要はありません。私は pyephem のドキュメントを調べましたが、それを行うための明らかなものは何も見当たりませんでした。ありがとう。

私は PyEphem が初めてで、何ができるのか、どのように機能するのかを理解しようとしています。それをブラックボックスとして使用したくなく、得られた数値を盲目的に信頼したくないので、ここで説明されている例を再現したいと思いました。

この例では、1998 年 8 月 10 日 23:10 UT の特定のオブザーバーのオブジェクトの方位角と高度を計算します。次のパラメータが与えられます。

RA = 16時間41.7分、DEC = 36日28分

観測者の緯度は北緯 52 日 30 分、西経 1 日 55 分です。

例 (Excel で再現できます) によると、正解は AZ = 269.14634 度、ALT = 49.169122 度です。

同じ結果を得るために、pyephem を使用して次のコードを作成しました。

プログラムを実行すると、次の出力が得られます。

AZ + ALT の結果は明らかにこの例とほぼ同じですが、まったく同じではありません。また、RA と DEC が、入力したものと比較して印刷物でわずかに変更されているという事実にも困惑しています。

結果が異なる理由と、結果を再現するために私ができること、またはすべきことを誰かが明らかにするのを手伝ってくれるなら、私はそれを大いに感謝します. ありがとう。

編集：以下の回答で指摘されたタイプミスを修正しました。質問はまだ有効です。

EDIT2 ：OK、このリンクからPyEphemによって赤経と赤緯が調整される理由を読みました（そしてある程度理解しました）。私が理解していないのは、大気の屈折を無視させるのと同じように、PyEphem に光の相対論的たわみ、章動、収差の調整を無視させる方法があるかどうかということです。Azimuth の違いは RA と DEC の調整によるものだと推測していますが、確認していただければ幸いです。

次のコード フラグメントを実行します。

これにより、次の出力が得られます。

したがって、太陽オブジェクトの時間と太陽のオブザーバーについて計算された時間の間には 3 秒の差があります。2つのうち、どちらが最も信頼できますか? これらの値をローカルで計算された値と比較すると、オブザーバーの次の通過時間が最も近くなります。

火星についても同じ計算を行うと、次のようになります。

これにより、次の出力が得られます。

だから違いはありません。

よろしくお願いします、

大理石

観測者の地平線属性を調整しても、地球衛星の立ち上がり/設定時間には影響がないようです。さらに、私はこれをlibastroコードで見つけました：

ホライズンカットオフが実装されていない理由はありますか？たとえば、衛星の高度マスキングを行う場合に便利です。

太陽が次に沈む時刻を計算しようとすると、明日の値が得られますが、太陽が確実に地平線上にある間に計算しています!

m は私の観測サイトです。それが私が端末で得たものです：

それが明日の時間です。また：

したがって、previous_setting は未来です。

どこが間違っていますか？

次の TLE を計算しようとすると、計算結果が不完全になります。

この TLE は次を使用して読み取られます。

次のようにオブザーバーを設定しました。

私が使うとき

この関数は obj.ra や obj.elevation などに適切な値を入力できませんが、クラッシュしたりハングしたりしません (これもときどき発生します)。これはおそらく、オブジェクトが最近大気中に崩壊したという事実によるものですが、 、不完全な結果が私のプログラムの残りの部分をクラッシュさせます。この例外をキャッチする方法を見つけたいだけですが、失敗しました。

私は試した

も試した

これにより、~ 演算子が obj.elevation データ型 (またはそのようなもの) で使用できないことを示す別のエラーが発生しました。

Body.compute() に起因するこれらの例外をキャッチする方法はありますか?

オブザーバーが海面よりはるかに高くなる可能性があるという問題があります。RAとが与えられた場合Dec、この観測者とこのターゲットの高度と方位角を決定しようとしています。Observer()オブザーバーとFixedBody()ターゲットに使用しています。

高度を計算するために仮定する地平線はpyephem、実際の地球の地平線 (観測者にとって海面よりはるかに高い天頂にある) ですか? それとも、海面上の観測者の位置に描かれたローカル接線である地平線ですか? いくつかのテストから後者のようですが、これを確認したいと思います。

ephem初めて使用していますが、 oberver.sidereal_time() の出力を理解するのに苦労しています

時角から太陽時を判断するためのスクリプトをいくつか作成しました。1 つ目は、ephem を使用して赤経を計算し、Meeus の天文アルゴリズムの式を使用して、グリニッジ平均恒星時を取得します。これは、経度を使用してローカル平均恒星時に変換できます。

これにより、いくつかのデータをプラグインしたときに期待される出力が得られます。

私が書いた 2 番目のスクリプトは、同じ方法で赤経を計算しますが、ephem の sidereal_time() を使用してローカルの見かけの恒星時を取得します。

これは私が期待する出力を得られません。

私の知る限り、2つのスクリプトの唯一の違いは、最初の時角が地方の平均恒星時を基準にしているのに対し、2番目の時角は地方の見かけの恒星時を基準にしていることです。非常に小さな要因。代わりに、約 3 時間の違いが見られます。誰が私に何が起こっているのか説明できますか?

私は以前に pyephem を使用したことがなく、衛星測位の専門家でもありません。pyephem を利用して、TLE を使用して衛星の位置を計算したいと考えています。次のように、非常に簡単なことをしなければなりません。

-40:06:46.3 199:08:24.3 8834968.0

これらの 3 つの変数は、地平線基準系における衛星の位置を提供します。pyephem がこの値をどのように計算するかは明確ではありません。リファレンス ガイドを読みました: http://rhodesmill.org/pyephem/radec

ドキュメントを読むと、pyephem は歳差運動と章動を適用しているように見えますが、ドキュメントの最後の 2 行には次のように書かれています。

"最後の 2 つの座標セットのいずれにも歳差運動が適用されていないことに注意してください。最初の座標にのみ適用されます。これは、「アストロメトリック」位置のみが星図の線に対応することを意味します。他の位置は「 epoch-of-date」座標であり、観測自体のまさにその日の天の極と天の赤道の方向から測定されます。

az と alt に地球の歳差運動は適用されますか?

さらに、歳差運動と章動運動に pyephem が使用するモデルの種類を知りたいです (本当に参照が必要です)。Xephem と libastro へのリンクがありますが、アルゴリズムについては何も見つかりません。何か提案はありますか？

どうもありがとうございました！

私は python3 と pyephem を使用して天体航法を研究してきましたが、ほとんどの場合、手動で行う必要があるいくつかの計算を行っています。
視力低下の部分については、航海年鑑、主にグリニッジ時角と赤緯に匹敵するピエフェム出力データを取得しようとしています。牡羊座の最初の点の時角を取得するために、ephem.readdb を使用して 0.0 RA および 0.0 dec でボディを追加してみました。しかし、ドキュメントをさらに読んだ後、これは機能していると思います。

出力の最後の行を航海年鑑の日付の牡羊座の列と比較します

航海年鑑のオンライン版 www.tecepe.com.br