Coverage for tests/testApprox.py : 31%

Hot-keys on this page
r m x p toggle line displays
j k next/prev highlighted chunk
0 (zero) top of page
1 (one) first highlighted chunk
lsst.utils.tests.init()
""" Test the fast approximate ra,dec to alt,az transforms """
nside = 64 hpids = np.arange(hp.nside2npix(nside)) ra, dec = utils.hpid2RaDec(nside, hpids) mjd = 59852. obs = utils.ObservationMetaData(mjd=mjd)
alt1, az1, pa1 = utils.altAzPaFromRaDec(ra, dec, obs)
alt2, az2 = utils.approx_RaDec2AltAz(ra, dec, obs.site.latitude, obs.site.longitude, mjd)
# Check that the fast is similar to the more precice transform tol = 2 # Degrees tol_mean = 1. separations = utils.angularSeparation(az1, alt1, az2, alt2) self.assertLess(np.max(separations), tol) self.assertLess(np.mean(separations), tol_mean)
# Check that the fast can nearly round-trip ra_back, dec_back = utils.approx_altAz2RaDec(alt2, az2, obs.site.latitude, obs.site.longitude, mjd) separations = utils.angularSeparation(ra, dec, ra_back, dec_back) self.assertLess(np.max(separations), tol) self.assertLess(np.mean(separations), tol_mean)
nside = 64 hpids = np.arange(hp.nside2npix(nside)) ra, dec = utils._hpid2RaDec(nside, hpids) mjd = 59852. obs = utils.ObservationMetaData(mjd=mjd)
alt1, az1, pa1 = utils._altAzPaFromRaDec(ra, dec, obs)
alt2, az2 = utils._approx_RaDec2AltAz(ra, dec, obs.site.latitude_rad, obs.site.longitude_rad, mjd)
# Check that the fast is similar to the more precice transform tol = np.radians(2) tol_mean = np.radians(1.) separations = utils._angularSeparation(az1, alt1, az2, alt2)
self.assertLess(np.max(separations), tol) self.assertLess(np.mean(separations), tol_mean)
# Check that the fast can nearly round-trip ra_back, dec_back = utils._approx_altAz2RaDec(alt2, az2, obs.site.latitude_rad, obs.site.longitude_rad, mjd) separations = utils._angularSeparation(ra, dec, ra_back, dec_back) self.assertLess(np.max(separations), tol) self.assertLess(np.mean(separations), tol_mean)
lsst.utils.tests.init() unittest.main() |