matching.py 文件源码

python
阅读 34 收藏 0 点赞 0 评论 0

项目:astronomy-utilities 作者: astronomeralex 项目源码 文件源码
def matching(objra,objdec,ras,decs,boxsize=1):
    """
    objra and objdec are in degrees
    ras and decs are numpy arrays of all the catalog values (also in degrees)
    """
    obj = coord.ICRSCoordinates(ra=objra, dec=objdec, unit=(u.degree, u.degree))
    distances = []
    for i in range(len(ras)):
        if abs(ras[i] - objra) > 0.1 or abs(decs[i] - objdec) > 0.1:
            distances.append(9999)
        else:
            catobj = coord.ICRSCoordinates(ra=ras[i], dec=decs[i], unit=(u.degree, u.degree))
            distances.append(obj.separation(catobj).arcsecs)
    distances = np.array(distances)
    minindex = np.argmin(distances)
    mindist = np.amin(distances)
    b = distances < boxsize
    boxnum = len(distances[b])
    return [minindex,mindist,boxnum]
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号