def _check_coords(self, position):
"""
Check that a source position is in the correct format.
If not, change to default coords.
Args:
position - source position. Must be 2-element list-like object with
position[0] = RA, position[1] = dec.
Acceptable formats for RA & dec are decimal degrees as
single floats or ints, or sexagisimal as 3-element
list-like object of (h,m,s), or strings with h:m:s.
Returns:
source_pos - Tuple of (RA,dec) in degrees.
"""
if len(position) != 2:
raise IOError('Position {:} has the wrong number of elements'.format(position))
else:
if not self._in_degrees(position):
if isinstance(position[0], (np.ndarray, list)):
# SkyCoord interprets lists as multiple sources, not hour/min/sec
position = (tuple(position[0]), tuple(position[1]))
try:
coords = SkyCoord(position[0], position[1], unit=(u.hourangle, u.deg))
new_pos = (coords.ra.deg, coords.dec.deg)
except ValueError as exc:
logger.critical(exc)
raise
except:
logger.critical('WARNING: RA and dec format cannot be read')
raise
else:
new_pos = (float(position[0]), float(position[1]))
if not (0 <= new_pos[0] <= 360):
raise ValueError('RA = {:.2f}, but should be in the range [0,360]'.format(new_pos[0]))
if not (-90 <= new_pos[1] <= 90):
raise ValueError('dec should be in the range [-90,90]')
return new_pos
评论列表
文章目录
