def test_frames(tmpdir):
frames = [
cf.CelestialFrame(reference_frame=coord.ICRS()),
cf.CelestialFrame(reference_frame=coord.FK5(
equinox=time.Time('2010-01-01'))),
cf.CelestialFrame(reference_frame=coord.FK4(
equinox=time.Time('2010-01-01'), obstime=time.Time('2015-01-01'))),
cf.CelestialFrame(reference_frame=coord.FK4NoETerms(
equinox=time.Time('2010-01-01'), obstime=time.Time('2015-01-01'))),
cf.CelestialFrame(reference_frame=coord.Galactic()),
cf.CelestialFrame(
reference_frame=coord.Galactocentric(galcen_distance=5.0 * u.m,
galcen_ra=45 * u.deg,
galcen_dec=1 * u.rad,
z_sun=3 * u.pc,
roll=3 * u.deg)),
cf.CelestialFrame(
reference_frame=coord.GCRS(obstime=time.Time('2010-01-01'),
obsgeoloc=[1, 3, 2000] * u.pc,
obsgeovel=[2, 1, 8] * (u.m / u.s))),
cf.CelestialFrame(reference_frame=coord.CIRS(
obstime=time.Time('2010-01-01'))),
cf.CelestialFrame(reference_frame=coord.ITRS(
obstime=time.Time('2022-01-03'))),
cf.CelestialFrame(reference_frame=coord.PrecessedGeocentric(
obstime=time.Time('2010-01-01'),
obsgeoloc=[1, 3, 2000] * u.pc,
obsgeovel=[2, 1, 8] * (u.m / u.s)))
]
tree = {'frames': frames}
helpers.assert_roundtrip_tree(tree, tmpdir)
def create_test_frames():
"""Creates an array of frames to be used for testing."""
# Suppress warnings from astropy that are caused by having 'dubious' dates
# that are too far in the future. It's not a concern for the purposes of
# unit tests. See issue #5809 on the astropy GitHub for discussion.
from astropy._erfa import ErfaWarning
warnings.simplefilter("ignore", ErfaWarning)
frames = [
cf.CelestialFrame(reference_frame=coord.ICRS()),
cf.CelestialFrame(reference_frame=coord.FK5(
equinox=time.Time('2010-01-01'))),
cf.CelestialFrame(reference_frame=coord.FK4(
equinox=time.Time('2010-01-01'), obstime=time.Time('2015-01-01'))),
cf.CelestialFrame(reference_frame=coord.FK4NoETerms(
equinox=time.Time('2010-01-01'), obstime=time.Time('2015-01-01'))),
cf.CelestialFrame(reference_frame=coord.Galactic()),
cf.CelestialFrame(reference_frame=coord.Galactocentric(
# A default galcen_coord is used since none is provided here
galcen_distance=5.0 * u.m,
z_sun=3 * u.pc,
roll=3 * u.deg)),
cf.CelestialFrame(
reference_frame=coord.GCRS(obstime=time.Time('2010-01-01'),
obsgeoloc=[1, 3, 2000] * u.pc,
obsgeovel=[2, 1, 8] * (u.m / u.s))),
cf.CelestialFrame(reference_frame=coord.CIRS(
obstime=time.Time('2010-01-01'))),
cf.CelestialFrame(reference_frame=coord.ITRS(
obstime=time.Time('2022-01-03'))),
cf.CelestialFrame(reference_frame=coord.PrecessedGeocentric(
obstime=time.Time('2010-01-01'),
obsgeoloc=[1, 3, 2000] * u.pc,
obsgeovel=[2, 1, 8] * (u.m / u.s))),
cf.StokesFrame(),
cf.TemporalFrame(time.Time("2011-01-01"))
]
return frames
*args, frame='icrs', **kwargs)
GalacticCoordGenerator = lambda *args, **kwargs: coord.SkyCoord(
*args, frame='galactic', **kwargs)
FK5CoordGenerator = lambda *args, **kwargs: coord.SkyCoord(
*args, frame='fk5', **kwargs)
FK4CoordGenerator = lambda *args, **kwargs: coord.SkyCoord(
*args, frame='fk4', **kwargs)
ICRSCoord = coord.SkyCoord
CoordClasses = (coord.SkyCoord, BaseCoordinateFrame)
except ImportError:
from astropy.coordinates import SphericalCoordinatesBase as BaseCoordinateFrame
ICRSCoordGenerator = lambda *args, **kwargs: coord.ICRS(*args, **kwargs)
GalacticCoordGenerator = lambda *args, **kwargs: coord.Galactic(
*args, **kwargs)
FK5CoordGenerator = lambda *args, **kwargs: coord.FK5(*args, **kwargs)
FK4CoordGenerator = lambda *args, **kwargs: coord.FK4(*args, **kwargs)
ICRSCoord = coord.ICRS
CoordClasses = (coord.SphericalCoordinatesBase, )
from ..exceptions import TimeoutError
from .. import version
__all__ = [
'send_request', 'parse_coordinates', 'parse_radius', 'TableList',
'suppress_vo_warnings', 'validate_email'
]
def send_request(url,
data,
timeout,
def test___init__(self) -> None:
"""Test method ``__init__``.
Copying from astropy docs
"""
# -----------
c = icoord.InterpolatedSkyCoord(
[10] * self.num,
[20] * self.num,
unit="deg",
affine=self.affine,
) # defaults to ICRS frame
self._test_isc(c)
# -----------
c = icoord.InterpolatedSkyCoord(
[1, 2, 3, 4],
[-30, 45, 8, 16],
frame="icrs",
unit="deg",
affine=self.affine[:4],
) # 4 coords
self._test_isc(c)
# -----------
coords = [
"1:12:43.2 +31:12:43",
"1:12:43.2 +31:12:43",
"1:12:43.2 +31:12:43",
"1 12 43.2 +31 12 43",
]
c = icoord.InterpolatedSkyCoord(
coords,
frame=coord.FK4,
unit=(u.hourangle, u.deg),
obstime="J1992.21",
affine=self.affine[:4],
)
self._test_isc(c)
# -----------
c = icoord.InterpolatedSkyCoord(
["1h12m43.2s +1d12m43s"] * self.num,
frame=coord.Galactic,
affine=self.affine,
) # Units from string
self._test_isc(c)
# # -----------
c = icoord.InterpolatedSkyCoord(
frame="galactic",
l=["1h12m43.2s"] * self.num,
b="+1d12m43s", # NOT BROADCASTING THIS ONE
affine=self.affine,
)
self._test_isc(c)
# -----------
ra = coord.Longitude([1, 2, 3, 4], unit=u.deg) # Could also use Angle
dec = np.array([4.5, 5.2, 6.3, 7.4]) * u.deg # Astropy Quantity
c = icoord.InterpolatedSkyCoord(
ra,
dec,
frame="icrs",
affine=self.affine[:4],
)
self._test_isc(c)
# -----------
c = icoord.InterpolatedSkyCoord(
frame=coord.ICRS,
ra=ra,
dec=dec,
obstime="2001-01-02T12:34:56",
affine=self.affine[:4],
)
self._test_isc(c)
# -----------
c = coord.FK4(
[1] * self.num * u.deg,
2 * u.deg,
) # Uses defaults for obstime, equinox
c = icoord.InterpolatedSkyCoord(
c,
obstime="J2010.11",
equinox="B1965",
affine=self.affine,
) # Override defaults
self._test_isc(c)
# -----------
c = icoord.InterpolatedSkyCoord(
w=[0] * self.num,
u=1,
v=2,
unit="kpc",
frame="galactic",
representation_type="cartesian",
affine=self.affine,
)
self._test_isc(c, representation_type=coord.CartesianRepresentation)
# -----------
c = icoord.InterpolatedSkyCoord(
[
coord.ICRS(ra=1 * u.deg, dec=2 * u.deg),
coord.ICRS(ra=3 * u.deg, dec=4 * u.deg),
] * (self.num // 2),
affine=self.affine,
)
self._test_isc(c)
from astropy.tests.helper import pytest
import astropy.coordinates as coord
import astropy.units as u
from astropy.table import Table
import sys
import os
import re
import requests
from ...exceptions import TableParseError
from distutils.version import LooseVersion
is_python3 = (sys.version_info >= (3, ))
GALACTIC_COORDS = coord.Galactic(l=-67.02084, b=-29.75447, unit=(u.deg, u.deg))
ICRS_COORDS = coord.ICRS("05h35m17.3s -05h23m28s")
FK4_COORDS = coord.FK4(ra=84.90759, dec=-80.89403, unit=(u.deg, u.deg))
FK5_COORDS = coord.FK5(ra=83.82207, dec=-80.86667, unit=(u.deg, u.deg))
DATA_FILES = {
'id': 'query_id.data',
'coo': 'query_coo.data',
'cat': 'query_cat.data',
'bibobj': 'query_bibobj.data',
'bibcode': 'query_bibcode.data',
'error': 'query_error.data',
'sample': 'query_sample.data',
'region': 'query_sample_region.data',
}
class MockResponseSimbad(MockResponse):
