カタログ確認

カタログを眺める

カタログデータの中身を確認するためにfvというソフトを使います。例えばvisit=950, ccd=49のSRCファイルを開くには次のようにします。

fv HSC/rerun/tutorial/00816/HSC-I/output/SRC-0000950-049.fits

すると、下のような画面が出てきます。

_images/fv-1.png

「1」をクリックするとカタログの内容を表示できます。

_images/fvb.png

「2」をクリックするとヘッダが確認できそこから各カラムの定義や単位などを知ることができます。

_images/fvc.png

テキストに書き出す

必要なスクリプトをダウンロードします。

cd /work/<User Name>
wget http://hsc.mtk.nao.ac.jp/HSC_Training_download/script/2018/fits2txt.py

このスクリプトはsingleFrameDriver.pyの結果のカタログのうち、天体の重心が問題なく決まったものの座標、PSF等級、PSF等級誤差をテキストに書き出します。 例えばvisit=950, ccd=49の結果について走らせるには次のようにします。

python fits2txt.py --data ./HSC/rerun/tutorial --visit 950 --ccd 49

# 他のカラムに興味がある場合はスクリプトを各自で変更して走らせてみてください。
# 日本語が文字化けしている場合は、
# メニューから「Terminal/Set Character Encoding/Unicord(UTF-8)」を選び、export LC_ALL=ja_JP.UTF-8 をしてください。

ds9上にオーバープロット

必要なスクリプトをダウンロードします。

cd /work/[User Name]
wget http://hsc.mtk.nao.ac.jp/HSC_Training_download/script/2018/overplot.py

このスクリプトはsingleFrameDriver.pyの結果のカタログのうち、天体の重心が問題なく決まった天体を画像上にオーバープロットします。 例えばvisit=950, ccd=49の結果について走らせるには次のようにします。

python overplot.py --data ./HSC/rerun/tutorial --visit 950 --ccd 49

または、測定された2次モーメントから楕円を表示したい場合は

python overplot.py --data ./HSC/rerun/tutorial --visit 950 --ccd 49 --ellip
表示される画像には色がついているはずです。色とflagの関係は以下に書いた通りです。
BAD=RED,
CR=MAGENTA,
EDGE=YELLOW,
INTERPOLATED=GREEN,
SATURATED=GREEN,
DETECTED=BLUE,
DETECTED_NEGATIVE=CYAN,
SUSPECT=YELLOW,
NO_DATA=ORANGE,
INTRP=GREEN,
SAT=GREEN

次に画像を見やすくします。

  1. メニューの「Scale/Zscale」を選びます。

  2. メニューの「Analysis/Mask Parameters...」を選びます。
    • マスクの透明度をあげてマスクの裏が透けて見えるようにします。

下のような図ができあがります。

_images/ds9.png