def printBaselines(deploy_index=default_deploy,
calculate_phase=False,
verbose=False):
'''
Currently only intended to plot the most recent station with the three pulsers that we used for it.
'''
try:
antennas_physical, antennas_phase_hpol, antennas_phase_vpol = info.loadAntennaLocationsENU(
deploy_index=deploy_index, check=False)
print(
info.loadAntennaLocationsENU(deploy_index=deploy_index,
check=False))
colors = ['b', 'g', 'r', 'c']
if calculate_phase:
modes = ['physical', 'hpol', 'vpol']
else:
modes = ['physical']
for mode in modes:
if mode == 'hpol':
enu = antennas_phase_hpol
elif mode == 'vpol':
enu = antennas_phase_vpol
else:
enu = antennas_physical
print('%s baselines:' % mode)
for pair in [[0, 1], [0, 2], [0, 3], [1, 2], [1, 3], [2, 3]]:
try:
baseline = numpy.sqrt(
(enu[pair[0]][0] - enu[pair[1]][0])**2 +
(enu[pair[0]][1] - enu[pair[1]][1])**2 +
(enu[pair[0]][2] - enu[pair[1]][2])**2)
print('\t%s : %0.3f m' % (str(pair), baseline))
except:
if verbose:
print('%s failed or not present.' % mode)
except Exception as e:
print('\nError in %s' % inspect.stack()[0][3])
print(e)
exc_type, exc_obj, exc_tb = sys.exc_info()
fname = os.path.split(exc_tb.tb_frame.f_code.co_filename)[1]
print(exc_type, fname, exc_tb.tb_lineno)
print(valley_source_key)
print('Signal Delays:')
print(numpy.mean(timeshifts, axis=1))
print('Cable Delays:')
print(cable_delays_differences)
plt.legend(fontsize=8)
#### AIRPLANES ####
if len(included_airplanes) > 0:
print('Loading known plane locations.')
known_planes, calibrated_trigtime, output_tracks = pt.getKnownPlaneTracks(
ignore_planes=[]
) # ['1728-62026','1773-14413','1773-63659','1774-88800','1783-28830','1784-7166']#'1774-88800','1728-62026'
origin = info.loadAntennaZeroLocation(deploy_index=deploy_index)
antennas_physical, antennas_phase_hpol, antennas_phase_vpol = info.loadAntennaLocationsENU(
deploy_index=deploy_index)
plane_polys = {}
interpolated_plane_locations = {}
colors = plt.rcParams["axes.prop_cycle"].by_key()["color"]
for index, key in enumerate(list(calibrated_trigtime.keys())):
if key not in included_airplanes:
continue
run = int(key.split('-')[0])
reader = Reader(datapath, run)
eventids = known_planes[key]['eventids'][:, 1]
enu = numpy.array(
from beacon.tools import info
from beacon_deployment_2021.scripts.plot_array import plotStation, printBaselines
default_deploy = info.returnDefaultDeploy()
if __name__ == '__main__':
try:
deploy_index_1 = os.path.join(os.environ['BEACON_DEPLOYMENT_DIR'],
'config/deploy_30.json')
deploy_index_2 = os.path.join(os.environ['BEACON_DEPLOYMENT_DIR'],
'config/rtk-gps-day1-june20-2021.json')
printBaselines(deploy_index=deploy_index_1, calculate_phase=True)
printBaselines(deploy_index=deploy_index_2, calculate_phase=False)
#fig, ax = plotStation(deploy_index=deploy_file,plot_phase=False)
antennas_physical_1, antennas_phase_hpol_1, antennas_phase_vpol_1 = info.loadAntennaLocationsENU(
deploy_index=deploy_index_1)
antennas_physical_2, antennas_phase_hpol_2, antennas_phase_vpol_2 = info.loadAntennaLocationsENU(
deploy_index=deploy_index_2)
colors = ['b', 'g', 'r', 'c']
fig = plt.figure()
plt.title('New Physical Compared to Old')
ax = fig.add_subplot(111, projection='3d')
for i, a in antennas_physical_1.items():
ax.scatter(a[0],
a[1],
a[2],
marker='o',
color=colors[i],
def compareStations(deploy_index_1, deploy_index_2, plot_phase=False):
'''
Currently only intended to plot the most recent station with the three pulsers that we used for it.
'''
try:
antennas_physical_1, antennas_phase_hpol_1, antennas_phase_vpol_1 = info.loadAntennaLocationsENU(
deploy_index=deploy_index_1)
antennas_physical_2, antennas_phase_hpol_2, antennas_phase_vpol_2 = info.loadAntennaLocationsENU(
deploy_index=deploy_index_2)
print(info.loadAntennaLocationsENU(deploy_index=deploy_index))
colors = ['b', 'g', 'r', 'c']
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
for i, a in antennas_physical_1.items():
ax.scatter(a[0],
a[1],
a[2],
marker='o',
color=colors[i],
label='Physical %i' % i,
alpha=0.8)
for i, a in antennas_physical_2.items():
ax.scatter(a[0],
a[1],
a[2],
marker='*',
color=colors[i],
label='Physical %i' % i,
alpha=0.8)
ax.plot([antennas_physical_1[i][0], antennas_phase_physical_1[i][0]],
[antennas_physical_1[i][1], antennas_phase_physical_1[i][1]],
[antennas_physical_1[i][2], antennas_phase_physical_1[i][2]],
color=colors[i],
linestyle='--',
alpha=0.5)
if plot_phase == True:
for i, a in antennas_phase_hpol_1.items():
try:
ax.plot(
[antennas_hpol_1[i][0], antennas_phase_hpol_1[i][0]],
[antennas_hpol_1[i][1], antennas_phase_hpol_1[i][1]],
[antennas_hpol_1[i][2], antennas_phase_hpol_1[i][2]],
color=colors[i],
linestyle='--',
alpha=0.5)
except:
pass
for i, a in antennas_phase_vpol_1.items():
try:
ax.plot(
[antennas_vpol_1[i][0], antennas_phase_vpol_1[i][0]],
[antennas_vpol_1[i][1], antennas_phase_vpol_1[i][1]],
[antennas_vpol_1[i][2], antennas_phase_vpol_1[i][2]],
color=colors[i],
linestyle='--',
alpha=0.5)
except:
pass
ax.set_xlabel('E (m)')
ax.set_ylabel('N (m)')
ax.set_zlabel('Relative Elevation (m)')
plt.legend()
return fig, ax
except Exception as e:
print('\nError in %s' % inspect.stack()[0][3])
print(e)
exc_type, exc_obj, exc_tb = sys.exc_info()
fname = os.path.split(exc_tb.tb_frame.f_code.co_filename)[1]
print(exc_type, fname, exc_tb.tb_lineno)
def plotStation(deploy_index=default_deploy, plot_phase=False):
'''
Currently only intended to plot the most recent station with the three pulsers that we used for it.
'''
try:
antennas_physical, antennas_phase_hpol, antennas_phase_vpol = info.loadAntennaLocationsENU(
deploy_index=deploy_index)
print(info.loadAntennaLocationsENU(deploy_index=deploy_index))
colors = ['b', 'g', 'r', 'c']
fig = plt.figure()
if str(deploy_index).isdigit():
fig.canvas.set_window_title('Antenna Locations: Calibration %i' %
deploy_index)
else:
fig.canvas.set_window_title(
'Antenna Locations: %s' %
os.path.split(deploy_index)[-1].replace('.json', '').replace(
'_', ' '))
ax = fig.add_subplot(111, projection='3d')
for i, a in antennas_physical.items():
ax.scatter(a[0],
a[1],
a[2],
marker='o',
color=colors[i],
label='Physical %i' % i,
alpha=0.8)
if plot_phase == True:
for i, a in antennas_phase_hpol.items():
ax.plot([antennas_physical[i][0], antennas_phase_hpol[i][0]],
[antennas_physical[i][1], antennas_phase_hpol[i][1]],
[antennas_physical[i][2], antennas_phase_hpol[i][2]],
color=colors[i],
linestyle='--',
alpha=0.5)
ax.scatter(a[0],
a[1],
a[2],
marker='*',
color=colors[i],
label='%s Phase Center %i' % ('Hpol', i),
alpha=0.8)
for i, a in antennas_phase_vpol.items():
ax.plot([antennas_physical[i][0], antennas_phase_vpol[i][0]],
[antennas_physical[i][1], antennas_phase_vpol[i][1]],
[antennas_physical[i][2], antennas_phase_vpol[i][2]],
color=colors[i],
linestyle='--',
alpha=0.5)
ax.scatter(a[0],
a[1],
a[2],
marker='^',
color=colors[i],
label='%s Phase Center %i' % ('Vpol', i),
alpha=0.8)
ax.set_xlabel('E (m)')
ax.set_ylabel('N (m)')
ax.set_zlabel('Relative Elevation (m)')
plt.legend()
return fig, ax
except Exception as e:
print('\nError in %s' % inspect.stack()[0][3])
print(e)
exc_type, exc_obj, exc_tb = sys.exc_info()
fname = os.path.split(exc_tb.tb_frame.f_code.co_filename)[1]
print(exc_type, fname, exc_tb.tb_lineno)