[发明专利]测量FGK型恒星金属丰度的方法及滤光片V有效
| 申请号: | 201810162482.9 | 申请日: | 2018-02-27 |
| 公开(公告)号: | CN108469418B | 公开(公告)日: | 2020-01-10 |
| 发明(设计)人: | 范舟 | 申请(专利权)人: | 中国科学院国家天文台 |
| 主分类号: | G01N21/31 | 分类号: | G01N21/31;G06T5/50 |
| 代理公司: | 11401 北京金智普华知识产权代理有限公司 | 代理人: | 巴晓艳 |
| 地址: | 100024 北京市*** | 国省代码: | 北京;11 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 恒星 丰度 金属 测量 滤光片 拟合 观测 滤光片带宽 大气消光 恒星观测 流量标准 天文测量 吸收谱线 有效测量 准确测量 数据处理 大样本 可测量 测光 定标 滤光 视场 覆盖 吸收 统计 | ||
1.一种测量FGK型恒星金属丰度的方法,其特征在于,包括如下步骤:
步骤一、对待测FGK型恒星进行观测和数据处理;对待测FGK型恒星进行观测时,采用对FGK型恒星金属丰度进行有效测量的滤光片vSAGE,该滤光片vSAGE中心波长位于395nm,滤光片带宽29nm;
步骤二、对流量标准星进行观测和处理;
步骤三、拟合大气消光曲线,并对待测恒星进行流量定标;
步骤四、对“(vSAGE-g)-(g-i)颜色~金属丰度关系[Fe/H]”进行拟合;
步骤五、通过测光方法求出待测恒星的金属丰度;
其中:
步骤一具体包括:
步骤1.1)获取望远镜的本底图像、多波段的平场图像、暗场图像;所述多波段为vSAGE/g/i波段;
步骤1.2)利用望远镜对待测FGK型恒星进行观测,观测时使用vSAGE波段滤光片以及g/i波段滤光片分别进行观测;控制曝光时间,实时检查星像的点扩散函数PSF轮廓,保障其不拉长、不出现虚焦;
步骤1.3)合并本底图像,对所有待测FGK型恒星的观测图像进行本底改正,消除探测器及光电观测设备的偏置电压及不均匀性影响;对于观测的待测恒星图像,要减去合并好的当晚的本底图像;
步骤1.4)进行暗场改正,消除探测器上暗电流的影响;
步骤1.5)进行平场改正,消除光学系统照明和CCD探测器像素响应差异的影响;
步骤1.6)对待测恒星进行vSAGE波段和g/i波段观测图像测光;
步骤二具体包括:
步骤2.1)对流量标准星进行vSAGE波段和g/i波段观测,观测中要选择测光夜,大气整晚稳定、无云的夜晚;控制曝光时间,实时检查星像的点扩散函数PSF轮廓,保障其不拉长、不出现虚焦;
步骤2.2)对流量标准星vSAGE波段和g/i波段的观测图像进行本底、暗场和平场的改正;
步骤2.3)对流量标准星进行vSAGE波段和g/i波段的图像进行测光,获得其仪器星等;
步骤三具体包括:
步骤3.1)计算流量标准星的大气质量,并拟合大气消光曲线;
大气质量的计算公式:
vcali-vinst=k0,v+k1,vχ+k2,v(vcali-gcali) (4)
gcali-ginst=k0,g+k1,gχ+k2,g(gcali-icali) (5)
icali-iinst=k0,i+k1,iχ+k2,i(gcali-icali) (6)
在上面的公式中,k0,v,k0,g和k0,i是仪器零点,k1,v,k1,g和k1,i是主消光系数,k2,v,k2,g和k2,i是二级消光系数和颜色相关,χ是大气质量,可用如下方法进行计算:
χ=sec z (7)
其中z是天顶角,该公式适用于天顶角z<60°-75°的范围;
通过测光夜多波段观测标准星和拟合该消光曲线对当晚的恒星观测数据进行流量定标;
步骤3.2)利用大气消光曲线对待测恒星进行流量定标:
将待测量的FGK型星的vSAGE/g/i波段的仪器星等带入步骤3.1)中拟合得到的大气消光曲线,得到大气外的定标星等;对于每个波段的仪器星等要带入对应波段的大气消光曲线;
步骤四具体包括:
步骤4.1)从Kurucz恒星大气模型选出FGK型恒星光谱样本卷积得到的vSAGE/g/i波段星等;
步骤4.2)拟合Kurucz恒星大气模型(vSAGE-g)-(g-i)颜色~金属丰度[Fe/H]的关系;
步骤五具体包括:
步骤5.1)将观测定标后的vSAGE/g/i波段FGK型星的星等带入拟合得到的(vSAGE-g)-(g-i)颜色~金属丰度[Fe/H]关系,从而得到其金属丰度;
步骤5.2)对得到待测FGK型星的金属丰度结果进行整理,按照需求生成固定格式的表格。
2.如权利要求1所述的方法,其特征在于,步骤4.1)中,通过如下公式卷积Kurucz恒星大气模型得到恒星的vSAGE/g/i波段星等:
对于Vega系统卷积公式如下:
MAGVega(Obj)=-2.5log(∫{Fl(Obj)Sλdλ}/∫{Fl(Vega)Sλdλ}) (8)
MAGVega(Obj)=-2.5log(∫{Fν(Obj)Sνdν}/∫{Fν(Vega)Sνdν}) (9)
对于AB系统卷积公式如下:
MAGAB(Obj)=-2.5log(∫{Fν(Obj)Sνdν}/∫{Sνdν})-48.6 (10)
其中λ是波长,ν是频率,Fν(Obj)是所述恒星的流量,单位是ergs s-1cm-2Hz-1,Sν是仪器的相应曲线,Fν(Vega)是织女星的流量,单位是ergs s-1cm-2Hz-1。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院国家天文台,未经中国科学院国家天文台许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810162482.9/1.html,转载请声明来源钻瓜专利网。





