Pythonで特定のピクセルを中心に画像を回転するにはどうすればよいですか? 夜空の一連の画像の回転を解除しようとしています。星は北極星の周りを回転するので、北極星を回転の中心として定義し、各画像を回転させて星を並べることができます。

ephem を使用して日の出と日没を計算する例をここで見ましたが、うまく機能しています。

この 2 つの時間の中間点を計算しようとすると、問題が発生します。ここに私が持っているものがあります:

私は得る：

日の出: 2010/11/2 12:47:40 UTC
日没: 2010/11/2 23:24:25 UTC
正午: 5:18:22.679044

それが私が立ち往生しているところです。私は Python の初心者であり、率直に言って、一般的にはあまりプログラマーではありません。

どんな提案でも大歓迎です！

特定の暦日の日の出と日の入りの時間を教えてくれる Objective-C (または C) ライブラリ (コアの場所と互換性がある) はありますか?

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

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

ご挨拶、

私は、観測された空間の大きさとビッグバンからの経過時間を決定するのに役立つ研究を行っています。うまくいけば、あなたが助けることができます!

2 つの線形回帰を実行したい双一次データがあります。勾配と切片が変化するポイントがあり、このポイントを見つける (プログラムを作成する) 必要があります。

考え？

~ おおよその太陽正午

更新として、混乱の一部は、太陽正午が紀元前 4713 年 1 月 1 日グリニッジ正午以降のすべての計算が開始された場所であることです。

Date.ordinal_to_jd を正しく使用しても、この事実は補償されません。したがって、次のように 12 時間を加算または減算すると、次のようになります。

エラーが少なくなるはずです。昨日の正午から計算を開始するので、どちらを使用しますか?

このコードは、このページSunrise_equationの 2 つの方程式から派生しています。

ここでユーザーから得た最初の回答は、0.0009 と lw / 360 の使用法を理解していないというものでした。lw / 360 は、本初子午線からの円弧の分数日であるように見えます。0.0009 に関しては、紀元前 4713 年 1 月 1 日グリニッジ正午からの秒単位のわずかな差異に違いありません。詳細については、IAU 規格を参照してください

このページによると、0.007776秒と計算されます。

メソッドの詳細を含まない Date クラスからの情報が少しあります。

余談ですが、Ruby にユリウス日を計算する方法があるのは素晴らしいことです。NOAAの Javascript コードを調べています。

これは、リンクからインスピレーションを得て作成したクラスです。

これでそこにたどり着きましたが、一番上の式で計算されたような数値に戻る方法をまだ調査中です。これを試すと、JDN で 0.5 が得られることがわかります。誰が正しいですか？ルビーかNOAAか？

NOAA は、2000 年 1 月 1 日の値 2451545.0 を使用して、jd から減算して、このように分数の世紀で時間を取得します。

私のタスクは比較的単純です。入力ファイルの各行について、その行が特定の条件セットを満たしているかどうかをテストし、満たしていれば、その行の特定の列を新しいファイルに書き込みます。私はこれを行う python スクリプトを作成しましたが、1) 速度の向上、2) 列名に関して作業するための最良の方法 (列番号はファイルごとに異なる可能性があるため)、および 3 について助けてください。 ) フィルター条件と目的の出力列を指定する最良の方法です。

1) 私が扱っているファイルには、天体画像の測光が含まれています。各ファイルは約 1e6 行×150 列の float であり、通常、サイズは 1GB を超えます。このようなファイルを約 1 分で処理する古い AWK スクリプトがあります。私の python スクリプトには 5 ～ 7 分かかります。出力ファイルが目的のものになるまで、フィルタリング条件を微調整して数回再実行する必要があることがよくあるため、速度は間違いなく望ましいものです。for ループは非常に高速であることがわかりました。ループを遅くするのは、ループ内での処理方法です。itemgetter を使用して必要な列だけを選択することは、行全体をメモリに読み込むよりも大幅に改善されましたが、さらに速度を上げるために何ができるかわかりません。これは AWK と同じくらい高速ですか?

2) 特定の量 (フォトン数、バックグラウンド、シグナル対ノイズなど) の列番号はファイル間で変わる可能性があるため、列番号ではなく列名で作業したいと思います。私の AWK スクリプトでは、フィルタリングと出力が同じ量に適用される場合でも、条件と出力列が指定されている場所で列番号が正しいことを常に確認する必要があります。Python での私の解決策は、各数量に列番号を割り当てる辞書を作成することでした。ファイルに異なる列がある場合、新しい辞書を指定するだけで済みます。おそらくこれを行うためのより良い方法がありますか？

3) 理想的には、入力ファイルと出力ファイルの名前、フィルタリング条件、および出力する目的の列を指定するだけでよく、それらはスクリプトの先頭にあるため、検索する必要はありません。何かを微調整するためだけのコード。私の主な問題は、未定義の変数にありました。たとえば、一般的な条件は「SNR > 4」ですが、「SNR」(信号対ノイズ) には、測光ファイルからラインの読み取りが開始されるまで、実際には値が割り当てられません。私の解決策は、文字列と eval/exec の組み合わせを使用することでした。繰り返しますが、もっと良い方法があるのではないでしょうか?

私はコンピューター サイエンスの訓練を受けていません (私は天文学の大学院生です)。通常、何かを一緒にハックして、それが機能するまでデバッグします。しかし、上記の 3 点に関する最適化は、私の研究にとって非常に重要になってきました。長文で申し訳ありませんが、詳しく教えていただけると助かります。クリーンアップ/コーディングスタイルに加えて、あなたが私に持っているすべてのアドバイスをいただければ幸いです。

どうもありがとう、ジェイク

与えられた時間で周回する物体の真の異常を見つけるためのステップとして、ケプラーの方程式を解こうとしています。しかし、ケプラーの方程式は解くのが難しいことが判明し、ウィキペディアのページでは微積分を使用したプロセスが説明されています。微積分はわかりませんが、方程式を解くには無限の数の集合が必要であり、正解にますます近づく近似が生成されることは理解しています。

数学を見てもこれを計算的に行う方法がわからないので、数学のバックグラウンドが優れている人が助けてくれることを望んでいました。この獣を計算的に解くにはどうすればよいでしょうか?

FWIW、私は F# を使用しています -- そして、この方程式に必要な他の要素を計算できます。問題を抱えているのはこの部分だけです。

私はまた、与えられた時間、近点距離、および離心率を考慮して真の異常を近似する方法にもオープンです

特定の緯度、経度、時間の天頂付近で星を表示する天文学アプリケーションを構築しようとしています。

私は位置天文学の基本を認識していますが、私を困惑させているのは、膨大なデータセットからデータベースを作成することです!

私のクエリは次のとおりです。

1. どのカタログを優先しますか? ヒッパルコスのカタログやヘンリー・ドレーパーのような？その長所と短所は何ですか？
2. どのフィールドが必須ですか? いくつか特定しましたが、他に何を追加する必要があり、その理由は何ですか? 私が特定したフィールド: 星の名前、一意の星 ID、RA、Dec、等級、星の色。
3. ヒッパルコスのカタログからCommon Star の名前を取得するテーブルを見つけようとしました。どのテーブルにも見つかりませんでした。

ありがとう、ラージ

PyEphemを使用していて、影の長さを計算したいと思います（単位長の棒が地面に植えられていると仮定します）。長さはcot（phi）で与えられます。ここで、phiは太陽高度です（間違っている場合は訂正してください）。太陽でどのフィールドを使うべきかわかりませんか？以下の例では、角度altを使用しています。

以下の私の解釈を確認してください：

1. 接線が無限大の場合、太陽が真上にあり、影がないことを示します。
2. 接線がゼロの場合、太陽が地平線にあり、影が無限に長いことを示します。
3. cot（phi）の否定的な結果を解釈する方法がわかりません。誰かが私を助けることができますか？

最後に、ephem.Observer（）を指定して、PyEphemを使用して、影の長さから次に太陽がその長さの影を落とすまでの逆方向の作業方法について混乱しています。

これについて助けていただければ幸いです。