1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#!/usr/bin/env python
"""
Example Kitt Peak
./demo_azel2radec.py 264.9183 37.911388 31.9583 -111.597 2014-12-25T22:00:00MST
"""
from argparse import ArgumentParser
from pymap3d import azel2radec
p = ArgumentParser(
description="convert azimuth and elevation to " "right ascension and declination"
)
p.add_argument("azimuth", help="azimuth [deg]", type=float)
p.add_argument("elevation", help="elevation [deg]", type=float)
p.add_argument("lat", help="WGS84 obs. lat [deg]", type=float)
p.add_argument("lon", help="WGS84 obs. lon [deg]", type=float)
p.add_argument("time", help="obs. time YYYY-mm-ddTHH:MM:SSZ")
P = p.parse_args()
ra, dec = azel2radec(P.azimuth, P.elevation, P.lat, P.lon, P.time)
print("ra [deg] ", ra, " dec [deg] ", dec)
|
