IDL memo
参考資料
コマンドメモ
- png 形式で書き出す
- ps 形式で書き出す
- Windowに描画
- 図にテキストを書く
- グラフのプロット
- ファイルへ出力(format文つき)
- ファイル読み込み(任意のフォーマットで)
- 分点変換 (B1950->J2000, in DEGREES)
- ガウシアンPSFを作る
- 画面左上を座標原点に取る
-
set_plot,'Z'
device, set_resolution=[600,400]
plot, (なんか適当に)
write_png, pngファイルの名前, tvrd()
-
set_plot,'ps'
device,filename='file.ps'
plot, (なんか適当に)
device,/close
-
device,decompose=0 ;(カラー表示の準備)
loadct, NUMBER; (blue-redは33番がいい感じ)
Window, xsize=.., ysize=.., retain=2; (retain=2で再描画する)
tvscl, ARRAY
-
xyouts, x, y, テキスト, charsize=1.0
-
plot, x, y, xrange=[xmin, xmax], yrange=[ymin, ymax],[/ynozero(縦軸範囲を自動調整)], xtitle='タイトル', ytitle='タイトル', psym=4 (マイナスの値なら線で連結)
-
printf, file_id, x,y, format='(A16, " ", f6.3," ", e11.4)'
-
template_name=ascii_template('filename');GUIが立ち上がってformatを指定
data=read_ascii('filename',template=template_name)
-
jprecess, ra, dec, ra_J2000, dec_J2000
-
psf=psf_gaussian(FWHM=FWHM,npix=Size of psf image)
-
!Order=1
FITS画像関連
- fits画像を開く
- fitsを読む
- fits header をもとに、astrometry情報を取り出す(構造体形式で出力)
- fits header をもとに、任意のピクセル(X,Y)の天体座標(RA,DEC)を出力する
- fits header をもとに、任意の天体座標(RA,DEC)からピクセル座標(X,Y)を出力する
- fits header を作成する
- fits header から項目を抽出する
- fits header に新しい項目の追加・変更する
- 座標を decimal から sexigesimal に変換する
-
fits_open, 'FILENAME', ID
-
fits_read,ID,data,header
-
extast,header,astrometry_info
-
xyad, header, X, Y, RA, DEC
-
adxy, header, RA, DEC, X, Y
-
mkhdr,newheader,image_array
-
result = sxpar(header, 'parameter')
-
sxaddpar, header, 'parameter', 'value'
-
radec, ra, dec, ihr, imin, xsec, ideg, imn, xsc