sphDist

lsst.validate.drp.util.sphDist(ra_mean, dec_mean, ra, dec)

Calculate distance on the surface of a unit sphere.

Parameters:

ra_mean : float

Mean RA in radians.

dec_mean : float

Mean Dec in radians.

ra : numpy.array [float]

Array of RA in radians.

dec : numpy.array [float]

Array of Dec in radians.

Notes

Uses the Haversine formula to preserve accuracy at small angles.

Law of cosines approach doesn’t work well for the typically very small differences that we’re looking at here.