animate_night.py 文件源码

def RaDec2region(ra, dec, nside):
    SCP_indx, NES_indx, GP_indx, WFD_indx = mutually_exclusive_regions(nside)

    indices = _raDec2Hpid(nside, np.radians(ra), np.radians(dec))
    result = np.empty(np.size(indices), dtype = object)
    SCP = np.in1d(indices, SCP_indx)
    NES = np.in1d(indices,NES_indx)
    GP  = np.in1d(indices,GP_indx)
    WFD = np.in1d(indices,WFD_indx)

    result[SCP] = 'SCP'
    result[NES] = 'NES'
    result[GP]  = 'GP'
    result[WFD] = 'WFD'

    return result