||
在时域天文研究中,我们经常会涉及到时间的计算和提取,其中,对于物理变星中的脉动变星的时序测光过程中,我们通常用到的时间应该是由观测时的UTC时间改正到日心的日心儒略日时间,即HJD时间。但是对于初次使用IRAF的用户来说,要想顺利的利用IRAF中的软件包进行HJD时间的转化和提取,还需要进行一些必要的设置,例如拍摄的CCD图片的台址信息等。
以下就以新疆天文台南山1米望远镜拍摄的盾牌座脉动变星YZ Boo为例,简要介绍下这个计算和提取的过程,希望能够对感兴趣的童鞋有所帮助。
第一步是要在IRAF的安装目录中找到相应的望远镜台址的信息的数据文件,一般情况下这个文件位于**/iraf/iraf/noao/lib这个目录下(注意,不同的用户安装的目录可能不尽相同,但是一般都是可以找到类似上面星号后面的那个目录),相应的望远镜的台址信息位于这个目录的obsdb.dat文件中,打开这个文件,会看到好多世界上不同的望远镜的信息,
原来没有添加新的望远镜台址信息的这个文件,其最开头的是kpno的信息,我们按照她的标准格式把新疆南山1米望远镜的信息添加了进去,注意,看原有的信息,我们发现,我们的经度是和美国的相反的,所以经度的前面是有负号的,而且时区也是要有负号的。如此,我们添加好,进行保存就好啦。
第二步,进行IRAF中的参数设置,在noao/astutil这个包里进行设置,epa setjd
其中第一个参数images设置为相应的fits图像,如images= YZ*.fits
第二个参数observa=XAONS 这个即是我们刚才添加进去的台址
下面的参数依次设置为相应的fits头文件中的关键词即可,注意有可能不同的望远镜拍摄的fits图里面的关键词不完全一致,此处设置时最好能够打开头文件,利用命令head YZ*fits即可。
第三步,计算HJD时间,设置好之后,利用命令setjd YZ*.fits 即可显示出相应的jd,hjd和ljd时间,如果想输出到某个文件,则可利用setjd YZ*.fits >time.dat 相应的时间信息均以被储存到了time.dat文件中了。
此后,结合孔径测光给出的较差星等便能够得出一条完整的光变曲线,方便后续的深入分析。
若有问题,随时欢迎交流。yangtaozhi2018@163.com
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-21 22:21
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社