image thumbnail

read_teqc_compact

by

 

17 Apr 2011 (Updated )

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

doit2.m
data = read_teqc_compact_each ('5145K79585201210170620C.sn5');

%%
prnsys = unique([data.prn double(data.sys)], 'rows');
prn = prnsys(:,1);
sys = char(prnsys(:,2));
figure
for i=1:size(prnsys,1)
  idx = ( (data.prn == prn(i)) & (data.sys == sys(i)) );
  plot(data.epoch(idx), data.obs(idx)-median(data.obs(idx)),  '.-k', 'LineWidth',1/2)
  ylabel('SNR var.'), xlabel('Time')
  title(sprintf('%c%d', sys(i), prn(i)))
  if (i ~= size(prnsys,1)),  pause;  end
end

Contact us