def _set_altaz_fields(self, leafsize=10):
"""
Have a fixed grid of alt,az pointings to use
"""
tmp = read_fields()
names = ['alt', 'az']
types = [float, float]
self.fields = np.zeros(tmp.size, dtype=list(zip(names, types)))
self.fields['alt'] = tmp['dec']
self.fields['az'] = tmp['RA']
x, y, z = treexyz(self.fields['az'], self.fields['alt'])
self.field_tree = kdtree(list(zip(x, y, z)), leafsize=leafsize,
balanced_tree=False, compact_nodes=False)
hpids = np.arange(hp.nside2npix(self.nside))
self.reward_ra, self.reward_dec = _hpid2RaDec(self.nside, hpids)
评论列表
文章目录