問題タブ [fits]

私はこの問題の解決策を探して探しましたが、何も見つかりません。

私は、matplotlib を介して長方形の FITS 画像を生成し、続いて AstroPy (または PyFITS) を使用してそれらに WCS 座標を適用しています。私の画像は銀河の緯度と経度であるため、私のマップに適したヘッダー キーワードはGLON-CARand GLAT-CAR(直交投影の場合) である必要があります。SAO DS9でこの同じ地図投影法を使用する他の地図を見てきましたが、座標はうまく機能します...グリッドは本来あるべき完全に直交しています。FITS の標準的な投影法は、ここで見つけることができます。

しかし、マップを生成すると、座標はまったくデカルトではありません。これは、私のマップ (左) とほぼ同じ地域の別の参照マップ (右) を並べて比較したものです。どちらも FITS ヘッダーにリストされていますGLON-CARがGLAT-CAR、私のものは SAO DS9 で見ると厄介です (座標グリッドは FITS ヘッダーのデータに基づいて SAO DS9 が生成するか、少なくとも FITS ファイルのどこかに保存されていることに注意してください)。

投影が間違っていると、座標割り当てアルゴリズムが各ピクセルに間違った座標を割り当てるため、これは問題です。

誰かがこれに遭遇したことがありますか、それとも何が問題なのか知っていますか?

他の投影法を適用してみましたが (SAO DS9 でどのように機能するかを確認するためだけに)、問題なく表示されます... しかし、デカルト投影とメルカトル図法は本来の直交グリッドでは表示されません。

これが AstroPy のバグであるとは信じられませんが、他の原因を見つけることができません...ヘッダーの引数の形式が正しくない場合を除きますが、それがどのように問題を引き起こす可能性があるのか​​ まだわかりません.私は経験しています。それとも、他のものを使用することをお勧めしますか? （私はmatplotlibベースマップを見てきましたが、それを自分のコンピューターで動作させるのに問題がありました）。

私のヘッダーコードは以下の通りです:

ご協力いただきありがとうございます。

私はPythonが初めてで、これにアプローチする方法がわかりません! FITS ファイルのディレクトリがあり、それぞれから特定のヘッダーを読み取り、テーブルとして出力できるようにしたいと考えています。

私は英語でアルゴリズムを構築する方法を知っており、個々の FITS ファイルからヘッダーを読み取る方法も知っています。ディレクトリ内の束全体からそれを行うには助けが必要です。

1. 最初lsにすべてのファイルを実行して表示します
2. どういうわけかfor、私のディレクトリの各ファイルを調べてそれをhdulist
3. のコマンドを与えるhdulist[0].header['name of the header I want'](プライマリのみを見て)
4. おそらくASCIIテーブルまたは通常のテーブル/テキストファイルで、それらすべてを印刷してください。

イベント データの fit ファイルがあり、同じテーブルの既存の列に格納されているデータから派生したデータの新しい列を追加して、テーブルの 1 つを変更する必要があります。私が抱えている問題は、変更されたファイルを閉じることです。これはコードです：

入力と同じファイルに書き込んでいるファイルを閉じるにはどうすればよいですか? 試してみるとtable.write(events, format='fits')、既存のファイルへの書き込みが原因でエラーが発生しますが、閉じようとするとdata、変更がファイルに書き込まれません。

はめあい形式の画像のコントロール ポイントの選択を見つけるのに問題cpselectがあります。このコマンドは、はめあい形式の画像を表示できません。一般的に、2つの画像を適切な形式でコレジストレーションしたいのですが、これまでのところ、その方法を理解できませんでしたがimtransform、「cp2tform」の使用方法は知っていました。

PyFITS を使用して、2 つの古い FITS ファイルから新しい FITS ファイルを作成しようとしています。

出力されるのは最初の値、つまりsci1[r][c]. これは、変数がまったく変更されていないことを意味します。どうすれば変更できますか？私はFITSの使用に非常に慣れていません。