1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
#!/usr/bin/env python
from argparse import ArgumentParser
from pymap3d.haversine import anglesep, anglesep_meeus
from pytest import approx
p = ArgumentParser(description="angular distance between two sky points")
p.add_argument("r0", help="right ascension: first point [deg]", type=float)
p.add_argument("d0", help="declination: first point [deg]", type=float)
p.add_argument("r1", help="right ascension: 2nd point [deg]", type=float)
p.add_argument("d1", help="declination: 2nd point [degrees]", type=float)
a = p.parse_args()
dist_deg = anglesep_meeus(a.r0, a.d0, a.r1, a.d1)
dist_deg_astropy = anglesep(a.r0, a.d0, a.r1, a.d1)
print(f"{dist_deg:.6f} deg sep")
assert dist_deg == approx(dist_deg_astropy)
|
