#!/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.
