yangtaozhi的个人博客分享 http://blog.sciencenet.cn/u/yangtaozhi

博文

利用IRAF中的软件包进行日心儒略日(HJD)时间的计算和提取

已有 4413 次阅读 2017-8-24 17:12 |系统分类:科研笔记| iraf, HJD计算和提取

在时域天文研究中,我们经常会涉及到时间的计算和提取,其中,对于物理变星中的脉动变星的时序测光过程中,我们通常用到的时间应该是由观测时的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



https://blog.sciencenet.cn/blog-807279-1072582.html

上一篇:利用IRAF对于盾牌座脉动变星YZ Boo的孔径测光笔记
下一篇:利用MESA进行星震学模型构建中铁氢比与金属丰度Z的转化关系
收藏 IP: 61.148.245.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-14 12:56

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部