def test_get_new_observer(aia171_test_map):
initial_obstime = aia171_test_map.date
rotation_interval = 2 * u.day
new_time = initial_obstime + rotation_interval
time_delta = new_time - initial_obstime
observer = get_earth(initial_obstime + rotation_interval)
# The observer time is set along with other definitions of time
for time in (rotation_interval, new_time, time_delta):
with pytest.raises(ValueError):
new_observer = _get_new_observer(initial_obstime, observer, time)
# Obstime property is present but the value is None
observer_obstime_is_none = SkyCoord(12*u.deg, 46*u.deg, frame=frames.HeliographicStonyhurst)
with pytest.raises(ValueError):
new_observer = _get_new_observer(None, observer_obstime_is_none, None)
# When the observer is set, it gets passed back out
new_observer = _get_new_observer(initial_obstime, observer, None)
assert isinstance(new_observer, SkyCoord)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value,
observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value, decimal=3)
# When the time is set, a coordinate for Earth comes back out
for time in (rotation_interval, new_time, time_delta):
with pytest.warns(UserWarning, match="Using 'time' assumes an Earth-based observer"):
new_observer = _get_new_observer(initial_obstime, None, time)
assert isinstance(new_observer, SkyCoord)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value,
observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value, decimal=3)
# The observer and the time cannot both be None
with pytest.raises(ValueError):
new_observer = _get_new_observer(initial_obstime, None, None)
def test_get_new_observer(aia171_test_map):
initial_obstime = aia171_test_map.date
rotation_interval = 2 * u.day
new_time = initial_obstime + rotation_interval
time_delta = new_time - initial_obstime
observer = get_earth(initial_obstime + rotation_interval)
# The observer time is set along with other definitions of time
for time in (rotation_interval, new_time, time_delta):
with pytest.raises(ValueError):
new_observer = _get_new_observer(initial_obstime, observer, time)
# Obstime property is present but the value is None
observer_obstime_is_none = SkyCoord(12*u.deg, 46*u.deg, frame=frames.HeliographicStonyhurst)
with pytest.raises(ValueError):
new_observer = _get_new_observer(None, observer_obstime_is_none, None)
# When the observer is set, it gets passed back out
new_observer = _get_new_observer(initial_obstime, observer, None)
assert isinstance(new_observer, SkyCoord)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value,
observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value, decimal=3)
# When the time is set, a coordinate for Earth comes back out
for time in (rotation_interval, new_time, time_delta):
new_observer = _get_new_observer(initial_obstime, None, time)
assert isinstance(new_observer, SkyCoord)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lon.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value,
observer.transform_to(frames.HeliographicStonyhurst).lat.to(u.deg).value, decimal=3)
np.testing.assert_almost_equal(new_observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value,
observer.transform_to(frames.HeliographicStonyhurst).radius.to(u.au).value, decimal=3)
# The observer and the time cannot both be None
with pytest.raises(ValueError):
new_observer = _get_new_observer(initial_obstime, None, None)
