云南天文台2.4 m望远镜几何扭曲变化初步探究.pdf

.天 文 研 究 与 技 术 第 卷 第 期 年 月:./.云南天文台 望远镜几何扭曲变化初步探究陈 筱 郭碧峰 郑中杰 彭青玉(.暨南大学计算机科学系 广东 广州 .暨南大学中法天体测量、动力学与空间科学联合实验室广东 广州 .广东海洋大学数学与计算机学院 广东 湛江)摘要:几何扭曲()是电荷耦合器件()高精度测量不可忽略的因素 通常需要细致的观测工作导出 这必然占用宝贵的观测时间利用丽江观测站 望远镜观测密集星场的 图像 初步发现了图像几何扭曲随天顶距()的变化规律 具体地 本文参考 星表提供的位置信息并进行数据归算 改正了较差颜色折射()效应 求解出不同观测时间 的几何扭曲模型 并使用四阶多项式拟合 结果表明 图像的几何扭曲与观测时的天顶距存在一定的相关性 在像素 方向 几何扭曲模型四次项 和、二次项 和 与天顶距存在较明显的线性关系 在 方向 与天顶距存在线性关系的系数项为四次项 和、三次项 和二次项 关键词:几何扭曲 天体测量 图像处理 天体测量定标中图分类号:文献标识码:文章编号:()望远镜光学成像系统的几何扭曲会导致观测的 图像中目标星像偏离原本的像素坐标位置这是一种普遍存在的现象 该现象不仅存在于空间望远镜 也存在于地面望远镜 对视场中星像位置测量产生不可忽略的影响 例如 在哈勃空间望远镜()的光学成像系统中 视场存在严重的几何扭曲 其中第 代广角行星照相机()的几何扭曲最大可达约 为了准确求解哈勃空间望远镜的几何扭曲效应 文提出了通过提高线性项精度来改正几何扭曲的方法 使得行星相机视场精确到 宽场()视场精确到 在地面望远镜的观测中 年 文将哈勃空间望远镜的几何扭曲求解方法应用于欧洲南方天文台 望远镜 取得了约 的位置测量精度 随后 文提出了一种基于参考星表位置测量信息的几何扭曲求解方法 提高了地面望远镜的观测准确度与精度 并成功应用于土卫九()的位置测量 此后 该方法应用于多项研究 文在校正几何扭曲后得到近地小行星 的位置信息 在 星表发布后 该方法成功应用于木卫六()与海卫一()的位置测量 文参考 星表 并从准确度、精度等方面系统比较了文与文的方法 随着最新的 参考星表发布 文所述方法的应用将更广泛 能够更精确地测量天体位置文和文这两种方法在求解几何扭曲时对所使用的观测资料有严格要求 观测资料需要采用不同取向或指向的观测模式 均匀抖动拍摄 这一工作通常需要充足的观测时间和一定的观测经验才能完成 然而 望远镜的观测机会宝贵 观测时间有限 为了能充分利用有限、宝贵的观测时间本文希望从天顶距方面探究云南天文台丽江观测站 望远镜 图像的几何扭曲变化规律 天基金项目:科技部国家重点研发计划政府间国际科技创新合作项目()中国载人航天工程巡天空间望远镜专项科学研究项目()国家自然科学基金()资助.收稿日期:修订日期:作者简介:陈 筱 女 硕士研究生 主要研究图像处理与精确测量.:.通信作者:彭青玉 男 教授 主要研究图形图像与精确测量.:.天 文研究与技术 卷顶距是天体与天顶之间的角距离 天顶距越大 大气折射越大 然而 大气折射造成天体的视方向与真实方向产生偏差 影响天体位置测量 并且由大气折射造成的扭曲项是光学扭曲的重要项 文求解了大气折射扭曲项对暗能量相机天体测量定标的影响 因此 本文认为天顶距可能在一定程度上影响几何扭曲模型的变化 本文采用文提出的方法求解几何扭曲模型 参考最新发布的 星表进行数据归算 并考虑较差颜色折射效应的影响 计算出不同时间点观测的几何扭曲模型 试图发现几何扭曲模型随观测天顶距变化的规律 本文希望通过该规律 能够为后续的几何扭曲求解提供一定的便利 即通过较少次数定标场的观测计算 图像在其他时间的几何扭曲 从而有效地节约望远镜观测定标场的时间 降低望远镜观测定标的频率 观测资料本文采用云南天文台丽江观测站 光学望远镜在 年 月 日 日对疏散星团 与近地小行星 观测的资料探究几何扭曲模型的变化规律 观测所用望远镜和 的详细参数见表 本文选取 组不同天顶距抖动观测的密集星场的资料(视场中大约能观测到 颗恒星)其中包括 帧 星团图像(原用于视场定标)帧近地小行星 图像(原用于高精度位置测量)具体的观测资料如表 表 望远镜和 参数说明 /()()/(/)表 观测资料概要 /()/()/()/()/()/()/()/()/()/()/()/()方 法 视场定标对于每幅 图像 采用如下步骤进行天体测量归算 首先 对图像进行平场与本底校正 并裁剪图像的边缘区域 裁剪后的图像像素大小为 然后通过二维高斯定心算法测量每个星像的像素坐标位置 并将像素坐标位置与 星表提供的位置信息进行匹配 在星表匹配过程中 为避免双星系统导致的位置误差 本文剔除了()值大于 期陈 筱等:云南天文台 望远镜几何扭曲变化初步探究 的恒星 将匹配到的每颗恒星转换到观测时刻的站心视位置 并考虑大气折射等天体测量效应再通过中心投影公式得到每个星像的标准坐标 最后 使用加权模型求解四阶多项式底片常数模型得到每幅 图像的定标参数 根据求解的底片常数模型计算每个星像的观测位置 本文归算了不同时刻的观测资料 对于每组资料 计算得到了每颗恒星的残差()和标准差 较差颜色折射效应校正由于光在空气中的折射率与波长有关 不同波长的光在天顶方向受到不同的折射 使得不同恒星的光穿过地球大气时折射为不同的光谱 因此 具有不同光谱的天体成像位置受到不同程度的大气折射影响 这种现象称为较差颜色折射 该现象可能导致目标天体的位置测量结果包含沿天顶方向的系统误差 本文使用的观测资料均采用 滤光片拍摄 该滤光片受较差颜色折射影响较为严重 为了求出更准确的几何扭曲模型 本文首先通过文提出的较差颜色折射改正方法消除每个星像受较差颜色折射的影响 首先 参考 星表恒星的位置与测光信息项(色指数)匹配图像中的每个星像 根据星表中的天体测量数据 转换到地平坐标系中的视位置(包含大气折射)得到每颗测量星的地平高度和方位角 然后 通过加权的四阶多项式底片常数模型归算得到每颗恒星的残差 使用一阶多项式()()()拟合色指数和地平高度方向的残差之间的关系 其中 表示参考星的色指数 为观测资料的天顶距 和 是拟合参数 通过加权最小二乘法进行拟合 求解较差颜色折射改正系数 和 改正恒星在高度方向的残差 最后 对每个星像的像素坐标位置进行校正本文对 组观测资料均进行了较差颜色折射改正 图 给出了较差颜色折射改正前后资料集 在天顶距方向的残差随色指数的变化情况 通过()式用加权最小二乘法拟合 得到拟合参数 蓝色虚线为拟合结果 从图()可以看出 改正较差颜色折射影响后 系统误差得到了明显改善 表明较差颜色折射对天顶方向的影响基本改正图 ()较差颜色折射改正前资料集 的残差随色指数的变化 蓝色实线表示()式的拟合结果()较差颜色折射改正后的结果()()()几何扭曲校正关于几何扭曲求解的方法 文通过自定标方法直接将图像的像素坐标位置构建为自洽的无扭天 文研究与技术 卷曲参考框架求解几何扭曲 文基于天体测量星表并消除相当的星表系统误差的方法求解几何扭曲 对于每组观测资料 通过本文第 节的方法改正较差颜色折射影响后 通过如下步骤求解几何扭曲()使用加权四常数底片模型进行归算 归算步骤见本文第 节 计算每幅图像中所有匹配星像的残差值()()对于出现在()幅不同的 图像中的星像 令 和 分别表示其中两幅不同的图像 用 ()表示该星像在第 幅和第 幅图像中像素位置处的几何扭曲()和()之间的数学关系其中 /为四常数模型求解的估计参数 和 分别为 的取向和比例尺 和 分别为星像的测量像素位置()与通过使用估计参数的四常数模型间接计算得到的该星像的像素位置()的差值 为天球切平面切点在赤纬方向的赤道坐标 通过求解()式可以推导该星像像素位置的几何扭曲()把 图像划分为 的方形区域 将像素坐标位于同一区域的所有星像的几何扭曲求平均值 并把该平均值作为区域中心像素坐标位置处的几何扭曲()通过双线性插值计算所有星像在像素位置处的几何扭曲值 然后对每幅图像所有星像的像素坐标进行校正 得到初步改正几何扭曲的观测数据结果()再次执行上述求解过程 得到几何扭曲增量 与第 次几何扭曲模型相加 得到第 次几何扭曲模型()重复迭代以上求解过程 直到几何扭曲增量达到预定的精度要求 结果分析 几何扭曲模型变化规律本文对 组资料进行了几何扭曲求解 图 分别给出了每晚典型的几何扭曲矢量图(从左到右依次为资料集 )我们使用四次多项式 ()表示几何扭曲模型 其中 和 为输入位置 和 归一化后的位置 范围在 内()/()/()由归一化位置可以更容易看出每一项对扭曲值影响的大小 表 和表 分别给出了在 方向和 方向不同天顶距观测资料的几何扭曲模型多项式的四次项、三次项和二次项的系数 表的第 行依次为天顶距(以()为单位)、四次项系数、三次项系数和二次项系数图 和图 分别给出了在 方向和 方向几何扭曲模型的四次项、三次项和二次项系数随天顶距的变化情况 并在每个子图的上方列出了系数项、系数的均值与标准差 我们初步发现三次项系数的标准差比较小 系数比较稳定 而四次项和二次项的某些系数随着天顶距的变化具有一定的规律 因此 我们通过线性拟合分析各项系数与天顶距的关系 图 和图 中蓝色实线为线性拟合的结果 表 和表 分别列出了拟合结果的皮尔森()相关系数与标准差()的值 本文根据皮尔森相关性分析各项系数与天顶距的相关强度 皮尔森相关系数 为 期陈 筱等:云南天文台 望远镜几何扭曲变化初步探究 图 几何扭曲模型矢量图 每个子图的上方列出了观测日期以及以像素为单位的扭曲最大值与扭曲中值 每个子图的矢量值放大 倍 .()()表 在 方向几何扭曲模型的四阶多项式系数 /()表 在 方向几何扭曲模型的四阶多项式系数 /()天 文研究与技术 卷图 在 方向几何扭曲模型系数随天顶距的变化情况 蓝色实线为线性拟合得到的结果 期陈 筱等:云南天文台 望远镜几何扭曲变化初步探究 图 在 方向几何扭曲模型系数随天顶距的变化情况 蓝色实线为线性拟合得到的结果 天 文研究与技术 卷 ()其中 为样本量()和()分别为两个变量的观测值和均值 和 为样本标准差 相关系数绝对值大于 表示两个变量的相关程度为中等及以上 本文取相关系数绝对值在 以上的为关系显著 在 方向 关系显著的系数项为四次项 和、二次项 和 在 方向关系显著的系数项为四次项 和、三次项 和二次项 表 线性拟合结果的皮尔森相关系数 表 线性拟合结果的标准差 讨 论从系统误差角度考虑 几何扭曲现象与望远镜的重力形变和大气折射误差有关 地面望远镜受重力影响 望远镜支撑结构、光学元器件发生弯曲或变形 导致地面观测系统的扭曲解不稳定 望远镜结构的重力形变产生的误差通常随着支撑结构的仰俯角变化 光学元器件的形变一定程度上也会引起光学系统的误差 地面望远镜在观测目标天体时还受到大气折射效应影响 从而导致观测到的目标高度比实际高度要高 这种效应随天顶距的增大而逐渐明显 虽然我们根据观测时的平均温度和气压考虑了大气折射模型的改正 但是该模型的改正未必准确 例如湿度、风力、观测方向、测站地理位置、观测目标光谱型等因素都影响大气折射误差的校正 这些偏差导致几何扭曲模型高阶项存在几何扭曲模型的常数项表示望远镜坐标轴的偏移程度 通常由于望远镜弯曲随观测时间和位置发生变化 一次项分别为望远镜的比例尺与取向 实际比例尺的变化通常由望远镜的热变化导致 模型的取向项用于对齐测量坐标系与标准坐标系 与底片测量仪器的定标和高阶项的建模不准确有关因此 本文没有探究几何扭曲模型的常数项、一次项与天顶距的关系 二次项用于修正底片倾斜导致的误差 通常由 种效应引入:切平面坐标系之间的转换 从相对于光轴倾斜的坐标系到切平面坐标系的转换 以及光学场角扭曲()定心不准确 三次项主要用于修正光学场角扭曲 光学场角扭曲定心不准确会造成像素坐标位置偏移 光学场角扭曲通常表示为与坐标有关的三次项形式 其微分为二次多项式 在光学场角扭曲系数求解过程中产生的误差在转换到切平面坐标的过程中引入二次项与三次项 由于光学场角扭曲与光线相对于光轴的角度有关 我们猜测天顶距不同的天体 其光线通过望远镜光学系统的折射角度不同 所以部分二次项和三次项的系数也存在与天顶距的相关性 四次项可以修正大气折射未完全校正导致的系统偏移 天顶距小的天体受大气折射的影响小 天顶距大的天体受大气折射的影响大 由图 可以看出 在 方向 随着天顶距变大 四次项 和 的系数对扭曲的影响逐渐增加 在 方向 四次项 和 的系数变化同理 期陈 筱等:云南天文台 望远镜几何扭曲变化初步探究 几何扭曲模型高阶项在 方向与 方向的系数项呈中心对称 在 方向 四次项 和 的系数与三次项 和 的系数对扭曲的影响较大 在 方向 四次项 和 的系数与三次项 和 的系数对扭曲的影响较大 本文初步发现的规律显示 天顶距越大 对天体的位置测量影响越大 总结与展望本文探究了丽江观测站 望远镜 几何扭曲模型的变化规律 我们选取了 年 月 日 日对密集星场抖动观测的 组资料 分别求解几何扭曲模型 并用四阶多项式拟合每个星像的几何扭曲与像素坐标位置的关系 本文初步发现几何扭曲模型的四次项和二次项的一些系数与天顶距有一定的相关性 在 方向 四次项 和、二次项 和 与天顶距存在明显的线性关系 在 方向 四次项 和、三次项 和二次项 与天顶距存在明显的线性关系 本文希望通过此规律 为今后丽江观测站 望远镜观测定标场节约一定的时间 通过较少次数的定标场观测推导其他时间观测的 图像的几何扭曲 本文的工作仅初步从天顶距方面探究几何扭曲模型的变化规律 在后续工作中 我们将从温度、气压等方面继续深入研究 望远镜的几何扭曲变化规律致谢:感谢丽江 望远镜全体工作人员的支持 该望远镜资金由中国科学院和云南省政府提供参考文献:.():.():.():.傅夏青 彭青玉 冯杰.基于几何扭曲模型的 图像扭曲校正残差及分析.天文研究与技术 ():.():.():.():.():.():.天 文研究与技术 卷 .():.():.付丽 董健 凌权宝 等.天马望远镜结构重力变形对面形和指向精度影响.电波科学学报 ():.():.(.:.):.().().: