問題タブ [fits]

私はpython2.5を実行していて、astLibライブラリを使用して天文画像のWCS情報を分析しようとしています。次のスケルトンコードでオブジェクトをインスタンス化しようとします。

ここで、file.fitsは、有効なフィットファイルを指す文字列です。

pyfitsヘッダーオブジェクトを渡す別の方法を使用しようとしましたが、これも失敗します。

エラーは次のとおりです。

ipythonで実行すると、pastebinで完全なエラーが発生します

astWCSモジュールがラップされたバージョンのWCStoolsであることは知っていますが、残りのコードはPythonであるため、Pythonモジュールを使用したいと思います。

誰かがこの問題を手伝うことができますか？

私が持っているいくつかのテスト データに合わせて、いくつかの FIT 画像ヘッダーのいくつかの値を変更する必要があります。したがって、アプリケーションで実行する分で FITs イメージ ヘッダーをハックしようとしています。

ただし、現時点では、ヘッダーも表示されません。ハックしてもかまいません。Ubuntuを実行しています。

FITを表示するようにソフトウェアにアドバイスできる人はいますか?おそらくそれをハッキングすることさえできますか?

FITSカスタム メタデータ (FITSヘッダー) 情報を含むファイルとしてデータをエクスポートした経験のある人はいますか? FITSこれまでのところ、標準の MathematicaFITSヘッダー テンプレートを使用したファイルしか生成できませんでした。ドキュメントには、カスタム メタデータ エクスポートがサポートされているかどうか、およびその方法についてのヒントはありません。

comp.soft-sys.math.mathematica からの次の提案は機能しません。

また

自分のメタデータをFITSファイルにエクスポートする適切な方法は何ですか?

乾杯、

マーカス

更新: Wolfram サポートからの応答: 「Mathematica は、FITS ファイルのメタデータのエクスポートをまだサポートしていません。この例では、このデータのインポートについて言及しています。将来的にはこれをサポートする予定です...」バイナリ テーブルを FITS インポート機能に追加します。」

私はいくつかの回避策を考え出そうとします。

私はPythonでのプログラミングが初めてで、Pythonで1つのコードで何千ものファイルを処理するセットアップを作成しようとしています。これを行うための練習用フォルダーを作成しました。その中には、2 つの FITS ファイル (FITS1.fits と FITS2.fits) があります。両方を .txt ファイルに入れるために次のことを行いました。

これが私が次にしたことです：

Enter キーを 2 回押した後、期待どおりに最初の図が表示されます。それから私はそれを閉じて、2番目を見ます。ただし、最初の図のみが保存されます。ループを変更して必要なことを行う方法について誰か提案がありますか?

先輩から教えてもらったデータがあります。データは何千もの .fits ファイルで構成されています。一部の .fits ファイルは他のファイルの古いバージョンであり、データ テーブルの作成方法が異なります。これが私が意味することです：

FITS1.fits と FITS2.fits の 2 つの .fits ファイルがあるとします。

データテーブルの列を切り替える方法を知っている人はいますか? FITS2.fits のフォーマットが FITS1.fits に似ているように？

FITSファイルヘッダーから情報を取得し、その情報をASCIIテーブルに転送する方法を学びたいです。たとえば、これは私が情報を取得する方法です。

RA_OBJとDEC_OBJを含むASCIIテーブルを作成するにはどうすればよいですか？

編集：ヘッダーから2つの列（RAとDEC）を含む.datファイルを作成したいと思います。これが私が試していることの例です：

次のようなものを含む.datファイルを作成したいと思います。

代わりに、私のコードは前のファイルのキーを上書きするだけです。何か案は？

FITS ファイルには、複数の画像 (より一般的には n 次元配列) とテーブルを含めることができます。Mac OS X Finder で FITS ファイルをダブルクリックすると、たとえばds9アプリケーションで画像を、 TOPCATアプリケーションでテーブルを自動的に開く方法はありますか?

AstroSoft の FITS Explorerのようなものを考えていますが、無料で Mac にうまく統合されています (つまり、Finder、Quicklook)。

そのようなものが存在しない場合、FITS ファイルの画像や表を参照して開くためのワークフローは何ですか?

2D numpy 配列の形式で 3 つのフィット画像があります。つまり、各ピクセルが 3 つの入力配列内の同じピクセルの中央値である出力配列を生成します。これは、imcombine を使用して IRAF で簡単に実行できます。配列全体をループして各ピクセルの中央値を取得せずに Python でこれを行う方法はありますか?

pyfitsOBJNO（オブジェクト番号）列で2つのデータオブジェクトを一致させようとしています。IDLでは、これはmatch, A.objno, B.objno, ii, jjプロシージャを使用して実行されます。iiこれにより、2つのインデックスと、jjどのインデックスAなどBが返されますA[ii].objno == B[ii].objno。

これを行うためのnumpy/pythonicの方法はありますか？Aどちらの順序やサイズについても、またはBそれらが大きなフィットのテーブルである可能性があるため、私は何の規定もしたくありません。これは私が今していることです：

より良いnumpy配列マッチングアルゴリズムはありますか？

通常の（512x512）2Dグリッドでの強度で構成されるデータキューブfitsfileがあります。ただし、3番目の軸もあります。この軸の各値は周波数に対応します。各周波数で、そのような2D画像が撮影されます。

私はfitsファイル内でこれらの頻度値を見つけようとしてきましたが、今のところ運がありません。私は140の異なる周波数チャネルがあることを知っています、私はフィットファイルのヘッダーを見ることによって単位などを知っています。問題は、それらの140チャネルの値を見つけることです。

ちなみにPYTHONでのプログラミング。