image thumbnail

read_teqc_compact

by

 

17 Apr 2011 (Updated )

Reads teqc's QC compact files more rapidly -- without for loops.

doit.m
base_dir = 'c:\work\other\fx\read_teqc_compact\';
addpath(base_dir)

base_name = 'Ball1830';
tic
data = read_teqc_compact (base_dir, base_name, 'mp1');
toc

prn = unique(data.prn);
figure
for i=1:numel(prn)
  idx = (data.prn == prn(i));
  subplot(3,1,1), plot(data.epoch(idx), data.elev(idx), '-k')
  subplot(3,1,2), plot(data.epoch(idx), data.azim(idx), '-k')
  subplot(3,1,3), plot(data.epoch(idx), data.obs(idx),  '-k', 'LineWidth',1/2)
  subplot(3,1,1), ylabel('Elev.')
  subplot(3,1,2), ylabel('Azim.')
  subplot(3,1,3), ylabel('MP1'), xlabel('Time')
  subplot(3,1,1), title(sprintf('PRN %d', prn(i)))
  pause
end

i = 1;
figure
idx = (data.prn == prn(i));
plot(data.elev(idx), data.obs(idx),  '-k', 'LineWidth',1/2)
ylabel('MP1 (m)')
xlabel('Elev. (degrees)')
title(sprintf('PRN %d', prn(i)))
xlim([32.5 50])

Contact us