目的:用IRAF软件进行找星测光,并给出误差-星等图

下载数据

从SDSS网站https://dr12.sdss.org/fields,下载r波段数据 *.fits;

输入Run、Camcol、Field(或者RA、Dec):

hw1

选择r波段的fits文件下载:

hw2

打开数据获取信息

在 xgterm 中,利用 IRAF 及 ds9 打开下载的数据,并用 imexam 命令,给出天光值,天光的标准偏差值及星象的半高全宽值(多次测量取平均值给出即可)

进入存放fits文件的目录,在终端输入

$ ur_setup
$ ds9&xgterm

此时会弹出xgterm终端和ds9图像窗口

在xgterm终端输入:

$ cl
$ tv
$ display frame-r-007773-1-0166.fits[0]
#此时会在ds9中显示图像

hw3

hw4

再输入:

imexam frame-r-007773-1-0166.fits[0]

将鼠标的光标放在目标源的中心,按下键盘的a键,在xgterm中会出现信息(如下图所示),ENCLOSED是半高全宽

hw5

这里选取了5个星点,算得ENCLOSED的平均值为:2.1

再将鼠标移动到任意背景处,按下m键,在xgterm中会出现信息(如下图所示),STDDEV是标准差

hw5

这里选取了6个背景点,算得STDDEV的平均值为:0.037196

获取信息、记录完成后按q退出imexam模式。

自动找星

利用以上测量值,确定daophot 包里面的daofind 任务参数(请截图展示参数设置),并完成对 r 波段图像的自动找星,要求找到高于天光背景 5 sigma 的星,并给出找到星的数目,显示图像并利用tvmark 命令标记找到的恒星;

找星需要用到daophot包里的daofind,daofind通过以下命令打开:

$ noao
$ digiphot
$ daophot
$ daofind

将FWHM和std设置为上一步中得到的2.1和0.037196,并要求找到高于天光背景 5 sigma 的星:

hw5

此时目录下会生成frame-r-007773-1-0166.fits0.coo.2文件,打开后翻到底端,最后一个序号就是找到的星的数量。

在终端输入

$ tvmark 1 frame-r-007773-1-0166.fits0.coo.2 col=204 #col=204会在找到的星中心做一个红点标记

hw5

hw5

孔径测光

利用 daophot 包里面的 phot 进行孔径测光,测光参数基于以上测量值自行确定

终端中输入:

$ phot frame-r-007773-1-0166.fits[0]

hw5

hw5

此时会自动生成一个frame-r-007773-1-0166.fits0.mag.2的文件

终端中输入:

$ txdump frame-r-007773-1-0166.fits0.mag.2 "mag,merr"

xgterm中会显示两列数字,第一列是mag、第二列是magerr:

hw15

画出mag-merr图:

hw14

从以上星等-误差图可以看到,误差最大0.2等,且随着星等减小(源变暗),上图的误差上升趋势迅速增加。