#!/usr/bin/env python import antenna_array import correlator import numpy as np import matplotlib.pyplot as plt BASELINE_LENGTH = 2.55 arr = antenna_array.AntennaArray(BASELINE_LENGTH * 2 * np.pi) ref = arr.each_pair_phase_difference_at_angle(0) corr = correlator.Correlator(ref, arr) response = corr.many_directions(-np.pi, np.pi, 900) x = response.keys() x.sort() y = [] for i_idx, i_val in enumerate(x): y.append(response[i_val]) plt.plot(x, y) plt.show()
npol = du.npol ant_data = du.data core_ind = NP.logical_and((NP.abs(antpos[:, 0]) < max_antenna_radius), (NP.abs(antpos[:, 1]) < max_antenna_radius)) # core_ind = NP.logical_and((NP.abs(antpos[:,0]) <= NP.max(NP.abs(antpos[:,0]))), (NP.abs(antpos[:,1]) < NP.max(NP.abs(antpos[:,1])))) antid = antid[core_ind] antpos = antpos[core_ind, :] ant_info = NP.hstack((antid.reshape(-1, 1), antpos)) n_antennas = ant_info.shape[0] ant_data = ant_data[:, core_ind, :, :] with PyCallGraph(output=graphviz, config=config): ants = [] aar = AA.AntennaArray() for i in xrange(n_antennas): ant = AA.Antenna('{0:0d}'.format(int(ant_info[i, 0])), lat, ant_info[i, 1:], f0, nsamples=nts) ant.f = ant.f0 + DSP.spectax(2 * nts, dt, shift=True) ants += [ant] aar = aar + ant aar.grid() antpos_info = aar.antenna_positions(sort=True) if max_n_timestamps is None:
def setUp(self): self.ant_array = antenna_array.AntennaArray(5.6) # spacing of 5.6 wavelengths apart.