diff --git a/galaxies/galaxies.py b/galaxies/galaxies.py index acdf1c5..6e612cb 100644 --- a/galaxies/galaxies.py +++ b/galaxies/galaxies.py @@ -2,7 +2,7 @@ from astropy.coordinates import SkyCoord, Angle, FK5 from astroquery.ned import Ned import astropy.units as u -from astropy.table import Table +from astropy.table import Table, Column from astropy.io import fits from astropy.wcs import WCS import warnings @@ -83,7 +83,8 @@ def __init__(self, name, params=None): SkyCoord(t['RA(deg)'][0], t['DEC(deg)'][0], frame='fk5', unit='degree') - self.redshift = t['Redshift'][0] + self.redshift = t['Redshift'][0] * \ + u.dimensionless_unscaled except: warnings.warn("Unsuccessful query to NED") pass @@ -176,6 +177,36 @@ def position_angles(self, skycoord=None, ra=None, dec=None, return Angle(np.arctan2(Y, X)) + def to_table(self): + ''' + Return an `~astropy.table.Table` with the galactic parameters. + ''' + + tab = Table() + + tab["Name"] = Column([self.name], unit=None) + tab["Center Position"] = Column([self.center_position]) + + tab["Distance"] = Column([self.distance.value], + unit=self.distance.unit) + tab["Inclination"] = Column([self.inclination.value], + unit=self.inclination.unit) + tab["PA"] = Column([self.position_angle.value], + unit=self.position_angle.unit) + tab["Vsys"] = Column([self.vsys.value], + unit=self.vsys.unit) + + if self.redshift is not None: + tab["Redshift"] = Column([self.redshift.value], + unit=self.redshift.unit) + self.canonical_name = None + + if self.canonical_name is not None: + tab["Canonical Name"] = Column([self.canonical_name], + unit=None) + + return tab + # push or pull override table using astropy.table # Check name equivalencies