def test_manual_wcs_angle(): # Checked with DS9 x, y = (11, 11) ra, dec = (10.0, 0.0) ps = 36 # arcsec/px north = 45 wcs = wcs_from_coords(x, y, ra, dec, ps, north) assert_almost_equal(wcs.all_pix2world(11, 11, 1), (ra, dec)) assert_almost_equal(wcs.all_pix2world(21, 21, 1), (10.0, 0.14142)) assert_almost_equal(wcs.all_pix2world(1, 1, 1), (10.0, -0.14142)) assert_almost_equal(wcs.all_pix2world(1, 21, 1), (10.14142, 0.0)) assert_almost_equal(wcs.all_pix2world(21, 1, 1), (9.858579, 0.0))
def test_manual_wcs_top_flip_all(): # Checked with DS9 x, y = (11, 11) ra, dec = (10.0, 0.0) ps = 36 # arcsec/px north = 'top' flip = 'all' wcs = wcs_from_coords(x, y, ra, dec, ps, north, flip=flip) assert_almost_equal(wcs.all_pix2world(11, 21, 1), (10.0, -0.1)) assert_almost_equal(wcs.all_pix2world(21, 11, 1), (10.1, 0.0)) assert_almost_equal(wcs.all_pix2world(1, 11, 1), (9.9, 0.0)) assert_almost_equal(wcs.all_pix2world(11, 1, 1), (10.0, 0.1)) assert_almost_equal(wcs.all_pix2world(21, 21, 1), (10.1, -0.1)) assert_almost_equal(wcs.all_pix2world(1, 21, 1), (9.9, -0.1)) assert_almost_equal(wcs.all_pix2world(1, 1, 1), (9.9, 0.1))
def test_manual_wcs_right(): # Checked with DS9 x, y = (11, 11) ra, dec = (10.0, 0.0) ps = 36 # arcsec/px north = 'right' wcs = wcs_from_coords(x, y, ra, dec, ps, north) assert_almost_equal(wcs.all_pix2world(11, 11, 1), (ra, dec)) assert_almost_equal(wcs.all_pix2world(11, 21, 1), (10.1, 0.0)) assert_almost_equal(wcs.all_pix2world(21, 11, 1), (10.0, 0.1)) assert_almost_equal(wcs.all_pix2world(1, 11, 1), (10.0, -0.1)) assert_almost_equal(wcs.all_pix2world(11, 1, 1), (9.9, 0.)) assert_almost_equal(wcs.all_pix2world(21, 21, 1), (10.1, 0.1)) assert_almost_equal(wcs.all_pix2world(1, 21, 1), (10.1, -0.1)) assert_almost_equal(wcs.all_pix2world(1, 1, 1), (9.9, -0.1))
def test_raise_north_angle(): with pytest.raises(ValueError) as exc: wcs_from_coords(0, 0, 0, 0, 0, 'not a direction') assert_in('invalid value for north', str(exc.value))