Example #1
0
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
Example #2
0
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
Example #3
0
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
Example #4
0
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