天体源星表的概率位置关联:Aspects 代码
摘要
我们描述了一种基于位置及其不确定度对两个星表中的天体源进行交叉证认的概率方法。在两种互斥假设下推导了某天体与另一星表中源相关联或其无对应体的概率:首先是多对一关联的经典情形,其次是更为现实但更困难的一对一关联问题。在这两种情形下,均计算了在有效位置观测到两个星表中天体的似然度,并构建了具有对应体的源的比例(计算关联概率所需的量)的最大似然估计量。当一个或两个星表的位置不确定度未知时,该方法可用于估计其典型值,甚至研究其对天体大小的依赖性。该方法也可应用于源的真实中心与其在另一波长的对应体中心不重合的情形。为了在不同假设下计算似然度和关联概率,我们开发了一个名为"Aspects"([asp{\epsilon}],法语"ASsociation PositionnellE/ProbabilistE de CaTalogues de Sources"的缩写)的 Fortran 95 代码;其源文件已公开免费提供。为测试 Aspects,我们创建了包含多达 个天体的全天模拟星表,强制设定为多对一或一对一关联。对这些模拟的分析证实,在这两种情况下,具有最高似然度的假设即为正确假设,且针对适当关联模型构建的未知参数估计量是可靠的。
引用
@article{arxiv.1209.5361,
title = {Probabilistic positional association of catalogs of astrophysical sources: the Aspects code},
author = {Michel Fioc},
journal= {arXiv preprint arXiv:1209.5361},
year = {2014}
}
备注
One typo corrected and links added. One-column format (A&A official version is in two-column format). 20 pages, 5 figures. Version 2 of code Aspects available at http://www2.iap.fr/users/fioc/Aspects/ . Detailed documentation and complements at arXiv:1404.4224. Numerical Recipes routines not needed anymore