def test_hgs_hcc_observer_mismatch(): # Test whether the transformation gives the same answer regardless of what obstime the observer # coordinate is represented in observer1 = HeliographicStonyhurst(0*u.deg, 0*u.deg, 1*u.AU, obstime='2001-01-01') observer2 = observer1.transform_to(HeliographicStonyhurst(obstime='2001-03-31')) hgs = HeliographicStonyhurst(20*u.deg, 40*u.deg, 0.5*u.AU, obstime=observer1.obstime) hcc1 = hgs.transform_to(Heliocentric(observer=observer1, obstime=hgs.obstime)) hcc2 = hgs.transform_to(Heliocentric(observer=observer2, obstime=hgs.obstime)) assert_quantity_allclose(hcc1.cartesian.xyz, hcc2.cartesian.xyz)
def test_hcc_hgs_observer_mismatch(): # Test whether the transformation gives the same answer regardless of what obstime the observer # coordinate is represented in observer1 = HeliographicStonyhurst(0 * u.deg, 0 * u.deg, 1 * u.AU, obstime='2001-01-01') observer2 = observer1.transform_to( HeliographicStonyhurst(obstime='2001-03-31')) hcc1 = Heliocentric(0.2 * u.AU, 0.3 * u.AU, 0.4 * u.AU, observer=observer1, obstime=observer1.obstime) hgs1 = hcc1.transform_to(HeliographicStonyhurst(obstime=hcc1.obstime)) hcc2 = Heliocentric(0.2 * u.AU, 0.3 * u.AU, 0.4 * u.AU, observer=observer2, obstime=observer1.obstime) hgs2 = hcc2.transform_to(HeliographicStonyhurst(obstime=hcc2.obstime)) assert_quantity_allclose(hgs1.lon, hgs2.lon) assert_quantity_allclose(hgs1.lat, hgs2.lat) assert_quantity_allclose(hgs1.radius, hgs2.radius)
def test_hgs_hgc_roundtrip(): obstime = "2011-01-01" hgsin = HeliographicStonyhurst(lat=10*u.deg, lon=20*u.deg, obstime=obstime) hgcout = hgsin.transform_to(HeliographicCarrington(obstime=obstime)) assert_quantity_allclose(hgsin.lat, hgcout.lat) assert_quantity_allclose(hgsin.lon + sun.L0(obstime), hgcout.lon) hgsout = hgcout.transform_to(HeliographicStonyhurst(obstime=obstime)) assert_quantity_allclose(hgsout.lat, hgsin.lat) assert_quantity_allclose(hgsout.lon, hgsin.lon)
def test_hgs_hgc_roundtrip(): obstime = "2011-01-01" hgsin = HeliographicStonyhurst(lat=0*u.deg, lon=0*u.deg, obstime=obstime) hgcout = hgsin.transform_to(HeliographicCarrington(obstime=obstime)) assert_quantity_allclose(hgsin.lat, hgcout.lat) assert_quantity_allclose(hgcout.lon, get_sun_L0(obstime)) hgsout = hgcout.transform_to(HeliographicStonyhurst(obstime=obstime)) assert_quantity_allclose(hgsout.lat, hgsin.lat) assert_quantity_allclose(hgsout.lon, hgsin.lon)