37 lines
1.0 KiB
Python
37 lines
1.0 KiB
Python
|
from astropy.coordinates import SkyCoord, ITRS, representation
|
||
|
|
from astropy.time import Time
|
||
|
|
import json
|
||
|
|
import sys
|
||
|
|
|
||
|
|
|
||
|
|
def convert_to_itrs(coord):
|
||
|
|
c = SkyCoord(coord[0], coord[1], unit='degree', frame='icrs')
|
||
|
|
c_itrs = c.transform_to(ITRS(obstime=Time.now()))
|
||
|
|
rep = c_itrs.represent_as(representation.UnitSphericalRepresentation)
|
||
|
|
return [rep.lat.deg, rep.lon.deg]
|
||
|
|
|
||
|
|
|
||
|
|
def points_to_lines(points):
|
||
|
|
lines = []
|
||
|
|
for x in range(len(points) - 1):
|
||
|
|
lines.append([points[x], points[x+1]])
|
||
|
|
return lines
|
||
|
|
|
||
|
|
|
||
|
|
with open(sys.argv[1]) as f:
|
||
|
|
constellations = json.load(f)['features']
|
||
|
|
|
||
|
|
output = []
|
||
|
|
|
||
|
|
for const in constellations:
|
||
|
|
for line in const['geometry']['coordinates']:
|
||
|
|
transformed_line = []
|
||
|
|
for point in line:
|
||
|
|
transformed_line.append(convert_to_itrs(point))
|
||
|
|
|
||
|
|
line_combined = points_to_lines(transformed_line)
|
||
|
|
for l in line_combined: # noqa
|
||
|
|
output.append(f'{l[0][0]} {l[0][1]} {l[1][0]} {l[1][1]} # {const["id"]}') # noqa
|
||
|
|
|
||
|
|
print('\n'.join(output))
|