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))
