def hor2cel(coord, time, site, copy=True): coord = np.array(coord, copy=copy) trepr = time[len(time)/2] info = iers.lookup(trepr) ao = slalib.sla_aoppa(trepr, info.dUT, site.lon*utils.degree, site.lat*utils.degree, site.alt, info.pmx*utils.arcsec, info.pmy*utils.arcsec, site.T, site.P, site.hum, 299792.458/site.freq, site.lapse) am = slalib.sla_mappa(2000.0, trepr) # This involves a transpose operation, which is not optimal pyfsla.aomulti(time, coord.T, ao, am) return coord
def hor2cel(coord, time, site, copy=True): coord = np.array(coord, copy=copy) trepr = time[len(time) / 2] info = iers.lookup(trepr) ao = pyfsla.sla_aoppa(trepr, info.dUT, site.lon * utils.degree, site.lat * utils.degree, site.alt, info.pmx * utils.arcsec, info.pmy * utils.arcsec, site.T, site.P, site.hum, 299792.458 / site.freq, site.lapse) am = pyfsla.sla_mappa(2000.0, trepr) # This involves a transpose operation, which is not optimal pyfsla.aomulti(time, coord.T, ao, am) return coord
def cel2hor(coord, time, site, copy=True): from . import pyfsla from enlib import iers # This is very slow for objects near the horizon! coord = np.array(coord, copy=copy) trepr = time[len(time)/2] info = iers.lookup(trepr) ao = pyfsla.sla_aoppa(trepr, info.dUT, site.lon*utils.degree, site.lat*utils.degree, site.alt, info.pmx*utils.arcsec, info.pmy*utils.arcsec, site.T, site.P, site.hum, 299792.458/site.freq, site.lapse) am = pyfsla.sla_mappa(2000.0, trepr) # This involves a transpose operation, which is not optimal pyfsla.oamulti(time, coord.T, ao, am) return coord
def cel2hor(coord, time, site, copy=True): from . import pyfsla from enlib import iers # This is very slow for objects near the horizon! coord = np.array(coord, copy=copy) trepr = time[len(time) / 2] info = iers.lookup(trepr) ao = pyfsla.sla_aoppa(trepr, info.dUT, site.lon * utils.degree, site.lat * utils.degree, site.alt, info.pmx * utils.arcsec, info.pmy * utils.arcsec, site.T, site.P, site.hum, 299792.458 / site.freq, site.lapse) am = pyfsla.sla_mappa(2000.0, trepr) # This involves a transpose operation, which is not optimal pyfsla.oamulti(time, coord.T, ao, am) return coord