barentsen · September 20, 2016 17:32
diff --git a/compute-galactic-plane-coordinates.py b/compute-galactic-plane-coordinates.py
 """Compute ra and dec for the Galactic Plane.

 This script computes the ICRS Right Ascension (ra) and Declination (dec)
 coordinates for the entire Galactic Plane and writes them to a csv file.
 """
 import numpy as np
 import pandas as pd
 from astropy.coordinates import SkyCoord

 galactic_longitudes = np.arange(start=0, stop=360, step=0.1)
 galactic_latitudes = [0] * len(galactic_longitudes)
 icrs = SkyCoord(galactic_longitudes, galactic_latitudes,
                unit="deg", frame="galactic").icrs
 df = pd.DataFrame({'glon': galactic_longitudes,
                   'glat': galactic_latitudes,
                   'ra': icrs.ra,
                   'dec': icrs.dec})
 df[['glon', 'glat', 'ra', 'dec']].to_csv('galactic-plane-coordinates-icrs.csv',
                                         index=False, float_format='%.14f')
diff --git a/galactic-plane-coordinates-icrs.csv b/galactic-plane-coordinates-icrs.csv
	"""Compute ra and dec for the Galactic Plane.

	This script computes the ICRS Right Ascension (ra) and Declination (dec)
	coordinates for the entire Galactic Plane and writes them to a csv file.
	"""
	import numpy as np
	import pandas as pd
	from astropy.coordinates import SkyCoord

	galactic_longitudes = np.arange(start=0, stop=360, step=0.1)
	galactic_latitudes = [0] * len(galactic_longitudes)
	icrs = SkyCoord(galactic_longitudes, galactic_latitudes,
	unit="deg", frame="galactic").icrs
	df = pd.DataFrame({'glon': galactic_longitudes,
	'glat': galactic_latitudes,
	'ra': icrs.ra,
	'dec': icrs.dec})
	df[['glon', 'glat', 'ra', 'dec']].to_csv('galactic-plane-coordinates-icrs.csv',
	index=False, float_format='%.14f')