-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathread-coordinate-reg.py
More file actions
40 lines (34 loc) · 999 Bytes
/
read-coordinate-reg.py
File metadata and controls
40 lines (34 loc) · 999 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
from __future__ import print_function
import numpy as np
import json
import os
from astropy.io import fits
from astropy import wcs
from astropy.wcs import WCS
from astropy import coordinates as coord
from astropy import units as u
import argparse
import sys
#position = "halpha-emitter.reg"
position = "HII-region.reg"
ra, dec = [], []
f = open(position, 'r')
header1 = f.readline()
header2 = f.readline()
header3 = f.readline()
for line in f:
line = line.strip()
columns = line.split()
coor = line.split("(")[-1].split(")")[0]
try:
ra1, dec1 = coor.split(",")[0:2]
except ValueError:
continue
coor_degree = coord.SkyCoord(ra1, dec1, unit=(u.hourangle, u.degree))
ra.append(coor_degree.ra.degree)
dec.append(coor_degree.dec.degree)
asciifile = "HII-region-galaxy.txt"
file=open(asciifile,'w') #create file
for x,y in zip(ra, dec):
file.write('%f %f\n'%(x,y)) #assume you separate columns by tabs
file.close() #close file