ftoolsの中でよく使ってたやつら
HEASoftをインストールするとデータの処理などでftoolsを使えるようになる。
ここでは私が使ってたftoolsのメモを公開する。
だいたいのことはターミナルで
$ fhelp ftoolsの各コマンド名
と入力すれば確認できる。
flcol
table のコラムを確認するときに使った。
$ flcol filename
例↓
$ flcol input.arf ___Column_Names_________Formats______Dims______Units___ ENERG_LO E keV ENERG_HI E keV SPECRESP E cm^2 //-- List FITS table column information.
表で例えると
順位|国名 |人口 ←これがコラム
1 |中国 |13億9,273万
2 |インド |13億5,261万7,328
3 |アメリカ |3億2,716万7,434
4 |インドネシア |2億6,766万3,435
5 |パキスタン |2億1,221万5,030
6 |ブラジル |2億946万9,333
7 |ナイジェリア |1億9,587万4,740
8 |バングラデシュ|1億6,135万6,039
9 |ロシア |1億4,447万8,050
10 |日本 |1億2,652万9,100
引用:(キッズ外務省)人口の多い国|外務省
fplot
コラムを選んで二次元プロットを書く。
Name of FITS file and [ext#][] input.arf Name of X Axis Parameter[error][ENERG_LO] Name of Y Axis Parameter[error] up to 8 allowed[SPECRESP] Lists of rows[-] Device: /XWindow, /XTerm, /TK, /PS, etc[/xw] Any legal PLT command[-] Illegal command. At token = - Type HELP to get command list.
fv
ファイルの中身の情報やtableなどが見れる、修正も出来る。
HEASoftのインストールでヘマすると、開くまで長い時間がかかる場合もある。
例↓
$ fv input.img
fdump
fvコマンドで見れるfileの情報をターミナル上で見れる。
例↓
fdump input.arf Name of optional output file[STDOUT] //outfileを設定する。”STDOUT”を入れるとターミナル上に表示される Names of columns[STDOUT] - //欲しいパラメータの名前をいれる “-“で全部表示される Lists of rows[-] //まだ理解してない
fdiff
2つのファイルを比べるときに値の異なるパラメータを確認できる。
fdiff inputA.img inputB.img FDIFF V1.1.4 inputA.img[0] < inputB.img[0] > < card 264: CRVAL2 = -1.0567010710700E+01 > card 264: CRVAL2 = -1.0570837710700E+01 --- < card 57: DEC_NOM = -1.0567010710700E+01 > card 57: DEC_NOM = -1.0570837710700E+01 ========================================= **** End of file comparison: 2 differences were found
ftimgcalc
画像ファイルに領域ファイルを読み込ませてピクセルの値を操作する感じ。
切り出したい領域以外のピクセルの値を0にするときに使った。
Output image name [inputC.img] Calculation expression [] regfilter("input.reg",A.P1,A.P2)? A:0 //A.P1,A.P2については理解してない(Aがimageのpixelごとのfluxに相当していて、region以外の部分を0にしろ、と言っているはず) Image letter A (filename, or NAME=filename) [NONE] a=output.img Image letter B (filename, or NAME=filename) [NONE]
f2dhisto
イベントファイルから画像ファイルを抽出する。
Name of FITS file and [ext#][] input1.evt Name of output FITS file[] output1.fits X BINSIZE (0.0:) [0.01] Y BINSIZE (0.0:) [0.01] Name of x column to generate histogram[finalxpos] Name of y column to generate histogram[finalypos] Enter the X-RANGE(Eg: 10.0,20.0)[-30,30] Enter the Y-RANGE(Eg: 10.0,20.0)[-30,30] input file: file://fti_20.evt(a.fits)[1][binj finalxpos=-30:30:0.01, finalypos=-30:30:0.01] . $ ds9 output1.fits
faddcol
infileの中にcolfileの中身のパラメータ(colnameで指定)を追加する。
faddcol infile=test.arf colfile=input.fits colname=fudge
fcalc
inputしたファイルの中でパラメータの計算して計算結果を新しいコラムとして追加、もしくは既存のコラムを上書きする。
fcalc test.arf output.arf clname=SPECRESP expr="SPECRESP/fudge”
fselect
eventファイルに条件を与えてファイルを出力する。
下の操作だとPHA幅や領域を切り出す。
$ fselect Name of FITS file and [ext#][] A.evts Name of output FITS file[] B.evts Selection Expression[] PHA>120&&PHA<200
領域で切り出すなら
Selection Expression[] ( RAWX-976.5)-(RAWY-1060.5)<10000
f2dhisto
evtファイルからimageファイルを表示する。
Name of FITS file and [ext#][] AAA.evts Name of output FITS file[] BBB.img X BINSIZE (0.0:) [16] 4 Y BINSIZE (0.0:) [16] 4 Name of x column to generate histogram[RAWX] Name of y column to generate histogram[RAWY] Enter the X-RANGE(Eg: 10.0,20.0)[1,2040] 1,1240 Enter the Y-RANGE(Eg: 10.0,20.0)[1,2048] 1,1156
ftlist
XMMのcalibulation fileの中身に欠陥があるかどうか確認するときに使っていた。
$ ftlist EMOS1_ADUCONV_0097.CCF Print options: H C K I T [H] Name Type Dimensions ---- ---- ---------- HDU 1 Primary Array Null Array HDU 2 ADUCOEFF BinTable 4 cols x 14 rows HDU 3 OFFSET_GAIN BinTable 13 cols x 7 rows Dumping CFITSIO error stack: -------------------------------------------------- Extension doesn't start with SIMPLE or XTENSION keyword. (ffrhdu) TENSION= 'BINTABLE' / binary table extension B Failed to move to HDU number 4 (ffmahd). -------------------------------------------------- CFITSIO error stack dump complete. CFITSIO ERROR UNKNOWN_REC: 1st key not SIMPLE or XTENSION Task ftlist 1.01 terminating with status 252