||
具体处理过程:
1.进入IRAF setjd设置:ecl>noao>astutil>epar setjd
注: (observa= )输入观测站的名称,名称包含在obsdb.dat中,在终端查找
astro@astro-desktop:~$ locate obsdb.dat
/iraf/iraf/noao/lib/obsdb.dat
astro@astro-desktop:~$ sudo gedit /iraf/iraf/noao/lib/obsdb.dat (字符编码选择: 西欧)
打开可以看到如下的信息(列出部分),查看你的观测站名称是否包含其中,如果包含可以直接在(observa= )中输入简称,没有可以手动添加,下面的第一条‘konkoly’的信息就是我手动加入的,然后保存。
---------------------------------------------------------------------
# Observatory Parameters.
# Some were taken from the Almanac and others have been submitted.
# Observatories wishing to be added or make changes in the default
# distributed database should send information to iraf@noao.edu.
observatory = "konkoly"
name = "Hungarian Academy of Sciences, Konkoly Observatory"
longitude = -19:00.7 #观测站赤经
latitude = 46:10.8 #观测站赤纬
altitude = 160. #观测站海拔高度
timezone = -1 #观测站所处时区(注:这里是匈牙利时区 +1,输入为 timezone= -1)
#timezone只是当地时间和世界时间的转换,如果前面的utday 和uttime为UT时间(yes),timezone不起转换的作用
# Submitted by Biwei JIANG 11/28/95
observatory = "BAO"
name = "Beijing XingLong Observatory"
longitude = 242:25.5
latitude = 40:23.6
altitude = 950.
timezone = -8 #观测站所处时区(注:这里是中国时区 +8,输入为 timezone= -8)
# Submitted by Pat van Heerden 11/20/06
observatory = "boyden"
name = Boyden Observatory
longitude = 332:35:40.0
latitude = -29:02:20
altitude = 1387
timezone = -2
------------------------------------------------------------------------------
2、运行计算
astutil> setjd kasc*.fits > time.dat #将数据输入time.dat文件中
打开time.dat得到如下结果:
# Image jd hjd ljd
# SETJD: Observatory parameters for Hungarian Academy of Sciences, Konkoly Observatory
# timezone = -1
kasc-f2_V_0001.fits 2454522.27453 2454522.27606 2454522
kasc-f2_V_0002.fits 2454522.27614 2454522.27767 2454522
kasc-f2_V_0003.fits 2454522.27775 2454522.27929 2454522
.....
参考IDL计算程序: http://idlastro.gsfc.nasa.gov/ftp/pro/astro/observatory.pro
检验计算的HJD: http://www.physics.sfasu.edu/astro/javascript/hjd.html
2013-12-05于BNU
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 23:00
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社