@@ -0,0 +1,539 @@
+from astropy.coordinates import SkyCoord%0Aimport astropy.units as u%0A%0A%0Adef angular_distance(ra1, dec1, ra2, dec2):%0A    %22%22%22%0A    Compute angular distance between pairs of coordinates%0A%0A    :param ra1:%0A    :param dec1:%0A    :param ra2:%0A    :param dec2:%0A    :return: angular distance%0A    %22%22%22%0A    point_1 = SkyCoord(ra=ra1 * u.degree, dec=dec1 * u.degree, frame='icrs')%0A%0A    point_2 = SkyCoord(ra=ra2 * u.degree, dec=dec2 * u.degree, frame='icrs')%0A%0A    angular_distance = point_1.separation(point_2).to(u.degree)%0A%0A    return angular_distance.value%0A

Compute angular distances between pairs of coordinates