def gal2hor(l, b, time, date_obs=QubicSampling.DEFAULT_DATE_OBS,
latitude=DOMECLAT, longitude=DOMECLON):
"""
gal2hor(l, b, time, [date_obs, [latitude, [longitude]]]) -> az, el
Galactic to horizontal spherical conversion. Angles are in degrees.
Parameters
----------
time : array-like
Elapsed time in seconds since date_obs.
date_obs : string
The starting date, UTC.
latitude : float
The observer's latitude geolocation. Default is Dome C.
longitude : float
The observer's longitude geolocation. Default is Dome C.
Example
-------
>>> gal2hor(0, 0, 0)
(array(50.35837815921487), array(39.212362279976155))
"""
incoords, time = _format_sphconv(l, b, date_obs, time)
g2e = SphericalGalactic2EquatorialOperator(degrees=True)
e2h = SphericalEquatorial2HorizontalOperator(
'NE', time, latitude, longitude, degrees=True)
outcoords = e2h(g2e(incoords))
return outcoords[..., 0], outcoords[..., 1]
def gal2equ(l, b):
"""
gal2equ(l, b) -> ra, dec
Galactic to equatorial spherical conversion. Angles are in degrees.
"""
incoords = _format_sphconv(l, b)
outcoords = SphericalGalactic2EquatorialOperator(degrees=True)(incoords)
return outcoords[..., 0], outcoords[..., 1]