def test_xy2sky_nord3(self):
x = 15000
y = 20000
crpix1 = -7535.57493517
crpix2 = 9808.40914361
crval1 = 176.486157083
crval2 = 8.03697351091
cd = [[5.115244026718E-05, 7.064503033578E-07],
[-1.280229655229E-07, -5.123112374523E-05]]
pv = [[
-7.030338745606E-03, 1.01755337222, 8.262429361142E-03,
0.00000000000, -5.910145454849E-04, -7.494178330178E-04,
-3.470178516657E-04, -2.331150605755E-02, -8.187062772669E-06,
-2.325429510806E-02, 1.135299506292E-04
],
[
-6.146513090656E-03, 1.01552885426, 8.259666421752E-03,
0.00000000000, -4.567030382243E-04, -6.978676921999E-04,
-3.732572951216E-04, -2.332572754467E-02,
-2.354317291723E-05, -2.329623852891E-02, 1.196394469003E-04
]]
nord = 3
ra, dec = wcs.xy2sky(x, y, crpix1, crpix2, crval1, crval2, cd, pv,
nord)
assert_that(ra, almost_equal(177.62041959006154, SIGFIGS))
assert_that(dec, almost_equal(7.5256066570082263, SIGFIGS))
def _update_ra_dec(self):
fits_header = self.fits_header
self._ra, self._dec = wcs.xy2sky(self.pixel_x, self.pixel_y,
float(fits_header[astrom.CRPIX1]),
float(fits_header[astrom.CRPIX2]),
float(fits_header[astrom.CRVAL1]),
float(fits_header[astrom.CRVAL2]),
wcs.parse_cd(fits_header),
wcs.parse_pv(fits_header),
wcs.parse_order_fit(fits_header))
def _update_ra_dec(self):
fits_header = self.fits_header
self._ra, self._dec = wcs.xy2sky(self.pixel_x, self.pixel_y,
float(fits_header['CRPIX1']),
float(fits_header['CRPIX2']),
float(fits_header['CRVAL1']),
float(fits_header['CRVAL2']),
wcs.parse_cd(fits_header),
wcs.parse_pv(fits_header),
wcs.parse_order_fit(fits_header))
def test_xy2sky_nord3(self):
x = 15000
y = 20000
crpix1 = -7535.57493517
crpix2 = 9808.40914361
crval1 = 176.486157083
crval2 = 8.03697351091
cd = [[5.115244026718E-05, 7.064503033578E-07],
[-1.280229655229E-07, -5.123112374523E-05]]
pv = [[-7.030338745606E-03, 1.01755337222, 8.262429361142E-03,
0.00000000000, -5.910145454849E-04, -7.494178330178E-04,
-3.470178516657E-04, -2.331150605755E-02, -8.187062772669E-06,
-2.325429510806E-02, 1.135299506292E-04],
[-6.146513090656E-03, 1.01552885426, 8.259666421752E-03,
0.00000000000, -4.567030382243E-04, -6.978676921999E-04,
-3.732572951216E-04, -2.332572754467E-02, -2.354317291723E-05,
-2.329623852891E-02, 1.196394469003E-04]]
nord = 3
ra, dec = wcs.xy2sky(x, y, crpix1, crpix2, crval1, crval2, cd, pv, nord)
assert_that(ra, almost_equal(177.62041959006154, SIGFIGS))
assert_that(dec, almost_equal(7.5256066570082263, SIGFIGS))
