BFOSC长缝光谱自动处理软件.pdf

.天 文 研 究 与 技 术 第 卷 第 期 年 月:./.长缝光谱自动处理软件宋得阳 张君波 王 靓(.中国科学院南京天文光学技术研究所 江苏 南京 .中国科学院天文光学技术重点实验室(南京天文光学技术研究所)江苏 南京 .中国科学院大学 北京 .中国科学院国家天文台 北京.中国科学院光学天文重点实验室(国家天文台)北京)摘要:在天文观测中 长缝光谱仪是一种常用的光谱观测设备 国家天文台兴隆观测基地 望远镜搭载的北京暗天体摄谱仪()的长缝光谱模式应用广泛 但没有配备专用的数据处理软件 通常用户观测后采用()软件利用交互方式处理 介绍了一款针对 长缝光谱的自动化数据处理软件 数据处理过程包括本底()和平场()改正、背景改正、抽取光谱和波长定标等 展示了对实测数据进行处理的结果 并与 软件处理的结果进行了对比 该软件在完成所有光谱处理必需步骤的基础上实现了自动化运行 并采用优化抽谱方法 能够自动扣除大部分宇宙线 保证了高精度的输出结果关键词:长缝光谱 优化抽谱方法 自动处理程序中图分类号:文献标识码:文章编号:()我国自主研制的第 台 级望远镜 国家天文台兴隆基地的 望远镜 配备了北京暗天体摄谱仪 这台仪器具备 种不同的观测模式 其中长缝光谱是应用最为广泛的一种模式 可以在不同的缝宽、光栅和滤光片组合之间进行切换 年更换 相机后最新的色散、波长范围等参数如表 据公开资料统计 年兴隆观测基地为 望远镜用户分配的 个观测夜中 使用 的有 夜 占比超过 测光模式下的数据可以使用文开发的程序处理 但长缝光谱部分尚没有专门的数据处理软件 观测者一般自行使用 等软件利用交互方式对原始光谱图像进行处理 获得一维谱 但这类软件并不是针对 特点设计的 没有针对长缝光谱的特点进行优化且美国国立光学天文台()已经宣布自 年起不再对 进行官方支持本文针对 长缝光谱 开发了一套光谱处理软件 实现了自动化处理光谱数据 该程序可以读取用户指定的光谱仪原始数据 自动进行图像合并、波长定标、目标源光谱抽取等步骤 在关键步骤中保存中间结果并进行可视化表 不同光栅及基本性能参数 ./()./()./基金项目:国家自然科学基金()资助.收稿日期:修订日期:作者简介:宋得阳 男 硕士 主要研究天文技术与方法.:.通信作者:王 靓 男 副研究员 主要研究天体物理.:.天 文研究与技术 卷程序引入了优化抽谱方法 自动屏蔽了宇宙线 提高了最终数据产品的精度 用户可以使用该程序快速实现批量数据处理 大大减轻了用户的工作量 有效提升了工作效率 系统整体框架 环境配置本程序基于 编程语言编写 依赖一些开源和免费的第 方软件包 包括(用于 文件的读取与存储)(用于科学计算)和(可视化)本软件在./.以开源形式发布 程序结构与流程程序主要包括数据读取、数据合并和光谱抽取 部分 其中数据读取包括扫描生成观测日志、读取各个步骤数据 数据合并包括本底合并、平场合并与改正以及定标灯红蓝端合并光谱抽取包括本底改正、平场改正、目标星定位、图像畸变改正、背景改正、波长定标和一维谱抽取等步骤 具体流程如图 用户只需要输入原始数据文件所在的路径 程序自动执行必要的步骤并存储各步骤产生的中间结果以便后续使用 对关键步骤进行可视化并得到最终数据产品其中本底()是指没有光子输入时 像素在通电时(零秒曝)的读出值 平场()是指输入光源为均匀光的情况下 获得的图像 反映了 的各个像素对光子的响应差异 定标灯是指使用空心阴极射线灯的光谱包含波长已知的发射谱线 用于标定 图像上波长与像素之间的对应关系 科学目标是光谱图 程序流程图 仪的使用者拍摄的天体目标 例如彗星、恒星和星系等 程序使用本文介绍的程序是针对 的长缝光谱设计 用于将光谱仪获取的原始观测数据(通常为 格式的光谱图像)转换为可供科学分析的一维光谱 即目标天体的波长与流量的对应关系 最终输出数据以 格式保存 用户需要在操作系统的终端运行本程序 并输入原始数据文件和观测日志所在的路径 程序自动执行必要的步骤 存储中间结果并进行可视化 得到最终数据产品 即目标天体的一维谱 这些结果和中间过程文件都存储在指定的计算机目录中 此外 用户也可以更改运行模式 手动对目标星进行标记 程序会根据标记位置寻找目标星准确轮廓 实际处理问题与解决方法 扫描生成观测日志 望远镜的观测者在观测同时经常会记录电子版的观测日志 但是为了避免格式不一致 期宋得阳等:长缝光谱自动处理软件 程序会扫描原始数据所在的计算机目录 读取其中 格式的文件 根据每个文件的头部()信息 结合用户记录的观测日志 重新生成可以机读的格式化观测日志 包含文件名、图像类型、观测目标源的名称、曝光时间、观测开始时间等信息 本底改正通常情况下 观测者在观测开始前或结束后在 相机快门关闭的情况下进行多次零秒曝光获得一组本底图像 通过图像合并获得本底的数值并抑制 读出噪声带来的误差 这个步骤提供了两种不同的合并模式 即均值模式和中值模式 通常情况下默认为均值模式 本底合并后生成本底图像 并将合并后的本底图像以 格式单独存储为一个文件 之后对需要操作的数据直接扣除保存的本底合并数据即完成本底改正 目标源定位本文介绍的程序采用以下方法确定目标天体的光谱在 图像上的位置:从光谱图像的左侧开始 沿 轴方向每隔 个像素取图像的截面最高处作为初始位置 取上下各 像素范围内的数据用高斯函数进行拟合 获取纵向(方向)的中心位置与半高全宽 并将中心位置当作新的初始位置重复上述过程 直到拟合参数收敛 接下来对所有 方向的中心位置与横坐标(值)的关系进行三阶多项式函数拟合 作为目标源光谱的中心位置曲线 光谱图像如图 此外程序还为用户提供了手动选取目标源的模式 在该模式下用户可以在弹出的图形界面上用鼠标标注目标源的大致位置 程序在其标记的位置附近按上述流程寻找目标天体的光谱图 ()中心位置与半高全宽图()平均轮廓图()真实目标源位置图()轮廓中心拟合图 ()()().()()天 文研究与技术 卷 平场改正由于 的各个像素对光子的响应略有差异 即便输入光源为均匀光 获得的图像也是不均匀的 在拍摄过程中一般进行多次曝光 将获得的多幅平场图像进行叠加 以获得较高的信噪比并抑制读出噪声 然后使用三次样条插值补齐数据中的坏像元数据 之后使用 滤波过程获得主平场图像 每一幅目标天体的光谱图像需要除以主平场图像进行平场改正 视场光学畸变改正由于 相机的光学畸变 图像上狭缝所成的像存在一定程度的弯曲 在波长定标灯的光谱图像上选取沿着 轴中间位置的一维谱作为基准光谱 在 轴的上下两个方向 每间隔 个像素进行一次同样的一维谱提取操作 将得到的光谱与基准光谱计算交叉相关函数 ()()()其中 和 分别为待改正的光谱与基准光谱 为光谱取点个数 为相对像素差值 和 分别为两光谱的标准差 交叉相关函数最大值对应的像素之差即为实际的像素偏移量 根据各行的像素差获得改正函数 并生成改正图像 如图 不同颜色的线代表 方向不同位置的光谱 经过视场畸变改正后能有效提高一维光谱的波长定标精度 避免系统误差图()畸变改正后的 定标灯谱的一维谱图像()畸变改正前的一维谱图像()()背景改正实际观测过程中 我们获得的光谱不仅来自于目标源(图()也有天光背景以及仪器内部的 期宋得阳等:长缝光谱自动处理软件 杂散光 目标星的光谱也包含这些成分 因此为了获取目标源的真实流量 我们需要对图像数据进行背景改正程序一般会控制目标星使其落在沿 轴中间附近的位置 程序自动选取 图像上距离目标星足够远、像素一致性比较好的行(例如图()中第 行以及 行)作为测定背景光大小的区域 在选取过程中 程序会计算该区域每一行沿着 轴(即色散方向)的平均值和标准差 对于强度超过均值以上 倍标准差的行 可能包含目标源附近其他天体的光谱 这种情况在密集星场中尤为常见 程序在计算背景值时对这些行进行了屏蔽 以免获得比真实值显著偏高的背景流量 目标源光谱附近的像素值减去背景光谱 即为目标天体的光谱图()目标源天体的原始光谱图像()本例中测定天空背景采用的像素区域()()波长定标波长定标的目的是在 图像上建立像素位置与波长的映射关系 本文介绍的程序包含自动对波长进行标定 具体处理流程如下:()读取该模式下提供的已经定标完成的模板光谱数据、模板特征发射线波长与像素位置对应数据 其次对要处理的光谱数据与模板光谱数据进行交叉相关函数计算获得两个光谱数据位置的偏移量 将读入的模板像素位置数据增加偏移量获得新的像素位置数据()依次取新的像素位置左右各 个像素的一维光谱 用广义高斯函数()/进行拟合 获得峰值位置对应的像素数值 结果如图 和图 其中每幅图中的圆点代表一维光谱数据 实线代表拟合曲线 虚线代表拟合得到的谱线中心位置()用五阶多项式函数进行拟合 获得波长与像素位置的对应关系函数 并保存全部波长与像素位置数据组 画出波长定标函数图像 如图 此外 由于 采用的定标灯本身的特性 其红端存在大量明亮的 元素的发射线 而蓝端以 线为主 它们的强度差异较大 如果曝光时间过长会导致红端饱和 而曝光时间过短会导致蓝端谱线没有足够的强度 因此 我们通常拍摄一次长曝光和一次短曝光 分别作为蓝端和红端的波长定标谱 在本文的例子中 红蓝两端的波长定标图像曝光时间分别为 和 因此 在进行定标之前本方法对两组光谱数据进行截取之后合并处理生成一个新的定标光谱图像 处理过程中计算交叉相关函数时对红端和蓝端分别求解计算天 文研究与技术 卷图 用于波长证认的 谱线峰值位置拟合 期宋得阳等:长缝光谱自动处理软件 图 用于波长证认的 谱线峰值位置拟合 图()多项式拟合得到的波长()与像素()的关系 蓝色点代表定标谱线()拟合残差().()优化抽谱优化抽谱方法采用加权求和 保证具有足够流量的情况下包含尽量少的噪声像素 对暗弱源尤其有效 能显著提高信噪比 优化抽谱过程主要分为两步:()获得目标星的精确空间轮廓 确定每个像素的加权值()根据空间轮廓的形状 排除明显高于轮廓的像素点 即扣除宇宙线 在我们开发的光谱数据处理程序中 第 步首先默认空间轮廓宽度为 个像素 轮廓方向每个像素具有相同的权重值 然后计算光谱数据方差 根据方差改正加权值与空间轮廓 重复此过程迭代 次 最终获得精确的空间轮廓 第 步 先根据上一步得到的轮廓提取科学目标的一维谱 并与普通抽谱方法得到的一维谱计算差值 对其中大于 倍标准偏差的像素点进行屏蔽 我们认为这些像素是受到宇宙线照射的像素 重复上述步骤直至收敛 至此我们获得了屏蔽宇宙线后的科学目标一维谱天 文研究与技术 卷 结果比较为了验证程序的可靠性 我们使用 对超新星 的同一组光谱数据进行了处理并绘制了一维谱图像 如图 其中 黄色线是本软件自动化处理获得的结果 蓝色线是 处理结果由图 可以看出 两种处理方法结果基本一致 并且本软件有效自动屏蔽了宇宙线 在流量比较低的蓝端 信噪比明显优于 的处理结果 需要在处理过程中对每条光谱手动屏蔽宇宙线 相比之下本文方法节省了处理时间且精度更高图 结果比较 蓝色曲线为 处理结果 黄色曲线为本软件处理结果 .总 结本文系统介绍了一款针对 望远镜 光谱仪的长缝光谱数据处理软件 主要流程包括生成日志、本底和平场的合并与改正、目标星定位、视场畸变改正、背景改正、波长定标、优化抽谱 本文的数据处理软件使用方便 相比 需要手动处理 本软件只需要用户输入原始数据所在的路径即可 处理过程透明可视 每个步骤完成都会生成相应的图片并保存数据以便重复使用 工作效率高 能批量化自动处理光谱数据且自动扣除宇宙线数据 与 处理结果相比效果更好 具有较高的实用价值致谢:感谢兴隆 望远镜全体工作人员的支持参考文献:.:.():.黄磷 吴宏 李宏斌.使用说明/.().:/././.().:/./.兴隆观测站 中国科学院国家天文台.终端介绍/.()期宋得阳等:长缝光谱自动处理软件 .:/././.().:/./.郑捷 江林巧.光变曲线自动抽取程序.天文研究与技术 ():.():.():./.张艳 林众 郭迪福.软件中 包的应用.天文研究与技术 ():.():.():.():.(.:.):.().: