IRAF手动测光

1、图像切边查看一下待处理的CCD图片存储位置：/home/yue/ct/2016.09.01CT图片包括：平场flat*、本底bias*以及目标星图VZ_Psc-Bessell-*.fits如图所示在终端中依次输入命令：lscdyuelscdctmkiraf。在EnterterminalType:中输入“xgterm”，然后发现ct中多了“login.ct”和“uparm”两项。在终端命令lhp@LHP:~/yue/ct&中输入xgterm&，启动xgterm，如下图所示。在Xgterm中依次输入cl、noao、imred、ccdred设置ccdred参数在ccdred后输入epaccdred，打开参数设置窗口，并修改相关参数。显示设置输入epadisplay,修改参数。删除提示设置epadel在terminal中输入ds9&,打开ds9窗口在xgterm中依次输入ls、cd2016.09.01CT输入命令“lsflat*”,显示平场打开任意一幅flat,如输入命令displayflat-Bessell-0004.fits，将指针放在图像的边缘上，确定大概区域范围X(3.5-1021）Y（3.5-1021）Implot命令打开对应display打开的图像，即输入“implotflat-Bessell-0004.fits”设置步长shift+:step1。再确定y的准确范围，输入y的单个坐标shift+:l1021,按J或K调节，空格键寻找临界点（突变点）对应的位置1，shift+:3.5相同方法得到位置2。同理做shift+:c3.5/1021，得到X方向区域，按q键退出。空格键可以读出value值设置切边参数输入命令：epaccdproc修改参数填写切边范围Shift+:go，系统自动切边2、本底处理求出该CCD的增益（gain）和读出噪声（readnoise）,display某个Bias，取没有坏点的小块区域【X1：X2，Y1：Y2】，两坐标相间50-100个像素，一般选四个顶角往里靠的四个方块区域，记下各方块的范围【X1：X2,Y1：Y2】在xgterm窗口依次输入noaoobsutilepafindgain，打开参数设置窗口。By可以用来返回上一级。设置参数设置完后，Shift+:go运行，记下运行结果中的Gain和Readnoise.同理再做（B/R/I）波段的Gain和Readnoise读值，得到16组数据，最后求平均值！合并bias设置参数ccdredepazerocombine，在rdnoise和gain处要填写相应的值Shift+:go运行平场减bias参数设置：epaccdproc补充：（zero=bias20160901.fit）Zerolevelcalibrationimage目标星减bias在上一步平场减bias设置基础上，只需改动images4、平场处理对各个波段的平场进行合并，目标星图对应波段减掉处理得到的平场1、对V波段平场合并，设置参数epaflatcombine,记得填写rdnoise和gain.Shift+:go目标星图归一化处理得到的对应V波段的平场参数设置：epaccdproc按键设置另一页参数5、去除宇宙线分波段做统计标准偏差，displayv波段的一张目标星图，调节ds9窗口中选择框形选区：Regionshapebox,选取均匀的没有背景星且离目标星较近的相差50-100像素的范围并记下坐标——[450:523,115:184]分波段统计标准偏差如V波段：crutilimstatisticVZ_Pse-*sell-V*[450:523,115:184]记下STDDEV值2、分波段去除宇宙线crutilepacosmicrays六、手动测光1、display显示4张目标星图，从中选取星相较大的。【ds9中zoominvertx或inverty使图像与参考图一致，zoomzooom8放大8倍，找到直径最大星的中心坐标并记下】（404,145）用implot打开以上对应图像，分别输入上面的y坐标值：激活irafterm窗口，shift+:l145(l对应y坐标），设置步长：shift+:step1(即步长为1）鼠标指针先放在图像的坐标下左面按下E键，再把鼠标指针放在图像的坐标上右面按E键，下左上右应等宽放置才保证放大后的图像放在坐标中间按J或K来调整最大峰值，鼠标指针放在波峰处按空格读取最大峰值value=10604.26读取最大峰值value=10604.26的一半左右对应的半宽值所对应的坐标X1，X2,两者相差的绝对值为星的半宽值在波的两边较平缓处按空格读取两坐标X3，X4,两者相差的绝对值为星的直径。同理做shift+:c404,做X坐标方向的半宽和直径，得到的两个半径和直径，选较大的。鼠标放在irafterm窗口，按Q键退出。（c对应x坐标）设置查看头文件信息进入digiphotapphotepaimheadShift+:go运行，下图为运行后的文件信息。将鼠标放在xgterm窗口，关闭数码锁定键，按下ctrl键的同时点击鼠标右键，可以打开xgterm窗口中的滚动条，来查看完整文件信息。路径noaodigiphotapphotepaphot，设置参数。然后经光标放在datapar处按Shifi+:e进入TASK：datapar参数设置进入datapars设置参数，如下图所示，设置完成后按shift+q退出同样按shift+e进入centerpars参数设置。设置后按shift:q退出进入fitskyp，设置参数进入photpar设置参数返回到phot参数设置窗口，将interac和radplot设置为yes，再按shift+:go进行测光测光参考目标星图，确定测光顺序，用implot工具来读取确定要测光的目标星、参考星和对比星（至少要4颗星），value值高于6万的星要舍掉一窗：apphotdisplayVZ_Psc-Bessell-000*.fits序号（显示9幅图像）二窗：apphotphotVZ_Psc-Bessell-000*.fits（每次自动从1—9）DS9中显示的图像应与xgterm窗口中提示要处理的图像名一致，如不一致在ds9中按Framefirst/previous/next/last来调节。激活ds9窗口，把鼠标放在要测光的星的中心同时按下空格键读取数据，按顺序测完后Q键退出，此时鼠标退到2窗，激活2窗（n代表下一幅，同时在ds9中也要点next进入下一幅，退出时要关闭大写锁定键）分波段读取数据apphotepapdump进行如下参数设置，shift+:q退出设置文件名apphotpdump*sell-V*mag.1VZ_Psc20160901v（星名+观测日期+波段）查看读取的数据apphot!geditVZ_Psc20160901v