def get_heliographic_coordinates(self):
if self.heliographic_coordinates is None:
hcc_x, hcc_y, hcc_z = self.get_heliocentric_coordinates()
self.heliographic_coordinates = wcs.convert_hcc_hg(
hcc_x, hcc_y, hcc_z)[0:2]
return self.heliographic_coordinates
def test_conv_hcc_hg():
coord = [13.0, 58.0]
result = wcs.convert_hcc_hg(img.rsun_meters,
img.heliographic_latitude, img.heliographic_longitude,
coord[0], coord[1])
known_answer = [1.0791282e-06, -7.0640732]
magnitude = np.floor(np.log10(np.abs(known_answer)))
assert_array_almost_equal(result*10**(-magnitude),
known_answer*10**(-magnitude), decimal=2)
def test_conv_hcc_hg():
coord = [13.0, 58.0]
result = wcs.convert_hcc_hg(img.rsun_meters, img.heliographic_latitude,
img.heliographic_longitude, coord[0], coord[1])
known_answer = [1.0791282e-06, -7.0640732]
magnitude = np.floor(np.log10(np.abs(known_answer)))
assert_array_almost_equal(result * 10**(-magnitude),
known_answer * 10**(-magnitude),
decimal=2)
def test_conv_hcc_hg():
coord = [13.0, 58.0]
result = wcs.convert_hcc_hg(img.rsun_arcseconds,
img.heliographic_latitude, img.heliographic_longitude,
coord[0], coord[1])
assert_array_almost_equal(result, [1.0791282e-06, -7.0640732], decimal=2)
def test_conv_hcc_hg():
coord = [13.0, 58.0]
result = wcs.convert_hcc_hg(header, coord[0], coord[1])
assert_array_almost_equal(result, [1.0791282e-06, -7.0640732], decimal=2)