def process_vis_header (hdr): """handler for the visheader"""; # phase center (ra0,dec0) = hdr.phase_ref; set_state('ra0',value=ra0); set_state('dec0',value=dec0); # antenna positions pos = hdr.antenna_pos; if pos.rank != 2 or pos.shape[0] != 3: raise ValueError,'incorrectly shaped antenna_pos'; nant = pos.shape[1]; coords = ('x','y','z'); for iant in range(14): # range(nant), temporary hack sn = str(iant+1); # since some antennas may be missing from the tree, # ignore errors try: for (j,label) in enumerate(coords): print label+':'+sn, 'value = ',pos[j,iant] set_state(label+':'+sn,value=pos[j,iant]); except: pass; # array reference position #for (j,label) in enumerate(coords): # set_state(label+'0',value=pos[j,0]); print 'END OF READ_MSVIS_HEADER'
def process_vis_header (hdr): """handler for the visheader"""; # phase center (ra0,dec0) = hdr.phase_ref; set_state('ra0',value=ra0); set_state('dec0',value=dec0); # antenna positions pos = hdr.antenna_pos; if pos.rank != 2 or pos.shape[0] != 3: raise ValueError('incorrectly shaped antenna_pos'); nant = pos.shape[1]; coords = ('x','y','z'); for iant in range(14): # range(nant), temporary hack sn = str(iant+1); # since some antennas may be missing from the tree, # ignore errors try: for (j,label) in enumerate(coords): print(label+':'+sn, 'value = ',pos[j,iant]) set_state(label+':'+sn,value=pos[j,iant]); except: pass; # array reference position #for (j,label) in enumerate(coords): # set_state(label+'0',value=pos[j,0]); print('END OF READ_MSVIS_HEADER')
def process_vis_header (hdr): """handler for the standard MS visheader"""; # phase center (ra0,dec0) = hdr.phase_ref; set_state('ra0',value=ra0); set_state('dec0',value=dec0); # antenna positions pos = hdr.antenna_pos; if pos.rank != 2 or pos.shape[0] != 3: raise ValueError,'incorrectly shaped antenna_pos'; nant = pos.shape[1]; coords = ('x','y','z'); for iant in range(nant): sn = str(iant+1); # since some antennas may be missing from the tree, # ignore errors try: for (j,label) in enumerate(coords): set_state(label+'.'+sn,value=pos[j,iant]); except: pass; # array reference position for (j,label) in enumerate(coords): set_state(label+'0',value=pos[j,0]);
def process_vis_header(hdr): """handler for the standard MS visheader""" # phase center (ra0, dec0) = hdr.phase_ref set_state('ra0', value=ra0) set_state('dec0', value=dec0) # antenna positions pos = hdr.antenna_pos if pos.rank != 2 or pos.shape[0] != 3: raise ValueError, 'incorrectly shaped antenna_pos' nant = pos.shape[1] coords = ('x', 'y', 'z') for iant in range(nant): sn = str(iant + 1) # since some antennas may be missing from the tree, # ignore errors try: for (j, label) in enumerate(coords): set_state(label + '.' + sn, value=pos[j, iant]) except: pass # array reference position for (j, label) in enumerate(coords): set_state(label + '0', value=pos[j, 0])
def process_vis_header(hdr): """handler for the visheader""" try: # phase center (ra0, dec0) = hdr.phase_ref # print '[ReadVisHeader] phase centre: ',ra0,dec0; try: set_state('ra', value=ra0) set_state('dec', value=dec0) except: pass # antenna positions pos = hdr.antenna_pos if pos.ndim != 2 or pos.shape[0] != 3: raise ValueError('incorrectly shaped antenna_pos') nant = pos.shape[1] coords = ('x', 'y', 'z') for iant in range(nant): sn = ':num' + str(iant) # since some antennas may be missing from the tree, ignore errors try: for (j, label) in enumerate(coords): # print '[ReadVisHeader] ',label+sn, 'value = ',pos[j,iant] set_state(label + sn, value=pos[j, iant]) except: pass # also set nodes for the old naming convention sn = ":" + str(iant + 1) try: for (j, label) in enumerate(coords): # print '[ReadVisHeader] ',label+sn, 'value = ',pos[j,iant] set_state(label + sn, value=pos[j, iant]) except: pass # array reference position try: for (j, label) in enumerate(coords): set_state(label + '0', value=pos[j, 0]) except: pass # time extent (t0, t1) = hdr.time_extent # print '[ReadVisHeader] time extent: ',t0,t1; try: set_state('time0', value=t0) set_state('time1', value=t1) except: pass # freq range if isinstance(hdr.channel_freq, (int, float)): f0 = f1 = hdr.channel_freq else: f0, f1 = hdr.channel_freq[0], hdr.channel_freq[-1] # print '[ReadVisHeader] freq range: ',f0,f1; try: set_state('freq0', value=f0) set_state('freq1', value=f1) except: pass except: print(traceback.print_exc()) raise
def process_vis_header (hdr): """handler for the visheader"""; try: # phase center (ra0,dec0) = hdr.phase_ref; # print '[ReadVisHeader] phase centre: ',ra0,dec0; try: set_state('ra',value=ra0); set_state('dec',value=dec0); except: pass; # antenna positions pos = hdr.antenna_pos; if pos.ndim != 2 or pos.shape[0] != 3: raise ValueError,'incorrectly shaped antenna_pos'; nant = pos.shape[1]; coords = ('x','y','z'); for iant in range(nant): sn = ':num'+str(iant); # since some antennas may be missing from the tree, ignore errors try: for (j,label) in enumerate(coords): # print '[ReadVisHeader] ',label+sn, 'value = ',pos[j,iant] set_state(label+sn,value=pos[j,iant]); except: pass; # also set nodes for the old naming convention sn = ":" + str(iant+1); try: for (j,label) in enumerate(coords): # print '[ReadVisHeader] ',label+sn, 'value = ',pos[j,iant] set_state(label+sn,value=pos[j,iant]); except: pass; # array reference position try: for (j,label) in enumerate(coords): set_state(label+'0',value=pos[j,0]); except: pass; # time extent (t0,t1) = hdr.time_extent; # print '[ReadVisHeader] time extent: ',t0,t1; try: set_state('time0',value=t0); set_state('time1',value=t1); except: pass; # freq range if isinstance(hdr.channel_freq,(int,float)): f0=f1 = hdr.channel_freq; else: f0,f1 = hdr.channel_freq[0],hdr.channel_freq[-1]; # print '[ReadVisHeader] freq range: ',f0,f1; try: set_state('freq0',value=f0); set_state('freq1',value=f1); except: pass; except: print traceback.print_exc(); raise;