問題タブ [fits]

から提供された緯度に基づいて、fits 配列をスプライスしようとしていHeaderます。ただし、Python の知識と のドキュメントでは、そうすることができないようですastropy。私が持っているコードは次のようなものです：

latitude1これは度単位の値で、ヘッダーから呼び出された後に認識されます。だから私が達成する必要がある2つのことがあります：

1. 銀河緯度を認識するためにヘッダーを呼び出す方法は?
2. 緯度の範囲の値のみが含まれ、それ以外はすべて 0 になるように、配列をつなぎ合わせます。

satrapy と次のコードを使用して、FITS ファイルから世界座標系 (WCS) を読み取ろうとしました。

エラーが発生します：

ヘッダー ファイルは次のとおりです。

これをヘッダー ファイル内の他の WCS と比較したところ、見た目は異なりますが、問題を修正する方法は確かです。RADESYSaまた、それよりも使い方がわかりませんRADECSYS。私が見つけることができるドキュメントはないようです。どんな助けでも大歓迎です。

以前の投稿dtypeで、 の変更はrecarrayを使用して実行できることを見てきましたastype。recarrayただし、列の1つに配列がある でそれを行うことはできません。

Myrecarrayは FITS ファイル レコードから取得されます。

AC列をからintに変換する必要があるfloatので、試しました：

そして、それdtypeは確かに変わったように見えますが、

ACのデータを見ると、それらがまだ整数値であることがわかります。たとえば、sum計算が変数の限界に達したとしint16ます (すべてのAC列の値が正です)。

配列のデータ型を効果的に変更する方法はありますか?

私が使用しようとしているスクリプトでは、少し不可解なエラーが返されているようです。つまり、古いファイル拡張子から新しいファイルを作成しようとすると、作成しようとしているファイルが存在しません？なぜこれが起こっているのか、そしておそらく修正の可能性を説明できる人はいますか?

私が使用しているスクリプトは次のとおりです。

実行時:

次のエラーが表示されます。

完全なトレースバック:

astropy パッケージのファイル操作に問題があり、助けが必要です。

私は基本的に、ファイル形式に適合する画像を取得し、補正係数の入力を開始する必要がある新しいファイルと、補正係数と元の画像と共に使用して補正画像を生成できる新しい画像を作成したいと考えています。これらはそれぞれ同じ寸法になります。

これから始めます：

これを行うために欠けている明らかな知識を誰か教えてもらえますか...最後の2行は、何が欠けているかを概説するためのものです。含めました?? おそらく何か他のものが必要であることを知らせるために、おそらくfits.writeto()または同様のもの...