def test_hpc_low_precision_float_warning():
hpc = Helioprojective(u.Quantity(0, u.deg, dtype=np.float32),
u.Quantity(0, u.arcsec, dtype=np.float16),
observer=HeliographicStonyhurst(0*u.deg, 0*u.deg, 1*u.AU))
with pytest.raises(SunpyUserWarning, match="Tx is float32, and Ty is float16"):
hpc.make_3d()
def test_hpc_distance_cartesian():
# Test detection of distance in other representations
hpc1 = Helioprojective(CartesianRepresentation(0 * u.km, 0 * u.km, 1 * u.Mm))
assert isinstance(hpc1, Helioprojective)
assert isinstance(hpc1._data, CartesianRepresentation)
assert hpc1.make_3d() is hpc1
def test_hpc_distance_3D():
hpc1 = Helioprojective(1500 * u.arcsec, 0 * u.arcsec, 100 * u.Mm)
assert isinstance(hpc1, Helioprojective)
# Check that we have a 2D wrap180 representation
assert isinstance(hpc1._data, SphericalRepresentation)
# Check the attrs are correct
assert hpc1.Tx == 1500 * u.arcsec
assert hpc1.Ty == 0 * u.arcsec
hpc2 = hpc1.make_3d()
assert hpc2 is hpc1
def test_hpc_distance_off_limb():
hpc1 = Helioprojective(1500 * u.arcsec, 0 * u.arcsec,
observer=HeliographicStonyhurst(0*u.deg, 0*u.deg, 1*u.AU))
assert isinstance(hpc1, Helioprojective)
# Check that we have a 2D wrap180 representation
assert isinstance(hpc1._data, UnitSphericalRepresentation)
# Check the attrs are correct
assert hpc1.Tx == 1500 * u.arcsec
assert hpc1.Ty == 0 * u.arcsec
with pytest.warns(SunpyUserWarning, match="is all NaNs"):
hpc2 = hpc1.make_3d()
assert isinstance(hpc2._data, SphericalRepresentation)
# Check the attrs are correct
assert hpc2.Tx == 1500 * u.arcsec
assert hpc2.Ty == 0 * u.arcsec
assert_quantity_allclose(hpc2.distance, u.Quantity(np.nan, u.km))
def test_hpc_distance():
hpc1 = Helioprojective(0 * u.deg, 0 * u.arcsec,
observer=HeliographicStonyhurst(0*u.deg, 0*u.deg, 1*u.AU))
assert isinstance(hpc1, Helioprojective)
# Check that we have a 2D wrap180 representation
assert isinstance(hpc1._data, UnitSphericalRepresentation)
# Check the attrs are correct
assert hpc1.Tx == 0 * u.arcsec
assert hpc1.Ty == 0 * u.arcsec
hpc2 = hpc1.make_3d()
assert isinstance(hpc2._data, SphericalRepresentation)
# Check the attrs are correct
assert hpc2.Tx == 0 * u.arcsec
assert hpc2.Ty == 0 * u.arcsec
assert_quantity_allclose(hpc2.distance, DSUN_METERS - RSUN_METERS)