問題タブ [fits]

バイナリ テーブルを含む 60GB の FITS ファイルがあります。このテーブルを一度に 1 つの行/エントリ/行/ブロック* で読み取り (および処理) したいと思います。

(※正確な命名法はわかりません)

私はpyfitsを使用していますが、やりたいことは簡単に言えば次のとおりです。

変数entryの型<class 'pyfits.fitsrec.FITS_record'>は、バイナリ テーブルの列数と同じ長さです。ただし、次の行でバイナリ テーブル全体がメモリに読み込まれるように見えますentry = hdulist[1].data[i]。

pyfits のドキュメントを調べましたが、テーブル エントリごとに (または一度に小さなエントリ セットで) バイナリ テーブル拡張からデータを読み取るようなメソッドが見つかりません。テーブルから特定のエントリを選択するのではなく、単に順番にスキャンするだけです。

私の質問は次のとおりです。

0)hdulist[1].data[i]ステップで何が起きているか? すべてがメモリに読み込まれるのはなぜですか? （これを回避する方法はありますか？）

1）私は何かを逃しましたか？実際にpyfitsは私が望むことをすることができますか？

2）別のpythonライブラリはありますか？(つまり、FITS 拡張でバイナリ テーブルを使用する)

3）そうでない場合、データを別のバイナリ（または他の圧縮/非ASCII）形式（FITSではない）で書き直して、他のpythonライブラリまたはモジュールを見つけて、やりたいことを実行できますか？

望遠鏡でいくつかの FITS 画像を撮影しましたが、画像を分析して各画像の星の数を取得するには、Python コードを使用する必要があります。私はPythonに非常に慣れていないため、まだコードを持っていません。どんな助けでも大歓迎です。

pyfits で FITS ファイルを開きました。HEADER ファイルは、、、、、などの36列のタグで読み取りXTENSION='BINTABLE'ます。DIMENSION= 52989R x 36C'ZBEST''ZQUALITY''M_B''UB''PGAL'

ここで、「ZQUALITY」が 2 より大きく、「PGAL」が 3 に等しいデータからオブジェクトを選択する必要があります。次に、上記の条件に従って、対応するオブジェクトの「ZBEST」のヒストグラムを作成する必要があります。また、これらのオブジェクトの「M_B」と「UB」をプロットする必要があります。

最後に、「ZBEST」を 3 つのスライス (zbest < 0.5)、(0.5 < zbest < 1.0)、(zbest > 1.0) にスライスし、それらのヒストグラムと「M_B」対「UB」図を別々にプロットしたいと考えています。

2 つの条件に従うデータの選択に行き詰まっています。両方の条件を満たすデータからオブジェクトを選択する方法を教えてください ('ZQUALITY' > 2 & 'PGAL' == 3 )? like: を使用していますdata.field[('zquality' > 2) & ('pgal'==3)]が、機能していません。

for ループで計算される一連の画像を作成する必要があります。を使用してパスを定義するとsprintf、エラーが発生します

これを回避する方法について何か提案はありますか? 画像が計算されるコード セクションは次のとおりです。

Pythonで.fitsファイルを読み取ることができるように、.mapファイルを.fitsファイルに変換する方法を誰かが知っているかどうか疑問に思っていましたか?

.csh でコードを書いています。FITSファイルの bunit ヘッダーをK (ケルビン) から km/s に変更しようとしています。どうやってやるの？

Python ではnew_fitsfile.header['BUNIT']='km/s' を使用することはわかっていますが、現在の .csh コードでは機能せず、Python コードに切り替えるオプションではありません。

多くの列を含む FITS ファイルがあります。ファイルの単純化された例の一部は次のようになります。

ここで気が付くと、列A と Cの最初の 4 行は同じですが、列Bは 1 から 4 まで変化します。次に、列A と Cの次の 5 つの値は同じですが、列Bは 40 から 44 まで変化します。

私がやりたいことは、次のようなファイルを作成するプログラムを書くことです:

つまり、列 B には、列 A と C が同じ値の数が含まれている必要があります。

Pythonでこれを行う方法を知りたいです。必ずしも FITS ファイルを扱う必要はありません。FITS ファイルで使用できないルーチンがいくつかある場合は、ASCII 形式に変換することもできます。

私がこれまでに試したこと：

リスト内の等しい値の数をカウントし、それらを返すCollectionsサブルーチンが呼び出されたルーチンにぶつかりました。Counter

私は試した：

しかし、これは列Aの値が等しいだけです。このルーチンを使用して、これらの値を列Cと比較する方法も教えてもらえますか?

SNOOPY (点広がり関数フィッター) と IRAF を使用して、フィット画像に PSF フィッティングを実行しようとしています。imexam を使用してこの画像を開くことはできますが、ポイント (星など) を選択するとエラーが発生します。

私が見ているものと IRAF が見ているもの (舞台裏) は同じではないようです。あたかもある種の調整オフシフトか何かがあるかのように。

これを修正するにはどうすればよいですか？

[Scientific Linux 6、16ビット、IRAF v2.16]