def parse_file(self):
"""
Read a source catalog data file.
"""
# open data file
fileName = os.path.join(self.get_directory(), '4c.dat')
with open(fileName, 'r') as catFile:
# read source info
lineNum = 1
line = catFile.readline()
while line:
try:
name = line[0:8]
raHours = int(line[10:12])
raMinutes = int(line[13:15])
raSeconds = float(line[16:20])
decSign = line[22]
decDegrees = int(line[23:25])
decMinutes = float(line[26:30])
alias = line[64:-1]
except ValueError:
raise RuntimeError(
"file %s, line %d incorectly formated [%s]" %
(fileName, lineNum, line))
name = name.strip()
name = ('4C' + name)
alias = alias.strip()
alias = alias.rstrip('*')
ra = astro.hms(raHours, raMinutes, raSeconds)
(decSeconds, decMinutes) = math.modf(decMinutes)
if decSign == '-':
sign = True
else:
sign = False
dec = astro.dms(sign, decDegrees, int(decMinutes),
(60 * decSeconds))
sourcePos = astro.equ_posn(ra, dec)
# precess coordinates from B1950
entry = CatalogEntry(
name,
transform.CelestialPosition(sourcePos,
epoch='B1950',
name=name))
self.source_map[name] = entry
if len(alias.strip()):
alias = alias.rstrip()
entry.alias_list.append(alias)
self.alias_map[alias] = entry
line = catFile.readline()
lineNum += 1
def parse_file(self):
"""
Read a source catalog data file.
"""
# open data file
fileName = os.path.join(self.get_directory(), 'PKS90.txt')
with open(fileName, 'r') as catFile:
# read source info
lineNum = 1
line = catFile.readline()
while line:
if line == '\n':
line = catFile.readline()
lineNum += 1
continue
try:
name = line[0:10]
alias0 = line[139:148]
alias1 = line[168:175]
raHours = int(line[10:12])
raMinutes = int(line[13:15])
raSeconds = float(line[16:20])
decSign = line[23]
decDegrees = int(line[24:26])
decMinutes = int(line[27:29])
decSeconds = float(line[30:34])
except ValueError:
raise RuntimeError(
"file %s, line %d incorectly formated [%s]" %
(fileName, lineNum, line))
ra = astro.hms(raHours, raMinutes, raSeconds)
if decSign == '-':
sign = True
else:
sign = False
dec = astro.dms(sign, decDegrees, decMinutes, decSeconds)
sourcePos = astro.equ_posn(ra, dec)
entry = CatalogEntry(
name, transform.CelestialPosition(sourcePos, name=name))
self.source_map[name] = entry
if len(alias0.strip()):
alias = alias0.rstrip()
entry.alias_list.append(alias)
self.alias_map[alias] = entry
if len(alias1.strip()):
alias = alias1.rstrip()
entry.alias_list.append(alias)
self.alias_map[alias] = entry
line = catFile.readline()
lineNum += 1
def parse_file(self):
"""
Read a source catalog data file.
"""
# open data file
fileName = os.path.join(self.get_directory(), 'pkscat.txt')
with open(fileName, 'r') as catFile:
# read source info
lineNum = 1
line = catFile.readline()
while line:
if line == '\n':
line = catFile.readline()
lineNum += 1
continue
try:
name = line[0:8]
alias = line[12:19]
raHours = int(line[22:24])
raMinutes = int(line[25:27])
raSeconds = float(line[28:32])
decSign = line[33]
decDegrees = int(line[34:36])
decMinutes = int(line[37:39])
decSeconds = int(line[40:42])
except ValueError:
raise RuntimeError(
"file %s, line %d incorectly formated [%s]" %
(fileName, lineNum, line))
ra = astro.hms(raHours, raMinutes, raSeconds)
if decSign == '-':
sign = True
else:
sign = False
dec = astro.dms(sign, decDegrees, decMinutes, decSeconds)
sourcePos = astro.equ_posn(ra, dec)
# precess coordinates from B1950
entry = CatalogEntry(
name,
transform.CelestialPosition(sourcePos,
epoch='B1950',
name=name))
self.source_map[name] = entry
if len(alias.strip()):
alias = alias.rstrip()
self.alias_map[alias] = entry
entry.alias_list.append(alias)
line = catFile.readline()
lineNum += 1
def parse_file(self):
"""
Read a source catalog data file.
"""
# open data file
fileName = os.path.join(self.get_directory(), '3c.dat')
with open(fileName, 'r') as catFile:
# read source info
lineNum = 1
line = catFile.readline()
while line:
try:
name = line[0:3]
raHours = int(line[12:14])
raMinutes = int(line[15:17])
raSeconds = float(line[18:22])
decSign = line[27]
decDegrees = int(line[28:30])
decMinutes = float(line[31:35])
except ValueError:
raise RuntimeError(
"file %s, line %d incorectly formated [%s]" %
(fileName, lineNum, line))
name = ('3C' + name.strip())
ra = astro.hms(raHours, raMinutes, raSeconds)
(decSeconds, decMinutes) = math.modf(decMinutes)
if decSign == '-':
sign = True
else:
sign = False
dec = astro.dms(sign, decDegrees, int(decMinutes),
(60 * decSeconds))
sourcePos = astro.equ_posn(ra, dec)
# precess coordinates from B1950
entry = CatalogEntry(
name,
transform.CelestialPosition(sourcePos,
epoch='B1950',
name=name))
self.source_map[name] = entry
line = catFile.readline()
lineNum += 1
def parse_file(self):
"""
Read a source catalog data file.
"""
# open data file
fileName = os.path.join(self.get_directory(), 'lwa_catalog.dat')
with open(fileName, 'r') as catFile:
# read source info
lineNum = 0
for line in catFile:
lineNum += 1
if line.startswith('#') or line.isspace():
continue
try:
name = line[0:8]
raHours = int(line[9:11], 10)
raMinutes = int(line[12:14], 10)
raSeconds = float(line[15:20])
decSign = line[21]
decDegrees = int(line[22:24], 10)
decMinutes = int(line[25:27], 10)
decSeconds = float(line[28:32])
except ValueError as err:
raise RuntimeError(
"file %s, line %d incorectly formated: \'%s\' : %s]" %
(fileName, lineNum, line, err))
name = name.rstrip()
ra = astro.hms(raHours, raMinutes, raSeconds)
if decSign == '-':
sign = True
else:
sign = False
dec = astro.dms(sign, decDegrees, decMinutes, decSeconds)
entry = CatalogEntry(
name, transform.CelestialPosition((ra, dec), name=name))
self.source_map[name] = entry
aliasList = line[34:].split()
if len(aliasList) > 0:
for alias in aliasList:
entry.alias_list.append(alias)
self.alias_map[alias] = entry
station = stations.lwasv
elif site == 'ovrolwa':
station = stations.lwa1
station.name = 'OVRO-LWA'
station.lat, station.lon, station.elev = ('37.2397808', '-118.2816819', 1183.4839)
else:
raise RuntimeError("Unknown site name: %s" % site)
site = transform.GeographicalPosition((station.long*180.0/math.pi, station.lat*180.0/math.pi), name=station.name)
sun_pos = transform.PlanetaryPosition('Sun')
sun_dir = transform.PointingDirection(sun_pos, site)
jup_pos = transform.PlanetaryPosition('Jupiter')
jup_dir = transform.PointingDirection(jup_pos, site)
gal_pos = transform.CelestialPosition((astro.hms(17, 42, 0.1),
astro.dms(True, 28, 55, 8)), name = 'SgrA')
gal_dir = transform.PointingDirection(gal_pos, site)
tau_pos = transform.CelestialPosition((astro.hms(5, 34, 0.5),
astro.dms(False, 22, 0, 52)), name = 'TauA')
tau_dir = transform.PointingDirection(tau_pos, site)
cyg_pos = transform.CelestialPosition((astro.hms(19, 59, 0.8),
astro.dms(False, 40, 44, 2)), name = 'CygA')
cyg_dir = transform.PointingDirection(cyg_pos, site)
cas_pos = transform.CelestialPosition((astro.hms(23, 23, 0.7),
astro.dms(False, 58, 49, 16)), name = 'CasA')
cas_dir = transform.PointingDirection(cas_pos, site)
def parse_file(self):
"""
Read a source catalog data file.
"""
debug = False
# open data file
fileName = os.path.join(self.get_directory(), 'psrcat.db')
with open(fileName, 'r') as catFile:
# read source info
psrb = None
psrj = None
ra = None
dec = None
bad = False
for line in catFile:
if line.startswith('PSRB'):
psrb = line.split()[1]
if line.startswith('PSRJ'):
psrj = line.split()[1]
if line.startswith('RAJ'):
rastr = line.split()[1]
sp = rastr.split(':')
if len(sp) == 3:
(raHours, raMinutes, raSeconds) = sp
elif len(sp) == 2:
(raHours, raMinutes) = sp
raSeconds = 0.0
else:
if debug:
print('Bad format for RAJ line : ' + line)
bad = True
raHours = int(raHours)
raMinutes = int(raMinutes)
raSeconds = float(raSeconds)
try:
ra = astro.hms(raHours, raMinutes, raSeconds)
except:
if debug:
print('PSRCAT: Bad RA for ', psrj, " : ", rastr)
bad = True
if line.startswith('DECJ'):
decstr = line.split()[1]
sp = decstr.split(':')
if len(sp) == 3:
(decDegrees, decMinutes, decSeconds) = sp
elif len(sp) == 2:
(decDegrees, decMinutes) = sp
decSeconds = 0.0
else:
if debug:
print('PSRCAT: Bad format for DECJ line : ' + line)
bad = True
continue
if decDegrees.startswith('-'):
decDegrees = decDegrees[1:]
sign = True
else:
sign = False
decDegrees = int(decDegrees)
decMinutes = int(decMinutes)
decSeconds = float(decSeconds)
try:
dec = astro.dms(sign, decDegrees, decMinutes,
decSeconds)
except:
if debug:
print('PSRCAT: Bad DEC for ', psrj, " : ", decstr)
bad = True
if line.startswith('@-'):
# New source, save current source info
if psrb is not None:
name = psrb
alias = psrj
else:
name = psrj
alias = None
if ra is None or dec is None:
# These pulsars don't have RAJ, DECJ
# I think they may all have ecliptic positions
# which should be converted to ra,dec but I'm
# going to ignore them for now. -- paulr
#print "PSRCAT: No position for pulsar ",name
bad = True
# Add source to list if good.
if not bad:
sourcePos = astro.equ_posn(ra, dec)
entry = CatalogEntry(
name,
transform.CelestialPosition(sourcePos, name=name))
self.source_map[name] = entry
if debug:
print('Added ', name)
if alias is not None:
alias = alias.rstrip()
self.alias_map[alias] = entry
entry.alias_list.append(alias)
if debug:
print('Alias : ', alias.rstrip())
# Clear out vars for next source
psrb = None
psrj = None
ra = None
dec = None
bad = False
print('---------------------------------------------------------------')
print('RA: %s (%0.3f)' % (saturn_ra, saturn_equ.ra))
print('DEC: %s (%0.3f)' % (saturn_dec, saturn_equ.dec))
print('Ecl longitude: %s (%0.3f)' % (saturn_lng, saturn_ecl.lng))
print('Ecl latitude: %s (%0.3f)' % (saturn_lat, saturn_ecl.lat))
print('Rise: %s (%0.3f) [%s]' % (saturn_utc_rise, saturn_rst.rise, saturn_lcl_rise))
print('Transit: %s (%0.3f) [%s]' % (saturn_utc_trans, saturn_rst.transit, saturn_lcl_trans))
print('Set: %s (%0.3f) [%s]' % (saturn_utc_set, saturn_rst.set, saturn_lcl_set))
print('Azimuth: %0.3f %s' % (saturn_hrz.az, astro.hrz_to_nswe(saturn_hrz)))
print('Altitude: %0.3f' % saturn_hrz.alt)
print('Zenith: %0.3f' % saturn_hrz.zen())
print('Sun angle: %0.3f' % saturn_sun_ang)
# calculate SgrA phenomena
sgra_j2000_equ = astro.equ_posn(astro.hms(17, 42, 48.1), astro.dms(True, 28, 55, 8))
sgra_equ = astro.get_apparent_posn(sgra_j2000_equ, utc)
sgra_rst = astro.get_object_rst(utc, lwa_lnlat, sgra_equ)
(sgra_utc_rise, sgra_utc_set, sgra_utc_trans) = sgra_rst.format()
sgra_lcl_rise = sgra_utc_rise.to_zone()
sgra_lcl_trans = sgra_utc_trans.to_zone()
sgra_lcl_set = sgra_utc_set.to_zone()
(sgra_ra, sgra_dec) = sgra_equ.format()
sgra_hrz = sgra_equ.to_hrz(lwa_lnlat, utc)
sgra_gal = sgra_equ.to_gal(utc)
(sgra_l, sgra_b) = sgra_gal.format()
sgra_sun_ang = sun_equ.angular_separation(sgra_equ)
print('---------------------------------------------------------------')
print('Ecl latitude: %s (%0.3f)' % (saturn_lat, saturn_ecl.lat))
print('Rise: %s (%0.3f) [%s]' %
(saturn_utc_rise, saturn_rst.rise, saturn_lcl_rise))
print('Transit: %s (%0.3f) [%s]' %
(saturn_utc_trans, saturn_rst.transit, saturn_lcl_trans))
print('Set: %s (%0.3f) [%s]' %
(saturn_utc_set, saturn_rst.set, saturn_lcl_set))
print('Azimuth: %0.3f %s' %
(saturn_hrz.az, astro.hrz_to_nswe(saturn_hrz)))
print('Altitude: %0.3f' % saturn_hrz.alt)
print('Zenith: %0.3f' % saturn_hrz.zen())
print('Sun angle: %0.3f' % saturn_sun_ang)
# calculate SgrA phenomena
sgra_j2000_equ = astro.equ_posn(astro.hms(17, 42, 48.1),
astro.dms(True, 28, 55, 8))
sgra_equ = astro.get_apparent_posn(sgra_j2000_equ, utc)
sgra_rst = astro.get_object_rst(utc, lwa_lnlat, sgra_equ)
(sgra_utc_rise, sgra_utc_set, sgra_utc_trans) = sgra_rst.format()
sgra_lcl_rise = sgra_utc_rise.to_zone()
sgra_lcl_trans = sgra_utc_trans.to_zone()
sgra_lcl_set = sgra_utc_set.to_zone()
(sgra_ra, sgra_dec) = sgra_equ.format()
sgra_hrz = sgra_equ.to_hrz(lwa_lnlat, utc)
sgra_gal = sgra_equ.to_gal(utc)
(sgra_l, sgra_b) = sgra_gal.format()
sgra_sun_ang = sun_equ.angular_separation(sgra_equ)